×

采用以单元为中心的有限体积法求解具有二阶精确梯度的非梯度四叉树上的泊松问题。 (英语) Zbl 1380.65399号

摘要:本文介绍了泊松问题在自适应笛卡尔四叉树网格上的二维单元中心有限体积离散化,该网格在解及其梯度方面均表现出二阶精度,并且不需要相邻单元之间的分级条件。在T型结配置中,只要相邻单元之间的分辨率不同,就会出现这种配置,使用标准中心差分梯度模板需要通过插值来构造重影值。为了在生成的数值梯度中正确恢复二阶精度,之前处理块结构网格和分级树的工作表明,需要二次插值,而不是线性插值;否则,梯度只表现出一阶收敛,这限制了流体流动等潜在应用。然而,在一般非分级四叉树的情况下,由于存在更复杂的T结几何结构,以前的方案无法或失去准确性,这对相邻单元的分辨率没有任何限制。因此,我们针对这种情况提出了新的二次插值构造,通过依赖对角方向的模板实现二阶收敛,并根据需要递归应用。该方法处理复杂的树拓扑和相邻单元之间的大分辨率跳跃,甚至沿域边界,并且支持Dirichlet和Neumann边界条件。数值实验证实了该方法在\(L^\infty\)范数下的总体二阶精度。

MSC公司:

65牛顿50 涉及偏微分方程的边值问题的网格生成、细化和自适应方法
65号08 含偏微分方程边值问题的有限体积法
35J05型 拉普拉斯算子、亥姆霍兹方程(约化波动方程)、泊松方程
PDF格式BibTeX公司 XML格式引用
全文: 内政部 链接

参考文献:

[1] 格罗斯曼,C。;Roos,H.-G。;Stynes,M.:偏微分方程的数值处理。(2007) ·Zbl 1180.65147号 ·doi:10.1007/978-3-540-71584-9
[2] J.Guermond。;Minev,P。;沈,J.:不可压缩流投影方法概述。计算。方法应用。机械。工程195,编号44-47,6011-6045(2006)·Zbl 1122.76072号 ·doi:10.1016/j.cma.2005.10.010
[3] 最小值C。;Gibou,F。;Ceniceros,H.D.:非梯度网格上变系数泊松方程的超收敛有限差分格式。J.计算。物理学。218,第1期,123-140(2006)·Zbl 1106.65093号 ·doi:10.1016/j.jcp.2006.01.046
[4] 科尔拉,P。;格雷夫斯,D.T。;约翰逊,J.N。;基恩,N.D。;Ligocki,T。;马丁·D·F。;Mccorquodale,P.W。;莫迪亚诺,D。;P.O.施瓦茨。;斯特恩伯格,T.D。;Van Straalen,B.:AMR应用程序设计文件的Chombo软件包。(2013)
[5] Minion,M.L.:研究局部精细网格上涡旋斑块演化的两种方法。(1994年)
[6] 马丁·D·F。;Cartwright,K.L.:使用自适应网格细化求解泊松方程。(1996年)
[7] 约翰森,H。;Colella,P.:不规则区域上泊松方程的笛卡尔网格嵌入边界法。J.计算。物理学。147,第1号,60-85(1998)·Zbl 0923.65079号 ·doi:10.1006/jcph.1998.5965
[8] 马丁·D·F。;科尔拉,P。;Graves,D.T.:三维不可压缩Navier-Stokes方程的以细胞为中心的自适应投影方法。J.计算。物理学。227,第3期,1863-1886(2008)·Zbl 1137.76040号 ·doi:10.1016/j.jcp.2007.09.032
[9] Pletzer,A。;Jamroz,B。;克罗克特,R。;Sides,S.:紧凑的以单元为中心的离散化模板,位于精细的块结构网格界面。J.计算。物理学。260, 25-36 (2014) ·Zbl 1349.65570号 ·doi:10.1016/j.jcp.2013.12.020
[10] Liu,Q.:在局部细化交错网格上的一种稳定而精确的投影方法。国际期刊数字。方法流体67,No.1,74-92(2011)·Zbl 1303.76107号 ·doi:10.1002/fld.2338
[11] Henshaw,W.:重叠网格上不可压缩Navier-Stokes方程的四阶精确方法。J.计算。物理学。113,第1期,第13-25页(1994年)·Zbl 0808.76059号 ·doi:10.1006/jcph.1994.1114
[12] 英语,R.E。;邱,L。;Yu,Y。;Fedkiw,R.:使用多个移动笛卡尔网格对不可压缩流进行自适应离散。J.计算。物理学。254, 107-154 (2013) ·Zbl 1349.76456号 ·doi:10.1016/j.jp.2013.07.032
[13] Popinet,S.:Gerris:复杂几何中不可压缩Euler方程的基于树的自适应求解器。J.计算。物理学。190,第2号,572-600(2003)·Zbl 1076.76002号 ·doi:10.1016/S0021-9991(03)00298-5
[14] Losasso,F。;Gibou,F。;Fedkiw,R.:用八叉树数据结构模拟水和烟雾。ACM变速器。图表。23,第3期,457-462(2004)
[15] Losasso,F。;Fedkiw,R。;Osher,S.:水平集方法和不可压缩流的空间自适应技术。计算。流体35,No.10,995-1010(2005)·Zbl 1177.76295号 ·doi:10.1016/j.compfluid.2005.01.006
[16] 芬克尔,R。;Bentley,J.L.:四叉树:用于复合键检索的数据结构。Acta通知。第4期,第1期,第1-9期(1974年)·Zbl 0278.68030号 ·doi:10.1007/BF00288933
[17] Meagher,D.:使用八叉树编码的几何建模。计算。图表。图像处理。19,第2期,129-147(1982)
[18] Minion,M.L.:局部细化网格的投影方法。J.计算。物理学。127,第1期,158-178(1996)·Zbl 0859.76047号 ·doi:10.1006/jcph.1996.0166
[19] 吉特·A。;米歇尔·泰拉德。;Gibou,F.:在任意几何体和自适应四/八叉树上求解不可压缩Navier-Stokes方程的稳定投影方法。J.计算。物理学。292, 215-238 (2015) ·Zbl 1349.76336号 ·doi:10.1016/j.jcp.2015.03.024
[20] 最小值C。;Gibou,F.:非梯度自适应网格上不可压缩Navier-Stokes方程的二阶精确投影方法。J.计算机。物理学。219,第2期,912-929(2006)·Zbl 1330.76096号 ·doi:10.1016/j.jp.2006.07.019
[21] 古普塔:一种用于从有限网格过渡到粗网格的有限元。国际期刊《数值》。方法工程12,No.1,35-45(1978)·兹比尔0369.73073 ·doi:10.1002/nme.1620120104
[22] Tabarraei,A。;Sukumar,N.:一致四叉树网格上的自适应计算。有限元。分析。目的。41,第7期,686-702(2005)
[23] Natarajan,S。;Ooi,E.T。;Man,H。;Song,C.:四叉树网格上的有限元计算:应变平滑和半分析公式。国际高级工程师科学杂志。应用。数学。7、3号、124-133(2015)·Zbl 1342.74168号 ·doi:10.1007/s12572-015-0134-1
[24] 巴拉德,M。;Colella,P.:泊松方程的四阶精确局部求精方法。J.计算。物理学。209,第1期,第1-18页(2005年)·Zbl 1073.65126号 ·doi:10.1016/j.jcp.2005.02.027
[25] Seibold,B.:泊松方程无网格有限差分方法中的最小正模板。计算。方法应用。机械。eng.198,No.3,592-601(2008)·Zbl 1228.65209号 ·doi:10.1016/j.cma.2008.09.001
[26] Olshanskii,医学硕士。;Terekhov,K.M。;Vassilevski,Y.V.:基于八叉树的不可压缩Navier-Stokes方程求解器,具有增强的稳定性和低耗散性。计算。流体84,231-246(2013)·Zbl 1290.76106号 ·doi:10.1016/j.compfluid.2013.04.027
[27] 最小值C。;Gibou,F.:非梯度自适应笛卡尔网格上的二阶精确水平集方法。J.计算。物理学。225,第1期,300-321(2007)·兹比尔1122.65077 ·doi:10.1016/j.jcp.2006.11.034
[28] 陈,H。;最小值C。;Gibou,F.:具有超线性收敛速度的自适应笛卡尔网格上Stefan问题的数值格式。J.计算。物理学。228,第16号,5803-5818(2009)·兹比尔1176.80059 ·doi:10.1016/j.jcp.2009.04.044
[29] 帕帕克,J。;Helgadottir,A。;Ratsch,C。;Gibou,F.:四叉树/八叉树自适应笛卡尔网格上Robin边界条件下扩散和Stefan型问题的水平集方法。J.计算。物理学。233, 241-261 (2013)
[30] 米歇尔·泰拉德。;Djodom,L.F。;维耶,J.-L。;Gibou,F.:一种求解不规则区域和自适应网格上线性弹性方程的二阶锐利数值方法-在形状优化中的应用。J.计算。物理学。233, 430-448 (2013) ·Zbl 1286.74109号 ·doi:10.1016/j.jcp.2012.09.002
[31] 米尔扎德,M。;Gibou,M.Theillard;Gibou,F.:使用非梯度自适应笛卡尔网格对不规则几何体上的非线性泊松-玻耳兹曼方程进行二阶离散。J.计算。物理学。230,第5期,2125-2140(2011)·Zbl 1390.82056号
[32] 阿尔卑斯省赫尔加多蒂尔。;Gibou,F.:非梯度自适应网格上具有Neumann或Robin边界条件的不规则区域上的泊松-玻尔兹曼解算器。J.comput。物理学。230,第10号,3830-3848(2011)·Zbl 1369.76033号
[33] 米尔扎德,M。;Gibou,F.:自适应笛卡尔网格上泊松-能斯特-普朗克方程的保守离散化。J.计算。物理学。274, 633-653 (2014) ·Zbl 1351.82082号 ·doi:10.1016/j.jcp.2014.06.039
[34] 马丁,D.F。;Colella,P.:不可压缩欧拉方程的以细胞为中心的自适应投影方法。J.计算。物理学。163,编号271-312(2000)·Zbl 0991.76052号 ·doi:10.1006/jcph.2000.6575
[35] Chesshire,G。;Henshaw,W.:用于解偏微分方程的复合重叠网格。J.计算。物理学。90, 1-64 (1990) ·Zbl 0709.65090号 ·doi:10.1016/0021-9991(90)90196-8
[36] 肖特利,G.H。;Weller,R.:拉普拉斯方程的数值解。J.应用。物理学。9,No.5,334(1938)·Zbl 0019.03801号 ·doi:10.1063/1.1710426
[37] Gibou,F。;Fedkiw,R。;Cheng,L.-T。;Kang,M.:不规则区域上泊松方程的二阶精确对称离散化。J.计算。物理学。176,第1号,205-227(2002)·Zbl 0996.65108号 ·doi:10.1006/jcph.2001.6977
[38] Ng,Y.T。;陈,H。;最小值C。;Gibou,F.:使用鬼影流体方法在具有Dirichlet边界条件的不规则区域上的泊松解算器指南。科学杂志。计算。41,第2期,300-320(2009)·Zbl 1203.65223号 ·doi:10.1007/s10915-009-9299-8
[39] Stanoyevitch,A.:使用Matlab介绍数值常微分方程和偏微分方程。(2011) ·Zbl 1071.65002号
[40] Guennebaud,G。;雅各布:其他,特征v3。(2010)
[41] Li,X.S.:超级lu概述:算法、实现和用户界面。ACM变速器。数学。柔软。第302-325号第31页(2005年)·Zbl 1136.65312号 ·doi:10.1145/1089014.1089017
[42] Ng,Y.T.先生。;最小值C。;Gibou,F.:单相流的高效流固耦合算法。J.计算机。物理学。228,第23号,8807-8829(2009)·Zbl 1245.76019号 ·doi:10.1016/j.jcp.2009.08.032
[43] Demmel,J.W.:应用数值线性代数。(1997) ·Zbl 0879.65017号
[44] 陈,H。;最小值C。;Gibou,F.:不规则区域和非梯度自适应笛卡尔网格上泊松和热方程的超收敛有限差分格式。科学杂志。计算。31,第1期,19-60(2007)·Zbl 1120.65113号 ·doi:10.1007/s10915-006-9122-8
[45] 布雷齐,F。;Lipnikov,K。;Simoncini,V.:多边形和多面体网格上的一系列模拟有限差分方法。数学。模型方法应用。科学。15,第10期,1533-1551(2005)·Zbl 1083.65099号 ·doi:10.1142/S0218205000832
[46] Da Veiga,L.B。;Manzini,G.:模拟有限差分方法的高阶公式。SIAM J.科学。计算机。31,第1期,732-760(2008)·Zbl 1185.65201号 ·doi:10.1137/080717894
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。