扩散界面形式下一种节约时间步的质量分数混合模型

doi:10.19596/j.cnki.1001-246x.8453

计算物理 ›› 2022, Vol. 39 ›› Issue (5): 510-520.DOI: 10.19596/j.cnki.1001-246x.8453

扩散界面形式下一种节约时间步的质量分数混合模型

吴宗铎¹(), 严谨¹, 宗智², 庞建华³^,*()

1. 广东海洋大学海洋工程学院, 广东湛江 524088
2. 大连理工大学船舶工程学院, 辽宁大连 116024
3. 广东海洋大学深圳研究院智能海洋工程中心, 广东深圳 518055

收稿日期:2021-09-22 出版日期:2022-09-25 发布日期:2023-01-07
通讯作者: 庞建华
作者简介:
吴宗铎(1984-), 男, 博士, 副教授, 主要研究激波间断数值方法、水下爆炸冲击波模拟及SPH算法数值模拟等, E-mail: wuzongduo0@aliyun.com
基金资助:
国家自然科学基金(11702066); 广东省自然科学基金(2022A1515011562); 广东省促进经济高质量发展专项(GDNRC[2021]56); 广东省促进经济高质量发展专项(GDNRC[2021]161)

A Diffused Interface Type Mass Fraction Model with Time Steps Saving Mixture Rules

Zongduo WU¹(), Jin YAN¹, Zhi ZONG², Jianhua PANG³^,*()

1. College of Ocean Engineering, Guangdong Ocean University, Zhanjiang, Guangdong 524088, China
2. Department of Naval Architecture, Dalian University of Technology, Dalian, Liaoning 116024, China
3. Ocean Intelligence Technology Center, Shenzhen Institute of Guangdong Ocean university, Shenzhen, Guangdong 518055, China

Received:2021-09-22 Online:2022-09-25 Published:2023-01-07
Contact: Jianhua PANG

摘要/Abstract

摘要：

研究扩散形式的界面处理中混合物模型的处理方式对计算的影响。为了让界面的处理既简便又具有较强的适应性, 选用Mie-Grüneisen混合模型, 对原有基于体积分数的混合物模型进行改进。改进后的混合物模型, 以质量分数为色函数, 并利用流场中介质的密度完成对声速计算的优化。对声速的优化后, 可以在收敛机制下缩短时间步的数量, 提高计算效率。同时, 限制器的引入有助于使压力曲线变得稳定。数值算例证明, 改进后的质量分数模型, 能在多种计算格式下实现时间步数的降低。时间步数减少的幅度取决于物理条件。

关键词: Mie-Grüneisen混合模型, 质量分数, 时间步, 限制器

Abstract:

We focused on effects of treatment of mixture model as the interface is in a diffused style. In order to seek a convenient way to treat the interface and keep strong adaptability, Mie-Grüneisen mixture model is used. The original Mie-Grüneisen mixture model in terms of volume fraction is modified. In the model, the mass fraction is employed as color function, and the optimization of sound speed is achieved with the help of the density of fluid components. The optimized sound speed reduces time steps under the convergence mechanism and promotes efficiency of the calculation. Meanwhile, the introduction of limiter make the pressure curves stable. Numerical tests show that the modified mass fraction model reduces the number of time steps under different numerical schemes. Effects of this reduction depend on physical conditions of case.

Key words: Mie-Grüneisen mixture model, mass fraction, time steps, limiter

吴宗铎, 严谨, 宗智, 庞建华. 扩散界面形式下一种节约时间步的质量分数混合模型[J]. 计算物理, 2022, 39(5): 510-520.

Zongduo WU, Jin YAN, Zhi ZONG, Jianhua PANG. A Diffused Interface Type Mass Fraction Model with Time Steps Saving Mixture Rules[J]. Chinese Journal of Computational Physics, 2022, 39(5): 510-520.

图/表 9

图1 参考压力pref曲线

Fig.1 Curves of reference pressure pref

图2 爆轰气体与铜作用下，73 μs时刻的ρ、p计算结果

Fig.2 ρ and p results of the interaction between detonation gaseous product and copper at 73 μs

图3 算例1中，73 μs时刻声速计算结果

Fig.3 Results of the sound velocity at 73 μs in Case 1

图4 铜冲击惰性炸药时，85 μs的ρ、u计算结果

Fig.4 ρ and u results at 85 μs as an inert explosive suffers from the impact of copper

图5 算例2中，85 μs时刻声速计算结果

Fig.5 Results of the sound velocity at 85 μs in Case 2

图6 冲击波位置随时间的变化

Fig.6 Time evolutions of the positions of shock wave

图7 200 μs和500 μs时刻，NM爆炸实验的数值解

Fig.7 200 μs和500 μs solutions for NM explosion problem

图8 140 μs时刻，二维NM爆炸的压力分布

Fig.8 Pressure distribution at 140 μs for 2D NM explosion problem

图9 二维NM爆炸算例的声速云图对比

Fig.9 Contrast of sound velocity contours for 2D NM explosion problem

参考文献 26

1	TORO E F. Riemann solvers and numerical methods for fluid dynamics[M]. Springer-Verlag Berlin Heidelberg, 2009.
2	高斯, 刘铁钢. 多介质黎曼问题精确解解算器软件包MultiRP开发[J]. 数值计算与计算机应用, 2016, 37 (4): 245- 256.
3	LIU T G , KHOO B C , YEO K S . Ghost fluid method for strong shock impacting on material interface[J]. Journal of Computational Physics, 2003, 190 (2): 651- 681. DOI
4	XIAO F , HONMA Y , KONO T . A simple algebraic interface capturing scheme using hyperbolic tangent function[J]. International Journal for Numerical Methods in Fluids, 2005, 48 (9): 1023- 1040. DOI
5	XU L , FENG C L , LIU T G . Analysis of impact load on elastic plate in underwater explosions with modified ghost fluid method[J]. Chinese Journal of Computational Physics, 2017, 34 (1): 1- 9. DOI
6	赵海洋, 明平剑, 张文平, 等. 基于THINC/QQ格式的两相界面流动数值模拟[J]. 哈尔滨工程大学学报, 2020, 41 (12): 1804- 1810.
7	ALLAIRE G , CLERC S , KOKH S . A five-equation model for the simulation of interfaces between compressible fluids[J]. Journal of Computational Physics, 2002, 181 (2): 577- 616. DOI
8	刘铁钢, 许亮. 模拟模拟多介质界面问题的虚拟流体方法综述[J]. 气体物理, 2019, 4 (2): 1- 16.
9	ABGRALL R . How to prevent pressure oscillations in multicomponent flow calculations: A quasi-conservative approach[J]. Journal of Computational Physics, 1999, 25 (1): 150- 160.
10	SAUREL R , ABGRALL R . A simple method for compressible multifluid flows[J]. SIAM Journal on Scientific Computing, 1999, 21 (3): 1115- 1145. DOI
11	SAUREL R , ABGRALL R . A multiphase Godunov method for compressible multifluid and multiphase flows[J]. Journal of Computational Physics, 1999, 150 (2): 425- 467. DOI
12	梁姗, 刘伟, 袁礼. 七方程可压缩多相流模型的HLLC格式及应用[J]. 力学学报, 2012, 44 (5): 884- 895.
13	XUE C , LI X D , SUN W J , et al. A multi-material five-equation-reduced model and artificial compression method for interface capture[J]. Chinese Journal of Computational Physics, 2021, 38 (03): 257- 268.
14	陈芳, 李平, 刘坤, 等. 基于PPM的界面压缩方法研究[J]. 高压物理学报, 2019, 33 (5): 052302.
15	SHYUE K M . A fluid-mixture type algorithm for compressible multicomponent flow with Mie-Grüneisen equation of state[J]. Journal of Computational Physics, 2001, 171 (2): 678- 707.
16	SHYUE K . A wave-propagation based volume tracking method for compressible multicomponent flow in two space dimensions[J]. Journal of Computational Physics, 2006, 215 (1): 219- 244.
17	张宝坪, 张庆民, 黄风雷. 爆轰物理学[M]. 第1版北京: 兵器工业出版社, 2001.
18	WANG Z , ZOU G Y , GONG J , et al. A pressure-based algorithm for numerical simulation of one-dimensional two-phase flow[J]. Chinese Journal of Computational Physics, 2019, 36 (4): 413- 420.
19	吴宗铎, 赵勇, 严谨, 等. 球坐标系下多介质混合物模型的数值模拟[J]. 爆炸与冲击, 2019, 39 (5): 054204.
20	任玉新, 陈海昕. 计算流体力学基础[M]. 北京: 清华大学出版社, 2006.
21	LIU Y , MAO D K . An ADRS approximated Riemann solver in compatible Lagrangian method[J]. Chinese Journal of Computational Physics, 2020, 37 (2): 140- 152.
22	WU Z D , YAN J , ZONG Z , et al. A multi-component Mie-Grüneisen mixture model based on HLLC algorithm[J]. Chinese Journal of Computational Physics, 2020, 37 (1): 55- 62.
23	SHREYAS B , SARMA L R . Investigation of numerical viscosities and dissipation rates of second-order TVD-MUSCL schemes for implicit large-eddy simulation[J]. Journal of Computational Physics, 2015, 281 (2): 1003- 1031.
24	LIU J J , ZENG X Y , NI G X . Euler numerical methods for reactive flow with general equation of states in two dimensions[J]. Chinese Journal of Computational Physics, 2018, 35 (1): 29- 38.
25	STEINBERG D J. Spherical explosion and the equation of state of water, UCID-20974[R]. Livermore: Lawrence Livermore National Laboratory, 1978.
26	LEE E, FINGER M, COLLINS W. JWL equation of state coefficients for high explosives, UCID-16189[R]. Livermore: Lawrence Livermore National Laboratory, 1973.

[1]	沈亚玲, 封建湖, 郑素佩, 李雅蓉. 一种基于新型斜率限制器的理想磁流体方程的高分辨率熵相容格式[J]. 计算物理, 2022, 39(3): 297-308.
[2]	艾邦成, 张亮, 陈智. 低耗散非结构网格多维限制器[J]. 计算物理, 2018, 35(5): 545-553.
[3]	任炯, 封建湖, 刘友琼, 梁楠. 求解双曲守恒律方程的高分辨率熵相容格式[J]. 计算物理, 2014, 31(5): 539-551.
[4]	唐维军, 蒋浪, 程军波. 多组份流动质量分数保极值原理算法[J]. 计算物理, 2014, 31(3): 292-306.
[5]	许仁义, 钟德钰, 吴保生. 随机选取法和多波近似在一维浅水方程大时间步长格式中的应用[J]. 计算物理, 2013, 30(5): 649-658.
[6]	唐维军, 袁光伟, 沈智军, 程军波. Lagrange方法基于保正性的时间步长[J]. 计算物理, 2012, 29(5): 633-640.
[7]	吴宗铎, 宗智. Mie-Grüneisen状态方程下多介质守恒型欧拉方程组的数值计算[J]. 计算物理, 2011, 28(6): 803-809.
[8]	吴迪, 蔚喜军, 徐云. 局部时间步长间断有限元方法求解三维欧拉方程[J]. 计算物理, 2011, 28(1): 1-9.
[9]	罗力, 封建湖, 唐小娟, 向量. 求解双曲守恒律方程的高分辨率熵稳定格式[J]. 计算物理, 2010, 27(5): 671-678.
[10]	汪洪波, 孙明波, 梁剑寒, 王振国. 含双时间步法的化学非平衡流解耦算法[J]. 计算物理, 2010, 27(5): 685-691.
[11]	王东红, 赵宁, 刘剑明. 界面追踪方法中的激波限制器研究[J]. 计算物理, 2009, 26(4): 510-516.
[12]	王汉奎, 张雄, 刘岩. 并行化光滑分子动力学方法及其与分子动力学的耦合[J]. 计算物理, 2008, 25(6): 718-724.
[13]	赵慧勇, 乐嘉陵. 双时间步方法的应用分析[J]. 计算物理, 2008, 25(3): 253-258.
[14]	赵海洋, 刘伟, 杨小亮, 任兵. 采用Power限制器的PNND和PWNND格式[J]. 计算物理, 2007, 24(6): 660-666.
[15]	符尚武, 戴自换, 邬吉明. 二维三温流体力学计算中时间步长的自动控制[J]. 计算物理, 2005, 22(6): 45-49.

扩散界面形式下一种节约时间步的质量分数混合模型

A Diffused Interface Type Mass Fraction Model with Time Steps Saving Mixture Rules

RichHTML

PDF

可视化

摘要/Abstract

引用本文

使用本文

图/表 9

参考文献 26

相关文章 15

编辑推荐

Metrics

作者中心

审稿中心

期刊浏览

期刊介绍