The modified ghost fluid method (MGFM) provides a robust and efficient interface treatment for various multi-medium flow simulations and some particular fluid-structure interaction (FSI) simulations. However, this methodology for one specific class of FSI problems, where the structure is plate, remains to be developed. This work is devoted to extending the MGFM to treat compressible fluid coupled with a thin elastic plate. In order to take into account the influence of simultaneous interaction at the interface, a fluid-plate coupling system is constructed at each time step and solved approximately to predict the interfacial states. Then, ghost fluid states and plate load can be defined by utilizing the obtained interfacial states. A type of acceleration strategy in the coupling process is presented to pursue higher efficiency. Several one-dimensional examples are used to highlight the utility of thismethod over looselycoupled method and validate the acceleration techniques. Especially, this method is applied to compute the underwater explosions (UNDEX) near thin elastic plates. Evolution of strong shock impacting on the thin elastic plate and dynamic response of the plate are investigated. Numerical results disclose that this methodology for treatment of the fluid-plate coupling indeed works conveniently and accurately for different structural flexibilities and is capable of efficiently simulating the processes of UNDEX with the employment of the acceleration strategy.

Numerical methods are developed to simulate the wave propagation in heterogeneous 2D fluid/poroelastic media. Wave propagation is described by the usual acoustics equations (in the fluid medium) and by the low-frequency Biot’s equations (in the porous medium). Interface conditions are introduced to model various hydraulic contacts between the two media: open pores, sealed pores, and imperfect pores. Well-posedness of the initial-boundary value problem is proven. Cartesian grid numerical methods previously developed in porous heterogeneous media are adapted to the present context: a fourth-order ADER scheme with Strang splitting for time- marching; a space-time mesh-refinement to capture the slow compressional wave predicted by Biot’s theory; and an immersed interface method to discretize the interface conditions and to introduce a subcell resolution. Numerical experiments and comparisons with exact solutions are proposed for the three types of interface conditions, demonstrating the accuracy of the approach.

In this work, the modified ghost fluid method is developed to deal with 2D compressible fluid interacting with elastic solid in an Euler-Lagrange coupled system. In applying the modified Ghost Fluid Method to treat the fluid-elastic solid coupling, the Navier equations for elastic solid are cast into a system similar to the Euler equations but in Lagrangian coordinates. Furthermore, to take into account the influence of material deformation and nonlinear wave interaction at the interface, an Euler-Lagrange Riemann problem is constructed and solved approximately along the normal direction of the interface to predict the interfacial status and then define the ghost fluid and ghost solid states. Numerical tests are presented to verify the resultant method.

In this article, we discuss a least-squares/fictitious domain method for the solution of linear elliptic boundary value problems with Robin boundary conditions. Let Ω and ω be two bounded domains of Rd such that ω̅⊂Ω. For a linear elliptic problem in Ω\ω̅ with Robin boundary condition on the boundary ϒ of ω, our goal here is to develop a fictitious domain method where one solves a variant of the original problem on the full Ω, followed by a well-chosen correction over ω. This method is of the virtual control type and relies on a least-squares formulation making the problem solvable by a conjugate gradient algorithm operating in a well chosen control space. Numerical results obtained when applying our method to the solution of two-dimensional elliptic and parabolic problems are given; they suggest optimal order of convergence.