摘要

热工水力模型（Thermal-Hydraulics Model）是核反应堆安全分析、设计优化和运行支持的基石，其准确性与可靠性直接关系到核电站的安全边界和经济效益。报告将系统性地梳理热工水力模型在核工业中的具体实现方式，追溯其从20世纪中期至今的发展历程与关键里程碑。同时，本报告将重点分析该领域存在的核心技术争议，如“保守估算法”与“最佳估算加不确定性（BEPU）”的路线之争、模型验证与确认（V&V）的复杂性、以及多尺度/多物理场耦合的挑战，并阐述监管机构、工业界和学术界在这些争议中的各自立场。

第一章：引言与核心概念

1.1 研究背景与重要性

核能作为一种清洁、高效的能源形式，在全球能源结构中扮演着日益重要的角色。然而，核能的利用始终伴随着对安全的极高要求。核反应堆是一个极端复杂的能量系统，其核心是通过核裂变产生巨大的热能，再通过冷却剂将热能导出，最终驱动汽轮机发电。在这个过程中，冷却剂的流动、传热和相变行为，即“热工水力”现象，是整个能量转换链条的物理基础 。

热工水力模型正是为了描述和预测这些复杂现象而构建的数学和物理模型集合。它不仅是反应堆设计阶段进行性能评估和优化的关键工具，更是运行阶段进行安全分析、事故模拟和制定应急预案的根本依据。无论是正常运行工况下的热功率分布，还是如丧失冷却剂事故（LOCA）等极端条件下的系统响应，都必须通过精确的热工水力模型进行模拟预测 。因此，热工水力模型的能力和精度，直接决定了我们对核反应堆安全状态的认知深度和信心，是核安全审评和许可过程中的核心审查内容 。

1.2 热工水力模型的核心定义

热工水力模型是一个复杂的计算框架，其核心任务是求解描述流体运动和热量传递的控制方程组（质量、动量、能量守恒方程）。具体而言，它负责计算冷却剂在反应堆一回路（Primary Loop）和二回路（Secondary Loop）中的宏观与微观行为 。其模拟范围覆盖了从燃料芯块内部的热量产生，到通过燃料包壳传递给冷却剂，再到冷却剂流经堆芯、管道、蒸汽发生器等关键设备，最终将热量传递给二回路工质的全过程 。

该模型必须能够处理：

•流动（Hydraulics）： 冷却剂的速度场、压力场分布，包括单相流（液体或气体）和复杂得多的两相流（如气液混合流动、沸腾与凝结）。

•传热（Thermal）： 固体（燃料、包壳、管道壁）与流体之间的热量交换，包括导热、对流和辐射（尤其在高温气冷堆中）。

•相变（Phase Change）： 在特定压力和温度下，冷却剂可能发生的沸腾（液相到气相）或凝结（气相到液相）过程，这对传热效率和系统压力有剧烈影响。

这些过程相互耦合、瞬息万变，使得热工水力建模成为一项极具挑战性的工作。

1.3 模型在核反应堆设计与安全分析中的作用

热工水力模型在核电行业的应用贯穿于核电站的全生命周期。

•设计阶段： 工程师利用热工水力模型优化堆芯结构、燃料组件排布、冷却剂流道设计，以确保在额定功率下有足够的热量导出能力和安全裕度。例如，预测燃料包壳的最高温度（Peak Cladding Temperature, PCT）是否低于材料极限。

•安全分析与许可申请： 这是热工水力模型最关键的应用领域。监管机构要求电站运营商对一系列假想事故（设计基准事故）进行模拟分析，以证明电站的固有安全特性和专设安全系统（如紧急堆芯冷却系统 ECCS）的有效性 。热工水力模型是这些分析的核心工具，其计算结果是判断是否满足监管要求（如10 CFR 50.46中关于PCT和氧化的规定）的直接依据 。模型的计算结果，特别是关键位置（如燃料包壳表面、压力容器壁面）的温度和压力随时间变化的历程，是进行设备结构完整性和疲劳寿命分析的直接输入。

•运行支持与优化： 在电站运行期间，模型可用于分析异常瞬态、评估操作策略变更的影响、优化燃料管理方案，甚至用于开发操作员培训模拟器。

•先进反应堆研发： 对于第四代核能系统（Gen-IV）等新型反应堆，由于缺乏运行经验和实验数据，热工水力模型更是探索其可行性、发现潜在安全问题、指导实验设计的不可或缺的工具 。

第二章：热工水力模型的历史演进与里程碑

核反应堆热工水力学的发展史，是一部与计算机技术进步、核安全理念深化以及重大核事故教训紧密交织的历史。

2.1 早期探索与理论奠基 (1950s - 1960s)

核能发展的初期，即20世纪50至60年代，热工水力分析主要依赖于简化的解析模型和经验公式 。当时的计算机能力极为有限，无法支持复杂的数值模拟。研究重点集中在单相流传热、临界热流密度（CHF）的实验研究和经验关系式的建立上。这一时期的代码，如1966年发布的RELAP1，是基于简单一维模型的早期尝试 。随着商用核电站的设计和建造，特别是轻水堆（LWR）中两相流现象的重要性日益凸显，对更复杂模型的理论需求开始萌芽 。1957年至1968年间，苏联物理与动力工程研究所（IPPE）等机构开展了大规模的热工水力基础研究，为后续发展奠定了坚实的物理和实验基础 。

2.2 安全需求的驱动与系统代码的诞生 (1970s)

20世纪70年代是热工水力模型发展的黄金时期，其主要驱动力来自于对核反应堆安全的日益关切，特别是对紧急堆芯冷却系统（ECCS）有效性的担忧 。美国核管会（USAEC，后来的NRC）要求对失水事故（LOCA）进行详尽分析，这极大地推动了能够模拟整个反应堆冷却剂系统复杂瞬态行为的“系统代码”（System Codes）的开发 。

这一时期诞生了至今仍有深远影响的两大代码家族：

•RELAP系列： 由爱达荷国家实验室（INL）为美国能源部（DOE）开发，RELAP5的第一个版本于1979年发布，旨在模拟压水堆（PWR）的大破口LOCA 。

•TRAC系列： 由洛斯阿拉莫斯国家实验室（LANL）为NRC开发，其开发始于1975年，旨在成为更先进、更贴近物理真实情况的最佳估算代码 。

这些代码虽然受限于当时的计算能力，在物理模型和数值方法上做出了妥协（例如，使用集总参数法和一维模型），但它们首次实现了对整个一回路系统的综合模拟，是热工水力分析从局部经验公式到全局系统模拟的革命性转变。

1979年的三哩岛（Three Mile Island, TMI）事故成为该领域发展的又一个分水岭。这次事故并非设计基准的大破口LOCA，而是一个由小故障引发的复杂瞬态序列。事故暴露出当时的代码对于小破口LOCA和非预期两相流现象的模拟能力严重不足 。TMI事故的深刻教训，促使全球核工业界和监管机构投入巨量资源，加强对复杂两相流现象的实验研究和模型开发，系统代码的验证和确认（V&V）工作被提升到前所未有的高度 。

2.3 最佳估算方法与不确定性分析的兴起 (1980s - 1990s)

TMI事故后，核安全理念发生重要转变。人们认识到，仅仅依赖于添加了大量保守假设的“保守估算法”进行安全分析，虽然能够确保结果“安全”，但可能掩盖真实的物理过程，不利于深入理解事故演变和优化安全措施。因此，“最佳估算加不确定性”（Best-Estimate Plus Uncertainty, BEPU）的方法论应运而生并逐渐成为主流 。

BEPU方法的核心思想是：

1.最佳估算（Best-Estimate）： 使用最接近物理真实的模型和参数来模拟事故过程，而不是人为引入保守假设。RELAP5和TRAC等代码正是实现最佳估算分析的工具。

2.不确定性量化（Uncertainty Quantification, UQ）： 系统性地评估和量化由于模型自身的简化、输入参数的不确定性、经验关系式的误差等所有不确定性来源，对最终计算结果（如PCT）造成的影响范围。

这种转变推动了热工水力模型验证工作的系统化和科学化。大量的独立效应实验（SET）和积分效应实验（IET），如LOFT、Semiscale、BETHSY等，被用于验证和改进代码中的物理模型和整体性能 。同时，统计学方法被引入不确定性分析，使得安全分析结果能够以概率分布的形式呈现，为风险指引决策提供了更科学的依据。美国NRC也逐步修订法规，允许在满足严格验证和不确定性分析要求的前提下，使用BEPU方法进行许可申请 。

2.4 多维与多物理场耦合的深化 (2000s - 2010s)

进入21世纪，随着计算机性能的爆炸式增长，热工水力模型开始向更高维度和更广阔的物理领域拓展 。

•从一维到三维： 传统系统代码主要基于一维模型，难以精确描述压力容器下降段、堆芯内部等区域的复杂三维流动现象（如硼稀释、压力热冲击下的分层流）。因此，在系统代码中集成三维模块（如RELAP5-3D, TRACE中的三维容器模块）成为重要发展方向 。同时，独立的计算流体动力学（CFD）代码（如ANSYS-CFX, FLUENT）开始被用于对关键部件进行精细化三维模拟 。

•多物理场耦合： 反应堆内的物理过程是紧密耦合的。热工水力状态的变化会影响中子通量分布（通过慢化剂密度和温度反馈），而中子通量又决定了热量产生。因此，将热工水力代码与中子动力学代码（如PARCS）进行耦合，实现更真实的堆芯动态模拟，成为研究热点 。此外，与燃料性能分析代码、结构力学代码的耦合也得到越来越多的关注。

这一时期，美国NRC启动了代码现代化计划，旨在整合RELAP5、TRAC等多个传统代码的优点，开发新一代的系统代码TRACE（TRAC/RELAP Advanced Computational Engine），它拥有更现代化的软件架构，并内置了三维中子动力学模块PARCS 。

3.面向先进反应堆与高性能计算的当代发展 (2020s - 至今)

当前，热工水力模型的发展呈现出两大趋势：一是面向第四代核反应堆等先进堆型的特殊需求，二是充分利用现代高性能计算（HPC）平台的能力。

•先进反应堆建模： 熔盐堆（MSR）、钠冷快堆（SFR）、高温气冷堆（HTGR）等先进堆型采用了与传统轻水堆截然不同的冷却剂和设计理念，对热工水力模型提出了全新挑战（详见第五章）。例如，MSR中的燃料与冷却剂混合，需要模拟裂变产物输运；SFR的液态金属冷却剂具有极高的热导率；HTGR则涉及复杂的气体动力学和辐射传热。为此，一系列新的专用代码或在现有代码基础上扩展的模块正在被开发，如基于MOOSE框架的SAM和Pronghorn 。

•高保真度模拟： 随着HPC算力的提升，直接数值模拟（DNS）和大规模涡模拟（LES）等高保真度CFD方法开始被用于研究湍流、两相流界面等基础物理现象，其结果可用于开发或验证传统模型中的经验关系式 。这种从实验数据驱动向第一性原理模拟驱动的转变，预示着热工水力模型未来的发展方向。

总的来说，热工水力模型的历史是一部不断追求更高真实性、更广适用性和更强预测能力的历史，每一次重大进步都与安全需求的提升、计算机技术的发展和对物理现象认识的深化密不可分。

第三章：热工水力模型的具体实现方式与关键技术

热工水力模型的实现是一个多层次、多尺度的复杂体系。根据建模的精细程度和应用范围，可以大致分为三个层级：系统分析代码、子通道分析代码和计算流体动力学（CFD）代码 。这三个层级相辅相成，共同构成了核反应堆热工水力分析的完整工具箱。

3.1 建模方法的分类

•系统分析代码（System Analysis Codes）：

￮特点： 这是应用最广泛的一类代码。它们将整个反应堆冷却剂系统（一回路甚至二回路）离散为一系列由一维“组件”（如管道、堆芯、泵、阀门）和零维“节点”（连接点）构成的网络 。其目标是模拟整个系统的宏观、整体行为，特别是在各种瞬态和事故工况下的响应。

￮优势： 计算速度相对较快，能够进行数小时甚至数天的事故过程模拟，是安全分析和许可申请的主力工具。

￮局限性： 基于一维或简化的多维假设，无法精确描述组件内部的详细三维流动和传热分布。其物理模型大量依赖于从实验数据中总结出的经验关系式（如壁面摩擦、相间传热、CHF等），这些关系式的适用范围有限。

￮代表代码： RELAP5, TRACE, CATHARE, ATHLET, MARS等 。

•子通道分析代码（Subchannel Analysis Codes）：

￮特点： 这类代码专注于反应堆堆芯区域的详细热工水力分析。它将燃料组件内的冷却剂流道划分为一系列相互连接的平行“子通道”，并求解子通道间的横向质量、动量和能量交换 。

￮优势： 相比系统代码，它能提供更精细的堆芯内部温度和空泡份额分布信息，是进行热设计、计算偏离泡核沸腾比（DNBR）或临界功率比（CPR）等热工裕度的关键工具。

￮局限性： 仍属于简化模型，横向动量方程通常经过高度简化，且同样依赖经验关系式。其计算范围通常局限于堆芯或单个燃料组件。

￮代表代码： COBRA系列（如COBRA-TF）, FLICA, DYN3D/CRONOS等 。

•计算流体动力学（CFD）代码（Computational Fluid Dynamics Codes）：

￮特点： CFD代码通过在三维空间中求解完整的Navier-Stokes方程组，来模拟流体的详细流动和传热行为。它可以捕捉到湍流、涡流、二次流、热分层等复杂的局部三维现象。

￮优势： 提供最高保真度的流动细节，能够减少对经验关系式的依赖，特别适用于分析系统代码和子通道代码无法处理的复杂三维流动问题，如压力容器下降段的冷却剂混合、燃料组件入口的流量分布等。

￮局限性： 计算成本极其高昂，通常只能用于模拟小尺度、短时间的稳态或瞬态过程。其在两相流模拟方面的准确性和可靠性仍是当前研究的重点和难点。商业CFD代码在核工领域的特定流型下验证不足，是其应用的一大挑战 。

￮代表代码： ANSYS CFX/Fluent, STAR-CCM+, 以及专门为核领域开发的NEPTUNE_CFD等 。

这三类代码形成了一个金字塔结构：系统代码位于塔顶，覆盖范围广但精度较低；CFD位于塔底，精度高但适用范围窄；子通道代码居中。在实际应用中，常常采用多尺度耦合策略，例如，用CFD的精细计算结果来标定或验证系统代码中的局部模型参数。

3.2 主流系统分析代码详解

系统分析代码是核安全分析的基石，以下对几个国际主流代码进行详细介绍。

3.2.1 RELAP5系列

•历史与发展： RELAP (Reactor Excursion and Leakage Analysis Program) 系列代码是核工业中最著名、应用最广泛的系统代码之一。其历史可追溯至1966年的RELAP1 。经过不断发展，RELAP5 成为该系列的第五代主力代码。其第一个版本RELAP5/MOD0于1979年完成开发，主要用于模拟PWR的LOCA 。后续版本（MOD1, MOD2, MOD3）不断完善物理模型，扩展应用范围至BWR和非能动系统。

•主要版本与特点：

￮RELAP5/MOD3.3： 这是NRC官方维护的最后一个主要版本，发布于2010年左右 。它是一个非常成熟和经过广泛验证的代码。

￮RELAP5-3D©： 在MOD3.3的基础上，由Innovativ Systems Software (ISS) 公司商业化开发，并由国际RELAP5用户组（IRUG）监督。其最显著的特点是增加了三维流体动力学模块，可以对反应堆压力容器等关键部件进行三维模拟 。最新版本（如Ver. 4.4.2）于2018年发布 。

•维护与许可： NRC已停止对RELAP5的积极开发，其维护主要通过代码评估与维护计划（CAMP）进行Bug修复 。而RELAP5-3D©作为一个商业软件，由IRUG负责维护和改进，并通过会员制进行许可 。

•应用： 广泛应用于全球范围内的轻水堆瞬态与事故分析，包括LOCA、丧失厂外电（SBO）、预期瞬态无停堆（ATWS）等 。它也被扩展应用于钠冷快堆（SFR）和高温气冷堆（HTGR）等先进堆型的初步分析 。

3.2.2 TRAC/TRACE系列

•历史与演进： TRAC (Transient Reactor Analysis Code) 系列由洛斯阿拉莫斯国家实验室于1975年开始为NRC开发，旨在成为一个“最佳估算”的系统代码 。它最初分为TRAC-PWR和TRAC-BWR两个版本。

•从TRAC到TRACE： 20世纪90年代末，NRC启动了代码现代化计划，旨在将旗下的RELAP5、TRAC-P、TRAC-B和RAMONA等多个代码整合到一个统一、现代化的软件平台中，这个平台的产物就是 TRACE (TRAC/RELAP Advanced Computational Engine) 。TRACE继承了其前辈代码的物理模型，并采用了更先进的软件架构和数值解法。

•主要版本与特点： TRACE V5.0是其一个重要版本，于2007年左右发布 。TRACE的核心优势在于其高度的集成性和先进性：

￮统一平台： 可处理PWR、BWR以及未来先进反应堆的分析。

￮多维能力： 内置一维和三维热工水力模块，以及三维中子动力学模块。

￮多物理场耦合： 与中子动力学代码PARCS (Purdue Advanced Reactor Core Simulator) 紧密耦合，能够进行精细的耦合堆芯动力学/热工水力分析 。

•维护与许可： TRACE由美国NRC作为其官方的、主要的热工水力安全分析代码进行开发和维护，主要通过CAMP项目进行国际合作与验证 。

•应用： TRACE被NRC用于新型反应堆（如AP1000, ESBWR）的许可审评和安全研究，被认为是未来核电厂热工水力安全分析的首选工具之一 。

3.2.3 CATHARE

•背景： CATHARE (Code for Analysis of THermal-hydraulics during an Accident of Reactor and safety Evaluation) 是由法国原子能和替代能源委员会（CEA）、法国电力公司（EDF）、法国核安全与辐射防护研究院（IRSN）以及法马通（Framatome）联合开发的系统代码 。它是欧洲最主流的热工水力系统代码，相当于法国版的RELAP5/TRACE。

•主要版本与特点： CATHARE-2是其广为人知的版本 。该代码基于两流体六方程模型，经过了大量法国和国际实验的验证。CATHARE-3是其最新一代版本，拥有更灵活的架构，支持与CFD、燃料性能等其他代码的耦合。

•应用： 主要用于法国及其他采用法国核电技术的国家的PWR和VVER反应堆的热工水力事故分析 。它也被扩展用于第四代反应堆的计算 。

3.3 子通道分析代码

3.3.1 COBRA-TF

•背景： COBRA-TF (Coolant Boiling in Rod Arrays - Two Fluid) 是COBRA系列代码的现代版本，是一款先进的子通道分析代码 。它由美国NRC资助，由宾夕法尼亚州立大学（原北卡罗来纳州立大学）和太平洋西北国家实验室（PNNL）等机构开发和维护 。

•特点： 与其他子通道代码不同，COBRA-TF采用了完整的三维、两流体、三场（液、气、液滴）模型，使其能够更真实地模拟沸腾、再淹没等复杂的两相流现象。

•应用： 主要用于轻水堆堆芯的详细热工水力分析，特别是在LOCA再淹没阶段的传热计算。由于其模型的先进性，它也常被用作开发和验证其他代码中物理模型的基准工具。

3.4 计算流体动力学 (CFD) 的应用

CFD在核工程中的角色越来越重要，它填补了系统代码和子通道代码在精细化模拟方面的空白。

•NEPTUNE平台： 由法国EDF和CEA共同开发，是一个专为核能两相流模拟设计的多尺度软件平台 。它包括系统代码（NEPTUNE/CATHARE）、子通道代码、以及核心的CFD代码NEPTUNE_CFD。这种一体化平台有助于在不同尺度模型之间传递信息，实现更一致的分析。

•商业CFD软件： ANSYS CFX, Fluent等通用商业CFD软件也广泛应用于核工业 。它们通常用于分析稳态或准稳态下的局部三维问题，如压力容器内的热混合、燃料组件的压力损失特性、非能动系统的自然循环流动等。然而，将这些通用软件应用于核安全分析时，其在特定两相流工况下的模型验证和确认是一个必须严肃对待的挑战 。

3.5 新一代与先进堆型专用代码

为应对第四代反应堆的独特挑战，基于现代软件框架的新一代代码正在涌现。

•SAM 和 Pronghorn： 这两款代码均基于爱达荷国家实验室开发的开源多物理场对象导向模拟环境（MOOSE）构建 。

￮SAM (System Analysis Module)： 作为一款先进的系统分析代码，它专为非轻水堆（如SFR, LFR, MSR）设计，采用了更灵活的组件化架构和适用于液态金属、熔盐等特殊流体的物性模型。

￮Pronghorn： 一款专注于多孔介质流动和传热模拟的应用，特别适用于高温气冷堆的球床或棱柱块堆芯分析。

￮优势： MOOSE框架使得SAM和Pronghorn等应用能够便捷地实现与其他物理模块（如中子学、燃料性能、结构力学）的紧密耦合，极大地促进了高保真度、多物理场耦合模拟的发展 。

3.6 多物理场耦合技术

热工水力模型并非孤立存在，与其它物理过程的耦合是实现高保真度模拟的关键。

•中子动力学/热工水力耦合： 这是最重要的一种耦合。热工水力计算出的冷却剂密度和温度分布，会作为输入传递给中子动力学代码，影响中子慢化和吸收，从而改变功率分布。反过来，新的功率分布又会改变热工水力状态。这种反馈机制在反应性引入事故等快速瞬态中至关重要。TRACE/PARCS的耦合就是这种应用的典范 。

•热工水力/燃料性能耦合： 燃料芯块的温度、肿胀、裂变气体释放等行为受热工水力条件影响，同时又反作用于燃料到包壳的传热。

•热工水力/结构力学耦合： 流体压力和温度是结构应力分析的载荷输入。反之，结构的变形（如流动引起的振动）也可能影响流道和传热。

这些实现方式和技术共同构成了一个从宏观到微观、从单一物理场到多物理场耦合的复杂技术生态，支撑着现代核反应堆的设计、运行和安全评估。

第四章：核心争议点、各方立场与监管框架

热工水力模型领域并非铁板一块，其发展和应用始终伴随着深刻的技术路线之争和方法论辩论。这些争议的核心，是如何在计算成本、模型精度和安全保证之间取得最佳平衡。

4.1 核心争议一：保守估算法 vs. 最佳估算加不确定性 (BEPU)

这是热工水力安全分析领域最根本、影响最深远的方法论之争。

•4.1.1 保守估算法 (Conservative Approach) 的立场与依据

￮核心思想： 为了确保安全分析结果绝对包络所有不确定性，在建模时人为地、有选择地引入一系列保守的假设和参数。例如，在模拟LOCA时，假设热源功率最大、冷却系统性能最差、传热模型最不利等。

￮法规基础： 这一方法的法律基石是美国NRC法规 10 CFR 50.46 附录K (Appendix K) 。该附录在TMI事故前制定，详细规定了ECCS评估模型中必须采用的一系列保守假设。

￮立场与优势： 监管机构和部分保守的工业界人士认为，这种方法简单、明确，无需进行复杂的不确定性分析，只要计算结果满足接受准则（如PCT < 2200°F），就能提供一个确定的、有巨大安全裕度的“保证”。对于监管者而言，审查过程相对标准化。

￮劣势： 保守估算得到的结果往往与真实的物理过程相去甚远，可能导致对安全系统性能的误判。更重要的是，巨大的保守性裕度限制了核电站的运行灵活性和经济性，例如，它可能导致不必要的功率限制或过早的燃料更换。

•4.1.2 BEPU 方法 (Best-Estimate Plus Uncertainty) 的立场与优势

￮核心思想： 如第二章所述，BEPU方法使用最接近物理真实情况的模型进行计算，然后系统、科学地量化所有不确定性来源对计算结果的影响，最终给出一个带有置信区间的预测结果（例如，“有95%的置信度，PCT不会超过1800°F”）。

￮立场与优势：

▪工业界： 是BEPU方法的主要推动者。通过去除不必要的保守性，BEPU可以释放运行裕度，带来显著的经济效益，如功率提升、燃料循环延长、运行灵活性增加等。同时，它能提供对事故过程更深刻、更真实的理解，有助于优化事故管理策略。

▪学术界与研究机构： 积极发展BEPU所需的先进模型和不确定性量化（UQ）方法。他们认为BEPU代表了科学认知和工程实践的进步。

￮劣势： BEPU方法的实施技术门槛极高。它要求开发和验证高质量的最佳估算代码，并掌握复杂的统计学方法进行UQ分析。整个过程复杂、耗时，对计算资源和人员技能的要求都远高于保守估算法。

•4.1.3 监管机构的立场演变

￮美国NRC作为全球核监管的风向标，其立场经历了从强制保守到逐步接纳BEPU的转变。NRC认识到BEPU的科学性和经济价值，于1988年修订了10 CFR 50.46，允许在满足严格要求的情况下使用BEPU方法进行ECCS分析。

￮为了指导BEPU方法的实施，NRC发布了监管指南RG 1.157（后被RG 1.203替代），详细阐述了最佳估算代码的验证要求和可接受的UQ方法。

￮当前立场： NRC及全球大多数国家的监管机构都接受BEPU方法，但对其审评极为严格。申请者必须提供详尽的证据，证明其代码经过了充分的验证、UQ方法是合理的、且最终的不确定性区间足够包络真实情况 。保守估算法（附录K方法）依然是合规的选项之一，对于一些小型变更或不愿投入巨大资源进行BEPU分析的运营商来说，仍有其应用价值 。

4.2 核心争议二：模型的验证、确认与不确定性量化 (VVUQ)

如果说BEPU是目标，那么VVUQ（Verification, Validation, and Uncertainty Quantification）就是通往目标的必经之路，也是争议和挑战最集中的领域。

•4.2.1 VVUQ的挑战

￮Verification (验证)： 旨在确认“我们是否正确地求解了方程？”。这主要是一个数学和计算机科学问题，检查代码中是否存在Bug，数值解法是否收敛和准确。这相对直接，但对于复杂的耦合代码，依然充满挑战 。

￮Validation (确认)： 旨在确认“我们是否求解了正确的方程？”。这是将模型预测结果与实验数据进行比较，以评估模型描述物理现实的能力。这是VVUQ中最困难、最具争议的部分 。主要挑战包括：

▪实验数据稀缺： 获得高质量、适用于模型确认的实验数据成本高昂。特别是对于先进反应堆的新现象或事故工况下的极端条件，相关实验数据极其匮乏。

▪尺度转换问题： 实验通常在缩比设施中进行，如何将模型在缩比实验上的确认结果外推到全尺寸的真实反应堆，即“尺度效应”，是一个巨大的科学难题 。

▪现象的复杂性： 两相流等现象内在的随机性和复杂性，使得“完美”匹配实验数据几乎不可能。如何定义可接受的匹配程度，本身就是一个争议点 。

￮Uncertainty Quantification (UQ)： 量化所有已知和未知的不确定性。挑战在于如何全面地识别所有不确定性来源（模型参数、经验关系式、数值误差、输入条件等），以及如何评估“模型缺陷”这种难以量化的不确定性 。

•4.2.2 各方立场

￮工业界： 希望监管机构能提供清晰、稳定且高效的VVUQ路径和接受标准。他们倾向于利用现有的、经过广泛验证的国际标准问题（ISP）和基准算例来完成大部分确认工作，以降低许可成本和时间。

￮监管机构： 坚持审慎原则，要求申请者提供“足够”的确认证据，尤其是在模型应用于新设计或超出已验证范围的工况时。对于“足够”的定义，往往成为双方博弈的焦点。NRC等机构大力支持通过国际合作（如CAMP项目）来共享验证数据和经验，提高VVUQ的效率和一致性 。

￮学术界与研究机构： 致力于开发更先进的VVUQ方法论，如基于贝叶斯统计的UQ方法、多保真度模型融合等。他们强调基础物理实验的重要性，并推动建设开放、高质量的实验数据库。OECD/NEA下的核安全委员会（CSNI）及其工作组（如WGAMA）在组织国际合作、开展基准测试项目（Benchmarking）和制定最佳实践指南方面扮演着核心角色 。

•4.2.3 国际合作与基准测试项目

￮为了应对VVUQ的挑战，国际合作至关重要。OECD/NEA/CSNI等组织发起了大量的国际标准问题（ISP）和基准测试项目 。这些项目通常基于一个特定的实验（如BETHSY的SBLOCA实验）或一个明确定义的虚拟问题，邀请全球各地的机构使用自己的代码进行“盲算”或“开卷”计算。

￮通过比较不同代码的计算结果以及与实验数据的差异，可以评估代码的优劣、发现模型中的缺陷、量化代码间的不确定性。例如，HTGR热工水力基准、CANDU热工水力盲算基准等活动，为特定堆型的模型确认提供了宝贵的数据和经验 。

￮然而，一个突出的问题是，尽管存在大量基准活动，但目前缺乏一个公开的、统一的、量化的性能指标体系来系统比较各大主流代码在不同堆型、不同工况下的准确率、计算时间以及在现代HPC平台上的扩展性 。这使得用户在选择和评估代码时，仍需依赖零散的报告和各自的经验。

4.3 核心争议三：不同尺度模型的应用与耦合

•系统代码与CFD的鸿沟： 系统代码快但不准（在局部），CFD准但慢（且范围小）。如何在两者之间取长补短，是当前的一个技术热点和争议点。

￮争议点： 何时需要用CFD？CFD的结果如何可靠地反馈给系统代码？如何保证耦合过程中的数值稳定性和物理一致性？

￮立场： 一些人主张尽可能用高保真CFD替代系统代码中的经验模型，实现“由下至上”的建模。另一些人则认为，考虑到计算成本和CFD自身的不确定性，更现实的策略是“由上至下”，即用系统代码进行全局分析，仅在识别出的关键局部区域（如热点区域）进行CFD精细分析。

•多尺度建模策略： 理想的策略是根据问题的需要，无缝地在不同尺度的模型间切换或耦合。例如，用CFD模拟下降段的硼混合，将其出口的浓度分布作为边界条件输入到子通道代码，再将子通道代码计算的堆芯出口参数反馈给系统代码。实现这种无缝耦合的技术仍在发展中。

4.监管框架与关键文件解读

监管框架为上述所有争议设定了边界和规则。

•美国NRC体系：

￮10 CFR 50.46: 规定了ECCS的接受准则（PCT、氧化层厚度等），是所有LOCA分析的最终审判标准 。

￮NUREG系列报告: 这是NRC发布的技术报告系列，内容包罗万象。例如， NUREG-0800 (Standard Review Plan, SRP) 提供了NRC审评员审查许可申请的具体指南，其中章节15详细说明了对瞬态和事故分析的要求 。其他NUREG报告则可能发布关于某个特定代码的确认研究或UQ方法论的综述 。

￮Regulatory Guides (RG): 监管指南提供了NRC认可的、用于满足法规要求的方法。如前述的RG 1.203指导了BEPU方法的应用 ；RG 1.174则阐述了如何使用概率风险评估（PRA）进行风险指引决策，虽然不直接针对热工水力模型，但BEPU分析的结果是PRA的重要输入 。

•国际组织的角色：

￮OECD/NEA CSNI: 作为核安全研究的国际合作中心，其发布的报告和指南虽不具法律强制力，但在全球范围内拥有极高的技术权威。CSNI代码验证矩阵 (Code Validation Matrix) 为各大代码的确认工作提供了路线图和数据集。其下属的事故分析与管理工作组 (WGAMA) 专注于协调事故分析方法、代码开发和确认方面的国际合作。

￮IAEA (国际原子能机构): IAEA通过其安全标准系列（Safety Standards Series）为成员国提供核安全法规的框架和指导。IAEA TECDOCs 是一系列技术文件，常常发布关于特定技术主题（如热工水力建模、UQ方法）的会议纪要、研究报告和最佳实践指南 。

总而言之，热工水力模型领域正处于一个由争议和合作共同推动的动态发展阶段。BEPU方法论的胜利带来了VVUQ的巨大挑战，而监管框架的演进则不断为技术发展设定新的目标和边界。

第五章：不同堆型的热工水力模型特点与挑战

热工水力模型并非一成不变，其具体形式和关注重点与反应堆的类型密切相关。冷却剂的物性、堆芯的结构、能量转换的方式，共同决定了建模的特殊性和挑战性。

5.1 压水堆 (PWR) 与沸水堆 (BWR)

作为全球核电的主力堆型，PWR和BWR的热工水力模型最为成熟，拥有最丰富的实验数据和运行经验支持 。

•共同特点： 都使用水作为冷却剂和慢化剂，因此两相流是其热工水力建模的核心和难点。沸腾、气泡输运、相间滑移、干涸等现象的精确模拟至关重要。RELAP5, TRACE, CATHARE等主流系统代码最初都是为轻水堆开发的 。

•PWR的特殊性：

￮一回路在高压下运行，正常工况下冷却剂为单相液态水。沸腾只可能发生在某些事故工况或堆芯局部热点区域。

￮关注点：偏离泡核沸腾（DNB）是主要的热工限制，需要通过子通道代码精确计算DNBR。LOCA期间的再淹没、压力容器的压力热冲击（PTS）下的热分层和混合现象是安全分析的重点。

•BWR的特殊性：

￮冷却剂在堆芯内直接沸腾产生蒸汽，驱动汽轮机。堆芯内存在大范围的气液两相流，空泡份额极高。

￮关注点：临界功率（干涸）是主要的热工限制，需要计算CPR。两相流不稳定性（密度波振荡）可能影响反应堆的稳定运行，是BWR特有的安全问题，需要通过耦合的热工水力/中子动力学模型进行分析 。

5.2 高温气冷堆 (HTGR)

HTGR采用氦气作为冷却剂，石墨作为慢化剂和结构材料，出口温度可达750-950℃，具有优异的非能动安全性 。

•模型特点与挑战：

￮气体动力学： 冷却剂为单相氦气，不存在相变问题，模型相对简单。但氦气密度低、热容小，需要很高的流速和巨大的流量来导出热量，因此压降计算和风机/压缩机建模非常重要。

￮高温传热： 在极高温度下，热辐射成为与对流同等重要的传热方式，尤其是在堆芯和容器空腔内。模型必须包含精确的表面发射率模型和视角因子计算，这对于传统轻水堆代码是一个新挑战。

￮堆芯结构： 存在两种主要堆芯形式——球床堆（Pebble Bed）和棱柱块堆（Prismatic Block）。两者都可被视为多孔介质。因此，需要采用多孔介质模型来描述冷却剂在燃料球或石墨块冷却通道中的流动和传热。Pronghorn等基于多孔介质理论的代码在此领域有良好应用前景 。

￮事故工况： 在失压失流事故（DLOFC）下，堆芯依靠石墨的巨大热容和向环境的自然散热（导热、辐射、自然对流）来维持低温，这是其非能动安全性的核心。热工水力模型必须能准确模拟这些长时间、低流速下的多模式传热过程。

5.3 熔盐堆 (MSR)

MSR以熔融的氟盐或氯盐作为冷却剂，其中一些设计中燃料直接溶解在盐中，形成了液态燃料 。

•模型特点与挑战：

￮复杂流体物性： 熔盐是高温、高粘度的非牛顿流体，其物性（粘度、导热系数等）随温度和组分（特别是裂变产物的积累）动态变化。模型必须能处理这些复杂的、非恒定的物性。

￮燃料与裂变产物输运： 在液态燃料堆中，热源本身在流动。模型不仅要计算盐的流动和传热，还必须耦合求解缓发中子先驱核、气态和固态裂变产物在整个一回路中的输运、沉积和析出过程。这对代码架构提出了极高要求。

￮化学腐蚀： 熔盐对结构材料有腐蚀作用，腐蚀产物会进入盐中。这种化学过程与热工水力条件（温度、流速）相互影响，需要多学科耦合模型来描述。

￮氚的产生与输运： 熔盐中的锂会产生氚，氚会渗透到二回路甚至环境中。模拟氚在系统中的产生、输运、渗透和捕获是MSR许可申请中的一个关键安全问题。SAM等新一代代码正致力于解决这些非轻水堆特有的挑战。

5.4 钠冷快堆 (SFR)

SFR使用液态金属钠作为冷却剂，采用快中子谱，可实现燃料增殖 。

•模型特点与挑战：

￮液态金属传热： 钠具有极高的导热系数和极低的普朗特数（Prandtl Number）。这意味着其动量边界层远厚于热边界层，传热主要由导热主导，对流作用相对较弱。这使得传统用于水的传热关系式完全不适用，需要专门针对液态金属的传热模型。

￮不透明性： 液态钠完全不透明，使得辐射传热仅在自由液面发生。

￮化学活性： 钠与空气和水会发生剧烈化学反应。因此，必须对涉及钠泄漏的事故进行专门分析，这可能需要耦合化学反应动力学模型。

￮系统设计： SFR通常采用池式设计，一回路所有设备都浸泡在巨大的钠池中。这带来了优异的热惯性和自然循环能力，但同时也使得池内的热分层、射流冲击等三维热工水力现象变得非常复杂和重要，往往需要CFD进行精细分析。

5.5 其他先进堆型

•铅冷快堆 (LFR): 与SFR类似，但铅的密度更大、腐蚀性更强。

•超临界水冷堆 (SCWR): 冷却剂处于超临界状态，不存在相变，但物性在拟临界点附近剧烈变化，对传热模型提出严峻挑战。

总之，从轻水堆到第四代反应堆，热工水力模型的复杂性和挑战性不断增加，驱动着模型理论、数值方法和软件工具的持续创新。

第六章：未来发展方向与展望

热工水力模型正站在一个由数据科学、高性能计算和先进传感技术共同定义的新十字路口。未来的发展将朝着更高保真度、更强预测能力和更深度的智能化融合方向前进。

6.1 高保真度建模 (High-Fidelity Modeling)

•从经验到第一性原理： 未来的趋势是用更高保真度的模拟方法，如直接数值模拟（DNS）和大规模涡模拟（LES），来替代传统模型中高度依赖经验的关系式 。例如，通过DNS模拟沸腾过程中气泡的生成、合并与破碎，可以获得对相间传热和动量交换更基本的物理认识，从而构建更普适、更精确的宏观模型。虽然目前DNS/LES的计算成本极高，但随着HPC算力的增长，它们将成为模型开发而非仅仅是基础研究的工具。

•全堆芯CFD模拟： 随着千万亿次甚至百亿亿次计算机的出现，实现对整个反应堆压力容器甚至整个一回路的瞬态CFD模拟正从梦想变为可能。这将消除传统系统代码中大量的简化假设，提供前所未有的时空分辨率，从而发现未知的安全现象。

6.2 人工智能与机器学习的应用 (Application of AI/ML)

AI/ML技术正在为传统的热工水力建模领域注入新的活力，其应用潜力巨大。

•代理模型/降阶模型 (Surrogate/Reduced-Order Models)： 复杂的BEPU分析需要运行数千次数万次的热工水力代码，计算成本高昂。可以利用机器学习算法（如神经网络、高斯过程回归）训练一个代理模型，使其能够以极快的速度（毫秒级）精确复现昂贵代码的计算结果。这将彻底改变UQ分析、实时优化和概率安全评估的效率。

•数据驱动的物理模型发现： 利用先进的机器学习技术，可以直接从高保真模拟数据或海量实验数据中“发现”新的、更优的物理模型或经验关系式，克服人类认知和建模能力的局限。

•智能数据融合： 将不同来源、不同保真度的数据（实验数据、低保真模拟、高保真模拟）通过AI方法进行融合，以更低的成本获得对物理现象更全面的认识。

6.3 数字孪生与实时监控 (Digital Twins and Real-time Monitoring)

•概念： 数字孪生是物理核电站的实时、高保真度虚拟映射。它将高精度的热工水力模型与从电站传感器获取的实时运行数据相结合，通过数据同化技术不断校准模型，使其始终与物理实体保持同步。

•应用：

￮预测性维护： 通过模拟未来的运行状态，预测设备的老化和潜在故障，实现从“故障后维修”到“预测性维护”的转变。

￮运行优化： 实时模拟不同操作策略的效果，找到最优的运行参数，提升电站效率和经济性。

￮应急响应增强： 在事故发生时，数字孪生可以比实时快数百倍的速度预测事故的未来演变，为操作员和决策者提供宝贵的决策窗口和行动建议。

6.4 实验数据库的扩展与共享

无论是传统模型的确认，还是AI模型的训练，高质量的实验数据都是不可或缺的“燃料”。

•面向先进反应堆的实验： 未来需要大力投入建设针对MSR、SFR、HTGR等先进堆型中特有热工水力现象的实验设施，获取关键的确认数据 。

•高分辨率测量技术： PIV（粒子图像测速）、LIF（激光诱导荧光）等先进测量技术的发展，能够提供前所未有的高时空分辨率的流场和温度场数据，为高保真度CFD的确认提供了可能。

•开放数据倡议： 建立全球共享、格式标准、质量受控的核热工水力实验数据库，将极大地促进模型开发和确认的效率，降低整个行业的研发成本。

6.5 监管科学的现代化 (Modernization of Regulatory Science)

新技术的出现也对传统的核安全监管模式提出了挑战。

•如何审评AI模型？ 机器学习模型（特别是深度神经网络）通常被视为“黑箱”，其决策过程不透明。监管机构需要发展新的方法和标准，来审评和许可这些基于AI的模型的可靠性和安全性。

•如何适应高保真模拟？ 当模拟结果的精度可能超过实验测量精度时，如何定义“确认”？如何处理模拟揭示出的、现有安全准则未覆盖的新现象？这些都是监管科学需要回答的前沿问题。

•风险指引与实时决策： 数字孪生等技术将使实时、动态的风险评估成为可能。监管框架需要从当前基于静态、离线分析的模式，向能够适应动态风险信息、支持实时决策的模式演进。

第七章：结论

热工水力模型是贯穿核反应堆设计、建造、运行和退役全生命周期的核心技术，是保障核安全、提升核能经济性的科学基石。其发展史深刻地反映了核安全理念的演进、计算能力的飞跃以及对复杂物理现象认识的深化。

面向未来，热工水力模型正处在一个激动人心的变革时代。高性能计算正在催生更高保真度的模拟能力，人工智能和机器学习技术则有望通过代理模型、数据驱动建模等方式颠覆传统的分析范式。数字孪生技术预示着模型将从离线分析工具走向在线实时决策支持系统。然而，这些新技术的应用也对实验数据获取、模型可解释性以及监管科学的现代化提出了新的、深刻的挑战。

核技术论坛

阅读 赞 分享 言