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.     

Multi-variable special polynomials using an operator ordering method

Using an operator ordering method for some commutative superposition operators, we introduce two new multi-variable special polynomials and their generating functions, and present some new operator identities and integral formulas involving the two special polynomials. Instead of calculating complicated partial differential, we use the special polynomials and their generating functions to concisely address the normalization, photocount distributions and Wigner distributions of several quantum states that can be realized physically, the results of which provide real convenience for further investigating the properties and applications of these states.     

Flow-induced noise analysis for 3D trash rack based on LES/Lighthill hybrid method

A trash rack is used to obstruct dirt in a flow. As it is on the outside of the hull of a ship, such flow-induced noise cannot be ignored. This paper focuses on a hydrodynamic noise study of the trash rack installed on the inlet of pipelines to estimate the noise produced by the turbulent flow. A hybrid method of combining Large-Eddy Simulation (LES) and the Lighthill’s acoustic analogy theory simulate the flow-induced noise for a three-dimensional pipeline. The results of the simulation agree well with the experimental results available in the open literature. Following this, the hydrodynamic noise of flow through a three-dimensional pipeline with different types of trash racks is studied in order to determine the best form. For this objective, the flow and sound properties in different simulation models are analyzed and a low noise type of trash rack is proposed based on some general indexes.     

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.     

Assessment of LES of intermittent soot production in an aero-engine model combustor using high-speed measurements

Soot production in turbulent flames is an extremely intermittent phenomenon since it is the result of specific thermochemical conditions occasionally occurring in space and time. In realistic configurations such as the swirling flames used in gas-turbines, the presence of large-scale flow motions can additionally affect soot formation processes, leading to even more pronounced intermittency. Classically, the validation of numerical simulations is performed by comparing time-averaged results with experimental data of the phenomenon under investigation. This comparison can be considered as rigorous only if a statistically converged numerical representation is obtained. In case of sporadic events such as intermittent soot formation in turbulent flames, this means to perform the simulation over thousands of milliseconds of physical time, which can have extremely high CPU demands when performing Large Eddy Simulation (LES). In this work, a possible strategy to overcome this issue is proposed based on the use of high-speed measurements and numerically synthesized signals from LES. To illustrate the approach, numerical and experimental soot light scattering signals are considered here by looking at the model aero-engine combustor developed at DLR for the study of pressurized swirled sooting flames. The light scattering signal is numerically synthesized from an LES. Experimental high-speed measurements are used to statistically account for the high temporal and spatial variability of soot when considering time intervals similar to what is today achievable with LES. The feasibility of this approach is finally demonstrated by comparing numerical results to the ensemble of possible soot production states observed experimentally in the DLR burner allowing to eventually validate the present LES results.     

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.     

A basic hybrid finite element formulation for mid-frequency analysis of beams connected at an arbitrary angle

When beams are connected at an arbitrary angle and subjected to an external excitation, both longitudinal and bending waves are generated in the system. Since longitudinal wavelengths are considerably longer than bending wavelengths in the mid-frequency region, the number of bending wavelengths in the beams is considerably larger than the number of longitudinal wavelengths. In this paper, plannar beams connected at arbitrary angles are considered. The energy finite element analysis (EFEA) is employed for modelling the bending behavior of the beams and the conventional finite element analysis (FEA) is utilized for modelling the longitudinal vibration in the beams. Thus, a basic hybrid FEA formulation is presented for mid-frequency analysis of systems that contain two types of energy. The bending vibration is associated with the long members in the system and the longitudinal vibration is associated with the short members. The long members are considered to have high modal overlap and to contain several wavelengths within their dimension, and uncertainty effects are present. The short members contain a small number of wavelengths, and exhibit a low modal overlap. Due to the low modal overlap the resonant frequencies are spaced far apart in the frequency domain, therefore the short members exhibit resonant or non-resonant behavior depending on the frequency of the excitation.In this work, the bending and the longitudinal vibration within the same beam member are treated as a long and as a short member, respectively. A hybrid joint formulation is developed between long and short members. Power reflection and transmission coefficients are derived for each joint. The distribution of the energy throughout the system demonstrates a strong dependency on the power transfer coefficients. Several systems are analyzed by the hybrid FEA and by analytical solutions, and good correlation between them is observed.     

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.