Wave simulation in 2D heterogeneous transversely isotropic porous media with fractional attenuation: a Cartesian grid approach

Abstract : 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/21/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.
Type de document :
Article dans une revue
Journal of Computational Physics, Elsevier, 2014, 275, pp.118-142
Liste complète des métadonnées

Littérature citée [56 références]  Voir  Masquer  Télécharger

Contributeur : Bruno Lombard <>
Soumis le : vendredi 4 juillet 2014 - 14:47:24
Dernière modification le : mardi 22 mars 2016 - 01:26:10
Document(s) archivé(s) le : samedi 4 octobre 2014 - 11:26:15


Fichiers produits par l'(les) auteur(s)


  • HAL Id : hal-00949686, version 2


Emilie Blanc, Guillaume Chiavassa, Bruno Lombard. Wave simulation in 2D heterogeneous transversely isotropic porous media with fractional attenuation: a Cartesian grid approach. Journal of Computational Physics, Elsevier, 2014, 275, pp.118-142. 〈hal-00949686v2〉



Consultations de
la notice


Téléchargements du document