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.     

基于基因表达式编程的数据驱动湍流建模

计算流体动力学是湍流研究的重要手段, 其中雷诺平均模拟在航空航天等实际工程中得到了广泛应用. 雷诺平均模拟的结果很大程度上依赖于湍流模型的预测精度, 而实际工程应用中常用的模型往往精度有限. 近年来, 数据驱动的湍流建模方法得到越来越多的关注. 本文介绍了基于基因表达式编程 (gene-expression programming, GEP) 方法的湍流建模相关进展. 本文首先讨论基因表达式编程应用于湍流建模的具体方法, 包括基本算法、显式代数应力模型和湍流传热两种建模框架、模型测试方法以及损失函数设置等. 在此基础上, 基因表达式编程方法被应用于涡轮叶栅尾流混合、竖直平板间自然对流、三维横向流中的射流等问题. 结果表明, GEP可以有效提升常用模型对于尾流混合损失、壁面热通量等关键参数的预测精度. 基因表达式编程方法可以显式给出模型方程, 因此模型具有可解释性强等特点. 基于双向耦合方法得到的模型还被证明具有较好的后验测试精度和鲁棒性. 基因表达式编程方法还被初步应用于大涡模拟亚格子应力和边界层转捩等问题的建模, 在不同湍流建模领域表现出很大的潜力.      

KOLMOGOROV BEHAVIOR OF NEAR-WALL TURBULENCE AND ITS APPLICATION IN TURBULENCE MODELING

The near-wall behavior of turbulence is re-examined in a way different from that proposed by Hanjalic and Launder¹ and followers^2,3,4,5. It is shown that at a certain distance from the wall, all energetic large eddies will reduce to Kolmogorov eddies (the smallest eddies in turbulence). All the important wall parameters, such as friction velocity, viscous length scale, and mean strain rate at the wall, are characterised by Kolmogorov microscales. According t o this Kolmogorov behavior of near-wall turbulence, the turbulence quantities, such as turbulent kinetic energy, dissipation rate, etc. at the location where the large eddies become “Kolmogorov” eddies, can be estimated by using both direct numerical simulation (DNS) data and asymptotic analysis of near-wall turbulence. This information will provide useful boundary conditions for the turbulent transport equations. As a n example, the concept is incorporated in the standard κ - εmodel which is then applied t o channel and boundary layer flows. Using appropriate boundary conditions (based on Kolmogorov behaviour of near-wall turbulence), there is no need for any wall-modification to the κ - ε equations (including model constants). Results compare very well with the DNS and experimental data.     

Toward low‐noise synthetic turbulent inflow conditions for aeroacoustic calculations

The accuracy of boundary conditions for computational aeroacoustics is a well‐known challenge, due in part to the necessity of truncating the flow domain and replacing the analytical boundary conditions at infinity with numerical boundary conditions. In particular, the inflow boundary condition involving turbulent velocity or scalar fields is likely to introduce spurious waves into the domain, therefore degrading the flow behavior and deteriorating the physical acoustic waves. In this work, a method to generate low‐noise, divergence‐free, synthetic turbulence for inflow boundary conditions is proposed. It relies on the classical view of turbulence as a superposition of random eddies convected with the mean flow. Within the proposed model, the vector potential and the requirement that the individual eddies must satisfy the linearized momentum equations about the mean flow are used. The model is tested using isolated eddies convected through the inflow boundary and an experimental benchmark data for spatially decaying isotropic turbulence. Copyright © 2013 John Wiley & Sons, Ltd.     

Incompressible Roe‐like scheme for turbulent flows

In the current study, numerical investigation of incompressible turbulent flow is presented. By the artificial compressibility method, momentum and continuity equations are coupled. Considering Reynolds averaged Navier–Stokes equations, the Spalart–Allmaras turbulence model, which has accurate results in two‐dimensional problems, is used to calculate Reynolds stresses. For convective fluxes a Roe‐like scheme is proposed for the steady Reynolds averaged Navier–Stokes equations. Also, Jameson averaging method was implemented. In comparison, the proposed characteristics‐based upwind incompressible turbulent Roe‐like scheme, demonstrated very accurate results, high stability, and fast convergence. The fifth‐order Runge–Kutta scheme is used for time discretization. The local time stepping and implicit residual smoothing were applied as the convergence acceleration techniques. Suitable boundary conditions have been implemented considering flow behavior. The problem has been studied at high Reynolds numbers for cross flow around the horizontal circular cylinder and NACA0012 hydrofoil. Results were compared with those of others and a good agreement has been observed. Copyright © 2012 John Wiley & Sons, Ltd.     

Reynolds number effects on a turbulent boundary layer with separation,reattachment, and recovery

The present paper addresses experimental studies of Reynolds number effects on a turbulent boundary layer with separation, reattachment, and recovery. A momentum thickness Reynolds number varies from 1,100 to 20,100 with a wind tunnel enclosed in a pressure vessel by varying the air density and wind tunnel speed. A custom-built, high-resolution laser Doppler anemometer provides fully resolved turbulence measurements over the full Reynolds number range. The experiments show that the mean flow is at most a very weak function of Reynolds number while turbulence quantities strongly depend on Reynolds number. Roller vortices are generated in the separated shear layer caused by the Kelvin–Helmholtz instability. Empirical Reynolds number scalings for the mean velocity and Reynolds stresses are proposed for the upstream boundary layer, the separated region, and the recovery region. The inflectional instability plays a critical role in the scaling in the separated region. The near-wall flow recovers quickly downstream of reattachment even if the outer layer is far from an equilibrium state. As a result, a stress equilibrium layer where a flat-plate boundary layer scaling is valid develops in the recovery region and grows outward moving downstream.     

On the construction of manufactured solutions for one and two‐equation eddy‐viscosity models

This paper presents manufactured solutions (MSs) for some well‐known eddy‐viscosity turbulence models, viz. the Spalart & Allmaras one‐equation model and the TNT and BSL versions of the two‐equation k–ω model. The manufactured flow solutions apply to two‐dimensional, steady, wall‐bounded, incompressible, turbulent flows. The two velocity components and the pressure are identical for all MSs, but various alternatives are considered for specifying the eddy‐viscosity and other turbulence quantities in the turbulence models. The results obtained for the proposed MSs with a second‐order accurate numerical method show that the MSs for turbulence quantities must be constructed carefully to avoid instabilities in the numerical solutions. This behaviour is model dependent: the performance of the Spalart & Allmaras and k–ω models is significantly affected by the type of MS. In one of the MSs tested, even the two versions of the k–ω model exhibit significant differences in the convergence properties. Copyright © 2006 John Wiley & Sons, Ltd.     

Mono‐block and non‐matching multi‐block structured mesh adaptation based on aerodynamic functional total derivatives for RANS flow

An enhanced goal‐oriented mesh adaptation method is presented based on aerodynamic functional total derivatives with respect to mesh nodes in a Reynolds‐Averaged Navier‐Stokes (RANS) finite‐volume mono‐block and non‐matching multi‐block‐structured grid framework. This method falls under the category of methods involving the adjoint vector of the function of interest. The contribution of a Spalart–Allmaras turbulence model is taken into account through its linearization. Meshes are adapted accordingly to the proposed indicator. Applications to 2D RANS flow about a RAE2822 airfoil in transonic, and detached subsonic conditions are presented for the drag coefficient estimation. The asset of the proposed method is patent. The obtained 2D anisotropic mono‐block mesh well captures flow features as well as global aerodynamic functionals. Interestingly, the constraints imposed by structured grids may be relaxed by the use of non‐matching multi‐block approach that limits the outward propagation of local mesh refinement through all of the computational domain. The proposed method also leads to accurate results for these multi‐block meshes but at a fraction of the cost. Finally, the method is also successfully applied to a more complex geometry, namely, a mono‐block mesh in a 3D RANS transonic flow about an M6 wing. Copyright © 2016 John Wiley & Sons, Ltd.     

Implicit application of non‐reflective boundary conditions for Navier–Stokes equations in generalized coordinates

The non‐reflective boundary conditions (NRBC) for Navier–Stokes equations originally suggested by Poinsot and Lele (J. Comput. Phys. 1992; 101 :104–129) in Cartesian coordinates are extended to generalized coordinates. The characteristic form Navier–Stokes equations in conservative variables are given. In this characteristic‐based method, the NRBC is implicitly coupled with the Navier–Stokes flow solver and are solved simultaneously with the flow solver. The calculations are conducted for a subsonic vortex propagating flow and the steady and unsteady transonic inlet‐diffuser flows. The results indicate that the present method is accurate and robust, and the NRBC are essential for unsteady flow calculations. Copyright © 2005 John Wiley & Sons, Ltd.     

A novel hybrid finite element with a hole for analysis of plane piezoelectric medium with defects

A method is presented to formulate a novel hybrid finite element to obtain accurate distributions of mechanical and electrical quantities around a hole in plane piezoelectric mediums. The complex variable method is used in conjunction with Reissner's variational principle to formulate a hybrid special element with an elliptical hole. Detailed derivations are given and numerical examples are performed to demonstrate the accuracy and efficiency of the novel special element. Accurate results around the hole boundary are obtained for infinite and finite piezoelectric medium by using proposed special element.     

Adaptive wall functions for the v2‐f turbulence model

Adaptive wall functions for the v²‐f turbulence model have been derived for the flow over a flat plate at zero pressure gradient. These wall functions were implemented via tables for the turbulence quantities and the friction velocity u_τ. A special treatment for the ε and f boundary conditions is proposed. On fine grids (y⁺<1) this approach yields results consistent with the wall integration solution. Detailed numerical results are presented for a zero pressure gradient boundary layer and separated flow over a ramp. Copyright © 2006 John Wiley & Sons, Ltd.     

Simulation of the Effect of the Freestream Turbulence Parameters on Heat Transfer in an Unsteady Boundary Layer

A variant of the two-parameter turbulence model which makes it possible continuously to calculate a flow region with laminar, transition and turbulent regimes is proposed for investigating the flow under conditions of high freestream turbulence intensity. It is shown that the properties of the thermal transition can be theoretically described using the quasi-steady turbulence model in the case of periodic freestream velocity distribution. The numerical results are compared with theoretical and experimental data. The approach proposed is developed for determining the combined effect of the parameters of harmonic fluctuations of the external velocity and freestream turbulence on the heat transfer characteristics on a flat plate with different boundary conditions for the enthalpy.     

SPH computation of plunging waves using a 2‐D sub‐particle scale (SPS) turbulence model

The paper presents a 2‐D large eddy simulation (LES) modelling approach to investigate the properties of the plunging waves. The numerical model is based on the smoothed particle hydrodynamics (SPH) method. SPH is a mesh‐free Lagrangian particle approach which is capable of tracking the free surfaces of large deformation in an easy and accurate way. The Smagorinsky model is used as the turbulence model due to its simplicity and effectiveness. The proposed 2‐D SPH–LES model is applied to a cnoidal wave breaking and plunging over a mild slope. The computations are in good agreement with the documented data. Especially the computed turbulence quantities under the breaking waves agree better with the experiments as compared with the numerical results obtained by using the k–ε model. The sensitivity analyses of the SPH–LES computations indicate that both the turbulence model and the spatial resolution play an important role in the model predictions and the contributions from the sub‐particle scale (SPS) turbulence decrease with the particle size refinement. Copyright © 2006 John Wiley & Sons, Ltd.     

SPH modelling of depth‐limited turbulent open channel flows over rough boundaries

A numerical model based on the smoothed particle hydrodynamics method is developed to simulate depth‐limited turbulent open channel flows over hydraulically rough beds. The 2D Lagrangian form of the Navier–Stokes equations is solved, in which a drag‐based formulation is used based on an effective roughness zone near the bed to account for the roughness effect of bed spheres and an improved sub‐particle‐scale model is applied to account for the effect of turbulence. The sub‐particle‐scale model is constructed based on the mixing‐length assumption rather than the standard Smagorinsky approach to compute the eddy‐viscosity. A robust in/out‐flow boundary technique is also proposed to achieve stable uniform flow conditions at the inlet and outlet boundaries where the flow characteristics are unknown. The model is applied to simulate uniform open channel flows over a rough bed composed of regular spheres and validated by experimental velocity data. To investigate the influence of the bed roughness on different flow conditions, data from 12 experimental tests with different bed slopes and uniform water depths are simulated, and a good agreement has been observed between the model and experimental results of the streamwise velocity and turbulent shear stress. This shows that both the roughness effect and flow turbulence should be addressed in order to simulate the correct mechanisms of turbulent flow over a rough bed boundary and that the presented smoothed particle hydrodynamics model accomplishes this successfully. © 2016 The Authors International Journal for Numerical Methods in Fluids Published by John Wiley & Sons Ltd     

基于雷诺应力模型的平衡大气边界层模拟

基于雷诺应力湍流模型（简称RSM模型）,研究了平衡大气边界层风场数值模拟问题．假设流体不可压,且不计雷诺应力输运方程中的对流项、浮力产生项、系统旋转产生项和扩散项,在准各向同性的条件下,推导出RSM模型湍动能k的表达式是标准k-ε模型k常数表达式的0.893倍．考虑k沿高度变化的修正,根据在标准k-ε模型中满足水平均匀性的湍流来流边界条件,提出在RSM模型中产生平衡大气边界层的湍流来流边界条件．基于空风洞的数值模拟结果表明,与工程上常用的湍流来流边界条件相比,基于本文提出的湍流来流边界条件得到的风场水平均匀性更优,且在整个流域内,得到的雷诺应力剖面更合适．从而验证了该湍流来流边界条件的适用性．     

海岸城市上空热力内边界层的基本特征

海岸热力内边界层内大气在热力驱动下产生垂向对流运动,对海岸带的大气运动规律产生深刻的影响.在海岸城市地区,受城市冠层复杂结构的影响,热力内边界层呈现诸多独特的性质.采用大涡模拟的方法,研究了海岸城市热力内边界层内的流动特征.将海岸城市建筑物拟形为一系列有序排列的立方体块,并基于浸入边界方法考虑其对大气流动的影响.和常见的区域尺度数值研究中将城市下垫面假设为均质属性平面的方法相比,模拟结果的精度明显提高,海岸热力内边界层内各种尺度的动力过程也能在一定程度上得以描述.模拟计算了海岸城市热力内边界层、中性海岸城市边界层和开阔的自然海岸热力内边界层共3种情况,对比研究了城市摩阻和海岸热力作用对海岸城市上空大气运动的影响.研究结果表明:在空间上,边界层中各紊流特征量均在街谷尺度呈现出有规律而非均匀分布的特点;在强度上,热力作用与城市构筑物的摩阻作用相互促进,使得紊流强度显著提升并大于两种作用线性叠加的结果.此外,还发现有城市冠层海岸的热力内边界层发展远快于相同热力条件下的自然海岸带,边界层的发展规律也有显著的不同.      

基于局部移动Kriging无网格方法的层合板自由振动分析

利用基于局部移动Kriging插值无网格法对层合板自由振动进行了数值分析,基于一阶剪切层合理论导出了层合板振动的控制方程和边界条件,进一步得到了自由振动的离散化特征方程。由于Kriging插值函数具有Kronecker delta函数性质,可以直接施加本质边界条件。通过本文给出的方法,对不同边界条件、不同跨厚比、不同材料参数和铺设角度的层合板的振动频率进行了计算,均得到满意结果。最后用该方法对层合板的铺设角度进行优化设计,得到了与已有文献完全一致的优化结果。数值结果充分表明了无网格Kriging方法分析层合板自由振动问题的有效性和高精确度。     

On accurate boundary conditions for a shape sensitivity equation method

This paper studies the application of the continuous sensitivity equation method (CSEM) for the Navier–Stokes equations in the particular case of shape parameters. Boundary conditions for shape parameters involve flow derivatives at the boundary. Thus, accurate flow gradients are critical to the success of the CSEM. A new approach is presented to extract accurate flow derivatives at the boundary. High order Taylor series expansions are used on layered patches in conjunction with a constrained least‐squares procedure to evaluate accurate first and second derivatives of the flow variables at the boundary, required for Dirichlet and Neumann sensitivity boundary conditions. The flow and sensitivity fields are solved using an adaptive finite‐element method. The proposed methodology is first verified on a problem with a closed form solution obtained by the Method of Manufactured Solutions. The ability of the proposed method to provide accurate sensitivity fields for realistic problems is then demonstrated. The flow and sensitivity fields for a NACA 0012 airfoil are used for fast evaluation of the nearby flow over an airfoil of different thickness (NACA 0015). Copyright © 2005 John Wiley & Sons, Ltd.     

Large eddy simulation of free surface shallow‐water flow

Shallow‐water flow with free surface frequently occurs in ambient water bodies, in which the horizontal scale of motion is generally two orders of magnitude greater than the water depth. To accurately predict this flow phenomenon in more detail, a three‐dimensional numerical model incorporating the method of large eddy simulation (LES) has been developed and assessed. The governing equations are split into three parts in the finite difference solution: advection, dispersion and propagation. The advection part is solved by the QUICKEST scheme. The dispersion part is solved by the central difference method and the propagation part is solved implicitly using the Gauss–Seidel iteration method. The model has been applied to free surface channel flow for which ample experimental data are available for verification. The inflow boundary condition for turbulence is generated by a spectral line processor. The computed results compare favourably with the experimental data and those results obtained by using a periodic boundary condition. The performance of the model is also assessed for the case in which anisotropic grids and filters with horizontal grid size of the order of the water depth are used for computational efficiency. The coarse horizontal grid was found to cause a significant reduction in the large‐scale turbulent motion generated by the bottom turbulence, and the turbulent motion is predominately described by the sub‐grid scale (SGS) terms. The use of the Smagorinsky model for SGS turbulence in this situation is found inappropriate. A parabolic mixing length model, which accounts for the filtered turbulence, is then proposed. The new model can reproduce more accurately the flow quantities. Copyright © 2000 John Wiley & Sons, Ltd.     

复合材料旋转壳自由振动分析的新方法

提出了一种半解析区域分解法来分析任意边界条件的复合材料层合旋转壳自由振动. 沿壳体旋转轴线将壳体分解为一些自由的层合壳段, 视位移边界界面为一种特殊的分区界面;采用分区广义变分和最小二乘加权残值法将壳体所有分区界面上的位移协调方程引入到壳体的能量泛函中, 使层合壳的振动分析问题归结为无约束泛函变分问题. 层合壳段位移变量采用Fourier 级数和Chebyshev 多项式展开. 以不同边界条件的层合圆柱壳、圆锥壳及球壳为例, 采用区域分解法分析了其自由振动, 并将计算结果与其他文献值进行了对比. 算例表明, 该方法具有高效率、高精度和收敛性好等优点.