×

三元流体流动的相场流体动力学模型的物理信息结构保护数值格式。 (英文) Zbl 07814745号

摘要:相场模型广泛应用于多相流动力学研究。考虑到这些模型的复杂性和强非线性,设计精确、高效和稳定的数值算法来求解这些模型几十年来一直是一个活跃的研究领域。本文提出了一种新的数值格式,用于求解一个被广泛引用和使用的相场流体动力学模型,以模拟三元相流体流动。主要的新颖之处是引入了一个补充变量,将原始问题重新转化为约束优化问题。与求解该模型的许多现有数值技术相比,这种重新计算使我们提出的数值算法具有一些优点。首先,通过求解几个解耦的亥姆霍兹或泊松型系统以及一个恒定的预计算系数矩阵,所开发的方案允许对流体动力相场模型进行更直接的计算,显著降低了计算成本。其次,数值格式可以在离散水平上保持质量守恒和能量耗散。此外,基于二阶后向差分公式开发的方案尊重原始的能量耗散规律,这与许多现有方案不同,例如IEQ、SAV和Lagrange乘数法,这些方案尊重修改后的能量耗竭规律。此外,还提供了能量稳定性的严格证明和实际实施策略。我们进行了充分的二维和三维数值试验,以证明所提方案的准确性和有效性。

MSC公司:

52号B10 三维多面体
65D18天 计算机图形、图像分析和计算几何的数值方面
68单位05 计算机图形学;计算几何(数字和算法方面)
68单位07 计算机辅助设计的计算机科学方面
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] S.BOYAVAL、T.LELIEVRE和C.MANGOUBI,Oldroyd-B模型的自由能耗散方案,ESAIM数学。模型。数字。分析。43(3) (2009), 523-561. ·Zbl 1167.76018号
[2] F.BOYER和C.LAPUERTA,三组分Cahn-Hilliard流动模型的研究,ESAIM数学。模型。数字。分析。40 (2006), 653-687. ·Zbl 1173.35527号
[3] F.BOYER和S.MINJEAUD,三分量Cahn-Hilliard模型的数值格式,ESAIM数学。模型。数字。分析。45 (2011), 697-738. ·Zbl 1267.76127号
[4] L.CHEN,微观结构演化的相场建模,年。修订版材料。科学。32 (2002), 113-140.
[5] W.CHEN,C.WANG,S.WANG,X.WANG,AND S.WISE,三元Cahn-Hilliard系统的能量稳定数值格式,J.Sci。计算。84(2) (2020), 27. ·Zbl 1447.65098号
[6] 陈寅和沈俊杰,不可压缩Cahn-Hilliard Navier-Stokes相场模型的高效自适应能量稳定格式,J.Compute。物理学。308 (2016), 40-56. ·Zbl 1352.65229号
[7] Q.CHENG、C.LIU和J.SHEN,梯度流的新拉格朗日乘子方法,计算。应用方法。机械。工程367(2020),113070·Zbl 1442.65211号
[8] L.DONG,C.WANG,S.WISE,AND Z.ZHANG,具有奇异界面参数的三元Cahn-Hilliard系统的正守恒能量稳定格式,J.Compute。物理学。442 (2021), 110451. ·Zbl 07513797号
[9] D.EYRE,无条件梯度稳定时间推进Cahn-Hilliard方程,MRS在线程序库529(1998),39-46。
[10] H.GARCKE、B.NESTLER和B.STOTH,多相场概念:移动相边界和多结的数值模拟,SIAM J.Appl。数学。60 (2000), 295-315. ·Zbl 0942.35095号
[11] Y.GONG,Q.HONG,AND Q.WANG,热力学一致性偏微分方程的补充变量法,计算。应用方法。机械。工程381(2021),113746·Zbl 1506.74495号
[12] Z.GUO,Q.CHENG,P.LIN,C.LIU,AND J.LOWENGRUB,变密度两相流准不可压Navier-Stokes-Cahn-Hilliard系统的二阶近似,J.Compute。物理学。448 (2022), 110727. ·兹比尔1529.76062
[13] Z.GUO,P.LIN,J.LOWENGRUB,AND S.WISE,准不可压缩Navier-Stokes-Cahn-Hilliard系统的质量守恒和能量稳定的现场差分方法:原始变量和投影型格式,计算。应用方法。机械。工程326(2017),144-174·Zbl 1439.76121号
[14] Q.HONG,Y.GONG,J.ZHAO,AND Q.WANG,带非局部约束的Allen-Cahn模型的任意高阶结构保护算法,应用。数字。数学。170 (2021), 321-339. ·Zbl 1501.65080号
[15] 洪振强,李俊杰,王振强,用辅助变量法对带推导方程的偏微分方程进行结构保持逼近,应用。数学。莱特。110 (2020), 116576. ·Zbl 1452.65160号
[16] 华建华,林平,刘春秋,王庆庆,两相流计算中相场模型的能量守恒C 0有限元格式,J.Compute。物理学。230 (2011), 7115-7131. ·Zbl 1408.76550号
[17] 姜明明,张志军,赵建杰,用松弛法提高标量辅助变量(SAV)方法的准确性和一致性,J.Compute。物理学。456 (2022), 110954. ·Zbl 07518095号
[18] J.KIM,多组分流体流动的相场模型,Commun。计算。物理学。12 (2012), 613-661. ·兹比尔1373.76030
[19] J.KIM、K.KANG和J.LOWENGRUB,三元Cahn-Hilliard系统的保守多重网格方法,Commun。数学。科学。2 (2004), 53-77. ·Zbl 1085.65093号
[20] J.KIM和J.LOWENGRUB,三相流的相场建模和模拟,In-terfaces Free。边界7(2005),435-466·Zbl 1100.35088号
[21] KOU和SUN,彭罗宾逊状态方程非等温扩散界面两相流的热力学一致性模拟,J.Compute。物理学。371 (2018), 581-605. ·兹比尔1415.76644
[22] J.LEE和J.KIM,N分量Cahn-Hilliard系统的二阶精确非线性差分格式,Phys。A 387(2008),4787-4799。
[23] 李德华,乔志明,唐太宗,表征相场方程半隐式Fourier谱方法的稳定度,SIAM J.Numer。分析。54 (2016), 1653-1681. ·Zbl 1346.35162号
[24] C.LIU和J.SHEN,两种不可压缩流体混合物的相场模型及其通过傅里叶谱方法的近似,Phys。D 179(2003),211-228·Zbl 1092.76069号
[25] J.SHEN,C.WANG,X.WANG,AND S.M.WISE,《Ehrlich-Schwoebel型能量梯度流的二阶凸分裂方案:在薄膜外延中的应用》,SIAM J.Numer。分析。50 (2012), 105-125. ·Zbl 1247.65088号
[26] J.SHEN和X.YANG,Allen-Cahn和Cahn Hilliard方程的数值近似,离散Contin。动态。系统。28 (2010), 1669-1691. ·Zbl 1201.65184号
[27] L.SHEN、Z.XU、P.LIN、H.HUANG和S.XU,囊泡运动和变形相场模型的能量稳定的C 0有限元格式,SIAM J.Sci。计算。44(2022年),B122-B145·Zbl 07459363号
[28] 谭振中,吴建军,杨建军,任意形状三维表面上不可压缩多组分流体的高效实用相场方法,J.Compute。物理学。467 (2022), 111444. ·Zbl 07568547号
[29] 王强,广义Onsager原理及其应用,载《Froniter与当前软物质研究进展》,Springer(2020),101-132。
[30] S.WISE,Cahn-Hilliard-Hele-Shaw方程组的无条件稳定有限差分非线性多重网格模拟,J.Sci。计算。44 (2010), 38-68. ·Zbl 1203.76153号
[31] S.WISE、J.KIM和J.LOWENGRUB,用自适应非线性多重网格方法求解正则化强各向异性Cahn-Hilliard方程,J.Compute。物理学。226 (2007), 414-446. ·Zbl 1310.82044号
[32] C.XU和T.TANG,外延生长模型大时间步进方法的稳定性分析,SIAM J.Numer。分析。44 (2006), 1759-1779. ·Zbl 1127.65069号
[33] J.XU,Y.LI,AND S.WU,关于相场建模的部分隐式和完全隐式格式的稳定性和准确性,Comput。应用方法。机械。工程345(2019),826-853·Zbl 1440.80003号
[34] 杨振中,《Cahn-Hilliard方程组的无条件稳定二阶精确方法》,Commun。非线性科学。数字。模拟。87 (2020), 105276. ·Zbl 1455.76138号
[35] 杨振中,利用有效的修正标量辅助变量方法对三元Cahn-Hilliard流体进行数值研究,Commun。非线性科学。数字。模拟。102 (2021), 105923. ·Zbl 1497.65157号
[36] X.YANG,使用稳定不变能量二次化(IEQ)的流体动力学耦合三分量Cahn-Hilliard相场模型的有效和能量稳定方案,J.Comput。物理学。438 (2021), 110342. ·Zbl 07505949号
[37] X.YANG,一种新的流体动力学耦合三分量体积守恒Allen-Cahn相场模型的高效、二阶时间和能量稳定格式,数学。Mod-els方法应用。科学。31 (2021), 753-787. ·Zbl 1473.65247号
[38] X.YANG,J.ZHAO,Q.WANG,AND J.SHEN,基于不变能量积分法的三分量Cahn-Hilliard相场模型的数值近似,数学。模型方法应用。科学。27 (2017), 1993-2030. ·Zbl 1393.80003号
[39] Z.YANG和S.DONG,基于隐式辅助能量变量的无条件能量稳定格式,用于不同密度的不可压缩两相流,仅涉及预计算系数矩阵,J.Compute。物理学。393 (2019), 229-257. ·Zbl 1452.76244号
[40] Z.YANG和S.DONG,基于广义辅助变量和保证正性的耗散系统离散能量稳定方案路线图,J.Compute。物理学。404 (2020), 109121. ·Zbl 1453.65276号
[41] M.YUAN,W.CHEN,C.WANG,S.WISE,AND Z.ZHANG,高分子微球复合水凝胶三组分Cahn-Hilliard型模型的能量稳定有限元格式,J.Sci。计算。87 (2021), 78. ·Zbl 1473.65219号
[42] J.ZHANG和X.YANG,三相Cahn-Hilliard相场模型的解耦、非迭代和无条件能量稳定大时间步进法,J.Compute。物理学。404 (2020), 109115. ·Zbl 1453.76223号
[43] J.ZHANG和X.YANG,具有精确非局部体积守恒的基于L2梯度流的三元相场模型的无条件能量稳定大时间步进方法,计算。应用方法。机械。工程361(2020),112743·Zbl 1442.65186号
[44] 赵军,用松弛技术重温能量求积法,应用。数学。莱特。120 (2021), 107331. ·Zbl 1498.35538号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。