D. B. Spalding(1923-)
数值传热学的开拓者和奠基人为当时任职于伦敦大学
帝国理工学院
(Imperial College,University of London)的S.V. Patankar
和D.B. Spalding
。
明尼苏达大学
的E.M. Sparrow教授和W.J. Minkowycz教授对数值传热学的发展起到了重要的促进作用。
中国在数值传热学方面的知名学者有
西安交通大学
的
陶文铨
教授
。
求解对流换热问题的关键是确定流场,为确定流场则需求解反映了
粘性流体
流动的基本力学规律的动量方程,即
纳维-斯托克斯方程
(Navier-Stokes equations, N-S方程)。但N-S是一个
非线性偏微分方程
,求解起来非常困难和复杂,只能得到少数简化情况下的精确解。为解决该问题,学者们在N-S方程的数值解方面开展了一些列研究。
1972年,伦敦大学帝国理工学院(Imperial College,University of London)的S.V. Patankar和D.B. Spalding在总结前人研究的基础上提出了求解N-S方程的“求解压力耦合方程的半隐式方法”,即
SIMPLE算法
。从此该算法在
计算流体力学
与
数值传热学
中得到了广泛应用,并发展出了一些列高精度高稳定性的数值算法
。
2001年,B.Yu、HOzoe、W.Q.Tao(陶文铨)提出了一种加速SIMPLER算法收敛的MSIMPLER方法
。
数值传热学常用的数值方法
历史上最早采用的数值方法,对简单几何形状中的流动与换热问题最容易实施的数值方法。
其基本点是:将求解区域中用于
坐标轴
平行的一系列网格的交点所组成的点的集合来代替,在每个节点上,将
控制方程
中每一个
导数
用相应的差分表达式来代替,从而在每个节点上,形成一个
代数方程
,每个方程中包括了本节点及其附近一些节点上的未知值,求解这些代数方程就获得了所需的
数值解
。
将所计算的区域划分成一系列控制容积,每个控制容积都有一个节点做代表。通过将守恒型的控制方程对控制容积坐积分导出离散方程。在导出过程中,需要对界面上的被求函数本身及其
一阶导数
的构成做出假定,是流动与换热问题的
数值计算
中应用最广的一种方法。
把计算区域划分为一系列原题(在二维情况下,元体多为三角形或四边形),由每个元体上去数个点作为节点,然后通过对
控制方程
做积分来获得离散方程。有限元法最大的优点是对不规则区域的适应性较好。但计算的工作量一般要比有限容积法大,而且在求解流动与换热问题是,
对流
项的离散处理方法及
不可压缩流体
原始变量法求解方面没有
有限容积法
成熟。
由
陈景仁
教授在1981年提出。在这种方法中,也像
有限差分法
那样,用一系列网格线将区域离散,所不同的是每一个节点与相邻4个网格(二维)问题组成计算单元,即一个计算单元由一个中心节点与8个l 邻点组成。在计算单元中把控制方程中的非线性项局部
线性化
,并对该单元上未知函数的变化型线作出假设,把所选定型线表达式中系数和
常数项
用单元边界节点上位置的变量值来表示,找出其分析解。然后利用其分析解,得到该单元中点及其边界上的位置值的
代数方程
,即单元中点的离散方程。
1、给出物理模型;
2、借助基本原理、定律给出数学模型。主要包括质量守恒、能量守恒、动量守恒、
傅立叶定律
、牛顿冷却公式等。
3、对数学模型进行简化和化简;
4、求解区域的离散化;
5、数学模型的离散化;
此步骤需要恰当的方法和建立结点处待求变量近似值之间的代数关系:离散化方程。这也是数值传热学的基本内容,是研究成败的关键;
6、合理假设。因为物理上的简化和数学上的化简有利于抓住主要矛盾;
7、量级分析。分析时需要忽略小量级的项,并且建立边界层方程和粘性耗散函数;
9、运用与分析解对比(简单问题)、实验结果或者前人结果对研究成果的可靠性进行检验;
10、用图线、
动画
等可视化方法对研究成果进行表达和分析。
1、具有成本较低等优势;
2、结果的可靠性取决于模型的正确性和物性数据的可靠性;
3、非常全面和详细的过程;
4、不受实验等条件限制,便于分析单个因素的影响;
5、可以研究其他条件下无法进行试验的研究高温、高压等危险环境和超常尺寸等;