Hybrid Asynchronous Perfectly Matched Layer for seismic wave propagation in unbounded domains

Perfectly Matched Layer (PML) is recognized as a very effective tool for modeling unbounded domains. Nonetheless, the computation time required by the PML may be large, especially when an explicit time integration scheme is adopted for dealing with the wave propagation problem both in the domain of interest and in the PML medium. In this paper, it is proposed to investigate subdomain strategies enabling the appropriate time integration scheme in the PML with its own time step to be chosen, independently of the choice of the time scheme in the domain of interest. We focus on explicit time integrator in the physical subdomain (Central Difference scheme) associated with a fine time step satisfying the CFL stability criterion. The PML formulation proposed by Basu and Chopra (2004) [1] for 2D transient dynamics, has been coupled with the interior physical subdomain using the dual Schur approach proposed by Gravouil and Combescure (2001) [2]. Hybrid (implicit time integrator for the PML) asynchronous (multi time steps) PMLs have been derived. Their very good accuracy has been shown by considering the following numerical examples: Lambs test, loaded rigid strip footing on an half space and a layered half space. HighlightsSubdomain decomposition framework applied to a problem coupling a physical medium with a Perfectly Matched Layer (PML) medium.Coupling method based on Dual Schur approach.Derivation of Hybrid (different time integrators) Asynchronous (different time steps) PML.General framework for coupling complex PML formulations while conserving classical time integrators in other partitions.

[1]  Alain Combescure,et al.  Two FETI-based heterogeneous time step coupling methods for Newmark and α-schemes derived from the energy method , 2015 .

[2]  A. Combescure,et al.  Coupling subdomains with heterogeneous time integrators and incompatible time steps , 2009 .

[3]  Qing Huo Liu,et al.  PERFECTLY MATCHED LAYERS FOR ELASTODYNAMICS: A NEW ABSORBING BOUNDARY CONDITION , 1996 .

[4]  Shuo Ma,et al.  Modeling of the Perfectly Matched Layer Absorbing Boundaries and Intrinsic Attenuation in Explicit Finite-Element Methods , 2006 .

[5]  Loukas F. Kallivokas,et al.  A symmetric hybrid formulation for transient wave simulations in PML-truncated heterogeneous media , 2013 .

[6]  Alain Combescure,et al.  A numerical scheme to couple subdomains with different time-steps for predominantly linear transient analysis , 2002 .

[7]  Alain Combescure,et al.  External coupling software based on macro- and micro-time scales for explicit/implicit multi-time-step co-computations in structural dynamics , 2014 .

[8]  Alain Combescure,et al.  Multi-time-step explicit–implicit method for non-linear structural dynamics , 2001 .

[9]  Finite Difference Modeling of Elastic Wave Propagation , 2012 .

[10]  A. Majda,et al.  Absorbing boundary conditions for the numerical simulation of waves , 1977 .

[11]  Irini Djeran-Maigre,et al.  Design of an efficient multi‐directional explicit/implicit Rayleigh absorbing layer for seismic wave propagation in unbounded domain using a strong form formulation , 2016 .

[12]  Anthony Gravouil,et al.  Explicit/implicit multi-time step co-computations for blast analyses on a reinforced concrete frame structure , 2012 .

[13]  O. C. Zienkiewicz,et al.  Diffraction and refraction of surface waves using finite and infinite elements , 1977 .

[14]  Jean-Pierre Berenger,et al.  A perfectly matched layer for the absorption of electromagnetic waves , 1994 .

[15]  Prabhu Rajagopal,et al.  On the use of absorbing layers to simulate the propagation of elastic waves in unbounded isotropic media using commercially available Finite Element packages , 2012 .

[16]  Alain Combescure,et al.  Implicit/explicit multi-time step co-computations for predicting reinforced concrete structure response under earthquake loading , 2012 .

[17]  Youqing Wang,et al.  Equivalent dynamic infinite element for soil-structure interaction , 2013 .

[18]  Tsili Wang,et al.  Finite-difference modeling of elastic wave propagation: A nonsplitting perfectly matched layer approach , 2003 .

[19]  Nathan M. Newmark,et al.  A Method of Computation for Structural Dynamics , 1959 .

[20]  A. Chopra,et al.  Perfectly matched layers for transient elastodynamics of unbounded domains , 2004 .

[21]  A. Chopra,et al.  Erratum to ‘Perfectly matched layers for transient elastodynamics of unbounded domains’ (Int. J. Numer. Meth. Engng 2004; 59:1039–1074) , 2004 .

[22]  Loukas F. Kallivokas,et al.  Time‐domain hybrid formulations for wave simulations in three‐dimensional PML‐truncated heterogeneous media , 2015 .

[23]  Irini Djeran-Maigre,et al.  Multi-directional and multi-time step absorbing layer for unbounded domain , 2014 .

[24]  L. Lenti,et al.  A simple multi‐directional absorbing layer method to simulate elastic wave propagation in unbounded domains , 2010, 1009.1002.

[25]  Charbel Farhat,et al.  A transient FETI methodology for large‐scale parallel implicit computations in structural mechanics , 1994 .

[26]  C. Tsogka,et al.  Application of the perfectly matched absorbing layer model to the linear elastodynamic problem in anisotropic heterogeneous media , 2001 .

[27]  Loukas F. Kallivokas,et al.  Mixed perfectly-matched-layers for direct transient analysis in 2D elastic heterogeneous media , 2011 .

[28]  Horace Lamb,et al.  On the propagation of tremors over the surface of an elastic solid , 1904, Proceedings of the Royal Society of London.

[29]  René Matzen.  An efficient finite element time‐domain formulation for the elastic second‐order wave equation: A non‐split complex frequency shifted convolutional PML , 2011 .

[30]  John B. Schneider,et al.  Application of the perfectly matched layer (PML) absorbing boundary condition to elastic wave propagation , 1996 .

[31]  A. Chopra,et al.  Perfectly matched layers for time-harmonic elastodynamics of unbounded domains : Theory and finite-element implementation , 2003 .

[32]  Anthony Gravouil,et al.  Heterogeneous asynchronous time integrators for computational structural dynamics , 2015 .

[33]  Ushnish Basu,et al.  Explicit finite element perfectly matched layer for transient three‐dimensional elastic waves , 2009 .