×

标准有限元速度场的后处理,用于地下水流动的精确粒子跟踪。 (英语) Zbl 1439.65193号

小结:粒子追踪是一种计算上有利且快速的方案,用于确定地下水文中的旅行时间和轨迹。精确的粒子跟踪需要元素级质量守恒、一致的速度场。标准线性伽辽金有限元法(FEM)无法满足此条件。我们提出了一个投影,它将一个基于三角形和四面体计算的非协调的元素给定速度场映射到一个低阶Raviart-Thomas-Nédélec(mathcal{RTN}_0)\)空间,满足精确粒子跟踪的要求。投影的基础是最小化标准FEM解和符合模型的水力梯度之间单元质心处的水力梯度差异{RTN}_0\)速度场使元素质量守恒。使用mathcal中的一致速度场{RTN}_0\)在三角形和四面体上的空间中,我们提出了发散流和非发散流的半解析粒子跟踪方法。我们将结果与为相同单元定义的以单元为中心的有限体积法以及考虑水力各向异性的测试案例的分析结果进行了比较。基于标准有限元解投影的速度场和相关粒子轨迹与有限体积法的结果相当,但投影场在分段均匀水力传导率区域内更平滑。当\(\mathcal{RTN}_0因此,投影的标准FEM解更准确,以单元为中心的有限体积方法的计算成本要小得多。

MSC公司:

65N30型 含偏微分方程边值问题的有限元、Rayleigh-Ritz和Galerkin方法
65号08 含偏微分方程边值问题的有限体积法
76M10个 有限元方法在流体力学问题中的应用
76个M12 有限体积法在流体力学问题中的应用
76年 强制对流
76S05号 多孔介质中的流动;过滤;渗流
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] Aquanty,Inc.水文圈。描述完全集成的地下和地面流动及溶质传输的三维数值模型:用户手册。安大略省安大略州Aquanty公司技术代表(2015年)
[2] 亚利桑那州阿奇利;麦克斯韦,RM;Navarre-Sitchler,AK,使用流线模拟非均质含水层中的随机反应迁移:受CO2影响的饮用水含水层中金属的动态释放和迁移,Adv.water Resour。,52, 93-106 (2013)
[3] 巴里亚瓦蒂,C。;Carstensen,C.,具有后验错误控制的最低阶Raviart-Tomas MFEM的三种Matlab实现,Compute。方法应用。数学。,5, 4, 333-361 (2005) ·兹比尔1086.65107
[4] 贝尔·J。;Cheng,AHD,地下水流动和污染物运移建模。《多孔介质传输理论与应用》(2010)第23卷,海德堡:施普林格·Zbl 1195.76002号
[5] 贝林,A。;Salandin,P。;Rinaldo,A.,《非均质多孔地层中分散的模拟:统计学、一阶理论、计算收敛性、水资源》。第28、9、2211-2227号决议(1992年)
[6] Bennett,JP;哈斯劳尔,CP;Cirpka,OA,沉积各向异性对堆积冲刷池结构中溶质混合的影响,水资源。研究,53,4,2813-2832(2017)
[7] 布雷齐,F。;道格拉斯,J。;杜兰,R。;Fortin,M.,三变量二阶椭圆问题的混合有限元,Numer。数学。,51237-250(1987年)·Zbl 0631.65107号
[8] 布雷齐,F。;道格拉斯,J。;Marini,LD,二阶椭圆问题的两类混合有限元,数值。数学。,47, 2, 217-235 (1985) ·Zbl 0599.65072号
[9] 西尔普卡,OA;Chiogna,G。;罗尔,M。;Bellin,A.,三维非稳态各向异性非均质多孔介质中的横向混合,水资源。研究,51,1,241-260(2015)
[10] 西尔普卡,OA;弗里德,EO;Helmig,R.,矩形和流线型网格上反应输运的数值方法,高级水资源。,22, 7, 711-728 (1999)
[11] 绳索,C。;Kinzelbach,W.,线性、双线性和三线性有限元中的连续地下水速度场和路径线,水资源。Res.,28,11,2903-2911(1992)
[12] 绳索,C。;Kinzelbach,W.,评论R.Mosé、P.Siegel、P.Ackerer和G.Chavent,Water Resour所著的“混合有限元近似在地下水流动模型中的应用:奢侈还是必要?”。1905-1909年第32、6号决议(1996年)
[13] 起重机,MJ;Blunt,MJ,基于流线的溶质运移模拟,水资源。研究,35,10,3061-3078(1999)
[14] Davis,TA,非对称模式多面方法的列预排序策略,ACM Trans。数学。柔和。(TOMS),30,2,165-195(2004)·Zbl 1072.65036号
[15] Durlowsky,LJ,混合和控制体积有限元逼近达西速度和相关量的精度,水资源。研究,30,4965-973(1994)
[16] Edwards,MG,《基于流动网格的非结构化、控制体积分布、全张量有限体积格式》,计算。地质科学。,6,3343-452(2002年)·Zbl 1036.76034号
[17] Forsyth,P.,《NAPL地下水污染的控制体积有限元方法》,SIAM J.Sci。统计计算。,12, 5, 1029-1057 (1991) ·Zbl 0725.76087号
[18] 弗里斯,H。;爱德华兹,M。;Mykkeltveit,J.,非结构化以单元为中心的三角形网格上的对称正定含通量全张量有限体积格式,SIAM J.Sci。计算。,31, 2, 1192-1220 (2009) ·Zbl 1190.65163号
[19] Ginn,TR,《非线性反应和混合的随机对流传输:具有流管内弥散的多组分反应系统的有限流管系综公式》,J.Contam。水文。,47, 1-2, 1-28 (2001)
[20] Hoteit,H。;埃赫尔,J。;莫塞,R。;菲利普,B。;Ackerer,P.,应用于多孔介质流动问题的混合方法的数值可靠性,计算。地质科学。,6, 2, 161-194 (2002) ·Zbl 1079.76581号
[21] 休斯,TJ;恩格尔,G。;Mazzei,L。;Larson,MG,连续Galerkin方法是局部保守的,J.Compute。物理。,163, 2, 467-488 (2000) ·Zbl 0969.65104号
[22] 胡亚科恩,PS;GF Pinder,《地下水流的计算方法》(1983),纽约:学术出版社,纽约·Zbl 0577.76001号
[23] Juanes,R。;Matringe,SF,《二维非结构网格上高阶流线追踪的统一公式》,J.Sci。计算。,38, 1, 50-73 (2009) ·Zbl 1203.76145号
[24] CE Kees;法瑞特,MW;Dawson,CN,变饱和流的局部保守稳定有限元方法,计算。方法应用。机械。工程,197,51,4610-4625(2008)·Zbl 1194.76123号
[25] Kinzelbach,W.,Grundwasser Schadstown运输模型的数字方法。Schriftenreihe GWF Wasser-Abwasser第21卷。R(1992),慕尼黑:慕尼黑奥尔登堡
[26] Kinzelbach,W。;Ackerer,P.,《非污染物扩散和生态迁移模型化》,水文地质学,2197-205(1986)
[27] Larson,MG公司;Niklasson,AJ,基于非连续富集的连续Galerkin方法的保守通量,CALCOLO,41,2,65-76(2004)·Zbl 1168.65415号
[28] Loschko,M。;Wöhling,T。;鲁道夫,DL;Cirpka,OA,《在含水层尺度反应迁移的随机框架中解释非均质含水层的反应潜力下降》,《水资源》。研究,54,1,442-463(2018)
[29] Matringe,平方英尺;Juanes,R。;Tchelepi,HA,《在普通三角形和四边形网格上模拟多孔介质流动的稳健流线追踪》,J.Compute。物理。,2192992-1012(2006年)·Zbl 1330.76075号
[30] Nédélec,JC,(\mathbb{R}^3)中的混合有限元3,数字。数学。,35, 3, 315-341 (1980) ·Zbl 0419.65069号
[31] Nédélec,JC,(\mathbb{R}^3\)中一类新的混合有限元3,数字。数学。,50, 1, 57-81 (1986) ·Zbl 0625.65107号
[32] Ngo,AQT;Bastian,P。;Ippisch,O.,稳态地下水流动和溶质运移问题的数值解:基于间断Galerkin的方法与流线扩散方法的比较,计算。方法应用。机械。工程,294,331-358(2015)·Zbl 1423.76433号
[33] 比值,LH;惠勒,MF;Kvamsdal,T。;Larson,MG,非保守通量与非均匀介质中传输兼容性的后处理,计算。方法应用。机械。工程,315799-830(2017)·兹比尔1439.76158
[34] Pollock,DW,有限差分模型路径线的半分析计算,地下水,26,6,743-750(1988)
[35] 波维奇,TJ;道森,中国;法瑞特,MW;Kees,CE,变密度流动和溶质运移的有限元方法,计算。地质科学。,17, 3, 529-549 (2013) ·Zbl 1392.76031号
[36] 普提,M。;Cordes,C.,四面体上扩散算子的有限元近似,SIAM J.Sci。计算。,19, 4, 1154-1168 (1998) ·Zbl 0915.65111号
[37] 普提,M。;Sartoretto,F.,《线性Galerkin与混合有限元二维流场》,《国际数值杂志》。方法流体,60,9,1011-1031(2009)·Zbl 1187.76692号
[38] Raviart,P.A.,Thomas,J.M.:二阶椭圆问题的混合有限元方法。收录于:Galligani,I.,Magenes,E.(编辑)《有限元方法的数学方面》,第292-315页。柏林施普林格(1977)·兹比尔0362.65089
[39] 罗杰斯,M。;柯比,R。;Logg,A.,《H(div)和H(curl)协调有限元的高效组装》,SIAM J.Sci。计算。,31, 6, 4130-4151 (2009) ·Zbl 1206.65248号
[40] Schiavazzi,D.,《冗余多分辨率不确定性传播》。帕多瓦大学博士论文(2013)
[41] Scudeler,C。;普提,M。;Paniconi,C.,运输模拟中伽辽金速度场的质量守恒重建,水资源高级研究所。,94, 470-485 (2016)
[42] Shewchuk,J.R.:http://www.cs.cmu.edu/地震/三角形.html。访问时间:2018年10月30日(2005)
[43] Sun,S。;惠勒,MF,《与交通兼容的速度数据投影》,计算。方法应用。机械。工程,195,7,653-673(2006)·Zbl 1091.76041号
[44] AFB汤普森;Gelhar,LW,三维随机非均质多孔介质中溶质运移的数值模拟,水资源。研究,26,10,2541-2562(1990)
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。