权重积分方法计算界面张力消除表面“寄生”流

• TRANSPORT PHENOMENA & FLUID MECHANICS • 上一篇下一篇

权重积分方法计算界面张力消除表面“寄生”流

王剑峰^a,b; 杨超^a; 毛在砂^a

^a Institute of Process Engineering, Chinese Academy of Sciences, Beijing 100080, China
^b Chemical Engineering Department, Zhengzhou University, Zhengzhou 450052, China

收稿日期:1900-01-01 修回日期:1900-01-01 出版日期:2006-12-28 发布日期:2006-12-28
通讯作者: 王剑峰

A simple weighted integration method for calculating surface tension force to suppress parasitic flow in the level set approach

WANG Jianfeng^a,b; YANG Chao^a; MAO Zaisha^a

^a Institute of Process Engineering, Chinese Academy of Sciences, Beijing 100080, China
^b Chemical Engineering Department, Zhengzhou University, Zhengzhou 450052, China

Received:1900-01-01 Revised:1900-01-01 Online:2006-12-28 Published:2006-12-28
Contact: WANG Jianfeng

摘要/Abstract

摘要： Parasitic flows may occur in the numerical simulation of incompressible multiphase flow due to errors in the calculation of surface tension terms, specifically for the curvature and unit normal vector. An improved method for calculating the surface tension based on the level set approach is proposed, in which the contribution of not only the center node but also the rest area of a control volume to the calculation of surface tension is considered in a balanced manner. The weighted integration method (WIM) is more consistent with the concept of a banded in-terface in the level set method. It is applied to the temporal evolution of a two-dimensional neutrally buoyant liquid drop and a buoyancy driven deformable bubble in an immiscible fluid for the validation of WIM. The results show that the parasitic flows are evidently suppressed by the weighted integration method. The weight factors for WIM in 3-D cases are also suggested.

关键词: parasitic flow;weighted integration method;surface tension;level set approach;numerical simulation

Abstract: Parasitic flows may occur in the numerical simulation of incompressible multiphase flow due to errors in the calculation of surface tension terms, specifically for the curvature and unit normal vector. An improved method for calculating the surface tension based on the level set approach is proposed, in which the contribution of not only the center node but also the rest area of a control volume to the calculation of surface tension is considered in a balanced manner. The weighted integration method (WIM) is more consistent with the concept of a banded in-terface in the level set method. It is applied to the temporal evolution of a two-dimensional neutrally buoyant liquid drop and a buoyancy driven deformable bubble in an immiscible fluid for the validation of WIM. The results show that the parasitic flows are evidently suppressed by the weighted integration method. The weight factors for WIM in 3-D cases are also suggested.

Key words: parasitic flow, weighted integration method, surface tension, level set approach, numerical simulation

王剑峰,杨超, 毛在砂. 权重积分方法计算界面张力消除表面“寄生”流[J]. , DOI: .

WANGJianfeng,YANGChao, MAOZaisha. A simple weighted integration method for calculating surface tension force to suppress parasitic flow in the level set approach[J]. , DOI: .

参考文献

1    Brackbill, J.U., Kothe, D.B., Zemach, C., “A continuum method for modeling surface tension”, J. Comput. Phys., 100, 335—354(1992).
2    Coward, A.V., Renardy, Y., Renardy, M.J., Richards, R., “Temporal evolution of periodic disturbances in two-layer Couette flow”, J. Comput. Phys., 132, 346—361(1997).
3    Rider, W.J., Kothe, D.B., “Reconstructing volume track-ing”, J. Comput. Phys., 141, 112—152(1998).
4    Lafaurie, B., Nardone, C., Scardovelli, R., Zaleski, S., Zanetti, G., “Modelling merging and fragmentation in multiphase flows with SURFER”, J. Comput. Phys., 113, 134—147(1994).
5    Osher, S., Sethian, J., “Fronts propagating with curvature dependent speed: Algorithms based on Hmilton-Jacobi formulations”, J. Comput. Phys., 79, 12—49(1988).
6    Hirt, C.W., Nichols, B.D., “Volume of fluid (VOF) method for the dynamics of free boundaries”, J. Comput. Phys., 39, 201—225(1981).
7    Kothe, D.B., Rider, W.J., Mosso, S.J., Brock, J.S., “Vol-ume tracking of interfaces having surface tension in two and three dimensions”, In: The 34th Aerospace Sciences Meeting and Exhitit, AIAA, Reno, Nevada (1996).
8    Popinet, S., Zaleski, S., “A front-tracking algorithm for accurate representation of   surface tension”, Int. J. Numer. Meth. Fluids, 30, 775—793(1999)
9    Jamet, D., Torres, D., Brackbill, J.U., “On the theory and computation of surface tension: The elimination of para-sitic currents through energy conservation in the sec-ond-gradient method”, J. Comput. Phys., 182, 262—276(2002).
10    Renardy, Y., Renardy, M., “Prost: A parabolic reconstruc-tion of surface tension for the volume-of-fluid method”, J. Comput. Phys., 183, 400—421(2002).
11    Kang, M., Fedkiw, R.P., Liu, X.D., “A boundary condi-tion capturing method for multiphase incompressible flow”, J. Sci. Comput., 15, 323—360(2000).
12    Sussman, M., “A second order coupled level set and volume-of-fluid method for computing growth and col-lapse of vapor bubbles”, J. Comput. Phys., 187, 110—136(2003).
13    Cummins, S.J., Francois, M.M., Kothe, D.B., “Estimat-ing curvature from volume fractions”, Comput. Struct., 83, 425—434(2005).
14    Shin, S., Abdel-Khalik, S.I., Daru, V., Juric, D., “Accu-rate representation of surface tension using the level contour reconstruction method”, J. Comput. Phys., 203, 493—516(2005).
15    Shirani, E., Ashgriz, N., Mostaghimi, J., “Interface pres-sure calculation based on conservation of momentum for front capturing methods”, J. Comput. Phys., 203, 154—175(2005).
16    Yang, C., Mao, Z.S., “An improved level set approach to the simulation drop and bubble motion”, Chin. J. Chem. Eng., 10, 263—272(2002).
17     Bussmann, M., Mostaghimi, J., Chandra, S., “On a three-dimensional volume model of drop impact”, Phys. Fluids, 11, 1406—1417(1999).
18    van Doormaal, J.P., Raithby, G.D., “Enhancements of the SIMPLE method for prediction incompressible fluid flows”, Num. Heat Transfer, 7, 147—163(1984).
19    Patankar, S. V., Numerical Heat Transfer and Fluid Flow, Hemisphere, Washington (1980).
20    Scardovelli, R., Zaleski, S., “Direct numerical simulation of free surface and interfacial flow”, Ann. Rev. Fluid Mech., 31, 567—603(1999).
21    Bhaga, D., Weber, M.E., “Bubbles in viscous liquids: Shapes, wakes and velocities”, J. Fluid Mech., 105, 61—85(1981).