多物理场耦合现象广泛存在于各种核反应堆中,尤其是先进的反应堆系统中。多物理场耦合以及大规模并行也是当前核能软件研发的核心课题。传统的基于各独立物理模块开发接口程序进行外耦合的方法不适用于强耦合、高精度、大规模并行的需求。多物理场耦合的本质是偏微分方程组的联立求解。基于统一编程框架进行多物理场耦合并对各个不同物理场进行建模、离散、数值求解也是开发数值反应堆的研究趋势。本文工作基于开源CFD平台OpenFOAM,通过对先进反应堆热工水力学与堆内其他物理场耦合现象的研究,探索先进核能系统中的多物理场耦合的方法以及关键技术难点。主要工作内容如下:1)基于OpenFOAM开发液态金属在强磁场中流动的热工水力学-磁场耦合程序,并应用于先进聚变反应堆液态包层热工水力学特性研究。采用电势法求解磁场-流体耦合的方程组,电势梯度和电流密度的求解采用相容守恒的离散格式,流体NS方程的求解采用PISO算法。对开发的磁流体动力学程序进行Shericliff、Hunt、ALEX方管和IEA圆管等国际基准算例验证。验证结果表明耦合程序具有非常好的准确性和可靠性,同时可以适用于非结构、非正交网格,具有较好的几何适应性。采用该耦合程序分析了聚变堆液态金属包层内冷却剂在磁场下的三维流动现象。对聚变堆液态金属包层中的关键突扩、突缩管道进行三维磁流体现象分析,从而提出结构优化,降低局部三维压降。另外当包层内液态金属流动速度比较低时,磁对流的现象会比较突出,针对所开发的磁场-流体动力学耦合程序进行功能扩展,并分析不同哈特曼数对三维竖直管道内的磁对流现象的影响。2)基于OpenFOAM平台开发点堆中子动力学模型与CFD热工水力学模型以及燃料棒瞬态导热模型之间的耦合程序,探索基于多套网格的采用常微分方程组和偏微分方程组描述的多物理耦合计算处理方法的实现。并分别对点堆中子动力学模型、CFD热工水力学模型中的多孔介质计算功能、自然对流功能和瞬态燃料棒导热模型三个模块进行分离验证以及三个模型耦合之后的整体验证,针对PDS-XADS的验证结果表明各个独立模块计算结果准确、耦合程序数据传递正确,可以用于反应堆的瞬态安全分析。3)为满足非对称冷却、局部反应性扰动等有功率空间变化的瞬态安全事故分析的需求,基于OpenFOAM平台,采用有限体积法开发扩散中子动力学程序,添加外推边界条件的实现,稳态特征值求解采用源迭代算法,瞬态求解时以稳态结果为初始状态,耦合缓发中子先驱核方程和温度反馈模型,进行时间步进求解,每个时间步更新各群中子宏观截面。针对所开发的程序采用临界球、临界柱和临界块算例以及IAEA、TWIGL、Dodd Benchmark、LRA和HTGR五个稳态特征值算例和瞬态基准算例来对扩散中子动力学程序进行验证。验证结果表明扩散中子动力学程序计算精度高且几何适应性强。随后将本节开发的扩散中子动力学与上一小节开发CFD热工水力学程序进行耦合,探索多物理场耦合之间的网格映射、迭代收敛等问题。4)针对开发的核-热耦合程序开展应用研究,探索基于面向对象编程框架开发的多物理场耦合程序的适应性。开展的应用研究有:一是小型自然循环铅冷快堆的不同寿期反应性引入事故的对比分析,观察不同寿期时由于燃耗等引起的安全系数变化对堆安全的影响;二是ADS次临界堆芯的典型束流瞬态事故分析,包括不同时长的失束、超束以及叠加事故情况下反应堆功率、反应性以及燃料温度等变化;三是全堆芯LRA非对称反应性引入事故和对称反应性引入事故的对比分析。本文基于面向对象的CFD程序框架开发了热工水力学-磁场耦合程序、基于点堆中子动力学和CFD热工水力学的耦合程序以及基于扩散中子动力学和CFD热工水力学的耦合程序,探索了不同网格之间物理场的映射方法,并将耦合程序作了一定的应用研究。掌握了基于统一编程框架进行多物理场之间通过方程源项或方程系数耦合以及不同网格之间物理场映射的关键技术。本文的相关研究成果有助于我国核能软件的发展,为基于统一编程框架的先进核能系统的多物理场耦合程序开发提供经验积累和技术支持,具有一定的研究意义和工程应用价值。