一种新型的管壳式换热器多尺度耦合协同换热模拟方法与流程

本发明涉及管壳式换热器设计计算技术领域，具体涉及一种新型的基于多孔介质模型的多尺度耦合求解模拟方法。

背景技术：

管壳式换热器广泛应用在核电、锅炉、电厂、石化等能源行业，管侧和壳侧分别通有各自的热源流动介质和传热流动介质，两侧的介质包含不同的流动行为和传热及相变行为。目前，针对管壳式换热器常用的设计计算方法主要有两种：一种是基于传热学数学模型的经验公式计算方法；另外一种是基于全几何模型的计算流体力学(cfd)分析的方法。基于传热学数学模型的经验公式计算方法是根据管壳式传热经验公式，进行管壳式换热器传热性能计算；基于全几何模型的计算流体力学方法(cfd)是利用设计的详细几何模型进行流体力学传热计算分析。基于传热学数学模型的经验公式计算方法是换热器设计分析的传统计算方法，该方法无法详细展示换热器各区域的换热特性和物理量分布数据，多采用经验假设，对设计人员的经验要求较高，精度不高，所以采用经验计算方法设计管壳式换热器是有较大的工程限制的；采用整体模型的cfd计算方法，该方法计算结果详实，但是整体模型计算效率较低，设计计算周期太长，特别是针对存在相变传热的复杂大模型，当前的计算机技术很难满足整体模型计算的需求，因此，该方法在大型换热器设计中也存在较大的局限性。为了提高管壳式换热器设计计算的数据详实性及计算的快速性等问题，本发明提供了新型的针对管壳式换热器设计计算的模拟方法。

技术实现要素：

针对于管壳式换热器设计中的经验公式计算方法换热参数不够详实及完整计算模型设计计算速度慢等问题，本发明提供了新型的管壳式换热器多尺度耦合协同换热模拟方法。

为了解决上述问题，本发明是通过以下技术方案实现的：

步骤1：基于传热关系式，建立管侧一维传热模拟自定义函数(udf，userdefinedfunction)，通过分配c_udm存储空间存储每个网格的管侧流体温度和壁面温度等重要变量。

步骤2：以进口条件分别初始化管侧/壳侧工质的温度和流动状态。初始化整个多孔介质计算区域，基于每个单元的初始值，对管侧/壳侧进行一维、三维耦合传热模拟；

步骤3：随着迭代计算的进行，管侧内部温度升高，存在沸腾传热的可能，基于管侧温度的变化情况判定各温度条件下所采用的传热关系式，若存在某个单元的温度值超过其临界温度值(tw>＝tonb)，过冷沸腾发生。此时相应调整该部分网格管侧换热模型，基于新模型计算该部分网格管侧/壳侧的换热量。随着传热的继续进行，管侧流体和传热管壁面温度也随之不断上升，在loop循环时可检测到某部分位置网格tf>＝ts，饱和沸腾换热区出现。此时相应调整该部分网格管侧换热模型，基于新模型计算该部分网格管侧/壳侧的换热量。

步骤4：依步骤3的方法不断迭代进行计算，不断根据判别条件更新局部网格的二次换热模型，直至最后所有分区段的出现。最终形成稳定的多区段换热模式，直至所有控制方程计算获得收敛，计算结束。

附图说明

图1单元网格求解示意图

图2传热模型选择示意图

图3管-壳侧耦合热量传递示意图

具体实施方式

为了解决管壳式换热器设计过程中的快速设计及设计准确性问题，结合图1-图3对本发明进行了详细说明，具体实施步骤如下：

计算步骤:

步骤1：对已建立的管壳式换热器模型，采用多孔介质方法对模型进行网格划分，根据划分的网格尺寸，如每个网格内只包含一根传热管单元，如图1，则管侧流体能量控制方程为：

as(vs,outρs,outcs,outts,out-vs,inρs,incs,ints,in)＝qap(1)

其中as为传热管内截面积，ap一次侧传热管表面面积。v、ρ、c、t分别为二次侧流体流速、密度、比热容、温度。in表示流进该网格控制体，out表示流出。以流进流出该单元网格温度的平均数作为该网格管侧流体温度，用以计算管-壳侧换热。可以获得每一个网格的管侧流体温度值，在udf中通过c_udmi(c，t，i)存储。

步骤2：以传热数学模型为基础，利用udf创建管侧的多孔介质区域管侧换热数学模型，利用c_udmi分配存储空间，用以存储每个网格的管侧流体温度、管道壁面温度、管侧流体含气率、管壳传热量等重要变量。

根据管侧含气量及温度的不同选定不同的传热关系式，具体传热关系式选择如图2所示，基于已有的传热相关数学模型描述，有以下经验公式：

(1)αg>＝0.99，管侧为单相蒸汽，管侧换热模型采用dittus-bolter关系。针对轴向高度z处的某一网格控制体，通过获取该网格壳侧流体温度，可计算出一次侧流体经过换热器入口流至此处时的焓差，该焓差△ht即是其向该网格包含的管侧流体所释放的热量。通过计算管侧流体从入口温度达到相应网格中管侧温度下汽相份额0.99时所需要的热量q，当△ht>＝q,即判定αg>＝0.99。传热系数为dittus-bolter关系

(2)tw>＝tmfb管侧工质为膜态沸腾。tmfb采用leidenfrost温度确定

管侧换热热流密度采用zuber和tribus提出的计算公式

(3)tw>＝tchf，恶化沸腾换热区。采用曹丹提出的公式：

nu＝0.017θrecm^0.8prcr^0.8

(4)tf>＝ts，饱和沸腾换热区。chen关系式

h＝hc+hncb

hc＝0.023f(gd(1-x)/μ1)^0.85·(μcp/k)1^0.4(d/dc)^0.1kl/d

(5)tw>＝tonb，过冷沸腾区。tonb的计算参考关系为tonb＝ts+0.0225·q^0.5e^-p/8.7，相应条件下的传热关系为修正的chen换热系数关系式

h＝hc+hncb(ti-ts)/(ti-tl)

hc＝0.023(gd/μ1)^0.85·(μcp/kl)^0.4(d/dc)^0.1kl/d

(6)在不满足以上五种条件时，其他情况下，判定管侧流动为单相液体，管侧换热关系采用dittus-bolter关系

以上分段中，涉及每个单元格处传热管内壁面温度的计算，由管壳侧的热流密度可以获得换热管内壁面温度：

步骤3：基于以上换热关系式，利用udf将管侧传热代码耦合到壳侧多孔介质模型中进行计算分析，具体实现过程如图3所示：

步骤3.1：初始化整个多孔介质模型，分别给定管侧及壳侧的初始温度场及流场，基于初始给定管侧各个单元的初始温度和壳侧各单元的初始温度，利用自定义函数进行传热模拟，得到下一步各单元的管侧及壳侧的温度信息，利用c_udmi进行数据储存。

步骤3.2：随着计算的进行，loop循环各单元的温度及含汽量，通过对温度值及含汽量判定采用哪种关系式对管-壳侧进行耦合热量传递，通过计算将计算结果储存到指定的储存单元c_udmi中；

步骤3.3：随着计算的继续进行，重复步骤3.2的工作，直到各单元对应的管-壳侧温度值计算收敛，得到换热器管-壳侧的温度、热量、流速等所有信息，计算结束。