第21次课-数值解法概述.doc

传热学数值方法概述 非常不幸的是：虽然我们能从数学给出描述流动及传热问题的一般的微分方程，但却不能一般地给出它的解。这就是造成我们理论和实践脱节的根本原因，也正因为如此，才导致我们在流体力学、传热学中所针对的研究对象常常是理想的或简单问题。因为只有引入理想流体或简单问题这样一些假定，可以使上述复杂的方程得到大大的简化，以致于我们在数学上能构求解它们。 一、传热数值计算概述 1、通用偏微分方程 基于质量守恒、动量守恒和能量守恒以及过程变量连续性假说的一组相互耦合的偏微分方程组可以写成如下通式： ，（） 我们可以说这样一组基于质量、动量和能量守恒的基本方程，在配以适当的边界条件即可适合于任何流动及传热问题。但是接下来的问题是，也是最为关键的问题是我们如何求解这组相互耦合的强非线性偏微分方程组。即如何由这样一组方程获得它的解： 在编制计算程序时，我们只需要写出一个求解方程的通用程序，对不同意义的φ就可以重复使用这个程序，当然对不同的φ需要对相应的和S分别赋以各自合适的表达式，如下表。 方程 φ Γ S 连续性方程 1 任意 0 X向动量方程 u μ 能量方程 T或h 或λ 或Sh 湍流-动能方程 K Γk G-ρε 化学组分方程 mi Γi Ri 上述通用微分方程中的四项分别是不稳态项、对流项、扩散项以及源项。因变量可以代表各种不同的物理量，如化学组外的质量分量、焓或温度、速度分量、湍流功能或湍流的长度尺度。与此相应，对于这些变量中的每一个都必须给相对应的扩散系数Γ以及源项S赋以适当的意义。 2、单向与双向的坐标 所谓双向坐标就是指如果在一个坐标上的一个给定位置处的条件，受该位置两侧条件变化影响的坐标。而如果在一个坐标上的一个给定位置处的条件，只受该位置一侧条件变化的影响，这样的坐标就是单向坐标。通常，空间坐标是双向坐标，时间坐标则是单向坐标。但空间坐标有时也可以作为单向坐标。对流是一种单向过程，而扩散则具有双向效应。当流量很大时，对流作用远远大于扩散作用，因而空间坐标就近于单向坐标了。 在一块固体的不稳态冷却期间，其某一给定瞬时的温度只能受在该瞬时以前所发生的那些条件变化的影响。这是一个常识问题。昨天的事情影响着今天发生的事情，但是明天的条件对于今天会发生什么事情则没有影响。 3、方程的类型 包括抛物型、椭圆型与双曲型。从方程类型讲，抛物型方程代表一种单向状态，而椭圆型方程代表一种双向状态。不稳态的传热问题实际上是时间坐标上的抛物型和空间坐标上的椭圆型问题。一个二维的边界层当以流动方向的坐标表示时是抛物型的，而当以横向坐标表示时就是椭圆型的。 这样的描述方法可能读者不习惯，或许用下面的规则可以将它与读者熟悉的分类方法联系起来：至少存在一个单向坐标的状态是抛物型状态；否则就是椭圆型状态。有时将具有一个单向空间坐标的流动称为边界层流动，而全部都是双向坐标的流动则称为回流流动。这里我们把双曲型问题当成是椭圆型问题这一总类中的一个组成部分，即这类问题所对应的全部是双向坐标。 二、数值计算的基本步骤 本世纪四十年代第一台计算机产生以来，特别是近十几年以来，计算机得到了迅速的发展和应用，它为我们求解这些完备的方程提供了数值计算的途径，即面对真实条件和复杂工况，用数值计算的方法求解描述该真实过程的基本方程，从而获得该过程变量的分布特征。 一个传热问题数值求解的总体步骤大致如图所示。 电子计算机中的一切计算都是通过加、减、乘、除四则运算来完成的。为了用计算机解出节点上未知量的近似值，首先需要从给定的微分方程或基本物理定律出发，建立起关于这些节点上来知量近似值之间的代数方程(称为离散方程)，然后对之进行求解。在传热学中所应用助数值计算方法很多，大多数方法的基本思想可以归结为：把原来在时间、空间坐标中连续的物理量的场(如速度场、温度场、浓度场等)，用有限个离散点上的值的集合来代替，按一定方式建立起关于这些值的代数方程并求解之，以获得物理量场的近似解。 数值解法是一种离散近似的计算方法。它所能获得的不象分析解那样是被研究区域中未知量的连续函数，而只是某些代表性地点(称为节点)上的近似值。 同一物理问题的不同数值解法间的主要区别，在于子区域的划分与节点的确定、离散方程的建立及其求解这几个步骤上。传热学所采用的一些数值求解方法，主要的是有限差分法(FDM)、有限元法（FEM）、边界元法（BEM）、有限容积法（FVM）及有限分析法（FAM），前三种方法的比较如下表。 方法 有限差方法 有限元法 边界元法 基本思想 用差商代替微商 变分原理，分段逼近 积分变换 应用领域 适用于形状较简单的场合 适用于物理条件和几何条件复杂的场合 适用于几何条件复杂的场合 特 点 发展相对成熟，物理意义明确，便于实施，对几何形状的适应性差。 计算