Conditional Moment Closure for Turbulent Premixed Flames

The conditional moment closure (CMC) method has been successfully applied to various non-premixed combustion systems in the past, but its application to premixed flames is not fully tested and validated. The main difficulty is associated with the modeling of conditional scalar dissipation rate of the conditioning scalar, the progress variable. A simple algebraic model for the conditional dissipation rate is validated using DNS results of a V-flame. This model along with the standard k- turbulence modeling is used in computations of stoichiometric pilot stabilized Bunsen flames using the RANS-CMC method. A first-order closure is used for the conditional mean reaction rate. The computed non reacting and reacting scalars are in reasonable agreement with the experimental measurements and are consistent with earlier computations using flamelets and transported PDF methods. Sensitivity to chemical kinetic mechanism is also assessed. The results suggest that the CMC may be applied across the regimes of premixed combustion.

Conditional Moment Closure for Turbulent Premixed Flames

The conditional moment closure (CMC) method has been successfully applied to various non-premixed combustion systems in the past, but its application to premixed flames is not fully tested and validated. The main difficulty is associated with the modeling of conditional scalar dissipation rate of the conditioning scalar, the progress variable. A simple algebraic model for the conditional dissipation rate is validated using DNS results of a V-flame. This model along with the standard k- turbulence modeling is used in computations of stoichiometric pilot stabilized Bunsen flames using the RANS-CMC method. A first-order closure is used for the conditional mean reaction rate. The computed non reacting and reacting scalars are in reasonable agreement with the experimental measurements and are consistent with earlier computations using flamelets and transported PDF methods. Sensitivity to chemical kinetic mechanism is also assessed. The results suggest that the CMC may be applied across the regimes of premixed combustion.

Computations of turbulent lean premixed combustion using conditional moment closure

Conditional Moment Closure (CMC) is a suitable method for predicting scalars such as carbon monoxide with slow chemical time scales in turbulent combustion. Although this method has been successfully applied to non-premixed combustion, its application to lean premixed combustion is rare. In this study the CMC method is used to compute piloted lean premixed combustion in a distributed combustion regime. The conditional scalar dissipation rate of the conditioning scalar, the progress variable, is closed using an algebraic model and turbulence is modelled using the standard k-e{open} model. The conditional mean reaction rate is closed using a first order CMC closure with the GRI-3.0 chemical mechanism to represent the chemical kinetics of methane oxidation. The PDF of the progress variable is obtained using a presumed shape with the Beta function. The computed results are compared with the experimental measurements and earlier computations using the transported PDF approach. The results show reasonable agreement with the experimental measurements and are consistent with the transported PDF computations. When the compounded effects of shear-turbulence and flame are strong, second order closures may be required for the CMC. © 2013 Copyright Taylor and Francis Group, LLC.

Influence of turbulence-chemistry interaction for n-heptane spray combustion under diesel engine conditions with emphasis on soot formation and oxidation

The influence of the turbulence-chemistry interaction (TCI) for n-heptane sprays under diesel engine conditions has been investigated by means of computational fluid dynamics (CFD) simulations. The conditional moment closure approach, which has been previously validated thoroughly for such flows, and the homogeneous reactor (i.e. no turbulent combustion model) approach have been compared, in view of the recent resurgence of the latter approaches for diesel engine CFD. Experimental data available from a constant-volume combustion chamber have been used for model validation purposes for a broad range of conditions including variations in ambient oxygen (8-21% by vol.), ambient temperature (900 and 1000 K) and ambient density (14.8 and 30 kg/m3). The results from both numerical approaches have been compared to the experimental values of ignition delay (ID), flame lift-off length (LOL), and soot volume fraction distributions. TCI was found to have a weak influence on ignition delay for the conditions simulated, attributed to the low values of the scalar dissipation relative to the critical value above which auto-ignition does not occur. In contrast, the flame LOL was considerably affected, in particular at low oxygen concentrations. Quasi-steady soot formation was similar; however, pronounced differences in soot oxidation behaviour are reported. The differences were further emphasised for a case with short injection duration: in such conditions, TCI was found to play a major role concerning the soot oxidation behaviour because of the importance of soot-oxidiser structure in mixture fraction space. Neglecting TCI leads to a strong over-estimation of soot oxidation after the end of injection. The results suggest that for some engines, and for some phenomena, the neglect of turbulent fluctuations may lead to predictions of acceptable engineering accuracy, but that a proper turbulent combustion model is needed for more reliable results. © 2014 Taylor & Francis.

Turbulence modeling for flows around convex features giving rapid eddy distortion

Reynolds averaged Navier-Stokes model performances in the stagnation and wake regions for turbulent flows with relatively large Lagrangian length scales (generally larger than the scale of geometrical features) approaching small cylinders (both square and circular) is explored. The effective cylinder (or wire) diameter based Reynolds number, ReW ≤ 2.5 × 103. The following turbulence models are considered: a mixing-length; standard Spalart and Allmaras (SA) and streamline curvature (and rotation) corrected SA (SARC); Secundov's νt-92; Secundov et al.'s two equation νt-L; Wolfshtein's k-l model; the Explicit Algebraic Stress Model (EASM) of Abid et al.; the cubic model of Craft et al.; various linear k-ε models including those with wall distance based damping functions; Menter SST, k-ω and Spalding's LVEL model. The use of differential equation distance functions (Poisson and Hamilton-Jacobi equation based) for palliative turbulence modeling purposes is explored. The performance of SA with these distance functions is also considered in the sharp convex geometry region of an airfoil trailing edge. For the cylinder, with ReW ≈ 2.5 × 103 the mixing length and k-l models give strong turbulence production in the wake region. However, in agreement with eddy viscosity estimates, the LVEL and Secundov νt-92 models show relatively little cylinder influence on turbulence. On the other hand, two equation models (as does the one equation SA) suggest the cylinder gives a strong turbulence deficit in the wake region. Also, for SA, an order or magnitude cylinder diameter decrease from ReW = 2500 to 250 surprisingly strengthens the cylinder's disruptive influence. Importantly, results for ReW ≪ 250 are virtually identical to those for ReW = 250 i.e. no matter how small the cylinder/wire its influence does not, as it should, vanish. Similar tests for the Launder-Sharma k-ε, Menter SST and k-ω show, in accordance with physical reality, the cylinder's influence diminishing albeit slowly with size. Results suggest distance functions palliate the SA model's erroneous trait and improve its predictive performance in wire wake regions. Also, results suggest that, along the stagnation line, such functions improve the SA, mixing length, k-l and LVEL results. For the airfoil, with SA, the larger Poisson distance function increases the wake region turbulence levels by just under 5%. © 2007 Elsevier Inc. All rights reserved.

Flux closure vortexlike domain structures in ferroelectric thin films

Enhanced piezoresponse force microscopy was used to study flux closure vortexlike structures of 90° ferroelastic domains at the nanoscale in thin ferroelectric lead zirconium titanate (PZT) films. Using an external electric field, a vortexlike structure was induced far away from a grain boundary, indicating that physical edges are not necessary for nucleation contrary to previous suggestions. We demonstrate two different configurations of vortexlike structures, one of which has not been observed before. The stability of these structures is found to be size dependent, supporting previous predictions. © 2010 The American Physical Society.