Fluent求解器——热分析

Ronko_G

已于 2023-08-23 17:32:58 修改

阅读量7.2k

点赞数 9

分类专栏： Fluent 文章标签：数学建模算法

于 2023-06-12 17:27:11 首次发布

本文链接：https://blog.csdn.net/Ronko_G/article/details/131170649

版权

Fluent 专栏收录该内容

10 篇文章 52 订阅

订阅专栏

一、热的传递方式

1、热传导

物体各部分之间不发生相对位移时，依靠分子、原子、自由电子等微观粒子的热运动而产生的热能传递称为热传导。导热的基本定律为傅立叶定律：

$\Phi =-\lambda A\frac{\Delta T}{\delta }$

其中， $\Phi$ 为热流量； $\lambda$ 为导热系数；A为面积；T为温度； $\delta$ 为壁面厚度。除水、水溶液和甘油外，其他大多数液体的导热系数随温度升高而下降。

2、热对流

由于流体的宏观运动而引起流体各部分之间发生相对位移，冷热流体相互掺混所导致的热量传递过程称为热对流。热对流只能发生在流体当中，而且由于流体中的分子同时在进行着不规则的热运动，因而热对流必然伴随着热传导。

工程中感兴趣的是流体流过一个物体表面时流体与物体表面之间的热量传递过程，称为对流传热，实际上，大多数情况下所说的热对流也指这种情况。根据引起流动的原因来划分，对流传热可以区分为自然对流和强制对流两大类。对流传热的基本公式为牛顿冷却公式：

$\Phi =hA\Delta T$

其中，h为对流换热系数（表面传热系数）。

3、热辐射

物体由于热的原因而发出辐射能的现象称为热辐射。理论上，只要物体的温度高于绝对零度，物体就会不断的把热能变为辐射能，向外发出热辐射。热辐射的基本公式为斯忒藩 - 玻耳兹曼定律，又称为四次方定律：

$\Phi =\varepsilon A\sigma T^{4}$

其中， $\varepsilon$ 为物体的辐射系数，其值总小于1； $\sigma$ 为斯特藩-玻尔茨曼常数。

二、导热问题的三类边界条件

1、规定了边界上的温度值，称为第一类边界条件，也称为 Dirichlet 条件。此类条件一般规定边界的温度为常数。

2、规定了边界上的热流密度值，称为第二类边界条件，也称为 Neumann 条件。此类条件一般规定边界上的热流密度为常数。

3、规定了边界上物体与周围流体间的表面传热系数h及周围流体的温度，称为第三类边界条件，也称为 Robin 条件。

此外，在处理实际工程问题时，还会遇到辐射边界条件，即物体表面与外界环境之间只发生辐射换热，如航天器上的发热元件向太空的散热。

三、Fluent中热边界的设置

Fluent仿真时，默认壁面为绝热，即通过壁面的热流量为 0。如果想要分析传热，首先打开能量方程，之后可以在 "Wall"-"Thermal" 选项卡中选择，前三种 Heat Flux，Temperature，Convection 分别对应前面所说的第二类、第一类及第三类边界条件。Radiation 为辐射边界条件，Mixed 为混合边界条件。除了正常的边界条件外，Fluent 还可以设定壁面的发热功率，如果壁面为普通的固体壁面，则此处为0。

需要注意，当计算当中存在共轭传热问题时，导入网格时，Fluent 会自动为共轭传热交界面生成 shadow 面，一般情况下，此类壁面不需要进行额外设置。

四、不同传热方式的计算方法

1、热传导的计算

在工程计算中，传热导的问题通常会以热阻设置的情况呈现，热阻的定义为：

$R=\frac{\delta }{\lambda A}$

各参数的含义已在前文说明。Fluent 一般有三种方法处理热阻：

(1) 设置一个有厚度的薄壁，并为其划分网格，设置材料通过求解器来计算热阻的值。这种方法可以考虑到各个方向上的热量传递过程，但往往会使网格数量大幅增加。

(2) 为壁面设置一个虚拟厚度。几何模型中不体现其厚度，因而划分网格时也仅以 0 厚度壁面存在，导入 Fluent 中后，在壁面设置中为其设置一个虚拟的厚度。这种做法可以在考虑热阻的同时大幅降低网格数量，但这种做法只能考虑到垂直于壁面方向的热量传递过程。在 Fluent 中，我们可以通过指定壁面材料、壁面厚度来考虑薄壁热阻对传热的影响。

(3) 设置 Shell Conduction，这种方法类似于 (2)，但可以考虑到各个方向上的热量传递。同时可以指定多层不同材料的薄壁。这种做法在一些工程应用上能够在保证精度的前提下，大大降低网格划分的难度及数量。

2、热对流的计算

(1) 在强制对流计算中，一般需要选取合适的湍流模型。通常推荐使用 Realizable k-e 或者 SST k-w 模型。

(2) 在自然对流计算中，一般以瑞利数的大小来判断流动是否为湍流，公式为：

$Ra=\frac{\beta gL^{3}\Delta T}{\nu \alpha }$

其中， $\beta$ 为膨胀系数，g为重力加速度，L为特征长度， $\nu$ 为运动粘度， $\alpha$ 为扩散率。当 $Ra> 10^{9}$ 时，流动为湍流，需要打开对应的湍流模型。

与一般的仿真计算不同，自然对流需要进行以下设置：

由于自然对流是重力影响下密度的变化所引起的，所以应打开重力项，设置重力的大小和方向；
自然对流的压力离散格式设置为 Body Force Weighted 或者 PRESTO!。使用默认二阶格式会得到错误的结果；
勾选Operating Density，设置参考密度，可以使计算更加稳定；
材料密度的设置：对于气体，密度设置可以选择 Boussinesq 假设或不可压缩的理想气体模型。对于液体，只能选择 Boussinesq 假设。对于封闭区域的自然对流计算，同样需要使用 Boussinesq 假设。Boussinesq 模型假设在动量方程中，除了体积力项之外，其它各项的密度为常数。需要注意的是 Boussinesq 假设只能用在密度变化小于 20%的情况下。打开 Boussinesq 需要在密度设置中选择 boussinesq，并设定参考密度，同时需要设置流体膨胀系数，一般气体的膨胀系数为其热力学温度的倒数。

3、热辐射的计算

首先提出光学厚度（Optical thickness）的概念。光学厚度是介质吸收辐射能力的量度。在 Fluent 中，光学厚度的公式为：

$Optical\; Thickness=(\alpha +\sigma )L$

其中， $\alpha$ 为吸收系数，即由于介质吸收而导致的辐射强度在经过每单位长度的介质后改变的量。由于空气一般不吸收辐射，因此流体介质为空气时该系数可设为 0。 $\sigma$ 为散射系数，即由于介质散射而导致的辐射强度在经过每单位长度的介质后改变的量。同样的，流体介质为空气时该系数可设为 0。L为特征长度。

一般而言，热辐射模型用在高温工况及仅有自然对流存在的工况中。若要在计算中考虑热辐射的影响，需要打开 "Radiation Model" 面板，在其中选取相应的热辐射模型。

热辐射模型的区别
模型	光学厚度	计算量
Rosseland	>5	计算量小，但精度一般，使用较少
P1	>1	计算量较小，在光学厚度较大的问题中计算效果较好
DTRM	全部	计算量小，但无法用于并行计算，使用很少
S2S	0	光学厚度为0时，S2S计算精度与DO模型相当，但计算量小于DO模型。
DO	全部	使用范围广，但计算量最大，计算精度最高