传热方向与热量符号的讨论

学术   2024-09-01 23:19   山东  

传热方向与热量符号的讨论

如果你正因为热量符号问题困惑不已,这篇文章大概对你有帮助。

判别法则

傅里叶定律如下:

这里面有个负号,初学者看到这个负号可能有些困惑。它表示传热方向与温度梯度方向相反,因为热量总是从高温指向低温,而温度梯度却是从高温指向低温的。同时注意到 自己也是有符号的,表示 的方向就是 轴正方向。为啥呢,假设 轴上两个点 ,且 。对傅里叶定律进行离散:

  • ,则热量从 传到,也就是沿 轴正向传递,同时注意到
  • ,则热量从 传到,也就是沿 轴反向传递,同时注意到

至此我们可以知道 的正负号与热量沿 轴正向/反向传递是一致的。

到这里我们思路还是很清晰的,现在我们引入另一个传热学定律——牛顿冷却定律:

这个地方规定热量进入物体为正。注意这边就没有什么沿 轴正向还是反向了,因为物体是被包围在环境中的,所以热量会从各个方向同时进入/流出物体。

这两种方向的描述好像不太一样,曾经就引起了我的混乱。能不能统一呢?当然可以。我们把傅里叶定律中的 和牛顿冷却定律里的 统一用 表示,显然有 。故两个定律可以统一写成:

显然, 的正负号与热量是沿 的正向还是反向是严格一致的:对于傅里叶定律, 的正向恰为 轴的正向;对于牛顿冷却定律,如果 在环境, 在物体,那么 的正向恰为环境到物体的方向。这样两个判别方法就统一起来了。下面我们用其处理几个一维稳态导热的边界条件,参考陶文铨《数值传热学》(第2版)P92-98。

一维稳态导热的边界条件

按内点法进行离散

控制方程为

其中 为源项。如图所示,使用【外点法】进行离散,右端点为 。第一类边界条件太简单,我们只讨论第二、三类边界条件。这本书都是以进入计算区域的热量为正的。

从定义出发

所以第二类边界条件应该怎么表示呢?重新对方程(1)离散:

此处下标 意为边界(boundary),注意此处是右侧边界,不是左侧边界,对于第二类BC, 为一恒定值。咦,这个等式和方程(2)不一样啊,负号竟然没了???这就是我困惑的开端了。解决这个困惑就要用到前面的判别法则,我们规定了进入为正,因为是右侧边界,所以从右往左为正,而方程右侧的分子刚好就是 方向,这就一致了!

现在看看第三类边界条件,这次 就不再是恒定值了,这个热量由自然对流来提供(牛顿冷却定律):

结合(6)与(7)可得

分析一下符号:左侧方向是 ,右侧方向是,一致的。

从定义出发得到的离散方程具有一阶截差。而且物理意义不是很明确,下面我们用控制容积平衡法,导出更具物理意义的式子。

控制容积平衡法

代表的控制容积(CV)是图中阴影部分,宽度为 。对该CV做热量衡算:

根据前面说的进入为正,我们不妨让所有热量都“名义”地进入CV,所以没有流出的热量,

热量衡算示意1

我们规定由相邻节点 传过来的热量为 ,源项传过来的为 ,环境传过来的为 ,则有:

所以

这个离散方程的物理意义就很明确了。并且书上证明了它是具有二阶截差的。

理解了这个,再回头去看方程(6)和方程(7)也许会有些迷惑,为什么两个含 的项具有相反的符号。其实看图就明白那里的 和这里的 方向就是反的。我们需要知道的是(6)和(7)是从定义出发的,只是对边界处的热量做了一个等价代换而已。切不可用守恒的思想去理解(6)和(7)。

另外如果我们调换一下 的方向(仅对于 ,违反进入为正的规定):

热量衡算示意2

这时就有 了,于是

所以

这样一来,(14)和(6-7)里面含 的项就一致了。并且(14)与(12)是完全等价的。

附加源项法

附加源项法仅适用于内点法的离散。

附加源项法图示

P点离散方程的一般形式为:

为了把未知的 消去,减少一个未知量,在两侧同时减去

由于 ,所以

此处也是规定进入为正, 的方向就是进入方向。

至此,可根据第二类或第三类边界条件对 做替换。


有限元语言与编程
面向科学计算,探索CAE,有限元,数值分析,高性能计算,数据可视化,以及 Fortran、C/C++、Python、Matlab、Mathematica 等语言编程。这里提供相关的技术文档和咨询服务,不定期分享学习心得。Enjoy!
 最新文章