An official website of the United States government
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
A Fluid-Structure Coupled Computational Model for the Certification of Shock-Resistant Elastomer Coatings
Abstract Shock waves from underwater and air explosions are significant threats to surface and underwater vehicles and structures. Recent studies on the mechanical and thermal properties of various phase-separated elastomers indicate the possibility of applying these materials as a coating to mitigate shock-induced structural failures. To demonstrate this approach and investigate its efficacy, this paper presents a fluid-structure coupled computational model capable of predicting the dynamic response of air-backed bilayer (i.e. elastomer coating – metal substrate) structures submerged in water to hydrostatic and underwater explosion loads. The model couples a three-dimensional multiphase finite volume computational fluid dynamics model with a nonlinear finite element computational solid dynamics model using the FIVER (FInite Volume method with Exact multi-material Riemann solvers) method. The kinematic boundary condition at the fluid-structure interface is enforced using an embedded boundary method that is capable of handling large structural deformation and topological changes. The dynamic interface condition is enforced by formulating and solving local, one-dimensional fluid-solid Riemann problems, which is well-suited for transferring shock and impulsive loads. The capability of this computational model is demonstrated through a numerical investigation of hydrostatic and shock-induced collapse of aluminum tubes with polyurea coating on its inner surface. The thickness of the structure is resolved explicitly by the finite element mesh. The nonlinear material behavior of polyurea is accounted for using a hyper-viscoelastic constitutive model featuring a modified Mooney-Rivlin equation and a stress relaxation function in the form of prony series. Three numerical experiments are conducted to simulate and compare the collapse of the structure in different loading conditions, including a constant pressure, a fluid environment initially in hydrostatic equilibrium, and a two-phase fluid flow created by a near-field underwater explosion.
more »
« less
- Award ID(s):
- 1751487
- PAR ID:
- 10211313
- Date Published:
- Journal Name:
- ASME 2020 39th International Conference on Ocean, Offshore and Arctic Engineering
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
De Lorenzis, Laura; Papadrakakis, Manolis; Zohdi, Tarek I. (Ed.)This paper presents a graph-manifold iterative algorithm to predict the configurations of geometrically exact shells subjected to external loading. The finite element solutions are first stored in a weighted graph where each graph node stores the nodal displacement and nodal director. This collection of solutions is embedded onto a low-dimensional latent space through a graph isomorphism encoder. This graph embedding step reduces the dimensionality of the nonlinear data and makes it easier for the response surface to be constructed. The decoder, in return, converts an element in the latent space back to a weighted graph that represents a finite element solution. As such, the deformed configuration of the shell can be obtained by decoding the predictions in the latent space without running extra finite element simulations. For engineering applications where the shell is often subjected to concentrated loads or a local portion of the shell structure is of particular interest, we use the solutions stored in a graph to reconstruct a smooth manifold where the balance laws are enforced to control the curvature of the shell. The resultant computer algorithm enjoys both the speed of the nonlinear dimensional reduced solver and the fidelity of the solutions at locations where it matters.more » « less
-
Summary We present a spatially varying Robin interface condition for solving fluid‐structure interaction problems involving incompressible fluid flows and nonuniform flexible structures. Recent studies have shown that for uniform structures with constant material and geometric properties, a constant one‐parameter Robin interface condition can improve the stability and accuracy of partitioned numerical solution procedures. In this work, we generalize the parameter to a spatially varying function that depends on the structure's local material and geometric properties, without varying the exact solution of the coupled fluid‐structure system. We present an algorithm to implement the Robin interface condition in an embedded boundary method for coupling a projection‐based incompressible viscous flow solver with a nonlinear finite element structural solver. We demonstrate the numerical effects of the spatially varying Robin interface condition using two example problems: a simplified model problem featuring a nonuniform Euler‐Bernoulli beam interacting with an inviscid flow and a generalized Turek‐Hron problem featuring a nonuniform, highly flexible beam interacting with a viscous laminar flow. Both cases show that a spatially varying Robin interface condition can clearly improve numerical accuracy (by up to two orders of magnitude in one instance) for the same computational cost. Using the second example problem, we also demonstrate and compare two models for determining the local value of the combination function in the Robin interface condition.more » « less
-
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
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Conference Paper:
- https://doi.org/10.1115/OMAE2020-18501
