×

具有奇异位势和动态边界条件的Caginalp系统的双分裂格式。 (英语) Zbl 1459.65184号

近年来,考虑到Caginalp相场模型(CPFM)及其变体在广泛领域的适用性,人们对其进行了大量研究。本文致力于研究具有动态边界条件和奇异势的Caginalp相场系统。对于这种典型的PFM,作者考虑了一个非线性抛物线耦合系统,该系统控制(相对)温度和序参数的演化,预测长时间行为。现在,从物理角度来看,奇异势也很重要;特别是,作者考虑到以下热力学相关对数势。为了从数学上证明,与正则势相反,这种奇异势具有分离性质,可以确保序参数严格保持在(-1)和1之间,正如从物理角度所期望的那样。本文作者对具有动态边界条件和奇异势,特别是对数势的Caginalp系统感兴趣。它们证明了解的存在唯一性及其正则性。本研究的主要内容在于证明序参数u与势的奇异值是严格分离的。在大多数与相场方法(PFM)相关的文献中,已经证明了各种类型模型的解的存在性和唯一性,这些模型在为实现特定目标而定制的假设中存在差异。现在,从物理角度来看,奇异势也很重要;因此,作者在这里关注奇异势。作者进行的文献调查表明,
深入分析了具有奇异位势和动态边界条件的Caginalp系统的适定性和长期行为。
从数值角度来看,对具有对数函数的PFM进行了正则解和奇异解的计算,但仅限于奇异时间。
据作者所知,直到撰写手稿时,涉及对数势和动态边界条件的Cahn-Hilliard型问题奇异解的计算似乎还没有得到解决。
因此,作者在本文中的研究方向是提出并分析一个方案,该方案允许我们计算问题的奇异解,即使在奇异性发生之后。他们方法的一个基本思想是利用与问题相关的能量。Caginalp系统的双重分裂方案是基于两个思想:时间分裂,它解耦了问题中每个时间步长的方程的求解,以及能量的凸分裂方案(CSS),为证明所提方案的无条件唯一可解性奠定了基础。如果时间步长极小,则可以获得其离散版本(计算模型)。本文的主要结果是以命题的形式表述的:通过让时间步长趋于0并使用单调性论证,时间半离散解收敛到问题的能量解。本文的一部分致力于严格证明和引入能量解概念所必需的假设、符号和函数设置。在适当的假设下,从理论上证明了与能量耗散解的存在性和收敛性以及格式的无条件唯一可解性有关的两个主要定理,并用数值方法证明了这两个定理是本文的主要结果。单独的章节专门关注CPFM一维平稳奇异解的分析和数值计算,以及Caginalp系统在二维空间中正则和奇异解的数值计算。通过大量的数值实验验证了该方法的准确性和有效性。该半离散方法可以以较低的实现成本方便地求解变系数的一维、二维和三维问题。新方法简单有效。给出了数值实验结果,并对其精度进行了讨论和比较。模拟是通过FreeFem++软件进行的。对于优化问题,使用了ff-NLopt包中的截断牛顿算法。意见和评论
这是一个重要的研究领域,直接应用于材料科学的许多领域。本文的目的是研究具有对数势和动态边界条件的Caginalp相场模型在序参量和温度下的长期行为。
许多数值分析人员开发了时间分裂方法,将复杂的含时偏微分方程分解为简单的方程组,然后可以通过数值方法在时间步长的分数上分别求解这些方程组。在本文中,时间分裂在每一步都解耦了问题中方程的求解。
仅仅开发速度更快的计算机是不够的,还需要数值分析(快速有效的算法)。双分裂格式可能为设计精确有效的相场建模数值格式带来一个很有前途的方向。
在计算偏微分方程的数值解时,经常使用一种或另一种形式的分裂方法。Caginalp相场模型是一个含时(演化)偏微分方程系统,如作者所解释的,双重分裂还有一个额外的优点。
该方法考虑了三个主要方面。首先是时间分割方法相对于非分割方法的准确性和效率。第二是分裂方法的稳定性。最后,该方法的鲁棒性。提议的方案符合所有标准。
作者对Caginalp系统进行了二维数值模拟。使用FreeFem++软件进行模拟。
误差可能会引发数值不稳定性,从而破坏收敛性。在理论和数值实验中处理了解的存在唯一性和正则性。表1显示了最终时间T=0:5时,双分裂方案的解与线性隐式方案的解之间的L2–误差。
指出对于正则解,线性隐式格式比双分裂格式快得多。这如表2所示。然而,在作者看来,使用二阶方法可以减少计算时间,该方法不仅涉及目标函数的Hessian,而且还涉及梯度
通常需要对长期解决方案的进化行为进行可视化,以便对模型进行预测。使用Matlab软件进行可视化。
通过与线性隐式格式计算的CPFM解的比较,对该格式进行了数值验证。

MSC公司:

65M60毫米 涉及偏微分方程初值和初边值问题的有限元、Rayleigh-Ritz和Galerkin方法
65N30型 含偏微分方程边值问题的有限元、Rayleigh-Ritz和Galerkin方法
65K10码 数值优化和变分技术
35K67型 奇异抛物方程
80A22型 Stefan问题、相变等。

软件:

自由Fem++
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] P.F.Antonietti;格拉塞利(M.Grasselli);S.Stangalino;M.Verani,具有动态边界条件的线性抛物问题的间断Galerkin近似,科学杂志。计算。,66, 1260-1280 (2016) ·Zbl 1376.65126号 ·doi:10.1007/s10915-015-0063-y
[2] H.Attouch,函数和算子的变分收敛性《应用数学系列》,皮特曼(高级出版计划),马萨诸塞州波士顿,1984年·Zbl 0561.49012号
[3] V.Barbu,Banach空间中的非线性半群和微分方程布加勒斯特罗姆尼亚共和国社会党编辑学院;诺德霍夫国际出版公司,莱顿,1976年·Zbl 0328.47035号
[4] S.Bartels;R.Müller,具有对数势的Allen-Cahn和Cahn-Hilliard方程近似的误差控制,数值。数学。,119, 409-435 (2011) ·兹比尔1241.65075 ·doi:10.1007/s00211-011-0389-9
[5] N.巴坦古纳;M.Pierre,Caginalp相场系统时间分裂近似下指数吸引子的收敛性,Commun。纯应用程序。分析。,17, 1-19 (2018) ·Zbl 1375.37167号 ·doi:10.3934/cpaa.2018001年
[6] F.Boyer;F.Nabet,带动态边界条件的Cahn-Hilliard/Stokes相场模型的DDFV方法,ESAIM数学。模型。数字。分析。,51, 1691-1731 (2017) ·Zbl 1391.35196号 ·doi:10.1051平方米/2016073
[7] H.布雷齐斯,最大单调算子和半群收缩算子荷兰阿姆斯特丹-朗顿North-Holland出版公司;美国爱思唯尔出版公司,纽约,1973年·Zbl 0252.47055号
[8] H.布雷齐斯,分析功能内尔《数学收藏》,《Appliques pour la Maêtrise》,巴黎马森,1983年·Zbl 0511.46001号
[9] G.Caginalp,自由边界相场模型分析,Arch。理性力学。分析。,92, 205-245 (1986) ·Zbl 0608.35080号 ·doi:10.1007/BF00254827
[10] W.Chen,C.Wang,X.Wang和S.M.Wise,具有对数势的Cahn-Hilliard方程的保正能量稳定数值格式,J.计算。物理学。X(X), 3 (2019).
[11] L.Cherfils;H.Fakih;A.Miranville,图像修复中Bertozzi Esedoglu-Gilette Cahn Hilliard方程的有限维吸引子,逆问题。成像,9105-125(2015)·Zbl 1515.35061号 ·doi:10.3934/ipi.2015.9.105
[12] L.Cherfils;H.Fakih;A.Miranville,关于带对数非线性项的Bertozzi-Esedoglu-Gillette-Cahn-Hilliard方程,SIAM J.Imaging Sci。,8, 1123-1140 (2015) ·Zbl 1381.35082号 ·数字对象标识代码:10.1137/140985627
[13] L.切尔菲尔斯;S.Gatti;A.Miranville,更正:“具有动态边界条件和奇异势的Caginalp相场系统整体解的存在性”,J.Math。分析。申请。,348, 1029-1030 (2008) ·Zbl 1160.35433号 ·doi:10.1016/j.jmaa.2008.07.058
[14] L.Cherfils,S.Gatti和A.Miranville,具有奇异势和动态边界条件的Caginalp系统的有限维吸引子,牛市。Transilv公司。布拉什大学。, 2 (2009), 25-34. ·Zbl 1212.35011号
[15] L.Cherfils;S.加蒂;A.Miranville,具有奇异势和动态边界条件的Caginalp系统的长时间行为,Commun。纯应用程序。分析。,11, 2261-2290 (2012) ·Zbl 1267.35032号 ·doi:10.3934/cpaa.2012.11.2261
[16] L.Cherfils;A.Miranville,关于奇异势Caginalp系统渐近行为的一些结果,高等数学。科学。申请。,17, 107-129 (2007) ·Zbl 1145.35042号
[17] L.Cherfils;A.Miranville,关于具有动态边界条件和奇异势的Caginalp系统,Appl。数学。,54, 89-115 (2009) ·Zbl 1212.35012号 ·doi:10.1007/s10492-009-0008-6
[18] L.Cherfils;A.米兰维尔;S.Zelik,具有对数势的Cahn-Hilliard方程,米兰数学杂志。,79561-596(2011年)·Zbl 1250.35129号 ·doi:10.1007/s00032-011-0165-4
[19] L.Cherfils;M.Petcu,带非渗透壁的Cahn-Hilliard方程的数值分析,Numer。数学。,128, 517-549 (2014) ·Zbl 1325.65133号 ·doi:10.1007/s00211-014-0618-0
[20] L.Cherfils;M.Petcu;M.Pierre,具有动态边界条件的Cahn-Hilliard方程的数值分析,离散Contin。动态。系统。,1511-1533年(2010年)·Zbl 1204.65113号 ·doi:10.3934/dcds.2010.27.1511
[21] M.I.M.科佩蒂;C.M.Elliott,具有对数自由能的Cahn-Hilliard方程的数值分析,Numer。数学。,63, 39-65 (1992) ·Zbl 0762.65074号 ·doi:10.1007/BF01385847
[22] I.Ekeland和R.Temam,凸分析与变分问题《应用数学经典》,28,工业和应用数学学会(SIAM),宾夕法尼亚州费城,1999年·Zbl 0939.49002号
[23] H.P.费舍尔;P.Maass;W.Dieterich,旋节分解中的新型表面模式,物理学。修订稿。,79, 893-896 (1997) ·doi:10.1103/PhysRevLett.79.893
[24] H.P.Fischer;P.Maass;W.Dieterich,薄流中旋节分解模式的发散时间和长度尺度,Europhys。莱特。,42, 49-54 (1998)
[25] H.P.Fischer;J.Reinhard;W.迪特里希;J.-F.古耶;P.Maass;A.马霍弗;D.Reinel,时间依赖密度泛函理论和与壁接触的晶格气体系统动力学,J.Chem。物理。,108, 3028-3037 (1998) ·doi:10.1063/1.475690
[26] 富考;S.吉川;S.Wada,一维动态边界条件下Cahn-Hilliard方程的保结构有限差分格式,Commun。纯应用程序。分析。,16, 1915-1938 (2017) ·Zbl 1364.65222号 ·doi:10.3934/cpaa.2017093
[27] 格拉塞利(M.Grasselli);A.米兰维尔;G.Schimperna,具有耦合动态边界条件和奇异势的Caginalp相场系统,离散Contin。动态。系统。,67-98年8月28日(2010年)·Zbl 1194.35074号 ·doi:10.3934/dcds.2010.28.67
[28] 格拉塞利(M.Grasselli);H.Petzeltová;G.Schimperna,奇异势Caginalp系统解的长时间行为,Z.Ana。安文德。,25, 51-72 (2006) ·Zbl 1128.35021号 ·doi:10.4171/ZAA/1277
[29] F.Hecht,freefem++的新发展,J.Numer。数学。,20, 251-265 (2012) ·Zbl 1266.68090号 ·doi:10.1515/jnum-2012-0013
[30] H.以色列;A.米兰维尔;M.Petcu,具有动态边界条件的Cahn-Hilliard型方程的数值分析,Ric。材料,64,25-50(2015)·Zbl 1320.65142号 ·doi:10.1007/s11587-014-0187-7
[31] B.科瓦奇;C.Lubich,动态边界条件抛物线问题的数值分析,IMA J.Numer。分析。,37, 1-39 (2017) ·Zbl 1433.65216号 ·doi:10.1093/imanum/drw015
[32] J.-L.狮子,非莱茵艾利斯问题解决方案Quelques Méthodes de Résolution des Problémes aux Limites,Dunod;1969年,巴黎,Gauthier-Villars·Zbl 0189.40603号
[33] A.Miranville,《Cahn-Hilliard方程及其变体》,AIMS数学,2479-544(2017)·Zbl 1425.35086号 ·doi:10.3934/数学.2017.2.479
[34] A.米兰维尔;S.Zelik,具有奇异势和动态边界条件的Cahn-Hilliard方程,离散Contin。动态。系统。,28, 275-310 (2010) ·Zbl 1203.35046号 ·doi:10.3934/dcds.2010.28.275
[35] D.Mugnolo;S.Romanelli,Dirichlet形式的一般Wentzell边界条件,解析半群,余弦算子函数,Electron。J.微分方程,2006,1-20(2006)·Zbl 1117.47036号
[36] F.Nabet,具有动态边界条件的Cahn Hilliard方程有限体积格式的收敛性,IMA J.Numer。分析。,36, 1898-1942 (2016) ·Zbl 1433.80007 ·doi:10.1093/imanum/drv057
[37] 皮埃尔先生;M.Pierre,Allen-Cahn-Gurtin方程通过多值算子的全局存在性,离散Contin。动态。系统。,33, 5347-5377 (2013) ·Zbl 1277.35227号 ·doi:10.3934/dcds.2013.33.5347
[38] J.Simon,空间中的紧集(L^p(0,T;B)),Ann.Mat.Pura应用。(4), 146 (1987), 65-96. ·Zbl 0629.46031号
[39] R.Temam,Navier-Stokes方程。理论与数值分析,《数学及其应用研究》,2,北荷兰出版公司,阿姆斯特丹,1984年·Zbl 0568.35002号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。