Large Eddy Simulation of Premixed Turbulent Combustion Using Ξ Flame Surface Wrinkling Model

One commonly-used method for deriving the RANS equations for multicomponent flow is the technique of conditional averaging. In this paper the concept is extended to LES, by introducing the operations of conditional filtering and surface filtering. Properties of the filtered indicator function are investigated mathematically and computationally. These techniques are then used to derive conditionally filtered versions of the Navier–Stokes equations which are appropriate for simulating multicomponent flow in LES. Transport equations for the favre-averaged indicator function and the unresolved interface properties (the wrinkling and the surface area per unit volume) are also derived. Since the paper is directed towards modelling premixed combustion in the flamelet regime, closure of the equations is achieved by introducing physical models based on the picture of the flame as a wrinkled surface separating burnt and unburnt components of the fluid. This leads to a set of models for premixed turbulent combustion of varying complexity. The results of applying one of this set of models to propagation of a spherical flame in isotropic homogeneous turbulence are analysed. This revised version was published online in July 2006 with corrections to the Cover Date.     

Cahn–Hilliard/Navier–Stokes Model for the Simulation of Three-Phase Flows

In this article, we describe some aspects of the diffuse interface modelling of incompressible flows, composed of three immiscible components, without phase change. In the diffuse interface methods, system evolution is driven by the minimisation of a free energy. The originality of our approach, derived from the Cahn–Hilliard model, comes from the particular form of energy we proposed in Boyer and Lapuerta (M2AN Math Model Numer Anal, 40:653–987,2006), which, among other interesting properties, ensures consistency with the two-phase model. The modelling of three-phase flows is further completed by coupling the Cahn–Hilliard system and the Navier–Stokes equations where surface tensions are taken into account through volume capillary forces. These equations are discretized in time and space paying attention to the fact that most of the main properties of the original model (volume conservation and energy estimate) have to be maintained at the discrete level. An adaptive refinement method is finally used to obtain an accurate resolution of very thin moving internal layers, while limiting the total number of cells in the grids all along the simulation. Different numerical results are given, from the validation case of the lens spreading between two phases (contact angles and pressure jumps), to the study of mass transfer through a liquid/liquid interface crossed by a single rising gas bubble. The numerical applications are performed with large ratio between densities and viscosities and three different surface tensions.     

Investigation of the “TECFLAM” Non-premixed Flame Using Large Eddy Simulation and Proper Orthogonal Decomposition

The confined turbulent swirling non-premixed TECFLAM S09c flame has been investigated using Large Eddy Simulation and a pre-calculated mixture-fraction-based flamelet model in a preliminary step to simulate this flame. The simplified formulation, despite the fact that it cannot capture localized extinctions, is found to reproduce the experimentally-observed vortex breakdown and the results agree well with available experimental data for velocity and mixture fraction. The unsteady flow features before the burner exit and inside the combustion chamber are analyzed with spectral analysis, correlations, and Proper Orthogonal Decomposition. The results show the presence of longitudinal vortices whose axes rotate with the swirl and which cause separation inside the inlet pipe. With combustion, some of the structures are damped and the spectral peaks shift to higher frequencies. Finally, simulations with the full 3D Conditional Moment Closure equation, which allows spatial and temporal fluctuations of the conditionally-filtered reacting scalars, reproduce successfully the flame lift-off observed experimentally.     

A Spalart–allmaras Turbulence Model Implementation for High-order Discontinuous Galerkin Solution of the Reynolds-averaged Navier-stokes Equations

A novel and robust approach has been proposed for the high-order discontinuous Galerkin (DG) discretization of the Reynolds-averaged Navier-Stokes (RANS) equations with the turbulence model of Spalart-Allmaras (SA). The solution polynomials of the SA equation are reconstructed by the Hermite weighted essentially non-oscillatory (HWENO) scheme. Several practical techniques are suggested to simplify and extend a positivity-preserving limiter to further guarantee the positivity of SA working variable. The resulting positivity-preserving HWENO limiting method is compact and easy to implement on arbitrary meshes. Typical turbulent flows are conducted to assess the accuracy and robustness of the present method. Numerical experiments demonstrate that with the increasing grid or order resolution, the limited results of the working variable are getting closer to the unlimited ones. And the most obvious improvement with proposed method is on the computation of the working variable field in wake regions.     

An Experimental Study of the Modulation of the Bubble Motion by Gas–Liquid-Phase Interaction in Oscillating-Grid Decaying Turbulence

The purpose of the present research is to understand dynamic bubble–liquid interaction in a bubbly flow based on the experimental results of the modulation of the bubble motion in oscillating-grid decaying turbulence. By comparing the experimental results obtained from stagnant water and those from oscillating-grid decaying turbulence, we discussed and described detailed process of the modulation of the bubble motion in a water vessel. We discussed the enhancement of the transition of the bubble motion from 2D to 3D by combining the liquid-phase motion obtained through particle imaging velocimetry/laser-induced fluorescence (PIV/LIF) measurement and the bubble wake motion captured through the LIF/HPTS (8-hydroxypyrene-1, 3, 6-trisulfonic acid) method, under both conditions (in the stagnant water and in the oscillating-grid decaying turbulence) in which the initial bubble formation and the bubble motion (gravity-center motion and surface oscillation) were considered to be the same. In addition, by using PIV/LIF measurement along with an infrared shadow technique, we simultaneously obtained the bubble motion (2D zigzagging motion in stagnant water, and 3D motion in the decaying turbulence) and the standard deviation of the liquid-phase motion (the bubble Reynolds number: 775; the turbulent Reynolds number: 62.2). Taking all of the results together, the modulation of the bubble motion in the decaying turbulence, and the dynamic interaction between the bubble and the liquid-phase motion were experimentally and carefully investigated. Consequently, the enhancement and the modulation of the bubble wake motion were considered to be triggered by the collapse of the symmetric property of the bubble–liquid (i.e. ambient liquid-phase turbulence) interaction.     

Simulation of Structure Change in Porous Media During Gas–Solid Reactions Using Cellular Automata Model

Transport in Porous Media - In some gas–solid reactions, a new solid substance is produced. The product acts as a shield and prevents the collision between gas and solid reactants which...     

Distinguishing Features of One Model of an Elastic Solid Associated with the Long–Range Interaction at the Molecular Level

A model of an elastic solid in the form of a system of elastically connected rigid elements is proposed. It is shown that the long–range interaction should be taken into account. The mathematical model proposed is, in essence, the physical model of a solid, which substantially broadens the range of its application.     

Numerical Simulation of the Physical–Chemical–Thermal Processes During Hydration Reaction of the Calcium Oxide/Calcium Hydroxide System in an Indirect Reactor

Transport in Porous Media - Solutal convection in a horizontal layer filled with porous media with horizontal seepage of a mixture is investigated considering the solute immobilization and...     

Using S-statistic for investigating the effect of temperature on hydrodynamics of gas–solid fluidization

The influence of temperature on fluidization was investigated by a statistical chaotic attractor comparison test known as S-statistic. After calibration of the variables used in this method, the S-test was applied to the radioactive particle tracking (RPT) data obtained from a lab-scale fluidized bed. Experiments were performed with sand as fluidized particles and in temperatures from ambient up to 600 °C with superficial gas velocities of 0.29, 0.38 and 0.52 m/s. Considering the behavior of bubbles and comparing with frequency domain analysis, it was concluded that S-statistic is a reliable method for characterization of fluidization process behavior at different temperatures.     

Digital Image Correlation Analysis and Numerical Simulation of Aluminum Alloys under Quasi-static Tension after Necking Using the Bridgman’s Correction Method

Experimental Techniques - Quasi-static tensile test is a common, yet fundamental, experiment in determining the mechanical properties of materials. Often, the determination of the equivalent...     

Correction to: A New “λ2” Term for the Spalart–Allmaras Turbulence Model,Active in Axisymmetric Flows

Flow, Turbulence and Combustion -     

Development of a Generalized Version of the Poisson– Nernst–Planck Equations Using the Hybrid Mixture Theory: Presentation of 2D Numerical Examples

A numerical scheme for the transient solution of a generalized version of the Poisson–Nernst–Planck (PNP) equations is presented. The finite element method is used to establish the coupled non-linear matrix system of equations capable of solving the present problem iteratively. The PNP equations represent a set of diffusion equations for charged species, i.e. dissolved ions, present in the pore solution of a rigid porous material in which the surface charge can be assumed neglectable. These equations are coupled to the ‘internally’ induced electrical field and to the velocity field of the fluid. The Nernst–Planck equations describing the diffusion of the ionic species and Gauss’ law in use are, however, coupled in both directions. The governing set of equations is derived from a simplified version of the so-called hybrid mixture theory (HMT). The simplifications used here mainly concerns ignoring the deformation and stresses in the porous material in which the ionic diffusion occurs. The HMT is a special version of the more ‘classical’ continuum mixture theories in the sense that it works with averaged equations at macroscale and that it includes the volume fractions of phases in its structure. The background to the PNP equations can by the HMT approach be described by using the postulates of mass conservation of constituents together with Gauss’ law used together with consistent constitutive laws. The HMT theory includes the constituent forms of the quasistatic version of Maxwell’s equations making it suitable for analyses of the kind addressed in this work. Within the framework of HTM, constitutive equations have been derived using the postulate of entropy inequality together with the technique of identifying properties by Lagrange multipliers. These results will be used in obtaining a closed set of equations for the present problem.     

The Multigrain Two-Scale Model for the Ziegler–Natta Polymerization Process with Fragmentation of the Catalytic Aggregate

In this paper we construct a mathematical model for the Ziegler–Natta polymerization process which includes the initial stage of fragmentation of the support. The problem consists in a system of partial differential equations with free boundaries and extends the previous model proposed in [2, 3]. Local existence and uniqueness of a regular solution are proved.     

Verification of Wide-Range Constitutive Relations for Elastic-Viscoplastic Materials Using the Taylor–Hopkinson Test

Journal of Applied Mechanics and Technical Physics - A mathematical model of a solid body with mesoscopic defects is presented and validated. The constitutive relations proposed earlier allow...