Coupling with the Embedded Boundary Method in a Runge-Kutta Discontinuous-Galerkin Direct Ghost-Fluid Method (RKDG-DGFM) Framework for Fluid-Structure Interaction Simulations of Underwater Explosions

Solution of near-field underwater explosion (UNDEX) problems frequently require the modeling of two-way coupled fluid-structure interaction (FSI). This paper describes the addition of an embedded boundary method to an UNDEX modeling framework for multiphase, compressible and inviscid fluid using the combined algorithms of Runge-Kutta, discontinuous-Galerkin, level-set and direct ghost-fluid methods. A computational fluid dynamics (CFD) solver based on these algorithms has been developed as described in previous work. A fluid-structure coupling approach was required to perform FSI simulation interfacing with an external structural mechanics solver. Large structural deformation and possible rupture and cracking characterize the FSI phenomenon in an UNDEX, so the embedded boundary method (EBM) is more appealing for this application in comparison to dynamic mesh methods such as the arbitrary Lagrangian-Eulerian (ALE) method to enable the fluid-structure coupling algorithm in the fluid. Its limitation requiring a closed interface that is fully submerged in the fluid domain is relaxed by an adjustment described in this paper so that its applicability is extended. Two methods of implementing the fluid-structure wall boundary condition are also compared. The first solves a local 1D fluid-structure Riemann problem at each intersecting point between the wetted elements and fluid mesh. In this method, iterations are required when the Tait equation of state is utilized. A second method that does not require the Riemann solution and iterations is also implemented and the results are compared.

Damped wave-equation-based first-arrival traveltime tomography using the embedded boundary method

First-arrival traveltime tomography (FATT) is used to delineate shallow velocity structures to identify static effects in oil exploration as well as to characterize the near surface for geotechnical purposes. Because FATT is generally used for land seismic data processing, it becomes necessary to consider irregular topography especially when performing wave-based tomography. However, the standard Cartesian finite-difference method cannot properly handle irregular topography. Hence, the embedded boundary method (EBM) is incorporated into the frequency-domain damped-wave equation in order to correctly describe irregular topography. The developed modeling algorithm is used to calculate first-arrival traveltimes and to perform FATT. The EBM-based modeling algorithm accurately describes the irregular surfaces of numerical velocity models on a regular mesh by exploiting the mirror image principle. The accuracy of the EBM-based traveltime calculation is validated by using two homogeneous velocity models with dipping and complex surfaces. The validation results demonstrate that the proposed algorithm is unaffected by the staircase approximation. The FATT is then applied to synthetic and real data to demonstrate the applicability of the developed algorithm to velocity models with complex topography. For the real data example, the inverted velocity model is used to apply static corrections. The processing results demonstrate an improvement in the continuity of seismic events, thus confirming the accuracy of the developed FATT method.