[发明专利]基于混合计算架构的用于核素迁移控制方程的求解方法

权利要求书

1.一种基于混合计算架构的用于核素迁移控制方程的求解方法，其特征在于，其包括以下步骤：

S1：读取参数阶段：读入计算参数、边界条件和计算网格文件；

S2：生成矩阵阶段：生成刚度稀疏矩阵、时间演化稀疏矩阵和载荷向量，并生成对应的稀疏矩阵，之后刚度稀疏矩阵、时间演化稀疏矩阵；

S3：循环计算求解阶段：使用有限差分法循环计算，循环过程中采用Jacobi迭代法求解AX＝B的矩阵方程，得到的解向量X即为该时间点的核素浓度场数据。

2.如权利要求1所述的基于混合计算架构的用于核素迁移控制方程的求解方法，其特征在于，所述核素在缓冲/回填材料中迁移控制方程为：

A∂C∂t=∂∂xi(Dij∂C∂xj)-Vi∂C∂xi-FC+Q]]>

其中C为核素浓度，D_ij为核素在多孔介质中的动力扩散系数：D_ij＝Di_j'+D_ij″，Dij′=αTV‾δij+(αL-αT)ViVjV‾,]]>

α_T，α_L分别为横向和纵向弥散系数，为渗流场速度大小，V_j为j方向的渗流速度分量，D_ij″＝D_dTδ_ij，D_d为分子扩散系数，T为多孔介质的弯曲率；

这里D_d＝D₀exp(σΔT)，D₀为参考温度下的分子扩散系数，σ为温度修正系数，ΔT为相对于参考温度的温度差，F为等效衰减系数，Q为源汇项。