Direct numerical simulation and global stability analysis of three-dimensional instabilities in a lid-driven cavity

The first bifurcation in a lid-driven cavity characterized by three-dimensional Taylor–Görtler-Like instabilities is investigated for a cubical cavity with spanwise periodic boundary conditions at Re=1000. The modes predicted by a global linear stability analysis are compared to the results of a direct numerical simulation. The amplification rate, and the shape of the three-dimensional perturbation fields from the direct numerical simulation are in very good agreement with the characteristics of the steady S1 mode from the stability analysis, showing that this mode dominates the other unstable unsteady modes. To cite this article: J. Chicheportiche et al., C. R. Mecanique 336 (2008).     

Asymptotic modeling of linearly piezoelectric slender rods

The piezoelectric thin plate modeling already derived by the authors is extended to rod-like structures. Two models corresponding to sensor or actuator behavior are obtained. The conditions of existence of non-local terms in the limit models are discussed. To cite this article: T. Weller, C. Licht, C. R. Mecanique 336 (2008).     

A two grid algorithm based on perturbation and homotopy methods

In this paper, we propose a new class of bi-grid algorithm to solve large scale linear algebraic equations. This method is based on homotopy, perturbation technique and Padé approximants. To cite this article: R. El Mokhtari et al., C. R. Mecanique 330 (2002) 825–830.     

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).     

Homogenization of an elasticity problem in domains with a net of slender bars near surface

We consider an elasticity problem in a domain Ω⁽⁾=ΩF⁽⁾, where Ω is an open bounded domain in R³, F⁽⁾ is a connected nonperiodic set in Ω like a net of slender bars, and is a parameter characterizing the microstructure of the domain. We consider the case of a surface distribution of the set F⁽⁾, i.e., for sufficiently small , the set F⁽⁾ is concentrated in arbitrary small neighbourhood of a surface Γ. Under a hypothesis on the asymptotic behaviour of the energy functional, we obtain the macroscopic (homogenized) model. To cite this article: M. Goncharenko, L. Pankratov, C. R. Mecanique 331 (2003).     

Weighted Korn inequalities for thin-walled elastic structures

Appropriate weighted norms in H¹ are presented such that the Korn type inequality is asymptotically sharp with respect to relative thickness and stiffness of the elastic plates. The weights depend crucially on the geometric structure of the plates' junction. To cite this article: O.V. Izotova et al., C. R. Mecanique 334 (2006).     

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).     

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).     

Resolution of linear viscoelastic equations in the frequency domain using real Helmholtz boundary integral equations

Boundary integral equations are well suitable for the analysis of seismic waves propagation in unbounded domains. Formulations in elastodynamics are well developed. In contrast, for the dynamic analysis of viscoelastic media, there are very seldom formulations by boundary integral equations. In this Note, we propose a new and simple formulation of time harmonic viscoelasticity with the Zener model, which reduces to classical elastodynamics if a compatibility condition is satisfied by boundary conditions. Intermediate variables which satisfy the classical elastodynamic equations are introduced. It makes it possible to utilize existing numerical tools of time harmonic elastodynamics. To cite this article: S. Chaillat, H.D. Bui, C. R. Mecanique 335 (2007).     

Use of Faraday instabilities to enhance fuel pulverisation in air-blast atomisers

The atomization of liquids into a spray is an important process in many industrial applications and particularly in the aero-engine sector. Conventional air-blast injectors in aircraft engines today use aerodynamic shearing effects to atomize the liquid fuel. However, at operating conditions where the air velocity is below 30 m/s (such as ground start and high altitude restart) the atomization quality is poor. Consequently combustion is less efficient with high pollutant emissions. The objective of this study is to validate a new concept of injector which couples the shearing effects with the principle of ultrasonic atomization. The latter consists of using piezoelectric actuators to generate the oscillations of a wall in contact with the liquid film. This excitation perpendicular to the liquid film surface creates Faraday instabilities at the liquid/air interface. Amplitudes higher than a defined threshold value induce the break-up of ligaments and the formation of droplets. To cite this article: M. Boukra et al., C. R. Mecanique 337 (2009).     

A criterion of the continuous spectrum for elasticity and other self-adjoint systems on sharp peak-shaped domains

The spectra of the elasticity and piezo-electricity systems for a solid with a sharp peak point on the boundary, which is free of traction, are not discrete. An algebraic criterion of non-empty continuous spectrum is found for the Neumann problem for rather arbitrary formally self-adjoint elliptic systems of second-order differential equations on a sharp peak-shaped domain. To cite this article: S.A. Nazarov, C. R. Mecanique 335 (2007).     

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).     

Fluid structure interaction problems in large deformation

The present article deals with the simulation of fluid structure interaction problems in large deformation, and discusses two aspects of their numerical solution: (i) the derivation of energy conserving time integration schemes in presence of fluid structure coupling, moving grids, and nonlinear kinematic constraints such as incompressibility and contact, (ii) the introduction of adequate preconditioners efficiently chaining local fluid and structure solvers. Solutions are proposed, analyzed and tested using nonlinear energy correcting terms, and added mass based Dirichlet Neumann preconditioners. Numerical applications include nonlinear impact problems in elastodynamics and blood flows predictions within flexible arteries. To cite this article: P. Le Tallec et al., C. R. Mecanique 333 (2005).     

Analytical development of disturbed matrix eigenvalue problem applied to mixed convection stability analysis in Darcy media

This work consists in evaluating algebraically and numerically the influence of a disturbance on the spectral values of a diagonalizable matrix. Thus, two approaches will be possible; to use the theorem of disturbances of a matrix depending on a parameter, due to Lidskii and primarily based on the structure of Jordan of the no disturbed matrix. The second approach consists in factorizing the matrix system, and then carrying out a numerical calculation of the roots of the disturbances matrix characteristic polynomial. This problem can be a standard model in the equations of the continuous media mechanics. During this work, we chose to use the second approach and in order to illustrate the application, we choose the Rayleigh–Bénard problem in Darcy media, disturbed by a filtering through flow. The matrix form of the problem is calculated starting from a linear stability analysis by a finite elements method. We show that it is possible to break up the general phenomenon into other elementary ones described respectively by a disturbed matrix and a disturbance. A good agreement between the two methods was seen. To cite this article: H.B. Hamed, R. Bennacer, C. R. Mecanique 336 (2008).     

Modélisation du comportement thermique d'un outil de fraisage : approche par identification de système non entierThermal modelling of a milling tool: a noninteger system identification approach

In this work we present an experimental apparatus devoted to the thermal characterisation of a milling tool. The experimental device used thermistors, one for each insert. Each thermistor is located at a point in the tool close to the tip of the insert. The heat flux in each insert is expressed according to the temperature at the sensor from a non-integer model. The parameters of the model are identified from transient evolutions measurements of the temperature on the sensor and on the cutting edge. An application shows the difference in the behaviour of each insert during machining from the estimated heat fluxes. To cite this article: J.-L. Battaglia et al., C. R. Mecanique 330 (2002) 857–864.     

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.     

Asymptotics of eigenfrequencies of an elastic body with a heavy and hard peak-shaped inclusion

A heavy and hard peak-shaped inclusion in an elastic body provokes to concentration of eigenvalues in the low-frequency range of the spectrum and localization of the corresponding eigenmodes near the peak tip. To cite this article: S.A. Nazarov, C. R. Mecanique 335 (2007).     

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).     

Random-field model for the elasticity tensor of anisotropic random media

This Note deals with the construction of a non-Gaussian positive definite matrix-valued random field whose mathematical properties allow the fourth-order elasticity tensor of random non homogeneous anisotropic three dimensional elastic media to be modelled. If the usual parametric probabilistic approach was used, then 21 mutually dependent random fields should be modelled and identified by using experimental data. Such an approach would be very difficult because the systems of the marginal probability distributions of these random fields have to be identified due to the fact that, for a boundary value problem, the displacement field of the random medium is a non-linear mapping of the random elasticity tensor. The theory presented in this paper allows such a probabilistic model of the fourth-order elasticity tensor field to be constructed and depends only of four scalar parameters: three spatial correlation lengths and one parameter allowing the level of the random fluctuations to be controlled. To cite this article: C. Soize, C. R. Mecanique 332 (2004).

Résumé

On présente la construction d'un champ aléatoire à valeurs dans les matrices définies positives dont les propriétés mathématiques permettent de modéliser le tenseur d'élasticité du quatrième ordre des mileux élastiques anisotropes tridimensionnels aléatoires. Si l'approche probabiliste paramétrique usuelle était utilisée, alors il serait nécessaire de modéliser et d'identifier à l'aide de données expérimentales 21 champs aléatoires mutuellement dépendants. Une telle approche serait très difficile de part le fait que le système de lois marginales de ces champs aléatoires doit être identifié parce que, pour un problème aux limites, le champ de déplacement est une transformation non linéaire du tenseur d'élasticité. La théorie présentée dans ce papier permet de construire une modélisation probabiliste du champ de tenseur d'élasticité qui ne dépend que de quatre paramètres scalaires : trois échelles de corrélation spatiale et un paramètre permettant de contrôler le niveau des fluctuations aléatoires. Pour citer cet article : C. Soize, C. R. Mecanique 332 (2004).