面向反应堆燃料棒束热工水力仿真的流固耦合并行方法

1.一种面向反应堆燃料棒束热工水力仿真的流固耦合并行方法，其特征在于，所述方法包括：根据反应堆的燃料棒的空间位置将整体计算域划分为不同的格子，每个格子包含一根完整的燃料棒，并将每个格子分配给一个单独的计算核心，对每个计算核心分别加载固体域网格和流体域网格并初始化变量场；所述变量场包括温度、速度和压力；调用固体求解器在固体域网格上求解固体传热模型；对所述固体传热模型中热流密度和固体温度进行判断，根据判断结果传递热流密度到流体域网格；调用流体求解器计算流体域伤口交界面处的通量，根据所述通量计算得到流体域伤口交界面处网格单元上的速度、压力和温度并将速度、压力和温度发送给邻居计算核心，并从邻居计算核心接收伤口交界面另一侧网格单元上的变量场；对所述速度、压力和温度进行收敛判断，若收敛，则将所述温度通过流固交界面发送到流体域网格；对所述固体传热模型中热流密度和固体温度进行判断，根据判断结果传递热流密度到流体域网格，包括：对所述固体传热模型中热流密度和固体温度进行判断，判断热流密度和固体温度是否达到收敛精度，如果未收敛，则在固体域网格上求解固体传热模型；如果收敛，则将热流密度通过流固交界面传给流体域网格。 2.根据权利要求1所述的方法，其特征在于，在固体域网格上求解固体传热模型，包括：在固体域网格上求解固体传热模型为：其中， q 是热流密度， 为固体温度， 为固体的热导率， 为热源。 3.根据权利要求1至2任意一项所述的方法，其特征在于，计算流体域伤口交界面处的通量，包括：计算流体域伤口交界面处的通量为：其中， 为质量流量， 为扩散项系数， 和 分别为面矢量 非正交修正的两个分量， 为 的模； 为待求变量， 为计算单元与邻居形心连线的向量， 为它的模， 为当前计算单元体积， 为密度， 表示时间步长， 表示上一时间步密度， 表示上一时间步当前离散单元变量值， 表示方程源项， 表示对流项和扩散项计算得到的对当前离散单元产生影响的通量， 表示对流项和扩散项计算得到的对当前离散单元的邻居单元产生影响的通量， 表示对流项和扩散项计算得到的对矩阵源项产生影响的通量， 表示瞬态项计算得到的对当前离散单元产生影响的通量， 表示瞬态项和源项计算得到的对矩阵源项产生影响的通量。 4.根据权利要求3所述的方法，其特征在于，根据所述通量计算得到流体域伤口交界面处网格单元上的速度、压力和温度，包括：根据所述通量依次求解得到流体域伤口交界面处网格单元上的速度、压力和温度的表达式为：其中， 表示当前离散单元系数， 表示当前离散单元的邻居单元系数， 表示矩阵源项， 表示当前离散单元求解变量，变量包括速度、压力和温度， 表示当前离散单元的邻居单元求解变量， 表示邻居， 表示邻居数。