重力热管数学模型的建立及求解过程及程序设计

三章重力热管数学模型的建立及其模拟求解3.1重力式热管模型建立的分析思路对于重力式热管的凝结换热特性进行深入的分析。分析的基本思想如下:（1）对层流状态下的竖壁凝结换热建立模型，提出合理的假设条件；（2）取液膜微元控制体，对它进行受力分析和能量平衡分析，建立动量方程和能量方程，合理简化边界条件。简化方程，找出液膜厚度、切应力与液膜下降高度间的关系式；（3）通过能量平衡关系式，由己知的热管内压强（定性温度）和冷却介质带走的热热流，计算出热管的内壁温度；（4）离散液膜长度，假设在每一个微元液膜段内，液膜内的各个参数都是一样的（如果微元液膜足够小的话，这个假设是成立的），积分得到边界方程，并求得各个关系式；（5）在各个离散点处求出凝结换热系数，然后计算出所有离散点的凝结换热系数的算术平均值作为所要求的凝结换热系数的数值解。3.2物理模型建立和数学描述3.2.1模型合理的简化和假定图3-1重力式热管冷凝段模型热管内部的冷凝液膜和热管的内径比起来很薄，因此研究热管内部的凝结换热可以简化成研究大平板表面的凝结换热处理，这样可以使问题变得简单化。只考虑竖直方向，其他方向不考虑。在竖直的热管中，没有不凝结气体的影响，建立坐标系，它的流动模型如图3-1所示。在分析中，作若干合理的简化假定以忽略次要因素。除了已经明确的纯净饱和蒸汽层流液膜的假定外，还有:（1）蒸汽及凝结液的热物性是常数；（2）液膜的惯性力可以忽略（即控制方程中的对流项可以忽略不计）；（3）汽液界面上无温差，界面上液膜温度等于饱和温度tsat；（4）膜内的温度分布是线性的，即认为液膜内的热量转移只有导热，而无对流作用；（5）液膜表面平整无波动。（6）忽略不凝结气体的影响。3.2.2理论推导建立如图3-1所示的凝结液膜柱坐标系统，进行理论分析。应用N-S方程，列出液膜在竖直管内的连续性方程、动量方程和能量方程。（1）连续性方程0)(1rvrrxu（3.1）（2）轴向动量方程rurrrxuxgxprrvurxuu1)(1)(应用假设条件可以得到动量方程的简化形式：0)(1grurrrvl（3.2）动量方程的边界条件为0,uRrruRr,其中u为冷凝液膜的流动速度，μ为冷凝液膜的动力粘度。（3）能量方程rTrrrxTxrrvTcrxuTcpp1)(1)(应用假设条件可以得到能量方程的简化形式：01rTrrr（3.3）其对应的边界条件为wTTRr,vTTRr,（4）能量平衡方程图3-2控制体热量传递分析图在任意水平截面上，取高度为dx的单位宽度的液膜控制体，对该控制体进行能量平衡分析，如图所示，有000)()()(dxrdrTTciudxdrdrTTciudxrTRdrrTTciuimdvpfvpfRrvpfg通过初步处理可以得到dxrdrTTucdxdRdxTThmdvpwvfg0)()(（3.4）这就是控制体的能量守恒方程式。3.2.3对冷凝段模型的数学描述（1）冷凝段液膜速度及汽液界面处速度的数学描述方程。对动量方程（3.2）积分求解可以得到，液膜速度的分布式：RrRgRrRgulvllilvlln)(2)()(4222（3.5）汽液界面处液膜的流动速度即为Rr时的u值，有RRRgRRgulvllilvlliln)(2)()2(42（3.6）（2）冷凝段液膜内的温度分布数学描述方程。对能量方程（3.3）积分求解可以得到，液膜内的温度分布：)ln(ln)ln(lnln)ln(lnRRRTRTrRRTTTwvvw（3.7）（3）冷凝段液膜单位弧度质量流量的描述方程。0urdrm将速度分布式代入上式，化简可得24323246231RRRRgmlilll（3.8）（4）冷凝段的通用控制方程。将u、T、m的表达式（3.5）（3.7）（3.8）代入能量平衡方程，积分化解可以得到034CxBAi（3.9）其中，3/422323)(4llpwvfglgcTThA3/1232)/(9)(3lllpwvfglggcTThB3/122)/(lgC（5）冷凝段汽液界面切应力的数学描述。fmi其中，i为界面切应力，m为摩擦切应力，f为相变切应力。)()()(livfgwvllivmuUhTTuU2)(2livvffuUc其中cf为摩擦系数是Rev的函数[6]，其取值对应关系如下表表3-1Rev摩擦系数cf0Rev=2000cf=16/RevvlivvvuUD)(Re。根据质量平衡关系，来计算蒸汽平均流速vU。热管稳态工作时，近似认为在同一水平截面上汽液两相的质量流量相平衡，则有DmUvv4。2000Rev=4000cf=Rev0.33/15254000Rev=30000cf=0.079×Rev-0.2530000Rev=106cf=0.046×Rev-0.2