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

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

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

Numerical analysis of a quasistatic piezoelectric problem with damage

The quasistatic evolution of the mechanical state of a piezoelectric body with damage is numerically studied in this paper. Both damage and piezoelectric effects are included into the model. The variational formulation leads to a coupled system composed of two linear variational equations for the displacement field and the electric potential, and a nonlinear parabolic variational equation for the damage field. The existence of a unique weak solution is stated. Then, a fully discrete scheme is introduced by using a finite element method to approximate the spatial variable and an Euler scheme to discretize the time derivatives. Error estimates are derived on the approximate solutions, from which the linear convergence of the algorithm is deduced under suitable regularity conditions. Finally, a two-dimensional example is presented to demonstrate the behaviour of the solution. To cite this article: J.R. Fernández 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).     

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

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

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

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

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

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

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.     

Homogénéisation de problèmes d''évolution dans une structure renforcée de fibresHomogenization of evolution problems in a fiber reinforced structure

We study the homogenization of evolution equations such as:

where the coefficient a is -periodic and takes very high values on a subset TΩ (fibered structure) of very small measure. We find a non-local effective equation deduced from a homogenized system of several equations. To cite this article: M. Bellieud, C. R. Mecanique 330 (2002) 843–848.     

Shape optimization of an accommodative intra-ocular lens

Cataract surgery consists in replacing the clouded or opacified crystalline lens by an Intra-Ocular Lens (IOL) having the same mean dioptrical power. Clear vision is then achieved at a given distance and glasses are needed in many situations. A new kind of IOL, potentially accommodative, is proposed. Its design is based on the deep understanding of the accommodation mechanism and on the mathematical modeling and the numerical simulation of the IOL's comportment in vivo. A preliminary version of this IOL is now commercialized by the company HumanOptics under the name ‘1CU’. In a second phase, shape optimization techniques equipped with strong mechanical and physiological constraints, are used to enhance the IOL performance and build a new design. To cite this article: F. Jouve, K. Hanna, C. R. Mecanique 333 (2005).     

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

On the Lagrange problem about the optimal form for circular hollow columnsSur le problème de Lagrange de la forme optimale d'une colonne circulaire aux parois minces

A new approach to the problem of finding the form of the strongest clamped circular column with thin walls of fixed volume and height is proposed. The same model describes also the form of the horizontal beam with rectangular vertical sections of a fixed height and variable widths. The existence and the uniqueness of the optimal column are proved here for the first time and the optimal form is found. To cite this article: Yu.V. Egorov, C. R. Mecanique 331 (2003).     

On the confinement of a viscous fluid by means of a feedback external field

In this work we consider a planar stationary flow of an incompressible viscous fluid in a semi-infinite strip governed by the standard Stokes system. We show how this fluid can be stopped at a finite distance from the entrance of the semi-infinite strip by means of a feedback source depending in a sublinear way on the velocity field. This localization effect is proved by reducing the problem to a non-linear biharmonic type one for which the localization of solutions is obtained through the application of an energy method, in the spirit of the monograph by S.N. Antontsev, J.I. Díaz and S.I. Shmarev (Energy Methods for Free Boundary Problems: Applications to Non-Linear PDEs and Fluid Mechanics, Birkäuser, Boston, 2002). Since the presence of the non-linear terms defined by the source is not standard in fluid mechanics literature, we give also some results about the existence and uniqueness of weak solutions for this problem. To cite this article: S.N. Antontsev et al., C. R. Mecanique 330 (2002) 797–802.     

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

On the cavitation problem and the influence of the surface energy in a viscoplastic medium

In this Note, the solution for spherically symmetric cavitation in a viscoplastic material is analysed. To ensure of the reality of the physical behaviour of the material, the problem is studied by considering a hollow sphere whose matrix obeys to a modified Bodner and Partom model. This local phenomenon is understood in the sense of the rapid growth of a pre-existing void and a particular attention is made to understand the influence of the surface energy on the critical dilative stress. To cite this article: F. Zaïri et al., C. R. Mecanique 333 (2005).