A FEniCS-based programming framework for modeling turbulent flow by the Reynolds-averaged Navier-Stokes equations

Finding an appropriate turbulence model for a given flow case usually calls for extensive experimentation with both models and numerical solution methods. This work presents the design and implementation of a flexible, programmable software framework for assisting with numerical experiments in computational turbulence. The framework targets Reynolds-averaged Navier-Stokes models, discretized by finite element methods. The novel implementation makes use of Python and the FEniCS package, the combination of which leads to compact and reusable code, where model- and solver-specific code resemble closely the mathematical formulation of equations and algorithms. The presented ideas and programming techniques are also applicable to other fields that involve systems of nonlinear partial differential equations. We demonstrate the framework in two applications and investigate the impact of various linearizations on the convergence properties of nonlinear solvers for a Reynolds-averaged Navier-Stokes model. © 2011 Elsevier Ltd.

Geometrical properties and turbulent flame speed measurements in stationary premixed V-flames using direct numerical simulation

Three dimensional, fully compressible direct numerical simulations (DNS) of premixed turbulent flames are carried out in a V-flame configuration. The governing equations and the numerical implementation are described in detail, including modifications made to the Navier-Stokes Characteristic Boundary Conditions (NSCBC) to accommodate the steep transverse velocity and composition gradients generated when the flame crosses the boundary. Three cases, at turbulence intensities, u′/sL, of 1, 2, and 6 are considered. The influence of the flame holder on downstream flame properties is assessed through the distributions of the surface-conditioned displacement speed, curvature and tangential strain rates, and compared to data from similarly processed planar flames. The distributions are found to be indistinguishable from planar flames for distances greater than about 17δth downstream of the flame holder, where δth is the laminar flame thermal thickness. Favre mean fields are constructed, and the growth of the mean flame brush is found to be well described by simple Taylor type diffusion. The turbulent flame speed, sT is evaluated from an expression describing the propagation speed of an isosurface of the mean reaction progress variable c̃ in terms of the imbalance between the mean reactive, diffusive, and turbulent fluxes within the flame brush. The results are compared to the consumption speed, sC, calculated from the integral of the mean reaction rate, and to the predictions of a recently developed flame speed model (Kolla et al., Combust Sci Technol 181(3):518-535, 2009). The model predictions are improved in all cases by including the effects of mean molecular diffusion, and the overall agreement is good for the higher turbulence intensity cases once the tangential convective flux of c̃ is taken into account. © 2010 Springer Science+Business Media B.V.

Statistics and modelling of turbulent kinetic energy transport in different regimes of premixed combustion

The statistical behaviour of turbulent kinetic energy transport in turbulent premixed flames is analysed using data from three-dimensional Direct Numerical Simulation (DNS) of freely propagating turbulent premixed flames under decaying turbulence. For flames within the corrugated flamelets regime, it is observed that turbulent kinetic energy is generated within the flame brush. By contrast, for flames within the thin reaction zones regime it has been found that the turbulent kinetic energy decays monotonically through the flame brush. Similar trends are observed also for the dissipation rate of turbulent kinetic energy. Within the corrugated flamelets regime, it is demonstrated that the effects of the mean pressure gradient and pressure dilatation within the flame are sufficient to overcome the effects of viscous dissipation and are responsible for the observed augmentation of turbulent kinetic energy in the flame brush. In the thin reaction zones regime, the effects of the mean pressure gradient and pressure dilatation terms are relatively much weaker than those of viscous dissipation, resulting in a monotonic decay of turbulent kinetic energy across the flame brush. The modelling of the various unclosed terms of the turbulent kinetic energy transport equation has been analysed in detail. The predictions of existing models are compared with corresponding quantities extracted from DNS data. Based on this a-priori DNS assessment, either appropriate models are identified or new models are proposed where necessary. It is shown that the turbulent flux of turbulent kinetic energy exhibits counter-gradient (gradient) transport wherever the turbulent scalar flux is counter-gradient (gradient) in nature. A new model has been proposed for the turbulent flux of turbulent kinetic energy, and is found to capture the qualitative and quantitative behaviour obtained from DNS data for both the corrugated flamelets and thin reaction zones regimes without the need to adjust any of the model constants. © 2010 Springer Science+Business Media B.V.

Occurrence of Vibrio parahaemolyticus in marine prawns and environments

Qualitative studies on the microflora of slime and guts of prawns and of sea water off Nagapattinam showed the presence of Vibrio in the slime and sea water. They were further tested for Vibrio parahaemolyticus types and related bio-types. Evidence of its occurrence is given. This points to the need for further studies on the distribution of this organism in terms of public health significance.

Turbulent Premixed Flames

An improved understanding of lean fuel turbulent premixed flames must play a central role in the fundamental science of these new concepts.

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.

On the stabilisation of turbulent lifted flames

This study explores the stabilisation mechanisms of turbulent lifted flames by examining the scalar dissipation rate (SDR) of both passive and reactive scalars and their cross dissipation (CDR) in the stabilisation region. DNS results of a laboratory scale hydrogen turbulent lifted flame has been used for this analysis. Various definitions of the flame leading edge (FLE) has been compared and differences are illustrated. Time and spatial averaged statistic of SDR and CDR were examined. It was found that the averaged SDR for the mixture fraction at FLE was well below the reference quenching value for stoichiometric mixture. The averaged SDR for the progress variable is in the same order of the unstrained premixed laminar flame value. It was observed that the averaged CDR changed from negative to weakly positive at FLE. The change in sign was explained by a change in the relative alignment of the gradients of mixture fraction and progress variable. It was thus evident that the CDR was a good marker for stabilisation region and an important quantity in stabilisation mechanism.