计算流动与传热通用控制方程研究 李旺.doc

中国工程热物理学会 传热传质 学术会议论文 编号：113593 计算流动与传热通用控制方程研究( 李旺1 宇波 王情愿1 李瑞龙王欣然1 陶文铨2 ( 1. 中国石油大学工程学院，城市油气输配技术北京市重点实验室，北京 102249( 2. 西安交通大学能源与动力工程学院010 E-mail:yubobox@cup.edu.cn) 摘要: 在流动与传热数值计算中往往采用通用形式的控制方程来简化离散和编程，这种方法在数值计算中取得了巨大的成功。本文对现有的通用控制方程进行了研究，发现目前所采用的通用控制方程对能量方程的处理不尽合理，该处理在比热容不为常数时不能保证能量守恒，导致比热容变化时计算结果不正确。本文提出了广义密度的概念，给出了新的通用控制方程，保证了能量守恒性，并通过多个计算算例说明了本研究提出的方法的正确性和可行性。 关键词: SIMPLE算法；通用控制方程；守恒；广义变量；流固耦合 0 前言 在采用数值方法求解流动与传热问题时，各控制方程一般采用共同的形式表达，这种形式称为流体流动与传热的控制方程通用形式。各控制方程的不同则通过广义扩散系数和广义源项的不同来体现，这里的“广义”表示不必是原来物理意义上的量，而是数值计算模型中的一种定义。采用方程的通用形式具有明显的优点，可以对所有的控制方程用统一的形式离散和编程，显著地提高编程的效率，同时也使程序的通用性大大提升。 在流动与传热数值计算的过程中广为采用的控制方程的通用形式如（1）式所示[1,2,3]： （1） 方程左边第一项为非稳态项为对流项；方程右边前项为扩散项为源项，它是一个广义量，代表了那些不能包括到控制方程非稳态项、对流项与扩散项中的所有其它各项之和。式中通用变量，u，v，w，T等求解变量；和是与相对应的广义扩散系数及广义源项，具体含义见表1。表1中，为向坐标，为向坐标，分别为向速度为压力、分别为x方向和y方向上的体积力，为密度、和分别为动力粘性系数。表1的系数及源项方程 方程 1 0 0 动量方程 动量方程 能量方程 （2） 该方程在比热容cp为常数时，可以变换为通用控制方程中表1的形式，也就是 （3） 方程（3）与方程（2）在比热容cp不为常数时是不等价的。 将方程(3)×cp-(2)，可以得到 （4） 要保证方程（3）和方程（2）在比热容不为常数时等价，需满足、以及成立，显然这三个关系式是不可能恒成立的，因此，方程（3）是一个不守恒的方程，在比热容不为常数时应该慎用。 当比热容随温度变化较大时以及在计算区域中存在流固耦合区域时，方程（2）和方程（3）差异会很大，导致方程（3）计算的结果失真。为保证在任何比热容cp时都得到正确的结果，应该构建新的通用控制方程形式，我们提出如下的通用控制方程： （5） 方程（5）中，为广义密度项，和也更加符合物理意义。下面通过三个计算实例来对比分析两种通用控制方程计算结果的差异：1）强制对流换热，2）自然对流换热，3）流固耦合问题。 表控制方程中的系数及源项方程 方程 1 0 0 动量方程 动量方程 能量方程 ，其中Nu和Nu*分别为采用常用通用控制方程（1）和新提出的通用控制方程（5）计算得到的Nu数。 2.1 强制对流换热 本算例为边长为L的方腔顶盖驱动强制对流换热问题，L=0.01m。为了考察比热容cp不为常数时两种通用方程的不同，这里选取比热容cp随温度变化不能忽略的乙醇作为流动介质。它的主要物性参数为：ρ=797.9(kg/m3)， cp=0.0407T2+6.349T+2231.4(J/kg·℃)，λ=(0.185-2.66×10-4T)(W/·℃) , μ=(2×10-7T2-3.08×10-5T+1.7×10-3)(Pa·s)， β=1.1×10-3(℃-1)，这些物性参数的适用范围为0~60℃。计算中温度边界条件为：方腔底部：T=60℃；顶盖：T=0℃；方腔两侧壁均为绝热边界条件。速度边界条件除顶盖外均为无滑移边界条件，顶盖驱动速度取两种情况：u=0.08m/s，0.8m/s。 图1对比了温度等值线，图中点划线和实线分别表示采用常用通用控制方程（1）和新提出的通用控制方程（5）的计算结果，从图中可