Mathematical model of cavitation under the influence of a single stretching pulse Текст научной статьи по специальности «Физика»

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

1 Institute of Physical Research and Technology Peoples' Friendship University of Russia (RUDN University) 6, Miklukho-Maklaya str., Moscow, 117198, Russian Federation 2 Department of Applied Probability and Informatics Peoples' Friendship University of Russia (RUDN University) 6, Miklukho-Maklaya str., Moscow, 117198, Russian Federation

(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.

Key words and phrases: numerical simulation, cavitation, mathematical model of cavitation

1. Introduction

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.

2. Cavitation equation

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:

- • 1

RR + - R2 = -2 p

3 k

pb+Ro) (R) - 2R- po - p(t)

(i)

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.

3. Dynamics of a cavity under the action of single pulses of negative pressure

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 = —,

at at e

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)

T

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:

dN V-dt

B ■ exp

16na3

3(1 — pv / pi )2(pb — p)2 V =1 sm3,

. = 16na3 = 3(1 — pv/pi)2 , ( )

pb — p = at,

dN = B ■ exp ^--I dt,

N =1 = B JexJ — )dt. (10)

o ^

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,

a2

11

C

B • I exp ( - )dt = 1.

(11)

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:

3 •

RR + - R2 =

1

P

2a\

pb + Ro)

R

3k

2a

Hi - R - Pb - a (t +t1) exn - t

t +11

. (12)

The differential equation for the growth of an embryo vapor has the form:

Pi\ rr+- a

-*> - 2a (1 - pi

1

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

2a ~R

(13)

^ -

1

(14)

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

-1

(15)

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:

R

pi R

at exp( - I)- ï 0-

-1

3 R2

or

du(l)

pi u(2)

, A 2a

atexpl -r) - U(2)

l - p

Pv

1

2 Rpi

3 [u(l)]2

2 u(2)pi

(16)

(17)

The last equation is solved numerically by the Runge-Kutta method. The

-9

s.

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,

u

0.

The rate of pressure drop ap at the point with coordinate x = 0 was set by the boundary conditions

x

0, pi = pio - ®p • t,

du dx

After a bubble appeared, the following initial conditions were accepted for ordinary differential equations describing its development (12): R = Rc,

pv = pb, tv = Ti.

l

l

4. Method of numerical solution of the equation

We introduce the following notation, convenient for working with the program for the numerical solution of the basic cavitation equation in the form (1):

2a

uR0 = omr0, —- = sr0, Pb = pb, Ro

(uR0)2 = omr02, ( Pb + — ) = pbsr0, Pa = pa,

Ro

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:

3 1 u(1)

du(1) =----+

v ; 2 omr0 u(2)

1

+

rhomr u(2) du(2) = u(1).

pbsrO 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].

1

5. The program for the numerical solution

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,

- ky

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,

- 'm'

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

nk=1

6 continue pt(3)=q1 pt(4)=q2 u(1)=0.0 u(2)=1. du(1)=0.5 du(2)=0.5

call rkgs (pt,u,du,2,ih,rad,res,aux) if(ih-10)3,3,4

3 continue

print *, ' Error code ',ih goto 5

4 continue

print *, ' ATTENTION !!!'

print *, 'ERROR CODE SHOULD NOT EXCEED 10.'

print *, 'Correct, please, the parameters pt(3),pt(4)'

read(*,*)q1,q2

goto 6

5 continue stop end

Below are some results of a numerical simulation of the bubble behavior in a rarefaction wave.

6. Conclusions

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

1. N. Yu. Kravchenko, M. M. Martynyuk, Dynamics of homogeneous cavitation bubbles in water under large amplitude pressure oscillations, Bulletin of Peoples' Friendship University of Russia. Series: Physics 8 (2000) 118-121, in Russian.

2. A. A. Aganina, M. A. Ilgamov, D. Yu. Toporkova, Dependence of vapor compression in cavitation bubbles in water and benzol on liquid pressure, Proceedings of Kazan University. Physics and Mathematics Series [Uchenye zapiski kazanskogo universiteta. Seriya fiziko-matematicheskie nauki] 158 (2) (2016) 231-242.

URL http://mi.mathnet.ru/eng/uzkui365

3. M. M. Martynyuk, N. Yu. Kravchenko, Limit of thermodynamic stability of a liquid phase in the field of nefative pressure [Granitsa termodi-namicheskoy ustoychivosti zhidkoy fazy v oblasti otritsatel'nykh davleniy], Russian Journal of Physical Chemistry 72 (6) (1998) 998-1001, in Russian.

4. P. A. Tamanga, M. M. Martynyuk, N. Yu. Kravchenko, Spinodal of liquid phase on basis of generalized berthelo's equation, Bulletin of Peoples' Friendship University of Russia. Series: Physics (9) (2001) 56-58, in Russian.

5. N. Yu. Kravchenko, The numerical solution of the Rayleigh-Plisset equation for spark cavitation and calculation of the maximum temperature and pressure in a cavity, Journal of Mechanics of Continua and Mathematical Sciences (Special Issue-1) (2019) 465-473. doi:10.26782/jm-cms.2019.03.00046.

6. E. V. Butyugina, E. S. Nasibullaeva, Numerical study of the gas diffusion process between clustered bubbles and technical fluids [Issledovaniye protsessa diffuzii gaza mezhdu puzyr'kom v klastere i tekhnicheskoy zhidkost'yu], News of the Ufa Scientific Center [Izvestiya UNTS RAN] (2) (2016) 15-21, in Russian.

7. A. B. Kapranova, A. E. Lebedev, A. M. Meltser, S. V. Neklyudov, S. E. M., Methods of modeling the developmental stages of hydrodynamic cavitation, Fundamental Research [Fundamental'nyye issledovaniya] (3) (2016) 268-273, in Russian.

8. Modelling of cavitation and bubble growth during ultrasonic cleaning process, in: T. Sile, J. Virbulis, A. Timuhins, J. Sennikovs, U. Bethers (Eds.), Proc. of International Scientific Colloquium Modelling for Material Processing, September 16-17, Riga, Latvia, 2010, pp. 329-334.

9. E. V. Volkova, E. S. Nasibullaeva, N. A. Gumerov, Numerical simulations of soluble bubble dynamics in acoustic fields, in: Proc. of the ASME 2012 International Mechanical Engineering Congress & Exposition (IMECE 2012), Houston, Texas, USA, 2012, pp. 317-323, 1 CD ROM.

10. A. A. Gubaydullin, A. C. Gubkin, Behavior of bubbles in cluster with acoustic exposure [Povedeniye puzyr'kov v klastere pri akusticheskom vozdeystvii], Modern Science. Researches, Ideas, Results, Technologies (1 (12)) (2013) 363-367, in Russian.

11. E. V. Butyugina, E. S. Nasibullaeva, N. A. Gumerov, I. S. Akha-tov, Numerical simulation of gas microbubble dynamics in an acoustic field with account for rectified diffusion [Chislennoye modelirovaniye di-namiki gazovogo mikropuzyr'ka v akusticheskom pole s uchetom protsessa napravlennoy diffuzii], Vychislitelnaya mekhanika sploshnykh sred 7 (3) (2014) 234-244, in Russian.

12. R. P. Taleyarkhan, C. D. West, J. S. Cho, R. T. Lahey Jr., R. I. Nigmatulin, R. C. Block, Evidence for nuclear emissions during acoustic cavitation, Science 295 (2002) 1868-1873. doi:10.1126/science.1067589.

13. R. P. Taleyarkhan, C. D. West, R. T. Lahey Jr., R. I. Nigmatulin, R. C. Block, Y. Xu, Nuclear emissions during self-nucleated acoustic cavitation, Physical Review Letters 96 (2006). doi:10.1103/PhysRevLett.96.034301.

14. M. M. Martynyuk, N. Yu. Kravchenko, Nuclear fusion reaction in mat-aphase substance in the process of electrical explosion [Reaktsii yadernogo sinteza v mezofaznom veshchestve v protsesse elektricheskogo vzryva], Applied Physics (1) (2003) 79-90, in Russian.

15. R. I. Nigmatulin, I. S. Akhatov, A. S. Topolnikov, R. K. Bolotnova, N. K. Vakhitova, R. T. Lahey Jr., R. P. Taleyarkhan, The theory of supercompression of vapor bubbles and nano-scale thermonuclear fusion, Phys. Fluids 17 (10) (2005) 1-31, art. 107106. doi:10.1063/1.2104556.

16. M. M. Martynyuk, N. Yu. Kravchenko, Impact-cluster nuclear fusion. Conditions of excitation of the process, Bulletin of Peoples' Friendship University of Russia. Series: Physics (1) (2005) 118-128, in Russian.

17. Determination of temperature and pressure within the cavitational cavity [Opredeleniye temperatury i davleniya vnutri kavitatsionnoy polosti], in Russian.

18. M. M. Martynyuk, P. A. Tamanga, N. Yu. Kravchenko, The titanium phase diagram at the phase transition region liquid-vapor, Bulletin of Peoples' Friendship University of Russia. Series: Physics (10) (2002) 121-125, in Russian.

19. V. N. Khmelev, S. S. Khmelev, R. N. Golykh, G. A. Bobrova, O. N. Kra-sulja, V. I. Bogush, S. Anandan, Experimental determining of conditions of ultrasonic influence for providing maximum cavitation intensity in medium [Eksperimental'noye opredeleniye usloviy ul'trazvukovogo vozdeystviya dlya obespecheniya maksimal'noy intensivnosti kavitatsii v srede], Yuzhno-Sibirskiy nauchnyy vestnik (4 (12)) (2015) 50-55, in Russian.

20. A. A. Aganin, Dynamics of a small bubble in a compressible fluid, International Journal for Numerical Methods in Fluids 33 (2) (2000) 157-174. doi:10.1002/(SICI)1097-0363(20000530)33:2<157::AID-FLD6>3.0.C0;2-A.

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)