×

一种用于辐射传输问题的高效且精确的短特性求解器。 (英语) Zbl 1133.85004号

摘要:在S.Turek公司[计算54,No.1,27-38(1995;兹比尔0822.65129)]广义平均强度的概念被提出作为一种特殊的数值方法来处理(线性)辐射传输方程,它可以在不消除任何特定强度信息的情况下显著降低离散系统的维数。此外,结合Krylov空间方法(CG、Bi-CGSTAB等),开发了鲁棒且非常有效的解算器,作为经典近似(Lambda)迭代的扩展。本文的关键工具是将特殊的重新编号技术与有限差分离散化策略相结合,以产生的传输算子为基础,基于变阶短特征迎风技术,可应用于局部可变网格宽度的高度非结构化网格。我们演示了如何以一阶精度,特别是二阶精度构造这种特殊的上卷积格式,从而始终在一般网格上生成较低的三角形系统矩阵。因此,可以避免全局矩阵组装(“即时”),从而使存储成本几乎达到最优,并且可以非常有效地获得每个方向对应的对流-反应子问题的解。进一步的结果是,在无散射的情况下,这种方法可以直接求解,而在吸收和散射系数不为零的情况下整个系统的收敛速度仅取决于它们的比率和这些物理量的绝对大小,但仅在网格大小或网格拓扑上较弱。我们通过原型配置演示了这些结果,并针对不同的计算域、网格和问题参数检查了结果的准确性和效率。

MSC公司:

85A25型 天文学和天体物理学中的辐射传输
35A25型 适用于PDE的其他特殊方法

软件:

美食
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] Deo,N.:图论在工程和计算机科学中的应用。普伦蒂斯·霍尔(Prentice Hall),恩格伍德悬崖(Englewood Cliffs)(1974年),第222–333页·Zbl 0285.05102号
[2] Führer,C.:双曲问题有限元解算器与辐射传输应用的比较研究。技术报告SFB 359,65,海德堡大学(1993)
[3] Hübner,T.:Spezielle Diskretisierungs-und Lösungsmethoden für Integro-Differentialgleichungen am Beispiel der Strahlungstransportgleichung。毕业论文,多特蒙德(2005)。www.mathematik.uni-dortmund.de/lsiii/static/showpdffile_Huebner2005.pdf
[4] Johnson C.和Pitkäranta J.(1983年)。二维中子输运全离散格式的收敛性。SIAM J数字分析20:951–966·Zbl 0538.65097号 ·doi:10.1137/0720065
[5] Kalkofen W.(1987)。数值辐射传输。剑桥大学出版社,伦敦·Zbl 1170.85300号
[6] Klar,A.,Seaid,M.:辐射传输方程离散化引起的线性系统的有效预处理。《科学计算的挑战》,柏林,2002年。计算科学与工程讲稿(2003)
[7] Larsen E.W.(1982)。板几何离散坐标方程的无条件稳定扩散合成加速方法。核科学与工程82:47
[8] Larsen E.W.和Morel J.E.(1989)。光学厚扩散区数值输运问题的渐近解II。计算机物理杂志83:212·Zbl 0684.65118号 ·doi:10.1016/0021-9991(89)90229-5
[9] Meinköhn,E.:Modellierung dreidimensionaler Strahlungsfelder im frühen Universum。海德堡博士论文(2002)
[10] Mihalas D.和Weibel-Mihalas B.(1984年)。辐射流体力学基础。牛津大学出版社,纽约·Zbl 0651.76005号
[11] Richling S.、Meinköhn E.、Kryzhevoi N.和Kanschat G.(2001)。有限元辐射传输。天文天文图380:776–788·数字对象标识代码:10.1051/0004-6361:20011
[12] Schäfer,M.,Turek,S.(得到Durst,F.,Krause,E.,Rannacher,R.的支持):圆柱周围层流的基准计算。收录:Hirschel,E.H.(编辑)高性能计算机的流动模拟II。《数值流体力学注释》52、547–566。布伦瑞克·维埃格(1996)
[13] Sea id M.、Frank M.、Klar A.、Pinnau R.和Thömmes G.(2004年)。燃气轮机辐射的有效数值方法。《计算机应用数学杂志》170(1):217–239·Zbl 1221.80023号 ·doi:10.1016/j.cam.2004.01.003
[14] Steiner,O.:使用算符扰动的快速收敛温度校正程序。天文天体物理学231278-288。斯普林格,海德堡(1990)
[15] Turek S.(1993)。辐射传递方程的一种有效求解方法。影响计算科学工程5:201–214·兹比尔0786.65125 ·doi:10.1006/icse.1993.1009
[16] Turek S.(1995年)。辐射传递方程数值解的广义平均强度方法。计算54:27–38·兹比尔0822.65129 ·doi:10.1007/BF02238078
[17] Turek S.、Runge A.和Becker Ch.(2001)。FEAST指数——现代软件组件和处理器技术的现实评估。计算数学应用41:1431–1464·兹伯利0989.65051 ·doi:10.1016/S0898-1221(01)00108-0
[18] 范德福斯特H.(1992)。BI-CGSTAB:BI-CG的一种快速、平滑收敛的变体,用于求解非对称线性系统。SIAM科学统计杂志计算13:631–644·Zbl 0761.65023号 ·doi:10.1137/0913035
[19] Väth,H.M.:大规模并行计算机上的三维辐射传输。天文天体物理学284319。施普林格,海德堡(1994)
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。