Abstract

The integral-differential equation was obtained to simulate unsteady evolutions of the slender axisymmetric ventilated supercavity with the use of one-dimensional inviscid flow of the incompressible gas in the channel between the cavity surface and the body of revolution. For small ventilation rates, the solution of this equation was expressed as asymptotic series. In the steady case the nonlinear differential equation and its solutions were obtained. It was shown that the ventilation can increase and diminish the supercavity dimensions. Examples of calculations for different hull shapes are presented. At some critical values of the gas injection rate the cavity dimensions can become unbounded. Stability of steady and pulsating gas cavities was investigated in the case of the low gas injection rate.

1. Introduction

The drag of high speed underwater vehicles can be reduced with the use of supercavitation. This idea was developed in many theoretical, numerical, and experimental investigations in a lot of countries. In 2013 we celebrated the 100 anniversary of the famous Ukrainian academician G. V. Logvinovich who sufficiently contributed to both the theoretical and experimental research of supercavitation and practical applications of this phenomenon. His principle of independence (Logvinovich [1]) is still the very efficient tool for calculating the shape of long 3D cavities. In this paper we will discuss its accuracy and areas of application.

To obtain small cavitation numbers at small vehicle velocities or at large movement depths, a gas ventilation inside the cavity is used (see, e.g., Logvinovich [1], Epshtein [24], Knapp et al. [5], Yegorov et al. [6], Levkovsky [7], Kuklinski et al. [8], Spurk [9], Wosnik et al. [10], Zhuravlev and Varyukhin [11], Matveev et al. [12], and Vlasenko and Savchenko [13]). In most cases the gas and vapor densities inside the cavity are much smaller than the density of water (approximately 800 times). This fact allows neglecting the influence of the gas flow inside the cavity and assuming the constant pressure on its surface. As a result, in the case of small gravity effects, the cavitation number ( is the water density, is the velocity of movement, is the water vapor pressure at ambient temperature, and and are pressures measured in the cross section of the cavity origin far away in the water flow and in the injected gas, resp.) is supposed to be constant over all of the cavity surface and the ventilated cavity shape is assumed to be coinciding with the vapor one. This fact is usually used both in theories of the ventilated cavities (e.g., Logvinovich [1], Epshtein [24], and Spurk [9]) and in experiments. For example, in tests by Wosnik et al. [10] and Vlasenko and Savchenko [13] the pressure inside the ventilated cavity was not measured, but the shape of the corresponding vapor cavity was used to estimate the cavitation number. Even in this case, when the effects of the gas flow inside the cavity and gravity are negligible, there is no complete theory for the cavity shape as a function of the gas supply rate, cavitation number, and the shape of the body located inside the cavity.

To improve the efficiency of the supercavitation in comparison with the attached flow pattern, the cavitation number has to be reduced and the cavity volume has to be used completely to locate the vehicle hull (see, e.g., Nesteruk [14]). It means that an injected gas must flow in a narrow channel between the cavity surface and the vehicle hull; see Figure 1. If the ventilation rate is large enough, the pressure on the cavity surface is no more constant and changes the cavity shape in comparison with the case of vapor cavitation. This complicated phenomenon could be investigated numerically with the use of viscous fluid equations (e.g., Zhuravlev and Varyukhin [11]). The ideal fluid approach and the slender body theory allow one to obtain simple equations for the shape of axisymmetric ventilated supercavity provided that the gas flow between the cavity surface and the body of revolution is one-dimensional inviscid and incompressible. Some interesting results were obtained in Manova et al. [15] and Nesteruk and Shepetyuk [16, 17] for steady flow of liquid without gravity effects.

Here the results of these papers are generalized for the unsteady vertical flows in the gravity field. The integral-differential equation is obtained to simulate unsteady evolutions of the slender axisymmetric ventilated supercavity. For small ventilation rates, a solution is expressed as asymptotic series. Examples of the steady cavity shape calculations for different hulls are presented. Peculiarities of the cavity closing near the hull contour discontinuities were investigated.

Stability and pulsations of ventilated cavities were investigated experimentally in Silberman and Song [18], Song [19], Michel [20], and Zou et al. [21] and by many other authors. The theories of such empty cavities were proposed by Woods [22], Paryshev [23, 24], Nesteruk [25], and Semenenko [26]. Here we will try to construct simple stability theories for the steady and pulsating cavities with the use of the main term of the developed asymptotic solution. In comparison with the previous investigations we will take into account cases of the hull, located inside the cavity, different cavitator shapes, and the influence of the gravity.

2. Slender Body of Revolution with a Prescribed Pressure Distribution

The theory of supercavitation is very complicated due to the absence of exact solutions for 3D unsteady flows even in ideal fluid. But for long cavities similar to the slender bodies of revolution, the corresponding linear theories can be used. In monograph (Cole [27]) the potential of the steady flow of ideal incompressible liquid was obtained by the method of matched asymptotic expansions. In Nesteruk [28], this slender body theory was generalized for the unsteady case. The potential of the external water flow can be written as follows: We use the symbols with “prime” for physical (dimensional) parameters; is the typical time of unsteadiness, for example, the pulsation period; is the body or the cavity radius; is a small parameter, a ratio of the maximum radius of the body (or a cavitator-cavity system) to its length at a fixed moment of time ; are cylindrical coordinates shown in Figure 1; the integration in (5) must be carried out along the axis of symmetry from the nose till the tail of the body. If the body shape and the function are known, (2)–(5) allow one to calculate the flow field and the pressure distribution.

In supercavitating flows (or when the body shape is unknown), the radius may be calculated with the use of pressure distribution over the surface. Putting expression (2) into the Cauchy-Lagrange integral yields (Nesteruk [28, 29]) Here is the Froude number, is the gravity acceleration, and is pressure in liquid far from the body at the level . The sign “−” in the last term of (6) corresponds to the gravitational acceleration vector directed along the axis and the sign “+” corresponds to the case, when the directions of the body movement and gravity coincide (in order to have an axisymmetric cavity the direction of the gravity force is limited by these two cases). The dimensionless speeds are based on the current flow velocity .

Nonlinear integral-differential equation (6) relates the slender body (or the cavity) radius and the pressure distribution in liquid over its surface . If the functions , , and and the cavitator shape are given, the set of (3)–(8) allows one to calculate the unknown cavity radius . The asymptotic solution of the nonlinear integral-differential equation (6) was obtained in Nesteruk [30].

3. Equation of the First Approximation and Logvinovich’s Principle of Independence

The leading term of (6) (of order ) provides the equation of the first approximation (Nesteruk [28]): Thus the squared body (or cavity) radius can be found by successive integration of two linear differential equations (9) and (3) with partial derivatives of the first order. The corresponding general solution was obtained in Nesteruk [28].

Equation (9) can be rewritten in the absolute coordinates , (in which the liquid at infinity does not move and the body moves vertically). Let us suppose that the origin of this coordinate system is located on the undisturbed free water surface; the directions of the axis of symmetry and the cavitator movement are opposite and the pressure on the cavity surface depends on the time only. Then (9) attains the following form (Nesteruk [30]): where is some fixed moment of time and the atmospheric pressure on the undisturbed free water surface is assumed to be constant. The thickness parameter can be assumed to be coinciding with the ratio of the cavitator maximum radius to its length, for example, for slender conical cavitator ( is the cone angle).

According to (10), the evolution of every cross section of the cavity (at fixed value of ) is independent of the cavity behavior at other cross sections. This fact can be considered as one more theoretical support of Logvinovich’s principle of independence. An equation similar to (10) with an empirical constant instead of is used in the computer codes developed in the Institute of Hydromechanics NASU for the calculations of unsteady supercavities, for example, Semenenko and Naumova [31].

It must be noted that the accuracy of the first approximation (9) or (10) is limited by value , which yields only 44% for and 22% for . However these equations provide a nice qualitative description and were used to estimate the influence of many physical parameters such as gravity, capillarity, liquid compressibility, the cavitator shape, motion unsteadiness, and pressure gradients over the body surface; see Nesteruk [2830, 3234]. To improve the accuracy, the integral-differential equation (6) must be used, in which the independence of the cavity cross section is no more valid due to the integral term (5).

4. Shape of the Ventilated Axisymmetric Slender Cavity

In the case of ventilated cavity the total pressure of the water vapor and the gas must be equal to the pressure in water on the cavity surface . Here we neglect the capillarity forces and use the approach of one-dimensional flow of an ideal incompressible gas in the circular channel between the cavity surface and the hull located inside the cavity (see Figure 1). Then the continuity equation yields that the gas flux inside the cavity is only time dependent. Here the dimensionless gas velocity is based on ; the dimensionless cavity and body radii ( and ) are based on . By neglecting the longitudinal curvature of the circular channel (the body inside the cavity is assumed to be slender too), the potential of the one-dimensional gas flow can be determined from the relationship and (12) as follows: where is an arbitrary time dependent function.

The Cauchy-Lagrange integral for the gas flow has the following form: where is the gas density. Gravity and inertial forces are neglected in (15), since and they are very much smaller in comparison with the same forces in water flow, taken into account in (6). Putting expressions (12) and (14) into (15) yields Here and are the radii of the cavity and the body at the cross section of cavity origin (see Figure 1); the cavitator and body shapes can be changeable; is given by (1). Substitution of (16) into (6) yields the following nonlinear integral-differential equation for the ventilated cavity shape:

Even limited by the leading term of order , (17) remains integral-differential and the principle of independence is no more valid for unsteady ventilated cavities. In the case of a large velocity of the ventilated gas () the last term in (17) is proportional to and is much greater than two previous terms proportional to (see (12)). Thus, the shape of such slender intensively ventilated cavities can be estimated with the use of the following equation: which does not contain any integral terms and corresponds to the Logvinovich independence principle.

To solve both (17) and (18), the standard initial and boundary conditions in two domains can be used (Nesteruk [28]): where , , , and are given functions.

5. Steady Ventilated Supercavities

In steady case (17) limited by the leading term of order attains the following form: where all lengths are based on the cavity radius at its origin (see Figure 1). Relation (20) follows also from (18), but in the steady case (20) is valid at arbitrary values of ventilation rate. To characterize the influence of gas injection, other nondimensional parameters can be used instead of ; for example, where is the velocity of gas at . The dimensionless ventilation rate is the ratio of the pressure heads in the gas and water flows at . If the gravity forces can be neglected also in liquid flow (), nonlinear differential equation (20) coincides with the one obtained in Manova et al. [15].

With the use of parameter (20) can be rewritten as follows: The standard initial conditions at can be used to integrate (20) or (24), where is the derivative of the body radius at the cavity origin. It can be seen from (20) and (22) that the ventilation can strongly change the shape of the supercavity only at small values of or at the large ratio , since is very small.

Without ventilation () equation (20) or (23) coincides with the one proposed in Nesteruk [32]. Its solution in this case can be easily obtained with the use of initial conditions (26): For ventilated supercavities the nonlinear differential equation (20) or (23) must be solved. It can be seen that gas injection strongly changes the cavity shape not only at great values of the parameter but also at very small values of the circular channel width .

The term in brackets in (20) equals zero at and can be both positive and negative at different cavitator and hull shapes. It means that the corresponding ventilated cavity can, respectively, be both larger and smaller than the vapor one (27). Examples of calculations can be found in Manova et al. [15] and Nesteruk and Shepetyuk [16] for the first case and in Nesteruk and Shepetyuk [17] for the second one.

When the gravity forces can be neglected, (20) or (23) demonstrates that two independent parameters which are the pressure inside the cavity () and the gas injection rate ( or or ) influence the cavity shape. The gas leakage from the stable artificial cavity is equal to the gas injection rate and depends on the cavity dimensions, cavity closing behaviour, and so forth. Therefore, the gas leakage can also be different at the fixed cavitation number (as the equal value of the gas injection rate). In particular, it is written in Logvinovich [1] that “it is impossible to suggest a universal method for the theoretical determination of the gas loss.” This conclusion was supported by experiments (see, e.g., Wosnik et al. [10] and Vlasenko and Savchenko [13]), where different gas leakage/injection rates were obtained at the fixed cavitation number. Nevertheless, there are some attempts to present the gas leakage from the empty cavity () at high Froude numbers as a function of cavitation number (see, e.g., Logvinovich [1] and Spurk [9]). At small Froude numbers in horizontal flow (the cavity has no more the axis of symmetry) the formulas representing the gas leakage rate as function of the cavitation and Froude numbers were proposed in Epshtein [24].

6. Different Cavitator Shapes and Restrictions of the Flow Parameters

The developed theory and equations can be applied for a different slender cavitator—a part of a hull wetted by water (). For example, the slope of the cavitator radius at the cavity origin section can be both positive (; see Figure 1, a slender conical cavitator) and nonpositive (; see Figure 2, a case of base cavities). In many applications the disc or nonslender conical cavitators are used to create a long supercavity. Equations (17), (20), and (23) can be used in the case of the nonslender cavitator too, provided that the initial part of the long cavity is known. In particular, in the steady case the Logvinovich empirical formula (Logvinovich [1]) is valid. In (29) all lengths are based on the radius of the disc cavitator. In particular, from (29) it follows that at the corresponding cavity radius and its slope are 3 and , respectively. If at the cavity is slender enough, we can use (20) with all lengths based on . In the initial condition (26) must be replaced by . It is difficult to expect a good accuracy from such calculations, but an estimation of the supercavity shape is possible. In particular, at , the vapor cavity maximum radius can be determined from formula (27) as follows: Calculations of the maximum radius of the supercavity created by the disc cavitator with the use of the Logvinovich empirical formula (Logvinovich [1]) differ from (30) less than 10% at .

It is important to know which values of parameters , and are possible in real flows. To demonstrate the influence of the cavitator shape let us consider a vapor cavity in a liquid without gravity. The dependences (27) corresponding to are presented in Figure 3 and attached to the case in Figure 4. According to the stability principle, a small change in the parameters determining a solution ought to give rise to small changes in this solution.

At first let us consider the case . Since the parameter is negative at the positive cavitation numbers, then the branches of parabola (27) are directed to the bottom; accordingly at they are directed upward. For zero cavity number the parabola (27) degenerates into linear dependence (see Figure 3). Such character of solution testifies about its steadiness at . But already at zero cavitation number the stability disappears, because at formula (27) gives an infinite cavity, and at any positive cavitation number (which can be very near to zero) the cavity dimensions are bounded. The same character of instability takes place for . So, according to the stability principle only the flows with the positive cavitation numbers can be realized for cavitators with . The solution at for such cavitators is unstable and has the theoretical meaning only.

For the direction of parabola (27) branches does not change and it similarly degenerates into the linear dependence at zero cavitation number (see Figure 4). But due to the negative value of parameter , the zero cavitation number corresponds to a bounded cavity. The bounded cavities exist also for some negative cavitation numbers (small enough for module (see Figure 4)). Critical value of corresponds to the case when the parabola (27) touches the axis . At cavitation numbers, lesser than , cavity is unbounded. Application of the stability principle yields that at only flows with can be realized. At , the minimum cavitation number may be defined from (27) as follows: (Nesteruk [32]).

In the case of liquid with gravity the limitations of values of parameters follow from (27). Specifically, for the steady flow directed to the bottom, the coefficient of the polynomial (27) is negative; therefore a cavity always has the bounded dimensions, and steadiness principle is valid for arbitrary values of . For the opposite direction of the flow, this coefficient changes its sign, and the situation is similar to the case shown in Figure 4. The dependences of the minimum cavitation number versus the Froude number are calculated in Nesteruk [32].

The stability principle can be used to analyze the parameter restrictions for unsteady supercavity flows too. An example of a slender axisymmetric cavitator located in the steady flow of liquid with gravity directed upwards with the cavitation number, which varies according to the linear law , is presented in Nesteruk [35]. In this flow the pressure inside the cavity is only time dependent and we can use (6) with or (17) with . The solution of the equation of the first approximation (9) in this case can be written as follows (see Nesteruk [35]):

If the parameter satisfies the condition , then a moment of time comes, when the cavity becomes unbounded for . The critical cavitation number can be also calculated. For example, for the formula can be obtained from (32). Thus, the diminishing of the cavitation number to the value lesser than critical is impossible. Probably, the supercavity flow pattern is lost in the moment when .

It is interesting to emphasize that such flow of a liquid without gravity is always stable. Really, it follows from (32) that, when the cavitation number diminishes (), the cavity dimensions are finite both at the negative and at positive current values of the cavitation number . According to (32) the cavity dimensions diminish continuously when the cavitation number increases. In particular, there are possible unsteady supercavity flows with a conical cavitator at negative cavity numbers, while steady flows are impossible at .

In this example we considered the evolutions of the cavity shape forced by the changes of the pressure inside the cavity only (or by the variations of the cavitation number ). Solution (32) is independent of the shape of the hull located inside the cavity. At small values of the area of the circular channel between the cavity surface and the hull () equation (17) must be used. In particular, we show in the next sections that for the flows with ventilated supercavities their parameters restrictions depend on the pressure, ventilation rate, and the shape of the hull located inside the cavity.

7. Cavities on Different Hulls

Let us assume the shape of the hull located in the cavity to be a combination of functions as follows: Constants coincide with the values of the hull radius at ; that is, . Constants and can be chosen to ensure the continuity of the hull radius and its slope at ; . The hull shapes can be finite () and infinite (). Examples of such shapes are presented by the following formulas and shown in Figures 5 and 6 by solid lines:

By increasing the number of , arbitrary time independent hull shape can be interpolated by formulas (33) with a rather good accuracy. We will use these formulas in the next section to obtain a steady problem solution. Here we will try to answer the very important question: is it possible to cover any given hull by one steady cavity with the origin at the cross section ?

Let us consider the critical cavity shapes corresponding to the minimal possible values of the cavitation number . In particular, in the flow without gravity () the minimal possible cavitation number is zero at and the corresponding critical cavity shape can be determined from (27) as follows: In the case of the disc cavitator the critical cavity shape can be estimated with the use of (29) and (27) and continuity condition at : Examples of the critical cavity shape (37) are presented in Figure 5 by dashed line 1 () and in Figure 6 by dashed line 2 (). An example of the critical cavity shape (38) for the disc cavitator is presented in Figure 5 by dashed line 2.

For the base cavities at , the minimal possible cavitation number is and the corresponding conical critical cavity shape can be determined from (27) as follows: An example of the critical cavity shape (39) is presented in Figure 6 by dashed line 1 ().

If the hull pierces the corresponding critical cavity shape (e.g., dashed lines 1 in Figures 5 and 6), then it is impossible to cover such hull with one cavity starting at the cross section . Neither the cavity pressure increase (diminishing the cavitation number ) nor increasing the gas injection rate can create a cavity covering all the hull. It is possible only to cover the part of the hull located upstream to the piercing region (e.g., at for case shown in Figure 5) or to use other cavities, for example, the ones starting at the cross sections and on the hull shown in Figure 5.

For the vapor cavities on conical-cylindrical hulls ( in (30)), a similar limitation was revealed in Nesteruk [36]. It was shown that a cavity created by a slender cone can close on the cylindrical part only if (in accordance with (37)). If the radius of the cylindrical part is greater than value , all the cavities close on the conical part. For such cavities both the positive and the negative cavitation numbers are possible (the stability occurs in all these cases) and their shapes can be both convex and concave (Nesteruk [36]).

If the hull (35) touches the critical cavity shape only in one cross section , then the maximum length of the cavity closing on the conical part (or the minimal length of the cavity closing on the cylindrical part) is equal to (the cavity length is calculated along the axis of symmetry). Then the relationship between and the radius of the conical part coincides with (37) or (38) after replacing instead of and instead of . In Varghese et al. [37] the cavities created by the disk cavitator were calculated with the use of nonlinear approach and the dependence of the minimum cavity length versus the cylindrical part radius was revealed (no convergence was achieved at the values of cavity length smaller than ). The numerical results from the paper (Varghese et al. [37]) are shown in Figure 5 by markers. The rather good agreement between the data of Varghese et al. [37] and dashed line 2 supports the conclusions of the presented analysis and illustrates the accuracy of the Logvinovich formula (29) and the first approximation (27).

Any hull of bounded dimensions located under the corresponding critical cavity shape (like dashed lines 2 in Figures 5 and 6) can be entirely covered by a supercavity at zero gas injection rate by increasing the pressure inside the cavity or by reducing the cavitation number . If the cavitation number is not small enough, the corresponding supercavity closes on the hull.

In many cases the ventilation increases the supercavity dimensions at fixed (see Section 5, Manova et al. [15], and Nesteruk and Shepetyuk [16]). If the increase of the pressure at (or diminishing ) is impossible we could try to cover the hull by the cavity entirely only by increasing the gas injection rate. Equations (23) and (24) and a numerical solution of (23) obtained in Nesteruk and Shepetyuk [16] show that it is possible only for small enough values of . The ventilation decreases the effective cavity number (see (24)) and corresponding cavity must be larger (see (27) at ). On the other hand, according to the last term in (23) the ventilated cavity is smaller than the vapor one at . Thus, covering the hull by the cavity at small values of the gas pressure depends on the values of the ventilation rate and .

Let us illustrate this conclusion by an example of the numerical solution of (23) obtained in Nesteruk and Shepetyuk [16] for conical-cylindrical hulls ( in (30), , .8, and ; see Figure 7). The values of and correspond to the hull located under the corresponding critical cavity shape (in accordance with (37)). At the small fixed cavitation number it is possible to have a cavity which closes on the cylindrical part. In this case the cavity length can be rather good estimated with the use of linear theory—(27) at . In particular, the discontinuity of takes place (similar to the one revealed in Varghese et al. [37] and Nesteruk [36] for vapor cavities). If (or ) the cavity becomes unbounded. It means that ventilation rate is limited by the value

At the cavity closes only on conical part of the hull () and there is a large difference between the linear theory and nonlinear calculations based on (23). The gas injection rate is limited by some other critical value ( or ). Formally, a solution exists at greater values of ventilation, but the cavity becomes infinite at and such flow cannot be real.

Thus, covering the hulls located under the corresponding critical cavity shapes by ventilated cavities must be investigated with the use of solutions of nonlinear equation (20) or (23). In the next section it will be shown that in the case the resolving (20) can be reduced to the calculation of simple integrals.

8. Solution of the Steady Problem at High Froude Numbers

8.1. Integral Formula for the Solution

The nonlinear differential equation (20) or (23) can be solved numerically at any value of the Froude number and for any shape of the hull located inside the cavity. Let us consider the case and the hulls which can be presented or interpolated as (33). Substituting into (23) and using (33) yield The initial conditions (26) at can be rewritten as follows: Let us assume the function to be continuous at points ; . Its derivative can be discontinuous at these points if the slope of the hull contour is discontinuous (e.g., at for the hull shape shown in Figure 5). It means that, for every segment , , and the segment , the specific initial conditions should be used at the beginning of each segment , :

The order of (41) can be reduced by substitutions and . With the use of initial conditions (43) the following first-order equation can be obtained:

Equation (44) can be integrated with the use of initial conditions (43): Sign “+” or “−” must be used in (45) according to the value of .

The cavity closing corresponds to the zero value of , but from the physical point of view it is impossible to have the zero cross section area of gas flow in the steady ventilated cavity and from the mathematical point of view it is impossible to achieve the zero value of , since the differential equation (41) has a discontinuity. To avoid these difficulties we will suppose that the cavity “closure” corresponds to some small positive value of (e.g., ). Then the cavity length (along the axis of symmetry) can be also calculated from integral relations (45).

8.2. Maximum and Minimum Values of the Circular Channel Area

The inverse function may be many-valued (see curves 4 and 5 in Figure 8). In this case it is necessary to integrate in (45) from to corresponding to the maximum or the minimum value of the gas flow cross section area in the circular channel between the cavity surface and the hull and then to integrate from to with the opposite sign in (45). The values of can be easily calculated from the quadratic equation following from (44) at :

If , the maximum value of can be calculated from (47) as follows: An example of the cavity with the maximum value of (for , ) is presented in Figure 8 (curve 4). At the maximum of exists only for . In particular, at it follows from (22), (24), (46), and (47) that

The minimum of exists only for and and can be calculated from (47) as follows: An example of the cavity with the minimum value of is presented in Figure 8 (curve 6).

It must be noted that among the cavities presented in Figure 8 only curves 4 and 9 have classical elliptical shapes. Cavities 3 and 8 are parabolic, 2 and 7 are conical, 1 and 6 are hyperbolic, and 5 is cylindrical. If , it is possible to realize only cavities 4, 7, 8, and 9. Other cavities cannot exist in reality without conical part of the hull () according to the stability principle mentioned above. The cavity shapes corresponding to the zero cavitation number (dashed curves 3, 5, and 8) do not cover the hull shown in Figure 8. It means that at any value and at any ventilation rate all 3 cavitators (with , 0, and −0.1) cannot create a cavity large enough to cover the hull.

8.3. Degenerate Solutions, Neutral Hull Shapes, and Critical Values of the Ventilation Rate

Equation (41) can be rewritten as follows (with the use of (22), (24), and (43)): If then (51) has the obvious solution at any value of the ventilation rate. The inverse function is many-valued and such degenerate solution cannot be obtained with the use of (45). Relationship (33) in view of (43) and (52) allows calculating the corresponding body shape For example, for the part of the hull, corresponding to and shown in Figure 8, . Therefore (53) yields and . The degenerate solution is shown in Figure 8 by straight line 5. The hull (53) ensures the same cavity shape at any ventilation rate. This shape coincides with the shape of the vapor cavity at the same values of parameters and and can be refereed as neutral.

Equation (20) demonstrates that for the ventilation increases the cavity dimensions. For example, at the dimensions of cavities 1–4 shown in Figure 8 for will increase with increasing the ventilation rate. The maximum cavity radius can be also infinite, if the cavity covers all the hull or the hull is unbounded (e.g., ). Formulas (46) and (48) allow calculating the critical ventilation rate corresponding to zero value of If inequality (49) is valid, the cavity is bounded at the ventilation rates (54) but for the value of in (48) can be zero or negative. The equation yields the second critical value, which with the use of (46) can be represented as follows: If , the corresponding cavity is bounded, but unbounded at . It means that the ventilation rate is limited by the value . For unbounded cylindrical hulls the critical ventilation rates were calculated in Manova et al. [15] and coincide with (54) and (55).

On the other hand, for the ventilation decreases the cavity dimensions. For example, the dimensions of cavities 6–9 shown in Figure 8 for will decrease with increasing the ventilation rate at . In particular, according to the formula (50) the decreases from value at to at . The equation yields the critical value, presented by (55). If , the corresponding cavity has a minimum, but increasing the ventilation rate over the critical value diminishes the cavity dimensions and its shape has no minimum.

8.4. Examples of the Ventilated Cavity Shape Calculations

If the pressure in the cavity at some cross section is fixed (e.g., the cavitation number is fixed) and (the area of the circular gas channel diminishes), the increase of the gas injection rate yields the diminishing of the cavity dimensions (since local pressure in the cavity diminishes). This fact is illustrated by the calculation examples shown in Figures 9 and 10. In both cases and the increasing of the gas injection rate decreases the cavity dimensions at . To calculate the cavities at the values and were used according to formulas (43). The values are also negative for all cavities shown in Figure 10.

Curve 4 shown in Figure 10 corresponds to the critical value of gas injection rate calculated from (55). At greater values of (curves 5–7) cavity shape has no minimum. At smaller and zero values of the gas injection rate (curves 1–3) cavity shapes have a minimum point and some of them (curves 2 and 3) have also a maximum located near the closing on the conical part of the hull. For and the hull presented in Figure 8 corresponds to the neutral shape (53) at . It means that at the corresponding cylindrical cavity shape (shown in Figure 8 by curve 5) is the same at any values of ventilation. Since the value is negative, at , at this segment the cavity radius diminishes as the gas injection rate increases. The examples of calculations are shown in Figure 11.

Absolutely different behavior occurs in the case , when the area of the circular gas channel increases. The dimensions of the ventilated cavities increase as the gas injection rate increases. This fact is illustrated by the calculation examples shown in Figures 12 and 13. In both cases . To calculate the cavities at the values and were used according to formulas (43). At great enough values of the cavitation number , the cavity shapes transform from convex with the maximum thickness section at small values of the gas injection rates (see Figure 12, curves 1 and 2) to the concave at greater values of . At the cavity shapes are concave (see, e.g., Figure 13). All the cavities cannot be extended for (according to the stability principle); thus the values of gas injection rate are limited. For example, for the case shown in Figure 12 the critical value . At smaller values of the critical values of the gas injection rate diminish (e.g., at    (see Figure 13) and at   ). It must be noted that for the hull shape shown in Figures 813 the cavities corresponding to the critical values (54) and (55) are bounded (since is restricted). For example, the cavity shape corresponding to the critical value of the ventilation is shown in Figure 12 by curve 2.

8.5. Cavity Closing Near the Hull Contour Discontinuities

A special attention should be paid to the cases when the cavity end passes through the hull contour discontinuity as shown in Figures 14 and 15. If or (see (43)), the cavity length varies continuously. This case is shown in Figure 14. Absolutely different is the case shown in Figure 15. At some critical values of parameters and the cavity contour touches the hull at the point , and a “jump” in the cavity length is possible (see Figure 15).

These two cases of cavity behaviour are separated by the condition at or At formulas (56) and (57) allow obtaining the relationships between the critical values of parameters , , and . With the use of (27), equation (56) can be written as follows: where the critical value of the cavitation number may be determined from (57) as follows: Finally, from (58) and (59) the following simple relationship can be obtained:

In the particular case of a conic-cylindrical hull (), (60) coincides with the relationship proposed in Nesteruk [36]. When the value of is smaller than the critical one (60) and , the cavity length varies continuously as shown in Figure 14. Otherwise the cavity length has a discontinuity as shown in Figure 15. For the ventilated cavities the critical values of parameters , , and depend on and should be calculated with the use of (45). For the conic-cylindrical hull some examples can be found in Nesteruk and Shepetyuk [16] and are shown in Figure 7.

In the case a rather weak dependence of the cavity length on the ventilation rate should be expected. Even in the case when cavity dimensions increase with increasing at , the negative value of leads to a decrease of the cavity length at . This fact is illustrated in Figure 16 by the cavities on the conic-cylindrical hull (, , and ) created by a disc cavitator. The same cavitator and hull were tested in Vlasenko and Savchenko [13].

The value is smaller than the critical one based on (60); therefore no “jumps” in cavity length (as shown in Figure 15) are expected. At the greater value of the cavitation number the corresponding values of are negative and ventilation does not increase the cavity length significantly (see Figure 16). On the contrary, at the smaller cavitation number = 0.06 the values of become positive and the same values of the ventilation rate drastically increase the cavity length (see Figure 16). In the experiments performed in Vlasenko and Savchenko [13] the values of were smaller than 0.0001. The calculations show that the cavity shapes are very close to the vapor ones (shown in Figure 16 by dashed lines) at these small gas injection rates. Thus, the changes in the pressure inside the cavity (or different values of ) are the main reason of changes in the cavity length in the experiments (Vlasenko and Savchenko [13]). Unfortunately, the pressure inside the cavity was not measured in Vlasenko and Savchenko [13].

If at a zero or a small value of the cavity covers a hull, its length approaches the value and at (in particular, for a bounded hull ), then the value is negative (according to formula (43)). In this case the cavity length depends on the ventilation rate very slightly. This fact was revealed also in experiments (see Wosnik et al. [10] and Vlasenko and Savchenko [13]).

9. Asymptotic Solution at Small Values of the Ventilation Rate

9.1. Basic Equations

If the gas injection rate is small enough , the solution of the nonlinear equation (17) can be expressed as follows: where is the solution of (17) at zero ventilation rate , Function is a solution of the following equations: To solve the nonlinear integral-differential equation (62) the asymptotic series obtained in Nesteruk [30] and the standard initial and boundary conditions (19) can be used for the function . If the accuracy restrictions are not very severe, the first approximation equation (similar to (9) or (10) with main terms of order only) can be easily solved. Knowing the function the right hand part of (64) and (65) can be calculated and the function can be determined as a solution of the set of linear differential equations with the partial derivatives (64) and (66) and zero initial and boundary conditions. Thus the presentation (61) allows avoiding the nonlinear integral-differential equations.

9.2. Stability of the Steady Ventilated Cavities

Let us use the main term in (61) and the first approximation in (62) to analyze the stability of ventilated cavities at small values of the gas injection rates. Let us consider the steady cavity with a constant gas pressure inside the cavity . According to the first approximation of (62) the cavity radius is given by formula (27). Let us consider small changes in gas pressure , where is a new constant gas pressure inside the cavity. According to (27) and (1) the new cavity radius can be written as follows: Here we neglect the changes in the vapor pressure.

Integration of formulas (27) and (67) allows calculating the changes in the cavity volume : where is the initial cavity length, which is measured along the axis and can be found from (27) as follows: Here we suppose the body radius to be time independent; therefore yields also the changes in the volume of gas.

On the other hand for a fixed mass of gas its pressure and volume are related according to the polytropic relationship: where (in particular, for isothermal process ); the initial volume of gas can be calculated with the use of (27) as follows: Putting in (70) and and linearization yield Values and have opposite signs (see (68) and (72)) and can compensate each other, if . Otherwise the cavity becomes unstable. The critical value of the gas pressure can be obtained from , (68), and (72): At smaller values of the gas pressure the cavity is stable. If , the cavity becomes unstable. The existence of the critical gas pressure value was revealed in experiments (e.g., Michele [20]) and supported by the theory developed by Paryshev [23, 24].

Equation (73) can be rewritten in the dimensionless form taking into account the definition of the cavitation number (1): where the vapor cavitation number and the critical value of the cavitation number can be obtained from (69) as follows: It follows from (71), (73), and (76) that Formula (77) allows calculating the critical value of for every value of the cavity length. The corresponding critical cavitation number can be calculated with the use of (76). Thus the dependence () can be obtained. In particular, for an empty isothermal cavity at large Froude numbers (, , and ) it follows from (77) and (76) that

Examples of calculations based on (78) are presented in Figure 17. Stable cavities correspond to regions located above the corresponding curves. In the case of the slender cavitators, the function () depends on in comparison with the relationship obtained by Paryshev [23] for the isothermal cavities created by nonslender cavitators with the use of semiempirical formulas for the cavity shape. In the cases of long cavities () and for the base cavities with the relationships (78) yield the formulas respectively. The difference between the first relationships (80) and (79) can be explained by the limited accuracy of both the first approximation equation (27) and the semiempirical formulas used by Paryshev. In the experiments with 2D cavities the value was measured by Michel [20]. It can be seen from Figure 17 that the largest stability regions correspond to the base cavities (). In particular, stable cavities with negative cavitation number are possible at . It follows from (78) that the cavitation number is limited by value , obtained in Nesteruk [32].

The influence of the gravity forces is shown in Figure 18. For the conical cavitators the presence of gravity diminishes the stability regions, if its direction coincides with the direction of the ambient flow (see curves 1 and 2). In particular, there are no stable cavities at some small values of . Opposite behaviour is typical for the base cavities (see curves 9 and 10). If the directions of the ambient flow and the gravity are opposite the minimal cavitation number exists (this fact was already discussed in Section 6). The critical cavity shapes are similar to the solid line shown in Figure 4. For such lines ; therefore the minimum values of the cavitation number can be calculated from (76) with the use of the condition . These minimum possible values of the cavitation number correspond to the ends of the curves 4–7 (see Figure 18) and yield additional stability restrictions.

The cavities with the hulls located inside have smaller stability regions in comparison with the empty cavities (due to the smaller value of (see (73)). The same conclusion has been drawn out by Paryshev [23]. To illustrate this fact, let us take the neutral hull shape: For such hull formula (27) is a solution of (20) at any value of the gas injection rate; that is, (27) is not only the approximation of the steady cavity shape at a small ventilation but also the exact solution at any gas injection rate. If , (81) coincides with the relationship (53). The length of the neutral hull can be obtained from (81) as a solution of . Let us assume that at . Then and for the isothermal gas processes inside the cavity, (77) can be rewritten as follows:

Let us consider at first the case of very small area of the circular channel between the cavity and the hull surface . Then and formulas (82) and (76) yield Therefore, cavities with the very small area of the circular channel between the cavity and the hull can be stable at the cavitation numbers which are very close to maximum possible value, . Thus such cavities are mostly unstable.

Let us consider the case of the maximal area of the circular channel, that is, . For the base cavities ( and ) the hull volume tends to zero. The stability of the empty base cavities has been already investigated before and illustrated in Figures 17 and 18. It can be seen that the stability regions are rather large. In the case and , the corresponding cavities are not empty. Let us investigate their stability at . It follows from (76) and (81) that Substituting (84) into (82) yields that at fixed value of the corresponding value of is smaller by the value ; therefore the corresponding stability regions are smaller in comparison with the case of empty cavities (shown in Figure 17 by curves 1–3). The examples of calculations based on (78), (81), and (82) are shown in Figure 19. It can be seen that grater values of the volume of the hull part located inside the cavity (larger values of ) decrease the stability region. The most stable are the empty cavities presented in Figure 19 by dashed lines 7 and 8.

It must be noted that the experiments at standard room temperature and at values of the pressure close to the atmospheric one yield rather grate values of (e.g., at  m/s). At these values of the corresponding critical cavitation number is rather large and cavity length is rather small. Thus, it is difficult to obtain long stable empty cavities in the low speed experiments. Vlasenko and Savchenko [13] revealed stable cavities closing on the conical part of the hulls. Their length is restricted by the value , obtained at some gas injection rate for the hull shown in Figure 16. Further increasing of the ventilation leads to the drastic increase of the cavity length. Since the values of gas injection rate are too small to change the cavity shape significantly, such behaviour can be treated as a loose of stability.

Let us try to estimate corresponding to the experimental critical value with the use of (77) and compare with the experimental value (Vlasenko and Savchenko [13]). Since the experiments were performed for the disc cavitator, we need to use formula (29) for the initial part of the cavity at and take into account the corresponding volume of the gas . We assume that at , , , and . The calculations based on the modified formula (77) yield which is very close to the experimental value.

The obtained results can also be applied for the separated bubbles. In incompressible liquid and according to (72) such bubble is always unstable. In the case of a gas, the same gas is present in the external flow and inside the separated zone (there is no more a huge difference in densities as for a gas cavity located in the water flow). Nevertheless, the small velocities inside the separated zone allow supposing the pressure to be constant and using the developed theory. In this case the vapor cavitation number coincides with the Euler number as follows: According to the presented analysis, is difficult to obtain long enough stable separated bubbles for , but at smaller values of the Euler number the separated zones can be stable. Therefore, the Euler number must be taken into account in simulations of the gas flow separation (as mentioned in Nesteruk [25]).

9.3. Stability of the Pulsating Ventilated Cavities

In the previous section it was shown that steady stable cavities at low speeds ( m/s) have to be very short. In order to cover long hulls it is possible to use pulsating cavities which were investigated experimentally in Silberman and Song [18], Song [19], Michel [20], and Zou et al. [21] and by many other authors. The theories of such empty cavities were proposed by Woods [22], Paryshev [23, 24], Nesteruk [25], and Semenenko [26]. Here we consider the small oscillations of the pressure inside the cavity which covers the hull.

Let us suppose that the gas injection rates are small enough to use only the main term in (61) and to consider the pressure inside the cavity to be only time dependent. The corresponding changes of the cavitation number can be represented as follows: where is the average value of the cavitation number and the dimensionless time is based on the period of oscillations .

Let us suppose the velocity of the ambient flow and the hull shape to be time independent and use the radius of the cavity and the body at the cross section of cavity origin (see Figure 1) for dimensionless lengths and to calculate the parameters (4) and (7). Then Let us take only the main term of order in (62) and use the analytical solution obtained in Nesteruk [28]. In particular, for the harmonic low (86) the following formula was obtained (Nesteruk [33]): Formula (88) shows that for the pulsating cavity differs from corresponding to the steady cavity at the mean value of the cavitation number by function . Formula (89) corresponds to a propagating wave with the amplitude increasing downstream. The length of the wave—the minimal distance between maxima or minima of the function (89)—is equal to . The same behaviour was revealed in experiments (Song [19]).

Integration of (88) allows calculating the cavity volume: where the cavity length may be obtained from At small values of (or ) the difference between and the length of the steady cavity at the cavitation number must be small, since is a solution of (92) at . Here we do not take into account the cases of discontinuous changes in the cavity length discussed in Section 8.5 and shown in Figure 15.

By putting into (90) the pulsations of the cavity volume can be calculated: where is the average volume of the steady cavity at the cavitation number and can be obtained from (90) at . With the use of (93) the average value of the cavity volume pulsations can be calculated: Equations (94) and (95) show that the average value of the cavity volume pulsations is proportional to the amplitude of the gas pressure pulsations and can be diminished, when . It can be expected that the pulsating cavities with zero average value of the volume pulsations are stable; that is, they can exist during a long time. According to the formula (95) the zero values of are possible at Equation (96) shows that stable oscillating cavities occur only at integer values of the Strouhal number or at integer number of waves located on the average cavity length. This fact was revealed in experiments on empty 2D cavities (Song [19]) and supported by the theories developed by Woods [22] and Semenenko [26] for 2D cavities and by Paryshev [23, 24] and Nesteruk [25] for the axisymmetric empty cavities. Presented analysis illustrates that the same behaviour is typical also in the case of rigid hulls located inside the cavity. If the average cavity length is fixed, the corresponding frequencies of pulsations in the cavity and in its wake are limited by (96).

Formula (96) can also be applied for the separation bubbles both in gases and in liquids. In particular, let us compare (96) with the experimental Strouhal numbers for Karman vortex street, measured for circular cylinders of diameter at large Reynolds numbers (see, e.g., Theodore von Kármán [38]): Since the cylinder diameter is few times smaller than the average length of its separated zone, formulas (96) and (97) are in good agreement.

In order to investigate the influence of the ventilation rate on the pulsating cavities, (64) must be solved. But some conclusions can be drawn out with the use of simple physical considerations. At the local minima of the circular channel area the pressure in gas flow and cavity radius decrease and in maxima the pressure in gas flow and cavity radius increase. It leads to the increasing of the pulsation amplitude with the increasing of the ventilation rate. This fact was revealed by experiments (see, e.g., Silberman and Song [18]).

10. Conclusions

The integral-differential equation was obtained to simulate unsteady evolutions of the slender axisymmetric ventilated supercavity with the use of one-dimensional inviscid flow of the incompressible gas in the channel between the cavity surface and the body of revolution. For small ventilation rates, a solution was expressed as asymptotic series. In steady case the nonlinear differential equation and its solutions were obtained. It was shown that the cavity shape depends strongly on the values of the ventilation rate and the cavitation number.

The ventilation can strongly increase and decrease the cavity dimensions. In particular, for the cavities created by the slender conical cavitator the increase of the ventilation rate increases the cavity dimensions. The cavity could become unbounded at some critical values of the ventilation rate. On the other hand, the ventilation can strongly decrease the base cavity length. Some specific neutral hulls ensure no changes in the cavity shape at any value of the ventilation.

Examples of the steady cavity shape calculations for different hulls are presented. It was shown that the cavities can have classical elliptical shapes and be also conical, hyperbolic, and cylindrical. They may also have the minimum of the radius. Peculiarities of the cavity closing near the hull contour discontinuities were investigated.

Simple stability theories for the steady and pulsating cavities at small ventilation rates were proposed. The results are in good agreement with the known theoretical and experimental findings of other authors and generalize them for the cases of the hull, located inside the cavity, and different cavitators and take into account the influence of the gravity. In particular, it was shown that the cavities with a very small area of the circular gas channel between the hull and the cavity surface are unstable and that gravity and the cavitator shape strongly influence the stability. The fact that stable pulsating cavities are possible, when an integer number of waves are located on the average cavity length, is also valid for nonempty cavities.

More complicated models are necessary to simulate the gas flow inside the cavity. In particular, the viscosity of gas has to be taken into account.

Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.