matlab加速迭代法代码-dsdoe:18.086项目离散Ordinates方程的直接解决方案的论文和代码

matlab加速迭代法代码dsdoe 18.086项目离散Ordinates方程的直接解决方案的论文和代码 抽象的 本文的目的是研究和实现一维和二维多组中子输运的离散坐标方法，并直接求解所得方程。 从历史上看，离散纵坐标方程是通过幂迭代方法（PI）求解的。 通常实现的方法非常有效地提高内存效率，但是对于高度分散的问题，PI方法非常无效。 已经开发了几种加速PI方法的技术，但是即使对于简单的几何形状，它们也很复杂。 在这里，研究了一种替代方法，即直接求解离散坐标系。 一维和二维多维群离散坐标矩阵的一般形式已在MATLAB中使用其内置的稀疏矩阵实用程序导出并实现。 为了解决该系统，研究了直接消除的标准“反斜杠”以及内置Iry预处理的Krylov求解器（GMRES和BiCG-STAB）的组合。 将这些直接方法的计算费用与标准PI方法的计算费用进行了比较，标准PI方法也已在Matlab中实现，并已针对生产离散纵坐标代码进行了验证。 一维研究表明，对于低角度阶次（因此带宽较小），直接方法（使用消除器或Krylov求解器）比PI方法快大约两倍。 对于高度扩散的问题，直接方法的速度提高了一个数量级。