An FC-based spectral solver for elastodynamic problems in general three-dimensional domains

This paper presents a spectral numerical algorithm for the solution of elastodynamics problems in general three-dimensional domains. Based on a recently introduced "Fourier continuation" (FC) methodology for accurate Fourier expansion of non-periodic functions, the proposed approach possesses a number of appealing properties: it yields results that are essentially free of dispersion errors, it entails mild CFL constraints, it runs at a cost that scales linearly with the discretization sizes, and it lends itself easily to efficient parallelization in distributed-memory computing clusters. The proposed algorithm is demonstrated in this paper by means of a number of applications to problems of isotropic elastodynamics that arise in the fields of materials science and seismology. These examples suggest that the new approach can yield solutions within a prescribed error tolerance by means of significantly smaller discretizations and shorter computing times than those required by other methods.

[1]  P. Roache.  Code Verification by the Method of Manufactured Solutions , 2002 .

[2]  Mark Lyon,et al.  High-order unconditionally stable FC-AD solvers for general smooth domains I. Basic elements , 2010, J. Comput. Phys..

[3]  W. Press,et al.  Numerical Recipes: The Art of Scientific Computing , 1987 .

[4]  Oscar P. Bruno,et al.  Higher-Order Linear-Time Unconditionally Stable Alternating Direction Implicit Methods for Nonlinear Convection-Diffusion Partial Differential Equation Systems , 2014 .

[5]  S. Spekreijse.  Elliptic grid generation based on Laplace equations and algebraic transformations , 1995 .

[6]  David I. Gottlieb,et al.  On the Removal of Boundary Errors Caused by Runge-Kutta Integration of Nonlinear Partial Differential Equations , 1994, SIAM J. Sci. Comput..

[7]  Olivier Coutant,et al.  Seismic response of three-dimensional topographies using a time-domain boundary element method , 2000 .

[8]  Cristina Trillo,et al.  Transient elastic wave propagation and scattering in plates: comparison between pulsed TV-holography measurements and finite element method predictions , 2013 .

[9]  Nathan Albin,et al.  A spectral FC solver for the compressible Navier-Stokes equations in general domains I: Explicit time-stepping , 2011, J. Comput. Phys..

[10]  Tim Elling.  GPU-accelerated Fourier-continuation solvers and physically exact computational boundary conditions for wave scattering problems , 2013 .

[11]  Oscar P. Bruno,et al.  Spatially Dispersionless, Unconditionally Stable FC–AD Solvers for Variable-Coefficient PDEs , 2012, Journal of Scientific Computing.

[12]  K. Graff.  Wave Motion in Elastic Solids , 1975 .

[13]  Cristina Trillo,et al.  Video Ultrasonics by Pulsed TV Holography: A New Capability for Non-Destructive Testing of Shell Structures , 2007 .

[14]  M. Farrashkhalvat,et al.  Basic Structured Grid Generation: With an introduction to unstructured grid generation , 2003 .

[15]  David I. Gottlieb,et al.  The Theoretical Accuracy of Runge-Kutta Time Discretizations for the Initial Boundary Value Problem: A Study of the Boundary Error , 1995, SIAM J. Sci. Comput..

[16]  D. Komatitsch,et al.  The spectral element method: An efficient tool to simulate the seismic response of 2D and 3D geological structures , 1998, Bulletin of the Seismological Society of America.

[17]  William D. Henshaw,et al.  Parallel computation of three-dimensional flows using overlapping grids with adaptive mesh refinement , 2008, J. Comput. Phys..

[18]  Stefan A. Sauter,et al.  Is the Pollution Effect of the FEM Avoidable for the Helmholtz Equation Considering High Wave Numbers? , 1997, SIAM Rev..

[19]  Nathan Albin,et al.  Fourier continuation methods for high-fidelity simulation of nonlinear acoustic beams. , 2012, The Journal of the Acoustical Society of America.

[20]  Ibrahim Saidu,et al.  A Simplified Derivation and Analysis of Fourth Order Runge Kutta Method , 2010 .

[21]  Jan S. Hesthaven,et al.  Multi-dimensional hybrid Fourier continuation-WENO solvers for conservation laws , 2013, J. Comput. Phys..

[22]  Cristina Trillo,et al.  Numerical modeling and measurement by pulsed television holography of ultrasonic displacement maps in plates with through-thickness defects , 2010 .

[23]  Christopher J. Roy,et al.  Verification of a Compressible CFD Code Using the Method of Manufactured Solutions , 2002 .

[24]  Zhiyong Si,et al.  Modified Characteristics Gauge–Uzawa Finite Element Method for Time Dependent Conduction–Convection Problems , 2014, J. Sci. Comput..

[25]  D. Pathria,et al.  The Correct Formulation of Intermediate Boundary Conditions for Runge-Kutta Time Integration of Initial Boundary Value Problems , 1997, SIAM J. Sci. Comput..

[26]  Lars-Erik Eriksson,et al.  Practical Three-Dimensional Mesh Generation Using Transfinite Interpolation , 1985 .

[27]  João Marcelo Vedovoto,et al.  Application of the method of manufactured solutions to the verification of a pressure-based finite-volume numerical scheme , 2011 .

[28]  William H. Press,et al.  Numerical Recipes 3rd Edition: The Art of Scientific Computing , 2007 .

[29]  Mark Lyon,et al.  High-order unconditionally stable FC-AD solvers for general smooth domains II. Elliptic, parabolic and hyperbolic PDEs; theoretical considerations , 2010, J. Comput. Phys..

[30]  Nathan Albin,et al.  Multi-domain Fourier-continuation/WENO hybrid solver for conservation laws , 2011, J. Comput. Phys..

[31]  Craig A. Schultz,et al.  Effect of three‐dimensional topography on seismic motion , 1996 .

[32]  D. Appelö,et al.  A stable finite difference method for the elastic wave equation on complex geometries with free surfaces , 2007 .

[33]  Oscar P. Bruno,et al.  Higher-order in time "quasi-unconditionally stable" ADI solvers for the compressible Navier-Stokes equations in 2D and 3D curvilinear domains , 2015, J. Comput. Phys..

[34]  Faisal Amlani,et al.  A New High-Order Fourier Continuation-Based Elasticity Solver for Complex Three-Dimensional Geometries , 2014 .