An improved dynamic subgrid-scale model and its application to large eddy simulation of stratified channel flows

In the present paper, a new dynamic subgrid-scale (SGS) model of turbulent stress and heat flux for stratified shear flow is proposed. Based on our calculated results of stratified channel flow, the dynamic subgrid-scale model developed in this paper is shown to be effective for large eddy simulation (LES) of stratified turbulent shear flows. The new SGS model is then applied to the LES of the stratified turbulent channel flow to investigate the coupled shear and buoyancy effects on the behavior of turbulent statistics, turbulent heat transfer and flow structures at different Richardson numbers.     

A dynamic subgrid-scale model for the large eddy simulation of stratified flow

A new dynamic subgrid-scale (SGS) model, including subgrid turbulent stress and heat flux models for stratified shear flow is proposed by using Yoshizawa’s eddy viscosity model as a base model. Based on our calculated results, the dynamic subgrid-scale model developed here is effective for the large eddy simulation (LES) of stratified turbulent channel flows. The new SGS model is then applied to the large eddy simulation of stratified turbulent channel flow under gravity to investigate the coupled shear and buoyancy effects on the near-wall turbulent statistics and the turbulent heat transfer at different Richardson numbers. The critical Richardson number predicted by the present calculation is in good agreement with the value of theoretical analysis     

Stochastic parameterization for large eddy simulation of geophysical flows

Recently, stochastic, as opposed to deterministic, parameterizations are being investigated to model the effects of unresolved subgrid scales (SGS) in large eddy simulations (LES) of geophysical flows. We analyse such a stochastic approach in the barotropic vorticity equation to show that (i) if the stochastic parameterization approximates the actual SGS stresses, then the solution of the stochastic LES approximates the ``true" solution at appropriate scale sizes; and that (ii) when the filter scale size approaches zero, the solution of the stochastic LES approaches the true solution.

     

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.     

Direct Numerical Simulation of turbulent heat transfer behind a backward-facing step at low Prandtl number

The paper presents Direct Numerical Simulations of the turbulent flow of a low Prandtl number fluid over a backward-facing step with heat transfer. The backward-facing step flow is investigated as a generic configuration for sudden changes in cross section. Several simulations are reported: for isothermal conditions, for heat transfer with the Prandtl number of air, and for heat transfer with the Prandtl number of liquid sodium. The simulation for air is compared to results from literature. The differences induced by reduction of the Prandtl number are then assessed by comparison of the two cases. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Notes on modelling of environmental transport in wetland

This paper presents an effort towards a basic model for environmental transport of momentum, heat and mass transfer in the wetland. To smear out the discontinuity between the two phases of water and solid in the wetland, the continuum models distinctively applying for the water body and solid frame are transformed via the technique of phase average to give equations for a virtual single-phase flow in the entire domain of the wetland. Then to filter out the vortex and fluctuation common in the wetland flow, the operation of large eddy simulation (LES) is applied to yield a basic model for practical simulation. With reference to the modelling of flows in porous media and turbulent flows, closure relations are presented for the correlation terms due to the phase average and large eddy simulation.     

On large Eddy simulation and variational multiscale methods in the numerical simulation of turbulent incompressible flows

Numerical simulation of turbulent flows is one of the great challenges in Computational Fluid Dynamics (CFD). In general, Direct Numerical Simulation (DNS) is not feasible due to limited computer resources (performance and memory), and the use of a turbulence model becomes necessary. The paper will discuss several aspects of two approaches of turbulent modeling—Large Eddy Simulation (LES) and Variational Multiscale (VMS) models. Topics which will be addressed are the detailed derivation of these models, the analysis of commutation errors in LES models as well as other results from mathematical analysis.     

Numerical Study on Hypersonic Flow and Aerodynamic Heating北大核心CSCD

基于有限差分法开发了高超声速流动与换热问题气热耦合仿真求解器,运用该求解器对三种典型高超声速流动与换热问题开展了仿真研究,得到了相应的气动参数、热流密度分布。高超声速后台阶的存在使表面气动参数、热流分布不再连续;随着缝深的提高,缝隙局部流速迅速降低,对流换热效应减弱;高超声速无限长圆管绕流中,边界层外部区域气动参数随时间变化不大,边界层内存在较大的温度梯度,壁面温度随时间升高。三个算例的仿真结果均与试验测量值进行了对比,验证了所开发的求解器的计算能力。     

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.     

Comparison of isotropic and anisotropic subgrid scale models in Large Eddy Simulations of a backward-facing step and square duct flow

The present paper employs the anisotropic explicit algebraic subgrid-stress model (EASM), proposed by Marstorp et al. J. Fluid. Mech. 639, 403 - 432 (2009), developed in the spirit of an earlier explicit algebraic RANS-model for statistical closures. The EASM in its elementary, computationally efficient non-dynamic version is used for Large Eddy Simulation of three-dimensional flows over a backward-facing step at a bulk Reynolds number of 4805 and a square duct at 2205. Its performance is assessed by comparison with experimental data and an own Direct Numerical Simulation. Furthermore, a set of eddy-viscosity models, including the recent σ-model, is employed for comparison. Various statistical quantities are evaluated to assess the respective performance of the different models showing, that the anisotropic EASM compares favorably to the other models. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A rate-dependent algebraic stress model for turbulence

Based on the stress transport model, a rate-dependent algebraic expression for the Reynolds stress tensor is developed. It is shown that the new model includes the normal stress effects and exhibits viscoelastic behavior. Furthermore, it is compatible with recently developed improved models of turbulence. The model is also consistent with the limiting behavior of turbulence in the inertial sublayer and is capable of predicting secondary flows in noncircular ducts. The TEACH code is modified according to the requirements of the rate-dependent model and is used to predict turbulent flow fields in a channel and behind a backward-facing step. The predicted results are compared with the available experimental data and those obtained from the standard k-ε and algebraic stress models. It is shown that the predictions of the new model are in better agreements with the experimental data.     

CFD of turbulent transport of particles behind a backward-facing step using a new model—k-ε-Sp

The k-ε-S_p model, describing two-dimensional gas–solid two-phase turbulent flow, has been developed. In this model, the diffusion flux and slip velocity of solid particles are introduced to represent the particle motion in two-phase flow. Based on this model, the gas–solid two-phase turbulent flow behind a vertical backward-facing step is simulated numerically and the turbulent transport velocities of solid particles with high density behind the step are predicted. The numerical simulation is validated by comparing the results of the numerical calculation with two other two-phase turbulent flow models (k-ε-A_p, k-ε-k_p) by Laslandes and the experimental measurements. This model, not only has the same virtues of predicting the longitudinal transport of the solid particles as the present practical two-phase flow models, but also can predict the lateral transport of the solid particles correctly.     

湍流不同尺度间的近程和共振作用规律及应用

从Navier-Stokes方程出发,研究了湍流不同尺度间的相互作用规律,给出相近尺度间近程粘性应力的积分和微分表达式．引入极相近尺度之间共振相互作用的概念,得到共振粘性应力的微分表达式．利用共振粘性应力张量获得不含经验关系和常数、近似封闭的大涡模拟（LES）方程组．利用近程和共振粘性应力张量获得不含经验关系和常数、近似封闭的湍流多尺度方程组．讨论了湍流多尺度方程的性质及用于湍流计算的优点,尺度间相互作用的近程特性说明:多尺度模拟是湍流计算很有价值的方法,并列举了算例．     

A dynamical approach to large eddy simulation of turbulent flows: existence of weak solutions

We consider a system of equations coming from turbulence models using a large eddy simulation (LES) technique. The idea of this approach bases on decomposing the velocity into a part containing large flow structures and a part consisting of small scales. The equations for large‐scale quantities are derived from the Navier–Stokes equations with an additional constitutive relation for the contribution of small eddies. The mathematical difficulties in this paper focus on the non‐linear and non‐local turbulent term. Copyright © 2005 John Wiley & Sons, Ltd.     

The effect of modeling of velocity fluctuations on prediction of collection efficiency of cyclone separators

The effect of modeling of velocity fluctuations on the prediction of collection efficiency of cyclone separators has been numerically investigated using the Reynolds stress turbulence model (RSTM) and large eddy simulation (LES). The Eulerian–Lagrangian modeling approach of CFD code Fluent 6.3.26 has been employed to simulate the three dimensional, unsteady turbulent gas–solid flows in a Stairmand high efficiency cyclone. The simulated results have been compared with experimental observations available in the literature. The analysis of results shows that the RSTM and the LES have adequately predicted the mean flow field. Results of the present study demonstrate that the LES has good performance on prediction of fluctuating flow field and collection efficiency for each and every particle size. However, the performance of the RSTM is found poor in terms of prediction of velocity fluctuations and collection efficiency, especially for small particles. This relates to the precessing of the vortex core phenomenon, which is resolved more accurately by LES as compared to the RSTM simulation. The results suggest that the prediction of collection efficiency, especially for small particles is greatly influenced by the simulation of velocity fluctuations in cyclones.     

Computational fluid dynamics simulation of fluid particle fragmentation in turbulent flows

A simulation methodology is presented that allows detailed studies of the breakup mechanism of fluid particles in turbulent flows. The simulations, based on large eddy and volume of fluid simulations, agree very well with high-speed measurements of the breakup dynamics with respect to deformation time and length scales, and also the resulting size of the daughter fragments. The simulations reveal the size of the turbulent vortices that contribute to the breakup and how fast the interaction and energy transfer occurs. It is concluded that the axis of the deformed particle and the vortex core axis are aligned perpendicular to each other, and that breakup sometimes occurs due to interaction with two vortices at the same time. Analysis of the energy transfer from the continuous phase turbulence to the fluid particles reveals that the deformed particle attains it maximum in interfacial energy before the breakup is finalized. Similar to transition state theory in chemistry this implies that an activation barrier exists. Consequently, by considering the dynamics of the phenomenon, more energy than required at the final stage needs to be transferred from the turbulent vortices for breakup to occur. This knowledge helps developing new, more physical sound models for the breakup phenomenon required to solve scale separation problems in computational fluid dynamics simulations.     

Implicit subgrid‐scale modeling by adaptive local deconvolution

A class of implicit Subgrid‐Scale (SGS) models for Large‐Eddy Simulation (LES) is obtained from a new approach for the finite‐volume discretization of hyperbolic conservation laws. The extension of a standard deconvolution operator and the choice of a suitable numerical flux function result in a truncation error which can be forced to have physical significance. As determined by a modified‐differential‐equation analysis, the free parameters of this implicit SGS model can be adjusted to approximate a known explicit SGS model. Computational results for the viscous Burgers equation show that the model with parameters identified by evolutionary optimization gives significantly better predictions than other models. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

填充烧结铜球的T型管中流体混合的大涡模拟

在FLUENT软件平台上,运用大涡模拟湍流模型及Smagorinsky-Lilly亚格子尺度模型,对填充有烧结铜球多孔介质的T型管道内冷热流体混合过程的流动与传热情况进行了数值计算,与未填充多孔介质时混合区域内的平均温度和温度波动、平均速度和速度波动等数据进行了对比,并对温度波动进行了功率谱密度分析．数值结果表明,多孔介质可有效削弱T型通道流体混合区域内的温度和速度波动,有效降低1 Hz至10 Hz频域中的温度波动的功率谱密度．     

Artificial viscosity proper orthogonal decomposition

We introduce improved reduced-order models for turbulent flows. These models are inspired from successful methodologies used in large eddy simulation, such as artificial viscosity, applied to standard models created by proper orthogonal decomposition of flows coupled with Galerkin projection. As a first step in the analysis and testing of our new methodology, we use the Burgers equation with a small diffusion parameter. We present a thorough numerical analysis for the time discretization of the new models. We then test these models in two problems displaying shock-like phenomena. Of course, since the Burgers equation does not model turbulence, we next need to test our new models in realistic turbulent flow settings. This is the subject of a forthcoming report.