Effects of the Similarity Model in Finite-Difference LES of Isotropic Turbulence Using a Lagrangian Dynamic Mixed Model

A Lagrangian dynamic formulation of the mixed similarity subgrid (SGS) model for large-eddy simulation (LES) of turbulence is proposed. In this model, averaging is performed over fluid trajectories, which makes the model applicable to complex flows without directions of statistical homogeneity. An alternative version based on a Taylor series expansion (nonlinear mixed model) is also examined. The Lagrangian models are implemented in a finite difference code and tested in forced and decaying isotropic turbulence. As comparison, the dynamic Smagorinsky model and volume-averaged formulations of the mixed models are also tested. Good results are obtained, except in the case of low-resolution LES (32³) of decaying turbulence, where the similarity coefficient becomes negative due to the fact that the test-filter scale exceeds the integral scale of turbulence. At a higher resolution (64³), the dynamic similarity coefficient is positive and good agreement is found between predicted and measured kinetic energy evolution. Compared to the eddy viscosity term, the similarity or the nonlinear terms contribute significantly to both SGS dissipation of kinetic energy and SGS force. In order to dynamically test the accuracy of the modeling, the error incurred in satisfying the Germano identity is evaluated. It is found that the dynamic Smagorinksy model generates a very large error, only 3% lower than the worst-case scenario without model. Addition of the similarity or nonlinear term decreases the error by up to about 50%, confirming that it represents a more realistic parameterization than the Smagorinsky model alone.     

Resolution-optimised nonlinear scheme for secondary derivatives

A 5-point-stencil optimised nonlinear scheme with spectral-like resolution within the whole wave number range for secondary derivatives is devised. The proposed scheme can compensate for the dissipation deficiency of traditional linear schemes and suppress the spurious energy accumulation that occurs at high wave numbers, both of which are frequently encountered in large eddy simulation. The new scheme is composed of a linear fourth-order central scheme term and an artificial viscosity term. These two terms are connected by a nonlinear weight. The proposed nonlinear weight is designed based on Fourier analysis, rather than Taylor analysis, to guarantee a spectral-like resolution. Moreover, the accuracy is not affected by the optimisation, and the new scheme reaches fourth-order accuracy. The new scheme is tested numerically using the one-dimensional diffusion problem, one-dimensional steady viscous Burger’s shock, two-dimensional vortex decaying, three-dimensional isotropic decaying turbulence and fully developed turbulent channel flow. All the tests confirm that the new scheme has spectral-like resolution and can improve the accuracy of the energy spectrum, dissipation rate and high-order statistics of turbulent flows.     

Turbulence and cavitation models for time-dependent turbulent cavitating flows

Cavitation typically occurs when the fluid pressure is lower than the vapor pressure at a local thermodynamic state,and the flow is frequently unsteady and turbulent.To assess the state-of-the-art of computational capabilities for unsteady cavitating flows,different cavitation and turbulence model combinations are conducted.The selected cavitation models include several widely-used models including one based on phenomenological argument and the other utilizing interface dynamics.The kε turbulence model with additional implementation of the filter function and density correction function are considered to reduce the eddy viscosity according to the computed turbulence length scale and local fluid density respectively.We have also blended these alternative cavitation and turbulence treatments,to illustrate that the eddy viscosity near the closure region can significantly influence the capture of detached cavity.From the experimental validations regarding the force analysis,frequency,and the cavity visualization,no single model combination performs best in all aspects.Furthermore,the implications of parameters contained in different cavitation models are investigated.The phase change process is more pronounced around the detached cavity,which is better illus-trated by the interfacial dynamics model.Our study provides insight to aid further modeling development.     

可压缩多介质粘性流体和湍流的大涡模拟

在可压缩多介质粘性流体动力学高精度计算方法MVPPM(multi-viscous-fluid piecewise parabolicmethod)基础上,引入Smagorinsky和Vreman亚格子湍流模型,采用大涡数值模拟方法求解可压缩粘性流体NS(Navier-Stokes)方程,给出适用于可压缩多介质流体界面不稳定性发展演化至湍流阶段的计算方法和二维计算程序MVFT(multi-viscosity-fluid and turbulence)。在2种亚格子湍流模型下计算了LANL(Los Ala-mos National Laboratory)激波管单气柱RM不稳定性实验,分析了气柱的形状、流场速度以及涡的特征,通过与LANL实验和计算结果的比较可知,Vreman模型略优于Smagorinsky模型,MVFT方法和计算程序可用于对界面不稳定性发展演化至湍流阶段的数值模拟。     

Effective eddy viscosities in implicit modeling of decaying high Reynolds number turbulence with and without rotation

We assess the applicability of the numerical dissipation as an implicit turbulence model. The nonoscillatory finite volume numerical scheme MPDATA developed for simulations of geophysical flows is employed as an example of a scheme with an implicit turbulence model. A series of low resolution simulations of decaying homogeneous turbulence with and without Coriolis forces in the limit of zero molecular viscosity are performed. To assess the implicit model the long-time evolution of turbulence in the simulations is investigated and the numerical velocity fields are analyzed to determine the effective spectral eddy viscosity that is attributed to the numerical discretization. The detailed qualitative and quantitative comparisons are made between the numerical eddy viscosity and the theoretical results as well as the intrinsic eddy viscosity computed exactly from the velocity fields by introducing an artificial wave number cutoff. We find that the numerical dissipation depends on the time step and exhibits contradictory dependence on rotation: it is overestimated for rapid rotation cases and is underestimated for nonrotating cases. These results indicate that the numerical dissipation may fail to represent the effects of the physical subgrid scale processes unless the parameters of a numerical scheme are carefully chosen.     

A localised dynamic closure model for Euler turbulence

ABSTRACT

In this work, we present a localised form of the dynamic eddy viscosity model for computationally efficient and accurate simulation of the turbulent flows governed by Euler equations. In our framework, we determine the dynamic model coefficient locally using the information from neighbouring grid points through a test filtering process. We then develop an optimised Gaussian filtering kernel, using a consistent definition with respect to the test filtering ratio, which gives full attenuation at the grid cut-off wave number. A systematic a-posteriori analysis of our model is performed by solving two 3D test problems: (i) incompressible Taylor–Green vortex flow and (ii) compressible shear layer turbulence induced by Kelvin–Helmholtz instability to show the wide range of applicability of the proposed localised dynamic model. We demonstrate that the proposed dynamic model is robust and provides a better estimation of the inertial range turbulence dynamics than other numerical models tested in this study.     

Turbulence modeling for the axially rotating pipe from the viewpoint of analytical closures

A single-point model eddy viscosity model of rotation effects on the turbulent flow in an axially rotating pipe is developed based on two-point closure theories. Rotation is known to impede energy transfer in turbulence; this fact is reflected in the present model through a reduced eddy viscosity, leading to laminarization of the mean velocity profile and return to a laminar friction law in the rapid rotation limit. This model is compared with other proposals including linear redistribution effects through the rapid pressure-strain correlation, Richardson number modification of the eddy viscosity in a model of non-rotating turbulence, and the reduction of turbulence through the suppression of near-wall production mechanisms. PACS 47.27.Eq, 47.32.-y     

A dynamic eddy-viscosity closure model for large eddy simulations of two-dimensional decaying turbulence

This paper puts forth a simplified dynamic eddy-viscosity subgrid-scale model for the vorticity transport equation which is employed in a large eddy simulation study of freely evolving isotropic two-dimensional turbulent flows. The dynamic parameter is averaged in space, thereby retrieving a spatially constant value which only varies in time. The proposed dynamic model is applied to a two-dimensional decaying turbulence problem in a square periodic box, which is a standard prototype of more realistic turbulent flows in the atmosphere and oceans, in order to eliminate any possible errors associated with the boundary conditions or mesh non-uniformities. Compared with high-resolution direct numerical simulations, the performance of the dynamic model is systematically investigated considering various filtering strategies by means of test filters. The effects of the computational resolution and the filtering ratio between the test and the grid filters are also studied by using a huge set of parameters.     

On the eddy viscosity model of periodic turbulent shear flows

Physical argument shows that eddy viscosity is essentially different from molecular viscosity. By direct numerical simulation, it was shown that for periodic turbulent flows, there is phase difference between Reynolds stress and rate of strain. This finding posed great challenge to turbulence modeling, because most turbulence modeling, which use the idea of eddy viscosity, do not take this effect into account. The project supported by the National Natural Science Foundation of China (19732005) and Liu Hui Center for Applied Mathematics of Nankai & Tianjin University     

基于涡粘性模型的翼型湍流流场熵产率计算

利用有限体积法实现了基于非正交同位网格的SIMPLE算法。基于熵分析方法,采用涡粘性模型求解湍流熵产方程,系统研究了湍流模型对二维翼型绕流流场熵产率的影响。通过计算NACA0012翼型在来流雷诺数为2.88×10⁶时,0°攻角~16.5°攻角范围内的翼型表面压力系数分布和升阻力特性,验证了算法及程序的正确性。结果表明,选择不同湍流模型时,翼型流场熵产的计算结果存在差异,湍流耗散是引起流场熵产的主要原因;翼型流场的熵产主要发生在翼型前缘区、壁面边界层和翼型尾流区域,流场熵产率与翼型阻力系数线性相关;当产生分离涡时,粘性耗散引起的熵产下降。     

Residual‐based artificial viscosity for simulation of turbulent compressible flow using adaptive finite element methods

In this paper, we present a finite element method with a residual‐based artificial viscosity for simulation of turbulent compressible flow, with adaptive mesh refinement based on a posteriori error estimation with sensitivity information from an associated dual problem. The artificial viscosity acts as a numerical stabilization, as shock capturing, and as turbulence capturing for large eddy simulation of turbulent flow. The adaptive method resolves parts of the flow indicated by the a posteriori error estimates but leaves shocks and turbulence under‐resolved in a large eddy simulation. The method is tested for examples in 2D and 3D and is validated against experimental data. Copyright © 2012 John Wiley & Sons, Ltd.     

An eddy viscosity model with near‐wall modifications

An extended version of the isotropic k–ε model is proposed that accounts for the distinct effects of low‐Reynolds number (LRN) and wall proximity. It incorporates a near‐wall correction term to amplify the level of dissipation in nonequilibrium flow regions, thus reducing the kinetic energy and length scale magnitudes to improve prediction of adverse pressure gradient flows, involving flow separation and reattachment. The eddy viscosity formulation maintains the positivity of normal Reynolds stresses and the Schwarz' inequality for turbulent shear stresses. The model coefficients/functions preserve the anisotropic characteristics of turbulence. The model is validated against a few flow cases, yielding predictions in good agreement with the direct numerical simulation (DNS) and experimental data. Comparisons indicate that the present model is a significant improvement over the standard eddy viscosity formulation. Copyright © 2005 John Wiley & Sons, Ltd.     

Turbulence modulation for gas - particle flow in vertical tube and horizontal channel using four-way Eulerian-Lagrangian approach

Turbulence modulation of gas-solid flow in vertical tube and horizontal channel in dilute and moderately dense suspensions is investigated numerically using a four way Eulerian-Lagrangian approach. Low Reynolds number k-l model is used for analyzing the fluid phase motion. A new model is presented based on a source-term formulation, which can predict fluid phase turbulence augmentation due to the presence of large particles and damping of turbulence due to small particles in the core of the channel and tube. Particle-particle and particle-wall collisions are simulated based on a deterministic approach, and coupling terms representing the fluid-particle interactions are also taken into account. The predicted fluid mean velocity and turbulence intensity profiles are in good agreement with the available experimental data. Additional numerical simulation results for variation of the eddy viscosity, turbulence production and dissipation are also presented for different values of loading ratios.     

Assessment of Subgrid-Scale Models for Decaying Compressible MHD Turbulence

Large eddy simulation method is formulated for study of compressible magnetohydrodynamic turbulence and assessment of different subgrid-scale models as applied for decaying case is performed. The filtered equations of compressible magnetopause using the mass-weighted filtering operation are obtained. Mass-weighted filtered equations for large-scale turbulent component comprise subgrid-scale terms and five models for closure of the subgrid-scale terms are suggested. In present paper the obtained results of numerical computations for large eddy simulation are compared with the results of direct numerical simulation of three-dimensional compressible magnetohydrodynamic turbulence. Assessment of five subgrid-scale models of large eddy simulation for MHD flows is fulfilled. The comparisons between large eddy simulation and direct numerical simulation are carried out regarding the temporal evolution of the global quantities kinetic and magnetic energy, cross helicity and the spectra of kinetic and magnetic energy.     

Nonlinear turbulence models for predicting strong curvature effects

Prediction of the characteristics of turbulent flows with strong streamline curvature, such as flows in turbomachines, curved channel flows, flows around airfoils and buildings, is of great importance in engineering applications and poses a very practical challenge for turbulence modeling. In this paper, we analyze qualitatively the curvature effects on the structure of turbulence and conduct numerical simulations of a turbulent Uduct flow with a number of turbulence models in order to assess their overall performance. The models evaluated in this work are some typical linear eddy viscosity turbulence models, nonlinear eddy viscosity turbulence models （NLEVM） （quadratic and cubic）, a quadratic explicit algebraic stress model （EASM） and a Reynolds stress model （RSM） developed based on the second-moment closure. Our numerical results show that a cubic NLEVM that performs considerably well in other benchmark turbulent flows, such as the Craft, Launder and Suga model and the Huang and Ma model, is able to capture the major features of the highly curved turbulent U-duct flow, including the damping of turbulence near the convex wall, the enhancement of turbulence near the concave wall, and the subsequent turbulent flow separation. The predictions of the cubic models are quite close to that of the RSM, in relatively good agreement with the experimental data, which suggests that these models may be employed to simulate the turbulent curved flows in engineering applications.     

Effects of Turbulence on Self-sustained Combustion in Premixed Flame Kernels: A Direct Numerical Simulation (DNS) Study

The effects of mean flame radius and turbulence on self-sustained combustion of turbulent premixed spherical flames in decaying turbulence have been investigated using three-dimensional direct numerical simulations (DNS) with single step Arrhenius chemistry. Several flame kernels with different initial radius or initial turbulent field have been studied for identical conditions of thermo-chemistry. It has been found that for very small kernel radius the mean displacement speed may become negative leading ultimately to extinction of the flame kernel. A mean negative displacement speed is shown to signify a physical situation where heat transfer from the kernel overcomes the heat release due to combustion. This mechanism is further enhanced by turbulent transport and, based on simulations with different initial turbulent velocity fields, it has been found that self-sustained combustion is adversely affected by higher turbulent velocity fluctuation magnitude and integral length scale. A scaling analysis is performed to estimate the critical radius for self-sustained combustion in premixed flame kernels in a turbulent environment. The scaling analysis is found to be in good agreement with the results of the simulations.     

Modelling turbulence around and inside porous media based on the second moment closure

To predict turbulence in porous media, a new approach is discussed. By double (both volume and Reynolds) averaging Navier–Stokes equations, there appear three unknown covariant terms in the momentum equation. They are namely the dispersive covariance, the macro-scale and the micro-scale Reynolds stresses, in the present study. For the macro-scale Reynolds stress, the TCL (two-component-limit) second moment closure is applied whereas the eddy viscosity models are applied to the other covariant terms: the Smagorinsky model and the one-equation eddy viscosity model, respectively for the dispersive covariance and the micro-scale Reynolds stress. The presently proposed model is evaluated in square rib array flows and porous wall channel flows with reasonable accuracy though further development is required.     

剪切湍流的涡黏张量模式

对剪切湍流提出了涡黏系数为四阶张量的涡黏张量模式。引入近代数学中Ｍｏｏｒｅ－Ｐｅｎｒｏｓｅ广义逆矩阵的研究结果，给出了构造涡黏张量各分量的计算公式。用平面后台阶流动验算了剪切湍流的涡黏张量模式，比ＲＳＭ和ｋ－ε模式更接近实验结果。提出了剪切湍流涡黏张量模式的应用设想。     

Hybrid URANS/LES Turbulence Simulation of Vortex Shedding Behind a Triangular Flameholder

In this study a detached eddy simulation (DES) model, which belongs to the group of hybrid URANS/LES turbulence models, is used for the simulation of vortex shedding behind a triangular obstacle. In the near wall region or in regions where the grid resolution is not sufficiently fine to resolve smaller structures, the two-equation RANS shear-stress transport (SST) model is used. In the other regions with higher grid resolution a LES model, which uses a transport equation for the turbulent subgrid energy, is applied. The DES model is first investigated for two standard test cases, namely decaying homogeneous isotropic turbulence and the backward facing step, respectively. For the decaying homogeneous isotropic turbulence test case the evolution of the energy spectra in wavenumber space for different times are studied for both the DES and a Smagorinsky type LES model. Different grid resolutions are analyzed with a special emphasis on the modeling constant connecting the filter length scale to the grid size. The results are compared to experimental data. The backward facing step test case is used to study the model behavior for a case with a transition region between a RANS modeling approach close to the wall and LES based modeling in the intense shear flow region. The final application is the simulation of the vortex shedding behind a triangular obstacle. First, the influence of the inlet condition formulation is studied in detail as they can have a significant influence especially for LES based models. Detailed comparisons between simulation and experiment for the flow structure past the obstacle and statistical quantities such as the shedding frequency are shown. Finally the additional temporal and spatial information provided by the DES model is used to show the predicted anisotropy of turbulence.