Investigation of cloud cavitating flow in a venturi using Adaptive Mesh Refinement (AMR)
Dhruv Apte Mingming Ge Olivier Coutier-Delgosha


Investigation of cloud cavitating flow in a venturi using Adaptive Mesh Refinement (AMR)]Investigation of cloud cavitating flow in a venturi using Adaptive Mesh Refinement (AMR)










[1]Dhruv Aptedhruvga@vt.edu

2]Mingming Gemmge@vt.edu


1,3]Olivier Coutier-Delgoshaocoutier@vt.edu


*[1]Kevin T. Crofton Department of Aerospace and Ocean Engineering, Virginia Tech,
Blacksburg, 24060, VA, USA

[2]Macao Environmental Research Institute, Faculty of Innovation Engineering, Macau University of Science and Technology, Macao, 999078, Macao SAR, China

[3]Univ. Lille, CNRS, ONERA, Arts et Metiers ParisTech, Centrale Lille,  Lille, F-59000, France







Unsteady cloud cavitating flow is detrimental to the efficiency of hydraulic machinery like pumps and propellers due to the resulting side-effects of vibration, noise and erosion damage. Modelling such a unsteady and highly turbulent flow remains a challenging issue. In this paper, cloud cavitating flow in a venturi is calculated using the Detached Eddy Simulation (DES) model combined with the Merkle model. The Adaptive Mesh Refinement (AMR) method is employed to speed up the calculation and investigate the mechanisms for vortex development in the venturi. The results indicate the velocity gradients and the generalized fluid element strongly influence the formation of vortices throughout a cavitation cycle. In addition, the cavitation-turbulence coupling is investigated on the local scale by comparing with high-fidelity experimental data and using profile stations. While the AMR calculation is able to predict well the time-averaged velocities and turbulence-related aspects near the throat, it displays discrepancies further downstream owing to a coarser grid refinement downstream and under-performs compared to a traditional grid simulation . Additionally, the AMR calculations is unable to reproduce the cavity width as observed in the experiments. Therefore, while AMR promises to speed the process significantly by refining grid only in regions of interest, it is comparatively in line with a traditional calculation for cavitating flows. Thus, this study intends to provide a reference to employing AMR as a tool to speed up calculations and be able to simulate turbulence-cavitation interactions accurately.





















[
    [
    
=====





§ INTRODUCTION


Cavitation is a multiphase, highly unstable and turbulent phenomenon characterized by formation of clouds of bubbles when the ambient pressure drops below the vapour pressure. These clouds of bubbles travel with high velocity and upon exiting the low-pressure area burst, generating shock. This shock can cause erosion damage, vibration, load asymmetry and other detrimental effects that can drastically impact the performance of marine engineering devices like propellers, pumps and turbines. Additionally, the capability to cause erosion can also be utilized in jet drilling in the hydrocarbon and geothermal energy sector <cit.>. Therefore, it is important to investigate the phenomenon, especially the mechanism of periodic development of the main cavity followed by its shedding and collapse.

A primary component of cavitation is its interaction with turbulence and the resulting generated vortices. Gopalan & Katz <cit.> used PIV and high-speed photography to demonstrate that vortice generation was a result of the collapse of the vapor cavities in the cavity closure region. More recently, Arabnejad et al. <cit.> investigated the cavitation over the leading edge of a NACA0009 hydrofoil using High-Speed Visualization (HSV) to observe the transformation of the shedding of cavity structures into horse-show vortices. Ge et al. <cit.> conducted similar experiments in a venturi-type cavitation reactor and observed a similar  horseshoe structure. They further concluded that the horseshoe broke up into a long, thin chain of bubbles with the upper parts merging with the clouds downstream. 

Concurrently, studies using numerical simulations to model cavitating flows have been gaining more and more interest. A successful modelling of cavitating flows typically necessitates the combination of a cavitation model and a turbulence model. Regarding turbulence models, Direct Numerical Simulation (DNS) where all turbulence scales are resolved seems an evident approach but its demand for an extremely fine mesh and computational power hinders its application in practice. Large Eddy Simulation (LES) studies seem as a viable option and several studies have been conducted to investigate the cavity and vortex dynamics driving cavitating flows <cit.>. In fact, LES has also been coupled with multi-scale Euler-Lagrangian approaches <cit.> to investigate both the large vapor volumes from an Eulerian point-of-view and the micro-scale bubble dynamics using a Lagrangian reference. It was observed that the collision of the re-entrant jet with the cavity interface results in the formation of large number of bubbles. Tian et al. <cit.> further confirmed these findings with their multi-scale approach and concluded that the generation of the large number of bubbles is a result of Kelvin-Helmholtz instability, triggered by the shear flow between water and vapor. Ji et al. <cit.> have conducted a comprehensive review of multi-scale techniques for simulating cavitating flows. However, LES also requires high-resolution grids which results in higher computational costs <cit.>. Thus, Reynolds-averaged Navier-Stokes (RANS) models, where all turbulence scales are modelled seems as the practical option to model cavitating flows. 

Coutier-Delgosha et al. <cit.> conducted RANS calculations using the k-ϵ Re-Normalization Group (RNG) and k-ω models with and without an empirical correction suggested by Reboud et al. <cit.> to simulate cavitating flows. They concluded that the standard models over-predict the eddy-viscosity and thus are unable to predict the periodic shedding of the cavity as observed in experiments. The same model was also employed by Ji et al. <cit.> to simulate the unsteady cavitating flow around a twisted hydrofoil. They noted that the cavitating flow induced boundary layer separation and reported large increases of vorticity at the cavity interface. Long et al. <cit.> advanced this approach by combining the Reboud correction with the Filter-Based Method (FBM) <cit.>, a recent type of model, termed as hybrid RANS-LES models that behave as RANS model near the wall and as LES model away from the wall, therefore combining the accuracy of LES models and computational efficiency of RANS models. They conducted a vorticity budget analysis, based on the Vorticity Transport Equation  (VTE) in attached cavitation around a Clark-Y hydrofoil. They concluded that the baroclinic torque term in the VTE is the principal source of vorticity production during the collapse of the cavitation cloud. Sun et al. <cit.> conducted a similar calculation using Partially-Averaged Navier-Stokes model (PANS) proposed by Girimaji et al. <cit.> to simulate unsteady cavitation around a NACA0015 hydrofoil in thermo-sensitive fluid. They observed that the vortex structures and motion show considerable unsteady properties in the cavity shedding region at the trailing edge. Additionally, Apte et al. <cit.> conducted a systematic analysis of several RANS and hybrid RANS-LES models to simulate cloud cavitating flow inside a venturi nozzle and concluded that while these models are able to simulate the unsteady cavity shedding, they are unable to simulate the cavity dynamics on the local scale. They suggested grid refinement as a possible solution to ensure the LES zone occurs on the fine grid. 

To solve the problem of the resulting high requirement in computational resources, studies have proposed mesh refinement techniques that promote refinement in specific flow regions rather than refining the entire domain. Bai et al. <cit.> devised a Cartesian cut-cell method where the rotation factor, the ratio of the strain rate to the sum of the strain rate and vorticity tensor was applied as a marker for cell refinement while modelling a tip leakage cavitating flow generated by a hydrofoil. The methodology was successfully extended to Cheng et al. <cit.> who conducted a systematic analysis of various LES models to simulate tip-vortex cavitation. They concluded the method was able to simulate tip vortex cavitation accurately and were able to capture the influence of the vorticity dilatation term in decreasing the vorticity inside the tip-leakage vortex section.  Another technique, titled Adaptive Mesh Refinement (AMR) was proposed <cit.> that allocates high resolution grid based on the input property and can be a promising technique to understand better the vortex dynamics in cavitating flows. AMR has been used for various multiphase flow applications like primary atomization <cit.>, bubble dynamics <cit.> and cavitation as well. Li et al. <cit.> used AMR to aid LES for simulating cavitation over a cylindrical body. Wang et al. <cit.>, employed a similar approach to model cavitating flow around a Clark-Y hydrofoil. However, both these papers have leveraged LES, which continues to be prohibitively expensive in some case. In addition, while these studies delve into the aspects of cavitation generation mechanisms, there remains a gap to investigate the ability of AMR to capture the cavitation-turbulence interplay.  

This study aims to address this gap by simulating cloud cavitation using AMR coupled with the Detached Eddy Simulation (DES) model, a hybrid RANS-LES model. Thus, the study attempts to strike a delicate balance between computational time and accuracy. 



§ MATHEMATICAL FORMULATION AND NUMERICAL METHOD



 §.§ Adaptive Mesh Refinement

To overcome the issue of using highly-refined grids, the AMR is used to focus the refinement dynamically in the region most needed. Mathematically, AMR can be formulated as refining a domain comprising of non-overlapping rectangular grid G_l,k:

    G_l = ∪_m G_l,k

where l=0, 1, 2, ... l_max represents refinement levels. In OpenFOAM <cit.>, AMR is based upon three key drivers, the refinement engine, mesh cutter and the refinement tree. The refinement engine is responsible for ensuring the connectivity structures in the mesh while the mesh cutter refines the mesh as per the input. Finally, the refinement tree saves the refinement history and ensures a steady transition for mesh refinement. Formulated by Mitchell et al. <cit.>, the refinement tree is defined as:

    T (G) = {Ω, {C(ν_i)}}

which contains a set of nodes Ω = ν_i_i=0^M with each ν_i∈Ω containing a set of child C(ν_i) ⊂Ω where G represents the grid. Fig <ref> presents the workflow of the three key portions in AMR. The solver reads the refinement criterion as the solver reaches the refinement interval . This refinement criterion could be based on a numerical domain for turbulent kinetic energy or pressure field or as in this case, the void fraction. Next, the solver proceeds to the refinement engine which selects the appropriate zones for grid refinement. The mesh cutter then refines the specific areas while the refinement tree saves the entire refinement history. The refinement occurs in the form of an octree structure which means, on every single refinement level, the number of cells increases by a factor of eight. After refinement, a similar unrefinement process is conducted for cells that no longer fit the refinement criteria, with the solution being mapped to the new cells. The refinement criterion is based on the void fraction to ensure the grid refinement happens in the interface between the cavitating flow and single-phase, turbulent flow.  




 §.§ Numerical models



  §.§.§ Basic Governing Equations

As previously mentioned,modelling cavitating flows involves a coupling of a cavitation model and a turbulence model. The work here uses the Transport-Equation Model approach (TEM) where the two phases are considered to be strongly coupled and governed by the same momentum and mass transfer equations:

    ∂ (ρ_m u_i)/∂ t+∂ (ρ_m u_i u_j)/∂ x_j
    
    =-∂ p/∂ x_i+∂/∂ x_j((μ_m +μ_t)(∂ u_i/∂ x_j+∂ u_j/∂ x_i- 2/3∂ u_k/∂ x_kδ_ij))




    ∂ρ_lα_l/∂ t+∂(ρ_lα_l u_j)/∂ x_j= ṁ^- + ṁ^+


    ρ_m= ρ_lα_l+ρ_vα_v


    μ_m= μ_lα_l+μ_vα_v

where u_j is the velocity component in the jth direction, ρ_m and μ_m are respectively the density and viscosity of the mixture phase, u is the velocity, p is the pressure, ρ_l and ρ_v are respectively the liquid and vapor density, μ_l and μ_v are respectively the liquid and vapor dynamic viscosity while μ_t represents the turbulent viscosity. α_l and α_v are respectively the liquid and vapor void fraction. The source (ṁ^+) and sink (ṁ^-) terms represent the condensation (vapor destruction) and evaporation (vapor formation) terms respectively and will be discussed below.

The separation of condensation and evaporation terms in the above equation is the principal concept of the TEM methodology. The separation aids in treating the two processes separately if needed. The TEM model used in the study is the Merkle model <cit.> which defines the evaporation and condensation terms as:

    ṁ^-= C_dest min(p-p_sat,0)γρ_l/0.5 U_∞^2 t_∞ρ_v


    ṁ^+= C_prod max(p-p_sat,0)(1-γ)/0.5 U_∞^2 t_∞

where γ is the liquid volume fraction, ρ_v and ρ_l are the vapor density and the liquid density respectively, p and p_sat are the pressure and the saturation pressure respectively, t_∞ is the free stream time scale and U_∞ is the free stream velocity.  The empirical factors C_dest and C_prod are set as 80 and 1e-3 respectively, according to the study by Merkle et al. <cit.>.



  §.§.§ Turbulence model


The turbulence model used is the Detached Eddy Simulation (DES) model with the k-ω SST model as the baseline RANS model <cit.>. Here, the LES model calculates the turbulent separated zone while the RANS model is used to model the near-wall zone. The switch to from RANS to LES and vice versa is defined by the change in length scale:

    L_DES= min (L_t, C_DESΔ)

where L_T is the turbulent length scale defined as:

    L_t = √(k)/β * ω

and C_DESΔ is the DES filter length with C_DES = 0.61



§ COMPUTATIONAL SETUP





A converging-diverging (venturi) nozzle is selected as the geometry utilized for the study as shown in <ref>. The venturi nozzle has been widely used in jetting for geothermal reservoirs and hydrocarbon industry <cit.>. The geometrical setup here is identical to the experiment <cit.>.    The nozzle height at the inlet and outlet is originally 21 mm but steadily drops to 10 mm at the throat. The venturi has a converging angle of 18 degrees and divergent angle of 8 degrees (see Fig <ref>). The velocity inlet and pressure outlet boundary conditions are adopted. For consistency and validation with the experiments, the velocity at inlet is set to 8.38 m/s while the pressure at the outlet is adjusted to ensure both experiments and simulations have the same mean cavity length. The Re number is 1.8 × 10^5. To evaluate whether AMR is able to capture the cavity dynamics accurately at the local stage, a similar standard mesh calculation is conducted. The mesh consists of 8.4 million cells, since the cell limit of the AMR mesh is 8.4 million cells as well <cit.>. Table <ref> shows the mean cavity lengths of the two cases, as compared to the experimental case. All three cases have the same mean cavity length, thus validating the simulations and providing an opportunity to investigate the cavitation-vortex interaction and the cavitation-turbulence interplay.



The calculations are conducted using interPhaseChangeFoam and interPhaseChangeDyMFoam.
interPhaseChangeFoam is an unsteady, isothermal solver that couples a cavitation model with a turbulence model while the latter is its sister-solver that uses adaptive mesh refinement. It uses a dynamicMeshDict for the dynamic mesh properties as input. The process of refining the grid is dependant on the field of void fraction, ranging from 0.01 to 0.9. This ensures the grids near the vapor cavity-water interface are refined. For this study, the mesh properties are based on the void fraction field with a maximum cell limit of 8.4 million cells from a base size of 3.3 million cells (see Fig <ref>). Both solvers utilize the PIMPLE algorithm, a hybrid algorithm comprising of both the PISO and SIMPLE algorithms used in OpenFOAM. The PIMPLE algorithm consists of three parts:

    
  * momentum predictor
    
  * pressure solver
    
  * momentum corrector

In PIMPLE, each timestep has at least one outer corrector loop where the fields are solved using a number of iterations. The algorithm exhibits a more robust pressure-velocity coupling by coupling a SIMPLE outer corrector loop with a PISO inner corrector loop and has been shown to display better numerical stability. 
Regarding numerical schemes, the volume fraction interpolation uses the Gauss vanLeer scheme. To mitigate the unboundedness, a semi-implicit  multi-dimensional limiter for explicit solution (MULES) is implemented: an implicit corrector step is first implemented corresponding to the discretization schemes and then an explicit correction is applied with the MULES limiter. The calculations are run at a fixed timestep of 1e-5. However, to ensure stability of the calculation, the AMR mesh is run at a timestep subjected to the maximum Courant number, 1 in this study. The refineInterval utility is set to 10 thus refining the mesh at 10 timesteps. The refinement is based on the alpha.water field or the void fraction field and set between the limits of 0.01 and 0.9. This ensures the mesh would be refined on lower void fraction values, where the cavity is expected to develop.



The simulation is set in three consecutive stages. First, it is run for 0.03 s where the vaporization coefficients in the cavitation model are set to zero thereby modelling it as a single-phase, turbulent flow. The next stage, consisting of 0.03s is a sinusoidal ramp whereupon at the end, fully cavitation has been fully launched. At the third and final stage, the fully cavitating regime is activated.  The focus of the work is on the fully cavitating regime with AMR also initiated at this stage. To maintain lower computational costs, the full cavitation regime is launched for 0.02 s, sufficient enough to capture a complete shedding cycle.



§ RESULTS



 §.§ Cavitation flow structures



Fig <ref> shows the pressure and void fraction distributions at five successive time instants in the diverging section of the venturi. In the first column, the pressure distributions are shown, with the dark blue region denoting the low pressure regions. The second column represents the same snapshots with void fraction. Here, the blue region denotes vapor while the red region denotes water.  The distributions show the evolution of cavitation starting from the first one, where a cavity is initiated at the throat of the venturi. Near the throat, a sharp drop in pressure is observed, responsible for the inception of the cloud cavity. Downstream, a group of smaller cloud cavities is observed. These cavities are remnants of the previous cycle. The next figures Fig <ref> (c)-(d) demonstrate the growth of the cavity at the throat while the cavities downstream flow downstream and collapse as it leaves the low-pressure region. The next set of figures (e)-(f) show the cavity reaching its maximum size. In (g)-(h), it is observed that the large cavity at the throat is pinched-off due to a re-entrant jet rushing upstream, resulting in a secondary cloud cavity along with smaller vapor cavity clouds. The primary cavity continues to reduce and collapses while the secondary cavity continues to roll up and flow downstream. In (i)-(j), it is observed another small incipient cavity develops at the throat while the detached secondary cavity of this cycle continues to be shed downstream. Thus, the snapshots demonstrate the periodic nature of cloud cavitating flows and the presence of a re-entrant jet responsible for this periodic shedding of the primary cavity cloud.



 §.§ Adaptive Mesh Refinement inside the nozzle


As stated previously, the AMR method refines the grid only in the concerning region based on the input refinement criterion. To check the mesh refinement inside the venturi nozzle for unsteady cavitating flow, six snapshots are taken at the mid-plane of the venturi and focused on the diverging section of the nozzle. The grid is then visualized with the void fraction field. Fig <ref> represents the grid during the cavity initiation stage. A small cavity at the throat is observed followed by a highly refined grid , in sharp contrast to the structured grid at the top of the venturi nozzle. The grid refinement extends further beyond the initial cavity, possibly as a result of other smaller cavities throughout the region. Further downstream, the refinement exists in selected regions, where the remnants of the previous cavity shedding cycle are flowing downstream. As the incipient cavity enlarges, shown in Fig <ref>, the grid refinement closer to the throat is observed. However, further downstream where the detached cavities collapse after exiting the low-pressure region, the grid reverts back to the structured grid defined at the onset of the calculation. This hypothesis is substantiated in Fig <ref> where the cavity reaches its maximum length. The grid in this region is significantly refined, while the grid downstream reverts completely to the structured grid devised as the original setup. Thus, it can be noted that the AMR method works well in refining solely the low vapor regions and giving more insights about the flow physics driving the unsteady cavitating flow at a much reduced computational cost.      







Fig <ref> shows the cavity detachment process. The presence of both a reduced primary cavity at the throat and a secondary cavity downstream results in a high grid refinement throughout the lower portion of the venturi nozzle, with some refinement zones slightly upwards as well, as a result of smaller shed vapor cavities in other planes. The cavity collapse process, depicted in Fig <ref> shows the primary cavity collapsing and the detached cavity flowing downstream. Significant grid refinement continues to be conducted in the section accounted by the detached cavity downstream and the presence of smaller vapor cavities present closer to the throat. The presence of a refined grid in regions containing smaller vapor cavities illustrate that the AMR method is very efficient tool to simulate unsteady cavitating flows by dramatically increasing the grid refinement while striking a balance between grid refinement and computational cost.    



 §.§ Cavitation-vortex interaction


For a better understanding of the cavitating flow and the resulting vortex interaction, the Q-criterion is plotted for various stages in a cavitating cycle. The Q-criterion defined as the second invariant of the velocity gradient tensor is defined as

    Q = 1/2 (|Ω|^2 - |S|^2)

where Ω is the vorticity tensor and S is the rate of strain tensor. 
The Q-criterion and corresponding vapor structures are presented in the following figures. Here, a snapshot of the venturi from the front is taken of contours of the α = 0.6.


Fig <ref> shows the void fraction and Q-criterion during cavitation inception. Here, the cavity is just initiated at the cavity throat while the remnants of the previous cavity appear further downstream including, the entire detached cloud cavity of the previous cycle.  
To further understand the mechanism of cavitation-vortex interaction, the vorticity-transport equation is employed at the mid-plane in the z-direction to investigate the vorticity and the distribution of its contributing terms, shown in the following equation


    Dω_z/Dt = [(ω·∇) V]_z - [ω (∇· V)]_z +  ( ∇ρ_m×∇ p/ρ_m^2)_z + [(ν_m + ν_t) ∇^2ω]_z


Here,the Left Hand Side (LHS) denotes the rate of vorticity change while the Right Hand Side (RHS) indicate the vortex stretching, vortex dilatation, baroclinic torque and viscous diffusion of vorticity terms respectively. The vortex stretching term describes the stretching and tilting of a vortex due to the velocity gradients. The vortex dilatation describes the expansion and contraction of a fluid element. The baroclinic torque is a result of the misalignment between the pressure and velocity gradients. The viscous diffusion term can be ignored in high Reynolds number flows <cit.> and thus is omitted in the study. The equation clearly shows the effects of velocity and pressure gradients formed as a result of cavitating flow on the vorticity. The equation is further simplified as:

    ω_z = ∂ V_y/∂ x - ∂ V_x/∂ y


    [(ω·∇) V]_z = ω_x∂ V_z/∂ x + ω_y∂ V_z/∂ y + ω_z∂ V_z/∂ z


    [ω (∇· V)]_z = ω_z (∂ V_x/∂ x+∂ V_y/∂ y+ ∂ V_z/∂ z)


    ( ∇ρ_m×∇ p/ρ_m^2)_z = 1/ρ_m^2 (∂ρ_m/∂ x·∂ p/∂ y - ∂ρ_m/∂ y·∂ p/∂ x)




Fig <ref> shows the different terms of the vorticity transport equation at the cavity initiation stage. It is observed, that while the vortex stretching and dilatation terms show highly positive values, specially at the throat, thus dominating the vortex formation. The baroclinic torque does not influence the vorticity considerably and it is demonstrated that the velocity gradients and expansion of the fluid element influences the cavitation-vortex interaction. 


Fig <ref> shows the next stage of the cavitation shedding. The previously small incipient cavity grows larger, developing into a cloud shape. Here, the Q-criterion is negative across the primary cavity and the detached cavity of the previous cycle, indicating the vorticity tensor is smaller than the rate of strain tensor. However, it is observed that the vortex stretching and dilatation terms continue to dominate the vorticity tensor compared to the baroclinic torque. 






Fig <ref> shows the Q-criterion and the vorticity equation terms as the cavity reaches its maximum size. Here, it is observed that the secondary cavity of the previous cycle has collapsed upon exiting the low-pressure region. The Q-criterion is distributed messily throughout the cavity with the vortex stretching and vortex dilatation term also displaying similar behaviours: the terms are positive along the cavity interface but negative inside the cavity, showing the very drastic change in vorticity from the cavity interface to the center of the cavity. As the cavity reaches its maximum size, it encounters a re-entrant jet travelling in the direction opposite to the main flow direction that breaks the cavity into two, a reduced primary cavity at the throat and a detached secondary cavity downstream (Fig <ref>). Here, the Q-criterion is considerably high near the cavity interface but dominated by the strain rate tensor inside the cavity. While the vortex stretching term remains positive at the cavity-water interface at the point of cavity detachment, the vortex dilatation term gets extremely at the trailing edge of the primary cavity and leading edge of the secondary, rolled-up cavity. It can be observed, the shedding of the primary cavity leads to generation of vortices in both direction. The figures demonstrate the dramatic change of vortices caused by cloud cavitation and the resulting interplay between strain rate and the vorticity tensors and how the expansion and contraction of the fluid element overall dominates the process.  The shed cavity flows downstream and collapses, upon exiting the low pressure, as seen in Fig <ref>. At the same instant, a small incipient cavity is observed at the throat. Thus, the periodic shedding of cavity clouds is observed. With the collapse of the cloud cavity, the Q-criterion is positive in very few areas, demonstrating the role of strain rate tensor in the process. In the areas where vorticity is observed, near the cavity-water interface of the detached cloud cavity, it is observed the vortex shedding dominates the process. While studies <cit.> have postulated that baroclinic torque might be responsible for vorticity production during cloud cavity collapse, baroclinic torque has not been observed as extremely influential in vorticity generation in these numerical simulations. Throughout the stages of cavitation, it is observed that the velocity gradients and the growth  shrinking of the fluid element resulting due to cavitation influence the vortex formation thus providing insights into the cavitation-vortex interaction.    



 §.§ Cavity dynamics at the local scale

While the cavity dynamics at global scale is discussed, it is important to investigate if turbulence modelling with AMR is able to capture the cavity dynamics at the local scale.    



Fig <ref> shows the venturi nozzle at cut at mid-plane about the z-axis with the profile stations for local analysis. These stations are 1.5mm, 3 mm, 5mm, 10 mm, 15 mm and 20 mm from the throat respectively. The profile stations have been placed at varying distance from the throat to capture the cavity and investigate if the numerical results are able to simulate the turbulence properties as observed in experiments. Figs <ref> and <ref> shows the profiles for time-averaged velocities in both stream-wise and wall directions respectively. The black dots represent the experimental data, taken from Ge et al. <cit.>. Near the throat, the stream-wise velocity near the venturi's diverging wall is less than the flow velocity but jumps considerably between 1 and 2 mm from the wall. Downstream, the jump is much less pronounced, with the velocity increasing steadily. Both models are able to predict the stream-wise time-averaged velocity well, with slight differences in the 10 mm profile. Extending the discussion to the velocity in wall-direction, different results are obtained. The velocity jump away from the wall is less pronounced throughout the profiles as compared to the velocity in the stream-wise direction. The numerical simulations show considerable discrepancies near the throat as they predict a sharper velocity jump slightly closer to the wall than observed in the experiments. This indicates that while the simulations are able to reproduce the same cavity length as measured in the experiments, they are unable to reproduce the same cavity width at the throat. This discrepancy could stem from the cavitation modelling strategy. In addition, subtle differences appear in the downstream profiles, starting from the 10 mm profile. Away from the wall, the AMR mesh predicts higher velocity than the non-AMR mesh and experimental data. At the farthest profile, the 20 mm one, the AMR profile shows some minor pulsations close to the bottom wall of the venturi which are not captured by the non-AMR mesh. These pulsations could be attributed to the smaller cavity clouds formed after the primary cavity detaches due to the re-entrant jet rushing upwards.    







The analysis is extended to turbulence data like the Reynolds shear stress and the Turbulent Kinetic Energy (TKE) as well, in Figs <ref> and <ref>. It is observed at the throat, that the AMR calculation predicts the Reynolds stress closer to the experimental data as compared to the non-AMR mesh which over-predicts the Reynolds shear stress. This co-relation extends until the 5 mm profile. It is posited that since the AMR is a function of low void fraction, the grid near the throat becomes considerably more refined than the standard mesh. As a result, the area is modelled by LES rather than URANS, as is the case in DES models, thus leading to better results as compared to standard meshes. However, downstream predictions show considerable discrepancies as the AMR mesh under-predicts the Reynolds stress while the standard mesh performs comparatively better. Downstream, the absence of low void fraction results in a coarser mesh as compared to the region close to the throat. The coarser mesh results in the modelling of the region by URANS rather than LES, thus under-predicting the turbulence stress as previously observed in URANS simulations. The TKE plots show similar results. As we go downstream, the TKE increases continuously and farther away from the bottom wall. The AMR mesh under-predicts the TKE at the profiles near the throat but are able to predict the TKE much better away from the throat. The AMR mesh results are similar to the results displayed by the non-AMR mesh that contained 8.4 million cells throughout the calculation. Thus, while the AMR is reproducing the same cavity dynamics, it provides insights into cavitating flow with considerably much less computational time than standard refined mesh, at the local scale scale as well.  



§ CONCLUSIONS


This paper employs the Adaptive Mesh Refinement (AMR), coupled with DES model to simulate the unsteady cavitating flow in a converging-diverging (venturi) nozzle. The Merkle model is used to model for cavitation modelling. The vapor cavity structures, the mechanisms of cloud shedding and cavitation-vortex interaction are examined alongside the cavitation-turbulence coupling on the local level. The main conclusions are as follows:

    
  * The numerical simulations show that cloud cavitation is a periodic process where a cavity initiates at the throat, followed by its growth until it reaches a maximum size before being pinched-off by a re-entrant jet. Further analysis of the cavitation-vortex interaction using the terms comprising vorticity transport equation demonstrated the extensive domination of the vortex stretching and the vortex dilatation terms throughout the cavitation cycle. This indicates the  strong influence of the velocity gradients and the growth of the fluid element in the vorticity and therefore, the strong interplay between cavitation and vortex formation. The baroclinic torque is much smaller and less influential than other terms. 
    
  * Analysis of time-averaged velocities and turbulence-related aspects like, the Reynolds stress and TKE on local profile stations showed the influence of AMR as compared to a non-AMR mesh calculation. While the AMR mesh is able to predict well the velocity in stream-wise direction, it displays some differences in the wall direction velocity in the downstream profiles. However, while the numerical simulations are able to reproduce the cavity lengths as measured in experiments, there is significant difference in the cavity width at the throat, where the numerical simulations predict a thinner cavity. 
    
  * The AMR mesh calculation is able to predict the Reynolds stress data near the throat much better than the standard mesh calculation. However, the prediction worsens downstream, where the AMR calculation under-predicts the Reynolds stress. This is a consequence of the absence of low void fraction downstream, resulting in lower mesh refinement and the region being subsequently modelled by URANS. The TKE plots show both calculations display similar results and reproducing TKE downstream in a close agreement with experiments, thus highlighting the use of AMR as a tool to speed up calculations for unsteady cavitating flows. 
    
  * Future avenues include running the AMR simulation with more refinement levels to investigate the results in a better way at a reduced computational cost. However, as this study sheds more light on turbulence-cavitation interaction, future AMR calculations should focus on grid refinement based on turbulence data rather than void fraction in order to model the turbulence-cavitation interplay more accurately, throughout the domain.   





§ ACKNOWLEDGEMENTS

This work was supported by the Office of Naval Research, USA [grant number N00014-18-S-B001], the Macao Young Scholars Program (Project code: AM2022003) and Priority Postdoctoral Projects in Zhejiang Province (Project Code:341781). The authors would like to thank the ONR proposal manager Dr. Ki-Han Kim for his support. The authors would like to thank the anonymous reviewers for their feedback to enhance the quality of the manuscript. 



§ COMPLIANCE WITH ETHICAL STANDARDS


    
  * Conflict of interest: The authors declare they have no conflict of interest.
    
  * Ethical approval: The article does not contain any studies with human participants or animals performed by any of the authors.
    
  * Informed consent: Not applicable 
    
  * Funding:This work was supported by the Office of Naval Research, USA [grant number N00014-18-S-B001]. The authors would like to thank the ONR proposal manager Dr. Ki-Han Kim for his support. 









