Hybrid asynchronous isogeometric Perfectly Matched Layer for transient elastodynamics. (June 2023)

• Record Type: Journal Article
• Authors: Chau, Anh-Khoa; Brun, Michael; Elguedj, Thomas
• Is Part Of: Computers and geotechnics. Volume 158(2023)
• Publication Date: 2023-06
• Subjects: Isogeometric analysis -- Perfectly Matched Layer -- Mortar method -- Heterogenous Asynchronous Time Integrator -- Time-domain elastic wave propagation; Engineering geology -- Data processing -- Periodicals; Soil mechanics -- Data processing -- Periodicals; Rock mechanics -- Data processing -- Periodicals; 624.1510285
• DOI: 10.1016/j.compgeo.2023.105387
• Languages: English
• ISSNs: 0266-352X

Abstract

Abstract: This paper presents an isogeometric formulation (IGA) of the Perfectly Matched Layer (PML) for two-dimensional time-domain elastodynamics in unbounded media. The displacement-based unsplit-field PML formulation is first discretized in space using B-spline basis functions of increasing polynomial order p, from 2 to 10, able to deal with the strong solution gradients in the thickness of the PML thanks to the high order continuity provided by the B-spline basis functions. In order to keep the high accuracy of IGA, implicit time integration is chosen in the PML subdomains, whereas the physical domain is integrated with an explicit time integration scheme. The domain decomposition is based on a dual formulation and Heterogeneous Asynchronous Time Integrator (HATI) approach. The coupling between non-conforming IGA-patches for the PML subdomains and the physical domain is carried out using a standard mortar method. It enables us to select for the IGA-PML subdomains both the appropriate space discretization and the appropriate time discretization with an implicit time integrator and a larger time step size. The hybrid multi-time step IGA-PML performance in terms of computation time, is assessed for two-dimensional elastodynamic problems in single and multi-layer media depending on the chosen discretization parameters in space and time. Highlights: Perfectly Matched Layers (PML) with Isogeometric Analysis (IGA) for 2D elastodynamics High order continuity of IGA to optimize the PML performance Explicit/Implicit multi time step coupling between FE subdomain and IGA-PML Mortar coupling to deal with non-conforming IGA-PML Computational gain in comparison to classical FEM-PML at about 10