HIGH-ORDER WENO SIMULATIONS OF THREE-DIMENSIONAL RESHOCKED RICHTMYER-MESHKOV INSTABILITY TO LATE TIMES:DYNAMICS,DEPENDENCE ON INITIAL CONDITIONS,AND COMPARISONS TO EXPERIMENTAL DATA

The dynamics of the reshocked multi-mode Richtmyer-Meshkov instability is investigated using 513×257 2three-dimensional ninth-order weighted essentially nonoscillatory shock-capturing simulations.A two-mode initial perturbation with superposed random noise is used to model the Mach 1.5 air/SF6 Vetter-Sturtevant shock tube experiment. The mass fraction and enstrophy isosurfaces,and density cross-sections are utilized to show the detailed flow structure before,during,and after reshock.It is shown that the mixing layer growth agrees well with the experimentally measured growth rate before and after reshock.The post-reshock growth rate is also in good agreement with the prediction of the Mikaelian model.A parametric study of the sensitivity of the layer growth to the choice of amplitudes of the short and long wavelength initial interfacial perturbation is also presented.Finally,the amplification effects of reshock are quantified using the evolution of the turbulent kinetic energy and turbulent enstrophy spectra,as well as the evolution of the baroclinic enstrophy production,buoyancy production,and shear production terms in the enstrophy and turbulent kinetic transport equations.     

Investigation of physiological pulsatile flow in a model arterial stenosis using large-eddy and direct numerical simulations

Physiological pulsatile flow in a 3D model of arterial stenosis is investigated by using large eddy simulation (LES) technique. The computational domain chosen is a simple channel with a biological type stenosis formed eccentrically on the top wall. The physiological pulsation is generated at the inlet using the first harmonic of the Fourier series of pressure pulse. In LES, the large scale flows are resolved fully while the unresolved subgrid scale (SGS) motions are modelled using a localized dynamic model. Due to the narrowing of artery the pulsatile flow becomes transition-to-turbulent in the downstream region of the stenosis, where a high level of turbulent fluctuations is achieved, and some detailed information about the nature of these fluctuations are revealed through the investigation of the turbulent energy spectra. Transition-to-turbulent of the pulsatile flow in the post stenosis is examined through the various numerical results such as velocity, streamlines, velocity vectors, vortices, wall pressure and shear stresses, turbulent kinetic energy, and pressure gradient. A comparison of the LES results with the coarse DNS are given for the Reynolds number of 2000 in terms of the mean pressure, wall shear stress as well as the turbulent characteristics. The results show that the shear stress at the upper wall is low just prior to the centre of the stenosis, while it is maximum in the throat of the stenosis. But, at the immediate post stenotic region, the wall shear stress takes the oscillating form which is quite harmful to the blood cells and vessels. In addition, the pressure drops at the throat of the stenosis where the re-circulated flow region is created due to the adverse pressure gradient. The maximum turbulent kinetic energy is located at the post stenosis with the presence of the inertial sub-range region of slope −5/3.     

Chaotic responses on gear pair system equipped with journal bearings under turbulent flow

This study performs a systematic analysis of the dynamic behavior of a gear-bearing system with the turbulent flow effect, nonlinear suspension, nonlinear oil-film force and nonlinear gear mesh force. The dynamic orbits of the system are observed using bifurcation diagrams plotted using the dimensionless damping coefficient, unbalance coefficient and the dimensionless rotational speed ratio as control parameters. The onset of chaotic motion is identified from the phase diagrams, power spectra, Poincaré maps, Lyapunov exponents and fractal dimension of the gear system. The numerical results reveal that the system exhibits a diverse range of periodic, sub-harmonic and chaotic behaviors. The results presented in this study provide a detailed understanding of the operating conditions under which undesirable dynamic motion takes place in gear-bearing system and therefore offer a useful source of reference for engineers in designing and controlling such systems.     

Study of spectral characteristics of a homogeneous turbulent flow

A spectral representation of kinetic energy for a vortex cascade of instability in a compressible inviscid shear flow is considered, and the Rayleigh-Taylor instability is studied. A comparative analysis is given to the spectral decompositions of kinetic energy for both problems. The classical Kolmogorov −5/3 power law is proved to hold for developed turbulent flows.     

The growth of fractal dimension of an interface evolution from the interaction of a shock wave with a rectangular block of SF6

The interface between air and a rectangular block of sulphur hexafluoride (SF6), impulsively accelerated by the passage of a planar shock wave, undergoes Richtmyer–Meshkov instability and the flow becomes turbulent. The evolution of the interface was previously simulated using a multi-component model based on a thermodynamically consistent and fully conservative formulation and results were validated against available experimental data (Bates et al. Richtmyer–Meshkov instability induced by the interaction of a shock wave with a rectangular block of SF6, Phys Fluids, 2007; 19:036101). In this study, the CFD results are analyzed using the fractal theory approach and the evolution of fractal dimension of the interface during the transition of the flow into fully developed turbulence is measured using the standard box-counting method. It is shown that as the Richtmyer–Meshkov instability on the interface develops and the flow becomes turbulent, the fractal dimension of the interface increases asymptotically toward a value close to 1.39, which agrees well to those measured for classical shear and fully developed turbulences.     

Assessment of turbulence modeling for gas flow in two-dimensional convergent–divergent rocket nozzle

In the present study, the turbulent gas flow dynamics in a two-dimensional convergent–divergent rocket nozzle is numerically predicted and the associated physical phenomena are investigated for various operating conditions. The nozzle is assumed to have impermeable and adiabatic walls with a flow straightener in the upstream side and is connected to a plenum surrounding the nozzle geometry and extended in the downstream direction. In this integrated component model, the inlet flow is assumed a two-dimensional, steady, compressible, turbulent and subsonic. The physics based mathematical model of the considered flow consists of conservation of mass, momentum and energy equations subject to appropriate boundary conditions as defined by the physical problem stated above. The system of the governing equations with turbulent effects is solved numerically using different turbulence models to demonstrate their numerical accuracy in predicting the characteristics of turbulent gas flow in such complex geometry. The performance of the different turbulence models adopted has been assessed by comparing the obtained results of the static wall pressure and the shock position with the available experimental and numerical data. The dimensionless shear stress at the nozzle wall and the separation point are also computed and the flow field is illustrated. The various implemented turbulence models have shown different behavior of the turbulent characteristics. However, the shear-stress transport (SST) k–ω model exhibits the best overall agreement with the experimental measurements. In general, the proposed numerical procedure applied in the present paper shows good capability in predicting the physical phenomena and the flow characteristics encountered in such kinds of complex turbulent flow.     

A joint SSPDF–CMC method for simulating turbulent methane–air jet combustion

A joint single scalar probability density function and conditional moment closure (SSPDF–CMC) method is proposed for modeling a turbulent methane–air jet flame. In general, the probability density function (PDF) of passive scalar (such as mixture fraction) is non-Gaussian and not fully determined by the advecting velocity field, therefore the presumed shape of PDF of mixture fraction assumed as clipped Gaussian distribution or beta function in normal conditional moment closure (CMC) method is incorrect. In SSPDF–CMC method, the PDF of mixture fraction is obtained using a Monte-Carlo method to solve a PDF transport equation. An assumption that the averaged scalar advection is approximately equal to the averaged scalar dissipation in the wake of a grid-generated turbulence flow is adopted to model the averaged scalar dissipation. The predictions using the proposed method are compared with those using the conventional CMC method and the experimental data. It is seen that the predicted Favre conditional averaged statistics and Favre unconditional averaged statistics using the proposed method are in better agreement with the measurement data than those using the conventional CMC method. The predicted conditional or unconditional mean NO even using the SSPDF model is only in fair agreement with the experiments. It shows that the first-order closure for the conditional reaction rate of NO should be improved.     

On the dynamics in a transitional boundary layer

IntroductionIll 1883 Professor Osborne Reynolds published in Philosopl1ical Transactions of the RoyalSociety the outcomes of his flow visua1ization at Manchester. These had shown that whetherthe flow in a pipe was direct to sinuous (or, as nowadays we would say, laminar to turbulent)depended on its Reynolds number. Transition from Iaminar to turbuIent flow becomes animportant probIem i1l fluid mechanics, which has attracted the interest of investigators fOrmore than l00 years. The partic…     

Exergetic aspects of high-speed boundary layers

Exergy is an established concept for evaluating heat cycles. It quantifies the available work that can be extracted from a system. The extension of the classical definition to non-uniform flow field properties is applied here to high-speed boundary layers. The destruction of exergy is quantified by the corresponding loss thickness and is evaluated locally. Turbulent fluctuations in this approach are described by their turbulent exergy, which is a generalisation of the concept of turbulent kinetic energy for incompressible flow. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A stochastic method to account for the ambient turbulence in Lagrangian Vortex computations

This paper describes a detailed implementation of the Synthetic Eddy Method (SEM) initially presented in Jarrin et al. (2006) applied to the Lagrangian Vortex simulation. While the treatment of turbulent diffusion is already extensively covered in scientific literature, this is one of the first attempts to represent ambient turbulence in a fully Lagrangian framework. This implementation is well suited to the integration of PSE (Particle Strength Exchange) or DVM (Diffusion Velocity Method), often used to account for molecular and turbulent diffusion in Lagrangian simulations. The adaptation and implementation of the SEM into a Lagrangian method using the PSE diffusion model is presented, and the turbulent velocity fields produced by this method are then analysed. In this adaptation, SEM turbulent structures are simply advected, without stretching or diffusion of their own, over the flow domain. This implementation proves its ability to produce turbulent velocity fields in accordance with any desired turbulent flow parameters. As the SEM is a purely mathematical and stochastic model, turbulent spectra and turbulent length scales are also investigated. With the addition of variation in the turbulent structures sizes, a satisfying representation of turbulent spectra is recovered, and a linear relation is obtained between the turbulent structures sizes and the Taylor macroscale. Lastly, the model is applied to the computation of a tidal turbine wake for different ambient turbulence levels, demonstrating the ability of this new implementation to emulate experimentally observed tendencies.     

湍流边界层中固体小颗粒湍流运动的Lagrangian模型

给出了固体小颗粒在边界层中的Lagrangian运动方程,方程中包括受壁面影响的粘性阻力,Saffman升力及Magus升力等.使用频谱法,得到了颗粒响应流体的Lagrangian能谱的表达式,使用这些结果研究了各种响应特性.本文的结果清楚地表明了固体个颗粒在湍流扩散过程中,其湍流扩散是可能大于流体的.     

Numerical study of the effect of changing the cyclone cone length on the gas flow field

Models of six different cyclone separator designs have been investigated numerically to study the effect of changing cone length on the gas flow field. The cone length was changed in two ways, one by fixing the cone diameter and the other by fixing the cone angle. The Reynolds stress model was applied to predict the turbulence of gas flow. The model was validated by comparing the numerical results to published experimental data. The results show that, for the cyclone without cylindrical part, the pressure drop and the value of tangential and axial velocity components increase with increasing cone length with fixed cone angle until they reach maximum value. It was found that increasing the cone length with the same cone diameter increases the back flow and the lip leakage below the vortex finder. The back flow and the lip leakage phenomena was observed to be reduced by increasing cone length with fixed cone angle. The effect of contraction ratio D_x/D, which is defined as the ratio between the vortex finder diameter and the cone diameter, on the turbulent kinetic energy and the velocity components also was investigated. As the contraction ratio decreases, the maximum value of the turbulent kinetic energy and the tangential and axial velocity components increases. The pressure drop was fitted into a correlation as a function of dimensionless quantities (Re, D_x/D, and l_c/D) for each design, so that the designers could select the appropriate dimensions according to their design limitations.     

Theoretical Considerations about Near-Wall Turbulence and Resulting Flow Control Schemes

Based on analytical considerations about how near-wall turbulence needs to be modified in order to reduce the momentum loss towards solid walls and to yield lower energy losses, a mechanism of turbulent drag reduction is proposed. This mechanism suggests that drag reducing flow control at high Reynolds numbers should be designed to minimize the turbulent dissipation rate. A previously published approach on how a reduction of the turbulent dissipation in the near-wall region can be achieved is analyzed further. The obtained results provide some new insight on the parameters that need to be considered when designing flow control schemes for skin friction drag reduction. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Asymptotic behaviour of commutation errors and the divergence of the Reynolds stress tensor near the wall in the turbulent channel flow

The derivation of the space averaged Navier–Stokes equations for the large eddy simulation (LES) of turbulent incompressible flows introduces two groups of terms which do not depend only on the space averaged flow field variables: the divergence of the Reynolds stress tensor and commutation errors. Whereas the former is studied intensively in the literature, the latter terms are usually neglected. This note studies the asymptotic behaviour of these terms for the turbulent channel flow at a wall in the case that the commutation errors arise from the application of a non‐uniform box filter. To perform analytical calculations, the unknown flow field is modelled by a wall law (Reichardt law and 1/αth power law) for the mean velocity profile and highly oscillating functions model the turbulent fluctuations. The asymptotics show that near the wall, the commutation errors are at least as important as the divergence of the Reynolds stress tensor. Copyright © 2006 John Wiley & Sons, Ltd.     

Heat Transfer in Rotor/Rotor and Rotor/Stator Cavity: Physics,Correlations and Numerical Modeling

The transitional and turbulent flow in the near wall sublayer is now mostly modeled based on the existing knowledge of simple 2D flows. To determine the effect of three dimensionality on the turbulent flow structures and turbulent heat transfer in the near wall areas the authors investigate numerically (SVV) turbulent flow in rotor/stator and rotor/rotor flows (with and without axial throughflow). These simple model flows contain most of the phenomena that are needed to understand more complex, 3D transitional and turbulent flows. Attention is focused on the turbulent characteristics which should have more universal character. To stabilize calculations for high Reynolds numbers (up to Re=800 000) the SVV operator is introduced into the Navier-Stokes and energy equations solver for cylindrical coordinate system without using complex numbers. Code optimization and parallelization have speeded up computations 20 times. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

旋转流动的低模分析及仿真研究

为了探讨Couette-Taylor流从层流到湍流过渡的方式以及流动发展到湍流之后混沌吸引子的某些特征等问题,采用低模分析方法研究了Couette-Taylor流的部分动力学行为及仿真问题,讨论了Couette-Taylor流三模态类Lorenz型方程组的动力学行为,包括定态的失稳、极限环的出现、分岔与混沌的演变和全局稳定性分析等。通过线性稳定性分析和数值模拟等方法给出了此三维模型分岔与混沌等动力学行为及其演化历程,并借此解释了Couette-Taylor流试验中观察到的部分涡流的演化过程.基于系统的分岔图、Lyapunov指数谱、功率谱、Poincaré(庞加莱)截面和返回映射等揭示了系统混沌行为的普适特征.     

A kinetic energy preserving DG scheme applied to mechanical FSI

In the context of mechanical fluid-structure interaction (FSI) comprising moving or deforming structures, fluid discretizations need to cope with time-dependent fluid domains and resulting grid deformations in addition to the general challenges regarding e.g. boundary layers and turbulent phenomena. Recent approaches in the simulation of compressible turbulent flow are based on so-called split forms of conservation laws to guarantee the preservation of secondary physical quantities such as kinetic energy. For the simulation of turbulent flows, this often leads to a better representation of the kinetic energy spectrum. Initially, kinetic energy preserving(KEP) DG schemes have been constructed on Gauss-Legendre-Lobatto(GLL) nodes containing the interval end points, however, KEP DG schemes based on the classical Gauss-Legendre(GL) nodes are potentially more accurate and may be also more efficient than its GLL variant for certain applications. In this work, the KEP-DG schemes both on GL and GLL nodes are applied to the classical moving piston test case via an ALE formulation on moving fluid grids showing a more accurate frequency representation of the structure displacement in case of GLL nodes. (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Conservation laws of turbulence models

The conservation of mass, momentum, energy, helicity, and enstrophy in fluid flow are important because these quantities organize a flow, and characterize change in the flow's structure over time. In turbulent flow, conservation laws remain important in the inertial range of wave numbers, where viscous effects are negligible. It is in the inertial range where energy, helicity (3d), and enstrophy (2d) must be accurately cascaded for a turbulence model to be qualitatively correct. A first and necessary step for an accurate cascade is conservation; however, many turbulent flow simulations are based on turbulence models whose conservation properties are little explored and might be very different from those of the Navier-Stokes equations.We explore conservation laws and approximate conservation laws satisfied by LES turbulence models. For the Leray, Leray deconvolution, Bardina, and Nth order deconvolution models, we give exact or approximate laws for a model mass, momentum, energy, enstrophy and helicity. The possibility of cascades for model quantities is also discussed.     

基于有限分析法的平面铅垂壁羽流的相似解

应用计及浮力对湍动能及其耗散率的影响的k-epsilon湍流模式,结合有限分析法对密度差引起的平面铅垂紊动壁羽流进行了数值分析．在均匀环境条件下控制该类流动的连续性方程、流动方向的动量方程、浓度扩散方程、湍动能及其耗散率方程存在相似解．考虑到浮力通量守恒条件,应用有限分析法给出了铅直壁羽流的速度、相对密度差、湍动能及耗散率的分布,进而给出了各物理量最大值沿主流方向变化的关系式．湍流Schmidt数为1.0时的计算结果与实验资料吻合较好,表明应用有限分析法分析铅垂平面羽流是有效的,即在分析壁羽流时浮力对湍动能及其耗散率的影响应该予以考虑．