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.

[1]  J. Alexander,et al.  Theory and Methods: Critical Essays in Human Geography , 2008 .

[2]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[3]  Joël Piraux,et al.  Numerical modeling of elastic waves across imperfect contacts , 2005, SIAM J. Sci. Comput..

[4]  Joël Piraux,et al.  A new interface method for hyperbolic problems with discontinuous coefficients: one-dimensional acoustic example , 2001 .

[5]  Johan O. A. Robertsson,et al.  A numerical free-surface condition for elastic/viscoelastic finite-difference modeling in the presence of topography , 1996 .

[6]  José M. Carcione,et al.  Seismic modeling in viscoelastic media , 1993 .

[7]  Chrysoula Tsogka,et al.  A TIME DOMAIN METHOD FOR MODELING VISCOACOUSTIC WAVE PROPAGATION , 2004 .

[8]  D. Komatitsch,et al.  Introduction to the spectral element method for three-dimensional seismic wave propagation , 1999 .

[9]  Ronald Fedkiw,et al.  The immersed interface method. Numerical solutions of PDEs involving interfaces and irregular domains , 2007, Math. Comput..

[10]  George A. McMechan,et al.  Composite memory variables for viscoelastic synthetic seismograms , 1995 .

[11]  Claus-Dieter Munz,et al.  Lax–Wendroff-type schemes of arbitrary order in several space dimensions , 2006 .

[12]  R. LeVeque.  Finite Volume Methods for Hyperbolic Problems: Characteristics and Riemann Problems for Linear Hyperbolic Equations , 2002 .

[13]  Einar Kjartansson,et al.  Attenuation of seismic waves in rocks and applications in energy exploration , 1979 .

[14]  E. Reiss,et al.  Reflection of Plane Viscoelastic Waves from Plane Boundaries , 1966 .

[15]  R. Ruth,et al.  Fourth-order symplectic integration , 1990 .

[16]  R. Borcherdt.  Reflection—refraction of general P-and type-I S-waves in elastic and anelastic solids , 1982 .

[17]  B. Gustafsson.  The convergence rate for difference approximations to mixed initial boundary value problems , 1975 .

[18]  Michael Dumbser,et al.  Fast high order ADER schemes for linear hyperbolic equations , 2004 .

[19]  M. Dumbser,et al.  An arbitrary high-order Discontinuous Galerkin method for elastic waves on unstructured meshes — III. Viscoelastic attenuation , 2007 .

[20]  G. M.,et al.  A Treatise on the Mathematical Theory of Elasticity , 1906, Nature.

[21]  P. Joly,et al.  A mixed finite element approach for viscoelastic wave propagation , 2005 .

[22]  Paul G. Richards,et al.  Quantitative Seismology: Theory and Methods , 1980 .

[23]  L. Knizhnerman,et al.  Finite-difference modeling of viscoelastic materials with quality factors of arbitrary magnitude , 2004 .

[24]  Zhilin Li,et al.  The Immersed Interface Method: Numerical Solutions of PDEs Involving Interfaces and Irregular Domains (Frontiers in Applied Mathematics) , 2006 .

[25]  M. Korn,et al.  Incorporation of attenuation into time-domain computations of seismic wave fields , 1987 .

[26]  José M. Carcione,et al.  3-D Wave simulation in anelastic media using the Kelvin-Voigt constitutive equation , 2004 .

[27]  L. Wennerberg.  Snell's law for viscoelastic materials , 1985 .

[28]  Bruno Lombard,et al.  Modélisation numérique de la propagation et de la diffraction d'ondes mécaniques , 2010 .

[29]  H. Cooper.  Reflection and Transmission of Oblique Plane Waves at a Plane Interface between Viscoelastic Media , 1967 .

[30]  Joël Piraux,et al.  Numerical treatment of two-dimensional interfaces for acoustic and elastic waves , 2004 .

[31]  Joakim O. Blanch,et al.  Modeling of a constant Q; methodology and algorithm for an efficient and optimally inexpensive viscoelastic technique , 1995 .

[32]  J. Bernard Minster,et al.  Numerical simulation of attenuated wavefields using a Padé approximant method , 1984 .

[33]  Joakim O. Blanch,et al.  Viscoelastic finite-difference modeling , 1994 .

[34]  José M. Carcione,et al.  Wave Fields in Real Media: Wave Propagation in Anisotropic, Anelastic and Porous Media , 2011 .

[35]  R. Borcherdt.  Reflection and refraction of type-II S waves in elastic and anelastic media , 1977, Bulletin of the Seismological Society of America.

[36]  Mathieu Chekroun,et al.  Comparison between a multiple scattering method and direct numerical simulations for elastic wave propagation in concrete , 2009 .

[37]  Peter Moczo,et al.  On the rheological models used for time‐domain methods of seismic wave propagation , 2005 .

[38]  F. J. Lockett.  The reflection and refraction of waves at an interface between viscoelastic materials , 1962 .