Time-domain numerical modeling of poroelastic waves : the Biot-JKD model with fractional derivatives

• Auteur / Autrice : Emilie Blanc
• Direction : Bruno Lombard, Guillaume Chiavassa
• Type : Thèse de doctorat
• Discipline(s) : Acoustique
• Date : Soutenance le 05/12/2013
• Etablissement(s) : Aix-Marseille
• Ecole(s) doctorale(s) : École doctorale Sciences pour l'Ingénieur : Mécanique, Physique, Micro et Nanoélectronique (Marseille)
• Jury : Président / Présidente : Vincent Pagneux
• Jury : Examinateurs / Examinatrices : Boris Gralak, Pierre Suquet
• Jury : Rapporteur / Rapporteuse : Andrew Norris, Thierry Goudon

Mots clés

Mots clés contrôlés

Modèles mathématiques

Éléments finis, Méthode des

Porosité

Ondes -- Propagation

Mots clés libres

Milieux poreux

Ondes élastiques

Modèle de Biot-JKD

Dérivées fractionnaires

Méthode de splitting

Méthodes de différences finies

Méthode d'interface immergée

Porous media

Elastic waves

Biot-JKD model

Fractional derivatives

Time splitting method

Finite difference methods

Immersed interface method

Résumé

Une modélisation numérique des ondes poroélastiques, décrites par le modèle de Biot, est proposée dans le domaine temporel. La dissipation visqueuse à l'intérieur des pores est décrite par le modèle de perméabilité dynamique de Johnson-Koplik-Dashen (JKD). Certains coefficients du modèle de Biot-JKD sont proportionnels à la racine carrée de la fréquence, introduisant dans le domaine temporel des dérivées fractionnaires décalées d'ordre 1/2, revenant à un produit de convolution. Basé sur une représentation diffusive, le produit de convolution est remplacé par un nombre fini de variables de mémoire satisfaisant une équation différentielle ordinaire locale en temps, menant au modèle de Biot-DA (diffusive approximation). Les propriétés des deux modèles sont analysées : hyperbolicité, décroissance de l'énergie, dispersion. On montre que la meilleure méthode de détermination des coefficients de l'approximation diffusive - quadratures de Gauss, optimisation linéaire ou non-linéaire sur la plage de fréquence d'intérêt - est l'optimisation non-linéaire. Une méthode de splitting est utilisée numériquement : la partie propagative est discrétisée par un schéma aux différences finies ADER d'ordre 4, et la partie diffusive est intégrée exactement. Les conditions de saut aux interfaces sont discrétisées avec une méthode d'interface immergée. Des simulations numériques sont présentées pour des milieux isotropes et isotropes transverses. Des comparaisons avec des solutions analytiques montrent l'efficacité et la précision de cette approche. Des simulations numériques en milieux complexes sont réalisées : influence de la porosité d'os spongieux, diffusion multiple en milieu aléatoire.

A time-domain numerical modeling of Biot poroelastic waves is proposed. The viscous dissipation in the pores is described using the dynamic permeability model of 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 product is replaced by a finite number of memory variables that satisfy local-in-time ordinary differential equations, resulting in the Biot-DA (diffusive approximation). The properties of the two models are analyzed: hyperbolicity, decrease of energy, dispersion. To determine the coefficients of the diffusive approximation, different methods of quadrature are analyzed: Gaussian quadratures, linear or nonlinear optimization procedures in the frequency range of interest. The nonlinear optimization is shown to be the best way of determination. A splitting strategy is applied numerically: the propagative part is discretized using a fourth-order ADER scheme on a Cartesian grid, and the diffusive part is solved exactly. An immersed interface method is implemented to discretize the jump conditions at interfaces. Numerical experiments are presented for isotropic and transversely isotropic media. Comparisons with analytical solutions show the efficiency and the accuracy of this approach. Some numerical experiments are performed in complex media: influence of the porosity of a cancellous bone, multiple scattering across a set of random scatterers.