Direct Numerical Simulation of Turbulence in a Stably Stratified Fluid and Wave-Shear Interaction

Turbulence decay in a strongly stratified medium is simulated by a direct pseudo-spectral code solving the three-dimensional equations of motion under the Boussinesq approximation. The results are compared to non-stratified simulations results. We focus on the production of mean shear energy observed in the stratified case. We then simulate the decay of stratified turbulence when affected by an initial horizontal mean flow and show that this mean flow is the major component remaining at large t. Next, we give some analytical elements on wave-shear interaction by using a simple refraction calculation with WKB hypothesis. This calculation is illustrated by simulating the interaction between one monochromatic internal wave and a vertical shear profile. We conclude that the existence of singularities in the mean shear production term in the presence of internal gravity waves may be one of the possible mechanisms involved within stratified turbulent shear flows.     

Direct Numerical Simulation of a Recirculating, Swirling Flow

A direct numerical simulation (DNS) of a recirculating, swirling flow is performed at a Reynolds number of 5000. Detailed one and two point statistics are presented in this paper. Flow visualization and frequency analysis are used to identify a precessing vortex core and to characterize its position, extent and influence on the flow field. The results are compared with laser Doppler velocimetry (LDV) measurements as well as large eddy simulation (LES) data reported in the literature. The present work constitutes a first step in setting up a DNS data base for complex flows.     

Direct Numerical Simulation of the Three-Dimensional Transition to Turbulence in the Transonic Flow around a Wing

The three-dimensional transition to turbulence in the transonic flow around a NACA0012 wing of constant spanwise section has been analysed at zero incidence and within the Reynolds number range [3000, 10000], by performing the direct numerical simulation. The successive stages of the 2D and 3D transition beyond the first bifurcation have been identified. A 2D study has been carried out near the threshold concerning the appearance of the first bifurcation that is the von-Kármán instability. The critical Mach number associated with this flow transition has been evaluated. Three other successive stages have been detected as the Mach number further increases in the range [0.3, 0.99]. Concerning the 3D transition of this nominally 2D flow configuration, the amplification of the secondary instability has been studied within the Reynolds number range of [3000, 5000]. The formation of counter-rotating longitudinal vorticity cells and the consequent appearance of a large-scale spanwise wavelength have been obtained downstream of the trailing-edge shock. A vortex dislocation pattern is developed as a consequence of the shock-vortex interaction near the trailing edge. The subcritical nature of the present 3D transition to turbulence has been proven by means of the DNS amplification signals and the Landau global oscillator model.     

Direct Numerical Simulation of the Three-Dimensional Transition to Turbulence in the Incompressible Flow around a Wing

The onset of the secondary instability and the successive steps ofthe 3D transition to turbulence are examined in the flow around a wing of NACA0012 section, constant along the spanwise direction. The wing is placed in a uniform flow upstream, at 20° of incidence and at the Reynolds number of 800. The spanwise length is equal to four chords. The objectives of this study are the identification of the three-dimensional transition mechanism and the development of the early stages of turbulence in the present class of unsteady aerodynamic flows. A detailed processing of the DNS signals carried out by an appropriate conditional sampling allows the identification of the physical mechanisms related to the birth of turbulence and to the non-linear interaction with the 3D coherent structures in the near region.     

Reacting, Circular Mixing Layers in Transition to Turbulence

The evolution of a reacting, circular mixing layer - a model of round-jet flow - in its transition to turbulence was studied by direct numerical simulation. An economical Fourier pseudospectral method was combined with the third-order Adams-Bashforth scheme to integrate Navier-Stokes and scalar transport equations. The Reynolds number based on initial mixing-layer diameter and velocity difference was 1600. The initially thin mixing layer encloses a cylindrical core of fuel that mixes and reacts with the surrounding oxidizer. Both fast and finite-rate reactions were examined. The stages in transition are characterized by roll-up of the mixing layer into a sequence of vortex rings, pairing of adjacent rings, azimuthal instability, and breakdown to a disordered (turbulent) state. Reaction surfaces in the fast reaction limit become extended, folded and pinched off at various times corresponding to the dynamics of the vortices observed in the simulations. When the equivalence ratio is O(1) or smaller,the progress of reaction is determined by the dynamics of vortex rings. For larger ratios there is a qualitative difference: Initially, the flame is located well outside the rings and is relatively unaffected. Following breakdown to turbulence, there is a steep increase in flame surface area resulting in a noticeable change in fuel consumption rate. At smaller reaction rates (small Damkohler numbers), the reaction zones are diffuse and fill the vortical (mixed) regions. Product accumulates in and its presence raises the temperature of vortex cores, but reaction rates remain low due to low reactant concentrations. Reaction rates are highest in the braids between vortex rings where scalar dissipation rates and compressive strain rates show the highest values. This revised version was published online in July 2006 with corrections to the Cover Date.     

弯道水流数值模拟研究

采用边界拟合坐标技术,生成边界处正交的曲线网格,在曲线坐标系中进行弯道水流数值计算。借鉴现有关于弯道水流流速分布的研究成果,对曲线坐标系平面二维浅水方程作了修正,修正后的方程组能够计入弯道环流引起的横向的动量交换,即考虑了二次流对流线弯曲的复杂水力特性的影响。通过连续弯道实例计算,验证了该模型的可靠性和实用性。     

Direct Numerical Simulation of Self-Similar Turbulent Boundary Layers in Adverse Pressure Gradients

Direct numerical simulations of the Navier–Stokes equations have been carried out with the objective of studying turbulent boundary layers in adverse pressure gradients. The boundary layer flows concerned are of the equilibrium type which makes the analysis simpler and the results can be compared with earlier experiments and simulations. This type of turbulent boundary layers also permits an analysis of the equation of motion to predict separation. The linear analysis based on the assumption of asymptotically high Reynolds number gives results that are not applicable to finite Reynolds number flows. A different non-linear approach is presented to obtain a useful relation between the freestream variation and other mean flow parameters. Comparison of turbulent statistics from the zero pressure gradient case and two adverse pressure gradient cases shows the development of an outer peak in the turbulent energy in agreement with experiment. The turbulent flows have also been investigated using a differential Reynolds stress model. Profiles for velocity and turbulence quantities obtained from the direct numerical simulations were used as initial data. The initial transients in the model predictions vanished rapidly. The model predictions are compared with the direct simulations and low Reynolds number effects are investigated.     

Numerical analysis of a 3D hydrodynamic contact

We study here the numerical analysis of a hydrodynamic contact in a particular configuration: the 3D incompressible viscous flow of a fluid dragged by a smooth plate over a rough surface. The mathematical model takes into account and discretizes the local topography of the rough profile. The simulation outcome will be the 3D velocity and pressure fields of the fluid film within the contact borders. This work is limited to the study of numerical resolution methods working solely in finite differences. The algorithms will be tested by analysing and comparing their results with analytically known flows. Copyright © 2006 John Wiley & Sons, Ltd.     

用混合网格及各向异性多重网格法求解三维可压紊流流动

采用混合网格求解紊流Navier Stokes方程。在物面附近采用柱状网格 ,其他区域则采用完全非结构网格。方程的求解采用Jamson的有限体积法 ,紊流模型采用两层Baldwin Lomax代数紊流模型。用各向异性多重网格法来加速解的收敛。数值算例表明 ,用混合网格及各向异性多重网格求解紊流流动是非常有效的     

液固两相槽道流中湍流调制的数值研究

结合雷诺应力模型 (Reynolds Stress Model, RSM) 和混合模型 (Mixture Model) 对槽道湍流向下流动中的颗粒调制湍流问题进行了研究．该模型考虑了颗粒流的动能理论和颗粒对湍流的反馈作用．着重分析了颗粒对湍流的调制作用,以及颗粒参数变量（如颗粒密度和质量载荷）对湍流调制的影响．结果表明：(1)在颗粒抑制湍流的范围内,当颗粒密度小于载流体密度时,湍流强度的改变量与颗粒密度成反比;当颗粒密度大于载流体密度时,湍流强度的改变量与颗粒密度成正比;(2)在一定范围内,颗粒抑制湍流的能力随颗粒质量载荷增加而变强．     

各向同性湍流内颗粒碰撞率的直接模拟研究

对 Re_{\lambda } 约为51均匀各向同性湍流内 St_{k}(=\tau_{p}/\tau_{k}) 为 0 ~10.0 的 有限惯性颗粒的碰撞行为进行了直接数值模拟，以研究湍流对有限惯性 颗粒碰撞的影响. 结果表明，具有一定惯性颗粒的湍流碰撞率完全不同于零惯性的轻颗粒 (St_{k}=0) 和可忽略湍流作用的重颗粒 (St{k} \to \infty) , 其变化趋势极其复杂： 在Stk为 0~1.0 之间，颗粒的碰撞率随 St 的增加而近乎线性地剧烈增长，在 Stk≈1.0 3.0(对应的StE=τp/Te≈0.5)附近，颗粒碰撞率出现两个峰值，在Stk>3.0以后，颗粒的碰撞率随惯性增 大而逐渐趋向于重颗粒极限；在峰值处，有限惯性颗粒的平均碰撞率的峰值较轻颗粒增强了 30倍左右. 为进一步分析湍流作用下颗粒碰撞率的影响因素，分别使用可能发生碰撞 的颗粒对的径向分布函数和径向相对速度来量化颗粒的局部富集效应和湍流掺混效应，表明 St_{k} \approx 1.0 时局部富集效应最为强烈，使得颗粒的碰撞率出现第1个峰值； 湍流掺混效应则随着颗粒Stk的增大而渐近增大；局部富集和湍流掺混联合作用的结果， 使得颗粒碰撞率在 St_{k} \approx 3.0 附近出现另一个峰值.     

Spatial Direct Numerical Simulation of the Large Vortical Structures in Forced Plumes

Direct numerical simulation (DNS) of forced plumes arising frominput of both momentum and buoyancy into an ambient fluid is presented.The large vortical structures in the near field of thermal and reactiveplumes are investigated. Boundary conditions associated with the spatialDNS of open-boundary buoyant flows that are compatible with the modernnon-dissipative, high-order, finite-difference schemes have beendeveloped. The governing equations for flow and combustion at the plumecenterline are put into a special form to circumvent the singularity atthe axis associated with the cylindrical coordinates. Mixing is found tobe stronger in the planar thermal plume than in the axisymmetric case.An explanation is provided based on the vorticity budget. Axisymmetricreactive plumes with a one-step reaction governed by the Arrheniuskinetics have also been studied. The unsteady effects of chemical heatrelease and combustion-induced buoyancy on the flow structures areinvestigated. Budgets of the vorticity transport are examined to revealthe mechanisms leading to the formation and evolution of large vorticalstructures in forced plumes. It is found that volumetric expansion dueto chemical heat release tends to destroy vorticity, whilecombustion-induced buoyancy under the gravitational effect generatesvorticity. The gravitational term in the vorticity transport equation isfound to be the main mechanism for the buoyant flow instability and thedevelopment of counter-rotating vortices in reactive plumes.     

采用NND方法计算三维喷管气流场

本文运用NND显式差分格式,计算了三维喷管气流场。气流场计算的基本方程为一般贴体坐标系下三维守恒型的欧拉方程。采用了时间分裂法和Steger-Warming矢通量分裂技术。在喷管内沿周向的每个由轴线和壁面构成的子午面上根据泊松方程生成贴体网格。本文运用三维程序计算了轴对称JPL喷管,同时与实验结果和前人采用轴对称二维程序所计算的结果做了对比。最后,本文还计算了三维矢量喷管,计算结果与现有的实验结果一致。通过轴对称JPL喷管和三维矢量喷管的计算考核,表明建立的算法和编写的计算程序是正确的。文中提出了采用子午面形式的贴体网格时奇性轴的处理方法。计算结果表明在喷管壁面处,马赫数与压强的计算结果与实验值吻合较好,而在喷管轴线处,只有当网格较密时,才能得出与实验结果接近的计算结果。     

Large-Eddy Simulation of Transition to Turbulence in Boundary Layers

Large-eddy simulation (LES) results for laminar-to-turbulent transition in a spatially developing boundary layer are presented. The disturbances are ingested into a laminar flow through an unsteady suction-and-blowing strip. The filtered, three-dimensional time-dependent Navier–Stokes equations are integrated numerically using spectral, high-order finite-differences, and a three-stage low-storage Runge–Kutta/Crank–Nicolson time-advancement method. The buffer-domain technique is used for the outflow boundary condition. The localized dynamic model used to parametrize the subgrid-scale (SGS) stresses begins to have a significant impact at the beginning of the nonlinear transition (or intermittency) region. The flow structures commonly found in experiments are also observed in the present simulation; the computed linear instability modes and secondary instability $\Lambda$-vortex structures are in agreement with the experiments, and the streak-like structures and turbulent statistics compare with both the experiments and the theory. The physics captured in the present LES are consistent with the experiments and the full Navier–Stokes simulation (DNS), at a significant fraction of the DNS cost. A comparison of the results obtained with several SGS models shows that the localized model gives accurate results both in a statistical sense and in terms of predicting the dynamics of the energy-carrying eddies, while requiring fewer ad hoc adjustments than the other models. Received: 5 April 1996 and accepted 27 March     

A Robust High-Resolution Split-Type Compact FD Scheme for Spatial Direct Numerical Simulation of Boundary-Layer Transition

In this paper an efficient split-type Finite-Difference (FD) scheme with high modal resolution – most important for the streamwise convection terms that cause wave transport and interaction – is derived for a mixed Fourier-spectral/FD method that is designed for the spatial direct numerical simulation (DNS) of boundary-layer transition and turbulence. Using a relatively simple but thorough and instructive modal analysis we discuss some principal trouble sources of the related FD discretization. The new scheme is based on a 6th-order compact FD discretization in streamwise and wall-normal direction and the classical 4th-order Runge–Kutta time-integration scheme with symmetrical final corrector step. Exemplary results of a fundamental-(K-) type breakdown simulation of a strongly decelerated Falkner–Skan boundary layer (Hartree parameter _H = – 0.18) using 70 mega grid points in space are presented up to the early turbulent regime (Re_,turb 820). The adverse pressure gradient gives rise to local separation zones during the breakdown stage and intensifies final breakdown by strong amplification of (background) disturbances thus enabling rapid transition at moderate Reynolds number. The appearance and dynamics of small-scale vortical structures in early turbulence basically similar to the large-scale structures at transition can be observed corroborating Kachanov's hypothesis on the importance of the K-regime of breakdown for coherent structures in turbulence.     

二层密度分层流体重力流的数值模拟

从RANS方程和RNGk-ε湍流模型出发,采用流体体积法(VOF)来模拟密度分层流动,对盐水和淡水因密度差异导致的分层重力流动现象进行了数值模拟.文中报道了平底水槽重力流、狭孔交换流的数值模拟结果,分层重力流锋面运动速度的计算值与现有的半理论半经验公式一致.为了揭示地形变化对分层重力流的影响,对设有缓变潜堤的水槽内分层重力流动的形成过程进行了数值模拟,给出了重力头推进速度和局部流场的计算结果,并讨论了分层流界面、流量和锋面附近的流速分布特征.     

机动弹头高超声速粘性绕流数值模拟

采用空间二阶精度的交替方向隐式分解的NND格式求解完全气体假定下的非定常薄层近似NaVierStokes方程,并采用抛物化的椭圆型方程生成复杂带翼弹头的空间网格.最后给出了机动弹头在M∞=7.3,α=20°下的流场计算结果.计算结果表明本文所采用的计算网格及计算方法可以适用于高超声速复杂流场的数值模拟.     

Numerical simulation of 3D free surface flows,with multiple incompressible immiscible phases. Applications to impulse waves

A numerical method for the solution to the density‐dependent incompressible Navier–Stokes equations modeling the flow of N immiscible incompressible liquid phases with a free surface is proposed. It allows to model the flow of an arbitrary number of liquid phases together with an additional vacuum phase separated with a free surface. It is based on a volume‐of‐fluid approach involving N indicator functions (one per phase, identified by its density) that guarantees mass conservation within each phase. An additional indicator function for the whole liquid domain allows to treat boundary conditions at the interface between the liquid domain and a vacuum. The system of partial differential equations is solved by implicit operator splitting at each time step: first, transport equations are solved by a forward characteristics method on a fine Cartesian grid to predict the new location of each liquid phase; second, a generalized Stokes problem with a density‐dependent viscosity is solved with a FEM on a coarser mesh of the liquid domain. A novel algorithm ensuring the maximum principle and limiting the numerical diffusion for the transport of the N phases is validated on benchmark flows. Then, we focus on a novel application and compare the numerical and physical simulations of impulse waves, that is, waves generated at the free surface of a water basin initially at rest after the impact of a denser phase. A particularly useful application in hydraulic engineering is to predict the effects of a landslide‐generated impulse wave in a reservoir. Copyright © 2014 John Wiley & Sons, Ltd.     

基于流动仿真大数据应对旋涡-湍流的研究进展

文章以流体科学进入二十一世纪后,在大规模超算、云存储、数据通信和人工智能为支撑的大数据时代背景下,结合目前在复旦大学航空航天系所构建的热流体湍流直接数值仿真数据库,以及复旦大学团队近期与美国德州大学刘超群教授、上海理工大学蔡小舒教授以及国内水动力学杂志编辑部所合作开展的第三代涡识别技术研究,初步概念性地展示旋涡和湍流,特别是针对有工程实际背景和直接应用价值的壁湍流,在这两个流体力学关键基础议题上的最新认知,和基于大数据深度学习的相关湍流工程模拟实践成果．这些成果包括：(1) 基于第三代涡识别技术的尾迹湍流中的涡运动学和动力学探索;(2) 流-热统一完整的类-1、类-2湍流边界层壁面律构建;(3) 基于第三代涡识别量对Kolmogorov 幂次律的再认知;(4) 基于DNS统计数据和神经网络深度学习构建新型湍流封闭模型及RANS计算实践．通过这些成果展示,论证解决这两个基础流体科学议题的技术路径,进而促进流体及相关学科研究在现代大数据背景下取得实质性进展和突破,并惠及现代流体、气动、水利、动力和化工等工程领域．