(19)中华 人民共和国 国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202111516242.2 (22)申请日 2021.12.08 (71)申请人 东北电力大 学 地址 132012 吉林省吉林市船 营区长春路 169号 (72)发明人 孙建闯　蔡伟华　李智明　张文超　 李伟超　侯延栋　杜立鹏　 (74)专利代理 机构 吉林市达利专利事务所 22102 代理人 陈传林 (51)Int.Cl. G06F 30/20(2020.01) G06F 119/08(2020.01) (54)发明名称 基于自适应迭代步的自然循环系统特性多 尺度模拟方法 (57)摘要 本发明是一种基于自适应迭代步的自然循 环系统特性多尺度模拟方法， 其特点是， 包括建 立自然循环系统、 设计交互参数自适应收敛判 据、 构建自适应收敛阈值函数。 采用自然循环系 统的驱动压头和总阻力压降间的定量关系更新 交互流速； 将一维和三维程序交互参数的方差作 为判断三维程序计算参数是否达到稳定的收敛 阈值； 基于历史交互流速构建收敛阈值的自适应 函数， 实现耦合程序快速收敛至精确值。 本发明 具有自适应能力强、 高效易行和不依赖经验等优 点。 特别适用于分析具有复杂结构热源的自然循 环系统流动与换 热特性。 权利要求书1页 说明书4页 附图5页 CN 114329923 A 2022.04.12 CN 114329923 A 1.一种基于自适应迭代步的自然循环系统特性多尺度模拟方法， 其特征是， 它包括以 下步骤： 1)建立自然循环式非能动安全系统的第一管路、 第二管路、 第三管路和第 四管路的一 维流动阻力与换 热计算程序和热源的三维流动与换 热计算程序； 2)设初速流速为vi， 并将vi作为三维计算程序的入口流速边界条件， 采用三维程序计算 热源的压降ΔP和出口温度Tout； 3)计算三维程序得到的出口温度Tout和压降ΔPre的收敛阈值δ， 将 三维程序计算的热源 出口温度Tout和热源压降ΔPre， 根据连续n 次迭代步的历史值求方差， 并作为热源出口温度 和压降的参数收敛判据， 具体为 且 I为三维程序 迭代步， n 为间隔迭步数； 4)当出口温度Tout和压降ΔPre的方差均小于收敛阈值δ， 将三维程序计算的出口温度 Tout和压降ΔPre传递给一维计算 程序， 并计算驱动压 头ΔPd和总阻力压降ΔPz； 5)根据驱动压头ΔPd和总阻力压降ΔPz的定量关系更新流速vi， 得到更新后的流速 vi+1， 具体为： vi+1＝vi[1+(ΔPd‑ΔPz)/(ΔPd+ΔPz)]； 6)根据历史交互流速vi和vi+1更新步骤4)中的收敛阈值δ， 具体为δ＝arctan{[|vi+1‑vi |/max(vi+1,vi)]C}， 其中， C为常数， 取值为0.3， i为数据交互 次数， max()为取最大值函数， δ 根据一维和三维程序的交互流速对收敛阈值进行自适应更新， 并实现耦合程序的快速、 准 确收敛； 7)将更新后的流速vi+1传递给三维计算程序， 并转到步骤3)， 直至vi与vi+1相同， 完成一 维程序和三维程序的自适应多尺度耦合。权　利　要　求　书 1/1 页 2 CN 114329923 A 2基于自适应迭代步的 自然循环系统特性多尺度模拟方 法 技术领域 [0001]本发明涉及 流动与换热数值计算领域， 是一种基于自适应迭代步的自然循环系统 特性多尺度模拟方法， 特别适用于具有一维和三 维流动特征的自然循环系统热工水力特性 的定量分析。 背景技术 [0002]由于自然循环系统具有结构简单， 固有安全性好等优点， 自然循环式非能动安全 系统在第三代核电厂中得到普遍应用。 因此， 研究基于 自然循环式非能动安全系统的热工 水力特性对核电厂的安全设计至关重要。 [0003]目前， 本领域技术人员针对自然循环系统热工水力特性的研究主要采用基于集总 参数法的一维分析程序， 该程序具有计算流速快、 参数调节灵活等优点。 然而， 对于具有复 杂结构的热源如反应堆热源和列管式换热器， 其内部存在显著的三 维流动特征， 如二次流、 涡流等。 因此， 采用传统一维分析程序无法模拟其内部复杂的三维流场。 同时， 如果仅采用 三维程序对整个自然循环系统模拟则需要庞大的网格量， 时间成本 极高。 因此， 本发明结合 一维程序和三维程序软件的优点， 开发了一维程序和 三维程序的自适应耦合计算方法， 实 现了对具有一维和三 维流动特征的自然循环系统热工水力特性的快速和准确模拟， 可为研 究人员提供极具价 值的研究工具。 [0004]在现有的文献中， 针对一维程序和三维程序自适应多尺度耦合的研究很少， 其中 张银星等提出一维用户程序和三维程序STAR ‑CCM+的多尺度耦合方法(张银星,高璞珍,何 晓强,等.STAR ‑CCM+与一维用户程序耦合方法[J].哈尔滨工程大学学报,2020,41(11): 1669‑1674.)。 该方法讲述了一维用户程序和三维程序STAR ‑CCM+的数据交互方法。 通过 STAR‑CCM+预设固定迭代步数， 使其计算参数稳定， 然后与一维用户程序进行数据交互。 然 而， 实际应用中发现， 当固定迭代步数选取过大时， 会导致耦合收敛的时间显著增加。 当固 定迭代步数选取过小时， 交互参数未达到稳态或接近稳态， 数据 交互会导致耦合程序计算 发散。 因此， 固定迭代步数的预设值依赖数值模拟经验， 通常采用试算法求得。 [0005]综上分析可知， 采用固定迭代步的一维和三维多尺度耦合方法具有自适应能力、 普适性和鲁棒 性较差的缺 点。 发明内容 [0006]鉴于现有技术基于固定迭代步的一维和三维程序耦合方法的局限性， 本发明的目 的是， 提出一种科学合理， 自适应能力强， 应用价值高的基于自适应迭代步的自然循环系统 特性多尺度模拟方法， 本发明的方法能够实现一 维和三维耦合程序的快速、 准确收敛， 同时 该方法也可用于计算分析含其它复杂热源的自然循环系统特性。 [0007]实现本发明目的采用的技术方案是， 一种基于自适应迭代步的自然循环系统特性 多尺度模拟方法， 其特 征是， 它包括以下步骤： [0008]1)建立自然循环系统的第一管路、 第二管路、 第三管路和第四管路的一维流动阻说　明　书 1/4 页 3 CN 114329923 A 3