Formation of Radially Expanding Liquid Sheet by Impinging Two Round Jets
-
摘要: 采用两股互相冲击的圆射流可以形成环形的液体薄膜,液膜在径向扩展到一定的临界半径距离会破碎.数值模拟了液膜在周围气体中形成和破碎的非定常过程.考虑了液体和气体都是不可压缩Newton流体的轴对称问题.液体和气体的界面采用Level set函数来跟踪,Navier-Stokes 控制方程和物理边界条件采用有限差分格式离散求解.计算结果给出了环形液体薄膜形成并在其环形边缘处破碎,并缓慢运动的过程.液膜的厚度随着液膜在轴向的扩展会逐渐变薄,因此定义的局部Weber数会在径向逐渐减小,这里的局部Weber数定义为ρu2h/σ,其中ρ和σ分别为液体的密度和界面的张力,u和h分别为在径向某个位置的液膜的平均径向速度和半液膜厚度.数值结果表明就像实验中所观察到的那样,液膜径向扩展的过程的确会在局部Weber数趋向于1的时候终结而停止扩张.根据空间-时间线性稳定性理论,液膜的破碎最初是由正弦模式在临界局部Weber数Wec=1引起的,在临界局部Weber数小于1时会发生绝对不稳定性.在线性理论中另一个独立的模式,所谓的余弦模式,则增长比正弦模式要慢,从而会推测到正弦模式主导破碎的结论.然而,这里的数值结果却表明,余弦模式在界面波的非线性发展阶段实质的超越了正弦模式的增长,并对液膜的最终阶段的破碎起主导作用.这验证了线性理论只能够对触发时扰动波的性质进行预测,而对失稳后情况和结果的预测则不一定正确.
-
关键词:
- 流动稳定性 /
- Level set 方法 /
- 界面流动 /
- 液体射流
Abstract: A thin circular liquid sheet can be formed by impinging two identicalround jets against each other. The liquid sheet expands to a certain critical radial distance and breaks. The unsteady process of the formation and breakup of the liquid sheet in the ambient gas was smiulated num erically. Both liquid and gas were treated as incompressible Newtonian fluids. The flow considered was axi-symm etric. The liquid-gas interface was modeled with a level set function. A finite difference scheme was used to solve the governing Navier-Stokes equations with physical boundary conditions. The numerical results show how a thin circular sheet can be formed and broken at its circular edge, in slow motion. The sheet continues to thin as it expands radially. Hence the Weber number decreases rad ially. The Weber number is defined asu2h/, where and are respectively the liquid density and the surface tension, and u and h are, respectively, the average velocity and the half sheet thickness at a local radial location in the liquid sheet. The num erical results show that the sheet indeed terminates at a radial location where the Weber num ber reaches one as observed in experiments. The spatio-tem poral linear theory predicts that the breakup is initiated by the sinuousm ode at the critical Weber number Wec=1 be low which absolute instability occurs. The other independentmode called varicose mode grow smore slowly than the sinuousmode according to the linear theory. However our numerical results show that the varicose mode actually overtakes the sinuous mode during the non linear evolution, and is respon sible for the final breakup. The linear theory predicts the nature of disturbance waves correctly only at the onset of instability, but cannotpredict the exact consequence of the instability.-
Key words:
- flow instability /
- level setmethod /
- interfacial flow /
- liquid jet
-
[1] Savart F. Memoire sure le choc d’une veine liquid lancee contre un plan circulair[J]. Ann Chim Phys, 1833, 54: 56-87. [2] Savart F. Suit du memoire sur le choc dune veine liquid lancee contre un plan circulaire[J]. Ann Chim Phys, 1833, 54: 113-145. [3] Savart F. Memoire sure le choc de deux veine liquids, animees de mouvemens directement opposes[J]. Ann Chim Phys, 1833, 55: 257-310. [4] Chubb D L, Calfo F D, McConley M W, McMaster M S, Afjeh A A. Geometry of thin liquid sheet flows[J]. AIAA J, 1994, 32(6): 1325-1328. doi: 10.2514/3.12139 [5] Soderberg L D, Alfredsson P H. Experimental and theoretical stability investigations of plane liquid jets[J]. Eur J Mech B/Fluids, 1998, 17(5): 689-737. doi: 10.1016/S0997-7546(98)80022-8 [6] Squire H B. Investigation of the instability of a moving liquid film[J]. Br J Appl Phys, 1953, 4(6): 167-169. doi: 10.1088/0508-3443/4/6/302 [7] Hagerty W W, Shea J F. A study of the stability of plane fluid sheets[J]. J Appl Mech, 1955, 22(4): 509-514. [8] Fraser R P, Eisenklam P, Dombrowski N, Hasson D. Drop formation from rapidly moving liquid sheets[J]. AIChE J, 1962, 8(5): 672-680. doi: 10.1002/aic.690080522 [9] Taylor G I. The dynamics of thin sheets of fluid Ⅱ: waves on fluid sheets [J]. Proc R Soc Lond, Ser A, 1959, 253: 296-312. doi: 10.1098/rspa.1959.0195 [10] Taylor G I. The dynamics of thin sheets of fluid Ⅲ: disintegration of fluid sheets[J]. Proc R Soc Lond, Ser A, 1959, 253: 313-321. doi: 10.1098/rspa.1959.0196 [11] Huang J C P. The break-up of axisymmetric liquid sheets[J]. J Fluid Mech, 1970, 43(2): 305-319. doi: 10.1017/S0022112070002392 [12] Clanet C, Villermaux E. Life of a smooth liquid sheet[J]. J Fluid Mech, 2002, 462: 307-340. [13] Villermaux E, Clanet C. Life of a flapping liquid sheet[J]. J Fluid Mech, 2002, 462: 341-363. [14] Crapper G D, Dombrowski N, Pyott G A D. Large amplitude Kelvin-Helmholtz waves on thin liquid sheets[J]. Proc R Soc Lond, Ser A, 1975, 342: 209-224. doi: 10.1098/rspa.1975.0021 [15] Clark C J, Dombrowski N. Aerodynamic instability and disintegration of inviscid liquid sheets[J]. Proc R Soc Lond, Ser A, 1972, 329: 467-478. doi: 10.1098/rspa.1972.0124 [16] Crapper G D, Dombrowski N, Jepson W P. Wave growth on thin sheets of non-Newtonian liquids[J]. Proc R Soc Lond, Ser A, 1975, 342: 225-236. doi: 10.1098/rspa.1975.0022 [17] Lin S P, Jiang W Y. Absolute and convective instability of a radially expanding liquid sheet[J]. Phys Fluids, 2003, 15(6): 1745-1754. doi: 10.1063/1.1570422 [18] Brackbill J U, Kothe D B, Zemach C. A continuum method for modeling surface tension[J]. J Comp Phys, 1992, 100(2): 335-354. doi: 10.1016/0021-9991(92)90240-Y [19] Osher S, Sethian J A. Fronts propagating with curvature dependent speed: algorithms based on Hamilton-Jacobi formulations[J]. J Comp Phys, 1988, 79(1): 12-49. doi: 10.1016/0021-9991(88)90002-2 [20] Sethian J A. Level Set Methods and Fast Marching Methods[M]. London: Cambridge University Press, 1999. [21] Sussman M, Smereka P, Osher S. A level set approach for computing solutions to incompressible 2-phase flow[J]. J Comp Phys, 1994, 114(1): 146-159. doi: 10.1006/jcph.1994.1155 [22] Chorin A J. Numerical solution of the Navier-Stokes equations[J]. Math Comp, 1968, 22(104): 745-762. doi: 10.1090/S0025-5718-1968-0242392-2 [23] Temam R. On an approximate solution of the Navier-Stokes equations by the method of fractional steps, part 1[J]. Archiv Ration Mech Anal, 1969, 32(2): 135-153. [24] Fortin M, Peyret R, Temam R. Résolution numérique des équations de Navier-Stokes pour un fluide visqueux incompressible[J]. J Mech, 1971,10(3): 357 -390. [25] Li J. Cacul d’interface affine par morceaux(piecewise linear interface calculation[J]. C R Acad Sci Paris, Ser ⅡB, 1995, 320(8) : 391-396. [26] Jiang G S, Shu C W. Efficient implementation of weighted ENO schemes[J]. J Comp Phys, 1996, 126(1): 202-228. doi: 10.1006/jcph.1996.0130
点击查看大图
计量
- 文章访问数: 1404
- HTML全文浏览量: 134
- PDF下载量: 867
- 被引次数: 0