Source author record

Joël Piraux

Joël Piraux appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

physics.class-ph math.NA physics.comp-ph physics.geo-ph

Catalog footprint

What is connected

5works

4topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2012arXiv

Time-domain numerical simulations of multiple scattering to extract elastic effective wavenumbers

Elastic wave propagation is studied in a heterogeneous 2-D medium consisting of an elastic matrix containing randomly distributed circular elastic inclusions. The aim of this study is to determine the effective wavenumbers when the incident wavelength is similar to the radius of the inclusions. A purely numerical methodology is presented, with which the limitations usually associated with low scatterer concentrations can be avoided. The elastodynamic equations are integrated by a fourth-order time-domain numerical scheme. An immersed interface method is used to accurately discretize the interfaces on a Cartesian grid. The effective field is extracted from the simulated data, and signal-processing tools are used to obtain the complex effective wavenumbers. The numerical reference solution thus-obtained can be used to check the validity of multiple scattering analytical models. The method is applied to the case of concrete. A parametric study is performed on longitudinal and transverse incident plane waves at various scatterers concentrations. The phase velocities and attenuations determined numerically are compared with predictions obtained with multiple scattering models, such as the Independent Scattering Approximation model, the Waterman-Truell model, and the more recent Conoir-Norris model.

preprint2011arXiv

Numerical modeling of transient two-dimensional viscoelastic waves

This paper deals with the numerical modeling of transient mechanical waves in linear viscoelastic solids. Dissipation mechanisms are described using the generalized Zener model. No time convolutions are required thanks to the introduction of memory variables that satisfy local-in-time differential equations. By appropriately choosing the relaxation parameters, it is possible to accurately describe a large range of materials, such as solids with constant quality factors. The evolution equations satisfied by the velocity, the stress, and the memory variables are written in the form of a first-order system of PDEs with a source term. This system is solved by splitting it into two parts: the propagative part is discretized explicitly, using a fourth-order ADER scheme on a Cartesian grid, and the diffusive part is then solved exactly. Jump conditions along the interfaces are discretized by applying an immersed interface method. Numerical experiments of wave propagation in viscoelastic and fluid media show the efficiency of this numerical modeling for dealing with challenging problems, such as multiple scattering configurations.

preprint2008arXiv

Numerical modeling of 1-D transient poroelastic waves in the low-frequency range

Propagation of transient mechanical waves in porous media is numerically investigated in 1D. The framework is the linear Biot's model with frequency-independant coefficients. The coexistence of a propagating fast wave and a diffusive slow wave makes numerical modeling tricky. A method combining three numerical tools is proposed: a fourth-order ADER scheme with time-splitting to deal with the time-marching, a space-time mesh refinement to account for the small-scale evolution of the slow wave, and an interface method to enforce the jump conditions at interfaces. Comparisons with analytical solutions confirm the validity of this approach.

preprint2006arXiv

Modeling 1-D elastic P-waves in a fractured rock with hyperbolic jump conditions

The propagation of elastic waves in a fractured rock is investigated, both theoretically and numerically. Outside the fractures, the propagation of compressional waves is described in the simple framework of one-dimensional linear elastodynamics. The focus here is on the interactions between the waves and fractures: for this purpose, the mechanical behavior of the fractures is modeled using nonlinear jump conditions deduced from the Bandis-Barton model classicaly used in geomechanics. Well-posedness of the initial-boundary value problem thus obtained is proved. Numerical modeling is performed by coupling a time-domain finite-difference scheme with an interface method accounting for the jump conditions. The numerical experiments show the effects of contact nonlinearities. The harmonics generated may provide a non-destructive means of evaluating the mechanical properties of fractures.

preprint2005arXiv

Numerical modeling of elastic waves across imperfect contacts

A numerical method is described for studying how elastic waves interact with imperfect contacts such as fractures or glue layers existing between elastic solids. These contacts have been classicaly modeled by interfaces, using a simple rheological model consisting of a combination of normal and tangential linear springs and masses. The jump conditions satisfied by the elastic fields along the interfaces are called the "spring-mass conditions". By tuning the stiffness and mass values, it is possible to model various degrees of contact, from perfect bonding to stress-free surfaces. The conservation laws satisfied outside the interfaces are integrated using classical finite-difference schemes. The key problem arising here is how to discretize the spring-mass conditions, and how to insert them into a finite-difference scheme: this was the aim of the present paper. For this purpose, we adapted an interface method previously developed for use with perfect contacts [J. Comput. Phys. 195 (2004) 90-116]. This numerical method also describes closely the geometry of arbitrarily-shaped interfaces on a uniform Cartesian grid, at negligible extra computational cost. Comparisons with original analytical solutions show the efficiency of this approach.