A Fluid-Dynamical Subgrid Scale Model for Highly Compressible Astrophysical Turbulence

We formulate and implement the Euler equations with SGS dynamics and provide numerical tests of an SGS turbulence energy model that predicts the turbulent pressure of unresolved velocity fluctuations and the rate of dissipation for highly compressible turbulence. We test closures for the turbulence energy cascade by filtering data from high-resolution simulations of forced isothermal and adiabatic turbulence. Optimal properties and an excellent correlation are found for a linear combination of the eddy-viscosity closure that is employed in LES of weakly compressible turbulence and a term that is non-linear in the Jacobian matrix of the velocity. Using this mixed closure, the SGS turbulence energy model is validated in LES of turbulence with stochastic forcing. It is found that the SGS model satisfies several important requirements: 1. The mean SGS turbulence energy follows a power law for varying grid scale. 2. The root mean square (RMS) Mach number of the unresolved velocity fluctuations is proportional to the RMS Mach number of the resolved turbulence, independent of the forcing. 3. The rate of dissipation and the turbulence energy flux are constant. Moreover, we discuss difficulties with direct estimates of the turbulent pressure and the dissipation rate on the basis of resolved flow quantities that have recently been proposed. In combination with the energy injection by stellar feedback and other unresolved processes, the proposed SGS model is applicable to a variety of problems in computational astrophysics. Computing the SGS turbulence energy, the treatment of star formation and stellar feedback in galaxy simulations can be improved. Further, we expect that the turbulent pressure on the grid scale affects the stability of gas against gravitational collapse.Comment: 19 pages, 16 figures, submitted to A&

Inlet conditions for large eddy simulation of gas-turbine swirl injectors

Copyright © 2008 American Institute of Aeronautics and AstronauticsIn this paper, we present a novel technique for generating swirl inlets for large eddy simulation. The velocity a short distance downstream of the inlet to the main domain is sampled and the flow velocity data are reintroduced back into the domain inlet, creating an inlet section integrated into the main domain in which turbulence can develop. Additionally, variable artificial body forces and velocity corrections are imposed in this inlet section, with feedback control to force the flow toward desired swirl, mean, and turbulent profiles. The method was applied to flow in an axisymmetric sudden expansion, with and without swirl at the inlet, and compared against experimental and literature large eddy simulation data and against similar results in the literature. The method generates excellent results for this case and is elegant and straightforward to implement

Investigations of microwave stimulation of a turbulent low-swirl flame

Irradiating a flame by microwave radiation is one of several plasma-assisted combustion (PAC) technologies that can be used to modify the combustion chemical kinetics in order to improve flame-stability and to delay lean blow-out. One practical implication is that engines may be able to operate with leaner fuel mixtures and have an improved fuel flexibility capability including biofuels. In addition, this technology may assist in reducing thermoacoustic instabilities that may severely damage the engine and increase emission production. To examine microwave-assisted combustion a combined experimental and computational study of microwave-assisted combustion is performed for a lean, turbulent, swirl-stabilized, stratified flame at atmospheric conditions. The objectives are to demonstrate that the technology increases both the laminar and turbulent flame speeds, and modifies the chemical kinetics, enhancing the flame-stability at lean mixtures. The study combines experimental investigations using hydroxyl (OH) and formaldehyde (CH2O) Planar Laser-Induced Fluorescence (PLIF) and numerical simulations using finite rate chemistry Large Eddy Simulations (LES). The reaction mechanism is based on a methane (CH4)-air skeletal mechanism expanded with sub-mechanisms for ozone, singlet oxygen, chemionization, electron impact dissociation, ionization and attachment. The experimental and computational results show similar trends, and are used to demonstrate and explain some significant aspects of microwave-enhanced combustion. Both simulation and experimental studies are performed close to lean blow off conditions. In the simulations, the flame is gradually subjected to increasing reduced electric field strengths, resulting in a wider flame that stabilizes nearer to the burner nozzle. Experiments are performed at two equivalence ratios, where the leaner case absorbs up to more than 5% of the total flame power. Data from experiments reveal trends similar to simulated results with increased microwave absorption

Using LES to Study Reacting Flows and Instabilities in Annular Combustion Chambers

Great prominence is put on the design of aeronautical gas turbines due to increasingly stringent regulations and the need to tackle rising fuel prices. This drive towards innovation has resulted sometimes in new concepts being prone to combustion instabilities. In the particular field of annular combustion chambers, these instabilities often take the form of azimuthal modes. To predict these modes, one must compute the full combustion chamber, which remained out of reach until very recently and the development of massively parallel computers. Since one of the most limiting factors in performing Large Eddy Simulation (LES) of real combustors is estimating the adequate grid, the effects of mesh resolution are investigated by computing full annular LES of a realistic helicopter combustion chamber on three grids, respectively made of 38, 93 and 336 million elements. Results are compared in terms of mean and fluctuating fields. LES captures self-established azimuthal modes. The presence and structure of the modes is discussed. This study therefore highlights the potential of LES for studying combustion instabilities in annular gas turbine combustors

Large-Eddy Simulation: Current Capabilities, Recommended Practices, and Future Research

This paper presents the results of an activity by the Large Eddy Simulation (LES) Working Group of the AIAA Fluid Dynamics Technical Committee to (1) address the current capabilities of LES, (2) outline recommended practices and key considerations for using LES, and (3) identify future research needs to advance the capabilities and reliability of LES for analysis of turbulent flows. To address the current capabilities and future needs, a survey comprised of eleven questions was posed to LES Working Group members to assemble a broad range of perspectives on important topics related to LES. The responses to these survey questions are summarized with the intent not to be a comprehensive dictate on LES, but rather the perspective of one group on some important issues. A list of recommended practices is also provided, which does not treat all aspects of a LES, but provides guidance on some of the key areas that should be considered

Mixing in Circular and Non-circular Jets in Crossflow

Coherent structures and mixing in the flow field of a jet in crossflow have been studied using computational (large eddy simulation) and experimental (particle image velocimetry and laser-induced fluorescence) techniques. The mean scalar fields and turbulence statistics as determined by both are compared for circular, elliptic, and square nozzles. For the latter configurations, effects of orientation are considered. The computations reveal that the distribution of a passive scalar in a cross-sectional plane can be single- or double-peaked, depending on the nozzle shape and orientation. A proper orthogonal decomposition of the transverse velocity indicates that coherent structures may be responsible for this phenomenon. Nozzles which have a single-peaked distribution have stronger modes in transverse direction. The global mixing performance is superior for these nozzle types. This is the case for the blunt square nozzle and for the elliptic nozzle with high aspect ratio. It is further demonstrated that the flow field contains large regions in which a passive scalar is transported up the mean gradient (counter-gradient transport) which implies failure of the gradient diffusion hypothesis