A consistent dual-mesh framework for hybrid LES/RANS modeling

In this work, we propose a consistent turbulence-modeling framework for hybrid LES/RANS modeling. In this framework, the filtered and Reynolds averaged Navier–Stokes (RANS) equations are solved simultaneously in the whole domain on their respective meshes. Consistency between the two solutions is achieved in terms of velocity, pressure, and turbulent quantities through additional drift terms in the corresponding equations. This approach leads to clean conditions at the LES/RANS interfaces. Note that this general framework does not depend on the specific choice of LES and RANS models. A hybrid LES/RANS solver is developed within this framework and used to simulate the flow in a plane channel and that in a channel with periodic hills. The results demonstrate that the hybrid solver leads to significantly improved results with moderate computational overhead compared to traditional LES, making it a promising candidate for industrial flow simulations.     

用RANS/LES混合方法研究超声速底部流动

在剪切应力输运和弱非线性k-ω模式中引入压缩性修正,并在此基础上构造RANS/LES混合方法,对包含丰富流动结构和复杂流动机理的导弹类超声速底部流场进行数值分析,发现压缩性修正可更好地描述流场特征,适于具有较强压缩性效应流动的数值分析;混合方法捕捉到丰富的流动现象及非定常特性.     

An algorithm for LES of premixed compressible flows using the Conditional Moment Closure model

A novel numerical method has been developed to couple a recent high order accurate fully compressible upwind method with the Conditional Moment Closure combustion model. The governing equations, turbulence modelling and numerical methods are presented in full. The new numerical method is validated against direct numerical simulation (DNS) data for a lean premixed methane slot burner. Although the modelling approaches are based on non-premixed flames and hence not expected to be valid for a wide range of premixed flames, the predicted flame is just 10% longer than that in the DNS and excellent agreement of mean mass fractions, conditional mass fractions and temperature is demonstrated. This new numerical method provides a very useful framework for future application of CMC to premixed as well as non-premixed combustion.     

The boson-fermion hybrid representation formulation,I

A boson-fermion hybrid representation is presented. In this framework, a fermion system is described concurrently by the bosonic and the fermionic degrees of freedom. A fermion pair in this representation can be treated as a boson without violating the Pauli principle. Furthermore the “bosonic interactions” are shown to originate from the exchange processes of the fermions and can be calculated from the original fermion interactions. Both the formulation of the BFH representations for the even and odd nuclear systems are given. We find that the basic equation of the nuclear field theory (NFT) is just the usual Schrödinger equation in such a representation with the empirical NFT diagrammatic rules emerging naturally. This theory was numerically checked in the case of four nucleons moving in a single-j shell and the exactness of the theory was established.     

On the simulation of trailing edge noise with a hybrid LES/APE method

A hybrid method is applied to predict trailing edge noise based on a large eddy simulation (LES) of the compressible flow problem and acoustic perturbation equations (APE) for the time-dependent simulation of the acoustic field. The acoustic simulation in general considers the mean flow convection and refraction effects such that the computational domain of the flow simulation has to comprise only the significant acoustic source region. Using a modified rescaling method for the prediction of the unsteady turbulent inflow boundary layer, the LES just resolves the flow field in the immediate vicinity of the trailing edge. The linearized APE completely prevent the unbounded growth of hydrodynamic instabilities in critical mean flows.     

A consistent LES/filtered-density function formulation for the simulation of turbulent flames with detailed chemistry

A hybrid large-Eddy simulation/filtered-density function (LES–FDF) methodology is formulated for simulating variable density turbulent reactive flows. An indirect feedback mechanism coupled with a consistency measure based on redundant density fields contained in the different solvers is used to construct a robust algorithm. Using this novel scheme, a partially premixed methane/air flame is simulated. To describe transport in composition space, a 16-species reduced chemistry mechanism is used along with the interaction-by-exchange with the mean (IEM) model. For the micro-mixing model, typically a constant ratio of scalar to mechanical time-scale is assumed. This parameter can have substantial variations and can strongly influence the combustion process. Here, a dynamic time-scale model is used to prescribe the mixing time-scale, which eliminates the time-scale ratio as a model constant. Two different flame configurations, namely, Sandia flames D and E are studied. Comparison of simulated radial profiles with experimental data show good agreement for both flames. The LES–FDF simulations accurately predict the increased extinction near the inlet and re-ignition further downstream. The conditional mean profiles show good agreement with experimental data for both flames.     

A hybrid finite/boundary element method for periodic structures on non-periodic meshes using an interior penalty formulation for Maxwell’s equations

This paper presents a hybrid finite element/boundary element (FEBE) method for periodic structures. Periodic structures have been efficiently analyzed by solving for a single unit cell utilizing Floquet’s theorem. However, most of the previous works require periodic meshes to properly impose the boundary conditions on the outer surfaces of the unit cell. To alleviate this restriction, the interior penalty method is adopted and implemented in this work. Also, the proper treatment of the boundary element part is addressed to account for the non-conformity of the boundary element mesh. Another ingredient of this work is the use of the efficient boundary element computation, accelerated by the Ewald transformation for the calculation of the periodic Green’s function. Finally, the method is validated through examples which are discretized without the constraint of a periodic mesh.     

Gauge covariant formulation of the generating operator. 1. The Zakharov-Shabat system

A gauge covariant formulation of the generating operator Λ, related to the Zakharov-Shabat system L is proposed. The operator Λ̃, corresponding to $\text{L}\text{?}\text{= ψ}{}_0{}^{\mathrm{?1}}\text{Lψ}{}_0$ in the pole gauge is explicitly calculated. Thus the unified approach to the higher nonlinear Schrödinger equations, based on Λ is automatically reformulated with Λ̃ for the higher Heisenberg ferromagnet equations.     

Evaluation of the flame propagation within an SI engine using flame imaging and LES

This work shows experiments and simulations of the fired operation of a spark ignition engine with port-fuelled injection. The test rig considered is an optically accessible single cylinder engine specifically designed at TU Darmstadt for the detailed investigation of in-cylinder processes and model validation. The engine was operated under lean conditions using iso-octane as a substitute for gasoline. Experiments have been conducted to provide a sound database of the combustion process. A planar flame imaging technique has been applied within the swirl- and tumble-planes to provide statistical information on the combustion process to complement a pressure-based comparison between simulation and experiments. This data is then analysed and used to assess the large eddy simulation performed within this work. For the simulation, the engine code KIVA has been extended by the dynamically thickened flame model combined with chemistry reduction by means of pressure dependent tabulation. Sixty cycles have been simulated to perform a statistical evaluation. Based on a detailed comparison with the experimental data, a systematic study has been conducted to obtain insight into the most crucial modelling uncertainties.     

Correspondence between the classical and quantum canonical transformation groups from an operator formulation of the wigner function

An explicit expression of the Wigner operator is derived, such that the Wigner function of a quantum state is equal to the expectation value of this operator with respect to the same state. This Wigner operator leads to a representation-independent procedure for establishing the correspondence between the inhomogeneous symplectic group applicable to linear canonical transformations in classical mechanics and the Weyl-metaplectic group governing the symmetry of unitary transformations in quantum mechanics.     

A high-order accurate hybrid scheme using a central flux scheme and a WENO scheme for compressible flowfield analysis

A high-order accurate hybrid central-WENO scheme is proposed. The fifth order WENO scheme [G.S. Jiang, C.W. Shu, Efficient implementation of weighted ENO schemes, J. Comput. Phys. 126 (1996) 202–228] is divided into two parts, a central flux part and a numerical dissipation part, and is coupled with a central flux scheme. Two sub-schemes, the WENO scheme and the central flux scheme, are hybridized by means of a weighting function that indicates the local smoothness of the flowfields. The derived hybrid central-WENO scheme is written as a combination of the central flux scheme and the numerical dissipation of the fifth order WENO scheme, which is controlled adaptively by a weighting function. The structure of the proposed hybrid central-WENO scheme is similar to that of the YSD-type filter scheme [H.C. Yee, N.D. Sandham, M.J. Djomehri, Low-dissipative high-order shock-capturing methods using characteristic-based filters, J. Comput. Phys. 150 (1999) 199–238]. Therefore, the proposed hybrid scheme has also certain merits that the YSD-type filter scheme has. The accuracy and efficiency of the developed hybrid central-WENO scheme are investigated through numerical experiments on inviscid and viscous problems. Numerical results show that the proposed hybrid central-WENO scheme can resolve flow features extremely well.     

Development of an algebraic wall heat transfer model for LES in IC engines using DNS data

This paper presents a novel algebraic wall-modeled large-eddy simulation (WMLES) approach for wall heat transfer (WHT) in internal combustion engines (ICE) using recent high-fidelity simulation data from two real ICE and a flame-wall interaction (FWI) setup. The model formulation is based on intuitive arguments rather than simplified forms of near-wall governing equations. Input information from the two wall-normal wall-adjacent nodes is used, facilitating the interpretation of the local state of the thermal boundary layer (TBL). With filtered direct numerical simulation (DNS) data (i.e. a ‘perfect’ LES), model performance is evaluated locally at different near-wall resolutions. The proposed model is compared to a widely used wall function (WF) approach and to a data-driven model.     

Studies of the flow and turbulence fields in a turbulent pulsed jet flame using LES/PDF

A turbulent piloted jet flame subject to a rapid velocity pulse in its fuel jet inflow is proposed as a new benchmark case for the study of turbulent combustion models. In this work, we perform modelling studies of this turbulent pulsed jet flame and focus on the predictions of its flow and turbulence fields. An advanced modelling strategy combining the large eddy simulation (LES) and the probability density function (PDF) methods is employed to model the turbulent pulsed jet flame. Characteristics of the velocity measurements are analysed to produce a time-dependent inflow condition that can be fed into the simulations. The effect of the uncertainty in the inflow turbulence intensity is investigated and is found to be very small. A method of specifying the inflow turbulence boundary condition for the simulations of the pulsed jet flame is assessed. The strategies for validating LES of statistically transient flames are discussed, and a new framework is developed consisting of different averaging strategies and a bootstrap method for constructing confidence intervals. Parametric studies are performed to examine the sensitivity of the predictions of the flow and turbulence fields to model and numerical parameters. A direct comparison of the predicted and measured time series of the axial velocity demonstrates a satisfactory prediction of the flow and turbulence fields of the pulsed jet flame by the employed modelling methods.     

The Schrödinger equation via an operator functional equation

In this paper we derive the Schrödinger equation by comparing quantum statistics with classical statistical mechanics, identifying similarities and differences, and developing an operator functional equation which is solved in a completely algebraic fashion with no appeal to spatial invariances or symmetries.     

The effect of fuel inlet turbulence intensity on H2/CH4 flame structure of MILD combustion using the LES method

Large eddy simulations (LES) are employed to investigate the effect of the inlet turbulence intensity on the H₂/CH₄ flame structure in a hot and diluted co-flow stream which emulates the (Moderate or Intense Low-oxygen Dilution) MILD combustion regime. In this regard, three fuel inlet turbulence intensity profiles with the values of 4%, 7% and 10% are superimposed on the annular mixing layer. The effects of these changes on the flame structure under the MILD condition are studied for two oxygen concentrations of 3% and 9% (by mass) in the oxidiser stream and three hot co-flow temperatures 1300, 1500 and 1750 K. The turbulence-chemistry interaction of the numerically unresolved scales is modelled using the (Partially Stirred Reactor) PaSR method, where the full mechanism of GRI-2.11 represents the chemical reactions. The influences of the turbulence intensity on the flame structure under the MILD condition are studied by using the profile of temperature, CO and OH mass fractions in both physical and mixture fraction spaces at two downstream locations. Also, the effects of this parameter are investigated by contours of OH, HCO and CH₂O radicals in an area near the nozzle exit zone. Results show that increasing the fuel inlet turbulence intensity has a profound effect on the flame structure particularly at low oxygen mass fraction. This increment weakens the combustion zone and results in a decrease in the peak values of the flame temperature and OH and CO mass fractions. Furthermore, increasing the inlet turbulence intensity decreases the flame thickness, and increases the MILD flame instability and diffusion of un-burnt fuel through the flame front. These effects are reduced by increasing the hot co-flow temperature which reinforces the reaction zone.     

A class of hybrid DG/FV methods for conservation laws I: Basic formulation and one-dimensional systems

By comparing the discontinuous Galerkin (DG) and the finite volume (FV) methods, a concept of ‘static reconstruction’ and ‘dynamic reconstruction’ is introduced for high-order numerical methods. Based on the new concept, a class of hybrid DG/FV schemes is presented for one-dimensional conservation law using a ‘hybrid reconstruction’ approach. In the hybrid DG/FV schemes, the lower-order derivatives of a piecewise polynomial solution are computed locally in a cell by the DG method based on Taylor basis functions (called as ‘dynamic reconstruction’), while the higher-order derivatives are re-constructed by the ‘static reconstruction’ of the FV method, using the known lower-order derivatives in the cell itself and its adjacent neighboring cells. The hybrid DG/FV methods can greatly reduce CPU time and memory required by the traditional DG methods with the same order of accuracy on the same mesh, and they can be extended directly to unstructured and hybrid grids in two and three dimensions similar to the DG and/or FV methods. The hybrid DG/FV methods are applied to one-dimensional conservation law, including linear and non-linear scalar equation and Euler equations. In order to capture the strong shock waves without spurious oscillations, a simple shock detection approach is developed to mark ‘trouble cells’, and a moment limiter is adopted for higher-order schemes. The numerical results demonstrate the accuracy, and the super-convergence property is shown for the third-order hybrid DG/FV schemes. In addition, by analyzing the eigenvalues of the semi-discretized system in one dimension, we discuss the spectral properties of the hybrid DG/FV schemes to explain the super-convergence phenomenon.     

A conservative coupling algorithm between a compressible flow and a rigid body using an Embedded Boundary method

This paper deals with a new solid–fluid coupling algorithm between a rigid body and an unsteady compressible fluid flow, using an Embedded Boundary method. The coupling with a rigid body is a first step towards the coupling with a Discrete Element method. The flow is computed using a finite volume approach on a Cartesian grid. The expression of numerical fluxes does not affect the general coupling algorithm and we use a one-step high-order scheme proposed by Daru and Tenaud [V. Daru, C. Tenaud, J. Comput. Phys. (2004)]. The Embedded Boundary method is used to integrate the presence of a solid boundary in the fluid. The coupling algorithm is totally explicit and ensures exact mass conservation and a balance of momentum and energy between the fluid and the solid. It is shown that the scheme preserves uniform movement of both fluid and solid and introduces no numerical boundary roughness. The efficiency of the method is demonstrated on challenging one- and two-dimensional benchmarks.