More Like this

1. A Computational Analysis of Bubble-Structure Interaction in Near-Field Underwater Explosion

   https://doi.org/10.1115/IMECE2021-72854  

   Ma, Wentao; Ozenkoski, Timothy; Wang, Kevin (November 2021, Proceedings of ASME 2021 International Mechanical Engineering Congress and Exposition)

   Underwater explosion poses a significant threat to the structural integrity of ocean vehicles and platforms. Accurate prediction of the dynamic loads from an explosion and the resulting structural response is crucial to ensuring safety without overconservative design. When the distance between the explosive charge and the structure is relatively small (i.e., near-field explosion), the dynamics of the gaseous explosion product, i.e., the “bubble”, comes into play, rendering a multiphysics problem that features the interaction of the bubble, the surrounding liquid water, and the solid structure. The problem is highly nonlinear, as it involves shock waves, large deformation, yielding, contact, and possibly fracture. This paper investigates the two-way interaction between the cyclic expansion and collapse of an explosion bubble and the deformation of a thin-walled elastoplastic cylindrical shell in its vicinity. Intuitively, when a shock wave impinges on a thin cylindrical shell, the shell would collapse in the direction of shock propagation. However, some recent laboratory experiments have shown that under certain conditions the shell collapsed in a counter-intuitive mode in which the direction of collapse is perpendicular to that of shock propagation. In other words, the nearest point on the structural surface moved towards the explosion charge, despite being impacted by a compressive shock. This paper focuses on replicating this phenomenon through numerical simulation and elucidating the underlying mechanisms. A recently developed computational framework (“FIVER”) coupling a nonlinear finite element structural dynamics solver and a finite volume compressible fluid dynamics solver is used to complete this study. The solver utilizes an embedded boundary method to track the wetted surface of the structure (i.e. the fluid-structure interface), which is capable of handling large structural deformation and topological changes (e.g., fracture). The solver also adopts the level set method for tracking the bubble surface (i.e. the liquid-gas interface). The fluid-structure and liquid-gas interface conditions are enforced by constructing and solving one-dimensional multi-material Riemann problems, which naturally accommodates the propagation of shock waves across the interfaces. In this paper, mesh refinement study is made to examine the sensitivity of the results to various meshing parameters. The results show that the intermediate level of refinement is appropriate in terms of both the accuracy and the computation costs. Next, the deformation history of both the bubble and the structure are presented and analyzed to provide a detailed view of the counter-intuitive collapse mode mentioned above. We show that timewise, the structural collapse spans multiple cycles of bubble oscillation. Additional details about the time-histories of fluid pressure, structure displacement, and bubble size are presented to elucidate this dynamic bubble-structure interaction and the resulting structural failure. 

   more » « less
2. Shock-induced bubble collapse near solid materials: effect of acoustic impedance

   https://doi.org/10.1017/jfm.2020.810  

   Cao, S.; Wang, G.; Coutier-Delgosha, O.; Wang, K. (January 2021, Journal of Fluid Mechanics)

   null (Ed.)

   The fluid dynamics of a bubble collapsing near an elastic or viscoelastic material is coupled with the mechanical response of the material. We apply a multiphase fluid–solid coupled computational model to simulate the collapse of an air bubble in water induced by an ultrasound shock wave, near different types of materials including metals (e.g. aluminium), polymers (e.g. polyurea), minerals (e.g. gypsum), glass and foams. We characterize the two-way fluid–material interaction by examining the fluid pressure and velocity fields, the time history of bubble shape and volume and the maximum tensile and shear stresses produced in the material. We show that the ratio of the longitudinal acoustic impedance of the material compared to that of the ambient fluid, $$Z/Z_0$$ , plays a significant role. When $$Z/Z_0<1$$ , the material reflects the compressive front of the incident shock into a tensile wave. The reflected tensile wave impinges on the bubble and decelerates its collapse. As a result, the collapse produces a liquid jet, but not necessarily a shock wave. When $$Z/Z_0>1$$ , the reflected wave is compressive and accelerates the bubble's collapse, leading to the emission of a shock wave whose amplitude increases linearly with $$\log (Z/Z_0)$$ , and can be much higher than the amplitude of the incident shock. The reflection of this emitted shock wave impinges on the bubble during its rebound. It reduces the speed of the bubble's rebound and the velocity of the liquid jet. Furthermore, we show that, for a set of materials with $$Z/Z_0\in [0.04, 10.8]$$ , the effect of acoustic impedance on the bubble's collapse time and minimum volume can be captured using phenomenological models constructed based on the solution of Rayleigh–Plesset equation. 

   more » « less
3. Introduction to the Distributional Finite Difference Method for 3D Seismic Wave Propagation and Comparison With the Spectral Element Method

   https://doi.org/10.1029/2023JB027576  

   Lyu, Chao; Masson, Yder; Romanowicz, Barbara; Zhao, Liang (April 2024, Journal of Geophysical Research: Solid Earth)

   Abstract We have extended the distributional finite difference method (DFDM) to simulate the seismic‐wave propagation in 3D regional earth models. DFDM shares similarities to the discontinuous finite element method on a global scale and to the finite difference method locally. Instead of using linear staggered finite‐difference operators, we employ DFDM operators based on B‐splines and a definition of derivatives in the sense of distributions, to obtain accurate spatial derivatives. The weighted residuals method used in DFDM's locally weak formalism of spatial derivatives accurately and naturally accounts for the free surface, curvilinear meshing, and solid‐fluid coupling, for which it only requires setting the shear modulus and the continuity condition to zero. The computational efficiency of DFDM is comparable to the spectral element method (SEM) due to the more accurate mass matrix and a new band‐diagonal mass factorization. Numerical examples show that the superior accuracy of the band‐diagonal mass and stiffness matrices in DFDM enables fewer points per wavelength, approaching the spectral limit, and compensates for the increased computational burden due to four Lebedev staggered grids. Specifically, DFDM needs 2.5 points per wavelength, compared to the five points per wavelength required in SEM for 0.5% waveform error in a homogeneous model. Notably, while maintaining the same accuracy in the solid domain, DFDM demonstrates superior accuracy in the fluid domain compared to SEM. To validate its accuracy and flexibility, we present various 3D benchmarks involving homogeneous and heterogeneous regional elastic models and solid‐fluid coupling in both Cartesian and spherical settings. 

   more » « less
4. Applicability of 3D Spectral Element Method for Computing Close-Range Underwater Piling Noises

   https://doi.org/10.1142/S2591728519500129  

   Jeong, C.; Manalaysay, A.; Gharti, H. N.; Guan, S.; Vignola, J. (December 2019, Journal of Theoretical and Computational Acoustics)

   null (Ed.)

   Pile driving is used for constructing foundation supports for offshore structures. Underwater noise, induced by in-water pile driving, could adversely impact marine life near the piling location. Many studies have computed this noise in close ranges by using semi-analytical models and Finite Element Method (FEM) models. This work presents a Spectral Element Method (SEM) wave simulator as an alternative simulation tool to obtain close-range underwater piling noise in complex, fully three-dimensional, axially-asymmetric settings in the time domain for impacting force signals with high-frequency contents (e.g., frequencies greater than 1000[Formula: see text]Hz). The presented numerical results show that the flexibility of SEM can accommodate the axially-asymmetric geometry of a model, its heterogeneity, and fluid-solid coupling. We showed that there are multiple Mach Cones of different angles in fluid and sediment caused by the difference in wave speeds in fluid, a pile, and sediment. The angles of Mach Cones in our numerical results match those that are theoretically evaluated. A previous work 18 had shown that Mach Cone waves lead to intense amplitudes of underwater piling noise via a FEM simulation in an axis-symmetric setting. Since it modeled sediment as fluid with a larger wave speed than that of water, we examined if our SEM simulation, using solid sediment–fluid coupling, leads to additional Mach Cones. Because this work computes the shear wave in sediment and the downward-propagating shear wave in a pile, we present six Mach Cones in fluid and sediment induced by downward-propagating P- and S-waves in a pile in lieu of two previously-reported Mach Cones in fluid and sediment (modeled as fluid) induced by a downward-propagating P-wave in a pile. We also showed that the amplitudes of the close-range underwater noise are dependent on the cross-sectional geometry of a pile. In addition, when a pile is surrounded by a solid of an axially-asymmetric geometry, waves are reflected from the surface of the surrounding solid back to the fluid so that constructive and destructive interferences of waves take place in the fluid and affect the amplitude of the underwater piling noise. 

   more » « less
5. A Spectral/hp-Based Stabilized Solver with Emphasis on the Euler Equations

   https://doi.org/10.3390/fluids9010018  

   Ranjan, Rakesh; Catabriga, Lucia; Araya, Guillermo (January 2024, Fluids)

   The solution of compressible flow equations is of interest with many aerospace engineering applications. Past literature has focused primarily on the solution of Computational Fluid Dynamics (CFD) problems with low-order finite element and finite volume methods. High-order methods are more the norm nowadays, in both a finite element and a finite volume setting. In this paper, inviscid compressible flow of an ideal gas is solved with high-order spectral/hp stabilized formulations using uniform high-order spectral element methods. The Euler equations are solved with high-order spectral element methods. Traditional definitions of stabilization parameters used in conjunction with traditional low-order bilinear Lagrange-based polynomials provide diffused results when applied to the high-order context. Thus, a revision of the definitions of the stabilization parameters was needed in a high-order spectral/hp framework. We introduce revised stabilization parameters, τsupg, with low-order finite element solutions. We also reexamine two standard definitions of the shock-capturing parameter, δ: the first is described with entropy variables, and the other is the YZβ parameter. We focus on applications with the above introduced stabilization parameters and analyze an array of problems in the high-speed flow regime. We demonstrate spectral convergence for the Kovasznay flow problem in both L1 and L2 norms. We numerically validate the revised definitions of the stabilization parameter with Sod’s shock and the oblique shock problems and compare the solutions with the exact solutions available in the literature. The high-order formulation is further extended to solve shock reflection and two-dimensional explosion problems. Following, we solve flow past a two-dimensional step at a Mach number of 3.0 and numerically validate the shock standoff distance with results obtained from NASA Overflow 2.2 code. Compressible flow computations with high-order spectral methods are found to perform satisfactorily for this supersonic inflow problem configuration. We extend the formulation to solve the implosion problem. Furthermore, we test the stabilization parameters on a complex flow configuration of AS-202 capsule analyzing the flight envelope. The proposed stabilization parameters have shown robustness, providing excellent results for both simple and complex geometries. 

   more » « less