×

(SP_N)方程的块Arnoldi方法。 (英文) Zbl 1480.65346号

摘要:简化的球谐方程是稳态中子输运方程的有效近似。从计算的角度来看,与它们相关的特征值问题是一个具有挑战性的问题。在这项工作中,我们利用所涉及矩阵的块结构,提出了无块逆预处理Arnoldi方法,作为解决该特征值问题的有效方法。对于空间离散化,使用了无矩阵技术实现的连续Galerkin有限元方法来保持合理的内存需求。为了提高有限元方法的收敛性,提出了一种使用线性形状函数的多级初始化方法。此初始化只占用总计算时间的一小部分。将所提出的特征值求解器与标准幂迭代法、Krylov-Schur方法和广义Davidson方法进行了比较。数值结果表明,它减少了求解特征值问题的计算时间。

MSC公司:

65N30型 含偏微分方程边值问题的有限元、Rayleigh-Ritz和Galerkin方法
65N25型 含偏微分方程边值问题特征值问题的数值方法
35立方英尺47英寸 二阶椭圆系统
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] 阿本茨,P。;赫特马纽克,美国。;Lehoucq,R.B。;Tuminaro,R.S.,《使用amg预处理迭代方法进行大规模三维模态分析的特征解算器比较》,国际期刊编号。《工程方法》,64204-236(2005)·Zbl 1093.74021号
[2] Baglama,J.,增广块householder-arnoldi方法,线性代数应用。,429, 2315-2334 (2008) ·Zbl 1153.65034号
[3] 班杰斯,W。;哈特曼,R。;Kanschat,G.,交易。II-通用面向对象有限元库,ACM Trans。数学。软质。,33, 1-27 (2007) ·Zbl 1365.65248号
[4] 绑定,P。;Najman,B。;Ye,Q.,厄米矩阵铅笔特征值的变分原理,积分。Equat公司。运营商。理论。,35, 398-422 (1999) ·Zbl 0971.47008号
[5] Brantley,P.S。;Larsen,E.W.,简化的P3近似,Nucl。科学。技术。,134, 1-21 (2000)
[6] 卡皮拉,M。;塔拉维拉,C。;Ginestar,D。;Verdü,G.,多维PL方程的节点配置近似在3D takeda基准问题中的应用,Ann.Nucl。能源,40,1-13(2012)
[7] 卡尼,S。;布朗,F。;基德罗夫斯基,B。;Martin,W.,《连续能量蒙特卡罗裂变矩阵法的理论和应用》,Ann.Nucl。《能源》,73423-431(2014)
[8] Carreño,A。;Vidal-Ferrándiz,A。;Ginestar,D。;Verdü,G.,计算中子扩散方程主要λ-模的块混合多级方法,Ann.Nucl。能源,121,513-524(2018)
[9] Gelbard,E.M.,球形谐波方法在反应堆问题中的应用,技术代表WAPD-BT-20,Bettis原子能实验室,1960年。
[10] Golub,G.H。;Ye,Q.,对称广义特征值问题的无逆预处理krylov子空间方法,SIAM J.Sci。计算。,24, 312-334 (2002) ·Zbl 1016.65017号
[11] 汉密尔顿,S.P。;Evans,T.M.,简化PN方程的有效解,J.Compute。物理。,284, 155-170 (2015) ·Zbl 1351.82123号
[12] 汉密尔顿,S.P。;埃文斯,T.M。;Davidson,G.G。;约翰逊·S·R。;潘迪亚,T.M。;Godfrey,A.T.,《使用非逻辑代码计算热零功率反应堆》,J.Compute。物理。,314, 700-711 (2016) ·兹比尔1349.82086
[13] 埃尔南德斯,V。;罗曼,J.E。;Vidal,V.,SLEPc:用于解决特征值问题的可扩展且灵活的工具包,ACM Trans。数学。软质。,31, 351-362 (2005) ·Zbl 1136.65315号
[14] Knyazev,A.V.,预处理特征解算器-矛盾修饰?,电子。事务处理。数字。Ana,7104-123(1998)·Zbl 1053.65513号
[15] Kronbichler,M。;Kormann,K.,基于并行单元的有限元算子应用程序的通用接口,Comput。流体,63,135-147(2012)·Zbl 1365.76121号
[16] 刘易斯,E.E。;Miller,W.F.J.,《中子输运计算方法》(1984),John Wiley&Sons,Ltd.:John Willey&Sons有限公司,美国纽约
[17] 摩根R.B。;Scott,D.S.,davidson稀疏对称矩阵特征值计算方法的推广,SIAM J.Sci。和统计计算。,7, 817-825 (1986) ·Zbl 0602.65020号
[18] 奎伦,P。;Ye,Q.,对称广义特征值问题的无块逆预条件krylov子空间方法,J.Comput。申请。数学。,233, 1298-1313 (2010) ·Zbl 1186.65044号
[19] Saad,Y.,《稀疏线性系统的迭代方法》(2003),工业和应用数学学会:美国费城工业与应用数学学会·Zbl 1031.65046号
[20] Smith,M.A.、Lewis,E.E.和Na,B.C.,《无空间均匀化确定性运输计算基准——二维/三维MOX燃料组件基准》(C5G7 MOX基准),技术代表NEA/NSC/DOC(2003)16,OECD/NEA,2003。
[21] Stacey,W.M.,《核反应堆物理》(2007),John Wiley&Sons:John Willey&Sons,Weinheim,德国
[22] Stewart,G.W.,大型特征问题的krylov-schur算法,SIAM J.矩阵分析。申请。,23, 601-614 (2002) ·Zbl 1003.65045号
[23] Verdú,G。;Ginestar,D。;维达尔,V。;Muñoz-Cobo,J.,《中子扩散方程的三维λ模》,Ann.Nucl。能源,21405-421(1994)
[24] Vidal-Ferrándiz,A。;Fayez,R。;Ginestar,D。;Verdü,G.,《使用h-p有限元法求解核动力反应堆的lambda模式问题》,Ann.Nucl。能源,72,338-349(2014)
[25] Vidal-Ferrándiz,A。;González-Pintor,S。;Ginestar,D。;Verdú,G。;Demazière,C.,中子扩散方程的Schwarz型预条件,J.Compute。申请。数学。,309, 563-574 (2017) ·Zbl 1468.65180号
[26] 瓦萨,J.S。;Wareing,T.A。;莫雷尔,J.E。;McGee,J.M。;Lehoucq,R.B.,确定性k特征值计算的Krylov子空间迭代,Nucl。科学。技术。,147, 26-42 (2004)
[27] Wu,K。;萨阿德,Y。;Stathopoulos,A.,大型对称特征值问题的非精确牛顿预处理技术,电子。事务处理。数字。分析。,7, 202-214 (1998) ·兹比尔0916.65035
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。