Bruno Lombard

Emilie Blanc, Guillaume Chiavassa, Bruno Lombard

A time-domain numerical modeling of Biot poroelastic waves is presented. The viscous dissipation occurring in the pores is described using the dynamic permeability model developed by Johnson-Koplik-Dashen (JKD). Some of the coefficients in the Biot-JKD model are proportional to the square root of the frequency: in the time-domain, these coefficients introduce order 1/2 shifted fractional derivatives involving a convolution product. Based on a diffusive representation, the convolution kernel is replaced by a finite number of memory variables that satisfy local-in-time ordinary differential equations. Thanks to the dispersion relation, the coefficients in the diffusive representation are obtained by performing an optimization procedure in the frequency range of interest. A splitting strategy is then applied numerically: the propagative part of Biot-JKD equations is discretized using a fourth-order ADER scheme on a Cartesian grid, whereas the diffusive part is solved exactly. Comparisons with analytical solutions show the efficiency and the accuracy of this approach.

Guillaume Chiavassa, Bruno Lombard

Numerical methods are developed to simulate the wave propagation in heterogeneous 2D fluid / poroelastic media. Wave propagation is described by the usual acoustics equations (in the fluid medium) and by the low-frequency Biot's equations (in the porous medium). Interface conditions are introduced to model various hydraulic contacts between the two media: open pores, sealed pores, and imperfect pores. Well-possedness of the initial-boundary value problem is proven. Cartesian grid numerical methods previously developed in porous heterogeneous media are adapted to the present context: a fourth-order ADER scheme with Strang splitting for time-marching; a space-time mesh-refinement to capture the slow compressional wave predicted by Biot's theory; and an immersed interface method to discretize the interface conditions and to introduce a subcell resolution. Numerical experiments and comparisons with exact solutions are proposed for the three types of interface conditions, demonstrating the accuracy of the approach.

Guillaume Chiavassa, Bruno Lombard

This paper deals with the numerical modeling of wave propagation in porous media described by Biot's theory. The viscous efforts between the fluid and the elastic skeleton are assumed to be a linear function of the relative velocity, which is valid in the low-frequency range. The coexistence of propagating fast compressional wave and shear wave, and of a diffusive slow compressional wave, makes numerical modeling tricky. To avoid restrictions on the time step, the Biot's system is splitted into two parts: the propagative part is discretized by a fourth-order ADER scheme, while the diffusive part is solved analytically. Near the material interfaces, a space-time mesh refinement is implemented to capture the small spatial scales related to the slow compressional wave. The jump conditions along the interfaces are discretized by an immersed interface method. Numerical experiments and comparisons with exact solutions confirm the accuracy of the numerical modeling. The efficiency of the approach is illustrated by simulations of multiple scattering.

Bruno Lombard, Joël Piraux

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.

Bruno Lombard, Denis Matignon

A fractional time derivative is introduced into the Burger's equation to model losses of nonlinear waves. This term amounts to a time convolution product, which greatly penalizes the numerical modeling. A diffusive representation of the fractional derivative is adopted here, replacing this nonlocal operator by a continuum of memory variables that satisfy local-in-time ordinary differential equations. Then a quadrature formula yields a system of local partial differential equations, well-suited to numerical integration. The determination of the quadrature coefficients is crucial to ensure both the well-posedness of the system and the computational efficiency of the diffusive approximation. For this purpose, optimization with constraint is shown to be a very efficient strategy. Strang splitting is used to solve successively the hyperbolic part by a shock-capturing scheme, and the diffusive part exactly. Numerical experiments are proposed to assess the efficiency of the numerical modeling, and to illustrate the effect of the fractional attenuation on the wave propagation.

Guillaume Chiavassa, Bruno Lombard, Joël Piraux

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.

Emilie Blanc, Guillaume Chiavassa, Bruno Lombard

An explicit finite-difference scheme is presented for solving the two-dimensional Biot equations of poroelasticity across the full range of frequencies. The key difficulty is to discretize the Johnson-Koplik-Dashen (JKD) model which describes the viscous dissipations in the pores. Indeed, the time-domain version of Biot-JKD model involves order 1/2 shifted fractional derivatives which amounts to a time convolution product. To avoid storing the past values of the solution, a diffusive representation of fractional derivatives is used: the convolution kernel is replaced by a finite number of memory variables that satisfy local-in-time ordinary differential equations. The coefficients of the diffusive representation follow from an optimization procedure of the dispersion relation. Then, various methods of scientific computing are applied: the propagative part of the equations is discretized using a fourth-order ADER scheme, whereas the diffusive part is solved exactly. An immersed interface method is implemented to discretize the geometry on a Cartesian grid, and also to enforce the jump conditions at interfaces. Numerical experiments are proposed in various realistic configurations.

Emilie Blanc, Guillaume Chiavassa, Bruno Lombard

A time-domain numerical modeling of transversely isotropic Biot poroelastic waves is proposed in two dimensions. The viscous dissipation occurring in the pores is described using the dynamic permeability model developed by Johnson-Koplik-Dashen (JKD). Some of the coefficients in the Biot-JKD model are proportional to the square root of the frequency. In the time-domain, these coefficients introduce shifted fractional derivatives of order 1/2, involving a convolution product. Based on a diffusive representation, the convolution kernel is replaced by a finite number of memory variables that satisfy local-in-time ordinary differential equations, resulting in the Biot-DA (diffusive approximation) model. The properties of both the Biot-JKD and the Biot-DA model are analyzed: hyperbolicity, decrease of energy, dispersion. To determine the coefficients of the diffusive approximation, two approaches are analyzed: Gaussian quadratures and optimization methods in the frequency range of interest. The nonlinear optimization is shown to be the better way of determination. A splitting strategy is then applied to approximate numerically the Biot-DA equations. The propagative part is discretized using a fourth-order ADER scheme on a Cartesian grid, whereas the diffusive part is solved exactly. An immersed interface method is implemented to take into account heterogeneous media on a Cartesian grid and to discretize the jump conditions at interfaces. Numerical experiments are presented. Comparisons with analytical solutions show the efficiency and the accuracy of the approach, and some numerical experiments are performed to investigate wave phenomena in complex media, such as multiple scattering across a set of random scatterers.

Nicolas Favrie, Bruno Lombard, Cédric Payan

Heterogeneous materials, such as rocks and concrete, have a complex dynamics including hysteresis, nonlinear elasticity and viscoelasticity. It is very sensitive to microstructural changes and damage. The goal of this paper is to propose a physical model describing the longitudinal vibrations of this class of material, and to develop a numerical strategy for solving the evolution equations. The theory relies on the coupling between two processes with radically-different time scales: a fast process at the frequency of the excitation, governed by nonlinear elasticity and viscoelasticity; a slow process, governed by the evolution of defects. The evolution equations are written as a nonlinear hyperbolic system with relaxation. A time-domain numerical scheme is developed, based on a splitting strategy. The numerical simulations show qualitative agreement with the features observed experimentally by Dynamic Acousto-Elastic Testing.

Bruno Lombard, Joël Piraux

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.