(四)、传热学-热传导问题的数值解法

卓越智工坊

已于 2024-06-19 11:37:05 修改

阅读量750

点赞数 5

分类专栏：传热学学习总结文章标签：算法人工智能机器学习

于 2024-06-19 10:36:49 首次发布

本文链接：https://blog.csdn.net/zcz_822/article/details/139794810

版权

传热学学习总结专栏收录该内容

3 篇文章 0 订阅

订阅专栏

1、基本思想

把原来在时间、空间坐标系中连续的物理的场，用有限个离散点上的值的集合来代替，通过求解按一定方法建立起来的关于这些值的代数方程，来获得离散点上被求物理量的值。

2、基本步骤

(1)、建立控制方程及定解条件；

(2)、区域离散化；

(3)、建立节点物理量的代数方程；

(4)、设立迭代初场；

(5)、求解代数方程组；

(6)、解的分析。

3、内节点离散方程的建立方法

(1)、泰勒级数展开法

以节点(m,n)处的二阶偏导数为例。

(2)、热平衡法

对每个节点所代表的元体用傅里叶导热定律写出其能量守恒表达式。把节点看成是元体的代表。

通过元体的界面所传导的热流量可以对有关的两个节点应用傅里叶定律写出。

$\varPhi _w=\lambda A\frac{\varDelta t}{\varDelta x} \\ \varPhi _w=\lambda \varDelta y\cdot 1\cdot \frac{t_{m-1,n}-t_{m,n}}{\varDelta x}$

类似的可以写出通过其它三个界面e、n、s，而传导给节点(m,n)的热量。对于所研究的问题，元体(m,n)的能量守恒方程为：

$\varPhi _w+\varPhi _e+\varPhi _n+\varPhi _s=0$

带入得：

$\lambda \varDelta y\frac{t_{m-1,n}-t_{m,n}}{\varDelta x}+\lambda \varDelta y\frac{t_{m+1,n}-t_{m,n}}{\varDelta x}+\lambda \varDelta x\frac{t_{m,n+1}-t_{m,n}}{\varDelta y}+\lambda \varDelta x\frac{t_{m,n-1}-t_{m,n}}{\varDelta y}=0$

4、边界节点离散方程的建立代数方程的求解

(1)、边界节点离散方程的建立

内部角点最具代表：

$\lambda \varDelta y\frac{t_{m-1,n}-t_{m,n}}{\varDelta x}+\lambda \frac{\varDelta y}{2}\frac{t_{m+1,n}-t_{m,n}}{\varDelta x}+\lambda \varDelta x\frac{t_{m,n+1}-t_{m,n}}{\varDelta y}+\lambda \frac{\varDelta x}{2}\frac{t_{m,n-1}-t_{m,n}}{\varDelta y}+ \\ \frac{3\varDelta x\varDelta y}{4}\dot{\varPhi}_{m,n}+\frac{\varDelta x+\varDelta y}{2}q_w=0$

当 $\varDelta x=\varDelta y$ 时有：

$t_{m,n}=\frac{1}{6}\left( 2t_{m-1,n}+2t_{m,n+1}+t_{m,n-1}+t_{m+1,n}+\frac{3\varDelta x^2}{2\lambda}\dot{\varPhi}_{m,n}+\frac{2\varDelta x}{\lambda}q_w \right)$

其中 $q_w$ 可以按照绝热边界、不为0、对流边界，三种情况进行代入。

(2)、处理不规则区域得阶梯形逼近

(3)、求解代数方程的迭代法

5、非稳态导热问题的数值解法

(1)、时间-空间区域的离散化(泰勒展开法)

将函数t在节点(n,i+1)对点(n,i)作泰勒展开，可有：

$t_{n}^{\left( i+1 \right)}=t_{n}^{\left( i \right)}+\left. \varDelta \tau \frac{\partial t}{\partial \tau} \right|_{n,i}+\frac{\varDelta \tau ^2}{2}\left. \frac{\partial ^2t}{\partial \tau ^2} \right|_{n,i}+...$

于是有：

$\left. \frac{\partial t}{\partial \tau} \right|_{n,i}=\frac{t_{n}^{\left( i+1 \right)}-t_{n}^{\left( i \right)}}{\varDelta \tau}+O\left( \varDelta \tau \right)$

式中，符号 $O\left( \varDelta \tau \right)$ 表示余项中 $\varDelta \tau$ 的最低阶为一次。由上式可在点(n,i)处一阶导数的一种差分表示式。

$\left. \frac{\partial t}{\partial \tau} \right|_{n,i}=\frac{t_{n}^{\left( i+1 \right)}-t_{n}^{\left( i \right)}}{\varDelta \tau}$ (向前差分)

另外可得向后差分和中心差分：

$\left. \frac{\partial t}{\partial \tau} \right|_{n,i}=\frac{t_{n}^{\left( i \right)}-t_{n}^{\left( i-1 \right)}}{\varDelta \tau}$ (向后差分)

$\left. \frac{\partial t}{\partial \tau} \right|_{n,i}=\frac{t_{n}^{\left( i+1 \right)}-t_{n}^{\left( i-1 \right)}}{2\varDelta \tau}$ (中心差分)

另外可得稳态和非稳态二阶中心差分：

$\frac{\partial ^2t}{\partial x^2}=\frac{t_{n+1}^{}-2t_{n}^{}+t_{n-1}^{}}{\varDelta x^2}$ (稳态二阶中心差分)

$\frac{\partial ^2t}{\partial x^2}=\frac{t_{n+1}^{i}-2t_{n}^{i}+t_{n-1}^{i}}{\varDelta x^2}$ (非稳态二阶中心差分)

(2)、显式格式

对于一维非稳态导热方程，如扩散项取中心差分(并取i时刻之值)，非稳态项取向前差分，则有：

$\frac{t_{n}^{\left( i+1 \right)}-t_{n}^{\left( i \right)}}{\varDelta \tau}=a\frac{t_{n+1}^{i}-2t_{n}^{i}+t_{n-1}^{i}}{\varDelta x^2}$

显示差分格式：求解非稳态导热方程就是从已知的初始温度分布出发，根据边界条件依次求得以后各个时间层上的温度值。一旦i时层上各个节点的温度已知，可立即算出i+1时层上各内点的温度，而不必求解联立的方程。

优缺点：计算工作量小，对时间步长及空间步长有一定的限制，否则会出现不合理的震荡的解，称为稳定性问题。

(3)、隐式格式

对于一维非稳态导热方程，如扩散项取中心差分(并取i+1时刻之值)，非稳态项取向前差分，则有：

$\frac{t_{n}^{\left( i+1 \right)}-t_{n}^{\left( i \right)}}{\varDelta \tau}=a\frac{t_{n+1}^{\left( i+1 \right)}-2t_{n}^{\left( i+1 \right)}+t_{n-1}^{\left( i+1 \right)}}{\varDelta x^2}$

隐式差分格式：已知的是i时层的值 $t_{n}^{\left( i \right)}$ ，而未知量有3个，因此不能直接算出 $t_{n}^{\left( i+1 \right)}$ 之值，而必须求解(i+1)时层的一个联立方程组才能得出(i+1)时层各节点的温度。

优缺点：计算工作量大，但对步长没有限制，不会出现解的振荡现象。

(4)、边界节点的离散方程(能量守恒法)

无限大平板得右边界部分，其表面流体冷却，表面传热系数为h。

对元体应用能量守恒定律得：

$\varPhi _{\text{right}}+\varPhi _{\text{left}}=\varPhi _{\text{chang}}$

$\lambda \varDelta y\frac{t_{N-1}^{i}-t_{N}^{i}}{\varDelta x}+h\varDelta y\left( t_f-t_{N}^{i} \right) =\rho cV\frac{\partial t}{\partial \tau}=\rho c\frac{\varDelta x}{2}\varDelta y\frac{t_{N}^{i+1}-t_{N}^{i}}{\varDelta \tau}$