Center for Modeling of Turbulence and Transition (CMOTT). Research briefs: 1990

Brief progress reports of the Center for Modeling of Turbulence and Transition (CMOTT) research staff from May 1990 to May 1991 are given. The objectives of the CMOTT are to develop, validate, and implement the models for turbulence and boundary layer transition in the practical engineering flows. The flows of interest are three dimensional, incompressible, and compressible flows with chemistry. The schemes being studied include the two-equation and algebraic Reynolds stress models, the full Reynolds stress (or second moment closure) models, the probability density function models, the Renormalization Group Theory (RNG) and Interaction Approximation (DIA), the Large Eddy Simulation (LES) and Direct Numerical Simulation (DNS)

Institute for Computational Mechanics in Propulsion (ICOMP)

The Institute for Computational Mechanics in Propulsion (ICOMP) is a combined activity of Case Western Reserve University, Ohio Aerospace Institute (OAI) and NASA Lewis. The purpose of ICOMP is to develop techniques to improve problem solving capabilities in all aspects of computational mechanics related to propulsion. The activities at ICOMP during 1991 are described

Hermite regularization of the Lattice Boltzmann Method for open source computational aeroacoustics

The lattice Boltzmann method (LBM) is emerging as a powerful engineering tool for aeroacoustic computations. However, the LBM has been shown to present accuracy and stability issues in the medium-low Mach number range, that is of interest for aeroacoustic applications. Several solutions have been proposed but often are too computationally expensive, do not retain the simplicity and the advantages typical of the LBM, or are not described well enough to be usable by the community due to proprietary software policies. We propose to use an original regularized collision operator, based on the expansion in Hermite polynomials, that greatly improves the accuracy and stability of the LBM without altering significantly its algorithm. The regularized LBM can be easily coupled with both non-reflective boundary conditions and a multi-level grid strategy, essential ingredients for aeroacoustic simulations. Excellent agreement was found between our approach and both experimental and numerical data on two different benchmarks: the laminar, unsteady flow past a 2D cylinder and the 3D turbulent jet. Finally, most of the aeroacoustic computations with LBM have been done with commercial softwares, while here the entire theoretical framework is implemented on top of an open source library (Palabos).Comment: 34 pages, 12 figures, The Journal of the Acoustical Society of America (in press

Large eddy simulation of highly turbulent under-expanded hydrogen and methane jets for gaseous-fuelled internal combustion engines

Burning hydrogen in conventional internal combustion (IC) engines is associated with zero carbon-based tailpipe exhaust emissions. In order to obtain high volumetric efficiency and eliminate abnormal combustion modes such as preignition and backfire, in-cylinder direct injection (DI) of hydrogen is considered preferable for a future generation of hydrogen IC engines. However, hydrogen's low density requires high injection pressures for fast hydrogen penetration and sufficient in-cylinder mixing. Such pressures lead to chocked flow conditions during the injection process which result in the formation of turbulent under-expanded hydrogen jets. In this context, fundamental understanding of the under-expansion process and turbulent mixing just after the nozzle exit is necessary for the successful design of an efficient hydrogen injection system and associated injection strategies. The current study used large eddy simulation (LES) to investigate the characteristics of hydrogen under-expanded jets with different nozzle pressure ratios (NPR), namely 8.5, 10, 30 and 70. A test case of methane injection with NPR = 8.5 was also simulated for direct comparison with the hydrogen jetting under the same NPR. The near-nozzle shock structure, the geometry of the Mach disk and reflected shock angle, as well as the turbulent shear layer were all captured in very good agreement with data available in the literature. Direct comparison between hydrogen and methane fuelling showed that the ratio of the specific heats had a noticeable effect on the near-nozzle shock structure and dimensions of the Mach disk. It was observed that with methane, mixing did not occur before the Mach disk, whereas with hydrogen high levels of momentum exchange and mixing appeared at the boundary of the intercepting shock. This was believed to be the effect of the high turbulence fluctuations at the nozzle exit of the hydrogen jet which triggered Gortler vortices. Generally, the primary mixing was observed to occur after the location of the Mach disk and particularly close to the jet boundaries where large-scale turbulence played a dominant role. It was also found that NPR had significant effect on the mixture's local fuel richness. Finally, it was noted that applying higher injection pressure did not essentially increase the penetration length of the hydrogen jets and that there could be an optimum NPR that would introduce more enhanced mixing whilst delivering sufficient fuel in less time. Such an optimum NPR could be in the region of 100 based on the geometry and observations of the current study

Computational fluid dynamics combustion analysis evaluation

This study involves the development of numerical modelling in spray combustion. These modelling efforts are mainly motivated to improve the computational efficiency in the stochastic particle tracking method as well as to incorporate the physical submodels of turbulence, combustion, vaporization, and dense spray effects. The present mathematical formulation and numerical methodologies can be casted in any time-marching pressure correction methodologies (PCM) such as FDNS code and MAST code. A sequence of validation cases involving steady burning sprays and transient evaporating sprays will be included

Astrophysical turbulence modeling

The role of turbulence in various astrophysical settings is reviewed. Among the differences to laboratory and atmospheric turbulence we highlight the ubiquitous presence of magnetic fields that are generally produced and maintained by dynamo action. The extreme temperature and density contrasts and stratifications are emphasized in connection with turbulence in the interstellar medium and in stars with outer convection zones, respectively. In many cases turbulence plays an essential role in facilitating enhanced transport of mass, momentum, energy, and magnetic fields in terms of the corresponding coarse-grained mean fields. Those transport properties are usually strongly modified by anisotropies and often completely new effects emerge in such a description that have no correspondence in terms of the original (non coarse-grained) fields.Comment: 88 pages, 26 figures, published in Reports on Progress in Physic

Fuel jet injection and supersonic mixing for scramjets

Imperial Users onl

A compressible Navier-Stokes solver with two-equation and Reynolds stress turbulence closure models

This report outlines the development of a general purpose aerodynamic solver for compressible turbulent flows. Turbulent closure is achieved using either two equation or Reynolds stress transportation equations. The applicable equation set consists of Favre-averaged conservation equations for the mass, momentum and total energy, and transport equations for the turbulent stresses and turbulent dissipation rate. In order to develop a scheme with good shock capturing capabilities, good accuracy and general geometric capabilities, a multi-block cell centered finite volume approach is used. Viscous fluxes are discretized using a finite volume representation of a central difference operator and the source terms are treated as an integral over the control volume. The methodology is validated by testing the algorithm on both two and three dimensional flows. Both the two equation and Reynolds stress models are used on a two dimensional 10 degree compression ramp at Mach 3, and the two equation model is used on the three dimensional flow over a cone at angle of attack at Mach 3.5. With the development of this algorithm, it is now possible to compute complex, compressible high speed flow fields using both two equation and Reynolds stress turbulent closure models, with the capability of eventually evaluating their predictive performance

ASHEE: a compressible, equilibrium-Eulerian model for volcanic ash plumes

A new fluid-dynamic model is developed to numerically simulate the non-equilibrium dynamics of polydisperse gas-particle mixtures forming volcanic plumes. Starting from the three-dimensional N-phase Eulerian transport equations for a mixture of gases and solid particles, we adopt an asymptotic expansion strategy to derive a compressible version of the first-order non-equilibrium model, valid for low concentration regimes and small particles Stokes $St<0.2$. When $St < 0.001$ the model reduces to the dusty-gas one. The new model is significantly faster than the Eulerian model while retaining the capability to describe gas-particle non-equilibrium. Direct numerical simulation accurately reproduce the dynamics of isotropic turbulence in subsonic regime. For gas-particle mixtures, it describes the main features of density fluctuations and the preferential concentration of particles by turbulence, verifying the model reliability and suitability for the simulation of high-Reynolds number and high-temperature regimes. On the other hand, Large-Eddy Numerical Simulations of forced plumes are able to reproduce their observed averaged and instantaneous properties. The self-similar radial profile and the development of large-scale structures are reproduced, including the rate of entrainment of atmospheric air. Application to the Large-Eddy Simulation of the injection of the eruptive mixture in a stratified atmosphere describes some of important features of turbulent volcanic plumes, including air entrainment, buoyancy reversal, and maximum plume height. Coarse particles partially decouple from the gas within eddies, modifying the turbulent structure, and preferentially concentrate at the eddy periphery, eventually being lost from the plume margins due to the gravity. By these mechanisms, gas-particle non-equilibrium is able to influence the large-scale behavior of volcanic plumes.Comment: 29 pages, 22 figure