PARALLEL ADAPTIVE SIMULATION OF A PLUNGING LIQUID JET

doi:10.1016/S0252-9602(10)60060-4

Abstract

Abstract:

This paper is concerned with three-dimensional numerical simulation of a plunging liquid jet. The transient processes of forming an air cavity around the jet, capturing an initially large air bubble, and the break-up of this large toroidal-shaped bubble into smaller bubbles were analyzed. A stabilized finite element method (FEM) was employed under parallel numerical simulations based on adaptive, unstructured grid and coupled with a level-set method to track the interface between air and liquid. These simulations show that the inertia of the liquid jet initially depresses the pool's surface, forming an annular air cavity which surrounds the liquid jet. A toroidal liquid eddy which is subsequently formed in the liquid pool results in air cavity collapse, and in turn entrains air into the liquid pool from the unstable annular air gap region around the liquid jet.

Key words: plunging liquid jet, air entrainment, two-phase flows, level set method, parallel adaptive simulation

CLC Number:

76B10

Azat Yu. Galimov, Onkar Sahni, Richard T. Lahey Jr., Donald A. Drew, Kenneth E. Jansen. PARALLEL ADAPTIVE SIMULATION OF A PLUNGING LIQUID JET[J].Acta mathematica scientia,Series B, 2010, 30(2): 522-538.

Trendmd

References

[1] McCurdy E. The Notebooks of Leonardo Da Vinci, Vol 2. six-th ed. London, UK: Jonathan Cape, 1956

[2] Rayleigh L. On the instability of jets. Proc Lond Math Soc, 1878, 10: 4--13

[3] Rayleigh L. On the capillary phenomena of jets. Proc Roy Math Soc London A, 1879, 29: 71--97

[4] Rayleigh L. On the instability of a cylinder of viscous liquid under capillary force. Phil Mag, 1892, 34: 145--154

[5] Ehrenberger R. Wasserbewegung in steilen Rinnen (Susstennen) mit besonderer Berucksichtigung der SelSelbstbeluftung (Flow of water in Steep Chutes with Special reference to Self-aeration). Zeittschrift des Osterreichischer Ingebieur und Architektverein. No 15/16 and 17/18 (in German), 1926

[6] Weber C. Zum Zerfall eines Flussigkeitsstrahles. Z Angew Math Mech, 1931, 11: 136--154 (in German)

[7] Tomotika S. On the instability of a cylindrical thread of a viscous liquid surrounded by another viscous fluid. Proc Roy Math Soc London A, 1935, 150: 322--337

[8] McKeogh E J, Ervine D A. Air entrainment rate and diffusion pattern of plunging liquid jets. Chem Eng Sci, 1981, 36: 1161--1172

[9] Bonetto F, Lahey R T Jr. An experimental study on air carryunder due to a plunging liquid jet. Int J Multiphase Flow 1993, 19(2): 281--294

[10] Chanson H. Air Bubble Entrainment in Free-Surface Turbulent Shear Flows. Academic Press, 1996

[11] Davoust L, Achard J L, Hammoumi M El. Air entrainment by a plunging jet: the dynamical roughness concept and its estimation by a light absorbtion technique. Int J Multiphase Flow, 2002, 28(9): 1541--1564

[12] Chanson, H., Manasseh R. Air entrainment processes in a circular plunging jet: Void-fraction and acoustic measurements. J Fluids Eng, 2003, 125(5): 910--921

[13] Bonetto F, Drew D A, Lahey R T Jr. The analysis of a plunging liquid jet-the air entrainment process. Chem Eng Commun, 1994, 130: 11--29

[14] Zhu Y, Oguz H N, Prosperetti A. On the mechanism of air entrainment by liquid jets at a free surface. J Fluid Mech, 2000, 404: 151--177

[15] Lorenceau E, Reyssat E, Restango F, Eggers J, Quere D. On liquid/liquid impect//Proceeding of 6th Int Conf on Multiphase Flow (ICMF 2007). Leipzig, Germany, July 9-13, 2007

[16] Jansen K E, Whiting C, Hulbert G M. A generalized-alpha method for integrating the filtered Navier--Stokes equations with a stabilized finite element method. Comput Methods Appl Mech Eng, 2000, 190: 305--319

[17] Muller J, Sahni O, Li X, Jansen K E, Shephard M S, Taylor C A. Anisotropic adaptive finite element method for modeling blood flow. Comput Methods in Biomech Biomed Eng, 2005, 8(5): 295--305

[18] Whiting C H, Jansen K E. A stabilized finite element method for the incompressible Navier--Stokes equations using hierarchical basis.
Int J Numer Methods Fluids, 2001, 35: 93--116

[19] Nagrath S, Jansen K E, Lahey R T Jr. Computation of incompressible bubble dynamics with a stabilized finite element level set method. Comp Methods Applied Mech Eng, 2005, 194(42/44): 4565--4587

[20] Nagrath S, Jansen K, Lahey R T Jr, Akhatov I. Hydrodynamic simulation of air bubble implosion using a level set approach. J Comput Phys, 2006, 215(1): 98--132

[21] Osher S, Fedkiw R. Level Set Methods and Dynamic Implicit Surfaces. New York: Springer-Verlag, 2003

[22] Brackbill J U, Kothe D B, Zemach C. A continuum method for modeling surface tension. J Comput Phys, 1992, 100: 335--354

[23] Sussman M, Almgren A S, Bell J B, Howell L H, Colella P, Welcome W L. An adaptive level set approach for incompressible two-phase flows. J Comput Phys, 1999, 148: 81--124

[24] Taylor C A, Hughes T J R, Zarins C K. Finite element modeling of blood flow in arteries. Comp Meth Appl Mech Eng, 1998, 158: 155--196

[25] www.simmetrix.com.

[26] De Cougny H L, Shephard M S. Parallel refinement and coarsening of tetrahedral meshes. Int J Numer Meth Eng, 1999, 46: 1101--1125

[27] Sahni O, Zhou M, Shephard M S, Jansen K E. Scalable implicit finite element solver for massively parallel processing with demonstration to 160 K cores//Proc 2009 ACM/IEEE Conference on Supercomputing. Portland, OR, November, 2009

[28] Ainsworth M, Oden J T. A Posteriori Error Estimation in Finite Element Analysis. New York: John Wiley \& Sons, 2000

[29] Saad Y, Schultz M H. GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems. SIAM J Sci Stat Comput, 1986, 7: 856--869

[30] Karypis G, Kumar V. Parallel multilevel k-way partitioning scheme for irregular graphs. SIAM Review, 1999, 41: 278--300

[31] Catalyurek U, Boman E, Devine K, Bozdag D, Heaphy R, Riesen L. Hypergraph-based dynamic load balancing for adaptive scientific computations//Proc of IEEE Intl Parallel and Distributed Processing Symp (IPDPS). Long Beach, CA, March, 2007

[32] Beall M W, Shephard M S. A general topology-based mesh data structure. Int J Numer Meth Engng, 1997, 40(9): 1573--1596

[33] Shephard M S, Flaherty J E, Bottasso C L, de Cougny H L, Ozturan C, Simone M L. Parallel automated adaptive analysis. Parallel Computing, 1997, 23: 1327--1347

[34] Flaherty J E, Loy R M, Ozturan C, Shephard M S, Szymanski B K, Terresco J D, Ziantz L H. Parallel structures and dynamic load balancing for adaptive finite element computations. Appl Numer Math, 1998, 26(1/2): 241--263

[35] Li X, Shephard M S, Beall M W. 3-D anisotropic mesh adaptation by mesh modifications. Comp Meth Appl Mech Eng, 2005, 194(48/49): 4915--4950

[36] Sahni O, Mueller J, Jansen K E, Shephard M S, Taylor C A. Efficient anisotropic adaptive discretization of the cardiovascular system. Comp Meth Appl Mech Eng, 2006, 195(41/43): 5634--5655

[37] Li X, Shephard, M S, Beall M W. Accouting for curved domains in mesh adaptation. Int J Numer Meth Eng, 2003, 58: 246--276

[38] Shephard M S, Flaherty J E, Jansen K E, Li X, Luo X -J, Chevaugeon N, Remacle J -F, Beall M W, O'Bara R M. Adaptive mesh generation for curved domains. J Appl Num Math, 2005, 53(2/3): 251--271