[发明专利]基于自适应迭代步的自然循环系统特性多尺度模拟方法

说明书

本发明是一种基于自适应迭代步的自然循环系统特性多尺度模拟方法，其特点是，包括建立自然循环系统、设计交互参数自适应收敛判据、构建自适应收敛阈值函数。采用自然循环系统的驱动压头和总阻力压降间的定量关系更新交互流速；将一维和三维程序交互参数的方差作为判断三维程序计算参数是否达到稳定的收敛阈值；基于历史交互流速构建收敛阈值的自适应函数，实现耦合程序快速收敛至精确值。本发明具有自适应能力强、高效易行和不依赖经验等优点。特别适用于分析具有复杂结构热源的自然循环系统流动与换热特性。

技术领域

本发明涉及流动与换热数值计算领域，是一种基于自适应迭代步的自然循环系统特性多尺度模拟方法，特别适用于具有一维和三维流动特征的自然循环系统热工水力特性的定量分析。

背景技术

由于自然循环系统具有结构简单，固有安全性好等优点，自然循环式非能动安全系统在第三代核电厂中得到普遍应用。因此，研究基于自然循环式非能动安全系统的热工水力特性对核电厂的安全设计至关重要。

目前，本领域技术人员针对自然循环系统热工水力特性的研究主要采用基于集总参数法的一维分析程序，该程序具有计算流速快、参数调节灵活等优点。然而，对于具有复杂结构的热源如反应堆热源和列管式换热器，其内部存在显著的三维流动特征，如二次流、涡流等。因此，采用传统一维分析程序无法模拟其内部复杂的三维流场。同时，如果仅采用三维程序对整个自然循环系统模拟则需要庞大的网格量，时间成本极高。因此，本发明结合一维程序和三维程序软件的优点，开发了一维程序和三维程序的自适应耦合计算方法，实现了对具有一维和三维流动特征的自然循环系统热工水力特性的快速和准确模拟，可为研究人员提供极具价值的研究工具。

在现有的文献中，针对一维程序和三维程序自适应多尺度耦合的研究很少，其中张银星等提出一维用户程序和三维程序STAR-CCM+的多尺度耦合方法(张银星,高璞珍,何晓强,等.STAR-CCM+与一维用户程序耦合方法[J].哈尔滨工程大学学报,2020,41(11):1669-1674.)。该方法讲述了一维用户程序和三维程序STAR-CCM+的数据交互方法。通过STAR-CCM+预设固定迭代步数，使其计算参数稳定，然后与一维用户程序进行数据交互。然而，实际应用中发现，当固定迭代步数选取过大时，会导致耦合收敛的时间显著增加。当固定迭代步数选取过小时，交互参数未达到稳态或接近稳态，数据交互会导致耦合程序计算发散。因此，固定迭代步数的预设值依赖数值模拟经验，通常采用试算法求得。

综上分析可知，采用固定迭代步的一维和三维多尺度耦合方法具有自适应能力、普适性和鲁棒性较差的缺点。

发明内容

鉴于现有技术基于固定迭代步的一维和三维程序耦合方法的局限性，本发明的目的是，提出一种科学合理，自适应能力强，应用价值高的基于自适应迭代步的自然循环系统特性多尺度模拟方法，本发明的方法能够实现一维和三维耦合程序的快速、准确收敛，同时该方法也可用于计算分析含其它复杂热源的自然循环系统特性。

实现本发明目的采用的技术方案是，一种基于自适应迭代步的自然循环系统特性多尺度模拟方法，其特征是，它包括以下步骤：

1)建立自然循环系统的第一管路、第二管路、第三管路和第四管路的一维流动阻力与换热计算程序和热源的三维流动与换热计算程序；

2)设初速流速为v_i，并将v_i作为三维计算程序的入口流速边界条件，采用三维程序计算热源的压降ΔP和出口温度T_out；

3)计算三维程序得到的出口温度T_out和压降ΔP_re的收敛阈值δ，将三维程序计算的热源出口温度T_out和热源压降ΔP_re，根据连续n次迭代步的历史值求方差，并作为热源出口温度和压降的参数收敛判据，具体为且I为三维程序迭代步，n为间隔迭步数；