隐式重启的Arnoldi方法及其在高阶谐波求解中的应用.PDFVIP

　第 28 卷 　第 2 期 核 科 学 与 工 程 Vol . 28 　No . 2 　2008 年 　6 月 Chine se J our nal of N uclear Science and Engineering J un . 　2008 隐式重启的 Arnol di 方法 及其在高阶谐波求解中的应用 黄义超 ,张少泓 (上海交通大学核科学与工程学院 ,上海 200240) 摘要 : Krylov 子空间方法的出现是近年来大型线性方程组和特征值问题求解领域的重大进展 ,介绍其中 ( ) 一类适用于求解反应堆 k本征值问题的隐式重启的 Arnoldi 方法 IRA M , 以及该方法在高阶谐波求解 中的应用 。研究结果表明 ,IRA M 方法求解高阶谐波具有和源修正法同样的精度 ,但计算速度更快 ,尤 其是当所求的谐波阶次较高时 , IRA M 方法可获得 10 倍以上的速度优势 。同时 , IRA M 方法还具备较 好的处理重特征值问题的能力 。 关键词 : Krylov 子空间方法 ;Arnoldi 方法 ;高阶谐波 ( ) 中图分类号 : TL 329 　　　文献标识码 :A 　　　文章编号 :025809 18 2008 020 10205 Impl icitly restarted Arnol di method and its appl ication to the sol ution of higher order harmonics HU AN G Yichao , ZHAN G Shaohong ( Shanghai J iaotong U niver sit y , Shanghai 200240 , China) Abstract :The recent develop ment of Krylov sub sp ace met ho d i s a maj or p ro gre ss in t he area of lar gescale linear equation and eigenvalue p roblem re search . One t yp e of Krylov sub sp ace met ho d suit able for eigenvalue p roblem i s Imp licit ly Re st art ed A r nol di Met ho d ( IRA M ) it i s int ro duced and app lied to t he solutio n of higher or der har mo nic s. N umerical re sult s demon st rat e t hat w hile it s accuracy i s co mp arable wit h t he conven tional mo dified source it eration met ho d (M SIM) , IRA M met ho d r un s much f a st er t han