第十三周：CFD基础理论+Fluent实例学习

目录

摘要

Abstract

Fluent部分

相关概念学习

热传导

热对流

自然对流

热辐射

案例部分

前处理

后处理

CFD理论部分

总结

---

摘要

在本周中，主要进行了对CFD中的有限差分法的补充学习，关注其中对激波进行捕捉的方法，包括：人工粘性法，TVD格式、群速度控制格式（GVC）等激波捕捉格式，其中常用的有TVD和WENO格式。在Fluent实例的学习中，进行了关于自然对流热传导和热辐射的实例演算，明白了如何对其进行设置。

Abstract

This week, I mainly studied about the finite difference method in CFD, focusing on the methods of shock wave capture, including: artificial viscosity method, TVD format, group velocity control format (GVC) and other shock wave capture formats, among which TVD and WENO formats were commonly used. In the study of Fluent examples, I performed an example calculation on natural convective heat conduction and heat radiation, and understood how to set them up.

Fluent部分

相关概念学习

热传导

热传导由傅里叶定律决定，当能量方程被激活时，在所有流体和固体区域中，Fluent计算热传导傅里叶定律指出，传热率与温度梯度成正比。

根据不同仿真精度的需求，可把其定义为温度相关的变量（热导率）一些极端的现象时也可定义为与速度、压力相关的变量，需通过UDF（用户自定义函数）实现。对于复合材料，其不同方向上的热导率是不一样的，可通过定义矩阵的方式来定义各项异性热导率（定义坐标系后在定义方向）。

热对流

随着流体的移动，它会带走热量，这就是对流现象，热传递与流体流动相耦合。能量 + 流动方程激活时可用于求解对流传热，能量方程激活时，流体中热传导也将求解。对流传热速率很大程度上取决于流体速度流体性质可能随温度(例如空气)，而显著变化壁面处传热系数通过湍流热壁函数计算。

对于对流换热系数，如果能做到给定换热系数， 且知道固体的温度，即可直接根据公式进行计算换热量（换热系数稳定）。实际情况很复杂，与流体性质、流动情况以及温度差有关。对稍微复杂的表面，很难给定一个稳定的数值，同时考虑流动和换热。对于fluent来说，需构建两套网格，第一套网格包括流体区域的速度情况的计算，第二套就是固体区域的网格，用于绘制计算换热量。

自然对流

当流体密度与温度有关时，被加热后密度变化，就会发生自然对流，流动是由密度差引起的重力驱动的;当在Fluent中激活重力项时，动量方程的压力梯度和体积力项重写为:

自然对流的流速流量与很多因素有关，使用情况：没有强制对流的动力源，不太好收敛。其中常用的有如下两种模型：

Boussinesq 密度模型：把密度当作一个不变的量，重力直接和温度相关，将密度和压力很好的分割开来，提高了线性程度（更多的适用于温度变化不太大，10以内的流体）。

Incompressible ldea Gas 密度模型：温度和密度成反比，和压力无关。（在温度比较小，压力变化小适用）。共轭换热为固体和流体的共同换热情况，在流体与固体区域之间的壁面或两侧有流体的壁面，在读取网格文件时，Fluent会自动创建 wall 与wall shadowwall shadow为流体边界，两者共享同样的数值，画好固体和流体的网格即可。

热辐射

可见光和大部分的红外线都包含热辐射，对于半透明物体(例如玻璃，燃烧产物气体)，辐射是体积现象:因为辐射可以从体内逸出对于不透明的物体来说，辐射本质上是一种表面现象:因为几乎所有的内部辐射都在体内被吸收。

热辐射通量计算公式如下所示：

σ为Stefan-Boltzmann常数，5.6704x10-8 W/(m2.K4)。

估计系统中传导或对流传热的大小:

与对流和传导的换热量相当时需要考虑热辐射，否则不考虑。在fluent中有六种热辐射模型，常用的有DO模型和S2S模型，DO模型适用于各种场景，计算量较大，S2S模型，适用于电子产品散热。因为空气中大部分是双原子分子，对热辐射的吸收很少，可以忽略，只考虑热表面和冷表面的辐射交换吸收，热辐射模型可极大降低计算量。

案例部分

前处理

电子芯片散热的案例基于官方教材使用SDM和SUMING视频工作流处理模型关注整个计算区域的最高温度包括发射体、机箱、散热器和外壳等组成部分的温度分布情况。几何模型如下：

因为不存在曲面，故选择proximity（临近度）函数，选择face 设定为1，即两个面之间至少有一个网格。选择发热体和热吸收体，进行局部加密，完成后划分面网格，最大尺寸定义为20mm，面网格划分如下：

要注意的是，proximity函数为间隙填充单元层数，即无论多窄的区域都必须有（待定义）的网格，在本案例中设定为0.5.选0.5的网格能够满足仿真需求，尤其是在固体区域上表面网格的质量非常好，虽然里面有差别，但足够了。其对比图如下所示：

设置体网格最大尺寸为12.5mm，其余设置保持默认，划分体网格：

后处理

首先打开能量模型，发热率在75W，但其粘性属于层流，故选择层流模型进行计算。材料定义:定义三种材料: 散热器、CPU和PCB板。其材料为铜、fr-4、coponent，参数如下：

进入边界条件定义，PCB板的材料为fr-4，只导热，不发热，散热片材料为铜，也不存在发热情况。热源的材料为coponent，既导热又发热，选择源项，输入的是单位体积的总量，但还需要除以CPU的体积，可以自定义一个计算函数来进行计算，定义的函数如下：

后设置进出口温度，都为45摄氏度。后对自然对流的条件进行设置，在求解自然对流（或任何重力开启时）时，必须使用“body force weighted”或“PRESTO！”进行压力离散化，这样做的原因是其他选项可以预测进入墙壁的非物理流动右侧的矢量图来自不同的情况，例如，不是本教程在大多数自然对流问题中，PRESTO！和body force weighted同样好。其设计参数如下：

工作密度通常是指当前环境下的气体密度，通过理想气体运算的出45摄氏度下的工作密度为1.11，创建一个点，采用该点的平均温度来确定是否收敛：

伪瞬态，也包括一些瞬态计算的过程当中，固体的时间尺度远远大于流体，因为流体的流动时间太长，而固体只算导热，线性强，所以选择流体的时间尺度选择固体的千分之一来计算。

报告图如下，可见并未收敛，调整残差重新计算：

调整残差至10-9后，再次进行200次迭代运算，发现各指标都趋于平稳，完成收敛：

为包括热辐射因素的温度云图如下：

选择观察未考虑热辐射时的壁面的温度分布，不合理，保持从前数据，开启辐射模型再次进行计算

选择S2S模型进行计算，每个流体边界都会有一个发射率，只需关注流体，固体无需关注，由于在选择辐射方程时，选择的其迭代次数为计算10次辐射方程来算一次辐射，故会出现曲线弹跳的现象：

两者的仿真结果对比如下图所示：

CFD理论部分

总结

本周主要尝试了热传导和对流方面的Fluent案例学习，明白了其常用模型的理论基础和实际运用间的差别，几乎无需进行多处修改，在fluent的默认参数中即可进行收敛计算。但在网格划分的部分中，网格质量小于案例中的质量，现在原因还未找出，但版本不同的概率较大。在CFD理论中，对有限差分法的理论和数学物理基础进行了更深入的学习，受益匪浅。