A Compatible Cell-Centered ALE Method Based on General Riemann Solver
-
摘要: 深入分析了精确Riemann解法器、MFCAV(Multi Fluid Channel on Averaged Volume)和Dukowicz近似Riemann解法器不能直接应用于相容拉氏方法的原因,通过引入更一般形式的角点算子, 成功将以上3种Riemann解法器应用于新相容拉氏方法中;进一步结合高效的网格重分、重映技术,建立了一种基于任意Riemann解法器的相容中心型显式两步任意拉格朗日-欧拉(ALE)方法。将新的ALE方法应用于数值算例中,结果表明,新ALE方法不仅具有新相容拉氏方法的优点,而且具有处理大变形流动问题的能力。
-
关键词:
- 中心型拉氏方法 /
- 交错型拉氏方法 /
- 近似Riemann解法器
Abstract: We first deeply analyzed the reason why the exact Riemann sovler, MFCAV (Multi fluid channel on averaged volume) Riemann solver and Dukowicz Riemann solver cannot be applied directly to the compatible Lagrangian method and then, by introducing a new general nodal solver, we successfully applied the 3 Riemann solvers to the Lagrangian method.Combing the new Lagrangian method with efficiency remesh and remapping technology, we developed a new cell-centered compatible ALE (Arbitrary Lagrangian and Eulerian) method based on arbitrary Riemann solvers.Several numerical examples show that the new ALE method not only has the advantages of the new compatible Lagrangian method, but also has the ability to deal with the large deformation problems. -
表 1 Taylor Green Vortex算例在0.75时刻的误差分析
Table 1. Error analysis of the 3 Riemann solvers for Taylor Green Vortex problem at t=0.75
Grid
numberExact MFCAV Dukowicz L1 error
Convergence
rateL1 error
Convergence
rateL1 error Convergence
rate40 0.357 465 380 0.214 033 305 0.357 462 330 80 0.211 823 462 0.754 7 0.116 933 955 0.872 1 0.211 822 688 0.754 9 160 0.115 855 157 0.870 6 0.061 305 690 0.931 6 0.115 854 999 0.870 5 320 0.060 805 808 0.929 5 0.031 570 952 0.957 4 0.060 805 780 0.930 0 -
[1] Loubere R, Shashkov M J. A subcell remapping method on staggered polygonal grids for arbitrary Lagrangian-Eulerian method[J]. J Comput Phys, 2005, 209(1): 105-138. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=yysxhlx-e201910009 [2] Duckowicz J K, Meltz B J A. Vorticity errors in multidimensional Lagrangian codes[J]. J Comput Phys, 1992, 99(1): 115-134. doi: 10.1007/3-540-54960-9_62 [3] Addession F L, Carroll D E, Dukowicz J K, et, al. CAVEAT: A computer code for fluid dynamics problem with large distortion and internal slip, LA-10613-MS[R]. New Mexico: Los Alamos National Lab, 1992. [4] 田保林, 申卫东, 刘妍, 等. ALE框架下几种不同Godunov型格式的数值比较[J].计算物理, 2007, 24(5): 537-542. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=jswl200705006Tian B L, Shen W D, Liu Y, et al. Numerical comparing of several Godunov type schemes in Lagrangian framework[J]. Chinese Journal of Computational Physics, 2007, 24(5): 537-542. (in Chinese) http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=jswl200705006 [5] Maire P H, Abgrall R, Breil J, et al. A cell-centered Lagrangian scheme for two-dimensional compressible flow problems[J]. SIAM J Sci Comput, 2007, 29(4): 1781-1824. doi: 10.1137/050633019 [6] 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(7): 2391-2425. doi: 10.1016/j.jcp.2008.12.007 [7] Liu Y, Tian B L, Shen W D, et al. Application of MFCAV Riemann solver to Maire's cell-centered Lagrangian method[J]. J Comput Math, 2015, 37(3): 286-298. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=jssx-e201502002 [8] 李德元, 徐国荣, 水鸿寿, 等.二维非定常流体力学数值方法[M].北京: 科学出版社, 1987.Li D Y, Xu G R, Shui H S, et al. Two-Dimensional Hydrodynamic Numerical Algorithm[M]. Beijing: Science Press, 1987. (in Chinese) [9] Winslow A M. Adaptive mesh zoning by the equipotential method, UCID-9062[R]. California: Lawrence Livermore National Laboratory, 1981. [10] Benson D J. Momentum advection on a staggered mesh[J]. J Comput Phys, 1992, 100(1): 143-162. http://www.sciencedirect.com/science/article/pii/002199919290299E [11] Noh W F. Errors for calculations of strong shocks using artifical viscosity and artifical heat flux[J]. J Comput Phys, 1987, 72(1): 78-120. http://www.sciencedirect.com/science/article/pii/002199918790074X [12] Tian B L, Shen W D, Jiang S, et al. An arbitrary Lagrangian-Eulerian method based on the apaptive Riemann solver for general equations of state[J]. Int J Numer Meth Fluids, 2009, 59(11): 1217-1240. doi: 10.1002/fld.1871 [13] Woodward P, Colella P. The numerical simulation of two-dimensional fluid flow with strong shocks[J]. J Comput Phys, 1984, 54(1): 115-173. http://www.ams.org/mathscinet-getitem?mr=748569 [14] 刘妍, 田保林, 申卫东, 等. MFCAV近似Riemann解法器在相容拉氏方法中的熵条件分析[J].计算数学学报, 2015, 33(2): 1-18. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=jssx201503006Liu Y, Tian B L, Shen W D, et al. The analysis of entropy condition for MFCAV Riemann solver in a compatible Lagrangian method[J]. Journal of Computational Mathematics, 2015, 33(2): 1-18. (in Chinese) http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=jssx201503006 [15] Brio M, Zakharian A R, Webb G M. Two-dimensional Riemann solver for Euler equations of gas dynamics[J]. J Comput Phys, 2001, 167(1): 177-195. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=6eg32jCmOcxOiUravPVT3BNlMlk7aK5ixuQJy9lQSmE=