×

相场晶体模型的线性二阶能量稳定方案。 (英语) Zbl 1432.82028号

摘要:本文基于凸分裂格式和相位场晶体方程的半隐式光谱延迟校正(SISDC)方法,提出了一种线性、无条件能量稳定的二阶(时间)数值格式。我们使用的凸分裂格式是线性的、唯一可解的和无条件能量稳定的,但它是一阶的,因此我们采用SISDC方法来提高收敛速度。该方案继承了凸分裂方案的优点,从而在每个时间步长上都得到了线性方程组,易于实现。我们还证明了该方案是无条件弱能量稳定的,并且在时间上具有二阶精度。数值实验验证了所提出的数值策略的准确性、有效性和能量稳定性。

MSC公司:

82M10个 有限元、伽辽金及相关方法在统计力学问题中的应用
65M60毫米 涉及偏微分方程初值和初边值问题的有限元、Rayleigh-Ritz和Galerkin方法
65个M12 含偏微分方程初值和初边值问题数值方法的稳定性和收敛性
82C26型 统计力学中的动态和非平衡相变(一般)
82D25个 晶体统计力学
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] Baskaran,A。;Lowengrub,J.S。;王,C。;Wise,S.M.,修正相场晶体方程二阶凸分裂格式的收敛性分析,SIAM J.Numer。分析。,51, 5, 2851-2873 (2013) ·Zbl 1401.82046号
[2] 洛杉矶卡法雷利。;Muler,N.E.,Cahn-Hilliard方程解的界限,Arch。定额。机械。分析。,133、2、129-144(1995年12月)·Zbl 0851.35010号
[3] Chen,L.,(i)FEM:MATLAB中的集成有限元方法包(2009),加州大学欧文分校,技术报告
[4] Cheng,M。;Warren,J.A.,求解相场晶体模型的有效算法,J.Compute。物理。,227, 12, 6241-6248 (2008) ·Zbl 1151.82411号
[5] Condette,N。;梅尔彻,C。;Süli,E.,具有二次增长双阱势的模式形成非线性演化方程的谱近似,数学。计算。,80, 273, 205-223 (2011) ·Zbl 1209.65066号
[6] Dong,L。;冯·W。;王,C。;怀斯,S.M。;Zhang,Z.,三维相场晶体方程二阶数值格式的收敛性分析和数值实现,计算。数学。申请。,1912年6月75日至1928年(2018年)·Zbl 1409.82014年8月
[7] Dutt,A。;格林加德。;Rokhlin,V.,《常微分方程的谱延迟校正方法》,BIT,40,2,241-266(2000)·Zbl 0959.65084号
[8] Elder,K.R。;Grant,M.,《使用相场晶体模拟非平衡处理中的弹性和塑性变形》,Phys。E版,70,第051605条pp.(2004)
[9] Elder,K.R。;Katakowski,M。;Haataja,M。;Grant,M.,晶体生长中的弹性建模,Phys。修订稿。,88,第245701条,pp.(2002)
[10] Evans,L.C.,偏微分方程(2010),美国数学学会·Zbl 1194.35001号
[11] Eyre,D.J.,无条件梯度稳定时间推进Cahn-Hilliard方程,MRS Proc。,529, 39 (1998)
[12] X·冯。;Tang,T。;Yang,J.,使用p-自适应谱延迟校正方法进行相场问题的长期数值模拟,SIAM J.Sci。计算。,37,1,A271-A294(2015)·兹比尔1327.65161
[13] 格拉斯纳,K。;Orizaga,S.,《提高梯度流方程凸度分裂方法的精度》,J.Compute。物理。,315, 52-64 (2016) ·Zbl 1349.65515号
[14] 戈麦斯,H。;Nogueira,X.,相场晶体方程的无条件能量稳定方法,计算。方法应用。机械。工程,249-252,52-61(2012)·Zbl 1348.74280号
[15] 关,Z。;王,C。;Wise,S.M.,周期非局部Cahn-Hilliard方程的收敛凸分裂格式,数值。数学。,128, 2, 377-406 (2014) ·Zbl 1304.65209号
[16] 郭,R。;Xu,Y.,相场晶体方程的局部间断Galerkin方法和高阶半隐式格式,SIAM J.Sci。计算。,38,1,A105-A127(2016)·Zbl 1330.65149号
[17] 海伍德,J.G。;Rannacher,R.,非平稳Navier-Stokes问题的有限元近似。第四部分:二阶时间离散化的误差分析,SIAM J.Numer。分析。,27, 2, 353-384 (1990) ·Zbl 0694.76014号
[18] 胡,Z。;怀斯,S。;王,C。;Lowengrub,J.,相场晶体方程的稳定有效的非线性有限差分多重网格格式,J.Compute。物理。,228, 15, 5323-5339 (2009) ·Zbl 1171.82015年
[19] Kessler博士。;诺切托,R.H。;Schmidt,A.,Allen-Cahn问题的后验误差控制:规避Gronwall不等式,ESAIM:数学。模型。数字。分析。,38, 1, 129-142 (2004) ·Zbl 1075.65117号
[20] 李,S。;郑浩。;Layton,W.J.,非平稳Stokes-Darcy模型的不同子域时间步长解耦方法,数值。方法部分差异。Equ.、。,29, 2, 549-583 (2013) ·Zbl 1364.76096号
[21] Minion,M.L.,基于光谱延迟修正的不可压缩流半隐式投影方法,应用。数字。数学。,48, 3, 369-387 (2004) ·Zbl 1035.76040号
[22] 北普罗瓦塔斯。;Dantzig,J.A。;Athreya,B。;Chan,P。;斯特凡诺维奇,P。;Goldenfeld,N。;Elder,K.R.,《在微观结构演化的多尺度建模中使用相场晶体方法》,JOM,59,7,83-90(2007)
[23] 沈杰。;Tang,T。;Wang,L.-L.,谱方法:算法、分析和应用,计算数学中的Springer系列,第41卷(2011),Springer:Springer Berlin,Heidelberg·兹比尔1227.65117
[24] 沈杰。;王,C。;王,X。;Wise,S.M.,《具有Ehrlich-Schwoebel型能量的梯度流的二阶凸分裂方案:在薄膜外延中的应用》,SIAM J.Numer。分析。,50, 1, 105-125 (2012) ·Zbl 1247.65088号
[25] J.Shen,J.Xu,J.Yang,梯度流的一类新型高效稳健能量稳定方案,arXiv e-prints,2017。;J.Shen,J.Xu,J.Yang,梯度流的一类新型高效稳健能量稳定方案,arXiv电子版,2017年。
[26] 沈杰。;徐,J。;Yang,J.,梯度流的标量辅助变量(SAV)方法,J.Compute。物理。,353, 407-416 (2018) ·Zbl 1380.65181号
[27] 沈杰。;Yang,X.,Allen-Cahn和Cahn-Hilliard方程的数值近似,离散Contin。动态。系统。,28, 1669-1691 (2010) ·Zbl 1201.65184号
[28] Shin,J。;Lee,H.G。;Lee,J.-Y.,基于相场晶体方程新凸分裂的一阶和二阶数值方法,J.Compute。物理。,327, 519-542 (2016) ·Zbl 1373.82097号
[29] 斯威夫特,J。;Hohenberg,P.C.,对流不稳定性的流体动力学波动,物理学。A版,第15页,第319-328页(1977年)
[30] 维格纳尔,P。;达尔星。;D.布朗。;科利尔,N。;Calo,V.,相场晶体方程的能量稳定凸分裂,计算。结构。,158, 355-368 (2015)
[31] 怀斯,S.M。;王,C。;Lowengrub,J.S.,相场晶体方程的能量稳定收敛有限差分格式,SIAM J.Numer。分析。,47, 3, 2269-2288 (2009) ·兹比尔1201.35027
[32] 徐,C。;Tang,T.,外延生长模型大时间步进方法的稳定性分析,SIAM J.Numer。分析。,44, 4, 1759-1779 (2006) ·Zbl 1127.65069号
[33] 杨,X。;Han,D.,相场晶体模型的线性一阶和二阶无条件能量稳定格式,J.Compute。物理。,330, 1116-1134 (2017) ·Zbl 1380.65209号
[34] 杨,X。;Ju,L.,相场弹性弯曲能量模型的具有无条件能量稳定性的有效线性格式,计算。方法应用。机械。工程,315691-712(2017)·Zbl 1439.74165号
[35] 张,Z。;马云(Ma,Y.)。;乔,Z.,求解相场晶体模型的自适应时间步进策略,J.Compute。物理。,249, 204-215 (2013) ·Zbl 1305.82009年
此参考列表基于出版商或数字数学图书馆提供的信息。它的项目与zbMATH标识符启发式匹配,并且可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。