×

具有Flory-Huggins-de Gennes自由能的Cahn-Hilliard相场模型的非迭代、无条件能量稳定和大时间步长方法。 (英语) Zbl 1442.65223号

小结:在本文中,我们考虑用Flory-Huggins-de Gennes自由能数值近似求解均聚物共混物的Cahn-Hilliard相场模型。我们开发了一个有效的、二阶精度的、无条件能量稳定的方案,该方案将SAV方法与稳定技术相结合,其中(H^1)从总自由能中分离范数,并添加两个额外的线性稳定项,以提高稳定性,并在使用大时间步长时保持所需的精度。该方案非常容易实现,并且不需要迭代,只需在每个时间步长求解两个解耦的四阶常系数双调和方程。我们进一步严格证明了该方案的无条件能量稳定性。通过对一些二维和三维基准数值算例与非稳定SAV和MSAV等其他流行格式的比较,我们从数值上证明了所开发格式的稳定性和准确性。

MSC公司:

65平方米2 含偏微分方程初值和初边值问题离散方程的数值解
6500万06 含偏微分方程初值和初边值问题的有限差分方法
65M70型 偏微分方程初值和初边值问题的谱、配置及相关方法
31A30型 二维双调和、多调和函数和方程、泊松方程
82年第35季度 与统计力学相关的PDE
82D60型 聚合物统计力学
PDF格式BibTeX公司 XML格式引用
全文: DOI程序

参考文献:

[1] 粘结剂,K.,聚合物混合物中的集体扩散、成核和旋节分解,化学杂志。物理。,79, 6387 (1983) ·数字对象标识代码:10.1063/1.445747
[2] 陈,C。;Yang,X.,带熔体对流的枝晶凝固相场模型的高效数值格式,J Compute。物理。,388, 41-62 (2019) ·Zbl 1452.65195号 ·doi:10.1016/j.jcp.2019.03.017
[3] 陈,C。;Yang,X.,各向异性Cahn-Hilliard模型的快速、可证明的无条件能量稳定和二阶精度算法,计算。方法。申请。机械。工程,351,35-59(2019)·Zbl 1441.74052号 ·doi:10.1016/j.cma.2019.03.030
[4] 陈,F。;Shen,J.,各向异性Cahn-Hilliard系统空间谱离散的高效能量稳定格式,Commun Compute。物理。,05, 1189-1208 (2013) ·Zbl 1388.65117号 ·doi:10.4208/cicp.101111.110512a
[5] Chen,L。;赵,J。;Yang,X.,带斜率选择的分子束外延模型的正则化线性方案,应用。数字数学。,128, 139-156 (2018) ·Zbl 1486.65191号 ·doi:10.1016/j.apnum.2018.02.004
[6] Cheng,Q。;Shen,J.,多标量辅助变量(MSAV)方法及其在相场囊泡膜模型中的应用,SIAM J.Sci。计算。,40,A3982-A4006(2018)·Zbl 1403.65045号 ·doi:10.1137/18M1166961
[7] Cheng,Q。;杨,X。;Shen,J.,流体动力学耦合相场二嵌段共聚物模型的高效准确数值格式,J.Comp。物理。,341, 44-60 (2017) ·兹比尔1380.65203 ·doi:10.1016/j.jcp.2017.04.010
[8] 库克,HE,旋节分解运动中的布朗运动,《金属学报》,18297(1970)·doi:10.1016/0001-6160(70)90144-6
[9] 科佩蒂,MIM;Elliott,CM,具有对数自由能的Cahn-Hilliard方程的数值分析,数值数学。,63, 4, 39-65 (1992) ·Zbl 0762.65074号 ·doi:10.1007/BF01385847
[10] de Gennes,PG,《聚合物物理中的尺度概念》(1979),伊萨卡:康奈尔大学出版社
[11] de Gennes,PG,聚合物共混物中涨落和旋节分解动力学,J.Chem。物理。,7, 4756 (1980) ·Zbl 1110.82310号 ·数字对象标识代码:10.1063/1.439809
[12] 德彪西,A。;Dettori,L.,关于具有对数自由能的Cahn-Hilliard方程。非线性分析:理论、方法与应用,24,10,1491-1514(1995)·Zbl 0831.35088号
[13] 埃利奥特,CM;Garcke,H.,关于简并迁移率的Cahn-Hilliard方程,SIAMJ数学。分析。,27, 404-423 (1996) ·Zbl 0856.35071号 ·doi:10.1137/S0036141094267662
[14] Eyre,D.J.:无条件梯度稳定时间推进Cahn-Hilliard方程。In:显微结构演化的计算和数学模型(加州旧金山,1998),《材料》第529卷。Res.Soc.交响乐。程序。,第39-46页。沃伦代尔夫人(1998)
[15] X·冯。;Prol,A.,Allen-Cahn方程的数值分析和平均曲率流的近似,Numer。数学。,94, 33-65 (2003) ·Zbl 1029.65093号 ·doi:10.1007/s00211-002-0413-1
[16] Fialkowski,M。;Holyst,R.,《聚合物共混物中相分离动力学:形态、旋节、噪声和成核,宏观摩尔》。理论模拟。,17, 263 (2008) ·doi:10.1002/mats.20080020
[17] 福斯特,S。;阿拉斯加州坎德布尔;赵,J。;Bates,FS,聚异戊二烯-聚苯乙烯二嵌段共聚物在有序-二元转变附近的复杂相行为,大分子,276922-6935(1994)·doi:10.1021/ma00101a033
[18] 高,Y。;何,X。;梅,L。;Yang,X.,Cahn-Hilliard-Navier-Stokes-Darcy相场模型的解耦、线性和能量稳定有限元方法,SIAM J.Sci。计算。,40,B110-B137(2018)·Zbl 1426.76261号 ·doi:10.1137/16M1100885
[19] 戈麦斯,H。;Calo,VM;Bazilevs,Y。;Hughes,TJR,Cahn-Hilliard相场模型的等几何分析,计算。方法应用。机械。工程,197,4333-4352(2008)·Zbl 1194.74524号 ·doi:10.1016/j.cma.2008.05.003
[20] Gomez,H.,Van der,Z.,Kristoffer,G.:计算相场建模。收录:计算力学百科全书,第二版。John Wiley&Sons有限公司,ISBN 978-1-119-00379-3(2017)
[21] 戈麦斯,H。;Hughes,TJR,《相场模型的可证明无条件稳定、二阶时间精度混合变分方法》,J.Compute。物理。,230, 5310-5327 (2011) ·Zbl 1419.76439号 ·doi:10.1016/j.jp.2011.03.033
[22] Han,D。;Brylev,A。;杨,X。;Tan,Z.,两相不可压缩流相场模型的二阶全离散能量稳定格式的数值分析,科学计算杂志。,70, 965-989 (2017) ·Zbl 1397.76070号 ·doi:10.1007/s10915-016-0279-5
[23] Han,D。;Wang,X.,Cahn-Hilliard-Navier-Stokes方程的一个二阶时间、唯一可解、无条件稳定的数值格式,J Compute。物理。,290, 139-156 (2015) ·Zbl 1349.76213号 ·doi:10.1016/j.jcp.2015.02.046
[24] 黄,Q。;杨,X。;He,X.,蒙脱石-液晶流动模型的一阶、线性、解耦和能量稳定格式的数值近似,Disc。Conti公司。动态。系统-B、 232177-2192(2018)·Zbl 1408.65094号
[25] Kim,J。;Lowengrub,J.,三相流的相场建模和模拟,界面和自由边界,7435-466(2005)·Zbl 1100.35088号 ·doi:10.4171/IFB/132
[26] Lowengrub,J。;Truskinovsky,L.,《准不可压缩Cahn-Hilliard流体和拓扑转变》,R.Soc.Lond。程序。序列号。数学。物理学。工程科学。,454, 1978, 2617-2654 (1998) ·Zbl 0927.76007号 ·doi:10.1098/rspa.1998.0273
[27] Romero,I.,《非线性热力系统的热力学一致时间步长算法》,《国际数值杂志》。方法。工程师。,79, 706-732 (2009) ·Zbl 1171.80305号 ·doi:10.1002/nme.2588
[28] Schimperna,G。;Pawlow,I.,关于一类具有非线性扩散的Cahn-Hilliard模型,SIAM J.Math。分析。,45, 4, 31-63 (2013) ·Zbl 1276.35101号 ·doi:10.1137/110835608
[29] 沈,J.,王,C.,王,S,王,X.:具有Ehrlich-Schwoebel型能量的梯度流的二阶凸分裂方案:在薄膜外延中的应用。出版中,SIAM J.Numer。分析(2011)·Zbl 1247.65088号
[30] 沈杰。;Xu,J.,梯度流标量辅助变量(SAV)格式的收敛性和误差分析,SIAM J.数值分析。,56, 2895-2912 (2019) ·Zbl 1403.65047号 ·doi:10.1137/17M1159968
[31] 沈杰。;薛,J。;Yang,J.,梯度流的标量辅助变量(SAV)方法,J.Compute。物理。,353, 407-416 (2018) ·Zbl 1380.65181号 ·doi:10.1016/j.jcp.2017.10.021
[32] 沈杰。;Yang,X.,Allen-Cahn和Cahn-Hilliard方程的数值近似,Disc。Conti公司。动态。系统-A、 2016年9月28日至1691年(2010年)·Zbl 1201.65184号 ·doi:10.3934/dcds.2010.28.1669
[33] O.沃多。;Ganapathysubramanian,B.,《Cahn-Hilliard方程自适应隐式时间格式的计算效率解》,网格敏感性分析和3D等周问题,J.Comp。物理。,230, 15, 6037-6060 (2011) ·Zbl 1416.65364号 ·doi:10.1016/j.jcp.2011.04.012
[34] 徐,Z。;杨,X。;张,H。;Xie,Z.,使用稳定不变能量求积(s-IEQ)方法的各向异性Cahn-Hilliard模型的高效线性格式,计算。物理学。社区。,238, 36-49 (2019) ·Zbl 07683954号 ·doi:10.1016/j.cpc.2018.12.019
[35] Yang,X.,均聚物共混物相场模型的线性、一阶和二阶无条件能量稳定数值格式,J.Compute。物理。,327, 294-316 (2016) ·Zbl 1373.82106号 ·doi:10.1016/j.jcp.2016.09.029
[36] Yang,X.,二元流体表面活性剂系统Cahn-Hilliard相场模型的数值近似,科学杂志。计算,741533-1553(2017)·Zbl 1456.65080号 ·doi:10.1007/s10915-017-0508-6
[37] Yang,X.,各向异性相场枝晶生长模型的高效线性稳定二阶时间推进方案,计算。方法。申请。机械。工程,347316-339(2019)·Zbl 1440.74483号 ·doi:10.1016/j.cma.2018.12.012
[38] Yang,X.,Zhang,G.-D.:使用提交的不变能量求积方法对具有一般非线性势的Cahn-Hilliard和Allen-Cahn方程进行数值逼近(2018)
[39] 杨,X。;赵,J。;He,X.,粘性Cahn-Hilliard方程双曲松弛的线性、二阶和无条件能量稳定格式,使用不变能量求积方法,J.Compute。申请。数学。,343, 80-97 (2018) ·Zbl 1462.65117号 ·doi:10.1016/j.cam.2018.04.027
[40] 杨,X。;赵,J。;王,Q。;沈,J.,基于不变能量求积方法的三分量Cahn-Hilliard相场模型的数值近似,M3AS:应用科学中的数学模型和方法,271993-2030(2017)·Zbl 1393.80003号
[41] 袁,C。;Zhang,H.,水凝胶的自洽平均场模型及其数值模拟,J.Theor。计算。化学。,12, 1350048 (2013) ·doi:10.1142/S02196336135048X
[42] 赵,J。;杨,X。;龚,Y。;Wang,Q.,液晶流体动力学Q张量模型的新型线性二阶无条件能量稳定格式,计算。方法。申请。机械工程,318,803-825(2017)·兹比尔1439.76124 ·doi:10.1016/j.cma.2017.01.031
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。