- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
附录三撰写论文具体要求及论文模板-核工程计算物理试验室
撰写论文具体要求(参考论文模板)
论文录入请用word文档,纸张为A4,页边距取默认上下2.54厘米,左右3.17厘米,页面设置无网格;
标题:小二号黑体,单倍行距,居中;
作者姓名:小四号仿宋,居中,单倍行距;通信地址:10磅仿宋,居中,单倍行距;
脚注:包括基金资助项目名称和项目号及第一作者简介。“基金项目:”、“作者简介:”小五号宋体加粗,脚注内容小五号宋体,左对齐,单倍行距。
“摘要”:小五号宋体加粗,摘要内容:小五号宋体,1.5倍行距,两端对齐;
“关键词”:小五号宋体加粗,关键词内容:小五号宋体,1.5倍行距;
论文分节编码:阿拉伯数字分级编号,最高3级。黑体,两端对齐,字号(1级四号、2级小四号、3级五号)。
论文正文:五号宋体,两端对齐,首行缩进2字符,单倍行距,不分栏;
图题、表题:图表依次编号,小五号宋体,居中对齐;
公式:所有公式以“()”依次编号,右对齐,公式主体位于中心,单倍行距;
“参考文献”:小四号黑体加粗; 参考文献:以“[]”依次编号,悬挂缩进2字符,10磅宋体,两端对齐。
全文所有英文字母采用 Time New Roma字体。
论文以“姓名_单位_论文名称”命名。
请将“保密审查证明”附在征文后一并提交。
为了统一格式,将论文汇编成册,请您在编写论文时务必遵守通知中的规定。谢谢您的合作!
以下为参考模板。
结合多群GMRES的IRAM算法
用于求解矩阵MOC方程(
吴文斌1,2 李庆1 王侃2
(1中国核动力研究设计院核反应堆系统设计技术重点实验室
摘要:矩阵MOC方法通过构造并求解线性方程组,代替传统MOC方法中反复地特征线扫描。求解中子输运方程临界问题时,通常采用幂迭代法求得keff。然而幂迭代法的收敛速度严重依赖于占优比,实际的较大规模的堆芯占优比接近于1,收敛很慢。本研究采用隐式再启动的Arnoldi算法(IRAM)求解keff,并应用多群耦合的GMRES算法直接求解含有上散射的多群问题,以避免能群间的散射迭代。采用C++语言编写了相关计算程序,对多个基准题如2D C5G7的数值结果表明,和幂迭代法相比,结合多群GMRES的IRAM算法具有良好的计算精度和更高的计算效率。
关键词:矩阵MOC;幂迭代法;IRAM算法;多群耦合GMRES
引言
由于具有较高的计算精度和较强的几何适应性,特征线方法(MOC)被认为是一种极具潜力的确定论方法,被广泛应用于中子输运计算。传统MOC方法求解中子输运方程需要反复地进行特征线扫描,消耗大量的计算时间。矩阵MOC方法仅通过一次扫描,即可构造出与特征线扫描等效的线性方程组;随后,求解该线性方程组可代替反复地特征线扫描。该线性方程组具有良好的数值特性,有很多高效求解器可供使用。
理论模型
矩阵MOC方法
通过一次特征线扫描,矩阵MOC方法可构造出与扫描完全等效的线性代数方程组
方程组的系数矩阵A和B具有如下的数值特性:
1)采用CSR格式存储系数矩阵时,不同能群的系数矩阵稀疏结构一致。
2)矩阵A和矩阵B具有比例关系。
3)系数矩阵A可由基本矩阵经初等变换得到
4)互易关系,矩阵a的对称性。
多群耦合GMRES方法
在处理多群中子输运问题时,如果存在上散射,通常采用Gauss-Seidel迭代处理能群间的上散射耦合,具体过程是:从高能群到低能群逐群求解群内问题(within-group problem),右端源项采用Gauss-Seidel迭代不断更新。
(以下内容省略)
数值验证
按照第2节的理论分析,采用C++语言编制了结合多群GMRES的IRAM算法求解矩阵MOC方程的程序。为了分析该程序的计算精度和计算效率,对BWR栅元问题、C5G7 UO2组件和C5G7 2D基准题进行了计算。所有的计算都是在个人电脑上完成的,该电脑采用Windows XP 32位系统,处理器是Intel Core 2 Quad Q8300 2.5GHz,内存空间为2G。尽管处理器包含4个计算核心,但是仅有一个核心参与计算。IRAM算法的收敛准则是10-5。矩阵MOC 的计算条件是:16个方位角,2个Leonard 最优极角,每个方位角拥有200条特征线,射线间距约为0.01cm。
BWR栅格基准题
该基准题是一个4×4的BWR栅格问题,包含两根Gd棒,如图 1所示。燃料棒1,2,3,4,5是3%w的UO2棒,毒物棒6包含3%w的UO2和3%w的Gd2O3。燃料棒的包壳为Zircaloy-2,慢化剂为水。外边界均采用全反射边界条件。
BWR栅格基准题几何布置和栅元网格划分
keff和归一化的燃料棒功率分布列于表 1,IRAM与幂迭代的计算结果保持一致,且和基准结果都符合得很好。
矩阵MOC求解BWR栅格基准题的计算结果
DRAGON (referen
文档评论(0)