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

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

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

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

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

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

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.     

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

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

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

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

Transition of unsteady flows of evaporation to steady state

We investigate the half-space problem of evaporation and condensation in the scope of discrete kinetic theory. Exact solutions are found to the boundary value problem and the initial boundary value problems of the flow in the half space for a discrete velocity model. The results are used to analyze the transition of the unsteady solutions towards steady states. To cite this article: A. d'Almeida, C. R. Mecanique 336 (2008).     

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

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

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

Linear bending of star-like pyramids

A family of polyhedra possessing unusual deformation properties is found. On one hand, models of these polyhedra admit free continuous large reversible bending without visible distortions of the material. On the other hand, the polyhedra themselves are rigid and do not admit continuous bending in the sense of the Cauchy definition. The found polyhedra are called model flexors in order to distiguish them from theoretical flexsors of Connelly. Bendings of the models are asymptotically exactly approximated by linear bendings of polyhedra. They represent a nonrigid, soft or slow, loss of stability that corresponds to the loss of stability in small accordingly to Euler. This new phenomenon in mechanics of deformable solid bodies may be considered as an original geometric machine of catastrophe. To cite this article: A.D. Milka, C. R. Mecanique 331 (2003).     

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

Missing boundary data recovering for the Helmholtz problem

This Note is dedicated to the numerical treatment of the ill-posed Cauchy–Helmholtz problem. Resorting to the domain decomposition tools, these missing boundary data are rephrased through an ‘interfacial’ equation. This equation is solved via a preconditioned Richardson algorithm with dynamic relaxation. The efficiency of the proposed method is illustrated by some numerical experiments. To cite this article: R. Ben Fatma et al., C. R. Mecanique 335 (2007).     

Extended three-dimensional digital image correlation (X3D-DIC)

A correlation algorithm is proposed to measure full three-dimensional displacement fields in a three-dimensional domain. The chosen kinematic basis for this measurement is based on continuous finite-element shape functions. It is furthermore proposed to account for the presence of strong discontinuities, similarly to extended finite element schemes, with a suited enrichment of the kinematics with discontinuities supported by a (crack) surface. An optimization of the surface geometry is proposed based on correlation residuals. The procedure is applied to analyze one loading cycle of a fatigue-cracked nodular graphite cast iron sample by using computed tomography pictures. Subvoxel crack openings are revealed and measured. To cite this article: J. Réthoré et al., C. R. Mecanique 336 (2008).     

Chaos synchronization on the N-torus and cryptography

A class of chaotic dynamical systems on the N-dimensional torus is proposed for masking some information in secure communications. The information is then recovered thanks to a chaos synchronization process. To cite this article: L. Rosier et al., C. R. Mecanique 332 (2004).

Résumé

Nous proposons une classe de systèmes chaotiques sur le tore N-dimensionnel pour masquer une information à transmettre dans une communication sécurisée. Cette information est ensuite reconstruite à l'aide d'un mécanisme de synchronisation du chaos. Pour citer cet article : L. Rosier et al., C. R. Mecanique 332 (2004).