Time-harmonic elasticity with controllability and higher-order discretization methods

The time-harmonic solution of the linear elastic wave equation is needed for a variety of applications. The typical procedure for solving the time-harmonic elastic wave equation leads to difficulties solving large-scale indefinite linear systems. To avoid these difficulties, we consider the original time dependent equation with a method based on an exact controllability formulation. The main idea of this approach is to find initial conditions such that after one time-period, the solution and its time derivative coincide with the initial conditions. The wave equation is discretized in the space domain with spectral elements. The degrees of freedom associated with the basis functions are situated at the Gauss-Lobatto quadrature points of the elements, and the Gauss-Lobatto quadrature rule is used so that the mass matrix becomes diagonal. This method is combined with the second-order central finite difference or the fourth-order Runge-Kutta time discretization. As a consequence of these choices, only matrix-vector products are needed in time dependent simulation. This makes the controllability method computationally efficient.

[1]  J. Desanto.  Mathematical and numerical aspects of wave propagation , 1998 .

[2]  Patrick Joly,et al.  FICTITIOUS DOMAINS, MIXED FINITE ELEMENTS AND PERFECTLY MATCHED LAYERS FOR 2-D ELASTIC WAVE PROPAGATION , 2001 .

[3]  Ferdinand Kickinger,et al.  Algebraic Multi-grid for Discrete Elliptic Second-Order Problems , 1998 .

[4]  Gary Cohen.  Higher-Order Numerical Methods for Transient Wave Equations , 2001 .

[5]  Bruno Després,et al.  Using Plane Waves as Base Functions for Solving Time Harmonic Equations with the Ultra Weak Variational Formulation , 2003 .

[6]  Isaac Harari,et al.  Studies of FE/PML for exterior problems of time-harmonic elastic waves , 2004 .

[7]  Xiaobing Feng,et al.  A domain decomposition method for solving a Helmholtz-like problem in elasticity based on the Wilson nonconforming element , 1997 .

[8]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[9]  C. Farhat,et al.  The discontinuous enrichment method for elastic wave propagation in the medium‐frequency regime , 2006 .

[10]  M. Dumbser,et al.  An arbitrary high-order discontinuous Galerkin method for elastic waves on unstructured meshes - I. The two-dimensional isotropic case with external source terms , 2006 .

[11]  Marcus J. Grote,et al.  Dirichlet-to-Neumann map for three-dimensional elastic waves , 2003 .

[12]  A. Majda,et al.  Radiation boundary conditions for acoustic and elastic wave calculations , 1979 .

[13]  Ezio Faccioli,et al.  Spectral-domain decomposition methods for the solution of acoustic and elastic wave equations , 1996 .

[14]  Patrick Joly,et al.  Domain Decomposition Method for Harmonic Wave Propagation : A General Presentation , 2000 .

[15]  R. Glowinski,et al.  Controllability Methods for the Computation of Time-Periodic Solutions; Application to Scattering , 1998 .

[16]  Tuomo Rossi,et al.  Controllability method for the Helmholtz equation with higher-order discretizations , 2007, J. Comput. Phys..

[17]  Thomas A. Manteuffel,et al.  Algebraic multigrid for higher-order finite elements , 2005 .

[18]  S. Mönkölä,et al.  Controllability method for acoustic scattering with spectral elements , 2007 .

[19]  Roland Martin,et al.  WAVE PROPAGATION IN 2-D ELASTIC MEDIA USING A SPECTRAL ELEMENT METHOD WITH TRIANGLES AND QUADRANGLES , 2001 .

[20]  Isaac Harari,et al.  Improved finite element methods for elastic waves , 1998 .

[21]  A. J. Berkhout,et al.  Elastic Wave Field Extrapolation: Redatuming of Single- and Multi-Component Seismic Data , 1989 .

[22]  Erkki Heikkola,et al.  An algebraic multigrid based shifted-Laplacian preconditioner for the Helmholtz equation , 2007, J. Comput. Phys..

[23]  Alfio Quarteroni,et al.  Generalized Galerkin approximations of elastic waves with absorbing boundary conditions , 1998 .

[24]  Jari P. Kaipio,et al.  The Ultra-Weak Variational Formulation for Elastic Wave Problems , 2004, SIAM J. Sci. Comput..

[25]  G. Mur,et al.  The finite-element modeling of three-dimensional electromagnetic fields using edge and nodal elements , 1993 .

[26]  Jean-Pierre Vilotte,et al.  The Spectral Element Method For Elastic Wa V Eequations: Application to 2D And 3D Seismic Problems , 1998 .

[27]  David Bindel,et al.  Elastic PMLs for resonator anchor loss simulation , 2005 .

[28]  Jing Li,et al.  An iterative domain decomposition method for the solution of a class of indefinite problems in computational structural dynamics , 2005 .

[29]  R. Glowinski,et al.  Exact and approximate controllability for distributed parameter systems , 1995, Acta Numerica.

[30]  S. Timoshenko,et al.  Theory of elasticity , 1975 .

[31]  D Komatitsch,et al.  CASTILLO-COVARRUBIAS JM, SANCHEZ-SESMA FJ. THE SPECTRAL ELEMENT METHOD FOR ELASTIC WAVE EQUATIONS-APPLICATION TO 2-D AND 3-D SEISMIC PROBLEMS , 1999 .

[32]  Alfio Quarteroni,et al.  Numerical solution of linear elastic problems by spectral collocation methods , 1993 .