Inflow boundary condition for DNS of turbulent boundary layers on supersonic blunt cones

For direct numerical simulation (DNS) of turbulent boundary layers, gen- eration of an appropriate inflow condition needs to be considered. This paper proposes a method, with which the inflow condition for spatial-mode DNS of turbulent boundary layers on supersonic blunt cones with different Mach numbers, Reynolds numbers and wall temperature conditions can be generated. This is based only on a given instant flow field obtained by a temporal-mode DNS of a turbulent boundary layer on a flat plate. Effectiveness of the method is shown in three typical examples by comparing the results with those obtained by other methods.     

Modeling transonic weakly underexpanded turbulent jets

The results of the calculations of model and actual turbulent jet flows with shock waves at low supersonic Mach numbers are presented. The gasdynamic flow features characterizing shock reflection from a mixing layer are analyzed. A possible version of the modified model for the turbulent viscosity is proposed; the model makes it possible to improve the prediction of the shock (rarefaction wave) intensity distribution along jet flows.     

Experimental Study of Flow in a Cavity on an Axisymmetric Body

The flow past a cylindrical cavity on an axisymmetric body in the range of Mach numbers from 0.6 to 1.18 and the effect of the Mach number in the transition from subsonic to supersonic flow velocities are studied experimentally. In addition, a broad, 5.3—11.3 range of relative elongations of the cavity which permits one to determine the influence of the elongation on flow regimes including flows with closed and open separation zones is studied.     

Improvements to compressible Euler methods for low‐Mach number flows

In the present study improvements to numerical algorithms for the solution of the compressible Euler equations at low Mach numbers are investigated. To solve flow problems for a wide range of Mach numbers, from the incompressible limit to supersonic speeds, preconditioning techniques are frequently employed. On the other hand, one can achieve the same aim by using a suitably modified acoustic damping method. The solution algorithm presently under consideration is based on Roe's approximate Riemann solver [Roe PL. Approximate Riemann solvers, parameter vectors and difference schemes. Journal of Computational Physics 1981; 43 : 357–372] for non‐structured meshes. The numerical flux functions are modified by using Turkel's preconditioning technique proposed by Viozat [Implicit upwind schemes for low Mach number compressible flows. INRIA, Rapport de Recherche No. 3084, January 1997] for compressible Euler equations and by using a modified acoustic damping of the stabilization term proposed in the present study. These methods allow the compressible Euler equations at low‐Mach number flows to be solved, and they are consistent in time. The efficiency and accuracy of the proposed modifications have been assessed by comparison with experimental data and other numerical results in the literature. Copyright © 2000 John Wiley & Sons, Ltd.     

Van Leer格式在全速域范围的推广

通过对格式耗散项的修正将Van Leer格式推广至全速域流场求解范围.对格式耗散项的分析表明,在低马赫数流动情况下格式耗散项中不应包含声速项,以此为依据对Van Leer迎风分裂格式提出了耗散项的修正方法.结合对控制方程时间导数项的预处理,修正后的格式能够成功地模拟低速流动问题,同时在其他马赫数范围内也不损失格式的收敛...     

参考焓值法在高速槽道湍流中的修正

利用高来流马赫数为3, 5, 6, 7, 10的槽道湍流直接数值模拟(direct numerical simulations, DNS)数据, 评估和修正经典的参考焓值法. 研究表明在高来流马赫数槽道湍流中, 经典参考焓值法预测的壁面热流与DNS结果相差很大, 需要作适当的修正.修正参考焓值法Ⅰ和Ⅱ的预测结果明显优于经典参考焓值法;并且修正参考焓值法Ⅱ更加适用于高马赫数流动, 其壁面热流与DNS结果的相对误差在10%以内. 同时, 修正参考焓值法Ⅱ的普适性在超声速燃烧室隔离段热环境试验中得到了验证.     

Development of an Intermittency Equation for the Modeling of the Supersonic/Hypersonic Boundary Layer Flow Transition

An intermittency transport equation is developed in this study to model the laminar-turbulence boundary layer transition at supersonic and hypersonic conditions. The model takes into account the effects of different instability modes associated with the variations in Mach numbers. The model equation is based on the intermittency factor γ concept and couples with the well-known SST k–ω eddy-viscosity model in the solution procedures. The particular features of the present model approach are that: (1) the fluctuating kinetic energy k includes the non-turbulent, as well as turbulent fluctuations; (2) the proposed transport equation for the intermittency factor γ triggers the transition onset through a source term; (3) through the introduction of a new length scale normal to wall, the present model employs the local variables only avoiding the use of the integral parameters, like the boundary layer thickness δ, which are often cost-ineffective with the modern CFD methods; (4) in the fully turbulent region, the model retreats to SST model. This model is validated with a number of available experiments on boundary layer transition including the incompressible, supersonic and hypersonic flows past flat plates, straight/flared cones at zero incidences, etc. It is demonstrated that the present model can be successfully applied to the engineering calculations of a variety of aerodynamic flow transition with a reasonably wide range of Mach numbers.     

Explicit algebraic Reynolds stress model (EARSM) for compressible shear flows

We develop an explicit algebraic Reynolds stress model (EARSM) for high-speed compressible shear flows and validate the model with direct numerical simulation (DNS) data of homogeneous shear flow and experimental data of high-speed mixing-layers. Starting from a pressure–strain correlation model that incorporates compressibility effects, the weak-equilibrium assumption is invoked to derive the EARSM closure expression. The resulting closure is fully explicit and physically realizable and is a function of mean flow strain rate, rotation rate, turbulent kinetic energy, dissipation rate, and gradient Mach number. Homogeneous shear flow calculations show that the model captures the asymptotic behavior of DNS quite well. Linear EARSM calculations of a plane supersonic mixing-layer are performed, and comparison with experimental data shows good agreement. Salient results are agreement of streamwise velocity similarity profiles, mixing-layer spreading rates, and capturing the Langley curve trend.     

A hybrid pressure‐based solver for nonideal single‐phase fluid flows at all speeds

This paper describes the implementation of a numerical solver that is capable of simulating compressible flows of nonideal single‐phase fluids. The proposed method can be applied to arbitrary equations of state and is suitable for all Mach numbers. The pressure‐based solver uses the operator‐splitting technique and is based on the PISO/SIMPLE algorithm: the density, velocity, and temperature fields are predicted by solving the linearized versions of the balance equations using the convective fluxes from the previous iteration or time step. The overall mass continuity is ensured by solving the pressure equation derived from the continuity equation, the momentum equation, and the equation of state. Nonphysical oscillations of the numerical solution near discontinuities are damped using the Kurganov‐Tadmor/Kurganov‐Noelle‐Petrova (KT/KNP) scheme for convective fluxes. The solver was validated using different test cases, where analytical and/or numerical solutions are present or can be derived: (1) A convergent‐divergent nozzle with three different operating conditions; (2) the Riemann problem for the Peng‐Robinson equation of state; (3) the Riemann problem for the covolume equation of state; (4) the development of a laminar velocity profile in a circular pipe (also known as Poiseuille flow); (5) a laminar flow over a circular cylinder; (6) a subsonic flow over a backward‐facing step at low Reynolds numbers; (7) a transonic flow over the RAE 2822 airfoil; and (8) a supersonic flow around a blunt cylinder‐flare model. The spatial approximation order of the scheme is second order. The mesh convergence of the numerical solution was achieved for all cases. The accuracy order for highly compressible flows with discontinuities is close to first order and, for incompressible viscous flows, it is close to second order. The proposed solver is named rhoPimpleCentralFoam and is implemented in the open‐source CFD library OpenFOAM^®. For high speed flows, it shows a similar behavior as the KT/KNP schemes (implemented as rhoCentralFoam‐solver, Int. J. Numer. Meth. Fluids 2010), and for flows with small Mach numbers, it behaves like solvers that are based on the PISO/SIMPLE algorithm.     

A Mach‐uniform unstructured staggered grid method

A novel Mach‐uniform method to compute flows using unstructured staggered grids is discussed. The Mach‐uniform method is a generalization of the pressure‐correction approach for incompressible flows, and is valid for Mach numbers ranging from 0 (incompressible) to > 1 (supersonic). The primary variables (ρ u ,p and ρ) are updated sequentially. The grid consists of triangles. A staggered positioning of the variables is employed: the scalar variables are located at the centroids of the triangles, whereas the normal momentum components are positioned at the midpoints of the faces of the triangles. Discretization of the two‐dimensional flow equations on unstructured staggered grids is discussed. For the cell face fluxes there is a choice between first‐order upwind and central approximation. Flows around the NACA 0012 airfoil with freestream Mach numbers ranging from 0 to 1.2 are computed to demonstrate the Mach‐uniform accuracy and efficiency of the proposed method. Copyright © 2002 John Wiley & Sons, Ltd.     

Simulation of supersonic and hypersonic flows

A simple convection algorithm for simulation of time-dependent supersonic and hypersonic flows of a perfect but viscous gas is described. The algorithm is based on conservation and convection of mass, momentum and energy in a grid of rectangular cells. Examples are given for starting flow in a shock tube and oblique shocks generated by a wedge at Mach numbers up to 30·4. Good comparisons are achieved with well-known perfect gas flows.     

Features of a Turbulent Jet at High Supersonic Velocities

The results of modeling a turbulent supersonic jet at M = 5 using large-eddy simulation (LES) are presented. The structural features of turbulence formed in this flow are analyzed. The possibilities of the large-eddy simulation method and the complexities of simulation of the compressibility effects in jet flows at high Mach numbers are considered. Such features of the supersonic jet as the local turbulent shocklets and Mach waves are reproduced numerically. It is shown that in the neighborhood of the jet the trajectories of ejection flow are located along the front of Mach waves. Anisotropic turbulent structures whose longitudinal scale is greater than the transverse scale by an order of magnitude are revealed in the jet. An estimate of the baroclinic effects shows their weak influence on the vorticity generation in the jet flow considered.     

基于直接数值模拟的可压缩湍流模型评估和改进

对来流Mach数2.25和6的平板边界层湍流进行了直接数值模拟, 并通过与理论、实验及他人计算结果的对比对数值结果进行了验证. 基于直接数值模拟得到的湍流数据库, 对常用的湍流模型进行了先验评估. 评估的湍流模型有k-εvarepsilon模型(包括标准k-εvarepsilon 模型、可实现的k-εvarepsilon模型及低Reynolds数k-εvarepsilon模型)、SA模型及BL模型. 结果显示, 对于Mach2.25的平板边界层, 可实现的k-εvarepsilon 模型及低Reynolds 数k-εvarepsilon模型具有较好的预测能力, 而标准k-εvarepsilon模型预测的湍流黏性系数偏高; SA模型在边界层内层预测准确度较高, 而在外层预测值偏高. 而对于Mach6的平板边界层, k-εvarepsilon模型及SA模型预测的湍流黏性系数均偏高, 尤其是标准k-εvarepsilon模型. 对于Mach6的平板边界层, BL模型低估了内-外层交界位置, 造成湍流黏性系数预测值严重偏低. 作者通过修改模型系数及内-外层交界位置对BL模型进行了修改, 修改后模型预测的湍流黏性系数与DNS给出的值吻合较好.     

Characteristics of Oscillations in Supersonic Open Cavity Flows

Characteristics of oscillations in supersonic open cavity flows are investigated numerically using hybrid RANS/LES (Reynolds-Averaged Navier-Stokes/Large Eddy Simulation) method. The oscillation regimes and feedback mechanisms for the supersonic cavity flows are identified and analyzed. The calculation captures a mixed shear-layer/wake oscillation mode in the flow of Ma = 1.75, where these two modes occur alternately. The shear-layer mode and wake mode are driven by vortex convection-acoustic feedback and absolute instability, respectively. In particular, the results indicate that the feedback-acoustic-wave in the shear-layer mode is probably generated by the reflection of the downstream-traveling pressure wave, associated with the shed vortex in the shear layer, on the aft wall. The cavity flow of Ma = 2.52 is then simulated to see the influence of Mach number. It is found that the increase of Mach number may decrease the amplitude of the fluctuations in the shear layer, inhibiting the transition to wake mode. Furthermore, the influence of upstream injection is also studied, where the results show that the injection only weakens the oscillations and faintly shifts the resonant frequencies.     

A ghost-cell immersed boundary method for large-eddy simulations of compressible turbulent flows

A methodology to perform a ghost-cell-based immersed boundary method (GCIBM) is presented for simulating compressible turbulent flows around complex geometries. In this method, the boundary condition on the immersed boundary is enforced through the use of ‘ghost cells’ that are located inside the solid body. The computations of variables on these ghost cells are achieved using linear interpolation schemes. The validity and applicability of the proposed method is verified using a three-dimensional (3D) flow over a circular cylinder, and a large-eddy simulation of fully developed 3D turbulent flow in a channel with a wavy surface. The results agree well with the previous numerical and experimental results, given that the grid resolution is reasonably fine. To demonstrate the capability of the method for higher Mach numbers, supersonic turbulent flow over a circular cylinder is presented. While more work still needs to be done to demonstrate higher robustness and accuracy, the present work provides interesting insights using the GCIBM for the compressible flows.     

Starting of an axisymmetric convergent-divergent nozzle in a hypersonic flow

The starting of an axisymmetric convergent-divergent nozzle, with the result that supersonic flow is formed within almost the entire channel, is modeled, as applied to the hypersonic aerodynamic setup of the Institute of Mechanics of Moscow State University. A successful starting is realized when the nozzle is thrown in a uniform supersonic air flow at a fairly high Mach number. The steady flow structure is studied. It is numerically shown that in the convergent section of the channel there arises an oblique shock wave whose interaction with the nozzle axis leads to the formation of a reflected shock and a curvilinear Mach disk with a region of unsteady subsonic flow in the vicinity of the throat. The mathematical model is based on the two-dimensional Euler equations for axisymmetric gas flows.     

On the use of a two-layer model for large-eddy simulations of supersonic boundary layers with separation

The two-layer modeling approach has become one of the most promising and successful methodology for simulating turbulent boundary layers in the past ten years. In the present study, a mixed wall model for large-eddy simulations (LES) of high-speed flows is proposed which combine two approaches; the thin-Boundary Layer Equations (TBLE) model of Kawai and Larsson (1994) and the analytical wall-layer model of Duprat et al. (2011) for streamwise pressure gradients. The new hybrid model has been efficiently implemented into a three-dimensional compressible LES solver and validated against DNS of a spatially-evolving supersonic boundary layer (BL) under moderate and strong pressure gradients, before being employed for the prediction of nozzle flow separations at different flow conditions, ranging from weakly to highly over-expanded regimes. A good agreement is obtained in terms of mean and fluctuating quantities compared to the DNS results. Particularly, the current wall-modeled LES results are found to perfectly match the DNS data of supersonic BL with/out pressure gradient. It is also shown that the model can account for the effect of the large-scale turbulent motions of the outer layer, indicating a good interaction between the inner and the outer part of the wall layer. In terms of simulations costs and improvements of computing power, the obtained results highlight the capability of the current wall-modeling LES strategy in saving a considerable amount of computational time compared to the wall-resolved LES counterpart, allowing to push further the simulations limits. Furthermore, the application of these computationally low-costly LES simulations to nozzle flow separation allows to clearly identify the origin of the shock unsteadiness, and the existence of broadband and energetically-significant low-frequency oscillations (LFO) in the vicinity of the separation region.     

Slender axisymmetric cavities in a supersonic flow

Slender axisymmetric cavities in a subsonic flow of compressible fluid were investigated in [1–4]. In [5] a finite-difference method was used to calculate the drag coefficient of a circular cone, near which the shape of the cavity was determined for subsonic, transonic, and supersonic water flows; however, in the supersonic case the entire shape of the cavity was not determined. Here, on the basis of slender body theory an integrodifferential equation is obtained for the profile of the cavity in a supersonic flow. The dependence of the cavity elongation on the cavitation number and the Mach number is determined.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 179–181, January–February, 1989.     

Hyperbolic Approximation of the Navier-Stokes Equations for Viscous Mixed Flows

Simplified two-dimensional Navier-Stokes equations of the hyperbolic type are derived for viscous mixed (with transition through the sonic velocity) internal and external flows as a result of a special splitting of the pressure gradient in the predominant flow direction into hyperbolic and elliptic components. The application of these equations is illustrated with reference to the calculation of Laval nozzle flows and the problem of supersonic flow past blunt bodies. The hyperbolic approximation obtained adequately describes the interaction between the stream and surfaces for internal and external flows and can be used over a wide Mach number range at moderate and high Reynolds numbers. Examples of the calculation of viscous mixed flows in a Laval nozzle with large longitudinal throat curvature and in a shock layer in the neighborhood of a sphere and a large-aspect-ratio hemisphere-cylinder are given. The problem of determining the drag coefficient of cold and hot spheres is solved in a new formulation for supersonic air flow over a wide range of Reynolds numbers. In the case of low and moderate Reynolds numbers a drag reduction effect is detected when the surface of the sphere is cooled.     

Prediction of asymmetric vortical flows around slender bodies using Navier-Stokes equations

Steady and unsteady asymmetric vortical flows around slender bodies at high angles of attack are solved using the unsteady, compressible, this-layer Navier-Stokes equations. An implicit, upwind-biased, flux-difference splitting, finite-volume scheme is used for the numerical computations. For supersonic flows past point cones, the locally conical flow assumption has been used for efficient computational studies of this phenomenon. Asymmetric flows past a 5° semiapex-angle circular cone at different angles of attack, free-stream Mach numbers, and Reynolds numbers has been studied in responses to different sources of disturbances. The effects of grid fineness and computational domain size have also been investigated. Next, the responses of three-dimensional supersonic asymmetric flow around a 5° circular cone at different angles of attack and Reynolds numbers to short-duration sideslip disturbances are presented. The results show that flow asymmetry becomes stronger as the Reynolds number and angles of attack are increased. The asymmetric solutions show spatial vortex shedding which is qualitatively similar to the temporal vortex shedding of the unsteady locally conical flow. A cylindrical afterbody is also added to the same cone to study the effect of a cylindrical part on the flow asymmetry. One of the cases of flow over a cone-cylinder configuration is validated fairly well by experimental data.