Study on the Influence of Temperature on the Temporal and Spatial Distribution Characteristics of Natural Cavitating Flow around a Vehicle

Cavitation involves complex multiphase turbulence and has important research significance. In this study, the Schnerr–Sauer cavitation model was used to model cavitation, and the detached-eddy simulation (DES) method was used to calculate the unsteady natural cavitating flow. The predicted results are in good agreement with experimentally measured cavity evolution and pressure values, demonstrating the effectiveness of this numerical method. Low temperature causes changes in the properties of water. The density of water at 0° is 999.84 kg/m3 and the density of water at 25° is 997.04. Cavitation evolution and shedding are analyzed at temperatures of 0 °C and 25 °C. The results showed that lower temperature increased the frequency of cavitation and enhanced pressure pulsation. At the same time, low temperature also increases the frequency of cavity shedding and shortens the cycle. In addition, based on the Ω method, the difference between vortex dynamics at various temperatures was studied, and it was found that different cavity stages showed different vortex structure characteristics, and lower temperature would aggravate the change of wake vortex structure. At the same time, the analysis of the turbulence characteristics in the downstream of the cavity shows that the lower temperature reduces the velocity pulsation and reduces the turbulence integral scale. At the end of the model, large-scale pulsations are transformed into small-scale pulsations.


Introduction
Cavitation usually occurs when the pressure in a liquid stream drops below the local saturated vapor pressure, causing cavities to be generated at the internal or liquid-solid interface [1,2]. Cavitating flows are typically unsteady and exhibit phase changes, turbulence, and multi-scale vortices. Such flows occur in various fluidic machines, including liquid rocket engines, turbines, marine propellers, and hydrofoils. Cavitation can cause serious problems in fluid machinery, such as pressure fluctuations, noise, vibration, and corrosion. Therefore, cavitating flows continue to be widely studied [3][4][5][6].
To date, many experimental and numerical studies have investigated the flow characteristics of cavitation [7][8][9][10]. Knapp [11] observed and analyzed cloud cavitation flow and pointed out that the re-entrant flow from the closed cavity position is the cause of cloud cavitation ruptures. Reissman et al. [12] studied cloud-like cavitation around an oscillating hydrofoil using high-speed cameras and pressure sensors and investigated the mechanism of transient pressure pulses near the detached cavitation. Passandideh-Fard and Roohi [13] used an improved volume-of-fluid (VOF) method to capture the gas-liquid interface of a two-dimensional cavitating flow field. It was found that the improved VOF method offers higher calculation efficiency and accuracy. Ganesh et al. [14] studied how wedge-shaped cavitation evolved into cloud cavitation through wedge geometric model experiments. They found two types of cavity shedding mechanisms in the cavitating flow field: the re-entrant jet mechanism and the shock wave mechanism. The occurrence and collapse of cavitation are periodic. This periodic pressure shock will change the dynamic characteristics of the flow field, also causing strong vibration and noise. In addition, the reentrant jets and shock waves formed during cavitation collapse can cause high-pressure pulsations on the surface of the material [15][16][17]. Wang et al. [18] performed numerical calculations of the transient characteristics of cloud cavitating flows and shock wave dynamics. They reported several valuable results regarding the influence of water/vapor compressibility in the cavitation instabilities and the relation between cavity development and pressure fluctuations. Ji et al. [19] conducted a numerical study of the cavitating flow structure around a twisting hydrofoil. Three flow characteristics on the suction side of the twisted hydrofoil were analyzed, and the vorticity transport equation in the variable density flow was used to illustrate the effect of cavitation on the vorticity distribution. De et al. [20] studied the structure of the unsteady cavitation flow on the hydrofoil through appropriate orthogonal decomposition (POD) and fast Fourier transform (FFT) techniques, focusing on analyzing the effects of DCM and PBE models on three different cavitation methods (bubble cavitation, cloud cavitation and supercavitation). Kadivar et al. [21] applied the passive flow control method to numerically and experimentally study the unsteady cavitation fluctuation around the semicircular front edge plate. The research found that the passive control method can better prevent the development of the attachment cavity's spanwise instability and mitigate the large cavity structure. Timoshevskiy et al. [22] used PIV technology to analyze the cavitation state from the velocity field of the 2D symmetric hydrofoil cavitation under four velocity conditions. The data showed that the turbulence distribution and turbulence distribution of turbulence changed significantly with the development of cavitation. In particular, the behavior of the average speed and the most probable speed unexpectedly seems to be different.
The physical properties of fluids are related to temperature, and changes in fluid properties further affect the evolution of the cavity and the flow field [23,24]. Typical thermo-fluids include cryogenic fluids, refrigerants, fluoroketone, and high-temperature water. Hall et al. [25] conducted comprehensive experiments with four cavitation test models in water and Freon 113 at different inlet speeds and temperatures. Cervone et al. [26] conducted experiments in water around NACA0015 hydrofoils at different cavitation numbers and free flow temperatures (298-343 K) to study the characteristics of cavitation instabilities. Franc et al. [27] studied the cavitation instability of the refrigerant R-114 at three different temperatures and proved that a higher reference temperature delayed the occurrence of blade cavitation. Petkovšek and Dular [28] conducted experimental research on the temperature field of cavitation in hot water and measured the temperature dynamics under various operating conditions. Chen et al. [29] studied the fluoroketone phase transfer process by evaluating the temperature and speed of various free flows. These numerical studies have provided accurate predictions and insights into the characteristics of cavitating flow in thermosensitive fluids. Sun et al. [30] studied an unsteady cavitating flow on the NACA0015 hydrofoil in a fluoroketone thermosensitive fluid, with particular emphasis on the cavitation shedding dynamics. Studies have shown that the temperature of the hydrofoil increases and decreases under certain thermal conditions. Thermal effects suppress cavitation and have the potential to reduce pressure pulsation peaks. The current research on temperature mainly focuses on low-temperature fluids. In a low-temperature water environment, the mechanism of changes in the cavity evolution and vortex structure is still unclear. Wang et al. [31] used OpenFOAM to perform a 3-D numerical simulation of the transient sheet/cloud cavitation flow around a Clark-Y hydrofoil. The results show that the compressible method can predict the reentrant jet dynamics well. The velocity divergence comes from the compressibility and mass transfer of the cavitation two-phase fluid, and also depends on the fluid density. At the same time, the temperature and density changes in different cavitation structures are given.
In cavitating flow simulations, the selection of the turbulence model plays an important role, and the numerical prediction results produced by different turbulence models will typically vary [32][33][34][35]. At present, the Reynolds-averaged Navier-Stokes (RANS) method is very popular in practical engineering, and most numerical simulations use the RANS method for cavitating flow calculations. Combining the VOF model and the k-ω SST turbulence model, An et al. [36] performed a numerical simulation of the supercavitating flow field formed on the submerged body of ground vehicles during the artificial ventilation process and analyzed the influence factors of resistance. Despite its good stability, reasonable computational cost, and acceptable accuracy, RANS overpredicts the turbulent eddy viscosity. In recent years, with the development of computer technology, the large-eddy simulation (LES) framework has increasingly been used to simulate cavitating flows [37][38][39][40]. Huang et al. [41] used the LES method to numerically study the cavitating flow on the flexible NACA66 hydrofoil and analyzed the vibration and cavitation vortex structure caused by the flow. The results show that the change trend of the cavity is consistent with the evolution of the vortex, and the vortex gradually increases in the wake. LES can predict the unstable characteristics of the flow [42,43], but requires a very fine mesh and small time steps. Therefore, the computational cost of the numerical simulations is very high. As a compromise, separation vortices have been widely used in large-separation unsteady flows in recent years. Kunz et al. [44] simulated an unsteady cavitating flow on a hydrofoil using the RANS and detached-eddy simulation (DES) methods, comparing the evolution of the cavity and the changes in lift and drag. The separation vortex approach combines the advantages of RANS and LES: the RANS method is used for simulating the near-wall area to reduce the grid size and calculation time, with the LES method used to simulate the area away from the wall to capture large-scale separated flows [45][46][47].
Although a lot of encouraging work has been conducted on cavitation, further research is still needed. Therefore, this article introduces a numerical study of rotating bodies at two temperatures. What is different from the past is that we pay attention to the influence of low temperature water environment on the temporal and spatial characteristics of cavitation flow. The main purpose is to study the effect of changes in fluid physical properties on cavitation. Section 2 introduces the numerical simulation models, including the basic governing equations, cavitation models, turbulence models, and verification of the numerical methods. Section 3 discusses and analyzes the influence of temperature on the cavity evolution, period, pressure fluctuations, vortex structure, and turbulence scale. Finally, Section 4 summarizes the results of this study and states the main conclusions.

Governing Equations
In this study, the Euler multiphase flow model was used to simulate natural cavitating flow. By solving a set of momentum equations over the entire watershed and tracking the volume fraction of each phase in the finite volume of each fluid, the two-phase immiscible flow was simulated. The mixing density is expressed as: where ρ is the mixture density and α l represents the water volume fraction. For the multiphase cavitating flow process, the mass and momentum of the movement of matter are conserved. This is expressed as follows: 4 of 20 where p is the pressure and u represents velocity; the subscripts i and j indicate the direction of the velocity component.

Detached-Eddy Simulation Model
The standard Spalart-Allmaras model defines the turbulence scale as the distance to the nearest wall, d. The original DES model proposed by Spalart et al. [48] was modified using the length scale: where C DES is a calibrated constant that functions similar to the Smagorinsky constant in the sub-grid stress model. ∆ is based on the largest dimension of the grid cell: Based on the shear-stress transport (SST) model, Menter et al. [49] developed a twoequation SST-DES method and modified the dissipation term of the k equation. The model equation is: where F DES is the scale factor of the vortex separation method, given by: and the calculated turbulence length scale L t is:

Schnerr-Sauer Cavitation Model
The Schnerr-Sauer cavitation model is composed of transport equations. This model is obtained by simplifying and improving the Rayleigh-Plesset cavity dynamics equation (Rayleigh,[50]; Plesset, [51]). The self-simplified Rayleigh-Plesset equation can be expressed as: where P V is the saturation pressure at the local temperature, P is the pressure around the cavity, R B denotes the bubble radius, and ξ is the surface tension between liquid and vapor. Ignoring second-order terms and surface tension, this equation can be simplified to: The processes of evaporation and condensation are expressed by inserting a source term into the transmission equation to describe the phase transition rate. The Schnerr-Sauer model is used as the cavitation model [52]. The transfer equation is: .
where . m + is the evaporation term, . m − is the condensation term, and P v is the saturation vapor pressure at the local temperature. R b is the cavity radius, which is given by: (15) where N b is the density of cavity number. According to Schnerr and Sauer [52], N b = 10 13 . Figure 1 shows a schematic diagram of the computational geometric model. It has a conical axisymmetric body with a diameter of D = 40 mm and a length of L = 335 mm. The length of the head is 19 mm, the length of the cylinder is 316 mm, and the angle between the cone section of the head and the cylinder is 45 • . Four pressure monitoring points are set up to verify the numerical method and analyze the results. P1 is 100 mm from the model head, P2 is 12 mm from P1, P3 is 8 mm from P2, and P4 is 22 mm from P3. The P2 and P4 points correspond to the experimental monitoring points to verify the validity of the numerical simulation, and the P1 and P3 points, respectively, detect the pressure pulsation inside the cavity and the drop position.

Numerical Setup
where ̇+ is the evaporation term, ̇− is the condensation term, and is the saturation vapor pressure at the local temperature.
is the cavity radius, which is given by: where Nb is the density of cavity number. According to Schnerr and Sauer [52], Nb = 10 13 . Figure 1 shows a schematic diagram of the computational geometric model. It has a conical axisymmetric body with a diameter of D = 40 mm and a length of L = 335 mm. The length of the head is 19 mm, the length of the cylinder is 316 mm, and the angle between the cone section of the head and the cylinder is 45°. Four pressure monitoring points are set up to verify the numerical method and analyze the results. P1 is 100 mm from the model head, P2 is 12 mm from P1, P3 is 8 mm from P2, and P4 is 22 mm from P3. The P2 and P4 points correspond to the experimental monitoring points to verify the validity of the numerical simulation, and the P1 and P3 points, respectively, detect the pressure pulsation inside the cavity and the drop position.  Figure 2 shows the size of the calculation domain around the calculation model. The relative position of the calculation model in the calculation domain is also shown. The distance between the entrance to the calculation domain and the front end of the model is 15D, the distance between the exit and the rear edge of the model is 20D, and the distance between the upper wall and the lower wall is 10D. The entrance to the calculation domain is set as a velocity inlet, the surroundings are set as symmetry planes, and the outlet is assigned a pressure outlet boundary condition. The entry to the calculation domain is set to an inflow velocity of V∞ = 30 m/s, giving a Reynolds number of Re = 8.9 × 10 6 . The cavitation number changes according to the temperature (Brennen, [53]), and is defined as:  Figure 2 shows the size of the calculation domain around the calculation model. The relative position of the calculation model in the calculation domain is also shown. The distance between the entrance to the calculation domain and the front end of the model is 15D, the distance between the exit and the rear edge of the model is 20D, and the distance between the upper wall and the lower wall is 10D. The entrance to the calculation domain is set as a velocity inlet, the surroundings are set as symmetry planes, and the outlet is assigned a pressure outlet boundary condition. The entry to the calculation domain is set to an inflow velocity of V ∞ = 30 m/s, giving a Reynolds number of Re = 8.9 × 10 6 . The cavitation number changes according to the temperature (Brennen, [53]), and is defined as:

Numerical Setup
To study the characteristics of turbulence at different temperatures, detection surfaces were set up inside the cavity and downstream of the trailing edge of the model, as shown in Figure 3a. The detection surface is composed of 13 evenly spaced concentric circles, with 90 monitoring points distributed on each circle. Figure 3b shows the rear view of the monitoring surface. Figure 4 shows the calculation grid. To capture the detailed features of unstable flow, the surface of the model has a prism layer and the front end of the model has a higher resolution to better capture the cavitation evolution. In each time step, performing too many iterations will increase the calculation time without noticeably improving the accuracy, whereas performing too few iterations may lead to poor convergence and inaccurate calculations. In the simulations, 10 iterations were performed in each time step to ensure a suitable balance between accuracy and efficiency. The grid contained approximately 8 × 10 7 cells and the calculated y+ was less than 1. To study the characteristics of turbulence at different temperatures, detection surfaces were set up inside the cavity and downstream of the trailing edge of the model, as shown in Figure 3a. The detection surface is composed of 13 evenly spaced concentric circles, with 90 monitoring points distributed on each circle. Figure 3b shows the rear view of the monitoring surface.  Figure 4 shows the calculation grid. To capture the detailed features of unstable flow, the surface of the model has a prism layer and the front end of the model has a higher resolution to better capture the cavitation evolution. In each time step, performing too many iterations will increase the calculation time without noticeably improving the accuracy, whereas performing too few iterations may lead to poor convergence and inaccurate calculations. In the simulations, 10 iterations were performed in each time step to ensure To study the characteristics of turbulence at different temperatures, detection surfaces were set up inside the cavity and downstream of the trailing edge of the model, as shown in Figure 3a. The detection surface is composed of 13 evenly spaced concentric circles, with 90 monitoring points distributed on each circle. Figure 3b shows the rear view of the monitoring surface.  Figure 4 shows the calculation grid. To capture the detailed features of unstable flow, the surface of the model has a prism layer and the front end of the model has a higher resolution to better capture the cavitation evolution. In each time step, performing too many iterations will increase the calculation time without noticeably improving the accuracy, whereas performing too few iterations may lead to poor convergence and inaccurate calculations. In the simulations, 10 iterations were performed in each time step to ensure a suitable balance between accuracy and efficiency. The grid contained approximately 8 × 10 7 cells and the calculated y+ was less than 1.

Numerical Method Validation
The experiments mainly investigated a circulating high-speed ventilated cavitation in the water tunnel of the Harbin Institute of Technology. The experimental system has been described in detail in a previous article (Zhang et al. [54]). In the process of numerical simulation, a steady field is given to the flow field. After the flow field is stable, the cavi-

Numerical Method Validation
The experiments mainly investigated a circulating high-speed ventilated cavitation in the water tunnel of the Harbin Institute of Technology. The experimental system has been described in detail in a previous article (Zhang et al. [54]). In the process of numerical simulation, a steady field is given to the flow field. After the flow field is stable, the cavitation model is turned on to simulate the evolution of the cavity. In the test conditions considered in this study, the working pressure was 68.4 kPa, the flow velocity of the water tunnel is 10.21 m/s, and the cavitation number was 0.62. Figure 5a,b compare the instantaneous cavity shape in the experiments and the numerical predictions. At the same time, Figure 5c illustrates the volume fraction distribution of the XOY plane. At t 0 + 6.7 ms, at the front end of the cavity, a gap appears at the shoulder of the model, as shown in Figure 6a,b. From t 0 + 6.7 ms to t 0 + 13.4 ms, large-scale shedding occurs over the entire cavity and moves downstream. This is the main shedding process and causes a large pressure pulsation downstream. At t 0 + 20.1 ms, the large-scale shedding ends, and the front-end cavitation develops again. At t 0 + 26.8 ms, the length of the cavity is the same as at t 0 . As shown in Figure 5c, there is a gap at the shoulder of the cavity at t 0 + 6.7 ms. At t 0 + 13.4 ms, under the influence of the re-entrant flow, a large cavity detaches from the entire cavity on the main body. At t 0 + 20.1 ms, the shed cavity moves downstream and gradually collapses. The results of the experiment are in good agreement with those from the simulations.

Numerical Method Validation
The experiments mainly investigated a circulating high-speed ventilated cavitation in the water tunnel of the Harbin Institute of Technology. The experimental system has been described in detail in a previous article (Zhang et al. [54]). In the process of numerical simulation, a steady field is given to the flow field. After the flow field is stable, the cavitation model is turned on to simulate the evolution of the cavity. In the test conditions considered in this study, the working pressure was 68.4 kPa, the flow velocity of the water tunnel is 10.21 m/s, and the cavitation number was 0.62. Figure 5a,b compare the instantaneous cavity shape in the experiments and the numerical predictions. At the same time, Figure 5c illustrates the volume fraction distribution of the XOY plane. At t0 + 6.7 ms, at the front end of the cavity, a gap appears at the shoulder of the model, as shown in Figure  6a,b. From t0 + 6.7 ms to t0 + 13.4 ms, large-scale shedding occurs over the entire cavity and moves downstream. This is the main shedding process and causes a large pressure pulsation downstream. At t0 + 20.1 ms, the large-scale shedding ends, and the front-end cavitation develops again. At t0 + 26.8 ms, the length of the cavity is the same as at t0. As shown in Figure 5c, there is a gap at the shoulder of the cavity at t0 + 6.7 ms. At t0 + 13.4 ms, under the influence of the re-entrant flow, a large cavity detaches from the entire cavity on the main body. At t0 + 20.1 ms, the shed cavity moves downstream and gradually collapses. The results of the experiment are in good agreement with those from the simulations.  Figure 6 compares the experimental and numerical pressure curves at monitoring points P2 and P4. The large-scale shedding of the cavitation bubbles moves downstream. The pressure fluctuations are obvious, first decreasing and then increasing as the cavitation bubbles pass. These results verify that the numerical method proposed in this paper can effectively predict the flow field evolution and pressure pulsation of cavitating flow.
In this study, we use STAR-CCM+ to conduct numerical simulation research; the time step was set to ∆t = 2.5 × 10 −5 s so that the Courant number was less than 1 at every position in the computational domain, thus ensuring good convergence. The numerical simulations were based on the finite volume method, and the VOF model adopted a second-order implicit scheme. The governing equations were subjected to second-order Figure 6 compares the experimental and numerical pressure curves at monitoring points P2 and P4. The large-scale shedding of the cavitation bubbles moves downstream. The pressure fluctuations are obvious, first decreasing and then increasing as the cavitation bubbles pass. These results verify that the numerical method proposed in this paper can effectively predict the flow field evolution and pressure pulsation of cavitating flow. In this study, we use STAR-CCM+ to conduct numerical simulation research; the time step was set to ∆ = 2.5 × 10 −5 s so that the Courant number was less than 1 at every position in the computational domain, thus ensuring good convergence. The numerical simulations were based on the finite volume method, and the VOF model adopted a second-order implicit scheme. The governing equations were subjected to second-order discretization in space. The split solver based on the SIMPLE algorithm was used to handle the coupling between the velocity and the pressure. Figure 7 shows the typical cavity morphology with a water volume fraction of 0.8 at 16 points during the evolution process at 0 °C and 25 °C. The initial position of natural cavity is at the position of the dotted line a. From t0 to t0 + 6.0 ms, the volume of the cavitation bubble increases to its maximum size. With the generation of the re-entrant jet, the cavity downstream gradually shedding. The volume of the cavitation bubble gradually decreases until it collapses to the position where the initial cavitation bubble is generated. At t0 + 6.0 ms, the cavitation length is greater in the 25 °C flow than in the 0 °C flow, and a vortex structure is generated at the tail. At t0 + 7.5 ms, the shedding of the vortex structure reduces the cavitation size, and we can determine that the collapse rate of the 0 °C cavitation bubble is faster than in the 25 °C case by comparing dotted lines b and c. At t0 + 21.0 ms, the 0 °C cavitation bubble has collapsed to its smallest volume, and the system has entered the next cycle. The 25 °C cavitation bubble reaches its minimum size at t0 + 22.5 ms, which is about 1.5 ms slower than at the lower temperature. The density of water at 0 °C is 999.84kg/m 3 and the density of water at 25 °C is 997.04 kg/m 3 . The saturated vapor pressure at 0 °C is 611 Pa, whereas that at 25 °C is 3169 Pa. The cavitation number of the 25 °C flow is 0.328, and that in the 0 °C case is 0.333. Therefore, the void volume at 25 °C is greater than that at 0 °C. This indicates that a lower temperature increases the instability of cavitation bubbles, accelerates the internal phase transition rate of the cavity, and shortens the evolution period.  Figure 7 shows the typical cavity morphology with a water volume fraction of 0.8 at 16 points during the evolution process at 0 • C and 25 • C. The initial position of natural cavity is at the position of the dotted line a. From t 0 to t 0 + 6.0 ms, the volume of the cavitation bubble increases to its maximum size. With the generation of the re-entrant jet, the cavity downstream gradually shedding. The volume of the cavitation bubble gradually decreases until it collapses to the position where the initial cavitation bubble is generated. At t 0 + 6.0 ms, the cavitation length is greater in the 25 • C flow than in the 0 • C flow, and a vortex structure is generated at the tail. At t 0 + 7.5 ms, the shedding of the vortex structure reduces the cavitation size, and we can determine that the collapse rate of the 0 • C cavitation bubble is faster than in the 25 • C case by comparing dotted lines b and c. At t 0 + 21.0 ms, the 0 • C cavitation bubble has collapsed to its smallest volume, and the system has entered the next cycle. The 25 • C cavitation bubble reaches its minimum size at t 0 + 22.5 ms, which is about 1.5 ms slower than at the lower temperature. The density of water at 0 • C is 999.84 kg/m 3 and the density of water at 25 • C is 997.04 kg/m 3 . The saturated vapor pressure at 0 • C is 611 Pa, whereas that at 25 • C is 3169 Pa. The cavitation number of the 25 • C flow is 0.328, and that in the 0 • C case is 0.333. Therefore, the void volume at 25 • C is greater than that at 0 • C. This indicates that a lower temperature increases the instability of cavitation bubbles, accelerates the internal phase transition rate of the cavity, and shortens the evolution period.

Cavity Evolution and Pressure Pulsation
To analyze the evolution of cavitation bubbles, Figure 8 shows the changes in cavity volume at the two temperatures. The volume of the cavity changes periodically and oscillates within a certain interval. As the cavity collapses and sheds, the volume of the cavity reaches its minimum value, from which point it begins to grow to its maximum volume once again. Interestingly, the volume change of cavitation bubbles at 25 • C is smoother than at 0 • C, and the volume change in each cycle is also larger. In some periods, the maximum volume of the 0 • C case exceeds the void volume of the 25 • C case. The average volume of the cavity at 0 • C is 1.41 × 10 −4 m 3 , compared with 1.43 × 10 −4 m 3 at 25 • C. The average hole volume at 25 • C is slightly larger than that at 0 • C and there are a greater number of holes. This shows that the lower temperature makes the cavity smaller and accelerates the collapse and shedding processes. To analyze the evolution of cavitation bubbles, Figure 8 shows the changes in cavity volume at the two temperatures. The volume of the cavity changes periodically and oscillates within a certain interval. As the cavity collapses and sheds, the volume of the cavity reaches its minimum value, from which point it begins to grow to its maximum volume once again. Interestingly, the volume change of cavitation bubbles at 25 °C is smoother than at 0 °C, and the volume change in each cycle is also larger. In some periods, the maximum volume of the 0 °C case exceeds the void volume of the 25 °C case. The average volume of the cavity at 0 °C is 1.41 × 10 −4 m 3 , compared with 1.43 × 10 −4 m 3 at 25 °C. The average hole volume at 25 °C is slightly larger than that at 0 °C and there are a greater number of holes. This shows that the lower temperature makes the cavity smaller and accelerates the collapse and shedding processes.  Figure 9a shows the pressure change at two positions at the two temperatures. The pressure change is also periodic. When the cavity reaches position P1, the minimum pressure at each temperature is the saturated vapor pressure. With the collapse and shedding of the cavity bubbles, the pressure fluctuates violently. When the shedding cavity moves to position P3, pressure pulsations can be detected, and the pulsation period is the same as the cavity collapse and shedding period. At the two points, the average pulsation pressure in the 0 °C case is higher than in the 25 °C case, indicating that the pressure pulsations To analyze the evolution of cavitation bubbles, Figure 8 shows the changes in cavity volume at the two temperatures. The volume of the cavity changes periodically and oscillates within a certain interval. As the cavity collapses and sheds, the volume of the cavity reaches its minimum value, from which point it begins to grow to its maximum volume once again. Interestingly, the volume change of cavitation bubbles at 25 °C is smoother than at 0 °C, and the volume change in each cycle is also larger. In some periods, the maximum volume of the 0 °C case exceeds the void volume of the 25 °C case. The average volume of the cavity at 0 °C is 1.41 × 10 −4 m 3 , compared with 1.43 × 10 −4 m 3 at 25 °C. The average hole volume at 25 °C is slightly larger than that at 0 °C and there are a greater number of holes. This shows that the lower temperature makes the cavity smaller and accelerates the collapse and shedding processes.  Figure 9a shows the pressure change at two positions at the two temperatures. The pressure change is also periodic. When the cavity reaches position P1, the minimum pressure at each temperature is the saturated vapor pressure. With the collapse and shedding of the cavity bubbles, the pressure fluctuates violently. When the shedding cavity moves to position P3, pressure pulsations can be detected, and the pulsation period is the same as the cavity collapse and shedding period. At the two points, the average pulsation pressure in the 0 °C case is higher than in the 25 °C case, indicating that the pressure pulsations  Figure 9a shows the pressure change at two positions at the two temperatures. The pressure change is also periodic. When the cavity reaches position P1, the minimum pressure at each temperature is the saturated vapor pressure. With the collapse and shedding of the cavity bubbles, the pressure fluctuates violently. When the shedding cavity moves to position P3, pressure pulsations can be detected, and the pulsation period is the same as the cavity collapse and shedding period. At the two points, the average pulsation pressure in the 0 • C case is higher than in the 25 • C case, indicating that the pressure pulsations caused by cavitation are more obvious at lower temperatures. This is because the saturated vapor pressure at 0 • C is 611 Pa and that at 25 • C is 3169 Pa. At the same flow rate, the 25 • C cavity is more stable, shedding and collapsing more smoothly, than in the 0 • C case, demonstrating the effect of the pressure pulsations on the shedding and collapse processes. Figure 10 shows the power spectral density (PSD) frequency distribution diagram at pressure detection points P1 and P3 (sampling frequency is 0-2.5 × 10 4 .). The main frequency of natural shedding at 0 • C is 50.0 Hz, compared with 45.8 Hz at 25 • C. It proves that low temperature increases the shedding frequency, that is, the cycle of cavitation evolution is shortened, and the cavitation evolution is intensified, which also corresponds to the result in Figure 8.
processes. Figure 10 shows the power spectral density (PSD) frequency distribution diagram at pressure detection points P1 and P3 (sampling frequency is 0-2.5 × 10 4 .). The main frequency of natural shedding at 0 °C is 50.0 Hz, compared with 45.8 Hz at 25 °C. It proves that low temperature increases the shedding frequency, that is, the cycle of cavitation evolution is shortened, and the cavitation evolution is intensified, which also corresponds to the result in Figure 8.  that low temperature increases the shedding frequency, that is, the cycle of cavitation evolution is shortened, and the cavitation evolution is intensified, which also corresponds to the result in Figure 8.

Vortex Structure Characteristics
Turbulence produces a large number of vortex structures with different scales and strengths. These vortex structures play a key role in the generation and maintenance of turbulence and are important in the analysis of cavitation problems. Therefore, it is of great significance to analyze the vortex structure. The third-generation Ω method (Liu [55]) for vortex identification offers enhanced vortex recognition without the need to adjust the threshold according to the flow. Here, the vorticity represents the curl of the velocity: We can decompose the velocity gradient tensor as follows: where A is the strain rate tensor, A = 1 2 ∇V + ∇V T , B is the rotation tensor, B = 1 2 ∇V − ∇V T , A F B F are the Frobenius norm of the matrix, and ε is the correction value: where ε c = 10 −7 is an empirical constant (Liu et al. [55]). Ω represents the ratio of the vorticity of the rotating part to the total vorticity, which is expressed as Figure 11a shows the vortex structure at the two temperatures during the development stage of the cavity. A strong vortex structure can be found on the surfaces of the two cavities. As the cavity grows, some vortex structures become attached to the surface of the cavity. At both temperatures, the Ω value inside the cavity zone is smaller than that outside the cavity. Inside the cavity, the shearing effect is equivalent to the rotation, while the rotation effect on the surface of the cavity indicates that there is strong rotational movement in this area. As shown by the black dotted line, a strong vortex structure is concentrated at the tail of the model, and the vortex structure is the same at the two temperatures. In the wake area, a large number of linear vortex structures are produced, and these are stronger in the 0 • C case than in the 25 • C case. The change period of 0 • C cavity is faster than that of 25 • C cavity, so the frequency of 0 • C cavity generation and shedding is greater, and a large number of vortex structures move to the wake region. As a result, the 0 • C vortex structure in the wake region is stronger than the 25 • C vortex structure. Figure 11b shows the vortex structure at the two temperatures at the stage of cavity stability. The vortex structure attached to the surface of the cavity has disappeared, and a strong vortex structure has appeared inside the cavity. The vortex structure at 0 • C is stronger than that at 25 • C. Due to the generation of the back jet, a large amount of vortex motion appears inside the cavity. The 0 • C cavity is more unstable, and its shedding frequency is faster, leading to the stronger vortex structure and the shedding of the next stage of a cavity. At the stage of cavity stability, no vortex shedding occurs and there is no disturbance in the wake region. Therefore, the linear vortex structures in the wake region are the same at both temperatures. Figure 11c shows the stage of cavity shedding. Several vortex structures are attached to the surface of the cavity, but the vortex structure is weaker than during the cavity development stage. The vortex structure attached to the rear wall of the model is stronger than in the first two stages. This vortex structure is larger in the 0 • C case than at 25 • C. At low temperatures, the vortex structure on the model wall is enhanced. The cavity breaks and sheds, and a large number of vortices move along the wall surface wake. The cavity shedding is faster at 0 • C than at 25 • C, so the wall vortex structure in the back section of the model is stronger at 0 • C.
Comparing the wake field of the three stages, it is found that the vortex structure changes significantly at 0 • C, especially during the development of cavity, and the tail vorticity exhibits a divergent trend. At 25 • C, the overall wake field vortex structure changes more smoothly, although the scale of the vortex structure is larger than at 0 • C. This is because the cavity shedding frequency is lower at 25 • C than at 0 • C, so the vortex structure has a large scale and a small overall strength.
Comparing the wake field of the three stages, it is found that the vortex structure changes significantly at 0 °C, especially during the development of cavity, and the tail vorticity exhibits a divergent trend. At 25 °C, the overall wake field vortex structure changes more smoothly, although the scale of the vortex structure is larger than at 0 °C. This is because the cavity shedding frequency is lower at 25 °C than at 0 °C, so the vortex structure has a large scale and a small overall strength.

Cavitation surface Cavitation surface
Wake field Wake field  Figure 12 compares the vortex structures produced in three-dimensional natural cavitation evolution at 0 • C and 25 • C. Due to the growth and shedding of the cavity, the vortex structure is always attached to the non-cavitated area of the rear wall of the model, and the scale varies with time. At t 0 + 0 ms, there are many linear vortex structures in the wake field. The vortex structure is stronger at 0 • C than at 25 • C, which is the same conclusion as for Figure 11a. The volume of the cavity has increased at t 0 + 10 ms, resulting in a smooth Ω isosurface in the cavity head and a small velocity gradient, indicating laminar flow. The re-entrant jet causes the flow to undergo a rapid laminar-turbulent transition in the closed position of the cavity, forming a series of vortex structures. These vortex structures move with the direction of the incoming flow and develop into a series of ring-like vortices. At this moment, due to the closure and shedding of the cavity, an annular vortex structure appears in the latter part of the model at 0 • C and at the position of the cavity shedding at 25 • C. As the vortex structure moves downstream, at t 0 + 15 ms, the tail end presents an obvious annular vortex structure. According to the cavity evolution process, the cavity shedding frequency is faster at 0 • C than at 25 • C. Therefore, there is a time difference between the vortex structures at the two temperatures. In general, the vortex structure scale is smaller at 0 • C than at 25 • C. However, there are more vortex structures at the lower temperature, especially in the wake field area of the model. Figure 12 compares the vortex structures produced in three-dimensional natural cavitation evolution at 0 °C and 25 °C. Due to the growth and shedding of the cavity, the vortex structure is always attached to the non-cavitated area of the rear wall of the model, and the scale varies with time. At t0 + 0 ms, there are many linear vortex structures in the wake field. The vortex structure is stronger at 0 °C than at 25 °C, which is the same conclusion as for Figure 11a. The volume of the cavity has increased at t0 + 10 ms, resulting in a smooth Ω isosurface in the cavity head and a small velocity gradient, indicating laminar flow. The re-entrant jet causes the flow to undergo a rapid laminar-turbulent transition in the closed position of the cavity, forming a series of vortex structures. These vortex structures move with the direction of the incoming flow and develop into a series of ring-like vortices. At this moment, due to the closure and shedding of the cavity, an annular vortex structure appears in the latter part of the model at 0 °C and at the position of the cavity shedding at 25 °C. As the vortex structure moves downstream, at t0 + 15 ms, the tail end presents an obvious annular vortex structure. According to the cavity evolution process, the cavity shedding frequency is faster at 0 °C than at 25 °C. Therefore, there is a time difference between the vortex structures at the two temperatures. In general, the vortex structure scale is smaller at 0 °C than at 25 °C. However, there are more vortex structures at the lower temperature, especially in the wake field area of the model.

Analysis of Turbulence Intensity and Integral Scale
Turbulence characteristics are important to the evolution of cavitation. To study the effect of temperature on the turbulence of natural cavitation and the turbulence integral scale, we established detection surfaces inside the cavity and downstream of the cavity shedding at the two temperatures. The turbulent flow characteristics considered in this paper are based on the assumption of random stationary flow. The time series of the x-axis velocity at time t is used to find the time average. The axial pulsation velocity of the corresponding point can be expressed as: The turbulence intensity is expressed as: where U ref is the flow speed at infinity and ε x is the turbulence intensity in the x-direction. Figure 13a shows a schematic diagram of the two-dimensional turbulence intensity inside the cavity at the two temperatures. The inside of the cavity exhibits a high degree of turbulence, while the boundary layer area near the model wall has a lower level of turbulence. Due to the violent phase change process in the cavity, the water is transformed into vapor and the velocity changes, which is characterized by higher turbulence intensity and strong instability. Figure 13b shows the three-dimensional turbulence intensity inside the cavity at the two temperatures. The inside of the cavity again exhibits a high degree of turbulence, while the boundary layer area near the model wall shows a lower degree of turbulence. A violent phase transition process occurs inside the cavity. Due to the cavitation effect, the water transforms into vapor and the velocity changes. Therefore, it is characterized by high turbulence intensity and strong instability. It can be seen that the turbulence intensity at 25 • C is 0.3 and that at 0 • C is 0.28, a decrease of approximately 6%. Furthermore, under low-temperature conditions, the internal phase change in the cavity is violent and the cavity becomes unstable. effect of temperature on the turbulence of natural cavitation and the turbulence integral scale, we established detection surfaces inside the cavity and downstream of the cavity shedding at the two temperatures. The turbulent flow characteristics considered in this paper are based on the assumption of random stationary flow. The time series of the xaxis velocity at time t is used to find the time average. The axial pulsation velocity of the corresponding point can be expressed as: The turbulence intensity is expressed as: where Uref is the flow speed at infinity and is the turbulence intensity in the x-direction. Figure 13a shows a schematic diagram of the two-dimensional turbulence intensity inside the cavity at the two temperatures. The inside of the cavity exhibits a high degree of turbulence, while the boundary layer area near the model wall has a lower level of turbulence. Due to the violent phase change process in the cavity, the water is transformed into vapor and the velocity changes, which is characterized by higher turbulence intensity and strong instability. Figure 13b shows the three-dimensional turbulence intensity inside the cavity at the two temperatures. The inside of the cavity again exhibits a high degree of turbulence, while the boundary layer area near the model wall shows a lower degree of turbulence. A violent phase transition process occurs inside the cavity. Due to the cavitation effect, the water transforms into vapor and the velocity changes. Therefore, it is characterized by high turbulence intensity and strong instability. It can be seen that the turbulence intensity at 25 °C is 0.3 and that at 0 °C is 0.28, a decrease of approximately 6%. Furthermore, under low-temperature conditions, the internal phase change in the cavity is violent and the cavity becomes unstable.  Figure 14a shows the two-dimensional turbulence intensity of the wake field at t two temperatures. The higher-turbulence area is mainly concentrated in the tail flow fiel and its radius is the same as the model radius. This flow phenomenon occurs because t velocity of the flow field in the tail region of the model is small, whereas the velocity o both sides is large. When the two velocities meet in the wake field, a large velocity puls tion occurs. Figure 14b shows the three-dimensional turbulence intensity at the dow  Figure 14a shows the two-dimensional turbulence intensity of the wake field at the two temperatures. The higher-turbulence area is mainly concentrated in the tail flow field, and its radius is the same as the model radius. This flow phenomenon occurs because the velocity of the flow field in the tail region of the model is small, whereas the velocity on both sides is large. When the two velocities meet in the wake field, a large velocity pulsation occurs. Figure 14b shows the three-dimensional turbulence intensity at the downstream cavity shedding position. The turbulence intensity at 25 • C is 0.175, whereas that at 0 • C is 0.15, a reduction of 14%. Compared with the inside of the cavity, the turbulence intensity is significantly reduced at both temperatures. This is because, as the cavity sheds, the vortex structure is gradually formed and the turbulent dissipation rate increases, causing the turbulent kinetic energy to be converted into heat energy. This is manifested as an overall decrease in turbulence intensity. The vortex structure sheds at both temperatures, and so the turbulence intensity decreases even more.  Figure 14a shows the two-dimensional turbulence intensity of the wake field at the two temperatures. The higher-turbulence area is mainly concentrated in the tail flow field, and its radius is the same as the model radius. This flow phenomenon occurs because the velocity of the flow field in the tail region of the model is small, whereas the velocity on both sides is large. When the two velocities meet in the wake field, a large velocity pulsation occurs. Figure 14b shows the three-dimensional turbulence intensity at the downstream cavity shedding position. The turbulence intensity at 25 °C is 0.175, whereas that at 0 °C is 0.15, a reduction of 14%. Compared with the inside of the cavity, the turbulence intensity is significantly reduced at both temperatures. This is because, as the cavity sheds, the vortex structure is gradually formed and the turbulent dissipation rate increases, causing the turbulent kinetic energy to be converted into heat energy. This is manifested as an overall decrease in turbulence intensity. The vortex structure sheds at both temperatures, and so the turbulence intensity decreases even more. The turbulence integral scale is the size of the vortex scale in the turbulent pulsation. According to Taylor's hypothesis (Piomelli et al. [56]), the axial turbulence integral scale can be expressed as: Figure 15a shows the two-dimensional turbulence integral scales inside the cavity. It can be seen that the larger turbulence integral scale is distributed around the model. The The turbulence integral scale is the size of the vortex scale in the turbulent pulsation. According to Taylor's hypothesis (Piomelli et al. [56]), the axial turbulence integral scale can be expressed as: Figure 15a shows the two-dimensional turbulence integral scales inside the cavity. It can be seen that the larger turbulence integral scale is distributed around the model. The turbulence integral scale is smaller inside the 0 • C than inside the 25 • C cavity, and the overall trend is gentler at the lower temperature. There are mainly small-scale vortices inside the 0 • C cavity, whereas those inside the 25 • C cavity are significantly larger. Figure 15b shows the three-dimensional turbulence integral scale inside the cavity. The 0 • C turbulence integral scale is slightly less than 0.02 m, and the 25 • C turbulence integral scale ranges from 0.02 to 0.03 m overall. This shows that, in the cavity, the integral scale of turbulence decreases at lower temperatures.  Figure 16a shows the two-and three-dimensional turbulence integral scales downstream of the cavity shedding point at the two temperatures. The shedding process moves the vortex structure downstream, and the turbulence integral scale increases significantly compared with the case inside the cavity. At 25 °C, the turbulence integral scale is greater than at 0 °C. This is because the overall scale of the vortex shedding is greater at 25 °C than at 0 °C, but the shedding frequency decreases, and the period becomes longer. Figure  16b shows the three-dimensional turbulence integral scale downstream of the cavity shedding point. The maximum value of the 0 °C turbulence integral scale ranges from 0.02 to 0.04 m, and the maximum value of the 25 °C turbulence integral scale is around 0.05 m. Larger turbulence integral scales are concentrated on both sides of the model tail. The overall vortex shedding scale is larger at 25 °C than at 0 °C, but the shedding frequency decreases, and the period becomes longer.  Figure 16a shows the two-and three-dimensional turbulence integral scales downstream of the cavity shedding point at the two temperatures. The shedding process moves the vortex structure downstream, and the turbulence integral scale increases significantly compared with the case inside the cavity. At 25 • C, the turbulence integral scale is greater than at 0 • C. This is because the overall scale of the vortex shedding is greater at 25 • C than at 0 • C, but the shedding frequency decreases, and the period becomes longer. Figure 16b shows the three-dimensional turbulence integral scale downstream of the cavity shedding point. The maximum value of the 0 • C turbulence integral scale ranges from 0.02 to 0.04 m, and the maximum value of the 25 • C turbulence integral scale is around 0.05 m. Larger turbulence integral scales are concentrated on both sides of the model tail. The overall vortex shedding scale is larger at 25 • C than at 0 • C, but the shedding frequency decreases, and the period becomes longer.

Conclusions
This paper has described how the DES method has been combined with the Schnerr-Sauer cavitation model to numerically simulate the unsteady cavitating flow of an experimental model at different temperatures. The effects of temperature on the cavity evolution, pressure fluctuations, vortex structure, and turbulence characteristics have been evaluated in detail. The main conclusions of this study are as follows: (1) The temperature affects the evolution of the cavity. Lower temperatures make the cavity unstable and increase the pressure pulsations inside the cavity, and also enhance the pressure pulsations where the cavity collapses. The cavity shedding frequency increases at lower temperatures, and the shedding cycle shortens. The cavity shedding size decreases, and the average volume of the cavity decreases. Additionally, the peak volume fluctuation increases, and the volume vibration amplitude becomes larger. This study also found that temperature changes have little effect on resistance.
(2) The temperature changes the vortex structure. At lower temperatures, during the cavity development stage, a larger number of linear vortex structures appear in the model wake area, and the vortex structure is significantly stronger than under

Conclusions
This paper has described how the DES method has been combined with the Schnerr-Sauer cavitation model to numerically simulate the unsteady cavitating flow of an experimental model at different temperatures. The effects of temperature on the cavity evolution, pressure fluctuations, vortex structure, and turbulence characteristics have been evaluated in detail. The main conclusions of this study are as follows: (1) The temperature affects the evolution of the cavity. Lower temperatures make the cavity unstable and increase the pressure pulsations inside the cavity, and also enhance the pressure pulsations where the cavity collapses. The cavity shedding frequency increases at lower temperatures, and the shedding cycle shortens. The cavity shedding size decreases, and the average volume of the cavity decreases. Additionally, the peak volume fluctuation increases, and the volume vibration amplitude becomes larger. This study also found that temperature changes have little effect on resistance. (2) The temperature changes the vortex structure. At lower temperatures, during the cavity development stage, a larger number of linear vortex structures appear in the model wake area, and the vortex structure is significantly stronger than under higher-temperature conditions. In the cavity stabilization stage, a lower temperature enlarges the vortex structure inside the cavity, thus accelerating the shedding of the cavity and shortening the cavity shedding cycle. Finally, in the cavity shedding stage, lower temperatures enhance the model wall vortex structure. In general, we found that the wake field vortex structure changes more drastically at the lower temperature. (3) The temperature affects the velocity pulsations around the model. At lower temperatures, the velocity pulsations are reduced, and this has a more significant impact on the wake flow field than the flow field inside the cavity. The lower temperature increases the turbulence dissipation rate and reduces the overall turbulence. At the same time, lower temperatures reduce the turbulence integral scale both outside and inside the cavity. The influence on the tail flow field is mainly concentrated on the two sides of the model tail, that is, the cavity moves downstream, and large-scale shedding becomes small-scale shedding.
In this paper, through numerical simulation of the axisymmetric body, the flow structure characteristics of the cavitation flow field are obtained. The focus is on the influence of temperature on the temporal and spatial distribution of natural cavitation flow.
Based on the axisymmetric body that is simulated, obtaining structural characteristics of flow cavitation flow field focuses on the effect of temperature on the spatial and temporal distribution of the natural air stream.
Although some research progress has been made, the modeling process of cavitation multiphase flow is still slightly inadequate. The sensitivity of grid independence is not compared, and the dependence of the number of generations at different time steps is not considered enough. Although some research progress has been made, the modeling process of cavitation multiphase flow is still slightly inadequate. The sensitivity of grid independence is not compared, and the dependence of the number of generations at different time steps is not considered enough.
Although research has made some progress, there are still some shortcomings in the process of multiphase flow modeling, there is no comparison for grid-independent sensitivity, and research is still lacking on the number of different time steps in generationdependent considerations.
In the follow-up research process, the grid sensitivity and the number of iterations will be analyzed more deeply.  Informed Consent Statement: Informed consent was obtained from all subjects involved in the study.
Data Availability Statement: Data sharing is not applicable to this article.

Conflicts of Interest:
The authors declare no conflict of interest.