(19)国家知识产权局 (12)发明 专利 (10)授权公告 号 (45)授权公告日 (21)申请 号 202210919848.9 (22)申请日 2022.08.02 (65)同一申请的已公布的文献号 申请公布号 CN 114970293 A (43)申请公布日 2022.08.30 (73)专利权人 西安交通大 学 地址 710049 陕西省西安市碑林区咸宁西 路28号 (72)发明人 刘宙宇　周欣宇　吴宏春　曹良志　 (74)专利代理 机构 西安智大知识产权代理事务 所 61215 专利代理师 何会侠 (51)Int.Cl. G06F 30/23(2020.01) G06F 30/25(2020.01) G06F 17/14(2006.01)G06F 111/10(2020.01) (56)对比文件 CN 114817834 A,202 2.07.29 CN 114387421 A,202 2.04.22 JP 2008051509 A,2008.03.06 CN 114528719 A,202 2.05.24 高明敏.二维MOC与一维NE M耦合堆芯计算方 法研究. 《中国优秀硕士学位 论文全文数据库 工 程科技II辑》 .2018, Pei Cao等.A 3-D Neutro n Distributi on Reconstructi on Method Based o n the Of f- Situ Measurement for Reactor. 《IE EE Transacti ons on Nuclear Science 》 .2021, 高智.对流扩散方程的绝对稳定高阶中心差 分格式. 《力学 学报》 .2010, 审查员 刘莹莹 (54)发明名称 基于二维一维耦合傅里叶分析的粗网扩散 系数计算方法 (57)摘要 本发明公开了基于二维一维耦合傅里叶分 析的粗网扩散系数计算方法， 将所要模拟的压水 堆三维堆芯划分成若干个粗网； 在每个粗网中划 分多个细网； 在每个细网上建立二维一维耦合输 运模型， 在每个粗网上建立粗网有限差分模型； 根据二维一维耦合输运模型和粗网有 限差分模 型建立二维一维傅里叶分析模型， 获得谱半径最 小时， 随粗网光学厚度、 散射率变化的粗网扩散 系数的修正因子； 基于修正因子计算修正粗网有 限差分模型的粗网扩散系数， 求解修正粗网有限 差分模型， 获得粗网中子标通量密度及特征值， 并更新细网中子标通量密度及特征值实现加速。 本发明方法使粗网有 限差分加速二维一维耦合 输运计算的过程更稳定， 可用于压水堆数值反应 堆模拟计算及工程应用。 权利要求书9页 说明书17页 附图2页 CN 114970293 B 2022.10.11 CN 114970293 B 1.一种基于二维一维耦合傅里叶分析的粗网扩散系数计算方法， 其特征在于： 包括如 下步骤： 步骤1： 读取所要模拟的压水堆三维堆芯的几何信息、 材料信 息、 边界条件信息、 细群及 宽群的能群结构信息； 步骤2： 根据步骤1获取的压水堆三维堆芯的几何信息和材料信息， 将压水堆三维堆芯 划分成若干个粗网； 在每 个粗网中划分多个细网； 步骤3： 根据步骤2划分的每个细网， 建立二维一维耦合输运模型， 如公式 （1） 表示二维 输运模型、 公式 （2） 表示一 维输运模 型， 二维输运模型和一 维输运模 型通过粗网的表面中子 角通量密度耦合； 求解二 维一维耦合输运模型获得细网的中子标通量密度和粗网表面的中 子角通量密度， 并将粗网表面的中子角通量密度投影到粗网表面的法线方向获得中子流密 度； （1） （2） 其中， i—细网沿 x轴方向的编号； j —细网沿 y轴方向的编号； k —细网沿 z轴方向的编号； ix —粗网沿 x轴方向的编号； iy—粗网沿 y轴方向的编号； iz—粗网沿 z轴方向的编号； h —细网沿 x,y轴方向的尺寸； hz—细网沿 z轴方向的尺寸； Δ—粗网沿 x,y轴方向的尺寸；权　利　要　求　书 1/9 页 2 CN 114970293 B 2Δz—粗网沿 z轴方向的尺寸； m2d—二维输运方程的内迭代次数； m1d—一维输运方程的内迭代次数； M2d—二维输运方程的最大内迭代次数； M1d—一维输运方程的最大内迭代次数； l —二维输运方程和一维输运方程的外迭代次数； p —中子角通 量密度的角度编号； —方向p与x轴正方向的夹角的余弦值； —方向p与y轴正方向的夹角的余弦值；  —方向p与z轴正方向的夹角的余弦值； —压水堆 三维堆芯的中子总截面， t为中子总截面符号； —压水堆 三维堆芯的中子 散射截面， s为中子 散射截面符号；  —压水堆 三维堆芯的中子裂 变截面， f为中子裂 变截面符号； —平均每次裂 变中子数； —第l迭代步的特 征值； —压水堆三维堆芯的第( l+1/2)次外迭代， 第 m2d次内迭代， 位于 z轴方向 第iz个粗网层内的 x轴方向第 i个、y轴方向第 j个细网的西面的沿方向 p的中子角通 量密度； —压水堆三 维堆芯的第( l+1/2)次外迭代， 第 m2d次内迭代， 位于 z轴方向 第iz个粗网层内的 x轴方向第 i个、y轴方向第 j个细网的东 面的沿方向 p的中子角通 量密度； —压水堆 三维堆芯的第( l+1/2)次外迭代， 第 m2d次内迭代， 位于 z轴方向第 iz个粗网层内的 x轴方向第 i个、y轴方向第 j个细网的南 面的沿方向 p的中子角通 量密度； —压水堆三维堆芯的第( l+1/2)次外迭代， 第 m2d次内迭代， 位于 z轴方向权　利　要　求　书 2/9 页 3 CN 114970293 B 3