Subgrid-modelling in LES of compressible flow

Subgrid-models for Large Eddy Simulation (LES) of compressible turbulent flow are tested for the three-dimensional mixing layer. For the turbulent stress tensor the recently developed dynamic mixed model yields reasonable results.A priori estimates of the subgrid-terms in the filtered energy equation show that the usually neglected pressure-dilatation and turbulent dissipation rate are as large as the commonly retained pressure-velocity subgrid-term. Models for all these terms are proposed: a similarity model for the pressure-dilatation, similarity andk-dependent models for the turbulent dissipation rate and a dynamic mixed model for the pressure-velocity subgrid-term. Actual LES demonstrates that for a low Mach number all subgrid-terms in the energy equation can be neglected, while for a moderate Mach number the effect of the modelled turbulent dissipation rate is larger than the combined effect of the other modelled subgrid-terms in the filtered energy equation.     

Large eddy simulation of fully developed turbulent flow in a curved channel

In this paper large eddy simulation of the fully developed turbulent flow in a curved channel is carried out. The computational results are presented and compared with the experimental results of Eskinazi and Yeh^[1]. It is shown that the numerical results of the present LES are reliable and the influence of the curvature on the turbulence feature is correctly revealed.     

An immersed boundary finite difference method for LES of flow around bluff shapes

A three‐dimensional numerical model using large eddy simulation (LES) technique and incorporating the immersed boundary (IMB) concept has been developed to compute flow around bluff shapes. A fractional step finite differences method with rectilinear non‐uniform collocated grid is employed to solve the governing equations. Bluff shapes are treated in the IMB method by introducing artificial force terms into the momentum equations. Second‐order accurate interpolation schemes for all sorts of grid points adjacent to the immersed boundary have been developed to determine the velocities and pressure at these points. To enforce continuity, the methods of imposition of pressure boundary condition and addition of mass source/sink terms are tested. It has been found that imposing suitable pressure boundary condition (zero normal gradient) can effectively reproduce the correct pressure distribution and enforce mass conservation around a bluff shape. The present model has been verified and applied to simulate flow around bluff shapes: (1) a square cylinder and (2) the Tsing Ma suspension bridge deck section model. Complex flow phenomena such as flow separation and vortex shedding are reproduced and the drag coefficient, lift coefficient, and pressure coefficient are calculated and analyzed. Good agreement between the numerical results and the experimental data are obtained. The model is proven to be an efficient tool for flow simulation around bluff bodies in time varying flows. Copyright © 2004 John Wiley & Sons, Ltd.     

Large-eddy simulation of stratified channel flow

An investigation of large-eddy simulation (LES) for turbulent channel flow with buoyancy effects was performed by solving the resolved incompressible Navier-Stokes equations under the Boussinesq approximation. The Smagorinsky eddy-viscosity model and Yoshizawa eddy-viscosity model were used to describe the unresolved subgrid scale (SGS) fluctuations respectively. After some numerical testing, the latter was further simplified so that it can be used in the dynamic model closure. A LES code was developed for parallel computations by using the parallel technique, and was run on the Dawn-1000 parallel computer. To demonstrate the viability and accuracy of the code, our results are compared with and found in good agreement with available LES results. The project supported by the National Natural Science Foundation of China and by the Youngster Funding of Academia Sinica     

A hybrid RANS/LES simulation of high-Reynolds-number channel flow using additional filtering at the interface

A hybrid method combining large eddy simulation (LES) with the Reynolds-averaged Navier-Stokes (RANS) equation is used to simulate a turbulent channel flow at high Reynolds number. It is known that the mean velocity profile has a mismatch between the RANS and LES regions in hybrid simulations of a channel flow. The velocity mismatch is reproduced and its dependence on the location of the RANS/LES interface and on the type of RANS model is examined in order to better understand its properties. To remove the mismatch and to obtain better velocity profiles, additional filtering is applied to the velocity components in the wall-parallel planes near the interface. The additional filtering was previously introduced to simulate a channel flow at low Reynolds number. It is shown that the filtering is effective in reducing the mismatch even at high Reynolds number. Profiles of the velocity fluctuations of runs with and without the additional filtering are examined to help understand the reason for the mismatch. Due to the additional filtering, the wall-normal velocity fluctuation increases at the bottom of the LES region. The resulting velocity field creates the grid-scale shear stress more efficiently, and an overestimate of the velocity gradient is removed. The dependence of the velocity profile on the grid point number is also investigated. It is found that the velocity gradient in the core region is underestimated in the case of a coarse grid. Attention should be paid not only to the velocity mismatch near the interface but also to the velocity profile in the core region in hybrid simulations of a channel flow at high Reynolds number. PACS47.27.Eq; 47.27.Nz; 47.60.+i     

Large eddy simulation in a turbulent channel flow using exit boundary conditions

The influence of the exit boundary conditions on the vanishing first derivative of the velocity components and constant pressure on the large eddy simulation of the fully developed turbulent channel flow has been investigated for equidistant and stretched grids at the channel exit. Results show that the chosen exit boundary conditions introduce some small disturbances that are mostly damped by the grid stretching. The difference of rms values between the fully developed turbulent channel flow with periodicity conditions and the fully developed channel flow using inlet and the exit boundary conditions is less than 10% for the equidistant grids and less than 5% for the stretched grids. The chosen boundary conditions are of interest because they may be used in complex problems with back flow at the exit. Copyright © 1999 John Wiley & Sons, Ltd.     

An immersed boundary method wall model for high‐Reynolds‐number channel flow over complex topography

High‐Reynolds‐number channel flows regularly encounter topographies composed of multiple length scales and that protrude into the boundary layer. Physically, the presence of immersed obstacles leads to increased velocity gradients, turbulence production, and manifestation of wakes. Considerable challenges are associated with numerically describing the presence of obstacles in channel flows. Common approaches include generation of a computational mesh that is uniquely designed for the flow and obstacle, the immersed boundary method, and terrain‐following coordinates. There are challenges and limitations associated with each of these techniques. Specification of boundary conditions representing the perimeter of solid obstacles is a primary challenge of the immersed boundary method. In this document, a simplistic canopy stress‐like wall model is used to impose boundary conditions. The model isolates aerodynamically relevant local frontal areas through evaluation of the gradient of the topographic height field. The gradient of the height field describes both the surface‐normal direction and the frontal area, making it ideal for detecting areas on which the flow impinges. The model is tested in numerical simulations of turbulent half‐channel flow over topographies with different obstacles affixed–right prisms, rectangular prisms, ellipsoidal mounds, and sinusoids. In all cases, the performance is strong relative to datasets presented in the literature. Results are finally presented for numerical simulation of flow over complex synthetic fractal‐like topography and a synthetic city. These results show interesting trends in how the turbulent multiscale flow field responds to multiscale topography. Copyright © 2012 John Wiley & Sons, Ltd.     

Grid-averaged Lagrangian LES model for multiphase turbulent flow

A grid-averaged Lagrangian (GAL) model for dispersed particle motion in multiphase turbulent flow is presented to provide a large eddy simulation (LES) model for multiphase turbulent flow in which a quite large number of particles are involved. The GAL model is based on an averaging operation for a Lagrangian-type equation of motion of a particle over a computational grid volume and a procedure of reallocation of a dispersed particle cloud with its centroid movement to each grid. The model is therefore a mixed Eulerian–Lagrangian model which can effectively reduce computational time compared with existing Lagrangian-type models, without losing the advantage of Lagrangian-type models that they can properly describe the dynamical evolution of particles. Since the GAL model adopts the grid-volume averaging operation it can easily provide an effective SGS model for LES modeling of multiphase turbulent flow. The validity of the multiphase LES model developed, which is named the GAL-LES model, is confirmed through its application to a particle plume, in which the present model is found to simulate large-eddy motion usually observed in a jet and plume, and to give good agreements with experimental data.     

Shear‐improved Smagorinsky modeling of turbulent channel flow using generalized Lattice Boltzmann equation

Generalized Lattice Boltzmann equation (GLBE) was used for computation of turbulent channel flow for which large eddy simulation (LES) was employed as a turbulence model. The subgrid‐scale turbulence effects were simulated through a shear‐improved Smagorinsky model (SISM), which is capable of predicting turbulent near wall region accurately without any wall function. Computations were done for a relatively coarse grid with shear Reynolds number of 180 in a parallelized code. Good numerical stability was observed for this computational framework. The results of mean velocity distribution across the channel showed good correspondence with direct numerical simulation (DNS) data. Negligible discrepancies were observed between the present computations and those reported from DNS for the computed turbulent statistics. Three‐dimensional instantaneous vorticity contours showed complex vortical structures that appeared in such flow geometries. It was concluded that such a framework is capable of predicting accurate results for turbulent channel flow without adding significant complications and the computational cost to the standard Smagorinsky model. As this modeling was entirely local in space it was therefore adapted for parallelization. Copyright © 2010 John Wiley & Sons, Ltd.     

基于非嵌入多项式混沌的爆轰不确定度量化

固体炸药爆轰是极为复杂的极端物理和化学过程,发生在极小的时间尺度和极小的空间尺度内,产生极高的爆压和极快的爆速,给理论研究和实验带来巨大的挑战,数值模拟成为极为重要的研究手段。验证、确认和不确定度量化综合了实验和数值模拟的优点。针对圆筒试验问题,使用基于全张量积的多元非嵌入多项式混沌方法,给出了初始密度和JWL参数不确定时、圆筒位置和速度的期望、方差以及置信区间。结果发现试验数据均落在置信区间内。最后,使用信息熵方法给出实验不确定度量化,从而确认了模型的有效性。     

LES-based filter-matrix lattice Boltzmann model for simulating fully developed turbulent channel flow

In this paper, a three-dimensional filter-matrix lattice Boltzmann (FMLB) model based on large eddy simulation (LES) was verified for simulating wall-bounded turbulent flows. The Vreman subgrid-scale model was employed in the present FMLB–LES framework, which had been proved to be capable of predicting turbulent near-wall region accurately. The fully developed turbulent channel flows were performed at a friction Reynolds number Re_τ of 180. The turbulence statistics computed from the present FMLB–LES simulations, including mean stream velocity profile, Reynolds stress profile and root-mean-square velocity fluctuations greed well with the LES results of multiple-relaxation-time (MRT) LB model, and some discrepancies in comparison with those direct numerical simulation (DNS) data of Kim et al. was also observed due to the relatively low grid resolution. Moreover, to investigate the influence of grid resolution on the present LES simulation, a DNS simulation on a finer gird was also implemented by present FMLB–D3Q19 model. Comparisons of detailed computed various turbulence statistics with available benchmark data of DNS showed quite well agreement.     

Assessment of stretched vortex subgrid‐scale models for LES of incompressible inhomogeneous turbulent flow

The physical space version of the stretched vortex subgrid scale model is tested in LES of the turbulent lid‐driven cubic cavity flow. LES is carried out by using a higher order finite‐difference method. The effects of different vortex orientation models and subgrid turbulence spectrums are assessed through comparisons of the LES predictions against DNS. Three Reynolds numbers 12000, 18000, and 22000 are studied. Good agreement with the DNS data for the mean and fluctuating quantities is observed. Copyright © 2013 John Wiley & Sons, Ltd.     

Large eddy simulation in a fully developed turbulent flow in a channel and comparison of subgrid eddy viscosity models

The accuracy and computational efficiency are compared for a number of models of subgrid eddy viscosity (Smagorinsky model, renormalization group model, and dynamic and one-parameter models). Space-filtered Navier-Stokes equations are solved numerically by the control-volume approach on a nonuniform grid with the use of high-resolution schemes in time and space. The numerical data are compared with the results of a physical experiment and direct numerical simulation. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 47, No. 3, pp. 31–42, May–June, 2006.     

超声速来流横向狭缝喷流标量输运的混合RANS/LES模拟

采用基于时均湍流模型的混合RANS/LES模拟方法对超声速来流狭缝横向喷流标量输运过程进行大涡模拟,控制方程对流项采用五阶精度的WENO格式求解,过滤后的组分方程中的亚格子组分对流通量采用梯度扩散模化.模拟得到了超声速来流狭缝横向喷流大尺度涡和组分的演化过程,研究表明喷流具有明显的非定常周期性特点,并且组分随时间的空间分布特征决定于流动的大尺度结构.模拟得到的压力的统计时均结果和实验结果定性一致,但喷流前后的分离区较实验值过大,这表明混合RANS/LES需要进一步改进.     

Non-Equilibrium Mixing of Turbulence Scales Using a Continuous Hybrid RANS/LES Approach: Application to the Shearless Mixing Layer

A new subgrid-scale (SGS) model based on partially integrated transport method (PITM) is applied to the case of a turbulent spectral non-equilibrium flow created by the mixing of two turbulence fields of differing scales: the shearless mixing layer. The method can be viewed as a continuous hybrid RANS/LES approach. In this model the SGS length scale is no longer given by the size of the discretization step, but is dynamically estimated using an additional transport equation for the dissipation rate. The results are compared to those corresponding to the classical model of Smagorinsky and to the experimental data of Veeravalli and Warhaft. A method for creating an anisotropic analytical pseudo-random field for inflow conditions is also proposed. This approach based on subgrid-scale transport modelling combined with anisotropic inlet conditions gives better results for the prediction of the shearless mixing layer.     

Investigation of Anisotropy-Resolving Turbulence Models by Reference to Highly-Resolved LES Data for Separated Flow

The predictive properties of several non-linear eddy-viscosity models are investigated by reference to highly-resolved LES data obtained by the authors for an internal flow featuring massive separation from a curved surface. The test geometry is a periodic segment of a channel constricted by two-dimensional (2D) `hills' on the lower wall. The mean-flow Reynolds number is 21560. Periodic boundary conditions are applied in the streamwise and spanwise directions. This makes the statistical properties of the simulated flow genuinely 2D and independent from boundary conditions, except at the walls. The simulation was performed on a high-quality, 5M-node grid. The focus of the study is on the exploitation of the LES data for the mean-flow, Reynolds stresses and macro-length-scale. Model solutions are first compared with the LES data, and selected models are then subjected to a-priori studies designed to elucidate the role of specific model fragments in the non-linear stress-strain/vorticity relation and their contribution to observed defects in the mean-flow and turbulence fields. The role of the equation governing the length-scale, via different surrogate variables, is also investigated. It is shown that, while most non-linear models overestimate the separation region, due mainly to model defects that result in insufficient shear stress in the separated shear layer, model forms can be derived which provide a satisfactory representation of the flow. One such model is identified. This combines a particular quadratic constitutive relation with a wall-anisotropy term, a high-normal-strain correction and a new form of the equation for the specific dissipation ω = ∈/k. This revised version was published online in July 2006 with corrections to the Cover Date.     

风口区空气流动的有限元模拟与实验研究

通风和空调出风口的射流特性主要取决于出风口截面之前的一段管道内的流动状况,研究送风口的射流偏转特性有重要实用价值.应用湍流大涡模拟技术结合Taylor-Galerkin有限单元法对工程中常见工况下的风口区管道内的三维流场进行了数值模拟和实验研究,数值计算结果与实验结果符合较好.表明湍流大涡模拟方法适合于边界形状复杂,存在各向异性的大尺度涡的内流情况,能可靠地预报风口区空气流动的射流偏转特性.