skip to main content

Title: A Fluid-Structure Coupled Computational Model for the Certification of Shock-Resistant Elastomer Coatings

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

« less
; ;
Award ID(s):
Publication Date:
Journal Name:
ASME 2020 39th International Conference on Ocean, Offshore and Arctic Engineering
Sponsoring Org:
National Science Foundation
More Like this
  1. 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 muchmore »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.« less
  2. As a step towards addressing a scarcity of references on this topic, we compared the Eulerian and Lagrangian Computational Fluid Dynamics (CFD) approaches for the solution of free-surface and Fluid–Solid Interaction (FSI) problems. The Eulerian approach uses the Finite Element Method (FEM) to spatially discretize the Navier–Stokes equations. The free surface is handled via the volume-of-fluid (VOF) and the level-set (LS) equations; an Immersed Boundary Method (IBM) in conjunction with the Nitsche’s technique were applied to resolve the fluid–solid coupling. For the Lagrangian approach, the smoothed particle hydrodynamics (SPH) method is the meshless discretization technique of choice; no additional equations are needed to handle free-surface or FSI coupling. We compared the two approaches for a flow around cylinder. The dam break test was used to gauge the performance for free-surface flows. Lastly, the two approaches were compared on two FSI problems—one with a floating rigid body dropped into the fluid and one with an elastic gate interacting with the flow. We conclude with a discussion of the robustness, ease of model setup, and versatility of the two approaches. The Eulerian and Lagrangian solvers used in this study are open-source and available in the public domain.
  3. Metal-ion batteries (e.g., lithium and sodium ion batteries) are the promising power sources for portable electronics, electric vehicles, and smart grids. Recent metal-ion batteries with organic liquid electrolytes still suffer from safety issues regarding inflammability and insufficient lifetime.1 As the next generation energy storage devices, all-solid-state batteries (ASSBs) have promising potentials for the improved safety, higher energy density, and longer cycle life than conventional Li-ion batteries.2 The nonflammable solid electrolytes (SEs), where only Li ions are mobile, could prevent battery combustion and explosion since the side reactions that cause safety issues as well as degradation of the battery performance are largely suppressed. However, their practical application is hampered by the high resistance arising at the solid–solid electrode–electrolyte interface (including cathode-electrolyte interface and anode-electrolyte interface).3 Several methods have been introduced to optimize the contact capability as well as the electrochemical/chemical stability between the metal anodes (i.e.: Li and Na) and the SEs, which exhibited decent results in decreasing the charge transfer resistance and broadening the range of the stable energy window (i.e., lowing the chemical potential of metal anode below the highest occupied molecular orbital of the SEs).4 Nevertheless, mitigation for the cathode in ASSB is tardily developed because: (1) themore »porous structure of the cathode is hard to be infiltrated by SEs;5 (2) SEs would be oxidized and decomposed by the high valence state elements at the surface of the cathode at high state of charge.5 Herein, we demonstrate a universal cathode design strategy to achieve superior contact capability and high electrochemical/chemical stability with SEs. Stereolithography is adopted as a manufacturing technique to realize a hierarchical three-dimensional (HTD) electrode architecture with micro-size channels, which is expected to provide larger contact areas with SEs. Then, the manufactured cathode is sintered at 700 °C in a reducing atmosphere (e.g.: H2) to accomplish the carbonization of the resin, delivering sufficiently high electronic conductivity for the cathode. To avoid the direct exposure of the cathode active materials to the SEs, oxidative chemical vapor deposition technique (oCVD) is leveraged to build conformal and highly conducting poly(3,4-ethylenedioxythiophene) (PEDOT) on the surface of the HTD cathode.6 To demonstrate our design strategy, both NCM811 and Na3V2(PO4)3 is selected as active materials in the HTD cathode, then each cathode is paired with organic (polyacrylonitrile-based) and inorganic (sulfur-based) SEs assembled into two batteries (total four batteries). SEM and TEM reveal the micro-size HTD structure with built-in channels. Featured by the HTD architecture, the intrinsic kinetic and thermodynamic conditions will be enhanced by larger surface contact areas, more active sites, improved infusion and electrolyte ion accessibility, and larger volume expansion capability. Disclosed by X-ray computed tomography, the interface between cathode and SEs in the four modified samples demonstrates higher homogeneity at the interface between the cathode and SEs than that of all other pristine samples. Atomic force microscopy is employed to measure the potential image of the cross-sectional interface by the peak force tapping mode. The average potential of modified samples is lower than that of pristine samples, which confirms a weakened space charge layer by the enhanced contact capability. In addition, through Electron Energy Loss Spectroscopy coupled with Scanning Transmission Electron Microscopy, the preserved interface between HTD cathode and SE is identified; however, the decomposing of the pristine cathode is clearly observed. In addition, Finite element method simulations validate that the diffusion dynamics of lithium ions is favored by HTD structure. Such a demonstrated universal strategy provides a new guideline to engineer cathode electrolyte interface by reconstructing electrode structures that can be applicable to all solid-state batteries in a wide range of chemical conditions.« less
  4. All-solid-state batteries (ASSBs) have garnered increasing attention due to the enhanced safety, featuring nonflammable solid electrolytes as well as the potential to achieve high energy density. 1 The advancement of the ASSBs is expected to provide, arguably, the most straightforward path towards practical, high-energy, and rechargeable batteries based on metallic anodes. 1 However, the sluggish ion transmission at the cathode-electrolyte (solid/solid) interface would result in the high resistant at the contact and limit the practical implementation of these all solid-state materials in real world batteries. 2 Several methods were suggested to enhance the kinetic condition of the ion migration between the cathode and the solid electrolyte (SE). 3 A composite strategy that mixes active materials and SEs for the cathode is a general way to decrease the ion transmission barrier at the cathode-electrolyte interface. 3 The active material concentration in the cathode is reduced as much as the SE portion increases by which the energy density of the ASSB is restricted. In addition, the mixing approach generally accompanies lattice mismatches between the cathode active materials and the SE, thus providing only limited improvements, which is imputed by random contacts between the cathode active materials and the SE during the mixingmore »process. Implementing high-pressure for the electrode and electrolyte of ASSB in the assembling process has been verified is a but effective way to boost the ion transmission ability between the cathode active materials and the SE by decreasing the grain boundary impedance. Whereas the short-circuit of the battery would be induced by the mechanical deformation of the electrolyte under high pressure. 4 Herein, we demonstrate a novel way to address the ion transmission problem at the cathode-electrolyte interface in ASSBs. Starting from the cathode configuration, the finite element method (FEM) was employed to evaluate the current concentration and the distribution of the space charge layer at the cathode-electrolyte interface. Hierarchical three-dimensional (HTD) structures are found to have a higher Li + transfer number (t Li+ ), fewer free anions, and the weaker space-charge layer at the cathode-electrolyte interface in the resulting FEM simulation. To take advantage of the HTD structure, stereolithography is adopted as a manufacturing technique and single-crystalline Ni-rich (SCN) materials are selected as the active materials. Next, the manufactured HTD cathode is sintered at 600 °C in an N 2 atmosphere for the carbonization of the resin, which induces sufficient electronic conductivity for the cathode. Then, the gel-like Li 1.4 Al 0.4 Ti 1.6 (PO 4 ) 3 (LATP) precursor is synthesized and filled into the voids of the HTD structure cathode sufficiently. And the filled HTD structure cathodes are sintered at 900 °C to achieve the crystallization of the LATP gel. Scanning transmission electron microscopy (STEM) is used to unveil the morphology of the cathode-electrolyte interface between the sintered HTD cathode and the in-situ generated electrolyte (LATP). A transient phase has been found generated at the interface and matched with both lattices of the SCN and the SE, accelerating the transmission of the Li-ions, which is further verified by density functional theory calculations. In addition, Electron Energy Loss Spectroscopy demonstrates the preserved interface between HTD cathode and SEs. Atomic force microscopy is employed to measure the potential image of the cross-sectional interface by the peak force tapping mode. The average potential of modified samples is lower than the sample that mix SCN and SEs simply in the 2D planar structure, which confirms a weakened space charge layer by the enhanced contact capability as well as the ion transmission ability. To see if the demonstrated method is universally applicable, LiNi 0.8 Co 0.1 Mn 0.1 O 2 (NCM811) is selected as the cathode active material and manufactured in the same way as the SCN. The HTD cathode based on NCM811 exhibits higher electrochemical performance compared with the reference sample based on the 2D planar mixing-type cathode. We believe such a demonstrated universal strategy provides a new guideline to engineer the cathode/electrolyte interface by revolutionizing electrode structures that can be applicable to all-solid-state batteries. Figure 1. Schematic of comparing of traditional 2D planar cathode and HTD cathode in ASSB Tikekar, M. D. , et al. , Nature Energy (2016) 1 (9), 16114 Banerjee, A. , et al. , Chem Rev (2020) 120 (14), 6878 Chen, R. , et al. , Chem Rev (2020) 120 (14), 6820 Cheng, X. , et al. , Advanced Energy Materials (2018) 8 (7) Figure 1« less
  5. Numerical modelling of deformation in hydromechanical systems can be time-consuming using fully coupled classical numerical methods for large representative porous media samples. In this paper, we present a new two-way coupled partitioned fluid–solid system. The coupled system is applied for simulating geomechanical processes at the pore-scale. We track the deformation of the solid resulting from the drainage of resident fluids in the pores, as well as the evolution of fluid properties from dynamic loading. The finite element method is responsible for capturing the structural deformation in the coupled system while the dynamic pore network is used for modelling multiphase flow in the fluid subsystem. A fictitious fluid–solid interface is created at each pore network-finite element node junction via convex hulling, followed by data exchange using linear interpolation. The results show good agreement with a pre-existing coupled finite volume model and the computations are completed in much less time.