三维全堆芯pin-by-pin中子输运模型的高效加速方法是核反应堆高精度计算的重点和难点。本文有效融合课题组开发的并行多维离散纵坐标(S_(N))中子输运程序comeSn和Jacobian-Free Newton Krylov(JFNK)通用求解框架comeJFNK的高效并行特性...三维全堆芯pin-by-pin中子输运模型的高效加速方法是核反应堆高精度计算的重点和难点。本文有效融合课题组开发的并行多维离散纵坐标(S_(N))中子输运程序comeSn和Jacobian-Free Newton Krylov(JFNK)通用求解框架comeJFNK的高效并行特性、鲁棒性和强收敛性,开发了一套三维稳态及瞬态中子输运模型的JFNK并行求解程序comeSn_JFNK。为了提高计算效率,选择中子标通量密度(而非中子角通量密度)作为JFNK全局求解变量,并利用基于空间区域并行的KBA输运扫描方法和物理预处理方法,分别构建了稳态及瞬态模型的JFNK统一残差计算模型。计算结果表明,comeSn_JFNK相比于comeSn,计算效率具有显著优势,对于三维pin-by-pin稳态KAIST-3A算例,加速比为10倍以上;对于栅元均匀化的二维七群瞬态C5G7-TD2系列基准算例,加速比约为30倍。展开更多
为了提高可压缩雷诺平均Navier-Stokes(RANS)方程的求解效率,基于多块对接结构化网格发展了求解RANS方程的Jacobian-Free Newton Krylov(JFNK)方法.JFNK方法将求解非线性方程的非精确Newton法和求解线性方程的Krylov子空间迭代法结合,...为了提高可压缩雷诺平均Navier-Stokes(RANS)方程的求解效率,基于多块对接结构化网格发展了求解RANS方程的Jacobian-Free Newton Krylov(JFNK)方法.JFNK方法将求解非线性方程的非精确Newton法和求解线性方程的Krylov子空间迭代法结合,通过非精确Newton方法中不精确条件控制不同阶段线性方程的求解精度,并利用无矩阵技术求解矩阵与向量的乘积;针对Krylov内迭代收敛停滞的问题,引入LU-SGS方法作为预处理器,降低线性系统的刚性从而大幅度提高了内迭代的计算效率.利用JFNK方法模拟NACA 0012翼型、带襟翼的NLR-7301两段翼与带发动机短舱的DLR-F6翼身组合体的绕流问题,研究不同参数对JFNK方法收敛特性的影响,对比LU-SGS研究JFNK方法的收敛速度,并对JFNK方法求解复杂绕流问题的RANS方程进行确认.结果表明,JFNK方法求解RANS方程具有良好的稳定性,相对于其他时间推进方法,JFNK方法具有更高的计算效率.展开更多
文摘三维全堆芯pin-by-pin中子输运模型的高效加速方法是核反应堆高精度计算的重点和难点。本文有效融合课题组开发的并行多维离散纵坐标(S_(N))中子输运程序comeSn和Jacobian-Free Newton Krylov(JFNK)通用求解框架comeJFNK的高效并行特性、鲁棒性和强收敛性,开发了一套三维稳态及瞬态中子输运模型的JFNK并行求解程序comeSn_JFNK。为了提高计算效率,选择中子标通量密度(而非中子角通量密度)作为JFNK全局求解变量,并利用基于空间区域并行的KBA输运扫描方法和物理预处理方法,分别构建了稳态及瞬态模型的JFNK统一残差计算模型。计算结果表明,comeSn_JFNK相比于comeSn,计算效率具有显著优势,对于三维pin-by-pin稳态KAIST-3A算例,加速比为10倍以上;对于栅元均匀化的二维七群瞬态C5G7-TD2系列基准算例,加速比约为30倍。
文摘为了提高可压缩雷诺平均Navier-Stokes(RANS)方程的求解效率,基于多块对接结构化网格发展了求解RANS方程的Jacobian-Free Newton Krylov(JFNK)方法.JFNK方法将求解非线性方程的非精确Newton法和求解线性方程的Krylov子空间迭代法结合,通过非精确Newton方法中不精确条件控制不同阶段线性方程的求解精度,并利用无矩阵技术求解矩阵与向量的乘积;针对Krylov内迭代收敛停滞的问题,引入LU-SGS方法作为预处理器,降低线性系统的刚性从而大幅度提高了内迭代的计算效率.利用JFNK方法模拟NACA 0012翼型、带襟翼的NLR-7301两段翼与带发动机短舱的DLR-F6翼身组合体的绕流问题,研究不同参数对JFNK方法收敛特性的影响,对比LU-SGS研究JFNK方法的收敛速度,并对JFNK方法求解复杂绕流问题的RANS方程进行确认.结果表明,JFNK方法求解RANS方程具有良好的稳定性,相对于其他时间推进方法,JFNK方法具有更高的计算效率.
