# Test case n°30: Unsteady cavitation in a Venturi type section (PN)

Olivier Coutier-Delghosa, Regiane . Fortes Patella, Jean-Luc Reboud, Benoît Stutz

► To cite this version:

Olivier Coutier-Delghosa, Regiane . Fortes Patella, Jean-Luc Reboud, Benoît Stutz. Test case n°30: Unsteady cavitation in a Venturi type section (PN). Multiphase Science and Technology, 2005, 16 (1-3), pp.207-218. 10.1615/MultScienTechn.v16.i1-3.290. hal-00211304

HAL Id: hal-00211304

https://hal.science/hal-00211304v1

Submitted on 3 Jan 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
# Test case number 30: Unsteady cavitation in a Venturi type section (PN)

By

Olivier Coutier-Delgosha, ENSTA - UER de Mécanique, 91761 Palaiseau, France
Phone: +33 (0)1 69 31 98 18, Fax: +33 (0)1 69 31 99 97, E-Mail: coutier@enstay.ensta.fr

Regiane Fortes-Patella, LEGI/INPG, BP 53, 38041 Grenoble cedex 9, France
Phone: +33 (0)4 76 82 50 81, Fax: +33 (0)4 76 82 50 01, E-Mail: fortes@hmg.inpg.fr

Jean-Luc Reboud, LEMD / UJF, BP 166, 38042 Grenoble cedex 9, France
Phone: +33 (0)4 76 88 10 78, Fax: +33 (0)4 76 88 79 45,
E-Mail: jean-luc.reboud@grenoble.cnrs.fr

Benoît Stutz, CETHIL/INSA, 9 rue de la Physique, 69621 Villeurbanne cedex, France
Phone: +33 (0)4 72 43 82 51, Fax: +33 (0)4 72 43 72 10,
E-Mail: stutz@genserver.insa-lyon.fr

## 1 Practical significance and interest of the test-case

A calculation of the cavitating flow field occurring in a Venturi type duct is presented in the present paper. Experimental results obtained in the same flow configuration by Stutz & Reboud (2000) are also reported. Convergent and divergent angles of the lower wall of the Venturi type section are respectively about 18° and 8° (see figure 2). According to experimental observations in this geometry the flow is characterized by unsteady cavitation behavior, with quasi-periodic fluctuations. Each cycle is composed of the following successive steps: the attached sheet cavity grows from the Venturi throat. A re-entrant jet is generated at the cavity closure and flows along the Venturi bottom toward the throat. Its interaction with the cavity surface results in the cavity break off. The generated vapor cloud is then convected by the main stream, until it collapses.

The challenge consists in simulating correctly this unsteady behavior. Two tests are proposed to evaluate the consistency of the numerical solution with the experiments:

* the evaluation of overall parameters (mean volume of vapor, standard deviation, frequency of the periodic fluctuations, phase average of the cavity shape evolution),

* the description of the flow inside the sheet of cavitation (time-averaged values and standard deviations of velocities and void fraction)

The numerical simulation of this problem requires a coupling between the Navier-Stokes equations, a model of turbulence (the Reynolds number in the proposed configuration equals 1.6 10⁶), and a physical model of cavitation to predict the inception of cavitation and the behavior of the liquid/vapour mixture.

## 2 Definitions and physical model description

### 2.1 Physical model of cavitation

The present work considers a single fluid model: the fluid density $ρ$ varies in the computational domain according to a barotropic equation of state, $ρ(P)$, that links the density to
the local static pressure (see figure 1). This equation considers that phase change occurs within a small range of pressure, $\Delta P_{vap}$, centered on the saturation pressure $P_{vap}$. When the pressure in a cell is larger than $P_{vap} + \Delta P_{vap}/2$, the fluid is supposed to be pure liquid, the entire cell is occupied by liquid, and its density $\rho_l$ is calculated by the Tait equation (Knapp et al., 1970).

$$ \frac{\rho}{\rho_{ref}} = \sqrt[n]{\frac{P + P_0}{P_{ref}^T + P_0}}, \qquad (1) $$

where $P_{ref}^T$ is the pressure at the domain outlet, $\rho_{ref}$ is the liquid density, and for water, $P_0 = 3 \ 10^8$ Pa and the exponent is $n = 7$. If the pressure is lower than $P_{vap} - \Delta P_{vap}/2$, the cell is full of vapor and its density $\rho_v$ is given by the perfect gas law (isotherm approach),

$$ \frac{P}{\rho} = RT, \qquad (2) $$

where $R = 462$ J/K/kg for water vapor. In the other situations, the cell is occupied by a liquid/vapor mixture, which is considered as one single fluid with a variable density $\rho$. This one is directly related to the void fraction $\alpha = (\rho(P) - \rho_l)/(\rho_v - \rho_l)$ corresponding to the local density of the fluid.

To model the mixture state, the barotropic equation of state presents a smooth transition in the vapor pressure value neighborhood, in the range $P_{vap} \pm \Delta P_{vap}/2$. In direct relation with the range $\Delta P_{vap}$, this equation is characterized mainly by its maximum slope $1/C_{min}^2$, where $C_{min}^2 = \partial P/\partial \rho$. $C_{min}$ can thus be interpreted as the minimum speed of sound in the mixture. Its calibration was done in previous studies (Coutier-Delgosha et al., 2003b). The optimal value was found to be independent of the hydrodynamic conditions, and is about 1.5 m/s for cold water (20°C), with $P_{vap} = 0.023$ bar, and corresponding to $\Delta P_{vap} \approx 0.06$ bar. These values are used here throughout the presented results.

Mass fluxes resulting from vaporization and condensation processes are treated implicitly by the barotropic state law, and no supplementary assumptions are required. Concerning the momentum fluxes, the model assumes that locally velocities are the same for liquid and for vapor: in the mixture regions vapor structures are supposed to be perfectly carried by the main flow. This hypothesis is often assessed to simulate sheet-cavity flows, in which the interface is considered to be in dynamic equilibrium (Merkle et al., 1998). The momentum transfer between the phases is thus strongly linked to phase change.

## 2.2 Numerical resolution

To solve the time-dependent Reynolds-averaged Navier-Stokes equations associated with the barotropic equation of state presented here above, the numerical code applies, on 2D structured curvilinear-orthogonal meshes, the SIMPLE algorithm (Patankar, 1981)), modified to take into account the cavitation process. It uses an implicit method for the time-discretization, and the HLPA non-oscillatory second order convection scheme proposed by Zhu (1991). The numerical model is detailed in Coutier-Delgosha et al. (2003b).
<figure><img src="image_2.png" /><figcaption>**Figure 1:** Barotropic state law $\rho(P)$. Water 20°C.</figcaption></figure>

## 2.3 Turbulence model

In our previous studies (Coutier-Delgosha et al., 2002, 2003b), either the k-ω model proposed by Wilcox (1998) or the k-ε RNG model presented by Yakhot et al. (1992) were applied to model cavitating flows. Results obtained have demonstrated that for both models, corrections of the influence of the vapor/liquid mixture compressibility on the turbulence should be taken into account to obtain the unsteady effects due to cavitation. For the present test case, the modified k-ε RNG model presented in Coutier-Delgosha et al. (2003a) is applied.

# 3 Geometry and boundary conditions

The main features of the geometry can be seen in figure 2. The precise description of the Venturi section is given as a list of coordinates for the lower wall of the Venturi (see table 1) and the upper wall of the Venturi (see table 2).

The velocity field is imposed at the computational domain inlet, and the static pressure is imposed at the outlet. Along the solid boundaries, the turbulence models are associated with laws of the wall. Details of the prescribed values are given in section 3.3.

## 3.1 Grid

The computational grid is composed of 160×50 orthogonal cells (figure 2). A special contraction of the mesh is applied in the main flow direction just after the throat, so that the two-phase flow area is efficiently simulated: about fifty grid points are used in this direction to model the 45 mm long mean cavity obtained by numerical calculations in the case $\sigma = 2.4$ where $\sigma$ is defined in (3) (results presented hereafter). In the other direction, a contraction is also applied close to the walls, to obtain at the first grid point the non-dimensional parameter $y^+$ of the boundary layer varying between 30 and 100, and to use standard laws of the wall. The grid is finer in the bottom part of the Venturi section than in its upper part, to enhance the accuracy in the cavitation domain: cavities obtained contain about thirty cells across their thickness.
<table><caption>**Table 1:** Description of the lower wall of the Venturi shown in figure 2. Coordinates must be read from top to bottom and left to right.</caption><thead><tr><th>x<br>(mm)</th><th>y<br>(mm)</th><th>x<br>(mm)</th><th>y<br>(mm)</th><th>x<br>(mm)</th><th>y<br>(mm)</th><th>x<br>(mm)</th><th>y<br>(mm)</th></tr></thead><tbody><tr><td>-152.093</td><td>0.000</td><td>13.361</td><td>14.980</td><td>61.921</td><td>8.665</td><td>151.062</td><td>-6.654</td></tr><tr><td>-51.706</td><td>0.000</td><td>14.508</td><td>14.833</td><td>64.027</td><td>8.376</td><td>154.938</td><td>-7.351</td></tr><tr><td>-42.924</td><td>2.163</td><td>15.676</td><td>14.683</td><td>66.178</td><td>8.079</td><td>158.903</td><td>-8.064</td></tr><tr><td>-35.553</td><td>4.558</td><td>16.869</td><td>14.531</td><td>68.396</td><td>7.770</td><td>162.936</td><td>-8.794</td></tr><tr><td>-29.280</td><td>6.594</td><td>18.085</td><td>14.377</td><td>70.636</td><td>7.450</td><td>167.058</td><td>-9.547</td></tr><tr><td>-23.971</td><td>8.324</td><td>19.326</td><td>14.218</td><td>72.921</td><td>7.117</td><td>171.269</td><td>-10.323</td></tr><tr><td>-19.437</td><td>9.795</td><td>20.591</td><td>14.054</td><td>75.251</td><td>6.771</td><td>175.548</td><td>-11.128</td></tr><tr><td>-15.589</td><td>11.045</td><td>21.881</td><td>13.889</td><td>77.626</td><td>6.411</td><td>179.894</td><td>-11.965</td></tr><tr><td>-12.318</td><td>12.108</td><td>23.187</td><td>13.719</td><td>80.045</td><td>6.037</td><td>184.353</td><td>-12.837</td></tr><tr><td>-9.538</td><td>13.011</td><td>24.531</td><td>13.540</td><td>82.510</td><td>5.647</td><td>188.878</td><td>-13.749</td></tr><tr><td>-7.175</td><td>13.779</td><td>25.898</td><td>13.357</td><td>85.041</td><td>5.240</td><td>193.493</td><td>-14.705</td></tr><tr><td>-5.166</td><td>14.432</td><td>27.309</td><td>13.180</td><td>87.595</td><td>4.818</td><td>198.198</td><td>-15.708</td></tr><tr><td>-3.459</td><td>14.986</td><td>28.720</td><td>12.998</td><td>90.216</td><td>4.378</td><td>202.992</td><td>-16.758</td></tr><tr><td>-2.007</td><td>15.458</td><td>30.177</td><td>12.805</td><td>92.882</td><td>3.923</td><td>207.876</td><td>-17.859</td></tr><tr><td>-0.774</td><td>15.859</td><td>31.655</td><td>12.613</td><td>95.615</td><td>3.450</td><td>212.827</td><td>-19.004</td></tr><tr><td>0.274</td><td>16.201</td><td>33.179</td><td>12.419</td><td>98.393</td><td>2.961</td><td>217.912</td><td>-20.188</td></tr><tr><td>0.509</td><td>16.298</td><td>34.724</td><td>12.220</td><td>101.216</td><td>2.456</td><td>223.065</td><td>-21.396</td></tr><tr><td>0.751</td><td>16.377</td><td>36.293</td><td>12.012</td><td>104.106</td><td>1.934</td><td>240.606</td><td>-24.016</td></tr><tr><td>0.992</td><td>16.455</td><td>37.906</td><td>11.799</td><td>107.041</td><td>1.398</td><td>258.752</td><td>-26.682</td></tr><tr><td>1.239</td><td>16.509</td><td>39.541</td><td>11.589</td><td>110.043</td><td>0.848</td><td>277.123</td><td>-29.392</td></tr><tr><td>2.143</td><td>16.395</td><td>41.199</td><td>11.379</td><td>113.089</td><td>0.283</td><td>295.717</td><td>-32.170</td></tr><tr><td>3.064</td><td>16.273</td><td>42.901</td><td>11.160</td><td>116.226</td><td>-0.294</td><td>314.983</td><td>-35.016</td></tr><tr><td>4.005</td><td>16.161</td><td>44.649</td><td>10.932</td><td>119.407</td><td>-0.881</td><td>334.474</td><td>-37.906</td></tr><tr><td>4.964</td><td>16.040</td><td>46.419</td><td>10.699</td><td>122.633</td><td>-1.479</td><td>354.188</td><td>-40.840</td></tr><tr><td>5.943</td><td>15.915</td><td>48.211</td><td>10.463</td><td>125.949</td><td>-2.085</td><td>374.575</td><td>-43.842</td></tr><tr><td>6.941</td><td>15.794</td><td>50.070</td><td>10.225</td><td>129.331</td><td>-2.699</td><td>395.185</td><td>-46.911</td></tr><tr><td>7.959</td><td>15.664</td><td>51.952</td><td>9.984</td><td>132.781</td><td>-3.328</td><td>416.244</td><td>-50.048</td></tr><tr><td>8.997</td><td>15.529</td><td>53.856</td><td>9.736</td><td>136.299</td><td>-3.975</td><td>437.751</td><td>-52.938</td></tr><tr><td>10.056</td><td>15.397</td><td>55.805</td><td>9.480</td><td>139.883</td><td>-4.632</td><td>459.929</td><td>-55.111</td></tr><tr><td>11.136</td><td>15.263</td><td>57.799</td><td>9.215</td><td>143.535</td><td>-5.296</td><td>1111.627</td><td>-119.228</td></tr><tr><td>12.238</td><td>15.123</td><td>59.838</td><td>8.945</td><td>147.254</td><td>-5.970</td><td></td><td></td></tr></tbody></table>

## 3.2 Initial conditions

To start unsteady calculations, the following numerical procedure is applied: first of all, a stationary step is carried out, with an outlet pressure value high enough to avoid any vapor in the whole computational domain. Then, this pressure is lowered slowly at each new time-step, down to the value corresponding to the desired cavitation number $\sigma$ defined by

$$ \sigma = \frac{P_{upstream} - P_{vap}}{\rho_{ref} V_{ref}^2 / 2}. \qquad (3) $$

Vapor appears during the pressure decrease. The cavitation number is then kept constant throughout the computation.
<table><caption>**Table 2:** Description of the upper wall of the Venturi shown in figure 2. Coordinates must be read from top to bottom and left to right.</caption><thead><tr><th>x<br>(mm)</th><th>y<br>(mm)</th><th>x<br>(mm)</th><th>y<br>(mm)</th></tr></thead><tbody><tr><td>-152.093</td><td>50.003</td><td>194.008</td><td>46.598</td></tr><tr><td>116.114</td><td>50.003</td><td>198.735</td><td>46.329</td></tr><tr><td>119.340</td><td>49.958</td><td>203.552</td><td>46.083</td></tr><tr><td>122.611</td><td>49.936</td><td>208.480</td><td>45.814</td></tr><tr><td>125.971</td><td>49.869</td><td>213.499</td><td>45.522</td></tr><tr><td>129.399</td><td>49.801</td><td>218.629</td><td>45.254</td></tr><tr><td>132.871</td><td>49.712</td><td>223.871</td><td>44.962</td></tr><tr><td>136.433</td><td>49.600</td><td>229.181</td><td>44.671</td></tr><tr><td>140.062</td><td>49.465</td><td>246.879</td><td>43.708</td></tr><tr><td>143.759</td><td>49.309</td><td>265.025</td><td>42.700</td></tr><tr><td>147.522</td><td>49.129</td><td>283.619</td><td>41.692</td></tr><tr><td>151.376</td><td>48.928</td><td>302.438</td><td>40.661</td></tr><tr><td>155.274</td><td>48.726</td><td>321.704</td><td>39.586</td></tr><tr><td>159.262</td><td>48.502</td><td>341.195</td><td>38.510</td></tr><tr><td>163.316</td><td>48.278</td><td>361.133</td><td>37.413</td></tr><tr><td>167.461</td><td>48.054</td><td>381.744</td><td>36.293</td></tr><tr><td>171.673</td><td>47.830</td><td>402.354</td><td>35.150</td></tr><tr><td>175.974</td><td>47.584</td><td>423.637</td><td>33.985</td></tr><tr><td>180.343</td><td>47.337</td><td>445.368</td><td>32.775</td></tr><tr><td>184.823</td><td>47.113</td><td>1120.140</td><td>-9.653</td></tr><tr><td>189.371</td><td>46.844</td><td></td><td></td></tr></tbody></table>

## 3.3 Calculations

Calculations are performed with non-dimensional variables based on the following reference parameters :

<table><tr><td>U<sub>a</sub> = U/U<sub>ref</sub></td><td>U<sub>ref</sub> = 7.2 m/s at the inlet</td></tr><tr><td>ρ<sub>a</sub> = ρ/ρ<sub>ref</sub></td><td>ρ<sub>ref</sub> = ρ<sub>liquid</sub></td></tr><tr><td>C<sub>p</sub> = (P - P<sub>ref</sub>) / (1/2ρ<sub>ref</sub>U<sub>ref</sub><sup>2</sup>)</td><td>P<sub>ref</sub> = P<sub>outlet of the domain</sub></td></tr><tr><td>σ = (P<sub>ref</sub> - P<sub>vap</sub>) / (1/2ρ<sub>ref</sub>U<sub>ref</sub><sup>2</sup>)</td><td>P<sub>vap</sub> = 2000 Pa</td></tr><tr><td>T<sub>a</sub> = T/T<sub>ref</sub></td><td>T<sub>ref</sub> = L<sub>ref</sub>/U<sub>ref</sub></td></tr><tr><td></td><td>L<sub>ref</sub> = chord of the Venturi = 0.224 m</td></tr></table>

Physical and numerical parameters applied:

<table><tr><td>C<sub>min</sub> (minimum speed of sound in the mixture)</td><td>1.5 m/s</td></tr><tr><td>Re (Reynolds number based on L<sub>ref</sub>, V<sub>ref</sub>, and water properties)</td><td>1.6 10<sup>6</sup></td></tr><tr><td>ρ<sub>v</sub>/ρ<sub>l</sub> (ratio of vapor to liquid density)</td><td>0.01</td></tr><tr><td>Inlet turbulence level</td><td>1%</td></tr><tr><td>Δt (non dimensional time-step)</td><td>0.005</td></tr><tr><td>time order discretization</td><td>1</td></tr></table>

It is proposed to consider six flow configurations corresponding to the six cavitation numbers, $\sigma$, reported in table 3. For each calculation, the total simulation time required to eliminate the initial transient effects is at least 40 $T_{ref}$. The sheet of cavitation systematically adopts the oscillating behavior observed in experiments. The oscillations
<figure><img src="image_3.png" /><figcaption>**Figure 2:** Curvilinear-orthogonal mesh of the Venturi type section.</figcaption></figure>

are almost periodic, and their frequency as well as the mean cavity length depends on the cavitation parameter, $\sigma$.

As the final cavity obtained is fundamentally unstable, it cannot be characterized by its final shape or the final void fraction distribution. The comparisons are thus based on the transient evolution of the cavitating flow. This evolution can be defined at each time by the vapor quantity present in the domain or by the cavity shape (length, volume). We propose to focus on the vapor volume oscillations. For each computation, the time-averaged vapor volume and its standard deviation are indicated in table 3. The cavitation cycle frequency can be calculated by using a FFT analysis of the inlet pressure signal and is also given in table 3.

<table><caption>**Table 3:** Values of the six cavitation numbers with the corresponding results.</caption><thead><tr><th colspan="4">RESULTS</th></tr><tr><th>σ</th><th>Mean vapor Volume /<br>L<sup>3</sup><sub>ref</sub> × 10<sup>4</sup></th><th>Standard deviation of the<br>vapor volume / L<sup>3</sup><sub>ref</sub> × 10<sup>4</sup></th><th>Oscillation<br>frequency (Hz)</th></tr></thead><tbody><tr><td>2.32</td><td>123</td><td>32.9</td><td>17</td></tr><tr><td>2.34</td><td>44.3</td><td>9.1</td><td>30</td></tr><tr><td>2.37</td><td>19.7</td><td>5.8</td><td>41</td></tr><tr><td>2.40</td><td>13.8</td><td>2.4</td><td>55</td></tr><tr><td>2.44</td><td>6.9</td><td>1.9</td><td>68</td></tr><tr><td>2.52</td><td>4.4</td><td>1.6</td><td>82</td></tr></tbody></table>

The transient evolution observed for $\sigma$=2.4 during the unsteady calculation is presented in figure 3. Figure 3(a) illustrates at a given time and for each cross section of the Venturi type duct the value of the minimal density present in the section. It gives information concerning the vapor cloud shedding process: the part of the cavity that breaks off clearly appears, and the fluctuation frequency can be easily calculated.
Moreover, it also shows the minimum density, i.e. the maximum void ratio in each section. The two other curves shown in figures 3(b) and 3(c) represent respectively the total vapor volume and the inlet pressure evolutions.

# 4 Comparison with experiments

## 4.1 Overall behavior

First, the evolution of the cavity shape at a given cavitation number is compared with pictures obtained from experiments. Two phase-averaged cavitation cycles are presented in figure 4. The right one results from experimental visualizations: Video frames acquired during a 100 ns exposure time under Laser sheet light are identified and digitized in 256 grey levels. A sampling technique is applied to classify them in nine sets corresponding to the different states of the recorded quasi-periodic pressure signal. Then, averaging the grey levels pixel per pixel for each set allows drawing a sequence of phase-averaged images, from an initial data set of 300 frames. The left part of the figure 4 corresponds to the same sequence obtained by numerical simulation (calculation duration equal to 20 $T_{ref}$, i.e. about 30 cycles). The same sampling technique is applied: the computational result is decomposed into 30×9 short sequences corresponding to the nine steps of the cavitation cycle and the phase-averaging process is applied.

From the results of table 3 concerning the effect of the cavitation number, a comparison can be proposed with results reported by Stutz & Reboud (2000): the frequency of the self-oscillation behavior is drawn with respect to the ratio $V_{ref}/L_{cav}$. In both the numerical simulation and the experiments $L_{cav}$ is chosen as the maximum length of the attached cavity.

## 4.2 Flow field inside the sheet of cavitation

Local comparisons are proposed in the case $\sigma = 2.4$ with experimental data obtained by double optical probes measurements. This technique and the results are presented in detail in Stutz & Reboud (1997, 2000). This is an intrusive technique, which allows measurements of the local void ratio and the velocities of the two-phase structures inside the cavitation sheet. Four data profiles located respectively at $x = 1.410^{-2}$ m, $3.110^{-2}$ m, $4.910^{-2}$ m, and $6.510^{-2}$ m, are available. The time-averaged and standard deviation values of the velocity $u$ and the void ratio $\alpha$ are compared along the four profiles in figure 6.

The three main features that should be obtained are:

* presence of the re-entrant jet, characterized by negative or zero mean values of the velocity $u$ close to the wall.

* the rather low mean void ratio observed in the main part of the cavitation sheet: it does not exceed 25%, excepted in the upstream end of the cavity.

* the general high level of velocity and void fraction fluctuations of the same order of magnitude than the mean values.

These three characteristics of the flow are strongly related to the overall unsteady oscillation of the cavitation sheet. So an accurate numerical simulation of these features is linked to the capability of the model to predict the unsteadiness of the cavitating flow.
<figure><img src="image_4.png" /><figcaption>(a) Temporal evolution (in abscissa) of the cavity length (graduated in ordinate) (Instantaneous attached and cloud cavities at $T = 12T_{ref}$ are given at left)</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>(b) Time evolution of the volume of vapor in the flow field</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>(c) Time evolution of the inlet pressure, $(P_{inlet} - P_{outlet})/\frac{1}{2}\rho lU_{ref}^2$</figcaption></figure>

**Figure 3:** Transient evolution of the unsteady cavitating flow in the Venturi type duct ($\sigma = 2.4$).
<figure><img src="image_7.png" /><figcaption>**Figure 4:** Numerical (on the left) and experimental (on the right) phase-averaged sequences of unsteady cavitation ($\sigma = 2.4$).</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 5:** Experimental and numerical oscillation frequency.</figcaption></figure>
<figure><img src="image_9.png" /><figcaption>**Figure 6:** Time-averaged and standard deviation values of void ratio $\alpha$ and velocity $u$. Comparison between numerical results (lines) and optical probes measurements (points) for $\sigma = 2.4$. Cavity external shape in dotted line from image processing - ratio 3 between vertical and horizontal scales.</figcaption></figure>

## Acknowledgements

The experimental and numerical works presented in this test-case have been obtained with the continuous support of the French Space Agency (CNES). The physical model and the numerical results presented for this test case have been published by Coutier-Delgosha *et al.* (2002, 2003a) and Coutier-Delgosha *et al.* (2003b).

## References

Coutier-Delgosha, O., Fortes-Patella, R., & Reboud, J. L. 2002. Simulation of unsteady cavitation with a two-equations turbulence model including compressibility effects. *J. of Turbulence*, **3**, 058. jot.iop.org.

Coutier-Delgosha, O., Fortes-Patella, R., & Reboud, J. L. 2003a. Evaluation of the turbulence model influence on the numerical simulations of unsteady cavitation. *J of Fluids Eng.*, **125**, 38-45.

Coutier-Delgosha, O., Fortes-Patella, R., & Delannoy, Y. 2003b. Numerical simulation of unsteady cavitating flow. *Int. J. for Numerical Methods in Fluids*, **42**(5), 527-548.

Knapp, R. T., Daily, J. T., & Hammit, F. G. 1970. *Cavitation*. Mc Graw Hill.

Merkle, C. L., Feng, J., & Buelow, P. E. O. 1998. Computational modeling of the dynamics of sheet cavitation. *Pages 307-314 of: Proc. of the 3rd Int. Symp. on Cavitation, Grenoble, France*, vol. 2.

Patankar, S. V. 1981. *Numerical heat transfer and fluid flow*. Hemisphere Publishing Corporation.

Stutz, B., & Reboud, J. L. 1997. Experiments on Unsteady Cavitation. *Experiments in Fluids*, **23**, 191-198.
Stutz, B., & Reboud, J.-L. 2000. Measurements within unsteady cavitation. *Experiments in Fluids*, **29**, 545–552.

Wilcox, D. 1998. *Turbulence modeling for CFD*. DCW Industries, Inc., La Canada, California, USA.

Yakhot, V., Orszag, S. A., Thangham, S., Gatski, T. B., & Speziale, C. G. 1992. Development of turbulence models for shear flows by a double expansion technique. *Phys. Fluids A*, **4**, 1510–1520.

Zhu, J. 1991. A low diffusive and oscillation-free convection scheme. *Comm. in Applied Num. Methods*, **7**, 225–232.