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. 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.

A hybrid FDTD-Rayleigh integral computational method for the simulation of the ultrasound measurement of proximal femur

“…The grid step ∆x should be sufficiently small to model the geometry accurately. Note that alternative FDTD schemes which have been proposed to model explicitly curved boundaries [29,30] could be used advantageously to model QUS experiments. Finally, the method of coupling between the Rayleigh integral formulation and the FDTD assumes that the diffracting object inside the FDTD box is in the far field.…”

Elastic surface waves in crystals – Part 2: Cross-check of two full-wave numerical modeling methods

“…A less straightforward issue using pseudospectral differential operators is to model the freesurface boundary condition. While in finite-element methods the implementation of traction-free boundary conditions is natural -simply do not impose any constraint at the surface nodes -finitedifference and pseudospectral methods require a particular boundary treatment [14,23,25,26].…”

High-order finite-volume methods for hyperbolic conservation laws on mapped multiblock grids

“…We therefore use least-squares interpolation to transfer cell-average information between the overlapping, curvilinear grids. Least-squares interpolation has also been used at block boundaries for embedded-boundary methods in [20,21] and for high-order coarse-fine mesh interpolation in AMR [25]. The methods of the present paper are applied in [26] to the solution of the shallow-water equations on the surface of a sphere, using AMR.…”