一种基于MOF界面重构的二维中心型MMALE方法

计算物理 ›› 2016, Vol. 33 ›› Issue (5): 523-538.

一种基于MOF界面重构的二维中心型MMALE方法

贾祖朋, 孙宇涛

北京应用物理与计算数学研究所, 北京 100094

收稿日期:2015-07-12 修回日期:2015-11-16 出版日期:2016-09-25 发布日期:2016-09-25
通讯作者: 孙宇涛,E-mail:syt02@mails.tsinghua.edu.cn
作者简介:贾祖朋(1967-),男,博士,研究员,从事计算数学及计算流体力学研究,E-mail:zpjia@iapcm.ac.cn
基金资助:
国家自然科学基金委员会与中国工程物理研究院联合基金资助项目

A 2D Cell-centered MMALE Method Based on MOF Interface Reconstruction

JIA Zupeng, SUN Yutao

Institute of Applied Physics and Computational Mathematics, Beijing 100094, China

Received:2015-07-12 Revised:2015-11-16 Online:2016-09-25 Published:2016-09-25

摘要/Abstract

摘要： 发展了一种基于MOF（Moment of Fluid）界面重构的二维中心型MMALE（Multi-Material Arbitrary Lagrangian-Eulerian）方法.其中，流体力学方程组采用中心型拉氏方法进行离散求解.混合网格的热力学封闭采用Tipton压力松弛模型.混合网格内的界面重构采用MOF方法，并对MOF方法作了简化和改进.重映步采用一种基于多边形剪裁算法的精确积分守恒重映方法.计算了若干数值例子，包括二维漩涡发展问题、Sedov问题、激波与氦气泡相互作用问题、水中强激波与空气泡相互作用问题、二维RT不稳定性问题等.数值算例表明，该方法具有二阶精度，能够计算界面两侧密度比和压力比很大的问题，并且其健壮性优于交错型MMALE方法，适合计算多介质复杂流体动力学问题.

关键词: MOF界面重构, 二维中心型MMALE方法

Abstract: A 2D cell-centered multi-material arbitrary Lagrangian-Eulerian(MMALE) method based on moment of fluid(MOF) interface reconstruction is developed. Hydrodynamic equations are discretized and solved by cell-centered Lagrangian scheme. Tipton's pressure relaxation model is adopted as closure model for mixed cells. MOF method is simplified and improved to reconstruct interface in mixed cells. A conservative integral remapping algorithm based on cell-intersection is used in remapping phase. Several numerical examples are given, such as 2D periodic vortex problem, Sedov problem, interaction of a shock wave with an helium bubble,a strong water shock impacting on a cylindrical air bubble in water,2D Rayleigh-Taylor instability, etc.It shows that the method is of second order accurate, and is capable of computing problems involving large density and/or pressure ratios across interface. Its robustness is better than that of staggered MMALE method and it is applicable to complex multi-material hydrodynamic problems.

Key words: MOF interface reconstruction, 2D cell-centered MMALE method

中图分类号:

O241.82

贾祖朋, 孙宇涛. 一种基于MOF界面重构的二维中心型MMALE方法[J]. 计算物理, 2016, 33(5): 523-538.

JIA Zupeng, SUN Yutao. A 2D Cell-centered MMALE Method Based on MOF Interface Reconstruction[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2016, 33(5): 523-538.

参考文献

[1] HIRT C W, AMSDEN A A, COOK H K. An arbitrary Lagrangian-Eulerian computing method for all flow speeds [J]. J Comput Phys,1974, 14:227-253.
[2] PEERY J S, CARROLL D E. Multi-material ALE methods in unstructured grids [J].Comput Methods Appl Mech Eng, 2000, 187: 591-619.
[3] HIRT C W, NICHOLS B D. Volume of fluid (VOF) method for the dynamics of free boundaries [J]. J Comput Phys, 1981, 39: 201-225.
[4] LIANG X H. Second order volume of fluid interface reconstruction method in three dimensions[J]. Chinese Journal of Computational Mechanics, 2013, 30(3): 353-360.
[5] 贾祖朋. 基于MOF界面重构的多物质ALE方法 [J]. 计算物理, 2010, 27(3): 353-360.
[6] JIA Z P, LIU J, ZHANG S D. An effective integration of methods for second-order three-dimensional multi-material ALE method on unstructured hexahedral meshes using MOF interface reconstruction [J]. J Comput Phys, 2013, 236: 513-562.
[7] 冯其金,刘军,王言金,等. 二维轴对称坐标下基于介质形心和体积的界面重构方法[J]. 中国科学:物理学力学天文学,2014,44(2):203-211.
[8] 曾清红,孙文俊. 多介质大变形流动的MOF-MMALE数值模拟研究 [J]. 应用数学和力学,2014,35(10):1163-1176.
[9] GODUNOV S K. Finite difference method for numerical computation of discontinuous solutions of the equations of fluid dynamics [J]. Matematicheskii Sbornik, 1959,47(3): 271-306.
[10] MAIRE P -H, ABGRALL R, BREIL J, et al. A cell-centered Lagrangian scheme for two-dimensional compressible flow problems [J]. SIAM Journal on Scientific Computing, 2007,29: 1781-1824.
[11] MAIRE P -H, BREIL J. A second-order cell-centered Lagrangian scheme for two-dimensional compressible flow problems [J]. Int J Numer Meth Fluids, 2007, 56: 1417-1423.
[12] MAIRE P -H. A high-order cell-centered Lagrangian scheme for two-dimensional compressible fluid flows on unstructured meshes [J]. J Comput Phys, 2009, 228: 2391-2425.
[13] DYADECHKO V, SHASHKOV M. Moment-of-fluid interface reconstruction[R/OL]. Technical Report LA-UR-05-7571,Los Alamos National Laboratory,Los Alamos, NM,Oct 2005. http://math.lanl.gov/~vdyadechko/doc/2005-mof.pdf.
[14] DYADECHKO V, SHASHKOV M. Reconstruction of multi-material interfaces from moment data [J]. J Comput Phys, 2008, 227: 5361-5384.
[15] GALERA S, MAIRE P H, BREIL J. A two-dimensional unstructured cell-centered multi-material ALE scheme using VOF interface reconstruction [J]. J Comput Phys, 2010, 229: 5755-5787.
[16] GALERA S, BREIL J, MAIRE P H. A 2D unstructured multi-material cell-centered arbitrary Lagrangian-Eulerian (CCALE) scheme using MOF interface reconstruction [J]. Comput Fluids, 2010, 229: 5755-5787.
[17] BENSON D J. A mixture theory for contact in multi-material Eulerian formulations [J]. Comput Methods Appl Mech Engrg, 1997,140: 59-86.
[18] VITALI E, BENSON D J. An extended finite element formulation for contact in multi-material arbitrary Lagrangian-Eulerian calculations [J].Int J Numer Meth Engng, 2006,67:1420-1444.
[19] SHASHKOV M. Closure models for multimaterial cells in arbitrary Lagrangian-Eulerian hydrocodes [J]. Int J Numer Meth Fluids, 2007, 56:1497-1504.
[20] KUCHARIK M, GARIMELLA R V, SCHOFIELD S P, et al. A comparative study of interface reconstruction methods for multi-material ALE simulations [J]. J Comput Phys, 2010, 229: 2432-2452.
[21] WINSLOW A M. Numerical solution of the quasilinear Poisson equation in a nonuniform triangle mesh [J]. J Comput Phys, 1996, 1(2): 149-172 (Reprinted in 1997,135(2): 128-138).
[22] KNUPP P. Achieving finite element mesh quality via optimization of the Jacobian matrix norm and associated quantities. Part I-A framework for surface mesh optimization [J]. Int J Numer Methods Eng, 2000, 48(3): 401-420.
[23] DUKOWICZ J, BAUMGARDNER J.Incremental remapping as a transport/advection algorithm [J].J Comput Phys,2000,160:318-335.
[24] 贾祖朋,蔚喜军. 基于level set的Eulerian-Lagrangian耦合方法及其应用 [J].力学学报, 2010,42(2):177-182.
[25] LIU Y K, YAN Y, SHI J Y. An efficient algorithm for line clipping against a polygon [J]. Chinese Journal of Computers, 1999,22(11):1209-1214.
[26] SEDOV L I. Similarity and dimensional methods in mechanics [M]. New York:Academic Press, 1959.
[27] BALL G J, HOWELL B P, LEIGHTON T G, et al. Shock-induced collapse of a cylindrical air cavity in water: a free-Lagrange simulation [J]. Shock Waves, 2000,10: 265-276.
[28] QUIRK J J, KAMI S. On the dynamics of a shock-bubble interaction [J]. J Fluid Mech, 1996, 318: 129-163.
[29] HAAS J F, STUTEVANT B. Interaction of weak-shock waves [J]. J Fluid Mech, 1987, 181: 41-76.
[30] LOUBèRE R, MAIRE P H, SJASHKOV M, et al. ReALE: A reconnection-based arbitrary-Lagrangian-Eulerian method [J]. J Comput Phys, 2010, 229: 4724-4761.
[31] WU J, XU S H, ZHAO N. A hybrid method for multiphase flows with large density ratio[J]. Chinese Journal of Computational Mechanics, 2013, 30(1): 1-10.
[32] GAO Y F, ZHANG X Y, LIU N. Particle finite element method for incompressible multi-fluid flows[J]. Chinese Journal of Computational Mechanics, 2013, 30(1): 35-43.