CC BY
11
Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki
[J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2020, vol. 24, no. 2, pp. 390-400 ISSN: 2310-7081 (online), 1991-8615 (print) d https://doi.org/10.14498/vsgtu1763
MSC: 76T10
An undamped oscillation model with two different contact angles for a spherical droplet impacting on solid surface
Shi Chen, Bozhong Cong, Dongqi Zhang, Xiaohui Liu, Shengqiang Shen
Dalian University of Technology, Dalian, 116024, China.
Abstract
In order to further elucidate the dynamic theory of droplet oscillating on solid surface, a new handling method of contact angle of the droplet during the process of the oscillation was founded, which is based on the spherical model. The influence of gravity on the contact angle and spreading radius was discussed. Thus, an equation between the spreading radius of the droplet and time flow was founded. The results of theoretical calculation were compared with smoothed numerical results.
Keywords: droplet oscillation, contact angle, theoretical analysis, spectrum analysis.
Received: 2nd January, 2020 / Revised: 12th March, 2020 / Accepted: 16th March, 2020 / First online: 30th June, 2020
Short Communication
3 ©® The content is published under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/) Please cite this article in press as:
Chen S., Cong B., Zhang D., Liu X., Shen S. An undamped oscillation model with two different contact angles for a spherical droplet impacting on solid surface, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2020, vol. 24, no. 2, pp. 390-400. doi: 10.14498/vsgtu1763. Authors' Details:
Shi Chen; Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education; e-mail: dlthermo@dlut.edu .cn
Bozhong Cong; Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education; e-mail: ttcbz@mail.dlut. edu.cn
Dongqi Zhang; Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education
Xiaohui Liu A; Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education; e-mail: lxh723@dlut.edu.cn
Shengqiang Shen; Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education; e-mail: zzbshen@dlut.edu.cn
1. Introduction
As a common phenomenon in nature, droplets impacting on solid surfaces exists in many fields, such as engineering, chemical industry, agriculture, aerospace and energy. For examples, the liquid fuel injection in the internal combustion engine [1], the ink jet of the printer [2], the process of seawater evaporation in desalination system [3], the spraying of the refrigerant in the evaporator in the refrigeration system [4], etc. The study of relationships among the dynamic parameters, initial parameters of droplets and initial can be beneficial to predict the processes of droplets impacting on solid flats.
For a droplet impacting on solid surfaces, according to the initial parameters, there may be oscillations, bounces, splashes, etc. Compared with other circumstances, the oscillations have more value and possibility to be studied. Since T. Young [5] presented that the contact angle between liquid and solid surface was constant, and P.S. Laplace [6] inferred the relationship between additional pressure and radius of curvature, many studies have been done by experiments, numerical simulations and theoretical analysis.
In the aspect of experimental studies, M. Marengo et al. [7] studied the effects of the impact parameters on the droplet impingement, and found that the surface wettability had a strong influence on the spreading of droplet in the later stages of the process. I.S. Bayer [8] studied the dynamic characters of contact angle between smooth surface and droplets with different wettability. M. Remer et al. [9] studied the variation of droplet dynamical contact angle after impacting on three kinds of coating surfaces.
In the aspect of numerical simulation, S.F. Lunkad et al. [10] simulated the drop impact and spreading process on horizontal and inclined surfaces using the volume of fluid (VOF) method, and investigated the effects of surface inclination, surface wetting characteristics, liquid properties and impact velocity on the droplet oscillation by using static contact angle (SCA) and dynamic contact angle (DCA) models. Y. Yao et al. [11] analyzed droplets oscillation with VOF method and introduced a model of dynamical contact angle to improve the accuracy of simulation. S. Sikalo et al. [12] carried out the numerical simulations of a single drop impacting onto a dry, partially wettable substratum.
In the aspect of theoretical analysis, S. Vafaei et al. [13,14] explained the dependence of contact angle on the size of liquid droplets on smooth solid substrates, and demonstrated that for sessile droplets on smooth surfaces, the contact angle can be uniquely determined for given droplet mass (or volume) and liquid/solid/gas properties. I.V. Roisman et al. [15] studied the normal impact of a liquid drop on a dry solid surface theoretically, and introduced a strictly theoretical model to predict the evolution of the drop diameter.
The purpose of this paper is to establish an amended theoretical model of undamped droplet oscillation.Compared with original model, in which droplet contact angle keeps constant value, the droplet contact angle in new model varies with spreading radius and droplet volume. The waveform of droplet spreading radius obtained from new model, original model and numerical simulation were compared.
2. Theoretical Model
2.1. The Oscillating Equation of Spherical Segment Droplet with Changing Contact Angle
For a droplet oscillating on a horizontal solid surface, when the size of droplet is short enough that the influence of gravity on the droplet shape can be ignored. The shape of the droplet can be considered as a sphere, shown in Fig. 1.
Figure 1. The geometrical model of droplet
At a certain moment of oscillation, spreading radius of droplet in horizontal direction is r, the radius of the sphere is R, the height is z and the contact angle between the droplet and the solid surface is . According to the geometrical relationship, the equation between r, R, and d can be written as:
3 V
—3 = (1 - cos d)2(2 + cos 0), (1)
nr 3
R = (2) * = r1-^, (3)
where V is the droplet volume.
When a droplet impacts on a solid surface, the shape of droplet is strongly irregular and hard to describe in beginning several oscillation periods after impact. But the oscillation will be gradually stabilized with close periods and amplitudes. Thus, it can be assumed that the droplet keeps the shape of spherical cap in the hereafter oscillation.
In previous study [16], the oscillation of droplet could be seen as the result of interaction of the surface tension, the internal pressure and the inertial force. A differential element is selected in the segment droplet of the angle, shown in Fig. 2.
And the oscillating equation of spherical segment droplet can be written as [16]:
ipRA0 ip(R — z) , 2 o. ' -- (2 r2 + 3 R2)
d2r M2
4 r 12
AlT(7 A / r3 \
— 2ia(l — v cos (p) +--—— +2npg[RA — — ) =0, (4)
R 3
<7 COS 6
COS 6
Figure 2. The force analysis of droplet differential. The surface tension and the internal pressure equally and vertically act on each vertical section of the droplet. The surface tension points to the inside of differential element and the internal pressure points to the outside of differential
element
where p — the intrinsic contact angle of the droplet; 9 — the dynamic contact angle of the droplet;1 a — the surface tension coefficient;
I = RA (5)
— arc length of the great circle of droplet;
A = 1 [Rl - r(R - z)] (6)
— the sectional area of the droplet in the vertical plane. Eq. (4) can be also written as:
rl2v
pw + Q + W = 0, (7)
P = - ^ + "ÉS-A (2,2 + 3R2),
4 r 12
^ ,, . 4" a A
Q = —2Tra(i — rcosp) +---—,
R
3,
W = 2irpg(^RA - y) ,
where P — the factor of inertia, Q — the factor of surface tension and internal pressure, W — the factor of gravity.
When substituted equations of geometrical relationship (Eq. (1)-(3), and Eq. (5), (6)), P, Q, and W can be written as:
3r cos 9 cos 9 9
P = "pr3 -- +-5---j- , (8
F V6sin0 4 sin3 9 4 sin4 9/ v
xThe angle ip is determined by the involved surface energies. When the gravity is absent and the droplet is static on the horizontal solid surface, ip = 0. The angle 0 can be influenced by inertial force and gravity of the droplet.
Q = 2n—r (cos p — cos (
W = npgt
f 9
cos 2
sin3 sin2 3
(9) (10)
In this model, the intrinsic contact angle is a constant valve, which determines direction of the contact force from the solid surface. And the dynamic contact angle determines the shape of droplet. The geometrical relationship between them is shown in Fig. 3.
Figure 3. The geometrical relationship between the intrinsic contact angle and the dynamic
contact angle
The deviation between the intrinsic contact angle and the dynamic contact angle increases with the increase of acceleration of the droplet. And the expression of the contact force from the solid surface in Fig. 2, "—cos d" can be replaced by "—cos (".
2.2. The Influence of Gravity
When the spreading radius reaches the balance value (or the droplet keeps static):
d2
pd2 + Q + W = °-
Substituted Eq. (8)-(10) and multiplied:
cos p = cos d —
cos
pgr2
2— Vsln3^ sln2^
Obviously, in the absence of gravity, d = p. In normal gravity condition that the droplet can still keep the shape of spherical segment (called micro gravity in followings), the dynamic contact angle d can be calculated by iteration from Eq. (1) and Eq. (5). The influence of gravity is shown in Fig. 4.
With the increase of droplet volume, the deviation between the micro gravity solutions and the absence gravity solutions on dynamic contact angle and spreading radius becomes larger, but that of spreading radius is quite small.
Droplet Volume, m3 xlO-9
Figure 4. (top) the influence of droplet volume on contact angle 0 (here, the gravity of droplet depends on the droplet volume); (bottom) the influence of gravity on balance radius
Figure 5. The frequency spectrum of numerical data
3. Numerical Simulation
3.1. Original Case
For the numerical solution,the processes of drop oscillation are simulated by using the VOF model of Fluent. The body force weighted model is used to calculate of pressure and gravity balance. The PISO algorithm is used to couple the droplet velocity and pressure in the momentum equation. The slip boundary condition is used and the shear stress is 0 Pa. The viscosity of water is 1.003 • 10-5 Pa ■ s which is multiplied by 0.01 in order to simulate the undamped oscillation of droplet. The time step is 4 ■ 10-6 s, and the residual error is 1 ■ 10-5. Others keep the default algorithm.
3.2. Spectral Analysis
The data from original case contains much ingredient of noise. This is because the oscillation is influenced by many minimal factors. After eliminating several data points with extreme deviation, the data was resampled by FFT operation. The spectral analysis of the original data is shown in Fig. 5.
In Fig. 5, the main frequency of the oscillation is concentrate upon [250 Hz, 550 Hz]. The corresponding frequency value of the peak value in figure is 426 Hz, which means the periods of the oscillation is about 2.35 ms with the same initial parameters shown in Tab. 1.
Table 1
The initial parameters of theoretical and numerical calculation
Initial Speed Droplet Volume Gravity Density Surface Tension
0.8 m/s 1.1310 • 10-10 m3 9.8 m/s 1003 kg/m3 0.073 N/m
3.3. Data Smoothing
To eliminate the noise, a Butterworth filter was applied. The response type is bandpass with the pass band [400 Hz, 500 Hz], the fluctuate is less than 1 dB. On the both side of pass band, the signal is decreased to 10 dB. The frequency spectrum after smoothing is shown in Fig. 6.
After smoothing, the data can be regarded as a superposition of several sine curves, with close frequency, amplitude and phase position.
4. Comparison
Eq. (7) is solved by using the fourth-order Runge-Kutta method in MATLAB R2016a.
The relationship between spreading radius and time flow of theoretical and numerical results (after smoothing) is shown in Fig. 7, with the same initial parameter shown in Tab. 1.
The period of theoretical results is 2.37 ms, which is almost equal to the average periods of numerical results. The amplitude of theoretical results is about 10 % larger than the average value of the numerical results. Compared with constant contact angle model [16], the average periods were more close to the numerical data (1.8 ms). The relationship between spreading radius and its accelerated velocity is shown in Fig. 8.
When ( is approaching or higher than the balance value (0), the deviation between theoretical and numerical results is relatively small, but when is far
Figure 6. The frequency spectrum after smoothing
Figure 7. The relationship between spreading radius and time flow (£ = (r — rb)/rb, rb is the
theoretical balance radius)
Figure 8. The relationship between spreading radius and its accelerated velocity
below the balance value, the deviation is considerable. This is because an obvious distortion would happen on the shape of droplet and the above-mentioned relationships (Eq. (8)-(10)) would lose efficacy. Besides, because the relationship between spreading radius and its acceleration in present model is not linear, the oscillation period is related to the initial velocity.
5. Discussion
The process of undamped oscillation for a droplet impacting on solid surface was theoretically described by an equation based on spherical segment model. Followings are several conclusions.
The contact of the droplet is changing in the oscillation and geometrically calculated. The influence of gravity on contact angle and spreading radius was analyzed quantitatively, the deviation was lower than the constant contact angle model. However, the result was based on the geometrical model of spherical segment of droplet, which might be inaccurate under the influence of gravity. More study was needed in this content.
After spectrum analyzing and smoothing by a filter, the numerical solutions can be regarded as a superposition of several sine curves with close frequency, amplitude and phase positions. The relationship between spreading radius and its accelerated velocity was founded by polynomial fit, and the curve was symmetrical about the balance point and strongly linear.
By comparing the results of theoretical and numerical calculation, when the spreading radius of the droplet approaches the balance value, the deviation between two solutions are small, but when the spreading radius is far from the balance value, the deviation can be considerable. And the deviation of far below from the balance value is larger than that of far above from the balance value. This is possibly because when the spreading radius is far from the balance value, the shape of the droplet would also deviate from spherical segment. More work about the amendment of the shape of the droplet in oscillation is needed in later study.
Competing interests. We declare that we have no conflicts of interest in the authorship or publication of this paper.
Authors' contributions and responsibilities. Each author has participated in the article concept development and in the manuscript writing. The authors are absolutely responsible for submitting the final manuscript in print. Each author has approved the final version of manuscript.
Funding. This research was funded by National Natural Science Foundation of China (grant no. 51936002).
References
1. Nakayama Y., Kidokoro T., Sakurai K. Fuel injection control system of an internal combustion engine, US Patent no. US9169758B2, 2015, https://patents.google.com/patent/ US9169758B2/.
2. Slater S.D., Clippingdale A.J., Newcombe G.C.F. Printing process and liquid ink jet ink, US Patent no. US9156256B2, 2015, https://patents.google.com/patent/US9156256B2/.
3. Qi C.H., Feng H.J., Lv H.Q., Miao C. Numerical and experimental research on the heat transfer of seawater desalination with liquid film outside elliptical tube, Int. J. Heat Mass Transfer, 2016, vol.93, pp. 207-216. doi: 10.1016/j.ijheatmasstransfer.2015.09.043.
4. Hartfield J.P., Sanborn D.F. Falling film evaporator with refrigerant distribution system, Canada Patent no. CA2219676A1, 1995, https://patents.google.com/patent/ CA2219676A1/.
5. Young T. An essay on the cohesion of fluids, Phil. Trans. Roy. Soc. London, 1805, vol. 95, pp. 65-87. doi:10.1098/rstl.1805.0005.
6. Laplace P.S. Sur l'action capillaire. Supplément à la théorie de l'action capillaire, In: Traité de mécanique céleste, vol. 4, Supplement 1, Livre X. Paris, Gauthier-Villars et fils, 1805, pp. 771-777.
7. Sikalo S., Marengo M , Tropea C., Ganic E.N. Analysis of impact of droplets on horizontal surfaces, Experimental Thermal and Fluid Science, 2002, vol. 25, no. 7, pp. 503-510. doi: 10. 1016/S0894-1777(01)00109-1.
8. Bayer I. S., Megaridis C. M. Contact angle dynamics in droplets impacting on flat surfaces with different wetting characteristics, J. Fluid Mechanics, 2006, vol. 558, pp. 415-449. doi: 10.1017/S0022112006000231.
9. Remer M., Psarski M., Gumowski K., Rokicki J., Sobieraj G., Kaliush M., Pawlak D., Celichowski G. Dynamic water contact angle during initial phases of droplet impingement, Colloids and Surfaces A: Physicochemical and Engineering Aspects, 2016, vol.508, pp. 5769. doi: 10.1016/j.colsurfa.2016.08.028.
10. Lunkad S. F., Buwa V. V., Nigam K.D.P. Numerical simulations of drop impact and spreading on horizontal and inclined surfaces, Chem. Eng. Sci., 2007, vol. 62, no. 24, pp. 7214-7224. doi: 10.1016/j.ces.2007.07.036.
11. Yao Y., Meng S., Li C., Chen X., Yang R. Droplet oscillation after impact on a solid surface, International Mechanical Engineering Congress and Exposition, 2016, vol. 7, Fluids Engineering, IMECE2016-66025. doi: 10.1115/IMECE2016-66025.
12. Sikalo S. , Wilhelm H.-D., Roisman I. V. , Jakirlic S., Tropea C. Dynamic contact angle of spreading droplets: Experiments and simulations, Phys. Fluids, 2005, vol. 17, no. 6, 062103. doi: 10.1063/1.1928828.
13. Vafaei S., Podowski M. Z. Theoretical analysis on the effect of liquid droplet geometry on contact angle, Nuclear Eng. Design, 2005, vol. 235, no. 10-12, pp. 1293-1301. doi: 10.1016/ j.nucengdes.2005.02.026.
14. Vafaei S., Podowski M. Z. Analysis of the relationship between liquid droplet size and contact angle, Adv. Colloid Interface Sci., 2005, vol. 113, no. 2-3, pp. 133-146. doi: 10.1016/j.cis. 2005.03.001.
15. Roisman I. V., Rioboo R., Tropea C. Normal impact of a liquid drop on a dry surface: model for spreading and receding, Proc. Royal. Soc. A, 2002, vol.458, pp. 1411-1430. doi: 10. 1098/rspa.2001.0923.
16. Chen S., Zhang D., Shen S., Liu X., Chen Y. Spherical drop impact on solid surfaces: Undamped oscillation theoretical model, AIP Conf. Proc., 2018, vol.1984, no. 1, 020032. doi:10.1063/1.5046616.
Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки. 2020. Т. 24, № 2. С. 390-400 ISSN: 2310-7081 (online), 1991-8615 (print) d https://doi.org/10.14498/vsgtu1763
УДК 517.958:531.32
Модель незатухающего колебания для сферической капли на твердой поверхности с двумя различными углами контакта
Shi Chen, Bozhong Cong, Dongqi Zhang,
Xiaohui Liu, Shengqiang Shen
Dalian University of Technology.
Dalian, 116024. China.
Аннотация
Предложен нелинейный подход описания колебания сферической капли на твердой поверхности. Интегрирование уравнений движений осуществляется без использования линеаризации тригонометрических функций, зависящих от угла контакта. Иными словами, угол контакта является произвольной конечной величиной. Проведено исследование влияние силы тяжести на угол контакта и радиус распространения капли по твердой поверхности. Таким образом, было найдено нелинейное уравнение, описывающее изменение радиуса распространения капли в зависимости от времени. Данное уравнение было численно проинтегрировано. Исследование численной сходимости осуществлялось посредством сравнения с известными модельными точными решениями и известными экспериментальными данными. На основании исследования методами численного интегрирования полученного в статье уравнения можно сделать вывод о целесообразности использования математической модели для описания и исследования новых физических эффектов при колебании капель.
Ключевые слова: колебание капли, угол контакта, теоретический анализ, спектральный анализ.
Получение: 2 января 2020 г. / Исправление: 12 марта 2020 г. / Принятие: 16 марта 2020 г. / Публикация онлайн: 30 июня 2020 г.
Краткое сообщение
3 ©® Контент публикуется на условиях лицензии Creative Commons Attribution 4.0 International (https://creativecommons.org/licenses/by/4.0/deed.ru) Образец для цитирования
Chen S., Cong B., Zhang D., Liu X., Shen S. An undamped oscillation model with two different contact angles for a spherical droplet impacting on solid surface, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2020, vol. 24, no. 2, pp. 390-400. doi: 10.14498/vsgtu1763.
Сведения об авторах
Shi Chen; Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education; e-mail: dlthermo@dlut.edu. cn
Bozhong Cong; Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education; e-mail: ttcbz@mail.dlut. edu.cn
Dongqi Zhang; Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education
Xiaohui Liu A; Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education; e-mail: lxh723@dlut.edu.cn
Shengqiang Shen; Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education; e-mail: zzbshen@dlut.edu.cn
