×

用移动网格混合有限元法模拟薄膜流动。 (英语) Zbl 1427.76153号

总结:我们提出了一种高效的混合有限元方法,用移动网格细化法求解四阶薄膜流动方程。移动网格策略基于由开发的调和映射R.李等[J.Comput.Phys.170,No.2,562-588(2001;Zbl 0986.65090号); 同上,177,第2号,365-393(2002年;Zbl 0998.65105号)]. 为了获得高质量的网格,我们采用了自适应监控功能,并基于扩散机制对其进行平滑处理。进行了各种数值测试,以证明该方法的准确性和有效性。与使用固定网格的数值模拟相比,移动网格细化可以准确地解决过冲和下冲结构,并降低计算成本。

MSC公司:

76M10个 有限元方法在流体力学问题中的应用
65M60毫米 涉及偏微分方程初值和初边值问题的有限元、Rayleigh-Ritz和Galerkin方法
76A20个 液体薄膜
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] Huppert,H.E.,《斜坡下粘性流的流动和不稳定性》,《自然》,300,5891,427-429(1982)
[2] 西尔维,N。;Dussan V,E.B.,液体对倾斜固体表面的重新润湿,物理学。流体,28,1,5-7(1985)
[3] Cazabat,A.M。;赫斯洛特,F。;Troian,S.M。;Carles,P.,温度梯度驱动的薄膜扩散的指状不稳定性,《自然》,346,6287,824(1990)
[4] Bernis,F。;Friedman,A.,高阶非线性退化抛物方程,J.Differ。Equ.、。,83179-206(1990年)·Zbl 0702.35143号
[5] 贝尔托齐,A.L。;Münch,A。;Shearer,M.,《薄膜流动中的欠压冲击》,Phys。D、 134、4、431-464(1999)·Zbl 1076.76509号
[6] 贝尔托齐,A.L。;Shearer,M.,薄膜方程中欠压行波的存在,SIAM J.Math。分析。,32, 1, 194-213 (2000) ·Zbl 0959.35122号
[7] Kondic,L。;Diez,J.,薄膜沿斜坡流动中的图案形成:恒定通量配置,Phys。流体,13,11,3168-3184(2001)·Zbl 1184.76297号
[8] Kondic,L.,薄膜重力驱动流动的不稳定性,SIAM Rev.,45,1,95-115(2003)·兹比尔1127.76313
[9] 利维,R。;希勒,M.,《驱动薄膜流动的动力学和成核》,物理学。D、 209、1、145-163(2005)·兹比尔1142.76335
[10] 贝尔托齐,A.L。;Bowen,M.,薄膜动力学:理论与应用,科学计算与应用中的现代方法,31-79(2002),Springer·Zbl 1198.76001号
[11] Ha,Y。;Kim,Y.-J。;Myers,T.G.,关于驱动薄膜方程的数值解,J.Compute。物理。,227, 15, 7246-7263 (2008) ·Zbl 1140.76027号
[12] 李毅。;Lee,H.G。;Yoon,D。;黄,W。;Shin,S。;Ha,Y。;Kim,J.,薄膜方程指进现象的数值研究,国际J。数值。《液体方法》,67,11,1358-1372(2011)·Zbl 1426.76047号
[13] Sun,P。;罗素·R·D。;Xu,J.,一种新的自适应局部网格细化算法及其在四阶薄膜流动问题中的应用,J.Compute。物理。,224, 2, 1021-1048 (2007) ·Zbl 1123.76343号
[14] 李毅。;Jeong,D。;Kim,J.,薄膜流动模拟的自适应网格细化,麦加尼卡,49,1,239-252(2014)·Zbl 1347.76034号
[15] Alharbi,A。;Naire,S.,具有表面张力的薄膜流动方程的自适应移动网格方法,J.Comput。申请。数学。,319, 365-384 (2017) ·Zbl 1457.76109号
[16] 黄,W。;Russell,R.D.,《自适应网格运动——MMPDE方法及其应用》,J.Compute。申请。数学。,128, 1, 383-398 (2001) ·Zbl 0971.65088号
[17] Alharbi,A.R.,使用自适应移动网格方法数值求解薄膜流动方程(2016),基尔大学(博士论文)
[18] 李,R。;Tang,T。;Zhang,P.,基于调和映射的多维移动网格方法,J.Compute。物理。,170, 2, 562-588 (2001) ·Zbl 0986.65090号
[19] 李,R。;Tang,T。;Zhang,P.,二维和三维奇异问题的移动网格有限元算法,J.Compute。物理。,177, 2, 365-393 (2002) ·兹比尔0998.65105
[20] Tang,H。;Tang,T.,一维和二维双曲守恒律的自适应网格方法,SIAM J.Numer。分析。,41887-515(2003年)·Zbl 1052.65079号
[21] 乔,Z。,Gierer-Meinhardt系统动力学行为和不稳定性的数值研究,Commun。计算。物理。,406-426年3月2日(2008年)·Zbl 1199.65338号
[22] Wang,H。;李,R。;Tang,T.,用r-自适应有限元方法有效计算枝晶生长,J.Compute。物理。,227, 12, 5984-6000 (2008) ·Zbl 1148.65052号
[23] Bao,G。;胡,G。;Liu,D.,采用自适应网格重分布技术的有限元方法对Kohn-Sham方程进行数值求解,J.Sci。计算。,55, 2, 372-391 (2013) ·Zbl 1272.82003年
[24] Tang,T.,计算流体动力学的移动网格方法,Contemp。数学。,383, 141-174 (2005) ·兹比尔1179.76065
[25] 格林斯潘,H.P.,《关于润湿表面的粘性小液滴的运动》,J.流体力学。,84, 01, 125-143 (1978) ·Zbl 0373.76040号
[26] Finn,R.,《毛细管表面界面》,通知AMS,46,7,770-781(1999)·Zbl 1194.76049号
[27] 马塔,M.R。;Bertozzi,A.L.,《二维颗粒层薄膜流动的数值格式》,J.Compute。物理。,230, 16, 6334-6353 (2011) ·Zbl 1419.76482号
[28] Witelski,T.P。;Bowen,M.,高阶非线性扩散方程的ADI格式,Appl。数字数学。,45, 2-3, 331-351 (2003) ·Zbl 1061.76051号
[29] Hundsdorfer,W。;Verwer,J.G.,时间相关对流扩散反应方程的数值解,3(2013),Springer Science&Business Media
[30] 阿舍尔,U.M。;Ruuth,S.J。;Spiteri,R.J.,含时偏微分方程的隐显Runge-Kutta方法,Appl。数字数学。,25, 2-3, 151-167 (1997) ·Zbl 0896.65061号
[31] Boyanova,P。;Do-Quang,M。;Neytcheva,M.,《大规模二进制Cahn-Hilliard模型的高效预条件》,计算。方法应用。数学。,12, 1, 1-22 (2012) ·Zbl 1284.65047号
[32] 博世,J。;斯托尔,M。;Benner,P.,使用隐式时间离散化和有限元快速求解Cahn-Hilliard变分不等式,J.Compute。物理。,262,38-57(2014)·Zbl 1349.65438号
[33] 博世,J。;Stoll,M.,向量值Cahn Hilliard方程的预处理,SIAM J.Sci。计算。,37、5、S216-S243(2015)·Zbl 1325.65043号
[34] 贝克特,G。;麦肯齐,J.A。;Ramage,A。;Sloan,D.M.,使用移动网格方法的二维非定常偏微分方程的计算解,J.Comput。物理。,182, 2, 478-495 (2002) ·Zbl 1016.65062号
[35] Zegeling,P.A。;De Boer,W.D。;Tang,H.Z.,二维欧拉方程的稳健高效自适应移动网格解,Contemp。数学。,383, 419-430 (2005) ·Zbl 1096.76035号
[36] van Dam,A。;Zegeling,P.A.,移动网格方法中流动现象的平衡监测,Commun。计算。物理。,7, 1, 138-170 (2010) ·Zbl 1364.76123号
[37] Huang,W.,移动网格偏微分方程的公式化和求解的实用方面,J.Compute。物理。,171, 2, 753-775 (2001) ·Zbl 0990.65107号
[38] Verwer,J.G。;布洛姆,J.G。;Furzeland,R.M。;Zegeling,P.A.,基于线方法的一维偏微分方程的移动网格方法(1989),SIAM,费城·Zbl 0697.65086号
[39] 黄,W。;Russell,R.D.,《带空间平滑的移动网格偏微分方程分析》,SIAM J.Numer。分析。,34, 3, 1106-1126 (1997) ·Zbl 0874.65071号
[40] 贝克特,G。;麦肯齐,J.A。;Ramage,A。;Sloan,D.M.,《使用基于均匀分布的自适应方法求解一维偏微分方程》,J.Compute。物理。,167, 2, 372-392 (2001) ·Zbl 0985.65097号
[41] 胡,G。;Zegeling,P.A.,用移动有限元法模拟多孔介质中的指状现象,J.Compute。物理。,230, 8, 3249-3263 (2011) ·Zbl 1316.76048号
[42] 张,H。;Zegeling,P.A.,多孔介质中非平衡方程非单调解的移动网格有限差分法,Commun。计算。物理。,935-964年4月22日(2017年)·Zbl 1488.74135号
[43] 张,H。;Zegeling,P.A.,通过移动网格法对多孔介质中具有动态毛细压力的两相流进行数值研究,J.Compute。物理。,345510-527(2017)·Zbl 1378.76077号
[44] Zegeling,P.A.,《自适应移动网格方法的理论与应用》,(Tao,T.;Xu,J.,《自适应计算:理论与算法》(2007),科学出版社,279-332
[45] R.Li,W.Liu,Afepack:自适应有限元软件包,可从以下网址获得:http://dsec.pku.edu.cn/rli/\(软件_e\).php;R.Li,W.Liu,Afepack:自适应有限元软件包,可从以下网址获得:http://dsec.pku.edu.cn/rli/\(软件\).php
[46] B.Niceno,Easymesh:基于Delaunay三角剖分的免费二维质量网格生成器,可在以下网站获得:http://www-dinma.univ.trieste.it/nirftc/research/easymesh/Default.htm; B.Niceno,Easymesh:基于Delaunay三角剖分的免费二维质量网格生成器,可在以下网站获得:http://www-dinma.univ.trieste.it/nirftc/research/easymesh/Default.htm
[47] 班杰斯,W。;哈特曼,R。;Kanschat,G.交易。II-通用面向对象有限元库,ACM Trans。数学。柔和。(TOMS),33,4,24(2007)·Zbl 1365.65248号
[48] Heroux,医学硕士。;Bartlett,R.A。;豪尔,V.E。;霍克斯特拉,R.J。;胡建杰。;科尔达·T·G。;Lehoucq,R.B。;Long,K.R。;Pawlowski,R.P。;Phipps,E.T.,《trilinos项目概述》,ACM Trans。数学。柔和。(TOMS),31,3,397-423(2005)·Zbl 1136.65354号
[49] The Mathworks Inc.,MATLAB版本8.3.0.532(R2014a)(2014),The Mathworks,Inc.:The Mathworks,Inc.,马萨诸塞州纳蒂克
[50] Getling,A.V.,Rayleigh-Bénard对流:结构和动力学,11(1998),《世界科学》·Zbl 0910.76001号
[51] 贝尔托齐,A.L。;Münch,A。;希勒,M。;Zumbrun,K.,《压缩和欠压缩薄膜行波的稳定性》,《欧洲应用杂志》。数学。,12, 03, 253-291 (2001) ·Zbl 1028.76012号
[52] Kang,M。;Kim,C.H。;Ha,Y.,WENO型方案与薄膜模拟的数值比较,J.Korea Soc.Ind.Appl。数学。,16, 3, 193-204 (2012) ·Zbl 1316.76011号
[53] Troian,S.M。;Herbolzheimer,E。;Safran,S.A。;Joanny,J.F.,驱动扩散膜的指状不稳定性,Europhys。莱特。,10, 1, 25 (1989)
[54] 贝尔托齐,A.L。;Brenner,M.P.,驱动接触线的线性稳定性和瞬态增长,Phys。流体,9,3,530-539(1997)·Zbl 1185.76627号
[55] 苏尔·J。;贝尔托齐,A.L。;Behringer,R.P.,驱动薄膜流中的反向欠压激波结构,物理。修订稿。,90, 12, 126105 (2003)
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。