Large eddy simulation of piloted pulverised coal combustion using extended flamelet/progress variable model

An extended flamelet/progress variable (EFPV) model for simulating pulverised coal combustion (PCC) in the context of large eddy simulation (LES) is proposed, in which devolatilisation, char surface reaction and radiation are all taken into account. The pulverised coal particles are tracked in the Lagrangian framework with various sub-models and the sub-grid scale (SGS) effects of turbulent velocity and scalar fluctuations on the coal particles are modelled by the velocity-scalar joint filtered density function (VSJFDF) model. The presented model is then evaluated by LES of an experimental piloted coal jet flame and comparing the numerical results with the experimental data and the results from the eddy break up (EBU) model. Detailed quantitative comparisons are carried out. It is found that the proposed model performs much better than the EBU model on radial velocity and species concentrations predictions. Comparing against the adiabatic counterpart, we find that the predicted temperature is evidently lowered and agrees well with the experimental data if the conditional sampling method is adopted.     

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.     

An analysis of subgrid-resolved scale interactions with use of results from direct numerical simulations

Subgrid nonlinear interaction and energy transfer are analyzed using direct numerical simulations of isotropic turbulence. Influences of cutoff wave number at different ranges of scale on the energetics and dynamics have been investigated. It is observed that subgrid-subgrid interaction dominates the turbulent dynamics when cut-off wave number locates in the energy-containing range while resolved-subgrid interaction dominates if it is in the dissipation range. By decomposing the subgrid energy transfer and nonlinear interaction into ‘forward’ and ‘backward’ groups according to the sign of triadic interaction, we find that individually each group has very large contribution, but the net of them is much smaller, implying that tremendous cancellation happens between these two groups.     

高亚声速湍流喷流气动噪声数值分析

为适应航空噪声管制规定要求,发动机喷流噪声控制成为目前气动声学研究中的重要课题,预测分析喷流噪声辐射并揭示其产生机理将为噪声控制奠定基础.采用高精度并行LES（large eddy simulation）方法计算分析马赫数0.9高亚声速喷流的湍流演化和气动噪声现象.首先,仔细验证喷流LES湍流场计算保真性,并分析流场中不同尺度涡结构的演化形态.其次,利用可穿透面FW-H（Ffowcs Williams and Hawkings）方法外推喷流近场声源数据获得精确声辐射远场,进而分析声场主导声模态特性.最后,通过分析声源机制、分离声模态等方法研究势流核末端大尺度拟序涡运动演化形成的低波数波包在噪声主导声模态产生中的重要作用.数值结果表明LES结合可穿透面FW-H方法可精确预测高亚声速喷流的流场及声场特征,且数值分析揭示涡环对并形成的大尺度拟序结构在喷流中心线上沿径向融合,产生了在远场低方位角占优的主导声模态,并构成强指向性声场,噪声峰值方位角约为30°.      

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.     

燃烧室两相流场亚网格燃烧模型的研究

在三维任意曲线坐标系下采用不同的亚网格燃烧模型对环形燃烧室火焰筒气液两相湍流瞬态反应流进行大涡模拟．计算中所采用的数学度模型有：k方程亚网格尺度模型估算亚网格湍流黏性;热通量辐射模型估算辐射换热,分别采用亚网格EBU燃烧模型（E-A model）、亚网格二阶矩输运方程模型（SOM）和亚网格二阶矩代数模型（SOM-A）估算化学反应速率．并在非交错网格系统下气相采用SIMPLE算法和混合差分格式求解,液相采用Lagrange处理,并用PSIC算法对其进行求解．通过实验结果和计算结果的比较,表明在三维任意曲线坐标系下对燃烧室火焰简两相湍流油雾燃烧流场进行大涡模拟,3种不同的亚网格燃烧模型都能真实反映两相湍流化学反应流流动及实际燃烧过程,而采用亚网格二阶矩输运方程模型稍优于其他两种亚网格燃烧模型．     

An investigation of thermally stratified turbulent channel flow with temperature oscillation on the bottom wall by large eddy simulation

Thermally stratified shear turbulent channel flow with temperature oscillation on the bottom wall of the channel is calculated to investigate the behavior of turbulent flow and heat transfer by use of large eddy simulation (LES) approach coupled with dynamic subgrid-scale (SGS) models. The objective of this study is to deal with the effect of the temperature oscillation on turbulent behavior of thermally stratified turbulent channel flow and to examine the effectiveness of the LES technique for predicting statistically unsteady turbulent flow driven by time-varying buoyancy force. To validate the present calculation, thermally stratified shear turbulent channel flow is computed and compared with available data obtained by direct numerical simulation (DNS), which confirm that the present approach can be used to predict thermally stratified turbulent channel flow satisfactorily. Further, to illustrate the effect of the temperature oscillation with different Richardson numbers and periods of the oscillation on turbulence characteristics, the phase-averaged mean value and fluctuation of the resolved velocities and temperature, and instantaneous velocity fluctuation structures are analyzed.     

Investigation of the sensitivity of turbulent closures and coupling of hybrid RANS‐LES models for predicting flow fields with separation and reattachment

Hybrid models have found widespread applications for simulation of wall‐bounded flows at high Reynolds numbers. Typically, these models employ Reynolds‐averaged Navier–Stokes (RANS) and large eddy simulation (LES) in the near‐body and off‐body regions, respectively. A number of coupling strategies between the RANS and LES regions have been proposed, tested, and applied in the literature with varying degree of success. Linear eddy‐viscosity models (LEVM) are often used for the closure of turbulent stress tensor in RANS and LES regions. LEVM incorrectly predicts the anisotropy of Reynolds normal stress at the RANS‐LES interface region. To overcome this issue, use of non‐linear eddy‐viscosity models (NLEVM) have started receiving attention. In this study, a generic non‐linear blended modeling framework for performing hybrid simulations is proposed. Flow over the periodic hills is used as the test case for model evaluation. This case is chosen due to complex flow physics with simplified geometry. Analysis of the simulations suggests that the non‐linear hybrid models show a better performance than linear hybrid models. It is also observed that the non‐linear closures are less sensitive to the RANS‐LES coupling and grid resolution. Copyright © 2016 John Wiley & Sons, Ltd.     

On stabilized space‐time FEM for anisotropic meshes: Incompressible Navier–Stokes equations and applications to blood flow in medical devices

In complex applications, such as the analysis of hydraulic performance of blood pumps (ventricular assist devices), the Navier–Stokes equations have to be discretized on very anisotropic meshes. If stabilized finite element formulations are applied, standard definitions of the stabilization parameter are usually not appropriate to handle elements with a high aspect ratio. If, in addition, rotating objects, moving meshes, or turbulence has to be considered in the simulation, further modifications of the stabilization procedure have to be applied. In this paper, we present stabilized space‐time finite element formulations of the incompressible Navier–Stokes equations that show very good convergence properties on complex anisotropic meshes and lead to reasonable numerical accuracy in complex flows when compared with experimental data. Copyright © 2017 John Wiley & Sons, Ltd.