ALBERT

Description: In premixed turbulent combustion, the modelling of the turbulent flux of the mean reaction progress variable remains somewhat controversial. Classical gradient transport assumptions based on the eddy viscosity concept are often used while both experimental data and theoretical analysis have pointed out the existence of countergradient turbulent diffusion. Direct numerical simulation (DNS) is used in this paper to provide basic information on the turbulent flux of and study the occurrence of counter-gradient transport. The numerical configuration corresponds to twoor three-dimensional premixed flames in isotropic turbulent flow. The simulations correspond to various flame and flow conditions that are representative of flamelet combustion. They reveal that different flames will feature different turbulent transport properties and that these differences can be related to basic dynamical differences in the flame-flow interactions: counter-gradient diffusion occurs when the flow field near the flame is dominated by thermal dilatation due to chemical reaction, whereas gradient diffusion occurs when the flow field near the flame is dominated by the turbulent motions. The DNS-based analysis leads to a simple expression to describe the turbulent flux of , which in turn leads to a simple criterion to delineate between the gradient and counter-gradient turbulent diffusion regimes. This criterion suggests that the occurrence of one regime or the other is determined primarily by the ratio of turbulence intensity divided by the laminar flame speed, and by the flame heat release factor, τ ≡ (Tb — Tu)/Tu, where Tu and Tb are respectively the temperature within unburnt and burnt gas. Consistent with the Bray-Moss-Libby theory, counter-gradient (gradient) diffusion is promoted by low (high) values and high (low) values of τ. DNS also shows that these results are not restricted to the turbulent transport of . Similar results are found for the turbulent transport of flame surface density, Σ. The turbulent fluxes of and Σ are strongly correlated in the simulated flames and counter-gradient (gradient) diffusion of always coincides with counter-gradient (gradient) diffusion of Σ.

Description: We consider a mathematical model which describes the quasistatic contact between two deformable bodies. The bodies are assumed to have a viscoelastic behavior that we model with Kelvin-Voigt constitutive law. The contact is frictionless and is modeled with the classical Signorini condition with zero-gap function. We derive a variational formulation of the problem and prove the existence of a unique weak solution to the model by using arguments of evolution equations with maximal monotone operators. We also prove that the solution converges to the solution of the corresponding elastic problem, as the viscosity tensors converge to zero. We then consider a fully discrete approximation of the problem, based on the augmented Lagrangian approach, and present numerical results of two-dimensional test problems.

Description: Group velocity of wave packet propagating radially across plasma immersed in static axial magnetic field, examining ringing phenomena due to pulsed transmission

Description: Although the different regimes of premixed combustion are not well defined, most of the recent developments in turbulent combustion modeling are led in the so-called flamelet regime. The goal of these models is to give a realistic expression to the mean reaction rate (w). Several methods can be used to estimate (w). Bray and coworkers (Libby & Bray 1980, Bray 1985, Bray & Libby 1986) express the instantaneous reaction rate by means of a flamelet library and a frequency which describes the local interaction between the laminar flamelets and the turbulent flowfield. In another way, the mean reaction rate can be directly connected to the flame surface density (Sigma). This quantity can be given by the transport equation of the coherent flame model initially proposed by Marble & Broadwell 1977 and developed elsewhere. The mean reaction rate, (w), can also be estimated thanks to the evolution of an arbitrary scalar field G(x, t) = G(sub O) which represents the flame sheet. G(x, t) is obtained from the G-equation proposed by Williams 1985, Kerstein et al. 1988 and Peters 1993. Another possibility proposed in a recent study by Mantel & Borghi 1991, where a transport equation for the mean dissipation rate (epsilon(sub c)) of the progress variable c is used to determine (w). In their model, Mantel & Borghi 1991 considered a medium with constant density and constant diffusivity in the determination of the transport equation for (epsilon(sub c)). A comparison of different flamelet models made by Duclos et al. 1993 shows the realistic behavior of this model even in the case of constant density. Our objective in this present report is to present preliminary results on the study of this equation in the case of variable density and variable diffusivity. Assumptions of constant pressure and a Lewis number equal to unity allow us to significantly simplify the equation. A systematic order of magnitude analysis based on adequate scale relations is performed on each term of the equation. As in the case of constant density and constant diffusivity, the effects of stretching of the scalar field by the turbulent strain field, of local curvature, and of chemical reactions are predominant. In this preliminary work, we suggest closure models for certain terms, which will be validated after comparisons with DNS data.

Description: In premixed turbulent flames the presence of intense mixing zones located in front of and behind the flame surface leads to a requirement to study the behavior of iso-concentration surfaces defined for all values of the progress variable (equal to unity in burnt gases and to zero in fresh mixtures). To support this study, some theoretical and mathematical tools devoted to level surfaces are first developed. Then a database of direct numerical simulations of turbulent premixed flames is generated and used to investigate the internal structure of the flame brush, and a new pdf model based on the properties of iso-surfaces is proposed.

Description: Current flamelot models based on a description of the flame surface dynamics require the closure of two inter-related equations: a transport equation for the mean reaction progress variable, (tilde)c, and a transport equation for the flame surface density, Sigma. The coupling between these two equations is investigated using direct numerical simulations (DNS) with emphasis on the correlation between the turbulent fluxes of (tilde)c, bar(pu''c''), and Sigma, (u'')(sub S)Sigma. Two different DNS databases are used in the present work: a database developed at CTR by A. Trouve and a database developed by C. J. Rutland using a different code. Both databases correspond to statistically one-dimensional premixed flames in isotropic turbulent flow. The run parameters, however, are significantly different, and the two databases correspond to different combustion regimes. It is found that in all simulated flames, the correlation between bar(pu''c'') and (u'')(sub S)Sigma is always strong. The sign, however, of the turbulent flux of (tilde)c or Sigma with respect to the mean gradients, delta(tilde)c/delta(x) or delta(Sigma)/delta(x), is case-dependent. The CTR database is found to exhibit gradient turbulent transport of (tilde)c and Sigma, whereas the Rutland DNS features counter-gradient diffusion. The two databases are analyzed and compared using various tools (a local analysis of the flow field near the flame, a classical analysis of the conservation equation for (tilde)(u''c''), and a thin flame theoretical analysis). A mechanism is then proposed to explain the discrepancies between the two databases and a preliminary simple criterion is derived to predict the occurrence of gradient/counter-gradient turbulent diffusion.

Description: In this paper, a new methodology to determine kinetic parameters for simple chemical models and simple transport properties classically used in DNS of premixed combustion is presented. First, a one-dimensional code is utilized to performed steady unstrained laminar methane-air flame in order to verify intrinsic features of laminar flames such as burning velocity and temperature and concentration profiles. Second, the flame response to steady and unsteady strain in the opposed jet configuration is numerically investigated. It appears that for a well determined set of parameters, one- and two-step mechanisms reproduce the extinction limit of a laminar flame submitted to a steady strain. Computations with the GRI-mech mechanism (177 reactions, 39 species) and multicomponent transport properties are used to validate these simplified models. A sensitivity analysis of the preferential diffusion of heat and reactants when the Lewis number is close to unity indicates that the response of the flame to an oscillating strain is very sensitive to this number. As an application of this methodology, the interaction between a two-dimensional vortex pair and a premixed laminar flame is performed by Direct Numerical Simulation (DNS) using the one- and two-step mechanisms. Comparison with the experimental results of Samaniego et al. (1994) shows a significant improvement in the description of the interaction when the two-step model is used.