Unfitted hybrid high-order methods for the wave equation

Erik Burman (1) , Omar Duran (2, 3, 4) , Alexandre Ern (2, 3)

1 UCL - University College of London [London]
2 CERMICS - Centre d'Enseignement et de Recherche en Mathématiques et Calcul Scientifique
3 SERENA - Simulation for the Environment: Reliable and Efficient Numerical Algorithms
4 UiB - University of Bergen

We design an unfitted hybrid high-order (HHO) method for the wave equation. The wave propagates in a domain where a curved interface separates subdomains with different material properties. The key feature of the discretization method is that the interface can cut more or less arbitrarily through the mesh cells. We address both the second-order formulation in time of the wave equation and its reformulation as a first-order system. We prove H 1-error estimates for a space semi-discretization in space of the second-order formulation, leading to optimal convergence rates for smooth solutions. Numerical experiments illustrate the theoretical findings and show that the proposed numerical schemes can be used to simulate accurately the propagation of acoustic waves in heterogeneous media, with meshes that are not fitted to the geometry. For the second-order formulation, the implicit, second-order accurate Newmark scheme is used for the time discretization, whereas (diagonally-implicit or explicit) Runge-Kutta schemes up to fourth-order accuracy are used for the first-order formulation. In the explicit case, we study the CFL condition on the time step and observe that the unfitted approach combined with local cell agglomeration leads to a comparable condition as when using fitted meshes.

Mots clés: Heterogeneous media; Wave equation; Error analysis; Unfitted mesh; Hybrid high-order methods

Domaines: Analyse numérique [math.NA]

Erik Burman, Omar Duran, Alexandre Ern. Unfitted hybrid high-order methods for the wave equation. Computer Methods in Applied Mechanics and Engineering, 2022, ⟨10.1016/j.cma.2021.114366⟩. ⟨hal-03086432v3⟩