Large-Eddy Simulation of transcritical flows

The present study uses the LES code AVBP, developed at CERFACS, to simulate transcritical flows. Real gas effects are accounted for by the use of a cubic equation of state, in conjunction with appropriate viscosity and thermal conductivity coefficients. First a single nitrogen round jet at supercritical pressure injected in a gaseous reservoir is simulated. Two cases are considered, one demonstrating a transcritical injection (high density injection), the other being directly injected at supercritical temperature (lower density injection). Comparison with available measurements shows good agreement. Finally, the simulation of a reacting case from the Mascotte bench (ONERA) is performed, consisting in a single coaxial injector injecting transcritical oxygen and supercritical hydrogen in a 60 bar chamber. Mean flow characteristics are in good agreement with the experimental observations of OH^* emission, whereas temperature comparisons are more difficult to interpret. To cite this article: T. Schmitt et al., C. R. Mecanique 337 (2009).     

The structure of multidimensional strained flames under transcritical conditions

Strained flames are commonly used to study the structure of reactive layers and describe the local properties of turbulent combustion. This model is attractive because constant strain rate flames only depend on a transverse coordinate and can be treated as a one-dimensional problem. This configuration is considered in a multidimensional context in which the strained flow is obtained by two counterflowing streams of reactants. It is used to examine the structure of transcritical strained flames in which one or two reactants are injected at a high pressure exceeding the critical value while their temperature is below the critical value. Calculations are carried out in a two-dimensional domain to test numerical models developed for multidimensional simulations and test thermodynamic and transport models devised to deal with high pressure real gas effects. Multidimensional strained flame calculations carried out in this study serve to check the validity of a new version of a Navier–Stokes flow solver (AVBP) conceived to deal with transcritical combustion of interest to liquid propellant rocket applications. This article describes the basic elements of such simulations and discusses results of calculations. It is shown that the calculated multidimensional strained flames have the expected features in terms of structure and response to the imposed strain rate. To cite this article: L. Pons et al., C. R. Mecanique 337 (2009).     

Comparison of numerical methods and combustion models for LES of a ramjet

Ramjets are very sensitive to instabilities and their numerical predictions can only be addressed adequately by Large Eddy Simulation (LES). With this technique, solvers can be implicit or explicit and handle structured, unstructured or hybrid meshes, etc. Turbulence and combustion models are other sources of differences. The impact of these options is here investigated for the ONERA ramjet burner. To do so, two LES codes developed by ONERA and CERFACS compute one stable operating condition. Preliminary LES results of the two codes underline the overall robustness of LES. Mean flow features at the various critical sections are reasonably well predicted by both codes. Disagreement mainly appear in the chamber where combustion positions differ pointing to the importance of the combustion and subgrid mixing models. The two LES produce different energy containing motions. With CEDRE, a low frequency dominates while AVBP produces different ranges of low frequencies that can be linked with acoustic modes of the configuration. To cite this article: A. Roux et al., C. R. Mecanique 337 (2009).     

Injection coupling with high amplitude transverse modes: Experimentation and simulation

High frequency combustion instabilities have technical importance in the design of liquid rocket engines. These phenomena involve a strong coupling between transverse acoustic modes and combustion. They are currently being investigated by combining experimentation and numerical simulations. On the experimental level, the coupling is examined in a model scale system featuring a multiple injector combustor (MIC) comprising five coaxial injectors fed with liquid oxygen and gaseous methane. This system is equipped with a novel VHAM actuator (Very High Amplitude Modulator) which comprises two nozzles and a rotating toothed wheel blocking the nozzles in an alternate fashion. This device was designed to obtain the highest possible levels of transverse oscillation in the MIC. After a brief review of the VHAM, this article reports cold flow experiments using this modulator. Velocity maps obtained under resonant conditions using the VHAM are examined at different instants during a cycle of oscillation. Experimental data are compared with numerical pressure and velocity fields obtained from an acoustic solver. The good agreement observed in the nozzle vicinity indicates that numerical simulations can be used to analyze the complex flow field generated by the VHAM. To cite this article: Y. Mery et al., C. R. Mecanique 337 (2009).     

Effect of multiperforated plates on the acoustic modes in combustors

The analytical model derived by Howe assessing the acoustic effect of perforated plates has been implemented in a 3D Helmholtz solver. This solver allows one to compute the acoustic modes of industrial chambers taking into account the multiperforated plates present for the cooling of the walls. An academic test case consisting of two coaxial cylinders, with the inner one being perforated is used to validate the implementation in the general purpose AVSP code. This case is also used to show the effects of the presence of the plates. In particular, the sensitivity of the acoustic damping to the bias flow speed will be studied. A maximum absorption speed is shown, and the behaviour towards an infinite speed will be illustrated by the academic case. Computations are also conducted in the case of an industrial helicopter chamber. The value of the maximum absorption speed is discussed to explain why the modes are in fact not much absorbed by the perforated plates, and that the frequencies are the same as for walls. To cite this article: E. Gullaud et al., C. R. Mecanique 337 (2009).     

Using data assimilation in numerical simulations of experimental geophysical flows

Data assimilation is used to couple numerical simulations and laboratory experiments of unsteady fluid flows in a stratified, rotating fluid. The experiments are performed on the large Coriolis turntable (Grenoble) and the simulations are performed with a multi-layer shallow water model. Sequential assimilation of high-resolution CIV (Correlation Image Velocimetry) measurements drives the numerical model close to the experimental flow and provides an estimation of all the flow variables at each time and each point. It is then possible (i) to analyse the flow dynamics in details, (ii) to determine the model errors starting from a realistic initial condition and (iii) to test the assimilation scheme when a reduced set of data is assimilated. To illustrate this, some results on the baroclinic instability of a two-layer vortex are presented. To cite this article: M. Galmiche et al., C. R. Mecanique 331 (2003).     

Numerical and analytical investigation of the indirect combustion noise in a nozzle

Analytical and numerical assessments of the indirect noise generated through a nozzle are presented. The configuration corresponds to an experimental setup operated at DLR by Bake et al. (2008) where an entropy wave is generated upstream of the nozzle by means of an electrical heating device. Both 3-D and 2-D axisymmetric simulations are performed to demonstrate that the experiment is mostly driven by linear acoustic phenomena, including pressure wave reflection at the outlet and entropy-to-acoustic conversion in the accelerated regions. Results show that the acoustic impedance downstream of the nozzle must be accounted for appropriately in order to recover the experimental pressure signal. A good agreement is also obtained with a purely analytical assessment based on the Marble and Candel compact nozzle approximation. To cite this article: M. Leyko et al., C. R. Mecanique 337 (2009).     

A DES method applied to a Backward Facing Step reactive flow

Hybrid Reynolds Averaged Navier Stokes–Large Eddy Simulation is a trend which is becoming of common use in aerodynamics but has seldom been employed to simulate reactive flows. Such methods, like the Delayed Detached Eddy Simulation (DDES) presented in this article, have been created to treat near wall flows with a RANS approach while switching to LES in the separated flow region. It is indeed an affordable solution to simulate complex and unsteady compressible flows and to have access to accurate skin friction and wall thermal fluxes. In order to validate this technique in combustion, we chose a simple and well documented Backward Facing Step combustor. To account for turbulent combustion a Dynamic Thickened Flame was used. The results obtained on this case show a good agreement with the experimental database and are of the same quality as LES in the separated region for both inert and reactive flows. To cite this article: B. Sainte-Rose et al., C. R. Mecanique 337 (2009).     

Conjugate heat transfer with Large Eddy Simulation for gas turbine components

CHT (Conjugate Heat Transfer) is a main design constraint for GT (gas turbines). Most existing CHT tools are developed for chained, steady phenomena. A fully parallel environment for CHT has been developed and applied to two configurations of interest for the design of GT. A reactive Large Eddy Simulations code and a solid conduction solver exchange data via a supervisor. A flame/wall interaction is used to assess the precision and the order of the coupled solutions. A film-cooled turbine vane is then studied. Thermal conduction in the blade implies lower wall temperature than adiabatic results and CHT reproduces the experimental cooling efficiency. To cite this article: F. Duchaine et al., C. R. Mecanique 337 (2009).     

Chemical kinetics modeling and LES combustion model effects on a perfectly premixed burner

Chemical kinetics modeling and coupling with turbulent combustion models for compressible Large Eddy Simulations (LES) is a critical issue. Accurate flow predictions can only be guaranteed if the coupling is well mastered. In a first attempt to qualify the effect of each model, the case of a lean premixed swirled combustor with comprehensive measures is targeted (species mass fractions and temperature fields). For the investigation, two turbulent combustion models are considered. The first model relies on a presumed PDF approach coupled to a look-up chemistry table obtained with a reduced chemical scheme. The second model makes use of the thickened flame approach using the same reduced chemical scheme but with reaction rates computed explicitly as the computation advances. Then, to estimate kinetic schemes reduction effects, the first model is compared to a third one, with the same PDF approach, but coupled to a look-up chemistry table obtained with a complete chemical scheme. All LES are very close to each other. The main difference between the different predictions relies on CO mass fractions. Although they are all able to return good outlet mass fractions, CO values inside the flame are different depending on the model used. To cite this article: G. Albouze et al., C. R. Mecanique 337 (2009).     

Fuel injection model for Euler–Euler and Euler–Lagrange large-eddy simulations of an evaporating spray inside an aeronautical combustor

Large-Eddy Simulations (LES) of an evaporating two-phase flow in an experimental burner are investigated. Two different numerical approaches for the simulation of the dispersed phase are coupled to the same gaseous solver: a mesoscopic Eulerian method and a Lagrangian particle tracking technique. The spray is represented by a single droplet size owing to the locally monodisperse formulation of the employed mesoscopic Eulerian approach. Both approaches use the same drag and evaporation models. They do not take into account the atomization process and a simplified injection model is applied instead. The presented methodology, referred as FIM-UR (Fuel Injection Method by Upstream Reconstruction) defines injection profiles for the monodisperse spray produced by a pressure-swirl atomizer. It is designed so as to ensure similar spray characteristics for both approaches and allows for a direct comparison between them. After a validation of the purely gaseous flow in the burner, liquid-phase dynamics and droplet dispersion are qualitatively and quantitatively evaluated for the Eulerian and Lagrangian simulations. Results obtained for both approaches are in very good agreement and compare reasonably with experiments, indicating that simplified injection methods are appropriate for the simulation of realistic combustor geometries.     

Stabilization of non-premixed flames in supersonic reactive flows

A Lagrangian framework is set out to describe turbulent non-premixed combustion in high speed coflowing jet flows. The final aim is to provide a robust computational methodology to simulate, in various conditions, the underexpanded GH2/GO2 torch jet that is used to initiate combustion in an expander cycle engine. The proposed approach relies on an early modelling proposal of Borghi and his coworkers. The model is well suited to describe finite rate chemistry effects and its recent extension to high speed flows allows one to take the influence of viscous dissipation phenomena into account. Indeed, since the chemical source terms are highly temperature sensitive, the influence of viscous phenomena on the thermal runaway is likely to be all the more pronounced since the Mach number values are high. The validation of the extended model has been recently performed through the numerical simulation of two distinct well-documented experimental databases. Only a brief summary of this preliminary validation step is provided here. The main purpose of the present work is to proceed with the numerical simulation of geometries that bring together the essential peculiarities of the underexpanded GH2/GO2 torch. The behavior of the corresponding supersonic coflowing jet flames for various conditions is discussed in the light of computational results. To cite this article: J.-F. Izard, A. Mura, C. R. Mecanique 337 (2009).     

Full coverage film cooling using compound angle

An experimental and numerical study is carried out on a cooling film issuing from a multiperforated wall of a simplified combustor. The objectives of this work are to achieve a better understanding of the dynamics of the film and to construct an experimental database on a simplified geometry in order to test numerical models. A parametric study of film cooling efficiency based on the direction of the cooling air injection is presented and shows that a swirling injection greatly enhances the cooling efficiency. As accounting for multiperforated walls in numerical simulations cannot be done at the jets scale because of computing resources, in this article are presented RANS computations performed using a uniform boundary condition to provide the injection of coolant. Two injection models are applied on this boundary and numerical results are compared to experimental data in the recovery region. The standard model is shown to be totally inappropriate while the multiperforation model delivers promising results although some weaknesses appear very close to the wall. To cite this article: B. Michel et al., C. R. Mecanique 337 (2009).     

Confined buckling of inextensible rods by convex difference algorithms

In this Note we present an approach to determine the local minima of a specific class of minimization problems. Attention is focused on the inextensibility condition of flexible rods expressed as a nonconvex constraint. Two algorithms are derived from a special splitting of the Lagrangian into the difference of two convex functions (DC). They are compared to the augmented Lagrangian methods used in this context. These DC formulations are easily extended to contact problems and applied to the determination of confined buckling shapes. To cite this article: P. Alart, S. Pagano, C. R. Mecanique 330 (2002) 819–824.     

A new formulation of immiscible compressible two-phase flow in porous media

A new formulation is proposed to describe immiscible compressible two-phase flow in porous media. The main feature of this formulation is the introduction of a global pressure. The resulting equations are written in a fractional flow formulation and lead to a coupled system which consists of a nonlinear parabolic (the global pressure equation) and a nonlinear diffusion–convection one (the saturation equation) which can be efficiently solved numerically. To cite this article: B. Amaziane, M. Jurak, C. R. Mecanique 336 (2008).     

Direct simulation of the motion of neutrally buoyant balls in a three-dimensional Poiseuille flow

In a previous article the authors introduced a Lagrange multiplier based fictitious domain method. Their goal in the present article is to apply a generalization of the above method to: (i) the numerical simulation of the motion of neutrally buoyant particles in a three-dimensional Poiseuille flow; (ii) study – via direct numerical simulations – the migration of neutrally buoyant balls in the tube Poiseuille flow of an incompressible Newtonian viscous fluid. Simulations made with one and several particles show that, as expected, the Segré–Silberberg effect takes place. To cite this article: T.-W. Pan, R. Glowinski, C. R. Mecanique 333 (2005).     

Vérification expérimentale de la relation de réciprocité d'Onsager pour l'électro-osmose et l'électro-filtration dans une argile naturelleExperimental verification of the Onsager's reciprocal relations for electro-osmosis and electro-filtration phenomena on a saturated clay

This Note is devoted to the experimental verification of the Onsager's reciprocal relations in the particular case of electro-osmosis and electro-filtration. A special set up has been designed to carry out the measurements of both the electro-osmotic permeability and the streaming potential. This has been performed by using a natural material i.e., saturated kaolinite. To cite this article: K. Beddiar et al., C. R. Mecanique 330 (2002) 893–898.     

Coupling tabulated chemistry with Large Eddy Simulation of turbulent reactive flows

A new modeling strategy is developed to introduce tabulated chemistry methods in the LES of turbulent premixed combustion. The objective is to recover the correct laminar flame propagation speed of the filtered flame front when the subgrid scale turbulence vanishes. The filtered flame structure is mapped by 1D filtered laminar premixed flames. Closure of the filtered progress variable and the energy balance equations are carefully addressed. The methodology is applied to 1D and 2D filtered laminar flames. These computations show the capability of the model to recover the laminar flame speed and the correct chemical structure when the flame wrinkling is completely resolved. The model is then extended to turbulent combustion regimes by introducing subgrid scale wrinkling effects on the flame front propagation. Finally, the LES of a 3D turbulent premixed flame is performed. To cite this article: R. Vicquelin et al., C. R. Mecanique 337 (2009).     

Modèles discrets de structures tissées : Analyse de stabilité et de drapéDiscrete models of woven structures: stability and draping analysis

A discrete model of a woven fabric structure is established, whereby nodes endowed with a mass and a rotational rigidity are connected by rigid bars to form a two-dimensional truss. The set of four bars that delineate a quadrilateral area is further endowed with a torsion deformation mode. The kinematics of the truss reproduces the large rotations and displacements encountered for real tissues. The equilibrium shape of such a structure is obtained as the minimum of its total potential energy versus the whole set of kinematic translational and rotational variables, accounting for eventual kinematic constraints due to contact with a rigid surface by the Lagrange multipliers method. A stability analysis is conducted, and the potentiality of the model is illustrated by fabric draping simulations. To cite this article: B. Ben Boubaker et al., C. R. Mecanique 330 (2002) 871–877.