8阶群速度控制格式及其应用

构造了8阶精度的群速度控制型差分格式.根据激波捕捉能力及对小尺度分辨能力对格式系数进行了优化.采用该格式对可压缩均匀各向同性湍流进行了直接数值模拟,所计算的最大湍流马赫数达到0.95.相同计算条件下的结果与他人结果吻合较好,说明了格式在捕捉激波的同时对湍流小尺度有较好的分辨能力.     

可压缩各向同性衰减湍流直接数值模拟研究

采用五阶有限差分WENO格式直接模拟了高初始湍流Mach数的可压缩均匀各向同性湍流,主要分析了湍流的统计特性 和压缩性的影响,包括能谱特征、激波串、耗散率、标度律等. 研究表明,湍动能主要来自于速度场螺旋分量的贡献;各向同性湍流的小尺度脉动对压缩性更为敏感,并且压缩性的增强加快了湍流大 尺度脉动向小尺度脉动的湍动能输运;随着湍流Mach数的升高,胀量(压缩)耗散率所占比率也显著增长. 标度律分析表明,强可压缩湍流的横向速度结构函数仍然具有扩展自相似性;当阶数较高(p ≥ 5)时,纵向速度结构函数的扩展自相似性则不再成立. 对于压缩性较弱的湍流,与不可压缩湍流一致,横向湍流脉动的间歇性要强于纵向湍流脉动;而对于强可压缩湍流,纵向湍流脉动的 间歇性要强于横向湍流脉动.     

GPU‐accelerated direct numerical simulations of decaying compressible turbulence employing a GKM‐based solver

Gas Kinetic Method‐based flow solvers have become popular in recent years owing to their robustness in simulating high Mach number compressible flows. We evaluate the performance of the newly developed analytical gas kinetic method (AGKM) by Xuan et al. in performing direct numerical simulation of canonical compressible turbulent flow on graphical processing unit (GPU)s. We find that for a range of turbulent Mach numbers, AGKM results shows excellent agreement with high order accurate results obtained with traditional Navier–Stokes solvers in terms of key turbulence statistics. Further, AGKM is found to be more efficient as compared with the traditional gas kinetic method for GPU implementation. We present a brief overview of the optimizations performed on NVIDIA K20 GPU and show that GPU optimizations boost the speedup up‐to 40x as compared with single core CPU computations. Hence, AGKM can be used as an efficient method for performing fast and accurate direct numerical simulations of compressible turbulent flows on simple GPU‐based workstations. Copyright © 2016 John Wiley & Sons, Ltd.     

On the nature of multitude scalings in decaying isotropic turbulence

We report multitude scaling laws for isotropic fully developed decaying turbulence through group theoretic method employing on the spectral equations both for modelling and without any modelling of nonlinear energy transfer. For modelling, besides the existence of classical power law scalings, an exponential decay of turbulent energy in time is obtained subject to exponentially decaying integral length scale at infinite Reynolds number limit. For the transfer without modelling, at finite Reynolds number, in addition to general power law decay of turbulence intensity with integral length scale growing as a square root of time, an exponential decay of energy in time is explored when integral length scale remains constant. Both the power and exponential decaying laws of energy agree to the theoretical results of George (1992), George and Wang (2009) and experimental results of fractal grid generated turbulence by Hurst and Vassilicos (2007). At infinite Reynolds number limit, a general power law scaling is obtained from which all classical scaling laws are recovered. Further, in this limit, turbulence exhibits a general exponential decaying law of energy with exponential decaying integral length scale depending on two scaling group parameters. The role of symmetry group parameters on turbulence dynamics is discussed in this study.     

Flow-thermodynamics interactions in decaying anisotropic compressible turbulence with imposed temperature fluctuations

The fundamental nature of the non-linear flow-thermodynamics interactions in a compressible turbulent flow with imposed temperature fluctuations is investigated. Direct numerical simulations (DNS) of decaying anisotropic compressible turbulence (turbulent Mach number 0.06–0.6) with imposed temperature fluctuations are performed to examine: (i) interactions between solenoidal and dilatational kinetic energy; (ii) partition between dilatational kinetic energy and thermodynamic potential energy; and (iii) redistribution of solenoidal and dilatational kinetic energy among the various Reynolds stress components. It is found that solenoidal kinetic energy levels and return-to-isotropy are weakly dependent on Mach number but independent of imposed temperature fluctuations in the parameter range studied. The dilatational kinetic energy generated is proportional to the square of the pressure fluctuations associated with the initial solenoidal and temperature fluctuations and thus a strong function of Mach number and heat release intensity. The energy exchange between dilatational kinetic and potential energy is driven by a strong proclivity toward equipartition. Consequently, the dynamics of pressure-dilatation ( ${\overline{pd}}$ ), which is the mechanism of this energy exchange between dilatational and potential energies, is dictated entirely by the requirement to impose energy equipartition. Based on the results, we provide a physical picture of the solenoidal–dilatational–potential energy interactions and the action of pressure-dilatation. The identification of the fundamental precepts underlying the various interactions is of great utility for turbulence closure model development.     

Analyzing the influence of compressibility on the rapid pressure–strain rate correlation in turbulent shear flows

The influence of compressibility on the rapid pressure–strain rate tensor is investigated using the Green’s function for the wave equation governing pressure fluctuations in compressible homogeneous shear flow. The solution for the Green’s function is obtained as a combination of parabolic cylinder functions; it is oscillatory with monotonically increasing frequency and decreasing amplitude at large times, and anisotropic in wave-vector space. The Green’s function depends explicitly on the turbulent Mach number M _t , given by the root mean square turbulent velocity fluctuations divided by the speed of sound, and the gradient Mach number M _g , which is the mean shear rate times the transverse integral scale of the turbulence divided by the speed of sound. Assuming a form for the temporal decorrelation of velocity fluctuations brought about by the turbulence, the rapid pressure–strain rate tensor is expressed exactly in terms of the energy (or Reynolds stress) spectrum tensor and the time integral of the Green’s function times a decaying exponential. A model for the energy spectrum tensor linear in Reynolds stress anisotropies and in mean shear is assumed for closure. The expression for the rapid pressure–strain correlation is evaluated using parameters applicable to a mixing layer and a boundary layer. It is found that for the same range of M _t there is a large reduction of the pressure–strain correlation in the mixing layer but not in the boundary layer. Implications for compressible turbulence modeling are also explored.      

L2Roe: a low dissipation version of Roe's approximate Riemann solver for low Mach numbers

A modification of the Roe scheme called L²Roe for low dissipation low Mach Roe is presented. It reduces the dissipation of kinetic energy at the highest resolved wave numbers in a low Mach number test case of decaying isotropic turbulence. This is achieved by scaling the jumps in all discrete velocity components within the numerical flux function. An asymptotic analysis is used to show the correct pressure scaling at low Mach numbers and to identify the reduced numerical dissipation in that regime. Furthermore, the analysis allows a comparison with two other schemes that employ different scaling of discrete velocity jumps, namely, LMRoe and a method of Thornber et al. To this end, we present for the first time an asymptotic analysis of the last method. Numerical tests on cases ranging from low Mach number (M_∞=0.001) to hypersonic (M_∞=5) viscous flows are used to illustrate the differences between the methods and to show the correct behavior of L²Roe. No conflict is observed between the reduced numerical dissipation and the accuracy or stability of the scheme in any of the investigated test cases. Copyright © 2015 John Wiley & Sons, Ltd.     

Comparison of magnetic resonance concentration measurements in water to temperature measurements in compressible air flows

Magnetic resonance imaging (MRI) measurements in liquid flows provide highly detailed 3D mean velocity and concentration data in complex turbulent mixing flow applications. The scalar transport analogy is applied to infer the mean temperature distribution in high speed gas flows directly from the MRI concentration measurements in liquid. Compressibility effects on turbulent mixing are known to be weak for simple flows at high subsonic Mach number, and it was not known if this would hold in more complex flows characteristic of practical applications. Furthermore, the MRI measurements are often done at lower Reynolds number than the compressible application, although both are generally done in fully turbulent flows. The hypothesis is that the conclusions from MRI measurements performed in water are transferable to high subsonic Mach number applications. The present experiment is designed to compare stagnation temperature measurements in high speed airflow (M = 0.7) to concentration measurements in an identical water flow apparatus. The flow configuration was a low aspect ratio wall jet with a thick splitter plate producing a 3D complex downstream flow mixing the wall-jet fluid with the mainstream flow. The three-dimensional velocity field is documented using magnetic resonance velocimetry in the water experiment, and the mixing is quantified by measuring the mean concentration distribution of wall-jet fluid marked with dissolved copper sulfate. The airflow experiments are operated with a temperature difference between the main stream and the wall jet. Profiles of the stagnation temperature are measured with a shielded thermocouple probe. The results show excellent agreement between normalized temperature and concentration profiles after correction of the temperature measurements for the effects of energy separation. The agreement is within 1 % near the edges of the mixing layer, which suggests that the mixing characteristics of the large scale turbulence structures are the same in the two flows.     

Direct numerical simulation of compressible turbulent flows

This paper reviews the authors＇ recent studies on compressible turbulence by using direct numerical simulation （DNS）,including DNS of isotropic（decaying） turbulence, turbulent mixing-layer,turbulent boundary-layer and shock/boundary-layer interaction.Turbulence statistics, compressibility effects,turbulent kinetic energy budget and coherent structures are studied based on the DNS data.The mechanism of sound source in turbulent flows is also analyzed. It shows that DNS is a powerful tool for the mechanistic study of compressible turbulence.     

Computational study of unsteady turbulent flows around oscillating and ramping aerofoils

The aim of this work is to computationally investigate subsonic and transonic turbulent flows around oscillating and ramping aerofoils under dynamic‐stall conditions. The investigation is based on a high‐resolution Godunov‐type method and several turbulence closures. The Navier–Stokes and turbulence transport equations are solved in a strongly coupled fashion via an implicit‐unfactored scheme. We present results from several computations of flows around oscillating and ramping aerofoils at various conditions in order to (i) assess the accuracy of different turbulence models and (ii) contribute towards a better understanding of dynamic‐stall flows. The results show that the employed non‐linear eddy‐viscosity model generally improves the accuracy of the computations compared to linear models, but at low incidence angles the Spalart–Allmaras one‐equation model was found to provide adequate results. Further, the computations reveal strong similarities between laminar and high‐Reynolds number dynamic‐stall flows as well as between ramping and oscillating aerofoil cases. Investigation of the Mach number effects on dynamic‐stall reveals a delay of the stall angle within a range of Mach numbers. Investigation of the reduced frequency effects suggests the existence of an (almost) linear variation between pitch rate and stall angle, with higher slope at lower pitch rates. The pitch rate affects both the onset of dynamic‐stall as well as the evolution of the associated vortical structures. Copyright © 2003 John Wiley & Sons, Ltd.     

A finite volume approach to large eddy simulation of compressible,homogeneous, isotropic,decaying turbulence

Large eddy simulation of compressible, homogeneous, isotropic, decaying turbulence in a rectangular box is performed using finite volume techniques. An analysis of the energy spectra obtained from the simulations shows that an agreement with the Kolmogorov law for the inertial range is found only when an appropriate spatial discretization method is used. This agreement is obtained both for a low (0·05) and a moderate (0·6) Mach number when Smagorinsky's subgrid model is employed.     

Interaction of isotropic turbulence with a normal shock wave

The evolution of incompressible and compressible isotropic 2-d turbulent fields interacting with a normal shock wave up to Mach numbers of 2.4 was investigated by means of direct numerical simulation using an ENO scheme. A comparison of statistics with linear analysis results is presented. Vorticity amplification in the DNS agrees well with the linear theory. Energy spectra are enhanced more in the small scales than in the large scales for incoming incompressible turbulence. The amplification rate for initially compressible turbulence is comparatively small.     

A high‐order discontinuous Galerkin method for all‐speed flows

In this article, we present a discontinuous Galerkin (DG) method designed to improve the accuracy and efficiency of steady solutions of the compressible fully coupled Reynolds‐averaged Navier–Stokes and k ? ω turbulence model equations for solving all‐speed flows. The system of equations is iterated to steady state by means of an implicit scheme. The DG solution is extended to the incompressible limit by implementing a low Mach number preconditioning technique. A full preconditioning approach is adopted, which modifies both the unsteady terms of the governing equations and the dissipative term of the numerical flux function by means of a new preconditioner, on the basis of a modified version of Turkel's preconditioning matrix. At sonic speed the preconditioner reduces to the identity matrix thus recovering the non‐preconditioned DG discretization. An artificial viscosity term is added to the DG discretized equations to stabilize the solution in the presence of shocks when piecewise approximations of order of accuracy higher than 1 are used. Moreover, several rescaling techniques are implemented in order to overcome ill‐conditioning problems that, in addition to the low Mach number stiffness, can limit the performance of the flow solver. These approaches, through a proper manipulation of the governing equations, reduce unbalances between residuals as a result of the dependence on the size of elements in the computational mesh and because of the inherent differences between turbulent and mean‐flow variables, influencing both the evolution of the Courant Friedrichs Lewy (CFL) number and the inexact solution of the linear systems. The performance of the method is demonstrated by solving three turbulent aerodynamic test cases: the flat plate, the L1T2 high‐lift configuration and the RAE2822 airfoil (Case 9). The computations are performed at different Mach numbers using various degrees of polynomial approximations to analyze the influence of the proposed numerical strategies on the accuracy, efficiency and robustness of a high‐order DG solver at different flow regimes. Copyright © 2014 John Wiley & Sons, Ltd.     

The Numerical Decomposition of Turbulent Fluctuations in a Compressible Boundary Layer

In many flows the turbulence is weakly compressible even at large Mach number. For example, in a compressible boundary layer Ma<5, the differences relative to an incompressible boundary layer understood as being caused by density variations that accompany variations temperature across the layer. Turbulent fluctuations in a boundary layer are therefore expected to be dominated by the effects nonconstant temperature, and low Mach number theories in which fluctuations are not dominant should be applicable to the fluctuating field. However, the analysis of compressible boundary layer DNS data reveals presence of significant acoustic fluctuations. To distinguish acoustic and thermal effects, a numerical decomposition procedure compressible boundary layer fluctuations is applied to determine the and nonacoustic fluctuations. Except for very near the wall, where decomposition procedure is not valid, it is found that the fluctuations are only weakly coupled to the acoustic fluctuations at numbers as high as 6. Received 13 March 2000 and accepted 21 May 2001     

Effects of Compressibility and Shock-Wave Interactions on Turbulent Shear Flows

Compressibility effects are present in many practical turbulent flows, ranging from shock-wave/boundary-layer interactions on the wings of aircraft operating in the transonic flight regime to supersonic and hypersonic engine intake flows. Besides shock wave interactions, compressible flows have additional dilatational effects and, due to the finite sound speed, pressure fluctuations are localized and modified relative to incompressible turbulent flows. Such changes can be highly significant, for example the growth rates of mixing layers and turbulent spots are reduced by factors of more than three at high Mach number. The present contribution contains a combination of review and original material. We first review some of the basic effects of compressibility on canonical turbulent flows and attempt to rationalise the differing effects of Mach number in different flows using a flow instability concept. We then turn our attention to shock-wave/boundary-layer interactions, reviewing recent progress for cases where strong interactions lead to separated flow zones and where a simplified spanwise-homogeneous problem is amenable to numerical simulation. This has led to improved understanding, in particular of the origin of low-frequency behaviour of the shock wave and shown how this is coupled to the separation bubble. Finally, we consider a class of problems including side walls that is becoming amenable to simulation. Direct effects of shock waves, due to their penetration into the outer part of the boundary layer, are observed, as well as indirect effects due to the high convective Mach number of the shock-induced separation zone. It is noted in particular how shock-induced turning of the detached shear layer results in strong localized damping of turbulence kinetic energy.     

Direct Numerical Simulation of Statistically Stationary One- and Two-Phase Turbulent Combustion: A Turbulent Injection Procedure

A new turbulent injection procedure dedicated to fully compressible direct numerical simulation (DNS) or large eddy simulation (LES) solvers is proposed. To avoid the appearance of spurious acoustic waves, this method is based on an accurate tracking of the turbulent structures crossing the boundary at the inlet of the domain. A finite difference DNS solver has been coupled with a spectral simulation in which a statistically stationary homogeneous turbulence evolves to provide fluctuating boundary conditions.A new turbulence forcing method, dedicated to spectral solvers, has been developed as well to control the major properties of the injected flow (turbulent kinetic energy, dissipation rate and integral length scale). One-dimensional Navier–Stokes characteristic boundary conditions extended to non-stationary flows are coupled with the injection procedure to evaluate is potential in four various configurations: spatially decaying turbulence, dispersion of vaporizing sprays, propagation of one- and two-phase V-shape turbulent flames.     

An immersed boundary/wall modeling method for RANS simulation of compressible turbulent flows

In this paper, an immersed boundary (IB) method is developed to simulate compressible turbulent flows governed by the Reynolds‐averaged Navier‐Stokes equations. The flow variables at the IB nodes (interior nodes in the immediate vicinity of the solid wall) are evaluated via linear interpolation in the normal direction to close the discrete form of the governing equations. An adaptive wall function and a 2‐layer wall model are introduced to reduce the near‐wall mesh density required by the high resolution of the turbulent boundary layers. The wall shear stress modified by the wall modeling technique and the no‐penetration condition are enforced to evaluate the velocity at an IB node. The pressure and temperature at an IB node are obtained via the local simplified momentum equation and the Crocco‐Busemann relation, respectively. The SST k ? ω and S‐A turbulence models are adopted in the framework of the present IB approach. For the Shear‐Stress Transport (SST) k ? ω model, analytical solutions in near‐wall region are utilized to enforce the boundary conditions of the turbulence equations and evaluate the turbulence variables at an IB node. For the S‐A model, the turbulence variable at an IB node is calculated by using the near‐wall profile of the eddy viscosity. In order to validate the present IB approach, numerical experiments for compressible turbulent flows over stationary and moving bodies have been performed. The predictions show good agreements with the referenced experimental data and numerical results.     

Modelling of pressure–strain correlation in compressible turbulent flow

Previous studies carried out in the early 1990s conjectured that the main compressible effects could be associated with the dilatational effects of velocity fluctuation. Later, it was shown that the main compressibility effect came from the reduced pressure-strain term due to reduced pressure fluctuations. Although better understanding of the compressible turbulence is generally achieved with the increased DNS and experimental research effort, there are still some discrepancies among these recent findings. Analysis of the DNS and experimental data suggests that some of the discrepancies are apparent if the compressible effect is related to the turbulent Mach number, Mt. From the comparison of two classes of compressible flow, homogenous shear flow and inhomogeneous shear flow （mixing layer）, we found that the effect of compressibility on both classes of shear flow can be characterized in three categories corresponding to three regions of turbulent Mach numbers： the low-Mr, the moderate-Mr and high-Mr regions. In these three regions the effect of compressibility on the growth rate of the turbulent mixing layer thickness is rather different. A simple approach to the reduced pressure-strain effect may not necessarily reduce the mixing-layer growth rate, and may even cause an increase in the growth rate. The present work develops a new second-moment model for the compressible turbulence through the introduction of some blending functions of Mt to account for the compressibility effects on the flow. The model has been successfully applied to the compressible mixing layers.