Discrete and Continuous Models and Applied Computational ScienceDiscrete and Continuous Models and Applied Computational Science2658-46702658-7149Peoples' Friendship University of Russia named after Patrice Lumumba (RUDN University)2219610.22363/2658-4670-2019-27-1-49-59Research ArticleMathematical model of cavitation under the influence of a single stretching pulseKravchenkoNikolay Yu<p>Deputy Director of Institute of Physical Research and Technology</p>kravchenko-nyu@rudn.ruKulyabovDmitry S<p>Doctor of Sciences in Physics and Mathematics, Full Professor, Department of Applied Probability and Informatics</p>kulyabov-ds@rudn.ruPeoples’ Friendship University of Russia (RUDN University)15122019271495920112019Copyright © 2019, Kravchenko N.Y., Kulyabov D.S.2019<p>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.</p>numerical simulationcavitationmathematical model of cavitationчисленное моделирование, кавитация, математическая модель кавитации<p>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. Cavitation equation The equation describing the dynamics of cavitation bubbles in an incom- pressible fluid without taking into account the vapor pressure, surface tension and viscosity of the fluid can be represented as: г1 RR + 3 R˙2 = 1 p 2 \ 1 Ro + \3k 2 - - p l - p(t) . (1) o 2 b R R R 0 Here: R0 is the radius of the nucleus at t = 0; R - radius of the nucleus at the next time instant t; - density of a liquid; 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; 2/R0 is the Laplace pressure; R0/R is the amplitude of oscillations of the cavity. Dynamics of a cavity under the action of single pulses of negative pressure A single impulse is presented in the form: p(t) = - (t + t1) exp 1 t + t1 \ - , (2) where t1 is the time of appearance of the first germ of homogeneous cavitation under the action of this pulse. In this case dp 1 1 \ - = dt (t + t1) exp 1 t + t1 \ + exp 1 t + t1 \ . (3) At t = 0, the rate of increase in pressure at the initial moment of time is dp = - dt . From the condition dp = 0, applied to (3) we find tm = , pm = dt , e where 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 1 = 0.95e(pb - ps). (5) Thus, the rate of pressure increase can be calculated for a given fluid temperature using the value of (pb - ps) for a given . The time of appearance of the first nucleus t1 is found from the condition (2), when at t = t1 the pressure is p = p1. p1 = -t1 exp 1 t1 \ - = -(pb - ps). (6) Here p1 is the pressure at which one germ appears in 1 cm3 at the moment of time t1, p = - b Substituting into (6), we get p1 - ps 1. (7) where x = t1/. = e t1 exp 1 t1 \ - = e x exp (-x) (8) Thus, in the equation (3), for a given value of , the value of ; is calculated by the formula (5) and the value of t1 is calculated by the formula (2) for a given value of (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 - 163 l , 3(1 - v /l)2(pb - p)2 V = 1 sm3, 163 A = 3(1 - v /l)2 pb - p = t, , (9) ( dN = B exp A - 2t2 \ dt, N = 1 = B ( t1 ( exp 0 A - 2t2 \ dt. (10) By calculating 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 A C = - 2 , t1 ( 1 C \ B exp t2 0 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: RR + 3 R˙2 = 2 1 г1 = pb + 2 \ 1 Ro \3k 2 - - pb - (t + t1) exp 1 t + t1 \l - . (12) Ro R R The differential equation for the growth of an embryo vapor has the form: 1 l RR + 3 R˙ 2\ 2 = -p(t) - 2 1 - 1 R l \-1 v , (13) where p(t) is defined by the formula (2). Substituting p(t) from (2) into (13), we get: 1 l RR + 3 R˙ 2\ 2 = (t + t1) exp 1 t + t1 \ - - 2 1 - 1 R l \-1 v . (14) For the reasons described above, we will assume t1 0, and then the equation (14) takes the form: 1 l RR + 3 R˙ 2\ 2 = t exp 1 t \ - - 2 1 - 1 R l \-1 v . (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 = 1 г 1 t \ t exp 2 1 1 l \ -1l 3 R˙ 2 (16) lR or - - R v 2 Rl du(1) = 1 г lu(2) t exp 1 t \ - - 2 1 - 1 u(2) l \-1l v - 3 [u(1)]2 2 u(2)l . (17) The last equation is solved numerically by the Runge-Kutta method. The calculations were carried out for = 10-9 s. 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, pl = pl0, Tl = Tl0, Nb = 0, u = 0. The rate of pressure drop p at the point with coordinate x = 0 was set by the boundary conditions x = 0, pl = pl0 - p t, u = 0. x After a bubble appeared, the following initial conditions were accepted for ordinary differential equations describing its development (12): R = Rc, pv = pb, tv = Tl. 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): R0 = omr0, 1 (R0)2 = omr02, 2 R0 Pb + = sr0, Pb = pb, 2 \ R = pbsr0, Pa = pa, 0 (R0)2 = rhomr, ϕo = 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) = - + 2 omr0 u(2) 1 1 Г pbsr0 + 3k 1 - sr0 l - pb + pa sin( + fio) , rhomr u(2) (u(2)) u(2) du(2) = u(1). 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]. 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. 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 t1 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, pbsr0, pk, sr0, pa, pb, pi, nk, c fa, r0, p, p1, 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 c bubble with acoustic effects and the frequency of Fa c according to Runge Kutta method.' print *, '-------------------------------------------' print *, 'time pressure dependence P = P (t / T) is c written to the file "p.dat", the time dependence of c the radius R / R0 = f (t / T) is written to the file c "r.dat".' print *, 'THE PROGRAM WILL NOT WORK without an auxiliary c file "bubble.txt" and also without the additional c program "rkgs.for"!' read (3, *) q1, q2, fa, pa, pb, p1, r0, sig, amu, ro, kt, c ky pi = 3.141592654 pk = 5 print *, 'The following parameters are entered here:' print *, '1) Time step integration (wt): pt (3) =', q1 print *, '2) Integration error: pt (4) =', q2 print *, '3) Acoustic frequency: Fa =', fa, 'Hz' print *, '4) Pressure amplitude at a given point: Pa =', c pa, 'Pa' print *, '5) Pressure on binodals at a given point: Pb =', c pb, 'Pa' print *, '6) Pressure at the point of emergence of the c unit: P1 =', p1, 'Pa' print *, '7) Starting radius of the nucleus: R0 =', r0, c '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 c =', 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 c the file' bubble.txt '' print *, 'To interrupt the program, press "Ctrl + C", for c 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 continue print *,' Error code ',ih goto 5 continue print *,' ATTENTION !!!' print *, 'ERROR CODE SHOULD NOT EXCEED 10.' print *, 'Correct, please, the parameters pt(3),pt(4)' read(*,*)q1,q2 goto 6 continue stop end Below are some results of a numerical simulation of the bubble behavior in a rarefaction wave. 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.</p>[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.][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/uzku1365][M. M. Martynyuk, N. Yu. Kravchenko, Limit of thermodynamic stability of a liquid phase in the field of nefative pressure [Granitsa termodinamicheskoy ustoychivosti zhidkoy fazy v oblasti otritsatel’nykh davleniy], Russian Journal of Physical Chemistry 72 (6) (1998) 998-1001, in Russian.][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.][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/jmcms.2019.03.00046.][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.][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.][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.][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.][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.][E. V. Butyugina, E. S. Nasibullaeva, N. A. Gumerov, I. S. Akhatov, Numerical simulation of gas microbubble dynamics in an acoustic field with account for rectified diffusion [Chislennoye modelirovaniye dinamiki gazovogo mikropuzyr’ka v akusticheskom pole s uchetom protsessa napravlennoy diffuzii], Vychislitelnaya mekhanika sploshnykh sred 7 (3) (2014) 234-244, in Russian.][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.][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.][M. M. Martynyuk, N. Yu. Kravchenko, Nuclear fusion reaction in mataphase 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.][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.][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.][Determination of temperature and pressure within the cavitational cavity [Opredeleniye temperatury i davleniya vnutri kavitatsionnoy polosti], in Russian.][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.][V. N. Khmelev, S. S. Khmelev, R. N. Golykh, G. A. Bobrova, O. N. Krasulja, 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.][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.CO;2-A.]