Numerical treatment of two-dimensional interfaces for acoustic and elastic waves

Abstract : We present a numerical method to take into account 2D arbitrary-shaped interfaces in classical finite-difference schemes, on a uniform Cartesian grid. This work extends the "Explicit Simplified Interface Method" (ESIM), previously proposed in 1D (2001, J. Comput. Phys., vol 168, pp. 227-248). The physical problem under study concerns the linear hyperbolic systems of acoustics and elastodynamics, with stationary interfaces. Our method maintains, near the interfaces, properties of the schemes in homogeneous medium, such as the order of accuracy and the stability limit. Moreover, it enforces the numerical solution to satisfy the exact interface conditions. Lastly, it provides subcell geometrical features of the interface inside the meshing. The ESIM can be coupled automatically with a wide class of numerical schemes (Lax-Wendroff, flux-limiter schemes,...) for a negligible additional computational cost. Throughout the paper, we focus on the challenging case of an interface between a fluid and an elastic solid. In numerical experiments, we provide comparisons between numerical solutions and original analytic solutions, showing the efficiency of the method.
Document type :
Journal articles
Complete list of metadatas

Cited literature [26 references]  Display  Hide  Download
Contributor : Bruno Lombard <>
Submitted on : Friday, April 29, 2005 - 2:24:07 PM
Last modification on : Monday, March 4, 2019 - 2:04:03 PM
Long-term archiving on : Thursday, April 1, 2010 - 9:23:50 PM



Bruno Lombard, Joël Piraux. Numerical treatment of two-dimensional interfaces for acoustic and elastic waves. Journal of Computational Physics, Elsevier, 2004, 195 (1), pp.90-116. ⟨10.1016/⟩. ⟨hal-00004813⟩



Record views


Files downloads