×

向列相液晶和胆甾相液晶的Oseen-Frank模型的增强拉格朗日预处理器。 (英语) Zbl 1482.76008号

作者摘要:我们提出了一个鲁棒有效的增广拉格朗日型预处理器,用于解决向列相和胆甾相液晶中Oseen-Frank模型的线性化问题。通过应用增广拉格朗日方法,可以通过拉格朗日乘法器的加权质量矩阵更好地近似指向矢块的Schur补,代价是使增广指向矢块更难求解。为了求解增广的方向图块,我们开发了一种鲁棒的多重网格算法,其中包括一种附加的Schwarz松弛,它捕获了半定项核的逐点版本。此外,我们证明了增广拉格朗日项改进了单位长度约束的离散执行。数值实验验证了该算法的有效性及其对问题相关参数(弗兰克常数和胆甾醇间距)和网格尺寸的鲁棒性。

MSC公司:

76甲15 液晶
65F08个 迭代方法的前置条件
65N30型 偏微分方程边值问题的有限元、Rayleigh-Riz和Galerkin方法
65号55 多重网格方法;含偏微分方程边值问题的区域分解
82天30分 随机介质、无序材料(包括液晶和自旋玻璃)的统计力学
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] 阿德勒,JH;阿瑟顿,TJ;Benson,T。;德意志银行艾默生;麦克拉克伦(MacLachlan),SP,《电和柔性电效应下液晶平衡的能量最小化》,SIAM J.Sci。计算。,37、5、S157-S176(2015)·Zbl 1446.76065号
[2] 阿德勒,JH;阿瑟顿,TJ;德意志银行艾默生;Maclachlan,SP,向列相液晶Frank-Oseen模型的能量最小化有限元方法,SIAM J.Numer。分析。,53, 5, 2226-2254 (2015) ·Zbl 1327.76015号
[3] 阿德勒,JH;阿瑟顿,TJ;艾默生,德意志银行;Maclachlan,SP,液晶平衡的约束优化,SIAM J.Sci。计算。,38, 1, 50-76 (2016) ·Zbl 1427.76028号
[4] Balay,S.、Abhyankar,S.,Adams,M.F.、Brown,J.、Brune,P.、Buschelman,K.、Dalcin,L.、Eijkhout,V.、Gropp,W.D.、Kaushik,D.、Knepley,M.、McInnes,L.C.、Rupp,K.,Smith,B.F.、Zhang,H.:PETSc用户手册。技术代表ANL-95/11-3.9版,阿贡国家实验室(2018)
[5] 鲍尔,吉咪,数学与液晶,分子晶体。液体结晶。,647,1,1-27(2017)·Zbl 1393.35157号
[6] 贝德福德,S.J.:变分法及其在液晶中的应用。牛津大学博士论文(2014)
[7] Beik,FPA;Benzi,M.,液晶控制器建模产生的鞍点系统的块预处理器,Calcolo,55,29,1-12(2018)·Zbl 1416.65079号
[8] Beik,FPA;Benzi,M.,双鞍点系统的迭代方法,SIAM J.矩阵分析。申请。,39, 902-921 (2018) ·Zbl 1391.65062号
[9] Benzi,M。;Beik,FPA;比尼,D。;Di Benedetto,F。;蒂尔蒂什尼科夫,E。;Van Barel,M.,双鞍点系统的Uzawa型和增广拉格朗日方法,数值线性代数中的结构矩阵(2019),Cham:Springer,Cham
[10] Benzi,M。;Golub,生长激素;Liesen,J.,鞍点问题的数值解,数值学报。,14, 1-137 (2005) ·Zbl 1115.65034号
[11] Benzi,M。;Olshanskii,MA,基于拉格朗日的Oseen问题增强方法,SIAM J.Sci。计算。,28, 6, 2095-2113 (2006) ·Zbl 1126.76028号
[12] Chandrasekhar,S.,Liq.克里斯特。(1992),剑桥:剑桥大学出版社,剑桥·Zbl 0824.00003号
[13] Ciarlet,PG,《椭圆型问题的有限元》(1978),阿姆斯特丹,纽约,牛津:北荷兰,阿姆斯特朗,纽约,剑桥·兹伯利0383.6058
[14] Clément,P.:使用局部正则化的有限元函数逼近。Française Automat版本。Informat公司。Recherche Opérationnelle Sér。胭脂分析。编号。9(R-2),77-84(1975)·Zbl 0368.65008号
[15] de Gennes,PG,《液晶物理学》(1974),牛津:牛津大学出版社,牛津·Zbl 0295.76005号
[16] 埃尔曼,HC;西尔维斯特,D。;Wathen,AJ,《有限元和快速迭代求解器:在不可压缩流体动力学中的应用》(2014),牛津:牛津大学出版社,牛津·兹比尔1304.76002
[17] Emerson,D.B.:液晶配置的高级离散化和多重网格方法。塔夫茨大学博士论文(2015)
[18] 德意志银行艾默生;法雷尔,体育;阿德勒,JH;麦克拉克伦,SP;Atherton,TJ,用放气算法计算椭圆通道中胆甾液晶的平衡态,Liq.Cryst。,45, 3, 341-350 (2018)
[19] 埃里克森,JL,《可变取向度液晶》,Arch。定额。机械。分析。,113, 2, 97-120 (1991) ·Zbl 0729.76008号
[20] Farrell,P.E.,Knepley,M.G.,Wechsung,F.,Mitchell,L.:Pcpatch:多重网格松弛方法拓扑构造软件。arXiv预印arXiv:1912.08516(2019)
[21] 法雷尔,体育;米切尔,L。;Wechsung,F.,高雷诺数下三维定常不可压缩Navier-Stokes方程的增广拉格朗日预处理,SIAM J.Sci。计算。,41,A3073-A3096(2019)·Zbl 1448.65261号
[22] Firedrake-Zenodo:用于“向列相液晶和胆甾相液晶的Oseen-Frank模型的增强拉格朗日预编码器”的软件(2020年)。doi:10.5281/zenodo.4249051
[23] Fortin,M。;Glowinski,R.,《增广拉格朗日方法:在边值问题数值解中的应用,数学及其应用研究》(1983年),阿姆斯特丹:爱思唯尔科学有限公司·Zbl 0525.65045号
[24] 弗兰克,FC,液晶,法拉第讨论。,25, 19-28 (1958)
[25] Friedel,V.,《马蒂岛的海洋生物》,Ann.Phys。,18, 273-474 (1922)
[26] Geuzaine,C。;Remacle,JF,Gmsh:一个内置预处理和后处理设施的三维有限元网格生成器,Int.J.Numer。方法工程,79,11,1309-1331(2009)·Zbl 1176.74181号
[27] Girault,V。;Raviart,PA,《Navier-Stokes方程的有限元方法:理论和算法》(2011),Cham:Springer,Cham·Zbl 0585.65077号
[28] Glowinski,R.,Le Tallec,P.:解变分问题的增广拉格朗日方法,第3章。摘自:《应用数学研究》,第45-121页。SIAM(1989)·Zbl 0698.73001号
[29] 格洛文斯基,R。;林,P。;Pan,XB,液晶模型的运算符拆分方法,计算。物理学。社区。,152, 3, 242-252 (2003)
[30] 何,X。;武伊克,C。;Klaij,CM,将增广拉格朗日预条件子与简单Schur补码近似相结合,SIAM J.Sci。计算。,40、3、A1362-A1385(2018)·Zbl 1392.65053号
[31] Heister,T。;Rapin,G.,《使用Grad-Div稳定性解决Oseen问题的高效增广拉格朗日型预处理》,国际期刊Numer。《液体方法》,71,1,118-134(2012)·Zbl 1430.76115号
[32] 胡,Q。;泰,X。;Winther,R.,调和映射计算的鞍点方法,SIAM J.Numer。分析。,47, 2, 1500-1523 (2009) ·Zbl 1219.49030号
[33] 约翰·V。;Linke,A。;Merdon,C。;Neilan,M。;Rebholz,L.,《关于不可压缩流动混合有限元方法中的发散约束》,SIAM Rev.,59,3,492-544(2017)·Zbl 1426.76275号
[34] Knabner,P。;Angermann,L.,Numerik partieller Differentialgleichungen(2000),柏林:施普林格出版社,柏林·Zbl 0953.65075号
[35] Lee,Y。;吴杰。;徐,J。;Zikatanov,L.,几乎奇异系统的鲁棒子空间校正方法,数学。模型方法应用。科学。,17, 11, 1937-1963 (2007) ·Zbl 1151.65096号
[36] Lin,FH,向列相液晶中缺陷的非线性理论;相变和流动现象,Commun。纯应用程序。数学。,42, 6, 789-814 (1989) ·兹比尔0703.35173
[37] 林,P。;Richter,T.,《高维液晶分子取向的自适应同伦多重网格法》,J.Compute。物理。,225, 2, 2069-2082 (2007) ·Zbl 1147.82361号
[38] 林,P。;泰,X。;Fitzgibbon,W。;库兹涅佐夫,Y。;Neitaanmaki,P。;Pironneau,O.,液晶模型微观结构的增强拉格朗日方法,第7章,科学技术建模、仿真和优化,123-137(2014),多德雷赫特:施普林格,多德雷赫特·Zbl 1319.82023号
[39] Nocedal,J。;Wright,SJ,《数值优化》(1999),柏林:施普林格出版社,柏林·Zbl 0930.65067号
[40] Olshanskii,MA,定常不可压缩流Navier-Stokes方程的低阶Galerkin有限元方法:稳定问题和迭代方法,Comput。方法应用。机械。工程,191,47,5515-5536(2002)·Zbl 1083.76553号
[41] Olshanskii,马萨诸塞州;润滑油,G。;Heister,T。;Löwe,J.,不可压缩Navier-Stokes方程的梯度-直径稳定和亚网格压力模型,计算。方法应用。机械。工程,198,49,3975-3988(2009)·Zbl 1231.76161号
[42] Oseen,CW,液晶理论,Trans。法拉第社,29,140,883-899(1933)·Zbl 0008.04203号
[43] 波利亚克,弗吉尼亚州;Tret'yakov,NV,条件极值问题的惩罚估计方法,苏联计算。数学。数学。物理。,13, 1, 42-58 (1974) ·Zbl 0273.90055号
[44] Ramage,A。;Gartland,E.,《用于液晶指向器建模的预处理零空间方法》,SIAM J.Sci。计算。,35、1、B226-B247(2013)·Zbl 1276.82061号
[45] Rathgeber,F。;哈姆,DA;米切尔,L。;Lange先生。;卢波里尼,F。;ATT麦克雷;GT Bercea;马卡尔,GR;Kelly,PHJ,Firedrake:通过编写抽象来实现有限元方法的自动化,ACM Trans。数学。软质。,43, 3, 1-27 (2017) ·Zbl 1396.65144号
[46] Reinitzer,F.,Beiträge zur Kenntnis des Cholesterins,Monatsh。化学。,9, 421-441 (1888)
[47] Saad,Y.,《一种灵活的内外预处理GMRES算法》,SIAM J.Sci。计算。,14, 2, 461-469 (1993) ·Zbl 0780.65022号
[48] 萨阿德,Y。;Schultz,M.,GMRES:求解非对称线性系统的广义最小残差算法,SIAM J.Sci。统计计算。,7, 3, 856-869 (1986) ·Zbl 0599.65018号
[49] Schöberl,J.,原始变量中参数相关问题的多重网格方法,Numer。数学。,84,197-119(1999)·Zbl 0957.74059号
[50] Schöberl,J.:参数相关问题的稳健多重网格方法。林茨约翰内斯·开普勒大学博士论文(1999)·Zbl 0957.74059号
[51] 西尔维斯特,D。;Wathen,AJ,稳定Stokes系统的快速迭代解。第二部分:使用通用块预条件,SIAM J.Numer。分析。,31, 1352-1367 (1994) ·Zbl 0810.76044号
[52] Stewart,IW,《液晶的静态和动态连续理论:数学导论》(2004),博卡拉顿:CRC出版社,博卡拉顿
[53] Vanka,SP,二维循环流的块隐式多重网格计算,计算。方法。申请。机械。工程师,59,1,29-48(1986)·Zbl 0604.76025号
[54] Wathen,AJ;Silvester,D.,稳定Stokes系统的快速迭代解。第一部分:使用简单的对角线预条件,SIAM J.Numer。分析。,30, 3, 630-649 (1991) ·Zbl 0776.76024号
[55] 夏,J.:ALpaper-numerics(2020)。doi:10.5281/zenodo.4257094
[56] Xu,J.,通过空间分解和子空间校正的迭代方法,SIAM Rev.,34,4,581-613(1992)·Zbl 0788.65037号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。