×

基于TT-M有限元方法的修正Cahn-Hilliard方程快速算法。 (英语) Zbl 1490.65213号

小结:我们考虑数值方法来求解含有强非线性的修正Cahn-Hilliard方程。为了克服非线性项计算耗时的问题,提出了一种基于时间双网格(TT-M)有限元格式的快速算法。TT-M有限元算法包括三个主要步骤:首先,在粗时间网格(tau_c)上求解非线性有限元格式。这里,空间离散采用有限元方法,时间离散采用隐式二阶(θ)格式(包含隐式Crank-Nicolson格式和二阶后向差分方法)。其次,利用拉格朗日插值在精细时间网格上得到插值结果。最后,在精细时间网格(τ)(τ<τ_c)上求解线性化有限元系统。文中详细给出了稳定性分析和先验误差估计。数值算例验证了该方案的有效性。将TT-M有限元方法与传统的Galerkin有限元方法进行了比较,结果表明TT-M方法可以节省计算时间。

MSC公司:

65M60毫米 涉及偏微分方程初值和初边值问题的有限元、Rayleigh-Ritz和Galerkin方法
PDF格式BibTeX公司 XML格式引用
全文: DOI程序

参考文献:

[1] A.C.Aristotelous、O.Karakashian和S.M.Wise,“修正Cahn-Hilliard方程的混合间断Galerkin凸分裂格式和有效的非线性多重网格解算器”,《离散Contin》。戴恩。系统。序列号。B18:9(2013),2211-2238·Zbl 1278.65149号
[2] M.Asle Zaeem、H.El Kadiri、M.F.Horstemeyer、M.Khafizov和Z.Utegulov,“内应力和中间相对相干沉淀粗化的影响:相场研究”,当前。申请。《物理》12:2(2012),570-580。
[3] B.Ayuso、B.García-Archilla和J.Novo,“Navier-Stokes方程的后处理混合有限元方法”,SIAM J.Numer。分析43:3(2005),1091-1111·Zbl 1094.76037号
[4] A.L.Bertozzi、S.Esedoíglu和A.Gillette,“使用Cahn-Hilliard方程绘制二值图像”,IEEE Trans。图像处理。16:1(2007),285-291·Zbl 1279.94008号
[5] J.W.Cahn,“非均匀系统的自由能。II:热力学基础”,J.Chem。Phys.30:5(1959),1121-1124。
[6] J.W.Cahn和J.E.Hilliard,“非均匀体系的自由能。I:界面自由能”,《化学杂志》。《物理学》28:2(1958),258-267·兹比尔1431.35066
[7] J.W.Cahn和J.E.Hilliard,“非均匀系统的自由能。III:双组分不可压缩流体中的成核”,J.Chem。《物理学》31:3(1959),688-699。
[8] D.Carolan、H.M.Chong、A.Ivankovic、A.J.Kinloch和A.C.Taylor,“含钴聚合物系统:数值研究”,计算。马特。《科学》第98卷(2015年),第24-33页。
[9] R.Choksi、M.A.Peletier和J.F.Williams,“关于二嵌段共聚物微相分离的相图:通过非局部Cahn-Hilliard泛函的方法”,SIAM J.Appl。数学69:6(2009),1712-1738·Zbl 1400.74089号
[10] C.M.Elliott和D.A.French,“二维Cahn-Hilliard方程的非协调有限元方法”,SIAM J.Numer。分析26:4(1989),884-903·Zbl 0686.65086号
[11] D.J.Eyre,“无条件梯度稳定时间推进Cahn-Hilliard方程”,第39-46页,微观结构演化的计算和数学模型(旧金山,1998年),由J.W.Bullard等人编辑,Mater。Res.Soc.交响乐。Proc.529,MRS,Warrendale,PA,1998年。
[12] G.-H.Gao、H.-W.Sun和Z.-Z.Sun,“基于某种超收敛性的一类时间分数次扩散方程的有限差分格式的稳定性和收敛性”,J.Compute。《物理学》第280卷(2015年),第510-528页·Zbl 1349.65295号
[13] F.Guillén-González和G.Tierra,“关于Cahn-Hilliard扩散界面模型的线性方案”,J.Compute。《物理学》234(2013),140-171·Zbl 1284.35025号
[14] F.Guillén-González和G.Tierra,“Allen-Cahn和Cahn-Hilliard模型的二阶格式和时间步长自适应性”,计算。数学。应用68:8(2014),821-846·Zbl 1362.65104号
[15] D.Han、A.Brylev、X.Yang和Z.Tan,“两相不可压缩流相场模型的二阶全离散能量稳定格式的数值分析”,《科学杂志》。计算70:3(2017),965-989·Zbl 1397.76070号
[16] Y.He和Y.Liu,“Cahn-Hilliard方程谱Galerkin方法的稳定性和收敛性”,数值。方法偏微分方程24:6(2008),1485-1500·Zbl 1157.65053号
[17] 何毅、刘毅和唐涛,“关于Cahn-Hilliard方程的大时间步长方法”,应用。数字。数学57:5-7(2007),616-628·Zbl 1118.65109号
[18] 何毅、刘毅和唐涛,“关于Cahn-Hilliard方程的大时间步长方法”,应用。数字。数学57:5-7(2007),616-628·Zbl 1118.65109号
[19] M.Heida,“关于Cahn-Hilliard-Navier-Stokes系统热力学一致边界条件的推导”,国际。《工程科学杂志》第62期(2013年),第126-156页·Zbl 1423.76094号
[20] Z.Hu、S.M.Wise、C.Wang和J.S.Lowengrub,“相场晶体方程的稳定有效的有限差分非线性多重网格格式”,J.Comput。《物理学》228:15(2009),5323-5339·Zbl 1171.82015年
[21] H.Jia、Y.Li、G.Feng和K.Li,“具有浓度相关迁移率和对数Flory-Huggins体势的Cahn-Hilliard方程的有效双网格方法”,应用。数学。计算387(2020)、124548、15·Zbl 1472.65121号
[22] L.Ju、J.Zhang和Q.Du,“模拟Cahn-Hilliard方程粗化动力学的快速准确算法”,计算。马特。《科学》第108卷(2015年),第272-282页。
[23] W.Layton和L.Tobiska,“Navier-Stokes方程的带回溯的两级方法”,SIAM J.Numer。分析。35:5 (1998), 2035-2054. ·Zbl 0913.76050号
[24] D.Li和Z.Qiao,“关于二维Cahn-Hilliard方程的二阶半隐式傅里叶谱方法”,科学杂志。计算70:1(2017),301-341·Zbl 1434.65206号
[25] Y.Li、J.Shin、Y.Choi和J.Kim,“使用相场模型从切片数据重建三维体积”,计算。视觉。《图像理解》137(2015),115-124。
[26] D.Li、Z.Qiao和T.Tang,“表征相位场方程半隐式Fourier谱方法的稳定尺寸”,SIAM J.Numer。分析54:3(2016),1653-1681·Zbl 1346.35162号
[27] Y.Li、J.-I.Choi和J.Kim,“复杂区域中具有不同边界条件的多组分Cahn-Hilliard系统”,J.Compute。《物理学》第323卷(2016年),第1-16页·兹比尔1415.65219
[28] Y.Li、J.-I.Choi和J.Kim,“带界面剖面校正项的相场流体建模和计算”,Commun。非线性科学。数字。模拟30:1-3(2016),84-100·Zbl 1461.76512号
[29] Q.Liu和D.Shi,“Cahn-Hilliard方程能量稳定的Ciarlet-Raviart型格式的双网格方法的超收敛分析”,Numer。《算法》85:2(2020),607-622·Zbl 1465.65101号
[30] Q.Liu、Y.Hou、Z.Wang和J.Zhao,“Cahn-Hilliard方程的两级方法”,《数学》。计算。模拟126(2016),89-103·Zbl 07313662号
[31] Y.Liu,Z.Yu,H.Li,F.Liu,J.Wang,“时间分数水波模型的时间双网格算法与有限元方法相结合”,《国际热质传输》120(2018),1132-1145。
[32] Y.Liu、Y.Du、H.Li、F.Liu和Y.Wang,“非线性分数阶拉索方程的一些二阶θ格式与有限元方法相结合”,数值。《算法》80:2(2019),533-555·兹比尔1433.65218
[33] M.Marion和J.Xu,“基于双网格有限元的新型非线性Galerkin方法的误差估计”,SIAM J.Numer。分析32:4(1995),1170-1184·兹比尔0853.65092
[34] Y.Shang,“稳态Navier-Stokes方程的二级子网格稳定Oseen迭代方法”,J.Compute。《物理学》第233卷(2013年),第210-226页·Zbl 1286.76035号
[35] Y.Wang、Y.Liu、H.Li和J.Wang,“非线性分数阶Cable方程的有限元法与二阶时间离散格式相结合”,《欧洲物理学》。J.Plus131(2016),第61条。
[36] D.Wang、Q.Du、J.Zhang和H.Jia,“Allen-Cahn方程的快速双网格算法”,Bull。马来人。数学。科学。Soc.43:3(2020),2417-2441·Zbl 1437.65147号
[37] G.N.Wells、E.Kuhl和K.Garikipati,“Cahn-Hilliard方程的间断Galerkin方法”,J.Compute。《物理学》218:2(2006),860-877·Zbl 1106.65086号
[38] S.M.Wise、J.S.Lowengrub、H.B.Friebes和V.Cristini,“三维多物种非线性肿瘤生长-I:模型和数值方法”,J.Theoret。《生物学》253:3(2008),524-543·Zbl 1398.92135号
[39] S.M.Wise、C.Wang和J.S.Lowengrub,“相场晶体方程的能量稳定收敛有限差分格式”,SIAM J.Numer。分析47:3(2009),2269-2288·Zbl 1201.35027号
[40] Xia Y.Xu和C.-W.Shu,“Cahn-Hilliard型方程的局部间断Galerkin方法”,J.Compute。《物理学》227:1(2007),472-491·Zbl 1131.65088号
[41] J.Xu,“线性和非线性偏微分方程的两种网格离散化技术”,SIAM J.Numer。分析33:5(1996),1759-1777·Zbl 0860.65119号
[42] X.Yang,“均聚物共混物相场模型的线性、一阶和二阶无条件能量稳定数值格式”,J.Compute。《物理学》第327卷(2016年),第294-316页·Zbl 1373.82106号
[43] B.Yin,Y.Liu,H.Li,and S.He,“基于TT-M有限元系统的空间分数阶Allen-Cahn方程光滑和非光滑解的快速算法”,J.Compute。《物理学》379(2019),351-372·Zbl 07581576号
[44] Z.Zhang和Z.Qiao,“Cahn-Hilliard方程的自适应时间步长策略”,Commun。计算。《物理学》11:4(2012),1261-1278·Zbl 1388.65072号
[45] S.Zhang和M.Wang,“Cahn-Hilliard方程的非协调有限元方法”,J.Compute。物理学。229:19 (2010), 7361-7372. ·兹比尔1197.65139
[46] J.Zhao、X.Yang、Y.Gong和Q.Wang,“液晶流体动力学Q张量模型的新型线性二阶无条件能量稳定格式”,计算。方法应用。机械。工程318(2017),803-825·兹比尔1439.76124
[47] S.Zhou和M.Y.Wang,“利用多相转变的广义Cahn-Hilliard模型进行多材料结构拓扑优化”,结构。多磁盘。优化33:2(2007),89-111·Zbl 1245.74077号
[48] J.Zhou、L.Chen、Y.Huang和W.Wang,“Cahn-Hilliard方程的有效双网格格式”,Commun。计算。《物理学》17:1(2015),127-145·Zbl 1373.76111号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。