Adaptive finite elements for the steady free fall of a body in a Newtonian fluid

The numerical simulation of the free fall of a solid body in a viscous fluid is a challenging task since it requires computational domains which usually need to be several order of magnitude larger than the solid body in order to avoid the influence of artificial boundaries. Toward an optimal mesh design in that context, we propose a method based on the weighted a posteriori error estimation of the finite element approximation of the fluid/body motion. A key ingredient for the proposed approach is the reformulation of the conservation and kinetic equations in the solid frame as well as the implicit treatment of the hydrodynamic forces and torque acting on the solid body in the weak formulation. Information given by the solution of an adequate dual problem allows one to control the discretization error of given functionals. The analysis encompasses the control of the free fall velocity, the orientation of the body, the hydrodynamic force and torque on the body. Numerical experiments for the two dimensional sedimentation problem validate the method. To cite this article: V. Heuveline, C. R. Mecanique 333 (2005).     

A two-dimensional effective model describing fluid–structure interaction in blood flow: analysis, simulation and experimental validation

We derive a closed system of effective equations describing a time-dependent flow of a viscous incompressible Newtonian fluid through a long and narrow elastic tube. The 3D axially symmetric incompressible Navier–Stokes equations are used to model the flow. Two models are used to describe the tube wall: the linear membrane shell model and the linearly elastic membrane and the curved, linearly elastic Koiter shell model. We study the behavior of the coupled fluid–structure interaction problem in the limit when the ratio between the radius and the length of the tube, , tends to zero. We obtain the reduced equations that are of Biot type with memory. An interesting feature of the reduced equations is that the memory term explicitly captures the viscoelastic nature of the coupled problem. Our model provides significant improvement over the standard 1D approximations of the fluid–structure interaction problem, all of which assume an ad hoc closure assumption for the velocity profile. We performed experimental validation of the reduced model using a mock circulatory flow loop assembled at the Cardiovascular Research Laboratory at the Texas Heart Institute. Experimental results show excellent agreement with the numerically calculated solution. Major applications include blood flow through large human arteries. To cite this article: S. Čanić et al., C. R. Mecanique 333 (2005).     

Experimental observation of frequency doubling in a viscoelastic mixing layer

An experimental mixing layer in water at Reynolds number 440 is investigated. A colored viscoelastic solution is introduced in the shear layer before the roll-up dynamics. On the basis of flow visualization and local velocity measurements, it is found that compared to the Newtonian case, the roll-up process is affected by the non-Newtonian behavior of the viscoelastic solution. The effect consists of the appearance of secondary eddies in the mixing layer corresponding to the production of higher harmonics in the vorticity distribution. Consequently, there is a frequency doubling of the local velocity oscillations in the mixing layer. To cite this article: F. Sausset et al., C. R. Mecanique 332 (2004).

Résumé

Une étude expérimentale de la couche de mélange à un nombre de Reynolds de 440 est réalisée. Une solution viscoélastique colorée est introduite dans la couche cisaillée stationnaire juste avant la formation des tourbillons. En utilisant un mesure locale de la vitesse et une technique de visualisation par fluorescence, il est montré que la formation tourbillonnaire est affectée par le comportement non-Newtonien de la solution viscoélastique par rapport au cas Newtonien. Cet effet correspond à la formation de tourbillons secondaires produisant un harmonique supérieur dans la répartition spatiale de vorticité. En conséquence, un doublement de la fréquence des oscillations de la vitesse locale dans la couche de mélange est observé. Pour citer cet article : F. Sausset et al., C. R. Mecanique 332 (2004).     

Effective flow of a viscous liquid through a helical pipe

We study the flow of a viscous fluid through a pipe with helical shape parameterized with , where the small parameter stands for the distance between two coils of the helix. The pipe has small cross-section of size . Using the asymptotic analysis of the microscopic flow described by the Navier–Stokes system, with respect to the small parameter that tends to zero, we find the effective fluid flow described by an explicit formula of the Poisseuile type including a small distorsion due to the particular geometry of the pipe. To cite this article: E. Marušić-Paloka, I. Pažanin, C. R. Mecanique 332 (2004).

Résumé

On considère un écoulement dans un tube de section circulaire et de forme hélicoïdale paramétré par , où est la distance entre deux tours de la spirale. Le rayon de la section du tube est lui aussi supposé égal à . A partir de l'écoulement microscopique décrit par le système de Navier–Stokes et en utilisant l'analyse asymptotique par rapport à ce petit paramètre on obtient l'écoulemment effectif décrit par une formule explicite de type Poiseuille associée à une petite déviation due à la géometrie du tube. Pour citer cet article : E. Marušić-Paloka, I. Pažanin, C. R. Mecanique 332 (2004).     

On the slow gravity-driven migration of arbitrary clusters of small solid particles

A new approach is advocated to compute at a low cpu time cost the rigid-body motions of settling solid particles when inertial effects are negligible. In addition to the relevant boundary-integral equations, the numerical implementation and a few convincing benchmark tests we address two configurations of equivalent spheres and spheroids, i.e. that exhibit when isolated the same settling velocity. To cite this article: A. Sellier, C. R. Mecanique 332 (2004).

Résumé

On propose une approche originale pour déterminer le mouvement d'une assemblée de particules solides et de formes arbitraires soumise à l'action de la pesanteur dans l'approximation de Stokes. Outre les intégrales de frontière et la méthode numérique associées on présente quelques comparaisons et examine le cas de deux configurations de sphères et ellipsoides de révolution équivalents, c'est-à-dire dotés lorsqu'ils sont seuls de la même vitesse de sédimentation. Pour citer cet article : A. Sellier, C. R. Mecanique 332 (2004).     

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

Thermodynamical functions for a gas of point vortices

We formulate nonlinear integro-differential equation for the averaged collective Hamiltonian of a gas of interacting two-dimensional vortices, derive its analytical solution, and discuss the equilibrium, axially-symmetrical, probability distributions that are possible for such a model. We also theoretically prove that the probability distribution for a system of 2D point vortices takes a form similar to the Gibbs distribution, but point out that the physical fundamentals of such a system differ from the standard theory of interacting particles. Furthermore, we find thermodynamical functions for positive and negative “temperature” of the system, and point out that the states with positive “temperature” correspond to stationary bell-shape vortex distributions, while the states with negative “temperature” correspond to distributions localized near container walls. To cite this article: E. Bécu et al., 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).     

Large Eddy Simulation of a plane impinging jet

Large Eddy Simulations of a plane turbulent impinging jet have been carried out using the dynamic Smagorinsky model. The statistical results are first validated with the measurements from the literature: mean and turbulent quantities along the jet axis and at different vertical locations are presented. This study is completed by the analysis of the wall shear stress at the impingement wall. The effect of the jet Reynolds number (3000Re13500) on the kinematic development of the jet is also discussed. To cite this article: F. Beaubert, S. Viazzo, C. R. Mecanique 330 (2002) 803–810.     

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

A simple multiaxial fatigue criterion for metals

A new simple multiaxial high-cycle fatigue endurance criterion, suitable for situations where the convex hull associated with the stress path approximates well an ellipsoid, is proposed. It considers, as measures of fatigue solicitation: (i) a new definition for the equivalent shear stress amplitude; and (ii) the maximum principal stress along the stress history. Assessment of the resulting criterion for a wide range of in-phase and out-phase cyclic loads shows that it compares very well with experimental data published in the literature. To cite this article: C.A. Gonçalves et al., C. R. Mecanique 332 (2004).

Résumé

On propose un critère simple d'endurance à la fatigue polycyclique, applicable à des situations pour lesquelles l'envelope convexe associé à l'histoire des contraintes s'approche bien d'une ellipse. Le critère considère, comme mesures de solicitation à la fatigue : (i) une nouvelle définition de l'amplitude de contrainte de cisaillement ; et (ii) la contrainte principale maximale au cours de l'histoire de chargement. Pour citer cet article : C.A. Gonçalves et al., C. R. Mecanique 332 (2004).     

Impact of a bead on a rotating wall

We study experimentally the impact of a plastic bead on a rotating wall made of steel (velocity Ω; radial position x₀). The results show that the restitution coefficient is directly function of the impact velocity x₀Ω and is invariant by changing frame reference. The influence of the height of release of the particle on its angular velocity after impact is also studied. We observe an increase of the angular velocity with height followed by a saturation. We propose an interpretation for this evolution considering that the particle may roll without sliding during all the impact. This physical feature is not always taken into account in existing models of impact between rigid bodies. To cite this article: F. Rioual et al., C. R. Mecanique 336 (2008).     

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

The influence of crack-imbalance orientation and orbital evolution for an extended cracked Jeffcott rotor

Vibration peaks occurring at rational fractions of the fundamental rotating critical speed, here named Local Resonances, facilitate cracked shaft detection during machine shut-down. A modified Jeffcott-rotor on journal bearings accounting for gravity effects and oscillating around nontrivial equilibrium points is employed. Modal parameter selection allows this linear model to represent first mode characteristics of real machines. Orbit evolution and vibration patterns are analyzed, yielding useful results. Crack detection results indicate that, instead of 1x and 2x components, analysis of the remaining local resonances should have priority; this is due to crack-residual imbalance interaction and to 2x multiple induced origins. Therefore, local resonances and orbital evolution around 1/2, 1/3 and 1/4 of the critical speed are emphasized for various crack-imbalance orientations. To cite this article: J. Gómez-Mancilla et al., C. R. Mecanique 332 (2004).

Résumé

Les pics de vibration apparaissant au passage des fractions de la vitesse de critique de rotation des systèmes tournants, appelées résonances locales, facilitent la détection de fissures sur les machines. Dans cette étude, un modèle de rotor Jeffcott modifié avec une fissure tournante, comportant des coussinets et prenant en compte les effets de pesanteur et de balourd est présenté. Le choix modal des paramètres permet de représenter les caractéristiques liées au premier mode des machines tournantes usuelles. Les évolutions des vibrations et des orbites du système comportant une fissure sont analysées et permettent d'obtenir des résultats utiles pour la détection des fissures sur les machines tournantes. Ainsi, ces résultats indiquent que, en plus des composants 1x et 2x, l'analyse des autres résonances locales restantes doivent être regardées avec attention du fait de l'interaction possible entre les différentes orientations de la fissure et du balourd, et des origines multiples pouvant engendrer la présence des résonances 2x. Par conséquent les résonances et l'évolution des orbites obtenus autour de 1/2, 1/3 et 1/4 de la vitesse critique sont étudiées pour différentes variations d'angle entre le balourd et l'orientation de la fissure. Pour citer cet article : J. Gómez-Mancilla et al., C. R. Mecanique 332 (2004).     

Edge singularities and structure of the 3-D Williams expansion

The elastic solution in a vicinity of a re-entrant wedge can be described by a Williams like expansion in terms of powers of the distance to a point on the edge. This expansion has a particular structure due to the invariance of the problem by translation parallel to the edge. We show here that some terms, so-called primary solutions, derive directly from solutions to the 2-D corner problem posed in the orthogonal cross section of the domain. The others, baptized shadow functions, derive of the primary solutions by integration along the axis parallel to the edge. This 3-D Williams expansion is shown to be equivalent to the edge expansion proposed by Costabel et al. [M. Costabel, M. Dauge, Z. Yosibash, A quasidual function method for extracting edge stress intensity functions, SIAM J. Math. Anal. 35 (5) (2004) 1177–1202]. To cite this article: T. Apel et al., C. R. Mecanique 336 (2008).     

Continuous thermodynamics for droplet vaporization: Comparison between Gamma-PDF model and QMoM

The Continuous Thermodynamics Model (CTM) (Cotterman et al., 1985) is a suitable method to reduce computational cost of multi-component vaporization models. The droplet composition is described by a probability density function (PDF) rather than tens of components in the classical Discrete Component Model (DCM). In the first CTM method developed for this application, the PDF was assumed to be a Γ-function (Hallett, 2000), but some problems had appeared in the case of vapor condensation at the droplet surface (Harstadt et al., 2003). The method put forward in this article, the Quadrature Method of Moments (QMoM), enables one to avoid any assumption on the PDF mathematical form. Following Lage who has developed this method for phase equilibria (Lage, 2007), this article widens the scope of QMoM to the modelling of multi-component droplet vaporization. The different CTM approaches are presented in the first part and the results obtained for a vapor condensation test case are then compared and analysed to illustrate improvements made by QMoM. To cite this article: C. Laurent et al., C. R. Mecanique 337 (2009).     

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

Computing Green's function of elasticity in a half-plane with impedance boundary condition

This Note presents an effective and accurate method for numerical calculation of the Green's function G associated with the time harmonic elasticity system in a half-plane, where an impedance boundary condition is considered. The need to compute this function arises when studying wave propagation in underground mining and seismological engineering. To theoretically obtain this Green's function, we have drawn our inspiration from the paper by Durán et al. (2005), where the Green's function for the Helmholtz equation has been computed. The method consists in applying a partial Fourier transform, which allows an explicit calculation of the so-called spectral Green's function. In order to compute its inverse Fourier transform, we separate as a sum of two terms. The first is associated with the whole plane, whereas the second takes into account the half-plane and the boundary conditions. The first term corresponds to the Green's function of the well known time-harmonic elasticity system in (cf. J. Dompierre, Thesis). The second term is separated as a sum of three terms, where two of them contain singularities in the spectral variable (pseudo-poles and poles) and the other is regular and decreasing at infinity. The inverse Fourier transform of the singular terms are analytically computed, whereas the regular one is numerically obtained via an FFT algorithm. We present a numerical result. Moreover, we show that, under some conditions, a fourth additional slowness appears and which could produce a new surface wave. To cite this article: M. Durán et al., C. R. Mecanique 334 (2006).     

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

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