杨晓峰；何晓明 针对三相流体流动系统的守恒Allen-Cahn型相场模型，提出了一种全离散解耦有限元方法。 （英语） Zbl 1507.76126号 计算。方法应用。机械。工程师。 389，文章ID 114376，31 p.（2022）. 小结：在本文中，我们开发并分析了一种新的全离散解耦有限元方法，用于求解由三种不混溶流体组成的系统的流动耦合三元相场模型。基于L^2梯度流方法，采用守恒Allen-Cahn型动力学描述自由界面运动，其中使用多个非局部型拉格朗日乘子精确地保持每个相的体积。由于结合了几种有效的数值技术，包括两步后向差分格式、有限元离散化、处理非线性的显式SAV（标量辅助变量）方法、，以及Navier-Stokes方程的投影方法。在每个时间步长，非局部分裂技术只需要解几个解耦的常系数椭圆方程。详细讨论了实现问题。严格证明了该格式的可解性和无条件能量稳定性。通过大量的二维和三维数值模拟，从数值上证明了该方案的准确性、能量稳定性和适用性。 引用于20文件 MSC公司： 76M10个 有限元方法在流体力学问题中的应用 76立方英尺 三个或更多组件流 关键词：守恒Allen-Cahn模型；相位场；三相流；有限元法；全离散格式；无条件稳定性 PDF格式BibTeX公司 XML格式引用 \textit{X.Yang}和\textit{X.He}，计算。方法应用。机械。工程389，文章ID 114376，31 p.（2022；Zbl 1507.76126） 全文： 内政部 参考文献： [1] Boyer，F。；Lapuerta，C.，三组分Cahn Hilliard流动模型的研究，ESAIM：数学。模型。数字。分析。，40, 4, 653-687 (2006) ·Zbl 1173.35527号 [2] Boyer，F。；Minjeaud，S.，三分量Cahn Hilliard模型的数值格式，ESAIM：数学。模型。数字。分析。，45, 04, 697-738 (2011) ·Zbl 1267.76127号 [3] Kim，J。；Kang，K。；Lowengrub，J.，三元Cahn-Hilliard系统的保守多重网格方法，Commun。数学。科学。，2, 53-77 (2004) ·Zbl 1085.65093号 [4] Wu，S。；Xu，J.，多相Allen-Cahn和Cahn-Hilliard模型及其在成对表面张力影响下的离散，J.Compute。物理。，343, 10-32 (2017) ·Zbl 1380.65208号 [5] 杨，X。；赵，J。；王，Q。；Shen，J.，基于不变能量四次化方法的三分量Cahn-Hilliard相场模型的数值近似，数学。模型方法应用。科学。，27, 1993-2030 (2017) ·Zbl 1393.80003号 [6] 张杰。；Yang，X.，三相Cahn-Hilliard相场模型的解耦、非迭代和无条件能量稳定的大时间步进方法，J.Comput。物理。，404，第109115条pp.（2020）·Zbl 1453.76223号 [7] 张杰。；Yang，X.，基于L2颗粒流的三元相场模型的无条件能量稳定大时间步长方法，具有精确的非局部体积守恒，计算。方法应用。机械。工程，361，第112743条pp.（2020）·Zbl 1442.65186号 [8] Chen，W。；王，C。；王，S。；王，X。；Wise，S.M.，三元Cahn-Hilliard系统的能量稳定数值格式，J.Sci。计算。，84, 27 (2020) ·Zbl 1447.65098号 [9] Guermond，J.L。；Minev，P。；Shen，J.，《不可压缩流投影方法概述》，计算。方法应用。机械。工程，195，6011-6045（2006）·兹比尔1122.76072 [10] Guermond，J.L。；Quartapelle，L.，变密度不可压缩流动的投影有限元法，J.计算。物理。，165, 1, 167-188 (2000) ·Zbl 0994.76051号 [11] Guermond，J.L。；Salgado，A.，基于压力泊松方程的变密度不可压缩流的分裂方法，J.Compute。物理。，228, 8, 2834-2846 (2009) ·Zbl 1159.76028号 [12] 吉蒙德，J.-L。；Salgado，A.，变密度不可压缩流分数时间步进技术的误差分析，SIAM J.Numer。分析。，49, 3, 917-944 (2011) ·Zbl 1241.76318号 [13] 诺切托，R。；Pyo，J.-H.，《Gauge-Uzawa有限元方法第一部分：Navier-Stokes方程》，SIAM J.Numer。分析。，43, 1043-1068 (2005) ·Zbl 1094.76041号 [14] Pyo，J。；沈，J.，变密度不可压缩流动的Gauge-Uzawa方法，J.计算。物理。，221, 181-197 (2007) ·Zbl 1109.76037号 [15] 沈杰。；Yang，X.，不同密度和粘度的两相不可压缩流动的相场模型及其数值近似，SIAM J.Sci。计算。，32, 1159-1179 (2010) ·Zbl 1410.76464号 [16] 迪格尔，A。；王，C。；王，X。；Wise，S.，Cahn-Hilliard-Navier-Stokes系统二阶精确有限元方法的收敛性分析和误差估计，编号。数学。，137, 495-534 (2017) ·Zbl 1523.65081号 [17] Feng，X.，两相流体流动Navier-Stokes-Cahn-Hilliard扩散界面模型的全离散有限元近似，ESAIM:Math。模型。数字。分析。，44, 1049-1072 (2006) ·Zbl 1344.76052号 [18] Han，D。；Brylev，A。；杨，X。；Tan，Z.，两相不可压缩流相场模型的二阶全离散能量稳定格式的数值分析，科学杂志。计算。，70, 965-989 (2017) ·Zbl 1397.76070号 [19] Han，D。；Wang，X.，Cahn-Hilliard-Navier-Stokes方程的一个二阶时间、唯一可解、无条件稳定的数值格式，J.Compute。物理。，290, 139-156 (2015) ·Zbl 1349.76213号 [20] Han，D。；Wang，X.，Cahn-Hilliard-Darcy系统的二阶时间解耦无条件稳定数值格式，J.Sci。计算。，14, 1210-1233 (2018) ·兹比尔1407.65158 [21] Yang，X.，使用稳定不变能量四边形（S-IEQ）方法的流体动力学耦合三分量Cahn-Hilliard相场模型的高效和能量稳定格式，J.Compute。物理。，438，第110342条pp.（2021）·兹伯利07505949 [22] Yang，X.，一种新的流体动力学耦合三组分体积守恒Allen-Cahn相场模型的高效、二阶时间和能量稳定方案，数学。模型方法应用。科学。，31，4753-787（2021）·Zbl 1473.65247号 [23] 凯·D·。；样式，V。；Welford，R.，Cahn-Hilliard-Navier-Stokes系统的有限元近似，界面自由边界。，10, 15-43 (2008) ·Zbl 1144.35043号 [24] 刘，C。；沈杰。；Yang，X.，变密度两相不可压缩流相场模型的解耦能量稳定格式，J.Sci。计算。，62, 601-622 (2015) ·Zbl 1326.76064号 [25] Minjeaud，S.，三相Cahn-Hilliard/Navier-Stokes模型的无条件稳定解耦方案，Numer。偏微分方程方法，29584-618（2013）·Zbl 1364.76091号 [26] 沈杰。；Yang，X.，两相复杂流体相场模型的解耦能量稳定格式，SIAM J.Sci。计算。，36、B122-B145（2014）·Zbl 1288.76057号 [27] 沈杰。；Yang，X.，两相不可压缩流相场模型的解耦能量稳定格式，SIAM J.Numer。分析。，53, 1, 279-296 (2015) ·Zbl 1327.65178号 [28] Yang，X.，三相不可压缩流Cahn-Hilliard相场模型的一种新的高效全解耦二阶时间精确格式，计算。方法应用。机械。工程，376，第113589条pp.（2021）·兹比尔1506.76193 [29] Rubinstein，J。；Sternberg，P.，非局部反应扩散方程和成核，IMA J.Appl。数学。，48, 249-264 (1992) ·Zbl 0763.35051号 [30] Antonopoulou，D.C。；贝茨，P.W。；Blomker，D。；Karali，G.D.，随机质量守恒Allen-Cahn方程的液滴运动，SIAM J.Math。分析。，48, 670-708 (2016) ·Zbl 1439.35232号 [31] 柴，Z。；Sun，D。；Wang，H。；Shi，B.，具有质量守恒的局部和非局部Allen-Cahn方程的比较研究，《国际传热杂志》，122，631-642（2018） [32] 陈，X。；Hilhorst，D。；Logak，E.，质量守恒Allen-Cahn方程和体积守恒平均曲率流，界面自由边界。，12, 527-549 (2010) ·Zbl 1219.35018号 [33] Chiu，P.-H。；Lin，Y.-T.，求解不可压缩两相流的保守相场方法，J.Comput。物理。，230, 1, 185-204 (2011) ·Zbl 1427.76201号 [34] Jeong，D。；Kim，J.，不可压缩两相流体流动的保守Allen-Cahn-Navier-Stokes系统，计算与流体，156，239-246（2017）·Zbl 1390.76577号 [35] Kim，J。；Lee，S。；Choi，Y.，带时空相关拉格朗日乘数的保守Allen-Cahn方程，Internat。工程科学杂志。，84, 11-17 (2014) ·Zbl 1425.65089号 [36] 科塔尼，体育。；Hilhorst，D。；Lee，K.，具有非线性扩散的随机质量守恒反应扩散方程，离散Contin。动态。系统。序列号。A、 385615-5648（2018）·Zbl 1401.60122号 [37] 林，F。；他，X.-M。；Wen，X.，新Allen-Cahn型方形相场晶体模型的快速、无条件能量稳定大时间步进法，应用。数学。莱特。，92, 248-255 (2019) ·Zbl 1433.65250号 [38] Vaibhav，J。；Rajeev，K.J.，保守和正保持Allen-Cahn相场模型的自适应变分程序，J.Compute。物理。，366, 478-504 (2018) ·Zbl 1406.76049号 [39] Vaibhav，J。；Rajeev，K.J.，两相流相场建模的保正保守变分格式，J.Comput。物理。，360137-166（2018）·Zbl 1391.76651号 [40] Baskaran，A。；周，P。；胡，Z。；Wang，C.，修正相场晶体方程的能量稳定且有效的非线性有限差分多重网格格式，J.Compute。物理。，250, 270-292 (2013) ·Zbl 1349.65265号 [41] 秦，Y。；王，C。；Zhang，“二元流体-表面活性剂体系的保正收敛数值格式”，国际J·数值。分析。型号。，18, 399-425 (2021) ·Zbl 1499.65424号 [42] 怀斯，S.M。；王，C。；Lowengrub，J.S.，相场晶体方程的能量稳定收敛有限差分格式，SIAM J.Numer。分析。，47, 3, 2269-2288 (2009) ·Zbl 1201.35027号 [43] 王，C。；Wise，S.M.，修正相场晶体方程的能量稳定收敛有限差分格式，SIAM J.Numer。分析。，49, 945-969 (2011) ·Zbl 1230.82005年 [44] 胡，Z。；怀斯，S.M。；王，C。；Lowengrub，J.S.，相场晶体方程的稳定且有效的有限差分非线性多重网格格式，J.Compute。物理。，228, 5323-5339 (2009) ·Zbl 1171.82015年 [45] 杨，J。；毛，S。；他，X.-M。；杨，X。；He，Y.，两相磁流体力学流动的扩散界面模型和半隐式能量稳定有限元方法，计算。方法应用。机械。工程，356435-464（2019）·Zbl 1441.76143号 [46] 戈麦斯，H。；Nogueira，X.，相场晶体方程的无条件能量稳定方法，计算。方法应用。机械。工程，249-252，52-61（2012）·Zbl 1348.74280号 [47] 黄，Q。；杨，X。；He，X.-M.，蒙脱石-a液晶流动模型的数值近似：一阶、线性、解耦和能量稳定格式，离散Contin。动态。系统。序列号。B、 23、6、2177-2192（2018）·Zbl 1408.65094号 [48] 徐，C。；Tang，T.，外延生长模型大时间步进方法的稳定性分析，SIAM J.Numer。分析。，44, 1759-1779 (2006) ·Zbl 1127.65069号 [49] 杨，X。；Feng，J.J。；刘，C。；Shen，J.，使用能量变分相场方法对射流夹断和液滴形成进行数值模拟，J.Compute。物理。，218, 417-428 (2006) ·兹比尔1158.76319 [50] 徐，C。；陈，C。；杨，X。；He，X.-M.，流体动力学耦合二元表面活性剂相场模型的数值近似：二阶线性无条件能量稳定格式，Commun。数学。科学。，17, 3, 835-858 (2019) ·Zbl 1450.65090号 [51] 杨，X。；Han，D.，相场晶体方程的线性一阶和二阶无条件能量稳定格式，J.Compute。物理。，330, 1116-1134 (2017) ·Zbl 1380.65209号 [52] 杨，X。；赵，J。；He，X.-M.，利用不变能量求积方法求解双曲松弛粘性Cahn-Hilliard方程的线性、二阶和无条件能量稳定格式，J.Compute。申请。数学。，343，1，80-97（2018）·Zbl 1462.65117号 [53] 张杰。；Yang，X.，关于具有面心立方（FCC）有序结构的双模相场晶体模型的有效数值格式，应用。数字。数学。，146, 13-37 (2019) ·Zbl 1428.65017号 [54] 张杰。；Yang，X.，具有精确非局部质量守恒的新型L2颗粒流相场晶体模型的数值近似，计算。物理。通信，24351-67（2019）·Zbl 07674815号 [55] 张杰。；Yang，X.，铁磁性固体材料磁耦合相场晶体模型的高效准确数值格式，计算。方法应用。机械。工程，371，第113110条pp.（2020）·Zbl 1506.74122号 [56] 迪格尔，A。；X·冯。；Wise，S.，《Cahn-Hilliard-Darcy-Stokes系统的混合有限元分析》，SIAM J.Numer。分析。，53, 1, 127-152 (2015) ·Zbl 1330.76065号 [57] Feng，X.，两相流体流动Navier-Stokes-Cahn-Hilliard扩散界面模型的全离散有限元近似，SIAM J.Numer。分析。，44, 3, 1049-1072 (2006) ·Zbl 1344.76052号 [58] X·冯。；何毅。；Liu，C.，两相流体相场模型的有限元近似分析，数学。公司。，76, 258, 539-571 (2007) ·Zbl 1111.76028号 [59] X·冯。；Li，Y.，Allen-Cahn方程和平均曲率流的对称内部罚分不连续伽辽金方法分析，IMA J.Numer。分析。，1622-1651年4月35日（2015年）·兹比尔1328.65205 [60] X·冯。；李毅。；Xing，Y.，Cahn-Hilliard方程和Hele-Shaw流的混合内罚间断Galerkin方法分析，SIAM J.Numer。分析。，54, 2, 825-847 (2016) ·Zbl 1382.65314号 [61] X·冯。；Wise，S.，Hele-Shaw流的Darcy-Cahn-Hilliard扩散界面模型分析及其全离散有限元近似，SIAM J.Numer。分析。，50, 3, 1320-1343 (2012) ·Zbl 1426.76258号 [62] 高，Y。；他，X.-M。；梅，L。；Yang，X.，Cahn-Hilliard-Navier-Stokes-Darcy相场模型的解耦、线性和能量稳定有限元方法，SIAM J.Sci。计算。，40、1、B110-B137（2018）·Zbl 1426.76261号 [63] 高，Y。；他，X.-M。；Nie，Y.，Cahn-Hilliard-Darcy系统的二阶完全解耦、线性化和无条件稳定SAV方案，Numer。偏微分方程方法（2021） [64] 李，R。；高，Y。；陈，J。；张，L。；他，X.-M。；Chen，Z.，耦合Navier-Stokes-Cahn-Hilliard相场模型的间断有限体积元方法，高级计算。数学。，46, #25 (2020) ·Zbl 1436.65125号 [65] 诺切托，R。；萨尔加多，A。；Tomas，I.，两相铁流体扩散界面模型，计算。方法应用。机械。工程，309，497-531（2016）·Zbl 1439.76088号 [66] 诺切托，R。；萨尔加多，A。；Walker，S.，《带有移动接触线的电润湿扩散界面模型》，数学。模型方法应用。科学。，24, 1, 67-111 (2014) ·兹比尔1280.35114 [67] 雷霍尔茨，L。；怀斯，S。；Xiao，M.，《Cahn-Hilliard Navier-Stokes两相流扩散界面模型的罚投影格式及其与无发散耦合格式的联系》，国际J.Numer。分析。型号。，15, 649-676 (2018) ·Zbl 1395.35109号 [68] 张，G。；他，X.-M。；Yang，X.，两相铁流体动力学模型的解耦、线性和无条件能量稳定的全离散有限元数值格式，SIAM J.Sci。计算。，43、1、B167-B193（2021）·Zbl 1469.65152号 [69] 张，G。；他，X.-M。；Yang，X.，不可压缩MHD方程的二阶时间精度和无条件能量稳定性的完全解耦线性化有限元方法，J.Compute。物理。(2021) [70] 诺切托，R.H。；萨尔加多，A.J。；Tomas，I.，两相铁流体扩散界面模型，计算。方法应用。机械。工程，309，497-531（2016）·Zbl 1439.76088号 [71] 列夫，B.I。；Nazarenko，V.G。；Nych，A.B。；Tomchuk，P.M.，外部电场中向列相液滴的变形和不稳定性，JETP Lett。，71226-265（2000年） [72] 诺切托，R.H。；萨尔加多，A.J。；Walker，S.W.，带移动接触线的电润湿扩散界面模型，数学。模型方法应用。科学。，24, 67-111 (2014) ·兹比尔1280.35114 [73] Bai，F。；Han，D。；他，X.-M。；Yang，X.，在均匀磁场下使用相场和改进的水平集方法在Rosensweig模型中的铁液滴变形和聚并，Commun。非线性科学。数字。模拟。，第85条，第105213页（2020年）·Zbl 1460.76773号 [74] Bai，F。；他，X.-M。；周，R。；杨，X。；Wang，C.，微流体流动聚焦装置中液滴形成的三维相场研究及实验验证，国际J.Multiph。流量。，93, 130-141 (2017) [75] Han，D。；他，X.-M。；王，Q。；Wu，Y.，叠合自由流和多孔介质中Cahn-Hilliard-Navier-Stokes-Darcy系统解的存在性和弱强唯一性，非线性分析。，211, #112411 (2021) ·Zbl 1475.35279号 [76] Han，D。；王，X。；Wu，H.，岩溶几何中两相不可压缩流动Cahn-Hilliard-Stokes-Darcy系统整体弱解的存在性和唯一性，J.微分方程，257，10，3887-3933（2014）·Zbl 1302.35217号 [77] 沈杰。；Yang，X.，一类高度非线性梯度流系统的IEQ和SAV方法及其扩展，Contemp。数学。，754, 217-245 (2020) ·Zbl 1484.65196号 [78] Girault，V。；Raviart，P.A.，《Navier-Stokes方程的有限元方法：理论和算法》，395-414（1987），Springer-Verlag:Springer-Verlag Berlin，海德堡 [79] Shen，J.，关于Navier-Stokes方程投影方法的误差估计：二阶格式，数学。公司。，65, 215, 1039-1065 (1996) ·Zbl 0855.76049号 [80] E、 W。；Liu，J.G.，投影法。I.收敛和数值边界层，SIAM J.Numer。，32, 1017-1057 (1995) ·Zbl 0842.76052号 [81] 陈，C。；Yang，X.，各向异性Cahn Hilliard模型的快速、可证明无条件能量稳定和二阶精确算法，Comput。方法应用。机械。工程，351，35-59（2019）·Zbl 1441.74052号 [82] 沈杰。；徐，J。；Yang，J.，梯度流的一类新型高效稳健能量稳定方案，SIAM Rev.，61，474-506（2019）·Zbl 1422.65080号 [83] 沈杰。；薛，J。；Yang，J.，梯度流的标量辅助变量（SAV）方法，J.Compute。物理。，353, 407-416 (2018) ·Zbl 1380.65181号 [84] 杨，Z。；Dong，S.，基于隐式辅助能量变量的无条件能量稳定格式，用于不同密度的不可压缩两相流，仅涉及预计算系数矩阵，J.Compute。物理。，393, 229-257 (2018) ·Zbl 1452.76244号 [85] 陈，C。；Yang，X.，带熔体对流的枝晶凝固相场模型的高效数值格式，J.Compute。物理。，388，41-62（2019）·Zbl 1452.65195号 [86] Cheng，K。；王，C。；Wise，S.M.，方相场晶体方程的能量加权BDF2傅里叶伪谱数值格式，通信计算。物理。，26, 1335-1364 (2019) ·Zbl 1518.65115号 [87] 沈杰。；Yang，X.，Allen-Cahn和Cahn-Hilliard方程的数值近似，离散Contin。动态。系统。A、 2016年9月28日至1691年（2010年）·Zbl 1201.65184号 [88] Yang，X.，均聚物共混物相场模型的线性、一阶和二阶无条件能量稳定数值格式，J.Compute。物理。，327, 294-316 (2016) ·Zbl 1373.82106号 [89] 杨，X。；Yu，H.，使用不变能量求积方法的相场移动接触线模型的有效二阶无条件稳定格式，SIAM J.Sci。计算。，40，B889-B914（2018）·Zbl 1395.65095号 [90] Ingram，R.，Navier-Stokes方程的一种新的线性外推曲柄-尼科尔森时间步进格式，数学。公司。，82, 284, 1953-1973 (2013) ·Zbl 1457.65115号 [91] van Kan，J.，粘性不可压缩流的二阶精确压力修正方案，SIAM J.Sci。统计计算。，7, 3, 870-891 (1986) ·兹比尔0594.76023 [92] Boyer，F。；拉普尔塔，C。；Minjeaud，S。；Piar，B。；Quintard，M.，Cahn-Hilliard/Navier-Stokes三相流模拟模型，交通。多孔介质，82463-483（2010） [93] Kim，J。；Lowengrub，J.，三相流的相场建模和模拟，界面自由边界。，7, 435-466 (2005) ·Zbl 1100.35088号 [94] 罗林森，J.S。；Widom，B.，《毛细管分子理论》（1989），克拉伦登出版社：牛津克拉伦登出版公司 [95] 秦，T。；拉加布，S。；Yue，P.，粘性流中上升气泡与刚性表面相互作用的轴对称模拟，国际期刊Multiph。流量。，52, 60-70 (2013) 此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配，可能包含数据转换错误。在某些情况下，zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献，而不要求完整或完全匹配。