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

JOURNAL OF COMPUTATIONAL PHYSICS

卷:195

Article

Lombard, B ; Piraux, J

关键词: acoustics and elastodynamics;    discontinuous coefficients;    jump conditions;    interface methods;    wave propagation algorithm;    singular value decomposition;    hyperbolic systems;   

DOI  :  10.1016/j.jcp.2003.09.024

来源: Elsevier

【 摘 要 】

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 ID [J. Comput. Phys. 168 (2001) 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, etc.) 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. (C) 2003 Elsevier Inc. All rights reserved.

【 授权许可】

Free