×

Cahn-Hilliard-Navier-Stokes方程的二阶解耦能量稳定格式。 (英语) Zbl 07515435号

概述:Cahn-Hilliard-Navier-Stokes(CHNS)方程代表了多相流体流动动力学的流体动力相场模型的基本构建块。由于Navier-Stokes方程和Cahn-Hilliard方程之间的耦合,CHNS系统的数值求解非常困难。传统上,使用耦合项的数值外推。然而,这种蛮力外推通常会破坏CHNS系统的固有热力学结构。本文提出了一种新的策略,将CHNS系统重新构造为约束梯度流公式,其中清楚地揭示了可逆和不可逆结构。这指导我们提出具有几个有利性质的算子分裂方案。首先,所提出的方案导致在每个时间步进中求解几个规模较小的解耦系统。这大大降低了计算成本。其次,所提出的方案仍然保证了CHNS系统在离散水平上的热力学定律。此外,与最近使用辅助变量的IEQ或SAV方法不同,我们得到的能量定律是在原始变量中制定的。这是一个显著的改进,因为带有辅助变量的修正能量定律有时会偏离原始能量定律。我们提出的框架为水动力相场模型的解耦和能量稳定数值算法的设计奠定了基础。此外,给定不同的分裂步骤,可以获得各种数值算法,这使得该框架相当通用。实现了所提出的数值算法。数值验证了它们的二阶时间和空间精度。通过算例和基准问题验证了所提方案的有效性。

理学硕士:

6500万 偏微分方程、初值和含时初边值问题的数值方法
7.6亿 流体力学基本方法
35季度xx 数学物理偏微分方程及其他应用领域
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] 艾伦,S.M。;Cahn,J.W.,《反相边界运动的微观理论及其在反相畴粗化中的应用》,《金属学报》。材料。,27, 1085-1095 (1979)
[2] Bao,Y。;Kim,J.,使用相位场模型的多相位图像分割,计算。数学。申请。,62, 737-745 (2011) ·邮编:1228.94009
[3] 贝尔托齐,A。;Esedoglu,S。;Gillette,A.,使用Cahn-Hilliard方程对二值图像进行修复,IEEE Trans。图像处理。,16, 1, 285-291 (2007) ·Zbl 1279.94008号
[4] Boyer,F。;拉佩尔塔,C。;明让,S。;Piar,B。;Quintard,M.,Cahn-Hilliard/Navier-Stokes三相流模拟模型,交通。多孔介质,82463-483(2010)
[5] 卡恩,J.W。;Hilliard,J.E.,非均匀系统的自由能。I.界面自由能,J.Chem。物理。,28, 258-267 (1958) ·Zbl 1431.35066号
[6] Chen,L。;Zhao,J.,Cahn-Hilliard-Navier-Stokes方程的新型二阶线性格式,J.Compute。物理。,423,第109782条pp.(2020)·Zbl 07508402号
[7] Chen,L.Q.,微观结构演化的相场建模,年。修订版材料。研究,32,113-140(2002)
[8] 陈,Y。;Shen,J.,《不可压缩Cahn-Hilliard Navier-Stokes相场模型的高效自适应能量稳定格式》,J.Compute。物理。,308, 40-56 (2016) ·Zbl 1352.65229号
[9] 丁·H。;斯佩尔特,P.D.M。;Shu,C.,大密度比不可压缩两相流扩散界面模型,J.Compute。物理。,226, 2, 2078-2095 (2007) ·Zbl 1388.76403号
[10] 杜琪。;刘,C。;Pyham,R。;Wang,X.,细胞膜自发曲率效应的相场模拟,Commun。纯应用程序。数学。,4, 3, 537-548 (2005) ·Zbl 1078.92004号
[11] Eyre,D.,《无条件梯度稳定时间推进Cahn-Hilliard方程》,(微观结构演化的计算和数学模型。微观结构演化计算和数学模式,加利福尼亚州旧金山,1998年,第529卷(1998年)),39-46
[12] 冯·W。;王,C。;怀斯,S。;Zhang,Z.,带斜率选择的外延薄膜方程的二阶能量稳定后向微分公式法,Numer。方法部分差异。Equ.、。,34, 6, 1975-2007 (2018) ·Zbl 1407.76095号
[13] Flory,P.J.,《高分子溶液的热力学》,J.Chem。物理。,1860年9月8日(1941年)
[14] 戈麦斯,H。;Hughes,T.J.R.,《相场模型的可证明无条件稳定的二阶时间精度混合变分方法》,J.Compute。物理。,230, 13, 5310-5327 (2011) ·Zbl 1419.76439号
[15] 龚,Y。;赵,J。;Wang,Q.,二元粘性流体流体动力相场模型交错网格上的二阶全离散能量稳定方法,SIAM J.Sci。计算。,40、2、B528-B553(2018)·Zbl 1393.65012号
[16] Griffith,B.E.,《利用投影法作为预条件,求解不可压缩Navier-Stokes方程的一种精确而有效的方法》,J.Compute。物理。,228, 20, 7565-7595 (2009) ·Zbl 1391.76474号
[17] Grmela,M。;Ottinger,H.C.,《复杂流体的动力学和热力学i一般形式主义的发展》,《物理学》。E版,56、6(1997年)
[18] 郭杰。;王,C。;怀斯,S。;Yue,X.,三维Cahn-Hilliard方程二阶凸分裂有限差分格式的收敛性,Commun。数学。科学。,14, 2, 489-515 (2015) ·Zbl 1338.65221号
[19] 郭,Z。;林,P。;Lowengrub,J.,《离散能量定律下变密度流动的准不可压缩Cahn-Hilliard-Navier-Stokes方程的数值方法》,J.Compute。物理。,276, 486-507 (2014) ·Zbl 1349.76057号
[20] 郭,Z。;林,P。;Lowengrub,J。;Wise,S.,拟不可压缩Navier-Stokes-Cahn-Hilliard系统的质量守恒和能量稳定有限差分方法:原始变量和投影型格式,计算。方法应用。机械。工程,326144-174(2017)·Zbl 1439.76121号
[21] Han,D。;Wang,X.,Cahn-Hilliard-Navier-Stokes方程的二阶时间唯一可解无条件稳定数值格式,J.Compute。物理。,第290页,第139-156页(2015年)·Zbl 1349.76213号
[22] Huggins,M.L.,《长链化合物的溶液》,J.Chem。物理。,9, 5, 440 (1941)
[23] 凯·D·。;Welford,R.,《二维Cahn-Hilliard-Navier-Stokes流体的高效数值解》,SIAM J.Sci。计算。,29, 6, 2241-2257 (2007) ·Zbl 1154.76033号
[24] Kim,J。;Lowengrub,J.,三相流的相场建模和模拟,界面自由边界。,7, 4, 435-466 (2005) ·Zbl 1100.35088号
[25] 寇,J。;Sun,S。;Wang,X.,多组分两相可压缩流动的线性解耦能量稳定数值方法,SIAM J.Numer。分析。,56, 6, 3219-3248 (2018) ·Zbl 1404.65226号
[26] Lee,H。;Shin,J。;Lee,J.,修正相场晶体方程的一阶和二阶能量稳定方法,计算。方法应用。机械。工程,321,1-17(2017)·Zbl 1439.74240号
[27] 李,X。;Shen,J.,关于两相不可压缩流Cahn-Hilliard-Navier-Stokes模型的完全解耦MSAV方案(2020)
[28] 刘,C。;弗兰克·F。;蒂勒,C。;Alpak,F.O。;Berg,S。;查普曼,W。;Riviere,B.,求解多孔介质中粘性对比度Cahn-Hilliard-Navier-Stokes系统的高效数值算法,J.Compute。物理。,400,第108948条pp.(2020)·Zbl 1453.76078号
[29] Lober,J。;齐伯特,F。;Aranson,I.S.,可变形细胞的碰撞导致集体迁移,科学。代表,5,第9172条pp.(2015)
[30] Lowengrub,J。;Truskinovsky,L.,《准不可压缩Cahn-Hilliard流体和拓扑转变》,Proc。R.Soc.A,数学。物理学。工程科学。,454, 1978, 2617-2654 (1998) ·Zbl 0927.76007号
[31] Minjeaud,Sebastian,三相Cahn-Hilliard Navier-Stokes模型的无条件稳定解耦方案,Numer。方法部分差异。Equ.、。,29, 2, 584-618 (2013) ·Zbl 1364.76091号
[32] Ottinger,H.C。;Grmela,M.,《复杂流体的动力学和热力学ii一般形式主义图解》,Phys。E版,56、6(1997年)
[33] Palmieri,B。;Bresler,Y。;维茨,D。;Grant,M.,单层细胞迁移的多尺度模型:细胞之间的弹性错配增强运动,科学。代表,5,第11745条pp.(2015)
[34] 邵,D。;Pappel,W。;Levine,H.,细胞形态动力学计算模型,物理学。修订稿。,105(2010年9月)
[35] Shen,J.,用相场方法对两相不可压缩流动进行建模和数值近似,(材料模拟的多尺度建模和分析。材料模拟的多级建模和分析,Lect.Notes Ser.Inst.Math.Sci.Natl.Univ.Singap.,第22卷(2012),《世界科学》。出版物:世界科学。出版物。新泽西州哈肯萨克),147-195年
[36] 沈杰。;王,C。;王,X。;Wise,S.,《具有Ehrlich-Schwoebel型能量的梯度流的二阶凸分裂方案:在薄膜外延中的应用》,SIAM J.Numer。分析。,50, 1, 105-125 (2012) ·Zbl 1247.65088号
[37] 沈杰。;徐,J。;Yang,J.,梯度流的标量辅助变量(SAV)方法,J.Compute。物理。,353, 407-416 (2018) ·Zbl 1380.65181号
[38] 沈杰。;Yang,X.,两相复杂流体相场模型的解耦能量稳定格式,SIAM J.Sci。计算。,36, 1, 122-145 (2014) ·Zbl 1288.76057号
[39] Wang,L。;Yu,H.,关于Cahn-Hilliard相场方程的有效二阶稳定半隐式格式,J.Sci。计算。,77, 2, 1185-1209 (2018) ·Zbl 1407.65163号
[40] 王,S。;Sekerka,R。;惠勒,A。;Murray,B.T。;Coriell,S.R。;Braun,R.J。;McFadden,G.B.,凝固热力学一致相场模型,《物理D》,69,189-200(1993)·Zbl 0791.35159号
[41] 王,X。;Du,Q.,《通过扩散界面方法对多组分脂质膜和开放膜进行建模和模拟》,J.Math。《生物学》,56,347-371(2008)·Zbl 1143.92001年
[42] Wise,S.,Cahn-Hilliard-Hele-Shaw方程组的无条件稳定有限差分非线性多重网格模拟,J.Sci。计算。,44, 38-68 (2010) ·Zbl 1203.76153号
[43] 怀斯,S。;Kim,J。;Lowengrub,J.,用自适应非线性多重网格法求解正则化强各向异性Cahn-Hilliard方程,J.Compute。物理。,226414-446(2007年)·Zbl 1310.82044号
[44] Yang,X.,Navier-Stokes方程与体积守恒多相场囊泡系统耦合的数值近似:完全解耦、线性、无条件能量稳定和二阶时间精确的数值格式,Comput。方法应用。机械。工程,375,第113600条pp.(2021)·Zbl 1506.76105号
[45] 杨,X。;赵,J。;Wang,Q.,基于不变能量求积方法的分子束外延生长模型的数值近似,J.Compute。物理。,333, 102-127 (2017) ·兹比尔1375.82121
[46] 赵,J。;李,H。;王,Q。;Yang,X.,三相不可压缩粘性流体流动相场模型的解耦能量稳定格式,J.Sci。计算。,70, 1367-1389 (2017) ·兹比尔1397.76098
[47] 赵,J。;沈毅。;Happasalo,M。;王振杰。;Wang,Q.,异质多物种生物膜中抗菌药物持久性的三维数值研究,J.Theor。生物学,39283-98(2016)·Zbl 1347.92048号
[48] 赵,J。;杨,X。;龚,Y。;X.赵。;杨,X。;李,J。;Wang,Q.,《非平衡模型数值近似的一般策略——第一部分:热力学系统》,国际期刊Numer。分析。型号。,15, 6, 884-918 (2018) ·Zbl 1412.65099号
[49] 赵,J。;杨,X。;沈杰。;Wang,Q.,向列相液晶和粘性流体混合物的流体动力相场模型的解耦能量稳定方案,J.Compute。物理。,305, 539-556 (2016) ·Zbl 1349.76019号
[50] 周,C。;岳,P。;Feng,J.J。;Ollivier-Gooch,C.F。;Hu,H.H.,牛顿流体和粘弹性流体界面动力学的三维相场模拟,J.Compute。物理。,229498-511(2010年)·Zbl 1355.76037号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。