鲁霍拉·塔瓦科利 Cahn-Morral方程的无条件能量稳定时间步长格式：应用于多分量调幅分解和最优空间拼接。 （英语） Zbl 1349.80047号 J.计算。物理学。 304, 441-464 (2016). 摘要：本文引入了一种无条件能量稳定的时间步长格式来求解Cahn-Morral-like方程。它是基于David Eyre的时间步进方案和Schur补码方法相结合而构造的。虽然所提出的方法是通用的，并且与齐次自由能密度函数项的选择无关，但本文特别考虑了对数和多项式能量函数。该方法用于研究多分量系统中的旋节分解和最优空间分块问题。针对后面的问题，制定了惩罚策略，以避免琐碎的解决方案。大量的数值实验证明了该方法的成功性和性能。数值结果表明，该方法收敛且能量稳定，与时间步长的选择无关。其MATLAB实现包含在附录中，用于算法的数值评估和所给出结果的再现。 引用于5文件 MSC公司： 80平方米 其他数值方法（热力学）（MSC2010） 65米70 偏微分方程初值和初边值问题的谱、配置及相关方法 80A22型 Stefan问题、相位变化等。 79年第35季度 PDE与经典热力学和传热 关键词：Cahn-Hilliard系统；MATLAB代码；多相场模型；空间划分；旋节分解；无条件能量稳定 软件：Matlab公司；拓扑_多重 PDF格式BibTeX公司 XML格式引用 \textit{R.Tavakoli}，J.计算。物理学。304、441--464（2016；Zbl 1349.80047） 全文： 内政部 参考文献： [1] Baldo，S.，《Cahn-Hilliard流体混合物相变的最小界面准则》，《Ann.Inst.Henri Poincaré》，Ana。Non Linéaire，第7、2、67-90页（1990年）·Zbl 0702.49009号 [2] 巴雷特，J.W。；Blowey，J.F.，多组分合金相分离模型的有限元近似误差界，IMA J.Numer。分析。，16, 2, 257-287 (1996) ·Zbl 0849.65069号 [3] 巴雷特，J.W。；Blowey，J.F.，具有非光滑自由能的多组分合金相分离模型的有限元近似，Numer。数学。，77, 1, 1-34 (1997) ·Zbl 0882.65129号 [4] 巴雷特，J.W。；Blowey，J.F.，具有浓度依赖迁移矩阵的多组分合金相分离模型的有限元近似，IMA J.Numer。分析。，18, 2, 287-328 (1998) ·Zbl 0899.76262号 [5] 巴雷特，J.W。；Blowey，J.F.，具有浓度相关迁移矩阵的多组分合金相分离模型的有限元近似的改进误差界，Numer。数学。，88, 2, 255-297 (2001) ·Zbl 0990.65105号 [6] 巴雷特，J.W。；布洛伊，J.F。；Garcke，H.，《关于退化Cahn-Hilliard系统的完全实用有限元逼近》，ESAIM:Math。模型。数字。分析。，35, 04, 713-748 (2001) ·Zbl 0987.35071号 [7] 贝尔托齐，A.L。；Esedoglu，S。；Gillette，A.，《使用卡恩-希利亚德方程对二值图像进行Inpainting of binary images using the Cahn-Hilliard equation》，IEEE Trans。图像处理。，16, 1, 285-291 (2007) ·Zbl 1279.94008号 [8] 巴塔查里亚，S。；Abinandanan，T.A.，三元合金中相分离的研究，Bull。马特。科学。，26, 1, 193-197 (2003) [9] 空白，L。；萨布，L。；Stoll，M.，带非局部约束的Allen-Cahn变分不等式的预处理，J.Compute。物理。，231, 16, 5406-5420 (2012) ·Zbl 1250.49007号 [10] 布洛伊，J.F。；科佩蒂，M.I.M。；Elliott，C.M.，多元合金相分离模型的数值分析，IMA J.Numer。分析。，16, 1, 111-139 (1996) ·Zbl 0857.65137号 [11] 博世，J。；斯托尔，M。；Benner，P.，使用隐式时间离散化和有限元快速求解Cahn-Hilliard变分不等式，J.Compute。物理。，262, 38-57 (2014) ·Zbl 1349.65438号 [12] Boyanova，P。；Neytcheva，M.，离散多分量Cahn-Hilliard系统的有效数值解，计算机。数学。申请。，67, 1, 106-121 (2014) ·Zbl 1381.76155号 [13] 博伊尔，F。；Minjeaud，S.，三分量Cahn-Hilliard模型的数值格式，ESAIM:Math。模型。数字。分析。，45, 04, 697-738 (2011) ·Zbl 1267.76127号 [14] Cahn，J.，关于旋节分解，金属学报。，9, 9, 795-801 (1961) [15] 卡恩，J。；Hilliard，J.，《Spinodal分解：再版》，《金属学报》。，19, 2, 151-161 (1971) [16] 卡恩，J.W。；Hilliard，J.E.，非均匀系统的自由能。I.界面自由能，J.Chem。物理。，28, 2, 258-267 (1958) ·Zbl 1431.35066号 [17] 陈，F。；Shen，J.，各向异性Cahn-Hilliard系统空间谱离散的高效能量稳定格式，Commun。计算。物理。，13, 1189-1208 (2013) ·Zbl 1388.65117号 [18] 陈，L.，三元系旋节分解的计算机模拟，金属学报。，42, 10, 3503-3513 (1994) [19] Cheng，M。；Warren，J.A.，求解相场晶体模型的有效算法，J.Compute。物理。，227, 12, 6241-6248 (2008) ·Zbl 1151.82411号 [20] Copetti，M.I.M.，三元混合物中相分离的数值实验，数学。计算。模拟。，52, 1, 41-51 (2000) [21] 考克斯·S·J。；Flikkema，E.，等面积n个受限可变形气泡的最小周长，电子。J.库姆。，17，1，第R45条pp.（2010）·Zbl 1225.05074号 [22] 考克斯，S.M。；Matthews，P.C.，刚性系统的指数时间差分，J.Compute。物理。，176, 2, 430-455 (2002) ·Zbl 1005.65069号 [23] Cueto-Felgueroso，L。；Peraire，J.，《Cahn-Hilliard和Kuramoto-Sivashinsky方程的时间自适应有限体积法》，J.Compute。物理。，227, 24, 9985-10017 (2008) ·Zbl 1153.76043号 [24] De Fontaine，D.，固体溶液中相干成分变化演变的计算机模拟（1967），西北大学，博士论文 [25] De Fontaine，D.，多组分固体溶液中的聚类和排序分析-I。稳定性标准，J.Phys。化学。固体，33，2，297-310（1972） [26] De Fontaine，D.，多组分固体溶液中的聚类和排序分析-II。波动与动力学，J.Phys。化学。固体，34，8，1285-1304（1973） [27] De Mello，E.V.L。；da Silveira Filno，Otton T.，《一维、二维和三维Cahn-Hilliard方程的数值研究》，《物理学A》，347429-443（2005） [28] 德西·J·P。；Edrees，H.H。；普莱斯，J.J。；桑德，E。；Wanner，T.，随机Cahn-Moral系统中的成核动力学，SIAM J.应用。动态。系统。，10, 2, 707-743 (2011) ·Zbl 1367.74041号 [29] Elliott，C.M。；Luckhaus，S.，具有界面自由能的多组分混合物相分离的广义扩散方程（1991），IMA预印本系列887 [30] Elsey，M。；Wirth，B.，一种简单有效的相场晶体模拟方案，ESAIM:Math。模型。数字。分析。，47, 05, 1413-1432 (2013) ·Zbl 1286.74118号 [32] Eyre，D.J.，《Cahn-Hilliard方程组》，SIAM J.Appl。数学。，53, 6, 1686-1712 (1993) ·Zbl 0853.73060号 [33] Eyre，D.J.，《无条件梯度稳定时间推进Cahn-Hilliard方程》（Mater.Res.Soc.Symp.Proc.，第529卷（1998），剑桥大学出版社），39 [34] Fife，P.C.，相分离模型及其数学，电子。J.差异。Equ.、。，2000，48，1-26（2000）·Zbl 0957.35062号 [35] Garcke，H。；雀巢公司，B。；斯汀纳，B。；Wendler，F.，Allen-Cahn体积约束系统，数学。模型方法应用。科学。，18, 08, 1347-1381 (2008) ·Zbl 1147.49036号 [36] Gómez，H。；Calo，V.M。；巴齐列夫斯，Y。；Hughes，T.，Cahn-Hilliard相场模型的等几何分析，计算。方法应用。机械。工程，197，49，4333-4352（2008）·Zbl 1194.74524号 [37] 戈麦斯，H。；Nogueira，X.，相场晶体方程的无条件能量稳定方法，计算。方法应用。机械。工程，249，52-61（2012）·Zbl 1348.74280号 [38] 戈麦斯，H。；Reali，A。；Sangalli，G.，《相场模型的精确、高效和（等）几何柔性配置方法》，J.Compute。物理。，262, 153-171 (2014) ·Zbl 1349.82084号 [39] Gräser，C。；Kornhuber，R。；Sack，U.，多组分Cahn-Hilliard系统的非光滑Schur-Newton方法，IMA J.Numer。分析。，35, 652-679 (2014) ·兹比尔1316.65086 [40] 关、珍；Lowengrub，John S。；王成；Wise，Steven M.，周期非局部Cahn-Hilliard和Allen-Cahn方程的二阶凸分裂格式，J.Compute。物理。，277, 48-71 (2014) ·Zbl 1349.65298号 [41] 关、珍；王成；Wise，Steven M.，周期非局部Cahn-Hilliard方程的收敛凸分裂格式，数值。数学。，128, 2, 377-406 (2014) ·Zbl 1304.65209号 [42] 吉伦·冈萨雷斯（Guillén-González），F。；Tierra，G.，《关于Cahn-Hilliard扩散界面模型的线性方案》，J.Compute。物理。，234, 140-171 (2013) ·Zbl 1284.35025号 [43] 美国赫赫特。；格拉纳西，L。；Pusztai，T。；Böttger，B。；阿佩尔，M。；威图谢维奇，V。；Ratke，L。；德怀尔德，J。；Froyen，L。；Camel，D.，《多元合金中的多相凝固》，马特出版社。科学。工程师R，46，1，1-49（2004） [44] M.Honjo。；Saito，Y.，使用Cahn-Hilliard方程对Fe-Cr二元和Fe-Cr-Mo三元合金中的相分离进行数值模拟，ISIJ Int.，40，9，914-919（2000） [45] Hoyt，J.，多组分系统成核的连续理论，金属学报。，38, 8, 1405-1412 (1990) [46] Hoyt，J.，规则三元合金中的线性旋节分解，金属学报。，38, 2, 227-231 (1990) [47] 胡正正；史蒂文·怀斯（Steven M.Wise）。；王成；Lowengrub，John S.，相场晶体方程的稳定有效的有限差分非线性多重网格格式，J.Compute。物理。，228, 15, 5323-5339 (2009) ·Zbl 1171.82015年 [48] Kassam，A.K。；Trefethen，L.N.，刚性偏微分方程的四阶时间步进，SIAM J.Sci。计算。，26, 4, 1214-1233 (2005) ·兹比尔1077.65105 [49] 凯·D·。；Welford，R.，《Cahn-Hilliard方程的多重网格有限元解算器》，J.Compute。物理。，212, 1, 288-304 (2006) ·Zbl 1081.65091号 [50] Kim，J。；Kang，K。；Lowengrub，J.，《Cahn-Hilliard流体的保守多重网格方法》，J.Compute。物理。，193, 2, 511-543 (2004) ·Zbl 1109.76348号 [51] Kirkaldy，J.S.，《三元旋节分解的Ginzburg-Landau处理》，J.Mater。科学。，35, 5, 1177-1180 (2000) [52] Kuwajima，T。；Saito，Y。；Suwa，Y.，铁基三元合金中的相分离动力学。二、。Fe-Cr-X（X=Mo，Cu）三元合金相分离的数值模拟，金属间化合物，11，11，1279-1285（2003） [53] Lamorgese，A.G。；莫林，D。；Mauri，R.，多相流建模的相场方法，米兰数学杂志。，79, 2, 597-642 (2011) ·Zbl 1237.82031号 [54] Lee，D。；哈，J。；Jeong，D。；Shin，J。；Yun，A。；Kim，J.，《Cahn-Hilliard方程的物理、数学和数值推导》，计算。马特。科学。，81, 216-225 (2014) [55] Lee，H。；Choi，J。；Kim，J.，《N分量Cahn-Hilliard系统的几乎无条件梯度稳定格式》，Physica A，391，4，1009-1019（2012） [56] Lee，H.G。；Kim，J.，N分量Cahn-Hilliard系统的二阶精确非线性差分格式，Physica A，387，19，4787-4799（2008） [57] 刘杰。；Dedè，L。；Evans，J.A。；Borden，M.J。；Hughes，T.，平流Cahn-Hilliard方程的等几何分析：剪切流下的旋节分解，J.Compute。物理。，242, 321-350 (2013) ·Zbl 1311.76069号 [58] Morral，J。；Cahn，J.，三元系统中的Spinodal分解，金属学报。，1037-1045年10月19日（1971年） [59] 雀巢公司，B。；Choudhury，A.，多组件系统的相场建模，Curr。操作。固态材料。科学。，15, 3, 93-105 (2011) [60] 雀巢公司，B。；温德勒，F。；塞尔泽，M。；斯汀纳，B。；Garcke，H.，具有保留体积分数的多相系统的相场模型，Phys。E版，78，1，第011604条，pp.（2008） [61] Oudet，E.，最小周长划分的（γ）收敛逼近：围绕开尔文猜想，Exp.Math。，20, 3, 260-270 (2011) ·Zbl 1261.49009号 [62] 北普罗瓦塔斯。；Elder，K.，《材料科学与工程中的相场方法》（2010），John Wiley&Sons [63] 乔，Z。；孙，Z。；Zhang，Z.，无斜率选择非线性外延生长模型二阶格式的稳定性和收敛性，数学。计算。，84, 292, 653-674 (2015) ·兹比尔1305.65186 [64] Raza，N。；Sial，S。；Butt，A.R.，使用加权Sobolev梯度对与Ginzburg-Landau泛函相关的时间演化进行数值近似，计算。数学。申请。，67, 1, 210-216 (2014) ·Zbl 1350.65069号 [65] Raza，N。；Sial，S。；Siddiqi，S.，在有限元环境中通过Sobolev梯度方法近似与Ginzburg-Landau泛函相关的时间演化，J.Compute。物理。，229, 5, 1621-1625 (2010) ·Zbl 1181.65124号 [66] Raza，N。；Sial，S。；Siddiqi，S.S.，Ginzburg-Landau泛函能量最小化相关时间演化的Sobolev梯度方法，J.Compute。物理。，228, 7, 2566-2571 (2009) ·Zbl 1166.65348号 [67] Saad，Y.，《稀疏线性系统的迭代方法》（2003），SIAM·Zbl 1002.65042号 [68] Schönlieb，C。；Bertozzi，A.，高阶修复的无条件稳定方案，Commun。数学。科学。，9413-457年9月2日（2011年）·Zbl 1216.94016号 [69] 沈杰。；王，C。；王，X。；Wise，S.M.，《具有Ehrlich-Schwoebel型能量的梯度流的二阶凸分裂方案：在薄膜外延中的应用》，SIAM J.Numer。分析。，50, 1, 105-125 (2012) ·Zbl 1247.65088号 [70] 沈杰。；Yang，X.，Allen-Cahn和Cahn-Hilliard方程的数值近似，离散Contin。动态。系统。，序列号。A、 2016年9月28日至1691年（2010年）·Zbl 1201.65184号 [71] 沈杰。；Yang，X.先生。；Yu，H.，相场移动接触线模型的高效能量稳定数值格式，J.Compute。物理。，284, 617-630 (2015) ·Zbl 1351.76184号 [72] Sial，S。；纽伯杰。；Lookman，T。；Saxena，A.，《使用Sobolev梯度实现能量最小化：相分离和排序的应用》，J.Compute。物理。，189, 1, 88-97 (2003) ·Zbl 1097.49002号 [73] 歌手罗格诺娃（Irina Loginova）；Singer，H.M.，《模拟多相材料的相场技术》，Rep.Prog。物理。，71006501（2008年） [74] Tavakoli，R.，体积约束Allen-Cahn系统和正则投影最速下降法的多材料拓扑优化，计算。方法应用。机械。工程，276，534-565（2014）·Zbl 1423.74764号 [75] Tavakoli，R.，体积约束向量值Ginzburg-Landau能量泛函最小化的计算效率方法，J.Compute。物理。，295, 355-378 (2015) ·Zbl 1349.65199号 [76] 塔瓦科利，R。；Mohseni，S.M.，多材料拓扑优化问题的交替有源相位算法：115线MATLAB实现，Struct。多磁盘。最佳。，49, 4, 621-642 (2014) [77] 蒂埃拉，G。；Guillén-González，F.，求解Cahn-Hilliard方程的数值方法及其对相关能量模型的适用性，Arch。计算。方法工程，22，269-289（2015）·Zbl 1348.82080号 [78] Trefethen，Lloyd N.，MATLAB中的谱方法，第10卷（2000年），SIAM·兹比尔0953.68643 [79] Vollmayr-Lee，B.P。；Rutenberg，A.D.，《具有无条件稳定时间步长的快速准确粗化模拟》，Phys。E版，68、6，第066703条，pp.（2003） [80] 怀斯，S。；Kim，J。；Lowengrub，J.，用自适应非线性多重网格法求解正则化强各向异性Cahn-Hilliard方程，J.Compute。物理。，226, 1, 414-446 (2007) ·Zbl 1310.82044号 [81] 怀斯，S.M。；王，C。；Lowengrub，J.S.，相场晶体方程的能量稳定收敛有限差分格式，SIAM J.Numer。分析。，47, 3, 2269-2288 (2009) ·兹比尔1201.35027 [82] O.沃多。；Ganapathysubramanian，B.，《Cahn-Hilliard方程的计算有效解：自适应隐式时间方案、网格敏感性分析和3D等周问题》，J.Compute。物理。，230, 15, 6037-6060 (2011) ·Zbl 1416.65364号 [83] 张，L。；Tonks，M.R。；D.加斯顿。；彼得森，J.W。；安德斯，D。；Millett，P.C。；Biner，B.S.，《用于求解卡恩-希利亚德方程的C0和C1元素之间的定量比较》，J.Compute。物理。，236, 74-80 (2013) [84] 张，Z。；乔，Z.，《卡恩-希利亚德方程的自适应时间步进策略》，Commun。计算。物理。，11, 1261-1278 (2012) ·Zbl 1388.65072号 此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配，可能包含数据转换错误。在某些情况下，zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献，而不要求完整或完全匹配。