裴帅超;侯燕仁;李，齐 修正相场晶体方程的线性二阶无条件能量稳定格式及其误差估计。 （英语） Zbl 1524.65384号 计算。数学。应用。 103, 104-126 (2021). 摘要：我们对一种线性、二阶、无条件能量稳定的半离散时间步进方案进行了误差估计，该方案基于拉格朗日乘子方法，也可以被视为不变能量平方方法的特例，用于修正的相场晶体方程。在每个时间步长，该方案只需要解几个线性方程组。给出了严格的证明，证明了该格式的唯一可解性、质量守恒性和无条件能量稳定性。在2D和3D中进行了各种数值实验，以验证所提出的数值策略的准确性、无条件的能量稳定性和质量守恒。 引用于8文件 理学硕士： 6500万06 含偏微分方程初值和初边值问题的有限差分方法 65个M12 含偏微分方程初值和初边值问题数值方法的稳定性和收敛性 65M60毫米 涉及偏微分方程初值和初边值问题的有限元、Rayleigh-Ritz和Galerkin方法 65M70型 偏微分方程初值和初边值问题的谱、配置及相关方法 35问题35 与流体力学相关的PDE 65岁15岁 涉及PDE的初值和初边值问题的误差界 65纳米35 偏微分方程边值问题的谱、配置及相关方法 80A22型 Stefan问题、相位变化等。 关键词：修正相场方程;无条件能量稳定性;误差估计 PDF格式BibTeX公司 XML格式引用 \textit{S.Pei}等人，计算。数学。申请。103、104-126（2021；Zbl 1524.65384） 全文： 内政部 参考文献： [1] Akrivis，G。；李，B。；Li，D.，Allen-Cahn和Cahn-Hilliard方程的能量衰减外推RK-SAV方法，SIAM J.Sci。计算。，41、6、A3703-A3727（2019）·Zbl 1435.65141号 [2] 巴迪亚，S。；吉伦·冈萨雷斯（Guillén-González），F。；Gutiérrez-Santacreu，J.V.，使用鞍点结构的向列相液晶流的有限元近似，J.Compute。物理。，230, 4, 1686-1706 (2011) ·Zbl 1211.82056号 [3] Baskaran，A。；胡，Z。；Lowengrub，J.S。；王，C。；怀斯，S.M。；周，P.，修正相场晶体方程的能量稳定且有效的非线性有限差分多重网格格式，J.Compute。物理。，250, 270-292 (2013) ·Zbl 1349.65265号 [4] Baskaran，A。；Lowengrub，J.S。；王，C。；Wise，S.M.，修正相场晶体方程二阶凸分裂格式的收敛性分析，SIAM J.Numer。分析。，51, 5, 2851-2873 (2013) ·Zbl 1401.82046号 [5] 布埃诺，J。；斯塔罗杜莫夫，I。；戈麦斯，H。；Galenko，P。；Alexandrov，D.，修正相场晶体方程预测的三维结构，计算。马特。科学。，111, 310-312 (2016) [6] 陈，H。；毛，J。；Shen，J.，梯度流标量辅助变量有限元格式的最佳误差估计，数值。数学。，145, 167-196 (2020) ·Zbl 1440.65192号 [7] 陈，N。；王，C。；Wise，S.，一类双稳态梯度流的全局时间Gevrey正则解，离散Contin。动态。系统。，序列号。B、 1689-1711年6月21日（2016）·Zbl 1350.35050号 [8] 陈，R。；杨，X。；Zhang，H.，蒙脱石-a液晶流体动力学模型的二阶线性无条件能量稳定格式，SIAM J.Sci。计算。，39，6，A2808-A2833（2017）·Zbl 1388.76201号 [9] Chen，W。；康德，S。；王，C。；王，X。；Wise，S.M.，无斜率选择的薄膜模型的线性能量稳定方案，科学杂志。计算。，52, 546-562 (2012) ·Zbl 1326.76007号 [10] Cheng，K。；王，C。；Wise，S.M.，方相场晶体方程的能量稳定的BDF2傅里叶伪谱数值格式，Commun。计算。物理。，26, 1335-1364 (2019) ·Zbl 1518.65115号 [11] Cheng，Q。；Shen，J.，多标量辅助变量（MSAV）方法及其在相场囊泡膜模型中的应用，SIAM J.Sci。计算。，40、6、A3982-A4006（2018）·Zbl 1403.65045号 [12] Cheng，Q。；沈，J。；Yang，X.，使用SAV方法的外延薄膜生长模型的高效和准确的数值格式，J.Sci。计算。，78, 1467-1487 (2019) ·Zbl 1444.76020号 [13] Dehghan，M。；Mohammadi，V.，通过全局和局部无网格方法对相场晶体（PFC）和修正相场晶体模型进行的数值模拟，计算。方法应用。机械。工程，298453-484（2016）·兹比尔1423.76321 [14] 迪格尔，A.E。；王，C。；Wise，S.M.，Cahn-Hilliard方程二阶混合有限元方法的稳定性和收敛性，IMA J.Numer。分析。，36, 1867-1897 (2016) ·Zbl 1433.80005 [15] 杜琪。；Feng，X.，几何运动界面的相场方法及其数值近似，（Bonito，A.；Nochetto，R.H.，《几何偏微分方程-第一部分几何偏微分方程式-第一部分，数值分析手册》，第21卷（2020年），Elsevier），425-508·兹比尔1455.35276 [16] Elder，K.R。；Grant，M.，《使用相场晶体模拟非平衡处理中的弹性和塑性变形》，Phys。E版，70，第051605条pp.（2004） [17] Elder，K.R。；Katakowski，M。；Haataja，M。；Grant，M.，《晶体生长弹性建模》，Phys。修订稿。，88，第245701条，pp.（2002） [18] Elder，K.R。；普罗瓦塔斯，N。；Berry，J。；斯特凡诺维奇，P。；Grant，M.，相场晶体建模和经典密度泛函冻结理论，Phys。B版，75，第064107条，pp.（2007） [19] Elliott，C.M。；Stuart，A.M.，离散双线性抛物方程的全局动力学，SIAM J.Numer。分析。，30, 6, 1622-1663 (1993) ·Zbl 0792.65066号 [20] Eyre，D.J.，无条件梯度稳定时间推进Cahn-Hilliard方程，MRS Proc。，529, 39-46 (1998) [21] 冯·W。；萨尔加多，A.J。；王，C。；Wise，S.M.，涉及p-Laplacian项的非线性椭圆方程的预处理最速下降法，J.Compute。物理。，334, 45-67 (2017) ·Zbl 1375.35149号 [22] 冯，X。；Tang，T。；Yang，J.，相场模型的稳定曲柄-尼科尔森/阿达姆-巴霍斯方案，东亚应用杂志。数学。，3, 1, 59-80 (2015) ·Zbl 1302.65224号 [23] Galenko，P.K。；戈麦斯，H。；Kropotin，N.V.公司。；Elder，K.R.，双曲相场晶体方程的无条件稳定方法和数值解，Phys。E版，88，第013310条，pp.（2013） [24] Golovin，A.A。；Nepomnyashchy，A.A.，正方形和六角形模式中的偏差，物理学。E版，67，第056202条，pp.（2003） [25] 戈麦斯，H。；Nogueira，X.，相场晶体方程的无条件能量稳定方法，计算。方法应用。机械。工程，249-252，52-61（2012）·Zbl 1348.74280号 [26] 龚，Y。；赵，J。；Wang，Q.，梯度流模型的任意高阶无条件能量稳定SAV格式，计算。物理学。社区。，249，第107033条pp.（2020）·Zbl 07678495号 [27] 龚，Y。；赵，J。；Wang，Q.，热力学一致梯度流模型的任意高阶无条件能量稳定格式，SIAM J.Sci。计算。，42、1、B135-B156（2020）·兹比尔1447.65099 [28] 格拉塞利，M。；Pierre，M.，修正相场晶体方程的能量稳定和收敛有限元格式，ESAIM:M2AN，50，5，1523-1560（2016）·Zbl 1358.82025号 [29] 关，Z。；王，C。；Wise，S.M.，周期非局部Cahn-Hilliard方程的收敛凸分裂格式，数值。数学。，128, 377-406 (2014) ·Zbl 1304.65209号 [30] 吉伦·冈萨雷斯（Guillén-González），F。；Tierra，G.，《关于Cahn-Hilliard扩散界面模型的线性方案》，J.Compute。物理。，234, 140-171 (2013) ·Zbl 1284.35025号 [31] 郭，R。；Xu，Y.，修正相场晶体方程的高阶自适应时间步长策略和局部间断Galerkin方法，Commun。计算。物理。，24, 1, 123-151 (2018) ·Zbl 1488.65420号 [32] Han，D。；Brylev，A。；杨，X。；Tan，Z.，两相不可压缩流相场模型的二阶全离散能量稳定格式的数值分析，科学杂志。计算。，70, 965-989 (2017) ·Zbl 1397.76070号 [33] 何毅。；刘，Y。；Tang，T.，关于Cahn-Hilliard方程的大时间步长方法，应用。数字。数学。，57, 5, 616-628 (2007) ·Zbl 1118.65109号 [34] 海伍德，J.G。；Rannacher，R.，非平稳Navier-Stokes问题的有限元近似。第四部分：二阶时间离散化的误差分析，SIAM J.Numer。分析。，27, 2, 353-384 (1990) ·Zbl 0694.76014号 [35] 胡，Z。；怀斯，S。；王，C。；Lowengrub，J.，相场晶体方程的稳定有效的非线性有限差分多重网格格式，J.Compute。物理。，228, 15, 5323-5339 (2009) ·兹比尔1171.82015 [36] 黄，Q.-A。；蒋伟（Jiang，W.）。；Yang，J.Z.，用于模拟柔性基底上弹性薄膜褶皱现象的无条件能量稳定方案，J.Compute。物理。，388, 123-143 (2019) ·Zbl 1452.74117号 [37] 寇，J。；Sun，S。；Wang，X.，多组分两相可压缩流动的线性解耦能量稳定数值方法，SIAM J.Numer。分析。，56, 6, 3219-3248 (2018) ·Zbl 1404.65226号 [38] Lee，H.G。；Shin，J。；Lee，J.-Y.，修正相场晶体方程的一阶和二阶能量稳定方法，计算。方法应用。机械。工程，321，1-17（2017）·Zbl 1439.74240号 [39] 李，D。；乔，Z。；Tang，T.，表征相场方程半隐式Fourier谱方法的稳定尺寸，SIAM J.Numer。分析。，54, 3, 1653-1681 (2016) ·Zbl 1346.35162号 [40] 李，H。；Ju，L。；张，C。；Peng，Q.，具有Peng-Robinson状态方程的扩散界面模型的无条件能量稳定线性格式，科学杂志。计算。，75, 993-1015 (2018) ·Zbl 1398.65305号 [41] 李强。；梅，L。；杨，X。；Li，Y.，修正相场晶体方程无条件能量稳定的高效数值格式，高级计算。数学。，45, 3, 1551-1580 (2019) ·Zbl 1420.82021 [42] 李，X。；Shen，J.，相场晶体方程SAV-Fourier谱方法的稳定性和误差估计，高级计算。数学。，46, 48 (2020) ·Zbl 1442.65294号 [43] 李，X。；Shen，J.，修正相场晶体方程的高效线性和无条件能量稳定格式，Sci。中国数学。，64 (2021) [44] 李，X。；沈，J。；Rui，H.，梯度流SAV块中心有限差分法的能量稳定性和收敛性，数学。计算。，88, 2047-2068 (2019) ·Zbl 1422.65165号 [45] 林，F。；何，X。；Wen，X.，新Allen-Cahn型方形相场晶体模型的快速、无条件能量稳定大时间步进法，应用。数学。莱特。，98, 248-255 (2019) ·Zbl 1433.65250号 [46] 林，L。；杨，Z。；Dong，S.，基于辅助能量变量的不可压缩Navier-Stokes方程的数值逼近，J.Compute。物理。，388, 1-22 (2019) ·Zbl 1452.76093号 [47] 刘，C。；沈，J。；Yang，X.，向列相液晶流中缺陷运动的动力学：建模和数值模拟，Commun。计算。物理。，2, 6, 1184-1198 (2007) ·Zbl 1164.76422号 [48] 刘，Z。；Li，X.，相场晶体方程的高效修正稳定不变能量求积方法，数值。算法，85107-132（2019）·Zbl 1446.65071号 [49] 刘，Z。；Li，X.，相场模型的指数标量辅助变量（E-SAV）方法及其显式计算，SIAM J.Sci。计算。，42、3、B630-B655（2020年）·Zbl 1447.65057号 [50] 劳埃德·D.J.B。；桑斯特德，B。；Avitabile，D。；Champneys，A.R.，平面Swift-Hohenberg方程的局部六边形图案，SIAM J.Appl。动态。系统。，7, 3, 1049-1100 (2008) ·Zbl 1168.35311号 [51] 普罗瓦塔斯，N。；Dantzig，J.A。；Athreya，B。；Chan，P。；斯特凡诺维奇，P。；Goldenfeld，N。；Elder，K.R.，《在微观结构演化的多尺度建模中使用相场晶体方法》，JOM，59，7，83-90（2007） [52] Shen，J.，通过相场方法对两相不可压缩流进行建模和数值近似，（（2011）），147-195 [53] 沈，J。；王，C。；王，X。；Wise，S.M.，《具有Ehrlich-Schwoebel型能量的梯度流的二阶凸分裂方案：在薄膜外延中的应用》，SIAM J.Numer。分析。，50, 1, 105-125 (2012) ·Zbl 1247.65088号 [54] 沈，J。；Xu，J.，梯度流标量辅助变量（SAV）格式的收敛性和误差分析，SIAM J.Numer。分析。，56, 5, 2895-2912 (2018) ·Zbl 1403.65047号 [55] 沈，J。；徐，J。；Yang，J.，梯度流的标量辅助变量（SAV）方法，J.Compute。物理。，353, 407-416 (2018) ·Zbl 1380.65181号 [56] 沈，J。；徐，J。；Yang，J.，梯度流的一类新型高效稳健能量稳定方案，SIAM Rev.，61，3，474-506（2019）·Zbl 1422.65080号 [57] 沈，J。；Yang，X.，Allen-Cahn和Cahn-Hilliard方程的数值近似，离散Contin。动态。系统。，28, 1669-1691 (2010) ·Zbl 1201.65184号 [58] Shin，J。；Lee，H.G。；Lee，J.-Y.，基于相场晶体方程新凸分裂的一阶和二阶数值方法，J.Compute。物理。，327, 519-542 (2016) ·Zbl 1373.82097号 [59] 斯特凡诺维奇，P。；Haataja，M。；Provatas，N.，弹性相互作用的相场晶体，物理。修订稿。，96，第225504条，pp.（2006） [60] 斯特凡诺维奇，P。；Haataja，M。；Provatas，N.，纳米晶材料变形和塑性的相场晶体研究，Phys。E版，80，第046107条，pp.（2009） [61] 斯威夫特，J。；Hohenberg，P.C.，对流不稳定性的流体动力学波动，物理学。A版，第15页，第319-328页（1977年） [62] Tang，T。；乔，Z.，相场方程的有效数值方法，科学。罪。，数学。，50775-794（2020），（中文）·Zbl 1499.65581号 [63] 蒂尔拉，G。；Guillén-González，F.，求解Cahn-Hilliard方程的数值方法及其对相关能量模型的适用性，Arch。计算。方法工程，22，269-289（2015）·兹比尔1348.82080 [64] 王，C。；Wise，S.M.，三维修正相场晶体方程的整体光滑解，方法应用。分析。，17, 2, 191-212 (2010) ·Zbl 1217.35040号 [65] 王，C。；Wise，S.M.，修正相场晶体方程的能量稳定收敛有限差分格式，SIAM J.Numer。分析。，49, 3, 945-969 (2011) ·Zbl 1230.82005年 [66] 王，M。；黄，Q。；Wang，C.，方相场晶体方程的二阶精确标量辅助变量（SAV）数值方法，J.Sci。计算。，88, 33 (2021) ·Zbl 1497.65196号 [67] 怀斯，S.M。；王，C。；Lowengrub，J.S.，相场晶体方程的能量稳定收敛有限差分格式，SIAM J.Numer。分析。，47, 3, 2269-2288 (2009) ·Zbl 1201.35027号 [68] 徐，C。；Tang，T.，外延生长模型大时间步进方法的稳定性分析，SIAM J.Numer。分析。，44, 4, 1759-1779 (2006) ·Zbl 1127.65069号 [69] 徐，Z。；杨，X。；张，H。；Xie，Z.，使用稳定不变能量二次化（S-IEQ）方法的各向异性Cahn Hilliard模型的有效和线性格式，计算。物理学。社区。，238, 36-49 (2019) ·Zbl 07683954号 [70] Yang，X.，均聚物共混物相场模型的线性、一阶和二阶无条件能量稳定数值格式，J.Compute。物理。，327, 294-316 (2016) ·Zbl 1373.82106号 [71] Yang，X.，二元流体-表面活性剂系统Cahn-Hilliard相场模型的数值近似，J.Sci。计算。，74, 1533-1553 (2018) ·Zbl 1456.65080号 [72] Yang，X.，各向异性相场枝晶生长模型的高效线性稳定二阶时间推进方案，计算。方法应用。机械。工程，347316-339（2019）·Zbl 1440.74483号 [73] 杨，X。；Han，D.，相场晶体模型的线性一阶和二阶无条件能量稳定格式，J.Compute。物理。，330, 1116-1134 (2017) ·Zbl 1380.65209号 [74] 杨，X。；Ju，L.，相场弹性弯曲能量模型的无条件能量稳定性高效线性格式，计算。方法应用。机械。工程，315691-712（2017）·Zbl 1439.74165号 [75] 杨，X。；Ju，L.，二元流体-表面活性剂相场模型的线性和无条件能量稳定格式，计算。方法应用。机械。工程，3181005-1029（2017）·Zbl 1439.76029号 [76] 杨，X。；Yu，H.，使用不变能量二次化方法的相场移动接触线模型的有效二阶无条件稳定方案，SIAM J.Sci。计算。，40、3、B889-B914（2018）·Zbl 1395.65095号 [77] 杨，X。；Zhang，G.-D.，求解具有一般非线性势的Cahn-Hilliard和Allen-Cahn方程的不变能量二次化（IEQ）格式的收敛性分析，J.Sci。计算。，82, 55 (2020) ·Zbl 1434.65222号 [78] 杨，X。；赵，J。；Wang，Q.，基于不变能量求积方法的分子束外延生长模型的数值近似，J.Compute。物理。，333, 104-127 (2017) ·Zbl 1375.82121号 [79] 张杰。；Yang，X.，具有强非线性空位势的修正相场晶体模型的高效二阶无条件稳定时间推进数值格式，计算。物理学。社区。，245，第106860条pp.（2019）·Zbl 07674875号 [80] 张杰。；赵，J。；Gong，Y.，Cahn-Hilliard型方程全离散不变能量求积方案的误差分析，J.Compute。申请。数学。，372，第112719条pp.（2020）·Zbl 1440.65257号 [81] 赵，J。；杨，X。；龚，Y。；Wang，Q.，液晶流体动力学Q张量模型的新型线性二阶无条件能量稳定格式，计算。方法应用。机械。工程，318，803-825（2017）·Zbl 1439.76124号 [82] 朱，G。；寇，J。；姚，J。；李，A。；Sun，S.，《含有可溶性表面活性剂的相场移动接触线模型》，J.Compute。物理。，405，第109170条pp.（2020）·Zbl 1453.76146号 [83] 朱，J。；陈立群。；沈，J。；Tilare，V.，变迁移率Cahn-Hilliard方程的粗化动力学：半隐式傅里叶谱方法的应用，Phys。版本E，60，4，3564-3572（1999） 此参考列表基于出版商或数字数学图书馆提供的信息。它的项目与zbMATH标识符启发式匹配，并且可能包含数据转换错误。在某些情况下，zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献，而不要求完整或完全匹配。