Discrete Continuous Models 2019 27 (1) 49—59
. „ . http://journals.rudn.ru/miph
Applied Computational Science
Mathematical modeling
Research article
UDC 519.63
DOI: 10.22363/2658-4670-2019-27-1-49-59
Mathematical model of cavitation under the influence of a single stretching pulse
Nikolay Yu. Kravchenko1, Dmitry S. Kulyabov2
(received: August 6, 2019; accepted: October 21, 2019)
This paper describes the created mathematical model that allows you to explore the dynamics of cavitation bubbles under the influence of a single negative pressure pulse. The time dependence and coordinates of the parameters of the carrier phase, the temperature and pressure of the vapor phase, the concentration and size of the bubbles are determined numerically. It is concluded that the model created gives a good agreement between the calculated and experimental data.
Cavitation in a fluid is a process in which a dynamic decrease in pressure inside a fluid occurs at a constant temperature, as a result of which vapor and vapor-gas bubbles increase in the fluid caused by the evaporation of fluid into these bubbles. Bubbles (caverns) are formed in those places where the pressure in the liquid becomes below a certain critical value [1]. In a real fluid, it is approximately equal to the saturated vapor pressure at a given temperature [2].
Some features of cavitation
— Cavitation is peculiar only to liquids and does not occur under normal conditions, either in solids or in gases.
© Kravchenko N. Yu., Kulyabov D. S., 2019
This work is licensed under a Creative Commons Attribution 4.0 International License http://creativecommons.org/licenses/by/4.0/
— Cavitation results from a decrease in pressure in a liquid, which means that it can be controlled by adjusting the pressure. If the pressure is maintained below a certain level for a long time, cavitation occurs.
— Cavitation is associated with the disappearance of cavities and their appearance in the liquid.
— Cavitation refers to non-stationary phenomena, since is a process of growth and collapse of cavities.
— Cavitation can occur both in the case of a moving fluid and in the case of a resting one.
— Cavitation can occur both in the volume of a liquid and on the solid boundary.
The equation describing the dynamics of cavitation bubbles in an incompressible fluid without taking into account the vapor pressure, surface tension and viscosity of the fluid can be represented as:
RR + - R2 = -2 p
pb+Ro) (R) - 2R- po - p(t)
Here: R0 is the radius of the nucleus at t = 0; R — radius of the nucleus at the next time instant t; p — density of a liquid; a is the surface tension of the fluid; k = 1 is the adiabatic index for steam in the bud; po is the hydrostatic pressure in a liquid (po = pb); R is the acceleration of the cavity wall; R is the speed of movement of the cavity wall; 2a/R0 is the Laplace pressure; R0/R is the amplitude of oscillations of the cavity.
A single impulse is presented in the form:
p(t) = -a (t + ti) exp( - tt+1) , (2)
where t\\ is the time of appearance of the first germ of homogeneous cavitation under the action of this pulse. In this case
| = (--) a (t + ti) exp (-+ a exp (-. (3) At t = 0, the rate of increase in pressure at the initial moment of time is
dp t-, , , dp , x ar
a = ——. From the condition — = 0, applied to (3) we find tm = t, pm = —,
where t is the rise time of the pressure pulse up to the maximum value pm.
Let us set the maximum pressure amplitude pm so that it approaches 95% to
the thermodynamic stability [3] of the fluid (spinodal, [4]), i.e.
pm = 0.95(pb — ps). (4)
Then from (4) it follows that
a = 10.95e(pb — Ps). (5)
Thus, the rate of pressure increase a can be calculated for a given fluid temperature using the value of (pb — ps) for a given t.
The time of appearance of the first nucleus t1 is found from the condition (2), when at t = ti the pressure is p = pi.
pi = —ati ■ exp ^—= —ft(pb — ps). (6)
Here p1 is the pressure at which one germ appears in 1 cm3 at the moment of time t1,
ft =--— < 1. (7)
pb — ps
Substituting a into (6), we get
ft = e ■ — ■ exp ^ ——^ = e ■ x ■ exp (—x) (8)
where x = -1/t.
Thus, in the equation (3), for a given value of t, the value of a; is calculated by the formula (5) and the value of t1 is calculated by the formula (2) for a given value of ft (formula (8)).
To determine the time t1, it is necessary to take into account the number of bubbles appearing per unit of time per unit volume of liquid:
B ■ exp
. = 16na3 = 3(1 — pv/pi)2 , ( )
pb — p = at,
dN = B ■ exp ^--I dt,
N =1 = B JexJ — )dt. (10)
By calculating a and A through the formulae (9) and (10) by integrating respectively and assuming that B ~ 1033 cm-3s-1 at N = 1 the time t1 is calculated.
Denoting by we get that
c = - A,
B • I exp ( - )dt = 1.
Taking into account the geometric meaning of a definite integral (11), one can determine the point t1, by numerical method knowing that the area of the figure bounded by the function f (x) on the interval [0, t1] should be equal to 1.
In the basic cavitation equation (1), we substitute the expression for pressure in the form (2). Then the basic cavitation equation (Rayleigh equation) takes the form:
RR + - R2 =
pb + Ro)
Hi - R - Pb - a (t +t1) exn - t
t +11
The differential equation for the growth of an embryo vapor has the form:
Pi\\ rr+- a
-*> - 2a (1 - pi
where p(t) is defined by the formula (2). Substituting p(t) from (2) into (13), we get:
Pi[ RR + - R2
a(t + t1) • ex^ t +11
For the reasons described above, we will assume t1 0, and then the equation (14) takes the form:
pi\\ rr + - r2
a •t •ex^- 7)- f i1 - fv
It is solved by us numerically. For this, it is necessary to present the last equation in the form of a system of two equations in which the following
change of variables is introduced:
R = u(2), R = du(2) = u(1), R = du(1).
The equation (15) by changing variables is converted to the form:
at exp( - I)- ï 0-1
atexpl -r) - U(2)
The last equation is solved numerically by the Runge-Kutta method. The
calculations were carried out for t =10
It was believed that at the initial time, a homogeneous liquid that does not contain vapor bubbles was at a given temperature and pressure
t = 0, pi = pio, Ti = Tio, Nb = 0,
The rate of pressure drop ap at the point with coordinate x = 0 was set by the boundary conditions
After a bubble appeared, the following initial conditions were accepted for ordinary differential equations describing its development (12): R = Rc,
pv = pb, tv = Ti.
We introduce the following notation, convenient for working with the program for the numerical solution of the basic cavitation equation in the form (1):
uR0 = omr0, —- = sr0, Pb = pb, Ro
(uR0)2 = omr02, ( Pb + — ) = pbsr0, Pa = pa,
p(uR0)2 = rhomr, = fio, R* = u(2), z * = u(1).
Then the original equation (1) takes the form of a system of two first order differential equations:
du(1) =----+
v ; 2 omr0 u(2)
rhomr u(2) du(2) = u(1).
— sr0—— — pb + pa • sin(r + jio)
(u(2))3k u(2)
The resulting system of differential equations is solved numerically by the Runge-Kutta method. It should be noted that numerical solutions were cited earlier for other cases [5-9].
We have created a program for the numerical solution of the cavitation equation in the Compaq Visual Fortran Professional programming language. Its work is based on the Runge-Kutta method.
Initially, the main program asks for the values of external parameters, such as fluid temperature, oscillation frequency, and others (see Fig. 1). Then the main program refers to an array of tabular data for the values of surface tension, fluid viscosity, fluid pressure, vapor pressure at a given temperature. These tabular data are discrete values and do not always correspond to a given temperature. Therefore, the main program refers to auxiliary subroutine 1, which approximates or extrapolates the table data to a given point.
Initial conditions j|
Ail ay of tabulai& data
Figure 1. Block diagram
To calculate parameters such as the initial radius of the cavity, the pressure at which the first cavitation nucleus appears, the initial phase of external oscillations, the main program refers to subroutine 2, which calculates these values based on the data already calculated by subroutine 1.
Subroutine 3 then receives from subprogram 2 a task to calculate the time t\\ during which the first cavitation nucleus appears in the fluid. The required
tabular data is requested from subroutine 1. The result of the calculation is reported to the main program.
Having collected all the necessary data, the main program calculates the basic cavitation equation for the maximum amplitude of oscillations of the cavity in case of acoustic cavitation [10,11]. It should be noted that this case is also interesting in the possibility of initiating a nuclear fusion reaction [12-14], which is confirmed by the theory [15,16].
Below is a part of the main program for the numerical solution of this system of equations, written in the programming language Fortran:
external rad,res
dimension pt(5), u(2), du(2), aux(8,2)
common omr0, rmr0, rhomr, pbsrO, pk, sr0, pa, pb, pi, nk,
- fa, r0, p, pi, fi0, ky open(1, file= &p.dat&) open(2, file= &r.dat&) open(3, file= &bubble.txt&)
print *, &The program for calculating the dynamics of the
- bubble with acoustic effects and the frequency of Fa
- according to Runge Kutta method.&
print *, &-------------------------------------------&
print *, &time pressure dependence P = P (t / T) is
- written to the file "p.dat", the time dependence of
- the radius R / R0 = f (t / T) is written to the file
- "r.dat".&
print *, &THE PROGRAM WILL NOT WORK without an auxiliary
- file "bubble.txt" and also without the additional
- program "rkgs.for"!&
read (3, *) qi, q2, fa, pa, pb, pi, r0, sig, amu, ro, kt,
pi = 3.141592654 pk = 5
print *, &The following parameters are entered here:& print *, &1) Time step integration (wt): pt (3) =&, qi print *, &2) Integration error: pt (4) =&, q2 print *, &3) Acoustic frequency: Fa =&, fa, &Hz& print *, &4) Pressure amplitude at a given point: Pa =&,
- pa, &Pa&
print *, &5) Pressure on binodals at a given point: Pb =&,
- pb, &Pa&
print *, &6) Pressure at the point of emergence of the
- unit: P1 =&, pi, &Pa&
print *, &7) Starting radius of the nucleus: R0 =&, r0,
print *, &8) Surface tension: sigma =&, sig, &N / m&
print *, &9) Molar mass: mu =&, amu, &kg / mol&
print *, &10) The density of the fluid at this point: po
- =&, ro, &kg / m * 3&
print *, &11) Number of periods studied: n =&, kt print *, &12) Number of points for one period: N =&, ky print *, &&
print *, &Changing these parameters is possible only in ^ the file& bubble.txt &&
print *, &To interrupt the program, press "Ctrl + C", for
continuation - "Enter"& pause
om=2.0*pi*fa
omr0=om*r0
omr02=omr0*omr0
rhomr=ro*omr02
rmr0=4.*amu/(r0*ro)
sr0=2.*sig/r0
fi0=asin((pb-p1)/pa)
pbsr0=pb+sr0
pt(1)=0.0
pt(2)=2.0*pi*kt
call rkgs (pt,u,du,2,ih,rad,res,aux) if(ih-10)3,3,4
print *, & Error code &,ih goto 5
print *, & ATTENTION !!!&
print *, &ERROR CODE SHOULD NOT EXCEED 10.&
print *, &Correct, please, the parameters pt(3),pt(4)&
read(*,*)q1,q2
Below are some results of a numerical simulation of the bubble behavior in a rarefaction wave.
The mathematical model created by us allows us to investigate the dynamics of cavitation bubbles with a change in fluid pressure. Numerical simulation made it possible to determine the time dependence and coordinates of the parameters of the carrier phase, the temperature and pressure of the vapor phase [17], the concentration and size of the bubbles. The proposed model is applicable not only to liquids, but also to metals in the liquid phase [18].
We concluded that the model created gives a good agreement between the calculated and experimental data [19], which demonstrates the applicability of the approach under consideration to the problem of rapid pressure drop.
In the studied problem, the following picture takes place: the pressure of the fluid drops in the rarefaction wave to a value below the saturation pressure, the fluid enters a metastable state [20], intense nucleation begins in the region of minimal pressure, after which the bubbles rapidly increase due to the interfacial mass transfer, which ultimately account stabilizes the pressure at a value close to the saturation pressure.
Acknowledgments
The publication has been prepared with the support of the "RUDN University Program 5-100".
References
URL http://mi.mathnet.ru/eng/uzkui365
For citation:
N. Yu. Kravchenko, D. S. Kulyabov, Mathematical model of cavitation under the influence of a single stretching pulse, Discrete and Continuous Models and Applied Computational Science 27 (1) (2019) 49-59. DOI: 10.22363/26584670-2019-27-1-49-59.
Information about the authors:
Nikolay Yu. Kravchenko (Russian Federation) — Deputy Director of Institute of Physical Research and Technology, Peoples& Friendship University of Russia (RUDN University) (e-mail: kravchenko-nyu@rudn.ru, phone: +7(495)9550839, ORCID: https://orcid.org/0000-0003-3397-1746, ResearcherID: E-6162-2018, Scopus Author ID: 14633789300)
Dmitry S. Kulyabov (Russian Federation) — Doctor of Sciences in Physics and Mathematics, Full Professor, Department of Applied Probability and Informatics, Peoples& Friendship University of Russia (RUDN University), (e-mail: kulyabov-ds@rudn.ru, phone: +7(495)9520250, ORCID: https://orcid.org/0000-0002-0877-7063, ResearcherID: I-3183-2013, Scopus Author ID: 35194130800)