王晓荣;张晓迪 无电感磁流体力学相场模型的解耦、线性、无条件能量稳定和电荷守恒有限元方法。 （英语） Zbl 07764085号 数学。计算。模拟。 215, 607-627 (2024). 摘要：在本文中，我们考虑了两相不可压缩无感应磁流体力学（MHD）问题扩散界面模型的数值近似。该模型由Cahn-Hilliard方程、Navier-Stokes方程和Poisson方程组成。我们提出了一种线性解耦有限元方法来求解这个高度非线性的多物理系统。对于时间变量，离散化是一阶Euler半隐式格式、几个一阶稳定项和耦合项的隐式显式处理的组合。对于空间变量，我们采用有限元离散化。特别地，我们通过inf-sup稳定的面-体积混合有限元对来近似电流密度和电势。利用这些技术，该方案只需要在每个时间步长求解一系列解耦线性方程组。我们证明了该方案是可证明的质量守恒、电荷守恒和无条件能量稳定的。通过数值实验验证了该方案的性能、准确性和有效性。 MSC公司： 76倍 流体力学 65-XX岁 数值分析 关键词：无电感磁流体动力学方程;Cahn-Hilliard方程;混合有限元法;解耦方案;能量稳定的;电荷守恒的 软件：自由女性++ PDF格式BibTeX公司 XML格式引用 \textit{X.Wang}和\textit{X.Zhang}，数学。计算。模拟。215607-627（2024；Zbl 07764085） 全文： 内政部 arXiv公司 参考文献： [1] Abdou，M。；Sze，D。；Wong，C。；Sawan，M。；Ying，A。；莫利，N.B。；Malang，S.，美国ITER覆盖测试计划和战略，Fusion Sci。技术。，47，3475-487（2005年） [2] Abdou，M。；APEX团队，D。；Ying，A。；莫利，N。；Gulec，K。；Smolentsev，S。；Kotschenreuther，M。；马朗，S。；Zinkle，S。；Rognlien，T。；Fogarty，P。；B.纳尔逊。；Nygren，R。；麦卡锡，K。；优素福，M。；古尼姆，N。；Sze，D。；Wong，C。；Sawan，M。；Khater，H。；Woolley，R。；马塔斯，R。；莫尔，R。；沙拉法特，S。；布鲁克斯，J。；Hassanein，A。；佩蒂，D。；蒂拉克，M。；Ulrickson，M。；Uchimoto，T.，《关于聚变室技术创新概念的探索》，聚变工程设计。，54, 2, 181-247 (2001) [3] 安德森，D.M。；麦克法登，G.B。；Wheeler，A.A.，《流体力学中的扩散界面方法》，（《流体力学年鉴》，第30卷（1998年），《年度评论》，加利福尼亚州帕洛阿尔托），139-165·Zbl 1398.76051号 [4] Arnold，D.N。；布雷齐，F。；Fortin，M.，Stokes方程的稳定有限元，Calcolo，21，4，337-344（1985）（1984）·兹比尔0593.76039 [5] 巴迪亚，S。；马汀，A.F。；Planas，R.，热耦合不可压缩无电感磁流体力学问题的块递归LU预条件，J.Compute。物理。，274, 562-591 (2014) ·Zbl 1351.76046号 [6] 洛杉矶卡法雷利。；Muler，N.E.，Cahn-Hilliard方程解的界，Arch。定额。机械。分析。，133, 2, 129-144 (1995) ·Zbl 0851.35010号 [7] Cahn，J.W.，关于旋节分解，金属学报。，9, 9, 795-801 (1961) [8] Cahn，S.M.A.J.W.，反相边界运动微观理论及其在反相畴粗化中的应用，金属学报。，27, 1085-1095 (1979) [9] 卡恩，J.W。；Hilliard，J.E.，非均匀系统的自由能。I.界面自由能，J.Chem。物理。，28, 2, 258 (1958) ·Zbl 1431.35066号 [10] Ceniceros，H.D。；Nós，R.L。；Roma，A.M.，《相场流体模型的三维完全自适应模拟》，J.Compute。物理。，229, 17, 6135-6155 (2010) ·Zbl 1425.76179号 [11] 陈，R。；Zhang，H.，Cahn-Hilliard MHD方程新模型的二阶能量稳定格式，高级计算。数学。，46, 6, 79 (2020) ·Zbl 1465.65082号 [12] 陈，C。；Zhang，T.，两相磁流体扩散界面模型的一阶半隐式稳定有限元方法的无条件稳定性和最佳误差估计，应用。数学。计算。，429（2022），论文编号127238，22·Zbl 1510.76084号 [13] Ciarlet，P.G.，（《椭圆问题的有限元方法》，《应用数学经典》，第40卷（2002年），工业和应用数学学会（SIAM），费城），xxviii+530·Zbl 0999.65129号 [14] Davidson，P.A.，（《磁流体动力学导论》，《磁流体力学导论》（2001），剑桥大学出版社：剑桥大学出版社，剑桥），xviii+431·Zbl 0974.76002号 [15] 丁·L。；曹毅。；李，H。；刘振华，铝电解槽磁流体力学数值模拟（中文），冶金。马特。金属。工程，42，4，8-13（2014） [16] 丁·H。；Spelt，医学博士。；Shu，C.，大密度比不可压缩两相流扩散界面模型，J.Compute。物理。，226, 2, 2078-2095 (2007) ·Zbl 1388.76403号 [17] Drazin，P.G。；Reid，W.H.，（《流体动力稳定性.流体动力稳定性》，剑桥数学图书馆（2004），剑桥大学出版社：剑桥大学出版社），xx+605·Zbl 1055.76001号 [18] 埃本贝克，M。；Garcke，H。；Nürnberg，R.，《肿瘤生长的Cahn-hilliard-brinkman系统》，离散Contin。动态。系统。序列号。S、 14、11、3989-4033（2021年）·Zbl 1480.35411号 [19] Eyre，D.J.，梯度系统的无条件稳定一步格式，15（1997），未出版 [20] Feng，J.J。；刘，C。；沈杰。；Yue，P.，复杂流体界面动力学的相场方法的能量变分公式：优势和挑战，（软物质建模。软物质建模，IMA Vol.Math.Appl.，Vol.141（2005），Springer，纽约），1-26·Zbl 1181.76019号 [21] Gerbeau，J.-F.，《不可压缩磁流体动力学方程的稳定有限元方法》，数值。数学。，87, 1, 83-111 (2000) ·Zbl 0988.76050号 [22] 杰布·J·F。；勒布里斯，C。；Lelièvre，T.，液态金属磁流体动力学的数学方法，（数值数学和科学计算（2006），牛津大学出版社：牛津大学出版社），xiv+310·Zbl 1107.76001号 [23] Hecht，F.，自由有限元++的新发展，J.Numer。数学。，2012年3月20日至4日，251-265·Zbl 1266.68090号 [24] Hua，J。；林，P。；刘，C。；Wang，Q.，两相流计算中相场模型的能量守恒C0有限元格式，J.Compute。物理。，230, 19, 7115-7131 (2011) ·Zbl 1408.76550号 [25] 黄，H。；Ying，A。；Abdou，M.，基于磁场感应方程的三维MHD自由表面流体流动模拟，Fusion Eng.Des。，63-64, 361-368 (2002) [26] 约翰斯顿，H。；Liu，J.-G.，基于压力项显式处理的精确、稳定和高效的Navier-Stokes解算器，J.Compute。物理。，199，1221-259（2004年）·Zbl 1127.76343号 [27] Ki，H.，磁场下两相不可压缩流动的水平集方法，计算。物理学。Comm.，181，6，999-1007（2010）·Zbl 1454.76063号 [28] Kim，J.，多组分流体流动的相场模型，Commun。计算。物理。，12, 3, 613-661 (2012) ·Zbl 1373.76030号 [29] Lee，H.G。；Kim，J.，多组分流体的二维Kelvin-Helmholtz不稳定性，《欧洲力学杂志》。B Fluids，49，A部分，77-88（2015）·Zbl 1408.76143号 [30] 李，L。；镍，M。；Zheng，W.，无电感MHD方程的电荷守恒有限元方法。第一部分：汇聚，SIAM J.Sci。计算。，41、4、B796-B815（2019）·兹比尔1447.65084 [31] 刘，C。；Shen，J.，两种不可压缩流体混合物的相场模型及其用傅里叶谱方法的近似，Physica D，179，3-4，211-228（2003）·Zbl 1092.76069号 [32] 毛，S。；Wang，X.，两相不可压缩无感应磁流体动力学流体三维扩散界面模型的全离散有限元近似（2021），提交出版 [33] Monk，P.，（麦克斯韦方程的有限元方法。麦克斯韦方程式的有限元法，数值数学和科学计算（2003），牛津大学出版社：牛津大学出版社，纽约），xiv+450·Zbl 1024.78009号 [34] Moreau，R.，（磁流体动力学。磁流体动力学、流体力学及其应用，第3卷（1990），Kluwer学术出版集团，多德雷赫特）·Zbl 0714.76003号 [35] Munger博士。；Vincent，A.，模拟铝电解槽磁流体动力学不稳定性的水平集方法，J.Compute。物理。，217, 2, 295-311 (2006) ·Zbl 1103.76046号 [36] Ni，M.-J。；Li，J.-F.，低磁雷诺数下不可压缩MHD流动的一致保守方案。第三部分：在交错网格上，J.Compute。物理。，231, 2, 281-298 (2012) ·Zbl 1316.76079号 [37] Ni，M.-J。；穆尼帕利，R。；莫利，N.B。；黄，P。；Abdou，M.A.，低磁雷诺数下不可压缩MHD流动的电流密度守恒方案。I.在矩形并置网格系统上，J.Comput。物理。，227, 1, 174-204 (2007) ·Zbl 1280.76045号 [38] 潘，Y。；张杰。；王，Z.-H。；Ni，M.-J.，基于自适应网格细化技术的MHD求解器开发，Fusion Eng.Des。，87, 5-6, 630-633 (2012) [39] 拉维亚特，P.A。；Thomas，J.M.，二阶椭圆问题的混合有限元方法，（Galligani，I.；Magenes，E.，有限元方法的数学方面（1977年），Springer Berlin Heidelberg:Springer Barlin Heitelberg Berlin，Heidelbrg），292-315·Zbl 0362.65089号 [40] Rayleigh，L.，《表面力理论》。二、。可压缩流体，菲洛斯。Mag.Ser.杂志。5, 33, 209-220 (1892) [41] 萨穆利亚克，R。；杜，J。；格利姆，J。；Xu，Z.，低磁雷诺数下自由表面流动MHD的数值算法，J.Compute。物理。，226, 2, 1532-1549 (2007) ·Zbl 1173.76391号 [42] 沈，J.，用相场方法对两相不可压缩流动进行建模和数值近似，（材料模拟的多尺度建模和分析。材料模拟的多级建模和分析，国立大学学报，第22卷（2012年），《世界科学》。出版物。，新泽西州哈肯萨克），147-195年 [43] 沈，J.，用相场方法对两相不可压缩流动进行建模和数值近似，（材料模拟的多尺度建模和分析。材料模拟的多级建模和分析，国立大学学报，第22卷（2012年），《世界科学》。出版物。，新泽西州哈肯萨克），147-195年 [44] 沈杰。；徐，J。；Yang，J.，梯度流的一类新型高效稳健能量稳定方案，SIAM Rev.，61，3，474-506（2019）·Zbl 1422.65080号 [45] 沈杰。；Yang，X.，Allen-Cahn和Cahn-Hilliard方程的数值近似，离散Contin。动态。系统。，28, 4, 1669-1691 (2010) ·Zbl 1201.65184号 [46] 沈杰。；Yang，X.，不同密度和粘度的两相不可压缩流动的相场模型及其数值近似，SIAM J.Sci。计算。，32, 3, 1159-1179 (2010) ·Zbl 1410.76464号 [47] 沈杰。；Yang，X.，两相不可压缩流相场模型的解耦能量稳定格式，SIAM J.Numer。分析。，53, 1, 279-296 (2015) ·Zbl 1327.65178号 [48] 苏，H。；Zhang，G.-D.，两相磁流体动力学2D/3D扩散界面模型的高效和能量稳定方案，J.Sci。计算。，90，1（2022），论文编号63，31·兹伯利07454946 [49] 苏，H。；Zhang，G.-D.，两相磁流体动力学扩散界面模型的二阶时间精度和解耦结构的能量稳定格式，Commun。非线性科学。数字。模拟。，119（2023），论文编号107126·Zbl 1508.65134号 [50] 苏，H。；Zhang，G.-D.，两相磁流体动力学扩散界面模型的二阶时间精度和解耦结构的能量稳定格式，Commun。非线性科学。数字。模拟。，第107126条pp.（2023）·Zbl 1508.65134号 [51] Takatani，K.，《带自由表面的不可压缩MHD流的数学建模》，ISIJ Int.，47，4，545-551（2007） [52] Temam，R.，Navier-STokes方程。《理论与数值分析》（1977年），North-Holland Publishing Co.，阿姆斯特丹-纽约-Oxford，x+500，数学及其应用研究，第2卷·Zbl 0383.35057号 [53] van der Waals，J.D.，《密度连续变化假设下的毛细现象热力学理论》，J.Stat.Phys。，20, 197-200 (1893) [54] 谢伟新。；蔡，L。；Feng，J.-H.，用加权鬼流体法追踪理想MHD方程中的熵波，应用。数学。型号。，31, 11, 2503-2514 (2007) ·Zbl 1117.76075号 [55] 杨，J。；毛，S。；何，X。；杨，X。；He，Y.，两相磁流体力学流动的扩散界面模型和半隐式能量稳定有限元方法，计算。方法应用。机械。工程，356435-464（2019）·兹比尔1441.76143 [56] 杨，J。；毛，S。；何，X。；杨，X。；He，Y.，两相磁流体力学流动的扩散界面模型和半隐式能量稳定有限元方法，计算。方法应用。机械。工程，356435-464（2019）·Zbl 1441.76143号 [57] 杨，X。；赵，J。；王，Q。；Shen，J.，基于不变能量求积方法的三分量Cahn-Hilliard相场模型的数值近似，数学。模型方法应用。科学。，27, 11, 1993-2030 (2017) ·Zbl 1393.80003号 [58] Zhang，X.，两相无感应磁流体动力学流体扩散界面模型的尖锐界面极限（2021） [59] 张，X。；Ding，Q.，无电感磁流体动力学方程的解耦、无条件能量稳定和电荷守恒有限元方法，计算。数学。申请。，127, 80-96 (2022) ·Zbl 1524.76533号 [60] 张杰。；Ni，M.，带相变的不可压缩多相磁流体动力学的直接数值模拟，J.Compute。物理。，375, 717-746 (2018) ·Zbl 1416.76237号 [61] 赵，J。；陈，R。；Su，H.，不可压缩磁流体动力学-Cahn-Hilliard耦合模型的能量稳定有限元方法，Adv.Appl。数学。机械。，13, 4, 761-790 (2021) ·Zbl 1488.65489号 [62] 周，J。；Chen，L。；黄，Y。；Wang，W.，Cahn-Hilliard方程的有效双网格格式，Commun。计算。物理。，17, 1, 127-145 (2015) ·兹比尔1373.76111 此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配，可能包含数据转换错误。在某些情况下，zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献，而不声称其完整性或完全匹配。