Prerpints.org logo
Preprint
Article

Application of CWENO Finite-Volume Interface Capturing Schemes for Modeling Cavitation Induced by an Underwater Explosion

Altmetrics

Downloads

126

Views

44

Comments

0

A peer-reviewed article of this preprint also exists.

Submitted:

20 December 2023

Posted:

21 December 2023

You are already at the latest version

Alerts
Abstract
Cavitation resulting from underwater explosions in compressible multiphase or multicomp-onent flows presents significant challenges due to the dynamic nature of shock-cavitation-structure interactions, as well as the complex and discontinuous nature of the involved interfaces. Achieving accurate resolution of interfaces between different phases or components, in the presence of shocks, cavitating regions, and structural interactions, is crucial for modeling such problems. Furthermore, pressure convergence in simulations involving shock-cavitation-structure interactions requires accurate algorithms. In this research paper, we employ the diffuse interface method, also known as the interface capturing scheme, to investigate cavitation in various underwater explosion test cases near different surfaces; free surface and a rigid surface. The simulations are conducted using the Unstructured Compressible Navier-Stokes (UCNS3D) finite volume framework employing central-weighted essentially non-oscillatory (CWENO) reconstruction schemes, utilising the five-equation diffuse interface family of methods. Quantitative comparisons are made between the performance of both models. Additionally, we examine the effects of cavitation as a secondary loading source on structures, and evaluate the ability of the CWENO schemes to accurately capture and resolve material interfaces between fluids with minimal numerical dissipation or smearing. The results are compared with existing high-order methods and experimental data, where possible, to demonstrate the robustness of the CWENO schemes in simulating cavitation bubble dynamics, as well as their limitations within the current implementation of interface capturing.
Keywords: 
Subject: Engineering  -   Marine Engineering

1. Introduction

Cavitation, which refers to the formation of vapor bubbles within a liquid due to high velocities and strong acceleration at low pressure, has diverse applications in various scientific and engineering fields. It can be beneficial in marine design engineering for reducing skin friction drag and achieving high speed for submerged projectiles, as well as in medicine for treating kidney stones and cancerous bodies. However, cavitation can also have adverse effects such as material erosion, noise, and vibrations in marine propellers and turbomachinery, that can cause significant damage. Therefore, understanding the dynamics of cavitation is crucial, but it poses challenges in numerical simulations of compressible multiphase problems due to the complex physics involved in phase change, high gradients of flow variables, numerical discretisation choices, and validation with experimental results. Numerical studies often focus on attached cavitations, as they are the primary subject of investigation in most research [1,2,3]. However, hull/bulk cavitation, which occurs during underwater explosions, is different as it evolves rapidly and unsteadily, changing dimensions before collapsing violently due to the surge in pressure, resulting in erratic and rapidly changing cavitation patterns. This makes numerical simulations of hull/bulk cavitation more challenging.
In hull cavitation, vapour bubbles are created on the surface of a submerged rigid surface, such as a projectile, ship hull, or submarine, due to rarefaction caused by high-pressure shock waves from an explosion. On the other hand, bulk cavitation occurs near a free surface between two interacting fluids due to the explosive shock waves, leading to widespread cavitation throughout the liquid. Hull cavitation is often more challenging to study using numerical methods due to its rapid and unsteady nature, while bulk cavitation is characterised by the formation of large-scale vapour bubbles throughout the liquid. Both forms of cavitation can have detrimental effects on marine structures and machinery.
In the context of hull or bulk cavitation, where regions experience extreme variations in pressure, negative pressure values are observed, which are physically unrealistic and can cause numerical codes to fail. To address this issue, two common approaches are used in diffuse interface models for simulating UNDEX (Underwater Explosion) phenomena. The first approach is to incorporate a single-fluid model, which has been used by several researchers [4,5,6,7,8,9] in similar UNDEX simulations. Various methods in this category, such as the cut-off model, vacuum model, isentropic model, Schmidt model, and modified Schmidt model, have been used to avoid negative pressure values during cavitation. The cut-off model assumes that cavitation bubbles do not form until the pressure drops below a critical value, known as the cut-off pressure, and is generally accurate in predicting peak pressure at the cavitation boundary, but may not accurately predict cavitation collapse periods. The Schmidt model is accurate when the void fraction (ratio of vapour bubble density to total fluid density) is less than 10 -5 [4], but has limitations in calculating the speed of sound, making it suitable for high-speed cavitating flows. The modified Schmidt model is an extension of the original Schmidt model and addresses its limitations, thereby expanding its range of applicability to larger vapor-to-liquid density ratios. The vacuum cavitation model assumes zero mass inside the cavitation bubble but is limited in capturing the complex shapes and dynamics of cavitation bubbles in multi-dimensional flows. The isentropic model, which has been used in the five-equation diffuse interface models [5,10], is robust as it considers the mixture as consisting of isentropic vapour and liquid phases, but it requires more computational resources and the predetermination of a model parameter k [4]. The cut-off model is preferred in this case as it is relatively simple to implement and does not alter peak pressure accuracy at the cavitation boundary, which is important for the UNDEX test cases considered in the studies. The one-fluid cavitation model aligns with the no-phase transition assumption of the five-equation interface capturing models used in this study.
The accuracy of the one-fluid model in predicting the details of the cavitation process, including thermodynamic aspects, in compressible two-phase flows, is limited due to its inability to account for mass and heat transfer between the liquid and vapour phases. To overcome this limitation, researchers have extended the simplified diffuse-interface models such as the four-, five-, or six-equation models by incorporating source terms in the energy and mass conservation equations to represent mass and heat transfer. This approach, known as the "two-fluid cavitation models," aims to prevent non-physical pressure occurrences in the cavitation region. The source terms are implicitly expressed using an algebraic relaxation solver for phase transition modelling, as proposed by Pelanti and Shyue [11,12]. However, this method is more complex and computationally expensive, especially for the seven-equation diffuse interface models such as the full non-equilibrium models of Baer-Nunziato [13] or Saurel-Abgrall [14]. A computational efficient and simple strategy in UNDEX is to employ a reduced diffuse interface models with phase transition, such as the four- or five-equation models, which enforce mechanical, thermal, or both equilibrium. Jun et al [15] extended the phase transition model used by Chiapolino et al. [16] for the four-equation model to study cavitation in underwater explosion test cases. Other existing diffuse interface models with phase transitions include those proposed by Pelanti and Shyue [9,11], Martelot et al. [17], Zein et al.[18], and others [14]. Ma et al. [19] utilised the Kapila [20] variant of the five-equation two-phase flow model with thermo-chemical relaxation terms to investigate liquid-vapour phase transition in cavitating flows for UNDEX test cases.
The present study focuses on hull and bulk cavitation in compressible multiphase or multi-component flows, specifically in the context of underwater explosions. Two interface capturing models, namely the ones proposed by Allaire et al. [21] and Kapila et al. [20], will be used to investigate the cavitation phenomenon in this study. The performance of these models will be evaluated using four underwater explosion (UNDEX) test cases in the open-source finite-volume UCNS3D solver [22].
To solve the interface capturing models, high-order central-weighted essentially non-oscillatory (CWENO) schemes for multi-component species within the finite volume framework, as implemented by Tsoutsanis and Dumbser[23], will be employed. The CWENO schemes offer computational advantages such as lower cost and faster computation compared to classical WENO schemes, due to the use of low-order polynomials associated with directional stencils and their reduced size compared to the high-order polynomials related to the central stencil. These advantages will be fully utilised in capturing the dynamics of cavitation in the test cases considered herein.
The paper proceeds to Section 2, where it presents the governing equations that describe the diffuse interface models of Allaire et al. [21] and Kapila et al. [20], which consist of five equations. It includes details about the additional K ∇ · v term that distinguishes these models, as well as the cut-off model used to prevent non-physical pressure values during cavitation. In Section 3, the implementation of the reconstruction process of the CWENO schemes is described, including the chosen fluxes and temporal discretisation, within the finite volume framework using the diffuse interface methods mentioned earlier. Section 4 is dedicated to the study of five test cases, where the proposed method is utilised to investigate cavitation occurrence in underwater explosions near different boundaries, such as a planar solid wall, a closed cylindrical container, and a free surface. The analysis focuses on parameters such as the influence of shock loading and cavitation collapse reloading on nearby structures like ship hulls or closed tanks. Additionally, the bubble dynamics of cavitating flows induced by underwater explosions near a free surface, as well as the impact of bubble separation distance on cavitation phenomena, are analysed. The numerical results obtained are compared with analytical, experimental, and reference solutions whenever possible. A comparison between the five-equation models of Allaire et al. [21] and Kapila et al. [20] is also presented. Finally, the last section provides concluding remarks and a summary drawn from this study.

2. The Five-equations Diffuse Interface Model

The five-equations model typically consists of two continuity equations, one momentum equation and one energy equation coupled with a non-conservative advection equation of the volume fraction equation. The common five-equations models are the Allaire et al. [21], Kapila et al. [20], Murrone et al.[24], and Wackers et al. [25] models, all these are simpler formulations or reductions to lesser unknowns of the full seven (7) Baer-Nunziato [13] equations, which are often too complex to solve numerically. In this paper, we shall only be considering the Allaire and Kapila five-equations models. The difference between Allaire’s and Kapila’s model is only in the non-conservative advection equation in the volume-fraction equation. The equations are presented below :
∂ ( α 1 ρ 1 ) ∂ t + ∇ · ( α 1 ρ 1 v ) = 0 ,
∂ ( α 2 ρ 2 ) ∂ t + ∇ · ( α 2 ρ 2 v ) = 0 ,
∂ ρ v ∂ t + ∇ · ( ρ v v + p I ) = 0 ,
∂ ρ E ∂ t + ∇ · ( ρ E + p ) v = 0 ,
∂ α 1 ∂ t + v · ∇ α 1 = 0 ( Allaire et al . ) ,
∂ α 1 ∂ t + v · ∇ α 1 = K ∇ · v ( Kapila et al . )
where the subscript 1 , 2 represents the fluid’s component, α is the volume fraction of each component, ρ is the density, v = ( u , v , w ) is the velocity, p is the pressure and E is the total energy.
This class of five-equations models (Allaire et al. [21] and Kapila et al. [20] model) are often referred to as the mechanical equilibrium models. We assume that the momentum, energy, and mass transfer between the phases reach equilibrium due to the thermodynamic difference between each component. This reasonable flow assumption involves a free surface like the one treated in this paper. Note that this is a mechanical equilibrium of the fluid; each phase remains in thermal disequilibrium. All the models in this category can be expressed into the following general form of a non-linear system of PDE in multiple space dimensions:
∂ Q ∂ t + ∇ · F ( Q ) + H ( Q ) ∇ · v = S ( Q )
where Q is the vector of evolution variables (conserved and not conserved), F is a flux function, v is the velocity field, and H and S are non-conservative quantities.
We then define Q, F and H in the form as presented in Eq.(8) ( Allaire et al. model [21]) and Eq. (9) (Kapila et al. model [21]) for a two-phase flow. The K∇· v term in Kapila’s model [20] describes the compressibility and thermodynamic properties of the mixture region of the two-phase fluid. It is derived from the asymptotic reduction of Baer-Nunziato’s seven-equation model[13,26,27]. From Eq.(5) and Eq.(6), K ∇ · v term differentiates Kapila’s model et al.[20] from Allaire’s et al. [21].
The Allaire’s et al. 5-equation model presented in a vector form:
Q = α 1 ρ 1 α 2 ρ 2 ρ v ρ E α 1 , F = α 1 ρ 1 v α 2 ρ 2 v ρ v 2 + p v ( ρ E + p ) α 1 v H = 0 0 0 0 − α 1 .
The Kapila’s et al. 5-equation model presented in a vector form:
Q = α 1 ρ 1 α 2 ρ 2 ρ v ρ E α 1 , F = α 1 ρ 1 v α 2 ρ 2 v ρ v 2 + p v ( ρ E + p ) α 1 v H = 0 0 0 0 − α 1 − k .
For each fluid’s component i, the mixture is assumed to be in mechanical equilibrium as stated previously, the velocity and pressure of each component is v = v 1 ≡ v 2 and p = p 1 ≡ p 2 (i.e. we maintain a single velocity and single pressure for the mixture).
The total density, momentum, the kinetic and internal energy of the mixture is stated below:
ρ = ( α ρ ) 1 + ( α ρ ) 2 = ∑ i = 1 N ( α ρ ) i : mixture density
ρ v = ( ( α ρ ) 1 + ( α ρ ) 2 ) ) v = ∑ i = 1 N ( α ρ ) i v : mixture momentum
ρ e K = ρ i v 2 2 : mixture kinetic energy .
ρ e I = ( α ρ ) 1 e 1 + ( α ρ ) 2 e 2 = ∑ i = 1 N ( α ρ e ) i : mixture internal energy .
So the total energy can be expressed as:
ρ E T = ( ( α ρ ) 1 e 1 + ( α ρ ) 2 e 2 + ρ i v 2 2 )
To determine the internal energy and put closure to the equations, it is convenient to use the stiffened EOS for both fluids components:
p i = γ i − 1 ρ i e i − γ i p c , i , i = 1 , 2
Using the stiffened EOS, the internal energy can now be expressed as:
ρ e i = p + γ i p c , i γ i − 1
and the total energy below:
ρ E = ∑ i = 1 N α i ( p + γ i p c , i γ i − 1 + ρ i . v 2 2 )
From the five-equations, as shown in Eq. (5) and Eq. (6), the advection equation of the volume fraction is in a non-conservative form which will present some difficulties when solving our five-model equations. We simplify it into a conservative term using Johnsen and Colonius [28] approach.
First, let f represent a vector in three-dimension and is defined as, f = α 1 v = α 1 u , α 1 v , α 1 w T .
Divergence of f then is:
∇ · f = ∂ ∂ x α 1 u + ∂ ∂ y α 1 v + ∂ ∂ z α 1 w = α 1 ∂ u ∂ x + ∂ v ∂ y + ∂ w ∂ z + u ∂ α 1 ∂ x + v ∂ α 1 ∂ y + w ∂ α 1 ∂ z = α 1 ∇ · v + v · ∇ α 1
And then,
v · ∇ α 1 = ∇ · f − α 1 ∇ · v .
Substituting v · ∇ α 1 into the non-conservative advection equation of the volume fraction in Eq. (5) and Eq. (6) to obtain the quasi-conservative form, the equation for the advection of the volume fraction can be written as Eq. (18) and Eq. (19). For Allaire’s et al:
∂ α 1 ∂ t + ∇ · f = α 1 ∇ · v
For Kapila’s et al.
∂ α 1 ∂ t + ∇ · f = α 1 + K ∇ · v
The K function in Kapila’s et al. can be determined from the Eq. (20) where the mixture’s speed of sound is gotten using the Wood’s speed of sound as seen in Eq. (22)
K = α 1 α 2 1 ρ 1 c 1 2 − 1 ρ 2 c 2 2 ρ c 2
where c 1 , c 2 and ρ are given as :
c 1 = γ 1 ρ 1 p + p c , 1 , c 2 = γ 2 ρ 2 p + p c , 2
1 ρ c 2 = α 1 ρ 1 c 1 2 + α 2 ρ 2 c 2 2 = ∑ i = 1 N α i ρ i c i 2
The “mixture-mixture” speed of sound derived from the mixture stiffened-gas SG-EoS is used for the Allaire et al. five -equations [21] and is given as follows:
c = γ ( p + p ∞ ) ρ

2.1. Cut-off Methods

The stiffened-gas equation of state (SG-EOS) introduced by Harlow and Amsden [29] has been commonly used by various authors [28,30,31,32,33] to achieve thermodynamic closure for multiphase compressible flow, including flows involving cavitation. However, in cases where the pressure difference between the two-phase medium is considerably high, such as in underwater explosions, the use of stiffened gas EOS alone may not be sufficient. This is because cavitation induced by the high-pressure surge is very fast, unsteady, and rapidly evolves into various dimensions before collapsing violently, resulting in non-physical negative pressures when using only the stiffened gas EOS. Similar findings have been reported by other authors in their numerical simulations.
To improve the accuracy of simulations for cavitating flows, modifications are needed. One common approach is to incorporate an additional single-fluid model to circumvent the issue, such as the cut-off model, isentropic model, or Schmidt model, which have been used by authors in previous studies. Among these models, the cut-off model is often preferred due to its simplicity of implementation and preservation of peak pressure accuracy during cavitation collapse, which is important in underwater explosion (UNDEX) test cases. Additionally, these models are pure phase models that do not consider phase exchange, as noted in previous studies [4,34]. Further details of these models can be found in the references cited. In the cut-off model, once the pressure in the cavitation region drops below a certain fixed saturation pressure psat, the non-physical pressure computed assumes a new value which is the value set for the psat.
p i = S G − E O S , if p i > p sat p sat , if p i < p sat
The physical saturation pressure ( p sat ) is chosen to be a comparatively small value in relation to the significantly higher ambient pressure exerted by the fluid, and p i refers to the static pressure at cell i. To avoid the negative physical pressure, most authors [4,5,35,36] assume a value between 2000-5000 Pa.
The mixture’s internal energy can be expressed as :
ρ i e i I s a t = ( p s a t + γ i p c , i γ i − 1 ) ,
and the total energy is final given by:
ρ e sat = ∑ i = 1 N α i ( p s a t + γ i p c , i γ i − 1 + ρ i . v 2 2 ) .

3. Numerical model

3.1. Finite volume framework

Lets consider a 3D domain Ω consisting of conforming tetrahedral, hexahedral, prism, and pyramid, triangles, or quadrilateral cells each one of them indexed by a unique mono-index i, and the governing equations of the five-equations model written in vector form as follows:
∂ ∂ t ∫ Ω i Q d Ω + ∮ ∂ Ω i F n ( ( Q ) d S = ∫ Ω i G d Ω
where Q = Q ( x , t ) is the vector of conserved variables and the volume fraction of one species, F n is the non-linear flux in the direction normal to the cells interface as given below:
∫ Ω i A l l a i r e α 1 ∇ · v d Ω ≈ ∫ Ω i α 1 d Ω · ∫ ∂ Ω i ( v n ) R i e m . d S .
∫ Ω i K a p i l a ( α 1 + K ) ∇ · v d Ω ≈ ∫ Ω i ( α 1 + K ) d Ω · ∫ ∂ Ω i ( v n ) R i e m . d S .
where ( n x , n y , n z ) is a unit normal vector across the boundary δ Ω , v n is the velocity normal to the bounded surface area, defined by v n = n x u + n y v + n z w . The source term G is with regards to the term α 1 ∇ · v of Eq. (18). Following the approach of Johnsen and Colonius [28] the source term is numerically approximated as a surface integral, rather than a volume one while using the same velocity estimate as the one used for the evaluation of the fluxes. The numerical flux F n ( Q ) can be computed considering the states to the left and the right of the high-order approximations of the solutions for the cell interface i and j respectively and expressed as:
F n ij = F n ij Q i j , L n ( x i j , α , t ) , Q i j , R n ( x i j , α , t )
Integrating Eq. (26) over the mesh element i using a high-order explicit finite-volume formulation, the following equation is obtained that incorporates the source term as previously defined:
d Q i d t = 1 V i ∑ j = 1 N f ∑ α = 1 N q p F n ij Q i j , L n x i j , α , t , Q i j , R n x i j , α , t − a i , 1 n · v n R i e m x i j , α , t ω α S i j ,
where Q i is the volume averaged vector of variables
Q i = 1 V i ∫ V i Q ( x , y , z ) d V ,
and F n ij is a numerical flux function in the direction normal to the cell interface between a considered cell i and one of its neighbouring cells j. N f is the number of faces per element, N q p is the number of quadrature points used for approximating the surface integrals, | S i j | is the surface area of the corresponding face, and Q i j , L n ( x i j , α , t ) and Q i j , R n ( x i j , α , t ) are the high-order approximations of the solutions for cell i and cell j respectively. α corresponds to different Gaussian integration points x α and weights ω α over each face. a i , 1 n corresponds to the volume-averaged volume fraction of cell i at time level n. A suitable Gauss-Legendre quadrature numerically approximates the volume, surface and line integrals.
The reconstruction process adopted in UCNS3D [22] follows the approaches of Tsoutsanis et al. [23,37,38,39,40,41,42,43], Titarev et al.[38] that have been previously applied to a wide-range of smooth and discontinuous flow compressible flow problems including laminar, transitional, and fully turbulent settings [44,45,46,47,48,49,50,51,52,53,54,55,56,57,58].

3.2. CWENO scheme

Tsoutsanis’s and Dumbser’s [23] CWENO and CWENOZ schemes for mixed-element unstructured meshes schemes were used to test compressible multiphase components [42,59,60]. The paper shows that the CWENO can produce lesser oscillations near the interface at a cheaper computational cost than the WENO schemes on the same grid. This is made possible by the fact that in the presence of discontinuity, CWENO, unlike the traditional WENO schemes, uses the lower-order polynomials from the directional stencils that might contain a smooth solution for its reconstruction, and in the presence of smooth use, the optimal polynomial in the desired order of accuracy. This slight but significant difference gives CWENO scheme advantages over the WENO. The CWENO will be discussed below briefly, and consequently, the same schemes will be used in capturing and resolving the material interface between the two-phase medium of UNDEX. More details of the CWENO schemes implemented by Tsoutsanis on an unstructured grid can be seen in the paper [23,32].
The optimal polynomial are defined as:
p o p t ( ξ , η , ζ ) = ∑ s = 1 s t λ s p s ( ξ , η , ζ ) ,
where s is the stencil index, with c = 1 being the central, c = 2 , 3 , . being the directional, s t being the total number of stencils, and λ s being the linear coefficients for each stencil, whose sum is equal to 1. The p 1 polynomial is not computed directly, but computed by subtracting the lower-order polynomials from the optimum polynomial as follows:
p 1 ( ξ , η , ζ ) = 1 λ 1 p o p t ( ξ , η , ζ ) − ∑ s = 2 s t λ s p s ( ξ , η , ζ ) .
The CWENO computes a reconstruction polynomial below as a non-linear combination of all the polynomials:
p ( ξ , η , ζ ) c w e n o = ∑ s = 1 s t ω s p s ( ξ , η , ζ ) ,
where ω s correspond to the non-linear weights assigned to each polynomial, and in regions with smooth data ω s ≈ λ s , hence obtaining the high-order approximation from the central stencil, and in regions of discontinuous solutions the reconstructed solution will be mostly influenced from the lower-order polynomials of the directional stencils using the Type 3 directional stencils. The reader is referred to the work of Tsoutsanis [61] for the stencil selection algorithms. a ˜ k are the reconstructed degrees of freedom; and the non-linear weight ω s is defined as: the nonlinear coefficients ω s are computed from the linear ones ω s as:
ω s = ω ˜ s ∑ s = 1 s t ω ˜ s where ω ˜ s = λ m ( ϵ + S I s ) b .
where ϵ is taken as 10-6 similarly to the WENO reconstructions, S I s denotes the smoothness indicator and b = 4 . The readers are reffered to [23] for a detailed description of the smoothness indicators.

3.3. Fluxes & Time advancement

For the numerical fluxes we employ the (Harten-Lax-van Leer-Contact) Riemann solver of Toro [62]. The solution is advanced in time by a third-order TVD Runge-Kutta method:
U i ( n + 1 / 3 ) = U i n + Δ t R i ( U ) ,
U i n + 2 / 3 = 3 4 U i n + 1 4 U i ( n + 1 / 3 ) + 1 4 Δ t R i ( U ( n + 1 / 3 ) ,
U i n + 1 = 1 3 U i n + 2 3 U i ( n + 2 / 3 ) + 2 3 Δ t R i ( U ( n + 2 / 3 ) ,
and the time step Δ t is selected according to:
Δ t = CFL 1 d min i ( h i S i ) ,
where S i is an estimate of the maximum in absolute value of the propagation speed in the cell i, d is the number of dimensions, h i is a characteristic length of the element i and CFL refers to the Courant-Friedrichs-Lewy condition, and a C F L = 0.6 is used for all tests unless otherwise specified. The reconstruction is performed with respect to the primitive variables since they were previously found [59] to be less oscillatory than the conserved variables.
All the volume/surface/line integrals are approximated by Gaussian quadrature rule suitable for the order of polynomial employed. All the schemes developed are implemented in the UCNS3D CFD code [22] which is written in object-oriented Fortran 2008, employing MPI message passing interface (MPI), and the Open Multi-Processing (OpenMP) application programming interface (API), and the reader is referred to [63] for more details on implementation and performance benchmarks.

4. Test Cases

Using the high-order CWENO schemes in conjunction with the 5-equation interface capturing schemes models [20,21], five underwater test problems were selected to study the cavitation effect of underwater shock near different boundaries.
1.
2D Underwater explosion near a planar wall.
2.
2D Underwater explosion near a free surface.
3.
3D Underwater explosions in a rigid cylinder.
4.
2D underwater explosion of two bubbles placed horizontally near a free surface.
5.
3D Underwater of two bubbles placed vertically in an enclosed region.

4.1. 2D underwater explosion near a free surface

To study the dynamics of bubbles at the interface between air and water in an underwater explosion, a two-dimensional model of an underwater explosion near the free surface of water is considered. Previous studies by various authors have used numerical [5,15,64,65,66,67,68,69] and experimental [69,70] methods to investigate this phenomenon. For the experimental setup, a spherical charge of silver azide with known initial energy is detonated in an open cylinder half-filled with water, and the process involves charge detonation, bubble pulsation, and jet formation [69].
For all of our UNDEX numerical studies in this paper, a pressurised bubble will be used as a substitute for explosives in underwater explosions. The pressurised gas bubble with a diameter of 0.24 m is placed at [0, -0.3] m in the two-phase medium, using a similar setup as seen in Shukla et al.’s work [68]. Our computational domain is a square area with x × y of dimensions ( [ − 2 , 2 ] × [ − 1.5 , 2.5 ] ) m as shown in Figure 1, with the air-water interface positioned at the centre [ 0 , 0 ] .
Transparent boundary conditions were set up at the top, left, and right surfaces, while the bottom is set to reflecting. The simulation is carried out until t = 1.2 ms as the final time, and the initial physical parameters of each fluid (density, velocities, pressure, specific ratio, and volume fraction) were set accordingly and presented below:
ρ , u , v , p , γ , α = 1.225 , 0 , 0 , 10 5 , 2.0 , 0 , Air 1250 , 0 , 0 , 10 9 , 2.0 , 1 , Explosive 1000 , 0 , 0 , 10 5 , 7.15 , 0 , Water .
We employ the high-order CWENO4 finite-volume method developed by Tsoutsanis et al. [23,32], on three grid resolutions (coarse, medium and fine) in our study to investigate the interaction between gas bubbles and the free surface of an air-water medium. Numerical Schlieren visualisations were generated to capture density and pressure gradients at various instants, as depicted in Figure 3 and Figure 4, respectively. The high-pressure air bubble undergoes an explosion, resulting in a shock wave that propagates outward towards the interface, reflects into the water upon hitting the air-water interface, and leaves behind cavitation bubbles, as evident from the pressure gradients in Figure 4. Subsequently, the bubble transforms from a spherical shape to an oval shape, as illustrated in the volume fraction contour in Figure 5.
We conducted a study to investigate the impact of grid resolution on the evolution of the interfacial region between the exploding gas bubble and the air-water free surface using the CWENO4 scheme, employing different grid sizes (400 × 400, 800 × 800, and 1600 × 1600). Changes in the CWENO4 scheme’s behaviour and its ability to capture the interfacial region between the gas bubble and the air-water free surface can be seen in Figure 6 as the grid size is refined. These include reductions in the smearing or numerical diffusion between the free surface and the expanding air bubble, resulting in a sharper description of the interface. The results are free from oscillations which also demonstrates the robustness of the numerical framework employed. However interface sharpening or anti-diffusion techniques could prove beneficial especially for up to 3rd-order methods at coarse grid resolutions since they are the ones that can benefit most from this strategy [35,59,68].
We have performed a quantitative investigation on the medium grid, for the obtained density and pressure profiles along the horizontal axis (x = 0) and compared them with published results that utilize sharpening techniques and anti-diffusive techniques in conjunction with high-order schemes [9,71]. The results, as shown in Figure 7 and Figure 8, are consistent with previously published findings[9,71,72].
Furthermore, we applied both the five-equation models proposed by Allaire et al. and Kapila et al. to investigate the influence of the additional term K ∇ · u in Kapila’s model, thereby differentiating between the two models. The density and pressure plots, depicted in Fig. Figure 9, demonstrate that both the five-equation diffuse-interface multi-component models of Allaire et al. and Kapila et al. yield comparable results, and are in good agreement with the published works of [73] and [71] . This shows that, for cavitation test cases such as the one under consideration in this paper, we can affirm that 5-equation model of Allaire et al. and Kapila et al. both accurately portrays cavitation formation and collapse. However, this accuracy does not extend to spherical bubble collapse in a liquid, as highlighted in the paper [26]. where the 5-equation model of Allaire et al. falls short in correctly representing the solution of spherical bubble collapse. In Figure 10, we observe that the maximum peak pressure on the finest mesh is approximately 8.5 MPa at x = -1.0 m, which is the location where the shock wave reflects upon hitting the expanding cavity from below. This finding is consistent with the observation made in [73] that traditional WENO schemes, although computationally more expensive than CWENO, fail to prevent the smearing of the air-water interfaces. Our results demonstrate that CWENO is capable of resolving material interfaces with reduced numerical smearing while maintaining their non-oscillatory nature, as reported in [59]. This makes CWENO schemes a favourable choice in terms of computational cost for simulating compressible multiphase or multicomponent flows involving cavitation.
Preprints 94023 g007
Figure 7. UNDEX near a free surface; plot of the density distribution along the horizontal axis x = 0 at two instants, specifically t = 0.4 ms and t = 0.8 ms. The obtained results are compared to those reported in previous studies of Kokh (2010) [71] and Shyue (2006)[74].
Figure 7. UNDEX near a free surface; plot of the density distribution along the horizontal axis x = 0 at two instants, specifically t = 0.4 ms and t = 0.8 ms. The obtained results are compared to those reported in previous studies of Kokh (2010) [71] and Shyue (2006)[74].
Preprints 94023 g007
Preprints 94023 g008
Figure 8. UNDEX near a free surface; pressure plot along horizontal axis x = 0 at instants t = (0.4; 0.8)ms. Comparison are made with the results of [71,74] and [72]
Figure 8. UNDEX near a free surface; pressure plot along horizontal axis x = 0 at instants t = (0.4; 0.8)ms. Comparison are made with the results of [71,74] and [72]
Preprints 94023 g008
Preprints 94023 g009
Figure 9. UNDEX near a free surface; density and pressure plot along horizontal axis x = 0 at instants t = 0.2ms using Allaire et al. [21] and Kapila et al. [20] models.
Figure 9. UNDEX near a free surface; density and pressure plot along horizontal axis x = 0 at instants t = 0.2ms using Allaire et al. [21] and Kapila et al. [20] models.
Preprints 94023 g009
Preprints 94023 g010
Figure 10. UNDEX near a free surface; density and pressure plot along horizontal axis x = 0 at instants t = 1.2 ms using different grid resolutions.
Figure 10. UNDEX near a free surface; density and pressure plot along horizontal axis x = 0 at instants t = 1.2 ms using different grid resolutions.
Preprints 94023 g010

4.2. 2D cylindrical underwater explosion near a planar rigid wall.

This test case is carried out to examine the impact of shock loading and cavitation loading or reloading, caused by an underwater explosion, on a rigid structure (e.g., a ship’s hull). Numerous studies, both numerical and experimental, have been conducted on this test case by various researchers [34,65,69,75,76]. The experiment involves detonating an explosive submerged in water underneath a plate.
The physical properties (density, velocities, pressure, specific ratio, and volume fraction) of each fluid involved in these test cases are presented below:
ρ , u , v , p , γ , α = 1270 , 0 , 0 , 8.29 × 10 8 , 2 , 1 , Explosive 1000 , 0 , 0 , 10 5 , 7.15 , 0 , Water .
The rectangular computational domain in two dimensions (as shown in Figure 11) has dimensions of ( [ − 6 , 6 ] × [ − 6 , 3 ] ) m. Similar to the previous test case, the explosive is modelled as a high-pressure air bubble with a spherical gas bubble of 1 m diameter placed in water at the origin [ 0 , 0 ] . A rigid wall is located at y = 3 m above the bubble, and a numerical sensor is placed at the centre point of the wall (denoted as P) to measure the load on the wall. The boundary conditions for the upper surface are set as reflecting, while the left, right, and bottom surfaces are set as non-reflecting. A coarse grid is used to discretise the domain, and the simulation is run for a final time of t = 100 ms.
The 4th-order CWENO scheme is used with the 5-equation diffuse interface model for this test problem. To address the issue of non-physical negative pressure values (i.e., p < 0) that can occur in cavitation regions during the explosion, a cut-off model, described in Section 2.1, is introduced to maintain the pressure at 5000 Pa.
Contour plots of pressure obtained using the 4th-order CWENO scheme at different times during the explosion are shown in Figure 12. The plots shows that the primary loading occurs when the pressurised air bubble explodes, generating shock waves that hit the wall above it at t = 1.5 ms and immediately reflect into the water, causing a rapid drop in density. Rarefaction can be observed in Figure 12 at t = 2 ms, as the pressure decreases around the wall, forming a low-pressure region at t = 3 ms. This pressure drop intensifies, resulting in a cavitation region around the wall at t = 4 ms and t = 5.3 ms. The cavitation collapses and generates a reload impact on the wall, known as cavitation collapse reloading (secondary reloading).
The pressure history of the wall along the horizontal axis is shown in Figure 13. The first peak pressure created by the shockwaves is 6200 MPa for the Allaire et al. model and 6400 MPa for the Kapila et al. model, indicating good agreement with the results of other authors [65,75]. As the bubbles continue to expand, the cavitation reloading impact on the wall occurs at a second peak pressure of 720 MPa for Allaire et al.’s model and 760 MPa for Kapila et al.’s model. It should be noted that although the impact of the shock loading is higher than that of the cavitation reloading, the duration of the cavitation reloading is longer, resulting in significant loading on the wall.

4.3. 3D underwater explosion in an enclosed, rigid cylindrical wall.

This test case has been investigated numerically by various authors [9,34,75,78,79,80]. It serves as an ideal problem for studying shock and cavitation responses in enclosed structures such as a cylindrical tank. The initial conditions for this simulation are based on Xie et al.’s work [4] and are presented below:
ρ , u , v , p , γ , α = 1770 , 0 , 0 , 2 × 10 9 , 2 , 1 , Explosive 1000 , 0 , 0 , 10 5 , 7.15 , 0 , Water .
The computational setup is illustrated in Figure 14, following that of [4,75].
The explosive, modeled as a gas bubble, has a diameter of 0.1 m and is placed within a rectangular cylindrical wall with dimensions of 0.2286 m in height and 0.0889 m in width. The water in the tank is partially aerated with 15% air. The test case can be simulated in two-dimensional axisymmetric form using cylindrical coordinates, or in three-dimensional form. In our study, the test is conducted in a three-dimensional setup.
Reflecting boundary conditions are imposed on all surfaces, and a hybrid hexa-dominant mesh is used to discretise the computational domain with approximately 0.5 million cells as shown in Figure 16. A CWENO4 scheme is utilized for numerical simulations, and the simulation is carried out until t = 120 µs. The time evolution of the gas bubble’s pressure is shown in the pressure contour plot in Figure 17. It can be observed that as the bubble explodes and expands, shock waves are generated and impact the side walls of the cylinder at t = 15 µs. These waves are immediately reflected back into the water, striking the expanding bubble at t = 30 µs. This results in rapid changes in density and rarefaction. At t = 45 µs, the shape of the bubble transforms into an oval. Over time, a low-pressure region is observed around the bubble and close to the center of the lateral walls at t = 30 µs, leading to cavitation at t = 45 µs, 60 µs, and 70 µs, as evident in Figure 17. The cavitation pocket is most prominent on the wall at t = 60 µs. The violent collapse of cavitation generates reloading pressure on the walls of the cylindrical structure, as depicted in Figure 15. This phenomenon of cavitation-induced reloading can have a significant impact on enclosed structures, with effects comparable to those of the primary shock-wave loading.
The pressure history plot in Figure 15 reveals that the initial peak pressure obtained in our study using the diffuse interface methods proposed by Allaire et al. is 6600 MPa. This value is in close agreement with the experimental result of Wardlaw et al. [83], where P = 6300 MPa. Other studies by Shyue et al. [9] and Ying et al. [82], which utilized the Allaire et al. diffuse interface methods in conjunction with cut-off cavitation models, reported numerical values in the range of 7000-7200 MPa, respectively.
The second impact of cavitation collapse occurs at a time of 90 µs, with a pressure value of p = 2100 MPa, as evident in our results. This value is lower than the numerical results reported by Shyue et al. [9] and Ying et al. [82], but closer to the experimental value than the latter. Notably, the study conducted by Chen [84] using Smoothed Particle Hydrodynamics (SPH) models exhibited a delay in the cavitation collapse and reload, as the Schmidt model [4,77,85] was employed instead of the cut-off model, as utilized in the present study.

4.4. 2D underwater explosion of two bubbles placed horizontally near a free surface

When two bubbles are positioned close to one another, the local pressure field may vary due to the bubble-bubble interactions. The growth, collapse, and oscillation periods of the bubbles—i.e., the time it takes for a bubble to complete one cycle of expansion and collapse can be affected by these changes in pressure. In a two-dimensional setup Figure 18, two bubbles are positioned horizontally in close proximity to each other near a free surface as an attempt to investigate the influence of the separation distance between the bubbles on cavitation. The separation distance, denoted as L o , is varied at four different values: ( 0.3 , 0.5 , 0.7 , 0.9 ) m, between two identical bubbles located 0.3 m below the free surface. This configuration closely is similar to the one utilized by Yu et al.[86]. Additionally, a single bubble with similar computational parameters and dimensions is also simulated for comparative analysis. The domain is discretized using a quadrilateral coarse mesh comprising approximately 0.025 million cells, with all surfaces imposing outflow boundary conditions.
Preprints 94023 g018
Figure 18. Computational domain setup for the two-dimensional underwater explosion of two bubbles placed horizontally near a free surface.
Figure 18. Computational domain setup for the two-dimensional underwater explosion of two bubbles placed horizontally near a free surface.
Preprints 94023 g018
Preprints 94023 g019
Figure 19. Underwater explosion of a single and two bubbles configurations placed horizontally near a free surface; plots of pressure history (at the centre location of the left wall ) for the underwater explosion in an enclosed cylindrical wall with CWENO4 scheme, different inter-bubble distances ( 0.3 , 0.5 , 0.7 , 0.9 ) m are considered. [9,81,87]
Figure 19. Underwater explosion of a single and two bubbles configurations placed horizontally near a free surface; plots of pressure history (at the centre location of the left wall ) for the underwater explosion in an enclosed cylindrical wall with CWENO4 scheme, different inter-bubble distances ( 0.3 , 0.5 , 0.7 , 0.9 ) m are considered. [9,81,87]
Preprints 94023 g019
The physical characteristics of the fluid involved in these test cases are described below, which includes parameters such as density, velocity, pressure, specific ratio (or specific heat capacity ratio), and volume fraction:
ρ , u , v , p , γ , α = 1270 , 0 , 0 , 8.29 × 10 8 , 2 , 1 , Explosive 1000 , 0 , 0 , 10 5 , 7.15 , 0 , Water .
The volume fraction contours presented in Figure 20 effectively illustrate the temporal evolution of the bubbles with the distances L o ranging from 0.3 m to 0.9 m. As the separation distance L o , increases from 0.3 m to 0.9 m, the height of the free surface center consistently decreases at each time instant. In the same figure, it is observed that for smallest inter-bubble distances L o = 0.3 m, the bubbles expand and eventually coalesce or merge.
As observed in previous UNDEX (Underwater Explosion) test cases in this paper, the two bubbles undergo rapid expansion, generating a shock wave that propagates through the water and interacts with the surrounding bubbles. The pressure-time history at the center of the explosion, depicted in the pressure contour plot Figure 19, reveals that the peak pressure caused by the initial shock wave hitting the free surface for a single bubble is 170 MPa, while for the two-bubble case with the largest inter-bubble distance of L o = 0.9 m, it reaches 290 MPa, which is 1.7 times higher than that of the single bubble. Furthermore, as evident in Figure 19, the initial peak pressure which is as a result of the intensity of the shock wave on the free surface increases with decreasing inter-bubble distance, from the maximum value at L o = 0.9 m to the smallest value at L o = 0.3 m. The close proximity of the two bubbles allows for stronger interaction, leading to a higher shock-intensity on the free surface.
The occurrence of cavitation beneath the free surface during a two-bubble explosion is influenced by the size and spacing of the bubbles. In same plot Figure 19, as the inter-distance bubbles increase, the cavitation bubble collapse that results in the cavitation reloading on the free surface is further delayed as seen in the second pressure peak. The second peak pressure associated with cavitation reloading is also more pronounced for the two charged bubbles compared to a single bubble. This observation further supports the conception that two charged bubbles have a higher potential for intensifying structural damage due to cavitation reloading compared to a single bubble. The ratio of the initial peak pressure on the free surface caused by the shock wave to that of the cavitation reloading is about four times for the smaller inter-bubble L o = ( 0.3 , 0.5 ) m, and twice for the larger inter-bubble distances L o = ( 0.7 , 0.9 ) m. This shows that for larger inter-bubble distances, the effect of structural damage caused by cavitation is more comparable to that of the primary shock-wave loading than for smaller inter-bubble distances.
In Figure 21, the pressure gradient contour is depicted for that of a single charge and the two bubbles with varying inter-bubble distances. The contour illustrates a region of higher pressure near bubble(s), while the pressure diminishes further away as the bubbles expand. The magnitude of the pressure gradients varies based on the distance between the bubbles.

4.5. 3D Underwater of two bubbles placed vertically in an enclosed region.

To investigate the interaction between two vertically aligned bubbles in water, enclosed in a region, it is best to consider this test case in either an axis-symmetric or three-dimensional setup, taking into account gravitational effects. The bubbles used in this test case are of smaller and identical diameter (d = 5.3 mm). They are initially positioned one above the other, with the first bubble placed at the bottom of the enclosed region and the second bubble above it.
Our numerical setup (Figure 22) uses the similar computational dimensions as the Test case Section 4.3, but with smaller bubble radii. The distance between the first and second bubbles is denoted by H b , and chosen to be 15, 30, and 60, all of which are given in millimetres. Reflecting boundary conditions are applied to all surfaces, and a quadrilateral mesh discretizes the domain with approximately 0.4 million cells. The initial physical properties of both fluid is given below:
ρ , u , v , p , γ , α = 1.225 , 0 , 0 , 1.57 × 10 7 , 2.0 , Air bubble 998.2 , 0 , 0 , 10 5 , 7.15 , 0 , Water .
The volume fraction contour in Figure 24 displays the shortest inter bubble distance H b = 15 mm, where the bubbles are in close proximity. As the explosion evolves with time, the water film between the interface of the two interacting gaseous bubbles thins out. This causes a surface tension force to pull the bubbles together until the thin film ruptures and the two bubbles merge together. The two bubbles combine to form a single, larger bubble after coalescing.
Preprints 94023 g023
Figure 23. Underwater explosion of two bubbles vertically placed in an enclosed region; plots of pressure history (at the centre location of the left wall ) for the underwater explosion in an enclosed cylindrical wall with CWENO4 schemes for inter-bubble lengths of (15, 30, and 60) mm
Figure 23. Underwater explosion of two bubbles vertically placed in an enclosed region; plots of pressure history (at the centre location of the left wall ) for the underwater explosion in an enclosed cylindrical wall with CWENO4 schemes for inter-bubble lengths of (15, 30, and 60) mm
Preprints 94023 g023
Preprints 94023 g024
Figure 24. Underwater explosion of two bubbles vertically placed in an enclosed region; volume fraction contours for the 3D underwater explosion in an enclosed for H b = 30 mm, rigid cylindrical wall.Captured using CWENO4. The green-yellow areas shows the interface between water and explosive. The two bubbles eventually merge into a single larger bubble.
Figure 24. Underwater explosion of two bubbles vertically placed in an enclosed region; volume fraction contours for the 3D underwater explosion in an enclosed for H b = 30 mm, rigid cylindrical wall.Captured using CWENO4. The green-yellow areas shows the interface between water and explosive. The two bubbles eventually merge into a single larger bubble.
Preprints 94023 g024
The pressure plot history at the centre of the left wall is displayed in Figure 23 for the three inter-bubble instances ( H b = 15, 30, and 60 mm). The pressure gradient contours in Figure 25 shows the cavitation regions are more prominent for the shortest inter bubble distance H b = 15 mm for varying time frames t = (1.54, 6, 30, 60 and 90) µs. The cavitation regions in same Figure 25 can be observed to be majorly concentrated at the walls of the cylinder and the surrounding areas close to the bubbles.

5. Conclusions

This research extends the CWENO high-order finite-volume numerical framework employed for simulating multicomponent and multiphase flows, as detailed in [59], to specifically investigate cavitation phenomena in underwater explosions (UNDEX). The application of the CWENO schemes proved successful across a range of intricate two- and three-dimensional test cases related to underwater explosions, yielding noteworthy observations.
Notably, the CWENO schemes showcase a remarkable capability to provide high-resolution and sharply defined representations of shock-cavitation structures, as well as the infinitesimal interfacial region. This is achieved by resolving material interfaces with minimized numerical smearing, eliminating the necessity for additional interface sharpening [68] or compression [71]. This takes away the increase in computational cost associated with incorporating interface sharpening or compression into the diffuse interface models.
The compressibility effect introduced by the K ∇ · u term in the Kapila et al. five-equation models [20], which distinguishes it from the Allaire five-equation models [21] from our results shows that this additional term does not significantly impact cavitation. The similarity in results between the two models underscores their agreement with other published papers, affirming their reliability in simulating cavitation phenomena.
The results also show that the proximity of the bubbles played a crucial role in intensifying the shock wave impacts and the occurrence of cavitation around free surfaces and rigid boundaries. All these findings contribute valuable insights into dynamics of cavitation within the context of underwater explosions.

Author Contributions

Conceptualization, E.A. and P.T.; methodology, E.A. and P.T.; software, E.A. and P.T.; validation, E.A. and P.T.; formal analysis, E.A. and P.T.; investigation, E.A. and P.T.; resources, E.A., P.T and K.J; data curation, E.A.; writing—original draft preparation, E.A. and P.T.; writing—review and editing, E.A. and P.T.; visualization, E.A. and P.T.; supervision, K.J. and P.T.; project administration, K.J. and P.T.; funding acquisition, E.M. All authors have read and agreed to the published version of the manuscript

Funding

This research is partly financially supported by the The Petroleum Trust Development Fund (PTDF) Nigeria through grant number PTDF/ED/OSS/PHD/DBB/1558/19 https://ptdf.gov.ng for the first author’s PhD studies.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Acknowledgments

The authors acknowledge the computing time at Cranfield University Delta2 HPC facility.

Conflicts of Interest

The authors declare that they have no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
UNDEX Underwater Explosion
DIM Diffuse Interface Models
UCNS3D Unstructured Compressible Navier Stokes 3D
CWENO Central Weighted Essentially Non-Oscillatory scheme
FV Finite Volume

References

  1. Petitpas, F.; Massoni, J.; Saurel, R.; Lapebie, E.; Munier, L. Diffuse interface model for high speed cavitating underwater systems. International Journal of Multiphase Flow 2009, 35, 747–759. [Google Scholar] [CrossRef]
  2. Petitpas, F.; Saurel, R.; Ahn, B.K.; Ko, S. Modelling cavitating flow around underwater missiles. International Journal of Naval Architecture and Ocean Engineering 2011, 3, 263–273. [Google Scholar] [CrossRef]
  3. Kwack, Y.K.; Ko, S.H. Numerical analysis for supercavitating flows around axisymmetric cavitators. International Journal of Naval Architecture and Ocean Engineering 2013, 5, 325–332. [Google Scholar] [CrossRef]
  4. Wenfeng, X. A numerical simulation of underwater shock-cavitation-structure interaction. PhD Thesis 2005, 53, 1–207. [Google Scholar]
  5. Daramizadeh, A.; Ansari, M.R. Numerical simulation of underwater explosion near air-water free surface using a five-equation reduced model. Ocean Engineering 2015, 110, 25–35. [Google Scholar] [CrossRef]
  6. Xie, W.F.; Liu, T.G.; Khoo, B.C. The simulation of cavitating flows induced by underwater shock and free surface interaction. Applied Numerical Mathematics 2007, 57, 734–745. [Google Scholar] [CrossRef]
  7. Saurel, R.; Lemetayer, O. A multiphase model for compressible flows with interfaces, shocks, detonation waves and cavitation. Journal of Fluid Mechanics 2001, 431, 239–271. [Google Scholar] [CrossRef]
  8. Ghidaglia, J.M.; Mrabet, A.A. A regularized stiffened-gas equation of state. Journal of Applied Analysis & Computation 2018, 8, 675–689. [Google Scholar] [CrossRef]
  9. Shyue, K.M. An Anti-Diffusion based Eulerian Interface-Sharpening Algorithm for Compressible Two-Phase Flow with Cavitation. Proceedings of the 8th International Symposium on Cavitation 2012, 268, 7–12. [Google Scholar] [CrossRef]
  10. Maltsev, V.; Skote, M.; Tsoutsanis, P. High-order methods for diffuse-interface models in compressible multi-medium flows: A review. Physics of Fluids 2022, 34. [Google Scholar] [CrossRef]
  11. Pelanti, M.; Shyue, K.M. A mixture-energy-consistent six-equation two-phase numerical model for fluids with interfaces, cavitation and evaporation waves. Journal of Computational Physics 2014, 259, 331–357. [Google Scholar] [CrossRef]
  12. Pelanti, M.; Shyue, K.M. A numerical model for multiphase liquid–vapor–gas flows with interfaces and cavitation. International Journal of Multiphase Flow 2019, 113, 208–230. [Google Scholar] [CrossRef]
  13. Baer, M.; Nunziato, J. A two-phase mixture theory for the deflagration-to-detonation transition (ddt) in reactive granular materials. International Journal of Multiphase Flow 1986, 12, 861–889. [Google Scholar] [CrossRef]
  14. Saurel, R.; Petitpas, F.; Abgrall, R. Modelling phase transition in metastable liquids: application to cavitating and flashing flows. Journal of Fluid Mechanics 2008, 607, 313–350. [Google Scholar] [CrossRef]
  15. Jun, Y.; Jian-qiang, P.; Hai-kun, W.; Hai-bin, M. Application of Compressible Multi-component Flow in Underwater Explosion Problems. Procedia Engineering 2015, 126, 339–343. [Google Scholar] [CrossRef]
  16. Chiapolino, A.; Saurel, R.; Nkonga, B. Sharpening diffuse interfaces with compressible fluids on unstructured meshes. Journal of Computational Physics 2017, 340, 389–417. [Google Scholar] [CrossRef]
  17. LeMartelot, S.; Nkonga, B.; Saurel, R. Liquid and liquid–gas flows at all speeds. Journal of Computational Physics 2013, 255, 53–82. [Google Scholar] [CrossRef]
  18. Zein, A. Numerical methods for multiphase mixture conservation laws with phase transition. PhD Thesis 2010, 120. [Google Scholar]
  19. Ma, Z.H.; Causon, D.M.; Qian, L.; Gu, H.B.; Mingham, C.G.; Martínez Ferrer, P. A GPU based compressible multiphase hydrocode for modelling violent hydrodynamic impact problems. Computers and Fluids 2015, 120, 1–23. [Google Scholar] [CrossRef]
  20. Kapila, A.K.; Menikoff, R.; Bdzil, J.B.; Son, S.F.; Stewart, D.S. Two-phase modeling of deflagration-to-detonation transition in granular materials: Reduced equations. Physics of Fluids 2001, 13, 3002–3024. [Google Scholar] [CrossRef]
  21. Allaire, G.; Clerc, S.; Kokh, S. A five-equation model for the simulation of interfaces between compressible fluids. Journal of Computational Physics 2002, 181, 577–616. [Google Scholar] [CrossRef]
  22. Antoniadis, A.; Drikakis, D.; Farmakis, P.; Fu, L.; Kokkinakis, I.; Nogueira, X.; Silva, P.; Skote, M.; Titarev, V.; Tsoutsanis, P. UCNS3D: An open-source high-order finite-volume unstructured CFD solver. Computer Physics Communications 2022, 279. [Google Scholar] [CrossRef]
  23. Tsoutsanis, P.; Dumbser, M. Arbitrary high order central non-oscillatory schemes on mixed-element unstructured meshes. Computers and Fluids 2021, 225. [Google Scholar] [CrossRef]
  24. Murrone, A.; Guillard, H. A five equation reduced model for compressible two phase flow problems. Journal of Computational Physics 2005, 202, 664–698. [Google Scholar] [CrossRef]
  25. Wackers, J.; Koren, B. A fully conservative model for compressible two-fluid flow. International Journal for Numerical Methods in Fluids 2005, 47, 1337–1343. [Google Scholar] [CrossRef]
  26. Schmidmayer, K.; Bryngelson, S.H.; Colonius, T. An assessment of multicomponent flow models and interface capturing schemes for spherical bubble dynamics. Journal of Computational Physics 2020, 402, 109080. [Google Scholar] [CrossRef]
  27. Zhang, J. A simple and effective five-equation two-phase numerical model for liquid-vapor phase transition in cavitating flows. International Journal of Multiphase Flow 2020, 132, 1–47. [Google Scholar] [CrossRef]
  28. Johnsen, E.; Colonius, T. Implementation of WENO schemes in compressible multicomponent flow problems. Journal of Computational Physics 2006, 219, 715–732. [Google Scholar] [CrossRef]
  29. Harlow, F.H.; Amsden, A.A. A numerical fluid dynamics calculation method for all flow speeds. Journal of Computational Physics 1971, 8, 197–213. [Google Scholar] [CrossRef]
  30. Saurel, R.; Abgrall, R. Simple method for compressible multifluid flows. SIAM Journal of Scientific Computing 1999, 21, 1115–1145. [Google Scholar] [CrossRef]
  31. Quirk, J.J.; Karni, S. On the dynamics of a shock-bubble interaction. Journal of Fluid Mechanics 1996, 318, 129–163. [Google Scholar] [CrossRef]
  32. Tsoutsanis, P.; Adebayo, E.; Merino, A.; Arjona, A.; Skote, M. CWENO Finite-Volume Interface Capturing Schemes for Multicomponent Flows Using Unstructured Meshes. Journal of Scientific Computing 2021, 89. [Google Scholar] [CrossRef]
  33. Schmidmayer, K.; Petitpas, F.; Daniel, E.; Favrie, N.; Gavrilyuk, S. A model and numerical method for compressible flows with capillary effects. Journal of Computational Physics 2017, 334, 468–496. [Google Scholar] [CrossRef]
  34. Liu, T.; Khoo, B.; Xie, W. Isentropic one-fluid modelling of unsteady cavitating flow. Journal of Computational Physics 2004, 201, 80–108. [Google Scholar] [CrossRef]
  35. Ansari, M.; Daramizadeh, A. Numerical simulation of compressible two-phase flow using a diffuse interface method. International Journal of Heat and Fluid Flow 2013, 42, 209–223. [Google Scholar] [CrossRef]
  36. Zhang, J. A simple and effective five-equation two-phase numerical model for liquid-vapor phase transition in cavitating flows. International Journal of Multiphase Flow 2020, 132, 1–47. [Google Scholar] [CrossRef]
  37. Tsoutsanis, P.; Titarev, V.; Drikakis, D. WENO schemes on arbitrary mixed-element unstructured meshes in three space dimensions. Journal of Computational Physics 2011, 230, 1585–1601. [Google Scholar] [CrossRef]
  38. Titarev, V.; Tsoutsanis, P.; Drikakis, D. WENO schemes for mixed-elementunstructured meshes. Communications in Computational Physics 2010, 8, 585–609. [Google Scholar] [CrossRef]
  39. Tsoutsanis, P. Extended bounds limiter for high-order finite-volume schemes on unstructured meshes. Journal of Computational Physics 2018, 362, 69–94. [Google Scholar] [CrossRef]
  40. Simmonds, N.; Tsoutsanis, P.; Antoniadis, A.; Jenkins, K.; Gaylard, A. Low-Mach number treatment for Finite-Volume schemes on unstructured meshes. Applied Mathematics and Computation 2018, 336, 368–393. [Google Scholar] [CrossRef]
  41. Maltsev, V.; Yuan, D.; Jenkins, K.; Skote, M.; Tsoutsanis, P. Hybrid discontinuous Galerkin-finite volume techniques for compressible flows on unstructured meshes. Journal of Computational Physics 2023, 473. [Google Scholar] [CrossRef]
  42. Tsoutsanis, P.; Pavan Kumar, M.; Farmakis, P. A relaxed a posteriori MOOD algorithm for multicomponent compressible flows using high-order finite-volume methods on unstructured meshes. Applied Mathematics and Computation 2023, 437. [Google Scholar] [CrossRef]
  43. Farmakis, P.; Tsoutsanis, P.; Nogueira, X. WENO schemes on unstructured meshes using a relaxed a posteriori MOOD limiting approach. Computer Methods in Applied Mechanics and Engineering 2020, 363. [Google Scholar] [CrossRef]
  44. Antoniadis, A.; Tsoutsanis, P.; Drikakis, D. High-order schemes on mixed-element unstructured grids for aerodynamic flows. 42nd AIAA Fluid Dynamics Conference and Exhibit 2012 2012. [Google Scholar] [CrossRef]
  45. Kokkinakis, I.; Antoniadis, A.; Drikakis, D.; Tsoutsanis, P.; Rana, Z. High-order methods for hypersonic shock wave turbulent boundary layer interaction. 20th AIAA International Space Planes and Hypersonic Systems and Technologies Conference, 2015.
  46. Antoniadis, A.; Tsoutsanis, P.; Drikakis, D. Numerical accuracy in RANS computations of high-lift multi-element airfoil and aircraft configurations. 53rd AIAA Aerospace Sciences Meeting 2015. [Google Scholar] [CrossRef]
  47. Antoniadis, A.; Tsoutsanis, P.; Rana, Z.; Kokkinakis, I.; Drikakis, D. Azure: An advanced CFD software suite based on high-resolution and high-order methods. 53rd AIAA Aerospace Sciences Meeting 2015. [Google Scholar] [CrossRef]
  48. Fernandez-Fidalgo, J.; Ramirez, L.; Tsoutsanis, P.; Colominas, I.; Nogueira, X. A reduced-dissipation WENO scheme with automatic dissipation adjustment. Journal of Computational Physics 2021, 425. [Google Scholar] [CrossRef]
  49. Tsoutsanis, P.; Kokkinakis, I.; Könözsy, L.; Drikakis, D.; Williams, R.; Youngs, D. Comparison of structured- and unstructured-grid, compressible and incompressible methods using the vortex pairing problem. Computer Methods in Applied Mechanics and Engineering 2015, 293, 207–231. [Google Scholar] [CrossRef]
  50. Srinivasan, H.; Tsoutsanis, P. Adaptive mesh refinement techniques for high-order finite-volume weno schemes. ECCOMAS Congress 2016 - Proceedings of the 7th European Congress on Computational Methods in Applied Sciences and Engineering 2016, 2, 2883–2899. [Google Scholar] [CrossRef]
  51. Tsoutsanis, P.; Drikakis, D. Addressing the challenges of implementation of high-order finite-volume schemes for atmospheric dynamics on unstructured meshes. ECCOMAS Congress 2016 - Proceedings of the 7th European Congress on Computational Methods in Applied Sciences and Engineering 2016, 1, 684–708. [Google Scholar] [CrossRef]
  52. Simmonds, N.; Tsoutsanis, P.; Gaylard, A. Implementation OFA low-Mach number modification for high-order finite-volume schemes for arbitrary hybrid unstructured meshes. ECCOMAS Congress 2016 - Proceedings of the 7th European Congress on Computational Methods in Applied Sciences and Engineering 2016, 2, 2900–2915. [Google Scholar] [CrossRef]
  53. Silva, P.; Tsoutsanis, P.; Antoniadis, A. Simple multiple reference frame for high-order solution of hovering rotors with and without ground effect. Aerospace Science and Technology 2021, 111. [Google Scholar] [CrossRef]
  54. A. S.F. Silva, P.; Tsoutsanis, P.; Antoniadis, A. Numerical investigation of full helicopter with and without the ground effect. Aerospace Science and Technology 2022, 122. [Google Scholar] [CrossRef]
  55. Ricci, F.; Silva, P.; Tsoutsanis, P.; Antoniadis, A. Hovering rotor solutions by high-order methods on unstructured grids. Aerospace Science and Technology 2020, 97. [Google Scholar] [CrossRef]
  56. Adebayo, E.; Tsoutsanis, P.; Jenkins, K. Implementation of CWENO schemes for compressible multicomponent/multiphase flow using interface capturing models. World Congress in Computational Mechanics and ECCOMAS Congress 2022. [Google Scholar] [CrossRef]
  57. Maltsev, V.; Tsoutsanis, P.; Skote, M. HHybrid finite-volume/discontinuous Galerkin framework for the solution of multiphysics problems using unstructured meshes. World Congress in Computational Mechanics and ECCOMAS Congress 2022. [Google Scholar] [CrossRef]
  58. Yuan, D.; Tsoutsanis, P.; Jenkins, K. Hybrid high-order finite volume discontinuous Galerkin methods for turbulent flows. World Congress in Computational Mechanics and ECCOMAS Congress 2022. [Google Scholar] [CrossRef]
  59. Tsoutsanis, P.; Adebayo, E.; Merino, A.; Arjona, A.; Skote, M. CWENO Finite-Volume Interface Capturing Schemes for Multicomponent Flows Using Unstructured Meshes. Journal of Scientific Computing 2021, 89. [Google Scholar] [CrossRef]
  60. Tsoutsanis, P.; Farmakis, P. A posteriori mood limiting approach for multicomponent flows on unstructured meshes. World Congress in Computational Mechanics and ECCOMAS Congress 2022. [Google Scholar] [CrossRef]
  61. Tsoutsanis, P. Stencil selection algorithms for WENO schemes on unstructured meshes. Journal of Computational Physics 2023, 475. [Google Scholar] [CrossRef]
  62. Toro, E.; Spruce, M.; Speares, W. Restoration of the contact surface in the HLL-Riemann solver. Shock Waves 1994, 4, 25–34. [Google Scholar] [CrossRef]
  63. Tsoutsanis, P.; Antoniadis, A.; Jenkins, K. Improvement of the computational performance of a parallel unstructured WENO finite volume CFD code for Implicit Large Eddy Simulation. Computers and Fluids 2018, 173, 157–170. [Google Scholar] [CrossRef]
  64. Petrov, N.V.; Schmidt, A.A. Multiphase phenomena in underwater explosion. Experimental Thermal and Fluid Science 2015, 60, 367–373. [Google Scholar] [CrossRef]
  65. Haimovich, O.; Frankel, S.H. Numerical simulations of compressible multicomponent and multiphase flow using a high-order targeted ENO (TENO) finite-volume method. Computers and Fluids 2017, 146, 105–116. [Google Scholar] [CrossRef]
  66. Yeom, G.S.; Chang, K.S. A modified HLLC-type Riemann solver for the compressible six-equation two-fluid model. Computers and Fluids 2013, 76, 86–104. [Google Scholar] [CrossRef]
  67. Shukla, R.K.; Pantano, C.; Freund, J.B. An interface capturing method for the simulation of multi-phase compressible flows. Journal of Computational Physics 2010, 229, 7411–7439. [Google Scholar] [CrossRef]
  68. Shukla, R.K. Nonlinear preconditioning for efficient and accurate interface capturing in simulation of multicomponent compressible flows. Journal of Computational Physics 2014, 276, 508–540. [Google Scholar] [CrossRef]
  69. Zhang, Z.F.; Wang, C.; Zhang, A.M.; Silberschmidt, V.V.; Wang, L.K. SPH-BEM simulation of underwater explosion and bubble dynamics near rigid wall. Science China Technological Sciences 2019, 62, 1082–1093. [Google Scholar] [CrossRef]
  70. Kleine, H.; Tepper, S.; Takehara, K.; Etoh, T.; Hiraki, K. Cavitation induced by low-speed underwater impact. Shockwaves 2009, 2, 895–900. [Google Scholar] [CrossRef]
  71. Kokh, S.; Lagoutière, F. An anti-diffusive numerical scheme for the simulation of interfaces between compressible fluids by means of a five-equation model. Journal of Computational Physics 2010, 229, 2773–2809. [Google Scholar] [CrossRef]
  72. Oomar, M.; Malan, A.; Jones, B.; Horwitz, R.; Langdon, G. An all-Mach number HLLC based scheme for Multi-phase Flow with Surface Tension. Applied Sciences 2021, 11. [Google Scholar] [CrossRef]
  73. Deng, X.; Inaba, S.; Xie, B.; Shyue, K.M.; Xiao, F. High fidelity discontinuity-resolving reconstruction for compressible multiphase flows with moving interfaces. Journal of Computational Physics 2018, 371, 945–966. [Google Scholar] [CrossRef]
  74. Shyue, K.M. A wave-propagation based volume tracking method for compressible multicomponent flow in two space dimensions. Journal of Computational Physics 2006, 215, 219–244. [Google Scholar] [CrossRef]
  75. Xie, W.F.; Young, Y.L.; Liu, T.G.; Khoo, B.C. Dynamic response of deformable structures subjected to shock load and cavitation reload. Computational Mechanics 2007, 40, 667–681. [Google Scholar] [CrossRef]
  76. Qin, Y.; Wang, Y.; Wang, Z.; Yao, X. The Influence of Various Structure Surface Boundary Conditions on Pressure Characteristics of Underwater Explosion. Computer Modeling in Engineering Sciences 2021, 126, 1093–1123. [Google Scholar] [CrossRef]
  77. Xie, W.; Liu, T.; Khoo, B. Application of a one-fluid model for large scale homogeneous unsteady cavitation: The modified Schmidt model. Computers & Fluids 2006, 35, 1177–1192. [Google Scholar] [CrossRef]
  78. Zhang, Z.F.; Wang, C.; Wang, L.K.; Zhang, A.M.; Silberschmidt, V.V. Underwater explosion of cylindrical charge near plates: Analysis of pressure characteristics and cavitation effects. International Journal of Impact Engineering 2018, 121, 91–105. [Google Scholar] [CrossRef]
  79. Hong, Y.; Wang, B.; Liu, H. Numerical simulation of compressible multiphase hydrodynamic problems using reduced five-equation model on body-fitted grids. Computers and Fluids 2022, 233, 105243. [Google Scholar] [CrossRef]
  80. Jafarian, A.; Pishevar, A. An exact multiphase Riemann solver for compressible cavitating flows. International Journal of Multiphase Flow 2017, 88, 152–166. [Google Scholar] [CrossRef]
  81. CHENG Juan, S.C.W. High Order Schemes for CFD: A Review. Chinese Journal of Computational Physics 2009, 26, 633–655. Available online: http://www.cjcp.org.cn/EN/abstract/article_825.shtml.
  82. Wang, Y.; Yao, X.; Qin, Y. Investigation on influence factors about damage characteristics of ice sheet subjected to explosion loads: Underwater explosion and air contact explosion. Ocean Engineering 2022, 260, 111828. [Google Scholar] [CrossRef]
  83. Wardlaw, A.B.; Luton, J.A. Fluid-structure interaction mechanisms for close-in explosions. Shock and Vibration 2000, 7, 265–275. [Google Scholar] [CrossRef]
  84. Cheng, L.; Deng, X.; Xie, B.; Jiang, Y.; Xiao, F. Low-dissipation BVD schemes for single and multi-phase compressible flows on unstructured grids. Journal of Computational Physics 2021, 428, 110088. [Google Scholar] [CrossRef]
  85. Schmidt, D.; Rutland, C.; Corradini, M. A fully compressible, two-dimensional model of small, high-speed, cavitating nozzles. Atomization and Sprays 1999, 9, 255–276. [Google Scholar] [CrossRef]
  86. Yu, J.; Liu, J.h.; Wang, H.k.; Wang, J.; Zhou, Z.t.; Mao, H.b. Application of two-phase transition model in underwater explosion cavitation based on compressible multiphase flows. AIP Advances 2022, 12. [Google Scholar] [CrossRef]
  87. Yingtang, D.; Lanhao, Z.; Jia, M.; Eldad, A. A resolved CFD-DEM-IBM algorithm for water entry problems. Ocean Engineering 2021, 240, 110014. [Google Scholar] [CrossRef]
Preprints 94023 g001
Figure 1. Computational domain setup for a two-dimensional underwater explosion near a free surface.
Figure 1. Computational domain setup for a two-dimensional underwater explosion near a free surface.
Preprints 94023 g001
Preprints 94023 g002
Figure 2. Underwater explosion near a free surface; 2D unstructured finest mesh with triangular elements used for this test case. The density of mesh elements in specific areas of the domain is increased to better capture the physical phenomena occurring in those regions. The zoomed-in area on the right provides a more detailed view of the mesh’s structure.
Figure 2. Underwater explosion near a free surface; 2D unstructured finest mesh with triangular elements used for this test case. The density of mesh elements in specific areas of the domain is increased to better capture the physical phenomena occurring in those regions. The zoomed-in area on the right provides a more detailed view of the mesh’s structure.
Preprints 94023 g002
Preprints 94023 g003
Figure 3. Underwater explosion near a free surface; Schileren visualization of density gradient using a CWENO4 scheme on a resolution of 1600 × 1600 cells (finest mesh) at t = ( 0.2 , 0.4 , 0.8 , 1.2 ) ms.
Figure 3. Underwater explosion near a free surface; Schileren visualization of density gradient using a CWENO4 scheme on a resolution of 1600 × 1600 cells (finest mesh) at t = ( 0.2 , 0.4 , 0.8 , 1.2 ) ms.
Preprints 94023 g003
Preprints 94023 g004
Figure 4. UNDEX near a free surface; Schileren visualization of pressure gradient using a CWENO4 scheme on the finest mesh at t = ( 0.2 , 0.4 , 0.8 , 1.2 ) ms.
Figure 4. UNDEX near a free surface; Schileren visualization of pressure gradient using a CWENO4 scheme on the finest mesh at t = ( 0.2 , 0.4 , 0.8 , 1.2 ) ms.
Preprints 94023 g004
Preprints 94023 g005
Figure 5. UNDEX near a free surface; contours of volume fraction for underwater explosion near a free surface using a CWENO4 scheme on finest grid at t = ( 0.2 , 0.4 , 0.8 , 1.2 ) ms.
Figure 5. UNDEX near a free surface; contours of volume fraction for underwater explosion near a free surface using a CWENO4 scheme on finest grid at t = ( 0.2 , 0.4 , 0.8 , 1.2 ) ms.
Preprints 94023 g005
Preprints 94023 g006
Figure 6. UNDEX near a free surface; density gradient contours at t = ( 0.2 , 0.4 , 0.8 , 1.2 ) ms. Captured using CWENO4 schemes with different grids resolution. More inter-facial structures are sharper as the grid resolution increases.
Figure 6. UNDEX near a free surface; density gradient contours at t = ( 0.2 , 0.4 , 0.8 , 1.2 ) ms. Captured using CWENO4 schemes with different grids resolution. More inter-facial structures are sharper as the grid resolution increases.
Preprints 94023 g006
Preprints 94023 g011
Figure 11. Computational domain setup for a two-dimensional underwater explosion near a planar rigid wall.
Figure 11. Computational domain setup for a two-dimensional underwater explosion near a planar rigid wall.
Preprints 94023 g011
Preprints 94023 g012
Figure 12. UNDEX near a rigid wall; the pressure contours at different instants. Captured using CWENO3 scheme with medium grid resolution. At very low pressure, cavitation can be noticed at time t = 4 ms and 5.3 ms near the rigid walls
Figure 12. UNDEX near a rigid wall; the pressure contours at different instants. Captured using CWENO3 scheme with medium grid resolution. At very low pressure, cavitation can be noticed at time t = 4 ms and 5.3 ms near the rigid walls
Preprints 94023 g012
Preprints 94023 g013
Figure 13. UNDEX near a rigid wall; plots of pressure history for the underwater explosion near a planar rigid wall at the center location of the upper wall obtained with CWENO4 schemes, and compared with other published results [65,77].
Figure 13. UNDEX near a rigid wall; plots of pressure history for the underwater explosion near a planar rigid wall at the center location of the upper wall obtained with CWENO4 schemes, and compared with other published results [65,77].
Preprints 94023 g013
Preprints 94023 g014
Figure 14. Computational domain setup for three-dimensional underwater explosion in a rigid cylinder.
Figure 14. Computational domain setup for three-dimensional underwater explosion in a rigid cylinder.
Preprints 94023 g014
Preprints 94023 g015
Figure 15. Underwater explosion in an enclosed cylindrical container; plots of pressure history (at the centre location of the left wall ) for the underwater explosion in an enclosed cylindrical wall with CWENO4 schemes, and compared with other published results. [9,81,82]
Figure 15. Underwater explosion in an enclosed cylindrical container; plots of pressure history (at the centre location of the left wall ) for the underwater explosion in an enclosed cylindrical wall with CWENO4 schemes, and compared with other published results. [9,81,82]
Preprints 94023 g015
Preprints 94023 g016
Figure 16. a. Mesh used for three-dimensional cylindrical underwater explosions. b. Corresponding mesh refinement for simulation zone.
Figure 16. a. Mesh used for three-dimensional cylindrical underwater explosions. b. Corresponding mesh refinement for simulation zone.
Preprints 94023 g016
Preprints 94023 g017
Figure 17. UNDEX in an enclosed cylindrical container; pressure gradient contours for the 3D underwater explosion in an enclosed, rigid cylindrical wall.
Figure 17. UNDEX in an enclosed cylindrical container; pressure gradient contours for the 3D underwater explosion in an enclosed, rigid cylindrical wall.
Preprints 94023 g017
Preprints 94023 g020
Figure 20. UNDEX of single bubble (top row) and two bubbles configurations placed horizontally near a free surface; contour of volume fraction using a CWENO4 scheme on a medium mesh at t=(0.2, 0.4, 0.8 and 1.2 ) ms (from left to right) and for inter-bubble lengths of (0.3, 0.5, 0.7 and 0.9) m from top to bottom respectively.
Figure 20. UNDEX of single bubble (top row) and two bubbles configurations placed horizontally near a free surface; contour of volume fraction using a CWENO4 scheme on a medium mesh at t=(0.2, 0.4, 0.8 and 1.2 ) ms (from left to right) and for inter-bubble lengths of (0.3, 0.5, 0.7 and 0.9) m from top to bottom respectively.
Preprints 94023 g020
Preprints 94023 g021
Figure 21. UNDEX of a single and two bubbles configurations placed horizontally near a free surface; contour of pressure gradient for the single bubble (Top) and the two-bubbles cases subsequently below for inter-bubble lengths of (0.3, 0.5, 0.7 and 0.9) m respectively using a CWENO4 scheme on a medium mesh at t=(0.4, 0.8 and 1.06 )ms. The red-yellow areas are cavitation zones.
Figure 21. UNDEX of a single and two bubbles configurations placed horizontally near a free surface; contour of pressure gradient for the single bubble (Top) and the two-bubbles cases subsequently below for inter-bubble lengths of (0.3, 0.5, 0.7 and 0.9) m respectively using a CWENO4 scheme on a medium mesh at t=(0.4, 0.8 and 1.06 )ms. The red-yellow areas are cavitation zones.
Preprints 94023 g021
Preprints 94023 g022
Figure 22. Computational domain setup for the three-dimensional underwater of two bubbles placed vertically in an enclosed region
Figure 22. Computational domain setup for the three-dimensional underwater of two bubbles placed vertically in an enclosed region
Preprints 94023 g022
Preprints 94023 g025
Figure 25. Underwater explosion of two bubbles vertically placed in an enclosed region; pressure gradient contours for the 3D underwater explosion in an enclosed rigid cylindrical wall with three different inter -bubble distances: H b = 15 mm(Top), 30 mm(Middle), 60 mm(Bottom) for varying time frames t = (1.54, 6, 30, 60 and 90) µs. Captured using CWENO4. The white-red-yellow areas are cavitation zones, the red regions are indicative of potential high cavitation activity.
Figure 25. Underwater explosion of two bubbles vertically placed in an enclosed region; pressure gradient contours for the 3D underwater explosion in an enclosed rigid cylindrical wall with three different inter -bubble distances: H b = 15 mm(Top), 30 mm(Middle), 60 mm(Bottom) for varying time frames t = (1.54, 6, 30, 60 and 90) µs. Captured using CWENO4. The white-red-yellow areas are cavitation zones, the red regions are indicative of potential high cavitation activity.
Preprints 94023 g025
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

facebook logotwitter logolinkedin logo
MDPI Initiatives
Important Links
Subscribe

Disclaimer

Terms of Use

Privacy Policy

© 2024 MDPI (Basel, Switzerland) unless otherwise stated