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

在可压缩多介质粘性流体动力学高精度计算方法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方法和计算程序可用于对界面不稳定性发展演化至湍流阶段的数值模拟。     

Assessment of LES Subgrid-scale Models and Investigation of Hydrodynamic Behaviour for an Axisymmetrical Bluff Body Flow

This work is concerned with the investigation of fluid-mechanical behaviour and the performance of different subgrid-scale models for LES in the numerical prediction of a confined axisymmetrical bluff-body flow. Four subgrid-scale turbulence models comprising the Smagorinsky model, Dynamic Smagorinsky model, WALE model and subgrid turbulent kinetic energy model, are validated and compared directly against the experimental data. Two different mesh counts are used for the LES studies, one with a higher mesh resolution in the shear layer than the other. It is found that increasing the mesh resolution improves the time-averaged fluctuating velocity profiles, but has less effect on the time-averaged filtered velocity profiles. A comparison against experiment shows that the recirculation zone length is well predicted using LES. The accuracy of the four different subgrid scale models is then assessed by comparing the LES results using the dense mesh with the experiment. Comparisons with the time-averaged axial and radial velocity profiles demonstrate that LES displays good agreement with the experimental data, with the essential flow features captured both qualitative and quantitatively. The subgrid velocity also matches well with the experimental results, but a slight underprediction of the inner shear layer is observed for all subgrid models. In general, it is found that the Smagorinsky and WALE models are more dissipative than the Dynamic Smagorinsky model and subgrid TKE model. Comparison of the spectra against the experiment shows that LES can capture dominant features of the turbulent flow with reasonable accuracy, and weak spectral peaks related to the Kevin-Helmholtz instability and helical vortex shedding are present.     

Large-Eddy Simulation of Lean Premixed Turbulent Flames of Three Different Combustion Configurations using a Novel Reaction Closure

This large eddy simulation (LES) study is applied to three different premixed turbulent flames under lean conditions at atmospheric pressure. The hierarchy of complexity of these flames in ascending order are a simple Bunsen-like burner, a sudden-expansion dump combustor, and a typical swirl-stabilized gas turbine burner–combustor. The purpose of this paper is to examine numerically whether the chosen combination of the Smagorinsky turbulence model for sgs fluxes and a novel turbulent premixed reaction closure is applicable over all the three combustion configurations with varied degree of flow and turbulence. A quality assessment method for the LES calculations is applied. The cold flow data obtained with the Smagorinsky closure on the dump combustor are in close proximity with the experiments. It moderately predicts the vortex breakdown and bubble shape, which control the flame position on the double-cone burner. Here, the jet break-up at the root of the burner is premature and differs with the experiments by as much as half the burner exit diameter, attributing the discrepancy to poor grid resolution. With the first two combustion configurations, the applied subgrid reaction model is in good correspondence with the experiments. For the third case, a complex swirl-stabilized burner–combustor configuration, although the flow field inside the burner is only modestly numerically explored, the level of flame stabilization at the junction of the burner–combustor has been rather well captured. Furthermore, the critical flame drift from the combustor into the burner was possible to capture in the LES context (which was not possible with the RANS plus k–ɛ model), however, requiring tuning of a prefactor in the reaction closure.     

Dynamic Model for LES Without Test Filtering: Quantifying the Accuracy of Taylor Series Approximations

The dynamic model for large-eddy simulation (LES) of turbulent flows requires test filtering the resolved velocity fields in order to determine model coefficients. However, test filtering is costly to perform in LES of complex geometry flows, especially on unstructured grids. The objective of this work is to develop and test an approximate but less costly dynamic procedure which does not require test filtering. The proposed method is based on Taylor series expansions of the resolved velocity fields. Accuracy is governed by the derivative schemes used in the calculation and the number of terms considered in the approximation to the test filtering operator. The expansion is developed up to fourth order, and results are tested a priori based on direct numerical simulation data of forced isotropic turbulence in the context of the dynamic Smagorinsky model. The tests compare the dynamic Smagorinsky coefficient obtained from filtering with those obtained from application of the Taylor series expansion. They show that the expansion up to second order provides a reasonable approximation to the true dynamic coefficient (with errors on the order of about 5% for c _s ²), but that including higher-order terms does not necessarily lead to improvements in the results due to inherent limitations in accurately evaluating high-order derivatives. A posteriori tests using the Taylor series approximation in LES of forced isotropic turbulence and channel flow confirm that the Taylor series approximation yields accurate results for the dynamic coefficient. Moreover, the simulations are stable and yield accurate resolved velocity statistics. Received 20 February 2001 and accepted 24 July 2001     

Subgrid Scale Modeling in Large-Eddy Simulation of Turbulent Combustion Using Premixed Flamelet Chemistry

Large-eddy simulation (LES) of turbulent combustion with premixed flamelets is investigated in this paper. The approach solves the filtered Navier–Stokes equations supplemented with two transport equations, one for the mixture fraction and another for a progress variable. The LES premixed flamelet approach is tested for two flows: a premixed preheated Bunsen flame and a partially premixed diffusion flame (Sandia Flame D). In the first case, we compare the LES with a direct numerical simulation (DNS). Four non-trivial models for the chemical source term are considered for the Bunsen flame: the standard presumed beta-pdf model, and three new propositions (simpler than the beta-pdf model): the filtered flamelet model, the shift-filter model and the shift-inversion model. A priori and a posteriori tests are performed for these subgrid reaction models. In the present preheated Bunsen flame, the filtered flamelet model gives the best results in a priori tests. The LES tests for the Bunsen flame are limited to a case in which the filter width is only slightly larger than the flame thickness. According to the a posteriori tests the three models (beta-pdf, filtered flamelet and shift-inversion) show more or less the same results as the trivial model, in which subgrid reaction effects are ignored, while the shift-filter model leads to worse results. Since LES needs to resolve the large turbulent eddies, the LES filter width is bounded by a maximum. For the present Bunsen flame this means that the filter width should be of the order of the flame thickness or smaller. In this regime, the effects of subgrid reaction and subgrid flame wrinkling turn out to be quite modest. The LES-results of the second case (Sandia Flame D) are compared to experimental data. Satisfactory agreement is obtained for the main species. Comparison is made between different eddy-viscosity models for the subgrid turbulence, and the Smagorinsky eddy-viscosity is found to give worse results than eddy-viscosities that are not dominated by the mean shear. Paper presented on the Eccomas Thematic Conference Computational Combustion 2007, submitted for a special issue of Flow, Turbulence and Combustion.     

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.     

On the Use of Embedded Meshes in the LES of External Flows

In the present study, the possibility of performing cheap and accurate LES computations using unstructured grids on complex geometries has been examined. The test case considered is a turbulent flow around a thin controlled diffusion airfoil at low speed, with angle of attack of 8° to the incoming free-stream flow. The flow velocity is set to satisfy the chord Reynolds number of 1.2 × 10⁵. Two subgrid models are tested; the classical Smagorinsky model and the WALE model. The predictions are found to be in good agreement with the experimental data both in terms of flow dynamics and pressure spectra. The results are found to be independent of the LES physical models tested herein, although this conclusion is issued with the caveat that a careful grid generation procedure using (2:3 and 3:4 type) refinements was required to reduce the jump size across the different grid resolutions.     

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.     

Large-eddy simulation of circular cylinder flow at subcritical Reynolds number: Turbulent wake and sound radiation

The flows past a circular cylinder at Reynolds number 3900 are simulated using large-eddy simulation(LES) and the far-field sound is calculated from the LES results. A low dissipation energy-conserving finite volume scheme is used to discretize the incompressible Navier–Stokes equations. The dynamic global coefficient version of the Vreman's subgrid scale(SGS) model is used to compute the sub-grid stresses. Curle's integral of Lighthill's acoustic analogy is used to extract the sound radiated from the cylinder. The profiles of mean velocity and turbulent fluctuations obtained are consistent with the previous experimental and computational results. The sound radiation at far field exhibits the characteristic of a dipole and directivity. The sound spectra display the-5/3 power law. It is shown that Vreman's SGS model in company with dynamic procedure is suitable for LES of turbulence generated noise.     

Grid Resolution Effects on VSFMDF/LES

The grid dependence of LES/VSFMDF is studied on a series of grids with progressively increased resolution reaching over 10 million grids for simulation of a turbulent piloted nonpremixed methane jet flame (Sandia D). In VSFMDF, the effects of the subgrid scale chemical reaction and convection appear in closed forms. The modeled transport equation for the VSFMDF is solved by a hybrid finite-difference/Monte Carlo scheme. A flamelet model is employed to relate the instantaneous composition to the mixture fraction. The simulated results are assessed via comparison with laboratory data. In addition, the dependence of predicted statistics on the grid size of the simulation is studied. The first order moments converge for the finest grid, but the higher order statistics including the PDFs are more sensitive to the grid resolution.     

Joint RANS/LES Approach to Premixed Flame Modelling in the Context of the TFC Combustion Model

We present an original timesaving joint RANS/LES approach to simulate turbulent premixed combustion. It is intended mainly for industrial applications where LES may not be practical. It is based on successive RANS/LES numerical modelling, where turbulent characteristics determined from RANS simulations are used in LES equations for estimation of the subgrid chemical source and viscosity. This approach has been developed using our TFC premixed combustion model, which is based on a generalization of the Kolmogorov’s ideas. We assume existence of small-scale statistically equilibrium structures not only of turbulence but also of the reaction zones. At the same time, non-equilibrium large-scale structures of reaction sheets and turbulent eddies are described statistically by model combustion and turbulence equations in RANS simulations or follow directly without modelling in LES. Assumption of small-scale equilibrium gives an opportunity to express the mean combustion rate (controlled by small-scale coupling of turbulence and chemistry) in the RANS and LES sub-problems in terms of integral or subgrid parameters of turbulence and the chemical time, i.e. the definition of the reaction rate is similar to that of the mean dissipation rate in turbulence models where it is expressed in terms of integral or subgrid turbulent parameters. Our approach therefore renders compatible the combustion and turbulent parts of the RANS and LES sub-problems and yields reasonable agreement between the RANS and averaged LES results. Combining RANS simulations of averaged fields with LES method (and especially coupled and acoustic codes) for simulation of corresponding nonstationary process (and unsteady combustion regimes) is a promising strategy for industrial applications. In this work we present results of simulations carried out employing the joint RANS/LES approach for three examples: High velocity premixed combustion in a channel, combustion in the shear flow behind an obstacle and the impinging flame (a premixed flame attached to an obstacle).     

An unsteady incompressible Navier–Stokes solver for large eddy simulation of turbulent flows

An unsteady incompressible Navier–Stokes solver that uses a dual time stepping method combined with spatially high‐order‐accurate finite differences, is developed for large eddy simulation (LES) of turbulent flows. The present solver uses a primitive variable formulation that is based on the artificial compressibility method and various convergence–acceleration techniques are incorporated to efficiently simulate unsteady flows. A localized dynamic subgrid model, which is formulated using the subgrid kinetic energy, is employed for subgrid turbulence modeling. To evaluate the accuracy and the efficiency of the new solver, a posteriori tests for various turbulent flows are carried out and the resulting turbulence statistics are compared with existing experimental and direct numerical simulation (DNS) data. Copyright © 1999 John Wiley & Sons, Ltd.     

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.     

Influence of mesh deformation on the quality of large eddy simulations

The influence of mesh motion on the quality of large eddy simulation (LES) was studied in the present article. A three‐dimensional, turbulent pipe flow (Re_τ=360) was considered as a test case. Simulations with both stretching and static meshes were carried out in order to understand how mesh motion affects the turbulence statistics. The spatial filtering of static and moving mesh direct numerical simulation (DNS) data showed how an ideal LES would perform, while the comparison of DNS cases with static and moving meshes revealed that no significant numerical errors arise from the mesh motion when the simulation is fully resolved. The comparison of the filtered fields of the DNS with a moving mesh with the corresponding LES fields revealed different responses to mesh motion from different numerical approaches. A straightforward test was applied in order to verify that the moving mesh works consistently in LES: when the mesh is stretched in the streamwise direction, the moving mesh results should be in between the two extremal resolutions between which the mesh is stretched. Numerical investigations using four different LES approaches were carried out. In addition to the Smagorinsky model, three implicit LES approaches were used: linear interpolation (non‐dissipative), the Gamma limiter (dissipative), and the scale‐selective discretisation (slightly dissipative). The results indicate that while the Smagorinsky and the scale‐selective discretisation approaches produce results consistent with the resolution of the non‐static mesh, the implicit LES with linear interpolation or the Gamma scheme do not. Copyright © 2016 John Wiley & Sons, Ltd.     

Coupling tabulated chemistry with Large Eddy Simulation of turbulent reactive flows

A new modeling strategy is developed to introduce tabulated chemistry methods in the LES of turbulent premixed combustion. The objective is to recover the correct laminar flame propagation speed of the filtered flame front when the subgrid scale turbulence vanishes. The filtered flame structure is mapped by 1D filtered laminar premixed flames. Closure of the filtered progress variable and the energy balance equations are carefully addressed. The methodology is applied to 1D and 2D filtered laminar flames. These computations show the capability of the model to recover the laminar flame speed and the correct chemical structure when the flame wrinkling is completely resolved. The model is then extended to turbulent combustion regimes by introducing subgrid scale wrinkling effects on the flame front propagation. Finally, the LES of a 3D turbulent premixed flame is performed. To cite this article: R. Vicquelin et al., C. R. Mecanique 337 (2009).     

A-priori testing of an eddy viscosity model for the density-weighted subgrid scale stress tensor in turbulent premixed flames

In this study, we report on the direct measurement of the density-weighted subgrid scale (SGS) stress tensor in turbulent premixed flames. In large-eddy simulations (LES), this unresolved tensor is typically modelled using eddy viscosity approaches. Additionally to the direct measurement, we provide a pure experimentally based a-priori test of the commonly used eddy viscosity model suggested by Smagorinsky. For two turbulent premixed V-shaped methane–air flames, a statistical analysis is presented where the correlation between the directly measured SGS stress tensor and the eddy viscosity model following Smagorinsky is tested. The measurement strategy is based on the application of a dual-plane stereo-PIV technique which enables the measurement of the 3D flow field in two parallel planes. This allows the determination of velocities as well as velocity gradients in all three directions. Here, a vector resolution of 118 μm was achieved. For a priori testing, the data are subjected to a spatial filtering procedure that reproduces the application of the filter function in LES. The calculation of velocity gradients is performed after the application of this spatial averaging. Additionally to the velocity field, the flame front position is deduced from the clearly observable step in the tracer particle number density between burnt and unburnt regions of the flame. This facilitates the direct single-shot-based evaluation of all components of the density-weighted SGS stress tensor. Additionally, the model expressions related to these terms can be determined, which is done in this first study for the static Smagorinsky model. With that, the instantaneous local comparison between directly measured stress terms and modelled terms is possible, based on the instantaneous local evaluation procedure. The measurement procedure is described, and first results are presented and discussed. They show a rather poor performance of the static form of the Smagorinsky model (with fixed Smagorinsky constant). Our future aims are to use the directly measured SGS data for the a-priori comparison with more advanced models.     

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.     

Fourier neural operator approach to large eddy simulation of three-dimensional turbulence

Fourier neural operator (FNO) model is developed for large eddy simulation (LES) of three-dimensional (3D) turbulence. Velocity fields of isotropic turbulence generated by direct numerical simulation (DNS) are used for training the FNO model to predict the filtered velocity field at a given time. The input of the FNO model is the filtered velocity fields at the previous several time-nodes with large time lag. In the a posteriori study of LES, the FNO model performs better than the dynamic Smagorinsky model (DSM) and the dynamic mixed model (DMM) in the prediction of the velocity spectrum, probability density functions (PDFs) of vorticity and velocity increments, and the instantaneous flow structures. Moreover, the proposed model can significantly reduce the computational cost, and can be well generalized to LES of turbulence at higher Taylor-Reynolds numbers.     

A computational study of supersonic combustion behind a wedge-shaped flameholder

In this study, large eddy simulation (LES) has been used to examine supersonic flow, mixing, self-ignition and combustion in a model scramjet combustor and has been compared against the experimental data. The LES model is based on an unstructured finite-volume discretization, using monotonicity-preserving flux reconstruction of the filtered mass, momentum, species and energy equations. Both a two-step and a seven-step hydrogen–air mechanism are used to describe the chemical reactions. Additional comparisons are made with results from a previously presented flamelet model. The subgrid flow terms are modeled using a mixed model, whereas the subgrid turbulence–chemistry interaction terms are modeled using the partially stirred reactor model. Simulations are carried out on a scramjet model experimentally studied at Deutsches Zentrum für Luft- und Raumfahrt consisting of a one-sided divergent channel with a wedge-shaped flame holder at the base of which hydrogen is injected. The LES predictions are compared with experimental data for velocity, temperature, wall pressure at different cross sections as well as schlieren images, showing good agreement for both first- and second-order statistics. In addition, the LES results are used to illustrate and explain the intrinsic flow, and mixing and combustion features of this combustor.     

A hybrid dynamic Smagorinsky model for large eddy simulation

A hybrid dynamic subgrid-scale model (HDSM) pertaining to Large-eddy simulation (LES) has been developed. The coefficient obtained by German dynamic Smagorinsky model (DSM) was integrated with a new dynamic coefficient, based on the dynamic subgrid characteristic length and controlled by the subgrid-scale (SGS) motions. In HDSM, the characteristic wave number determining the characteristic length of the dynamic subgrid is calculated from a new energy weighted mean method when the subgrid scale turbulent kinetic energy and the dissipation wave number are known. The dissipation wave-number is derived from the SGS turbulent kinetic energy spectrum equation. The total dissipation rate spectrum equation is based on the Pao energy spectrum and local equilibrium assumption. The dynamic subgrid characteristic length could take into account the rapidly fluctuating small scale behaviours and the spatial variation of turbulent characteristics. HDSM was used to simulate the fully developed channel and turbulent flow past a circular cylinder, and to determine the impact of the dam-break flow on downstream structure. The HDSM is robust in respect to anisotropic mesh and is less sensitive to grid resolution, and would accurately describe the energy transfer from large-scale to SGS fluctuations and capture more fluctuations of turbulence with same meshes compared to the DSM.