三维全堆芯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倍。展开更多
文摘三维全堆芯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倍。