一种新型的高分辨率通量差分裂格式

传统的Roe格式不满足熵条件并且在计算激波问题时会遭遇不同形式的不稳定现象,如慢行激波的波后振荡和红玉(carbuncle)现象.基于Zha-Bilgen对流-压力通量分裂方法,构造一种新型的通量差分裂格式.利用约旦标准型理论,通过添加广义特征向量构造通量差分裂方法来计算对流子系统.压力子系统具有一组完备的线性无关特征向量,因此可以构造传统的通量差分裂格式进行计算.为了提高接触间断的分辨率,利用界面变差下降(BVD)算法来重构对流通量耗散项中的密度差.激波稳定性分析表明,新格式可以有效地衰减数值误差,从而抑制不稳定现象的发生.一系列数值实验证明了本文构造的新型通量差分裂格式比Roe格式具有更高的分辨率和更好的鲁棒性.     

Numerical solutions of Euler equations by using a new flux vector splitting scheme

A new flux vector splitting scheme has been suggested in this paper. This scheme uses the velocity component normal to the volume interface as the characteristic speed and yields the vanishing individual mass flux at the stagnation. The numerical dissipation for the mass and momentum equations also vanishes with the Mach number approaching zero. One of the diffusive terms of the energy equation does not vanish. But the low numerical diffusion for viscous flows may be ensured by using higher-order differencing. The scheme is very simple and easy to be implemented. The scheme has been applied to solve the one dimensional (1D) and multidimensional Euler equations. The solutions are monotone and the normal shock wave profiles are crisp. For a 1D shock tube problem with the shock and the contact discontinuities, the present scheme and Roe scheme give very similar results, which are the best compared with those from Van Leer scheme and Liou–Steffen's advection upstream splitting method (AUSM) scheme. For the multidimensional transonic flows, the sharp monotone normal shock wave profiles with mostly one transition zone are obtained. The results are compared with those from Van Leer scheme, AUSM and also with the experiment.     

N-S方程在非结构网格下的求解

在Ｒｏｅ的矢通量差分分裂的基础上，吸收了ＮＮＤ格式的优点，提出了一种非结构网格下求解Ｅｕｌｅｒ方程和Ｎ－Ｓ方程的高分辨率高精度迎风格式．这种格式具有捕捉强激波和滑移线的良好性能．在时间方向上采用了显式和隐式两种解法．文中还给出了自适应技术．最后，成功地完成了ＧＡＭＭ超音速前台阶绕流、二维平板无粘激波反射、三维Ｈｏｂｓｏｎ叶栅流动、ＶＫＩ叶栅流动、Ｃ３Ｘ叶栅流动的数值模拟，得到了满意的结果     

Realization of contact resolving approximate Riemann solvers for strong shock and expansion flows

The Harten–Lax–van Leer contact (HLLC) and Roe schemes are good approximate Riemann solvers that have the ability to resolve shock, contact, and rarefaction waves. However, they can produce spurious solutions, called shock instabilities, in the vicinity of strong shock. In strong expansion flows, the Roe scheme can admit nonphysical solutions such as expansion shock, and it sometimes fails. We carefully examined both schemes and propose simple methods to prevent such problems. High‐order accuracy is achieved using the weighted average flux (WAF) and MUSCL‐Hancock schemes. Using the WAF scheme, the HLLC and Roe schemes can be expressed in similar form. The HLLC and Roe schemes are tested against Quirk's test problems, and shock instability appears in both schemes. To remedy shock instability, we propose a control method of flux difference across the contact and shear waves. To catch shock waves, an appropriate pressure sensing function is defined. Using the proposed method, shock instabilities are successfully controlled. For the Roe scheme, a modified Harten–Hyman entropy fix method using Harten–Lax–van Leer‐type switching is suggested. A suitable criterion for switching is established, and the modified Roe scheme works successfully with the suggested method. Copyright © 2009 John Wiley & Sons, Ltd.     

一种健壮的混合Roe黎曼求解器

使用Roe格式计算多维流动问题时,在强激波附近会出现数值激波不稳定现象。带有剪切粘性的HLLEC格式不仅可以捕捉接触间断,而且表现出很好的稳定性。混合Roe格式和HLLEC格式来消除数值激波不稳定性。在强激波附近,通过激波面法向和网格界面法向的夹角来定义开关函数,使得数值通量在激波面横向切换成HLLEC格式。在其余地方,数值通量依然使用Roe格式来计算。数值试验表明,混合格式不仅消除了Roe格式的数值激波不稳定性,还最大程度地减少了HLLEC格式所带来的剪切耗散,保留了Roe格式高分辨率的优点。     

变高度圆柱诱导的激波边界层干扰

采用分区方法及Ｒｏｅ三阶流通量差分分裂格式求解雷诺平均Ｎ－Ｓ方程,湍流附加黏性系数用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ模型计算,数值模拟了高超声速条件下变高度圆柱诱导的激波边界层层干扰,其流场的主要特性均与实验结果一致或规律相同,结果清晰地展示了流场结构以及气动载荷分布随柱高度的变化特征,产说明激波碰撞和旋涡运动都可能导致飞行器表面局部气动载荷的增加。     

Performance Improvement of a Supercritical Airfoil by a Multi-point Optimized Shock Control Channel

A shock control channel (SCC) is a flow control method introduced here to control the shock wave/boundarylayer interaction (SWBLI) in order to reduce the resulting wave drag in transonic flows. An SCC transfers an appropriate amount of mass and momentum from downstream of the shock wave location to its upstream to decrease the pressure gradient across the shock wave and as a result the shock-wave strength is reduced. Here, a multi-point optimization method under a constant-lift-coefficient constraint is used to find the optimum design of the SCC. This flow control method is implemented on a RAE-2822 supercritical airfoil for a wide range of off-design transonic Mach numbers. The RANS flow equations are solved using Roe’s averages scheme and a gradient-based adjoint algorithm is used to find the optimum location and shape of the SCC. The solver is validated against experimental works studying effect of cavities in transonic airfoils. It is shown that the application of an SCC improves the average aerodynamic efficiency in a range of off-design conditions by 13.2% in comparison with the original airfoil. The SCC is shown to be an effective tool also for higher angle of attack at transonic flows. We have also studied the SWBLI and how the optimization algorithm makes the flow wave structure and interactions of the shock wave with the boundary layer favorable.     

Performance of finite volume solutions to the shallow water equations with shock‐capturing schemes

Numerical methods have become well established as tools for solving problems in hydraulic engineering. In recent years the finite volume method (FVM) with shock capturing capabilities has come to the fore because of its suitability for modelling a variety of types of flow; subcritical and supercritical; steady and unsteady; continuous and discontinuous and its ability to handle complex topography easily. This paper is an assessment and comparison of the performance of finite volume solutions to the shallow water equations with the Riemann solvers; the Osher, HLL, HLLC, flux difference splitting (Roe) and flux vector splitting. In this paper implementation of the FVM including the Riemann solvers, slope limiters and methods used for achieving second order accuracy are described explicitly step by step. The performance of the numerical methods has been investigated by applying them to a number of examples from the literature, providing both comparison of the schemes with each other and with published results. The assessment of each method is based on five criteria; ease of implementation, accuracy, applicability, numerical stability and simulation time. Finally, results, discussion, conclusions and recommendations for further work are presented. Copyright © 2002 John Wiley & Sons, Ltd.     

Evaluation of Euler fluxes by a high-order CFD scheme: shock instability

The construction of Euler fluxes is an important step in shock-capturing/upwind schemes. It is well known that unsuitable fluxes are responsible for many shock anomalies, such as the carbuncle phenomenon. Three kinds of flux vector splittings (FVSs) as well as three kinds of flux difference splittings (FDSs) are evaluated for the shock instability by a fifth-order weighted compact nonlinear scheme. The three FVSs are Steger–Warming splitting, van Leer splitting and kinetic flux vector splitting (KFVS). The three FDSs are Roe's splitting, advection upstream splitting method (AUSM) type splitting and Harten–Lax–van Leer (HLL) type splitting. Numerical results indicate that FVSs and high dissipative FDSs undergo a relative lower risk on the shock instability than that of low dissipative FDSs. However, none of the fluxes evaluated in the present study can entirely avoid the shock instability. Generally, the shock instability may be caused by any of the following factors: low dissipation, high Mach number, unsuitable grid distribution, large grid aspect ratio, and the relative shock-internal flow state (or position) between upstream and downstream shock waves. It comes out that the most important factor is the relative shock-internal state. If the shock-internal state is closer to the downstream state, the computation is at higher susceptibility to the shock instability. Wall-normal grid distribution has a greater influence on the shock instability than wall-azimuthal grid distribution because wall-normal grids directly impact on the shock-internal position. High shock intensity poses a high risk on the shock instability, but its influence is not as much as the shock-internal state. Large grid aspect ratio is also a source of the shock instability. Some results of a second-order scheme and a first-order scheme are also given. The comparison between the high-order scheme and the two low-order schemes indicates that high-order schemes are at a higher risk of the shock instability. Adding an entropy fix is very helpful in suppressing the shock instability for the two low-order schemes. When the high-order scheme is used, the entropy fix still works well for Roe's flux, but its effect on the Steger–Warming flux is trivial and not much clear.     

UNSTEADY／STEADY NUMERICAL SIMULATION OF THREE-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS ON ARTIFICIAL COMPRESSIBILITY

A mixed algorithm of central and upwind difference scheme for the solution of steady/unsteady incompressible Navier-Stokes equations is presented. The algorithm is based on the method of artificial compressibility and uses a third-order flux-difference splitting technique for the convective terms and the second-order central difference for the viscous terms. The numerical flux of semi-discrete equations is computed by using the Roe approximation. Time accuracy is obtained in the numerical solutions by subiterating the equations in pseudotime for each physical time step. The algebraic turbulence model of Baldwin-Lomax is ulsed in this work. As examples, the solutions of flow through two dimensional flat, airfoil, prolate spheroid and cerebral aneurysm are computed and the results are compared with experimental data. The results show that the coefficient of pressure and skin friction are agreement with experimental data, the largest discrepancy occur in the separation region where the lagebraic turbulence model of Baldwin-Lomax could not exactly predict the flow.     

The Multi-point Optimization of Shock Control Bump with Constant-Lift Constraint Enhanced with Suction and Blowing for a Supercritical Airfoil

Both shock control bump (SCB) and suction and blowing are flow control methods used to control the shock wave/boundary layer interaction (SWBLI) in order to reduce the resulting wave drag in transonic flows. A SCB uses a small local surface deformation to reduce the shock-wave strength, while suction decreases the boundary-layer thickness and blowing delays the flow separation. Here a multi-point optimization method under a constant-lift-coefficient constraint is used to find the optimum design of SCB and suction and blowing. These flow control methods are used separately or together on a RAE-2822 supercritical airfoil for a wide range of off-design transonic Mach numbers. The RANS flow equations are solved using Roe’s averages scheme and a gradient-based adjoint algorithm is used to find the optimum location and shape of all devices. It is shown that the simultaneous application of blowing and SCB (hybrid blowing/SCB) improves the average aerodynamic efficiency at off-design conditions by 18.2 % in comparison with the clean airfoil, while this increase is only 16.9 % for the hybrid suction/SCB. We have also studied the SWBLI and how the optimization algorithm makes the flow wave structure and interactions of the shock wave with the boundary layer favorable.     

用隐式方法和WENO格式计算悬停旋翼跨声速无粘流场

寻找一种能够准确计算以涡为主要特征的复杂流场和克服尾迹耗散问题的数值方法,一直是旋翼空气动力学研究的热点和难点。本文发展了一种基于高阶迎风格式计算悬停旋翼无粘流场的隐式数值方法。无粘通量采用Roe通量差分分裂格式,为提高精度,使用五阶WENO格式进行左右状态插值,并与MUSCL插值进行比较。为提高收敛到定常解的效率,时间推进采用LU-SGS隐式方法。用该方法对一跨声速悬停旋翼无粘流场进行了数值计算,数值结果表明WENO-Roe的激波分辨率高于MUSCL-Roe,体现出了格式精度的提高对计算结果的改善,LU-SGS隐式方法的计算效率比5步Runge-Kutta显式方法的高。     

一种激波稳定的对流-压力通量分裂格式

针对欧拉方程三种流行的对流-压力通量分裂方法(Liou-Steffen,Zha-Bilgen和Toro-Vázquez)进行特征分析,进而提出一种新的对流-压力通量分裂格式。采用Zha-Bilgen分裂方法将欧拉方程的通量分裂成对流项和压力项两部分,使用TV格式来计算这两部分的数值通量。利用压力比构造激波探测函数,并且在强激波附近的亚声速区域增加TV格式的剪切粘性来克服数值模拟中的激波不稳定性。数值算例的计算结果表明,新的对流-压力通量分裂格式不仅保留了原始TV格式精确分辨接触间断的优点,而且具有更好的鲁棒性,在数值模拟多维强激波问题时不会出现不稳定现象。因此,该格式是一种精确并且具有强鲁棒性的数值方法,可以广泛地应用于可压缩流体的数值计算中。     

Solution to Euler equations by high‐resolution upwind compact scheme based on flux splitting

A high‐resolution upwind compact method based on flux splitting is developed for solving the compressive Euler equations. The convective flux terms are discretized by using the modified advection upstream splitting method (AUSM). The developed scheme is used to compute the one‐dimensional Burgers equation and four different example problems of supersonic compressible flows, respectively. The results show that the high‐resolution upwind compact scheme based on modified AUSM+ flux splitting can capture shock wave and other discontinuities, obtain higher resolution and restrain numerical oscillation. Copyright © 2007 John Wiley & Sons, Ltd.     

求解可压缩流的二维通量分裂格式

传统的一维通量分裂格式在计算界面数值通量时,只考虑网格界面法向的波系。采用传统的TV格式分别求解对流通量和压力通量。通过求解考虑了横向波系影响的角点数值通量来构造一种真正二维的TV通量分裂格式。在计算一维数值算例时,该格式与传统的TV格式具有相同的数值通量计算公式,因此其保留了传统的TV格式精确捕捉接触间断和膨胀激波的优点。在计算二维算例时,该格式比传统的TV格式具有更高的分辨率;在计算二维强激波问题时,消除了传统TV格式的非物理现象,表现出更好的鲁棒性;此外,该格式大大提高了稳定性CFL数,从而具有更高的计算效率。因此,本文方法是一种精确、高效并且具有强鲁棒性的数值方法,在可压缩流的数值模拟中具有广阔的应用前景。     

The region of nonlinear effects for intensive sound pulses in the ocean

The formation of a stationary shock wave is studied in media with an arbitrary power dependence of the damping coefficient on the frequency. The conditions for existence of a stationary shock wave are defined and it is shown that when acoustic signals propagate in the ocean the region of nonlinear effects is limited. For acoustic waves generated by explosive sources a calculation is given of the location of the transition point of the nonlinear wave into a linear one, and the dependence of this point on the charge weight is defined.     

Flux vector splitting solutions for coupling hydraulic transient of gas-liquid-solid three-phase flow in pipelines

The gas-liquid-solid three-phase mixed flow is the most general in multiphase mixed transportation. It is significant to exactly solve the coupling hydraulic transient problems of this type of multiphase mixed flow in pipelines. Presently, the method of characteristics is widely used to solve classical hydraulic transient problems. However, when it is used to solve coupling hydraulic transient problems, excessive interpolation errors may be introduced into the results due to unavoidable multiwave interpolated calculations. To deal with the problem, a finite difference scheme based on the Steger-Warming flux vector splitting is proposed. A flux vector splitting scheme is established for the coupling hydraulic transient model of gas-liquid-solid three-phase mixed flow in the pipelines. The flux subvectors are then discretized by the Lax-Wendroff central difference scheme and the Warming-Beam upwind difference scheme with second-order precision in both time and space. Under the Rankine-Hugoniot conditions and the corresponding boundary conditions, an effective solution to those points located at the boundaries is developed, which can avoid the problem beyond the calculation region directly induced by the second-order discrete technique. Numerical and experimental verifications indicate that the proposed scheme has several desirable advantages including high calculation precision, excellent shock wave capture capability without false numerical oscillation, low sensitivity to the Courant number, and good stability.     

A holographic interferometric study of shock wave focusing in a circular reflector

A holographic interferometric study was made of the focusing of reflected shock waves from a circular reflector. A diaphragmless shock tube was used for incident shock Mach numbers ranging from 1.03 to 1.74. Hence, the process of reflected shock wave focusing was quantitatively observed. It is found that a converging shock wave along the curved wall undergoes an unsteady evolution of mach reflection and its focusing is, therefore, subject to the evolution of the process of shock wave reflections. The collision of triple points terminates the focusing process at the geometrical focus. In order to interprete quantitatively these interferograms, a numerical simulation using an Eulerian solver combined with adaptive unstructured grids was carried out. It is found numerically that the highest density appears immediately after the triple point collision. This implies that the final stage of focusing is mainly determined by the interaction between shock waves and vortices. The interaction of finite strength shock waves, hence, prevents a curved shock wave from creating the infinite increase of density or pressure at a focal point which is otherwise predicted by the linear acoustic theory.     

An upstream flux‐splitting finite‐volume scheme for 2D shallow water equations

An upstream flux‐splitting finite‐volume (UFF) scheme is proposed for the solutions of the 2D shallow water equations. In the framework of the finite‐volume method, the artificially upstream flux vector splitting method is employed to establish the numerical flux function for the local Riemann problem. Based on this algorithm, an UFF scheme without Jacobian matrix operation is developed. The proposed scheme satisfying entropy condition is extended to be second‐order‐accurate using the MUSCL approach. The proposed UFF scheme and its second‐order extension are verified through the simulations of four shallow water problems, including the 1D idealized dam breaking, the oblique hydraulic jump, the circular dam breaking, and the dam‐break experiment with 45° bend channel. Meanwhile, the numerical performance of the UFF scheme is compared with those of three well‐known upwind schemes, namely the Osher, Roe, and HLL schemes. It is demonstrated that the proposed scheme performs remarkably well for shallow water flows. The simulated results also show that the UFF scheme has superior overall numerical performances among the schemes tested. Copyright © 2005 John Wiley & Sons, Ltd.     

Modelling of turbulent energy flux in canonical shock-turbulence interaction

High-speed turbulent flows often encounter high heat loads due to the presence of shock waves. The turbulent energy flux correlation in the mean energy conservation equation is a key unclosed term that determines the heat transfer rate. In this work, we employ existing turbulence models to predict the turbulent energy flux in canonical shock-turbulence interaction. The shortcomings of these models are highlighted, and a new heat-flux limiter model is proposed with the aid of linear theory results. We also write the transport equation for the turbulent energy flux across a shock wave and use it to develop a physics-based model for the same. It is found to predict the peak energy flux at the shock wave and its variation in the acoustic-adjustment region behind the shock. Numerical error incurred while solving the model equations at a shock wave are analyzed and a numerically robust model is obtained by eliminating the nonconservative source terms. The model predictions are compared with available direct numerical simulation data and a good match is obtained for a range of Mach numbers.