传热与流体流动的数值计算45章.pptVIP

传热与流体流动的数值计算;第四章 热传导;基本方程 稳态一维问题的控制微分方程： 推导出离散化方程 ;-网格间距;-界面导热系数ke;讨论这样一种情况：围绕着网格点P的控制容积由具有均匀导热系数kP的材料填满，围绕着E点的控制容积由导热系数kE的材料填满，对于P点和E点之间的组合板，根据稳态无内热源一维导热的分析，有：;应用于系数的定义式，得到aE：;-非线性;-源项的线性化;4. SC =4+20Tp*3，Sp= -25Tp*2。这一线性化比已知的S~T曲线陡，使收敛速度降低。 四种可能的线性化与实际曲线比较如图：;-边界条件;Boundary conditions; 如果边界上的热流密度qB已知，则要求的对TB的方程变成： 3. 通过放热系数和周围流体的温度来规定边界的热流密度。如果热流密度qB是放热系数h以及环境流体温度Tt规定，那么，方程TB方程变为： ;-线性代数方程的解;TDMA;TDMA;TDMA;TDMA 算法;通用的离散化方程 时间是一个单向坐标，由一已知的初始温度分布开始，沿着时间坐标逐步向前求解：已知t时刻T在网格点上的值，求得t+Δt时刻值。 对整个控制容积积分方程： 得到：;假设用下式归纳一般化有关TP、TE和TW如何随时间由t到t+Δt而变化的关系： 其中f是在0和1之间变化的加权因子。于是：;-显式，克兰克-尼科尔森模式，以及全隐式模式;显式格式Explicit scheme;Crank-Nicolson格式;隐式格式Implicit scheme;-全隐式离散化方程;4.4 Unsteady2-D heat conduction;Discretized unsteady 2-D heat conduction equation;Unsteady 3-D heat conduction equation;Discretized unsteady 3-D heat conduction equation;-三维问题的离散化方程;-代数方程的解;逐行法 把TDMA和高斯-赛德尔法结合起来。选择一条网格行（设在y方向选取这样的网格行），假定沿相邻的行上的T值批最新值构成。用TDMA法求得所选行上的T值。将在同一方向的所有行进行这种计算。如果想做的话，再按相同的方法在其他方向重复上述程序。 以二维为例，如图所示的情况需要注意：;依前后二次迭代之间因变量的变化究竟是被加速还是被减慢的过程称为超松弛或欠松弛。超松弛常用于和高斯-赛德尔法相结合，叫做持续超松弛（SOR）；欠松弛在强烈非线性方程组的迭代求解中用来避免发散。 取T*p作为前一次迭代所得Tp值。 引进松弛因子，得到： 可以根据经验以及对所给定的问题所作的试探性计算求得一个合适的?值。;通用惯量进行松弛。用下面公式代替离散化方程： 式中i是所谓的惯量。对于正的i值，方程具有欠松弛作用；对于负的i则产生超松弛。 ;控制容积面的位置 讨论控制容积面构成的两种不同的替代形式，并讨论它们各自有关的优点。为方便起见，描述针对二维问题。 方法A：控制容积面放在两个网格之间的中点。;方法B： 网格点放在控制容积的中心： 克服了A的缺点。 具有方便性。;我们所提出的这种方法不只限于直角坐标系，还可以用于任意一种正交坐标系。以二维极坐标问题为例，与方程 对应的r?形式是： 其中的网格与控制容积如图示： 设控制容积在z方向厚度为1，方程 两边乘以r，并在整个控制容积范围 内对r和进行积分，得到下面的离散 化方程：;由一个新的坐标系引入的补充特征主要是几何上的特征。;5-1 任务 在通用微分方程中将对流项考虑进去，只要对流项的加入不改变离散化的形式，同样的处理方法仍然适用。 本章任务是：在已知的流场（即速度分量和密度）的情况下，求得对φ的解。 已知流场必须满足连续性方程：;讨论只有对流与扩散这两项存在的情况下的一维稳态问题。 控制微分方程： 应用图示三网点群：;- 预备性的推导;定义两个新的符号： 两者具有相同因次， F表示对流或流动的强度；D是扩散传导性。（注意，D永远为正，而F不同） 离散化方程变为:;讨 论;- 上风方案;不会产生负的系数； 可以把这个方案说成是建立在“槽与管”的模型基础上，管内的流体不会“知道”将要流入那个槽内的任何情况，但它却携带了它所来自那个槽内的全部信息。这就是上风方案的本质。;精确解(Exact solution);不同的贝克列数时的φ~x变化如图;- 指数方案;- 混合方案;用特殊符号 代表其中包含的所有量的最大值。于是;- 幂函数方案;- 一个通用化的公式;如果φi和φi