第08章fluent基本物理模型.docx
《第08章fluent基本物理模型.docx》由会员分享，可在线阅读，更多相关《第08章fluent基本物理模型.docx（53页珍藏版）》请在冰豆网上搜索。
第08章fluent基本物理模型
第08章fluent基本物理模型
本章介绍了FLUENT所提供的差不多物理模型以及相关的定义和使用。
差不多物理模型概述
FLUENT提供了从不可压到可压、层流、湍流等专门大范畴模拟能力。
在FLUENT中,输运现象的数学模型与所模拟的几何图形的复杂情形是结合在一起的。
FLUENT应用的例子包括层流非牛顿流的模拟，涡轮机和汽车引擎的湍流热传导，锅炉内煤炭粉碎机的燃烧，可压射流，空气动力外流，以及固体火箭发动机的可压化学反应流。
为了与工业应用相结合，FLUENT提供了专门多有用的功能。
如多孔介质，块参数（风扇和热交换），周期性流淌和热传导，涡流，以及移动坐标系模型。
移动参考系模型能够模拟单一或者多个参考系。
FLUENT还提供了时刻精度滑动网格方法以及运算时刻平均流淌流场的混合平面模型，滑动网格方法在模拟涡轮机多重过程中专门有用。
FLUENT中另一个专门有用的模型是离散相模型，那个模型何以用于分析喷雾和粒子流。
，多项流模型能够用于推测射流的破散以及大坝塌陷之后流体的运动，气穴现象，沉淀和分离。
湍流模型是FLUENT中专门重要的一部分，湍流会阻碍到其它的物理现象如浮力和可压缩性。
湍流模型提供了专门大的应用范畴，而不需要对特定的应用做出适当的调剂，而且它涵括了其它物理现象的阻碍，如浮力和可压缩性。
通过使用扩展壁面函数和区域模型，它能够对近壁面的精度问题有专门好的考虑。
各种热传导模式能够被模拟，其中包括具有或不具有其它复杂性如变化热传导的，多孔介质的自然的、受迫的以及混合的对流。
模拟相应介质的辐射模型及子模型的设定通常能够将燃烧的复杂性考虑进来。
FLUENT一个最强大的功能确实是它能够通过耗散模型或者和概率密度函数模型来模拟燃烧现象。
关于燃烧应用十分有用的其它模型也能够在FLUENT中使用，其中包括碳和液滴的燃烧以及污染形成模型。
连续性和动量方程
关于所有的流淌，FLUENT差不多上解质量和动量守恒方程。
关于包括热传导或可压性的流淌，需要解能量守恒的附加方程。
关于包括组分混合和反应的流淌，需要解组分守恒方程或者使用PDF模型来解混合分数的守恒方程以及其方差。
当流淌是湍流时，还要解附加的输运方程。
本节所介绍的是层流流淌的守恒方程（在惯性（无加速度）的坐标系中）。
后面几节将会讨论热传导、湍流模拟以及组分输运的守恒方程。
关于旋转坐标系中的方程将在移动区域的流淌中介绍。
欧拉方程用于解决无粘流淌，将在无粘流淌一节中介绍
质量守恒方程
质量守恒方程又称连续性方程：
该方程是质量守恒方程的一样形式，它适用于可压流淌和不可压流淌。
源项S_m是从分散的二级相中加入到连续相的质量（比方说由于液滴的蒸发），源项也能够是任何的自定义源项。
二维轴对称问题的连续性方程为：
具体各个变量的意义能够参阅相关的流体力学书籍，其中有具体而详细地介绍。
动量守恒方程
在惯性（非加速）坐标系中i方向上的动量守恒方程为[8]：
其中p是静压，t_ij是下面将会介绍的应力张量，rg_i和F_i分别为i方向上的重力体积力和外部体积力（如离散相相互作用产生的升力）。
F_i包含了其它的模型相关源项，如多孔介质和自定义源项。
应力张量由下式给出：
上式的物理意义能够参阅流体力学教科书，其中会讲得专门清晰。
关于二维轴对称几何外形，轴向和径向的动量守恒方程分别为：
以及
其中：
w是漩涡速度（具体能够参阅模拟轴对称涡流中漩涡和旋转流淌的信息）
热传导
FLUENT承诺在你模型的流体和/或固体区域包含热传导。
本节中所介绍的物理模型和相关输入能够处理从流体内热混合到复合固体的热传导等问题。
自然对流问题会在浮力驱动流淌一节介绍，自然对流与辐射模型将在辐射模拟一节介绍
FLUENT能够推测周期性几何外形的热传导，如密集的热交换器，它只需要考虑单个的周期性模块进行分析。
关于如此流淌的处理，需要使用周期性边界条件，具体能够参阅周期性流淌和热传导一节。
在两个分离的流淌区域解决热传导问题
假如所模拟的流淌包括了两个流体区域，其中被固体区域或者壁面分离开，如下图所示，你需要更细心的定义问题。
要紧需要指定：
●两个流体区域都不能够使用质量出口边界条件
●每一个流体区域能够选择不同的流体材料。
（然而关于组分运算，你只能在整个区域选择唯独一种混合材料）
Figure1:
典型的逆流热交换，在两个流体区域包括了热传导
理论
能量方程
FLUENT所解的能量方程的形式为
其中k_eff是有效热传导系数（k+k_t,其中k_t是湍流热传导系数，依照所使用的湍流模型来定义），J_j^'是组分j^'的扩散流量。
上面方程右手边的前三项分别描述了热传导、组分扩散和粘性耗散带来的能量输运。
S_h包括了化学反应热以及其它用户定义的体积热源项。
在上面的方程中：
其中，理想气体的显焓定义为：
关于可压流为：
在方程5和7中，m_j^'是组分j^'的质量分数，而且
其中T_ref为298.15K.
PDF模型的能量方程
当激活非绝热PDF燃烧模型时，FLUENT解总焓形式的能量方程：
假定Lewis数为1，右手边第一项表示传导和组分扩散项。
非守恒形式的粘性扩散项的奉献由第二项描述。
总焓H定义为：
其中m_j^'为组分j^'的质量分数，而且
是在参考温度T_ref,j^'下组分j^'的生成焓
包括压力作用和动能项
能量方程中的方程1包含了不可压流淌中常常忽略的压力作用和动能项。
因此，在默认的情形下，分离解算器在解不可压流淌时不考虑压力作用和动能项。
假如你期望考虑这些作用，能够使用define/models/energy?
文本命令将所需的项激活。
模拟可压流或者使用耦合解算器时，压力作用和动能项总是压考虑的。
包括粘性耗散项
能量方程中的方程1和PDF模型的能量方程中的方程1包括了粘性耗散项，该项所描述的是粘性剪切所产生的热能。
使用分离解算器时，FLUENT默认的能量方程不包括它（因为粘性热能够忽略）。
当Brinkman数Br接近或者大于一，粘性热将会专门重要。
其中：
DT为系统温度的差分。
你需要考虑粘性耗散项同时使用分离解算器，你需要在粘性模型面板激活粘性热项。
关于可压流淌一样有Br
1。
然而需要注意的是，当使用分离解算器时，假如你定义了可压流淌模型，FLUENT并不自动激活粘性耗散项。
当使用耦合解算器时，所解的能量方程总会包含粘性耗散项。
包括组分扩散项
能量方程一节中的方程1和PDF模型的能量方程一节中的方程1包括了由于组分扩散而导致的焓的输运的阻碍。
当使用分离解算器时，在默认情形下，
会包含在能量方程一节的方程1中。
假如你不想包括它，你能够在组分模型面板中关闭扩散能量源项的选项。
当使用非绝热PDF燃烧模型时，该项并不是显式的显现在能量方程中，因为关于PDF模型的能量方程一节中的方程1来说，该方程右手边的第一项差不多包含了它。
当使用耦合解算器时，该项总是包含在能量方程中。
由于化学反应产生的能量源项
能量方程一节中的方程1的能量源项S_h包括了由于化学反应而产生的能量源项：
其中h^o_j^'是组分j^'的生成焓，R_j^'是组分j^'的体积生成速度。
非绝热PDF燃烧模型的能量方程中，焓的定义差不多包括了能量的生成（见PDF模型的能量方程一节中的方程5，因此能量的反应源项不包括在S_h中。
.
由于辐射产生的能量源项
当使用某一辐射模型时，能量方程一节中的方程1和PDF模型的能量方程一节中的方程1的S_h也包括了辐射源项。
详情参阅辐射模型一节。
相间的能量源项
需要注意的是，能量源项S_h还包括连续和离散相之间的热传导。
在后面的离散与连续相耦合一节将会详细讨论。
壁面处热传导的边界条件
壁面处热传导边界条件在标准壁面函数一节中讨论。
固体区域的能量方程
FLUENT所用的固体区域的能量输运方程的形式为：
其中r=密度
h=显焓（integral_T_ref^Tc_pdT）
k=传导系数
T=温度
q（dot）^"'=体积热源
方程1左手边的第二项表达了由于固体的平移和旋转而导致的能量对流热传导。
速度场u_i由指定固体区域的运动运算出来（见固体条件一节）。
方程1右手边的项分别是固体内部热传导流量和体积热源的热流量。
固体的各向异性热传导
当使用分离解算器时，FLUENT承诺你制定固体材料的各向异性热传导系数。
固体的各向异性传导项形式为：
其中k_ij是热传导系数矩阵。
关于固体材料的各向异性热传导系数的制定能够参阅固体的各向异性热传导系数一节。
入口处的扩散
入口处能量的净输入既包括对流部分也包括扩散部分。
对流部分由你所指定的入口温度确定。
扩散部分依靠于运算出温度场的梯度。
因此扩散部分（相应的净入口输运）不是提早指定的。
在某些情形下，你可能期望指定入口处的能量净输运而不是入口温度。
假如你使用分离解算器，你能够通过取消入口能量扩散来实现这一目标。
在默认的情形下，FLUENT在入口处会考虑能量的扩散流量。
要关闭入口扩散，能够使用文本命令：
define/models/energy?
。
假如你使用耦合解算器，入口扩散选项无法关闭。
热传导所需的用户输入
当FLUENT模型包含了热传导，你需要激活相关的模型，提供热边界条件，并输入操纵热传导和/或随温度变化的材料属性。
本节将会介绍这些输入。
下面将会介绍热传导问题的设定步骤。
（注意：
本步骤只包括热传导模型设定的必须步骤，你还要设定其它的模型，边界条件等。
）
1.要激活热传导的运算，请在能量面板中打开激活能量方程选项。
菜单：
Define/Models?
Energy...。
Figure1:
能量面板
2.（可选，只用于分离解算器）假如你模拟粘性流淌，而且期望在能量方程中包括粘性热传导项，请在粘性模型面板中打开粘性热传导项。
如包含粘性耗散一节中所述，当使用分离解算器时，FLUENT在默认的情形下会忽略能量方程中的粘性热传导项（假如使用耦合解散器，则会一直包含粘性热传导项。
当流体中的剪切应力较大（如：
润滑问题）和/或速度较高、可压流淌，就应该激活粘性耗散项（见包含粘性耗散项一节中的方程1）。
菜单Define/Models/Viscous...
3.在流淌入口、出口和壁面处定义热边界条件。
菜单：
Define/BoundaryConditions...。
在流淌的出入口你需要设定温度，在壁面处你可能需要设定下面的某一热条件：
●指定热流量
●指定温度
●对流热传导
●外部辐射
●外部辐射和外部对流热传导的结合
定义壁面处热边界条件一节详细地介绍了操纵热边界条件的模型输入。
入口处默认的热边界条件为指定的温度300K；壁面处默认的条件为零热流量（绝热）。
关于边界条件的输入请参阅边界条件一章。
4.定义适合于热传导的材料属性。
菜单：
Define/Materials...
如物理属性一节所述，必须定义热容和热传导系数，而且你能够指定专门多属性为温度的函数。
温度的上下限
出于稳固性考虑，FLUENT包括了推测温度范畴的限制。
设定温度上下限的目的是为了提高运算的稳固性，从物理意义上说，温度应该处于已知极限的范畴之内。
有时候方程中间解会导致温度超出这些极限，现在就无法专门好的定义属性。
温度极限保证你的问题的温度在期待的范畴之内。
假如运算的温度超出最大极限，那么所储备的温度就会固定在最大值处。
默认的温度上限是5000K。
假如运算的温度低于最小极限，那么储备的温度就会固定在最小值处。
默认的温度下限是1K。
假如你所预期的温度超过5000K，你应该使用解限制面板来增加最大温度。
菜单：
Solve/Controls/Limits...。
热传导的解过程
尽管使用Fluent默认的解参数能够成功的解决专门多简单的热传导问题，你依旧能够使用本节所提供的指导方针来加速收敛速度和解的稳固性。
能量方程的亚松驰
使用分离解算器时，FLUENT能够使用你在解操纵面板所定义的亚松驰参数来处理亚松驰能量方程，具体能够参阅设定放松因子一节所介绍的内容。
菜单：
Solve/Controls/Solution...。
假如使用非绝热PDF模型，你需要像通常一样设定能量亚放松因子，然而你也能够设定温度的亚放松因子，其用法和解焓方程时温度的亚松驰一节所介绍的一样。
FLUENT可不能管所解能量方程是温度依旧焓形式，它都会设定默认的亚放松因子为1.0。
在能量场阻碍流体流淌（通过温度相关属性或者焓）的问题中，你应该是用较小的亚放松因子，一样在0.8到1.0之间。
当流场和温度场解耦时（没有温度相关属性或者浮力），你能够保留放松因子的默认值1.0。
解焓方程时温度的亚松驰
当解焓形式的能量方程时（即当你使用非绝热PDF燃烧模型时），FLUENT也对温度进行亚松驰，也确实是说，只是用焓（亚松驰）变化对应的温度变化的某一分数来更新温度场。
当你期望焓场变化较快时，二层的亚松驰专门有用，只是温度响应比较之后，相应的温度对流场的阻碍也会滞后。
FLUENT关于温度的亚松驰默认设定为1.0，此设定使用解操纵面板来实现。
屏蔽组分扩散项
假如使用分离解算器来解决组分输运，而且遇到了收敛困难，你应该考虑在组分模型面板中关闭扩散能量源项。
菜单：
Define/Models/Species...。
当改选项关闭时，FLUENT会忽略能量方程的组分扩散阻碍。
注意：
当使用耦合解算器时组分扩散阻碍总会被考虑到的。
步进解
最为有效的推测热传导策略是先运算等温流淌然后加入能量方程的运算。
步骤稍有不同，要紧取决于流淌和热传导是否耦合。
假如流淌和热传导是解耦的（没有温度相关属性或浮力），你能够第一解等温流淌（关闭能量方程）来产生收敛的流场解，然后单独解能量输运方程。
注意：
因为耦合解算器总是一起解流淌和能量方程，因此单独解能量方程只应用于分离解算器。
你能够在解操纵面板中的方程列表中取消能量选项来临时关闭流淌方程或者能量方程（请参阅步进解一节）。
菜单：
Solve/Controls/Solution...。
假如流淌和热传导是耦合的（也确实是模型中包括温度相关属性或浮力），你能够在打开能量方程之前第一解流淌方程。
一旦你有了收敛的流场解，你就能够打开能量选项然后同时解流淌和能量方程完成热传导的模拟。
热传导的报告
FLUENT为热传导模拟提供了附加的报告选项。
你能够生成图形或者报告下面的变量或函数：
●静温
●总温
●焓
●相对总温
●壁面温度（内部表面）
●壁面温度（外部表面）
●总焓
●总焓误差
●熵
●总能量
●内能
●表面热流量
●表面热传导系数
●表面努塞尔（Nusselt）数
●表面斯坦顿（Stanton）数
上面所示的前11个变量包含在后处理面板中的变量选择下拉列表的温度类别中，剩下的变量在壁面流量类别中。
关于它们的定义能够参阅流场函数定义一节。
在报告和显示中焓与能量的定义
焓与能量报告值的定义是不同的，它取决于流淌可压与否。
完全的定义请参阅流场变量及其定义的列表。
报告通过边界的热传导
你能够使用流量报告面板来运算通过每一个边界的热传导或者将通过所有边界的热流量加起来来检查热平稳。
菜单：
Report/Fluxes...。
举荐检查热平稳以确认你的解是收敛的。
关于流量报告的生成请参阅通过边界的流量一节。
报告通过表面的热传导
你能够使用曲面积分面板（在曲面积分一节介绍）来运算通过任何边界的热传导或者运算通过曲面的热传导，那个曲面能够在显示和报告曲面数据一节中介绍的方法来创建。
菜单：
Report/SurfaceIntegrals...。
要报告焓的流速
在曲面积分面板选择流淌速度选项，选择焓（在温度类别中）作为流场变量，然后选择需要积分的一个或多个曲面。
报告平均热传导系数
曲面积分面板还能够报告在曲面上的平均热传导系数h，菜单：
Report/SurfaceIntegrals...。
在曲面积分面板中选择平均选项，选择曲面热传导系数（在壁面流量类别中）作为流场变量然后点击相应的曲面。
浮力驱动流淌和自然对流
当加热流体，而且流体密度随温度变化是，流体会由于重力缘故的而导致密度的变化。
这种流淌现象被称为自然对流（或者混合对流），Fluent能够模拟这种流淌。
理论
能够用Grashof数Reynolds雷诺数的比值来度量浮力在混合对流中的作用：
当那个数接近或者超过一，你应该考虑浮力关于流淌的奉献。
反之，你就能够忽略浮力的阻碍。
在纯粹的自然对流中，浮力诱导流淌由瑞利数（Rayleigh）度量：
其中热膨胀系数为：
热扩散系数为：
Rayleigh数小于10^8说明浮力诱导为层流流淌，当瑞利数在10^8到10^10之间就开始过渡到湍流了。
Boussinesq模型
关于专门多自然对流流淌，你能够用Boussinesq模型来得到更好的收敛速度，它要比设定密度为温度的函数来解决问题收敛得快。
除了动量方程的浮力项之外，该模型在所有解决的方程中将密度看成常数。
动量方程为：
其中r_0是流淌的常数密度，T_0是操作温度，b是热扩散系数。
上面的方程是通过Boussinesq近似等于r_0（1-bDT）来排除浮力项中的r得到的。
只要真实密度变化专门小，该近似是专门精确的.
使用Boussinesq模型的时机
在封闭区域使用Boussinesq模型来运算时刻相关自然对流是专门必要的。
假如温度变化专门小，该模型也能够用于定常问题。
Boussinesq模型不能用于组分，燃烧和反应流淌的运算。
浮力驱动流淌的用户输入
在混合或自然对流中，你必须提供下面的输入来考虑浮力问题：
1.在能量面板中打开能量方程选项。
菜单：
Define/Models/Energy...。
2.在操作条件面板（下图）中打开重力选项，并在每一个方向上输入相应的重力加速度数值。
菜单：
Define/OperatingConditions
Figure1:
操作条件面板
注意，FLUENT中默认的重力加速度为零
3.假如使用不可压理想气体定律，要在操作条件面板中检查操作压力的数值（非零值）。
4.下面的选项取决于你是否使用Boussinesq近似：
●假如不使用Boussinesq模型，输入如下：
1.必要的话在操作条件面板中激活操作密度选项，然后指定操作密度，详细设置能够参阅定义操作密度一节。
2.定义流体密度为温度的函数，具体能够参阅使用温度相关函数和密度定义属性一节。
菜单：
Define/Materials...。
●假如使用Boussinesq模型，输入如下：
1.在操作条件面板中指定操作温度（Boussinesq模型一节中方程1的T_0）
2.选择Boussines方法来运算在使用材料面板中的密度（具体能够参阅使用温度相关函数和密度定义属性一节）。
3.依旧在材料面板中，设定热扩散系数并指定常数密度。
注意：
假如模型包括多种材料，关于每一个材料你都能够选择是否使用Boussinesq模型。
因此你能够对某些材料使用Boussinesq模型其它的能够不使用。
关于每一个材料的设定步骤都和上面所介绍的一样。
5.在压力入口和出口边界处的你所输入的边界压力是重新定义的压力，该压力由操作密度的定义一节中的方程3给出。
一样说来，假如没有外部强加的压力梯度，FLUENT模型在入口和出口边界处的压力p^'应该是相等的。
菜单：
Define/BoundaryConditions...。
6.在解操纵面板中，选择加权的体积力或者二阶方法作为压力的离散方法。
菜单：
Solve/Controls/Solution...。
你需要在近壁面增加单元以解决边界层问题。
假如你使用四边形或六面体网格并使用分离解算器，举荐选择PRESTO!
作为压力的离散方法。
也能够参阅热传导运算设定所需的用户输入。
操作密度的定义
当不使用Boussinesq近似时，操作密度r_0在动量方程中显现在体积力一项中：
该种形式的体积力项遵从FLUEＮT中压力的重定义：
如此，静止流体能够保证静压平稳
变成：
因此，在所有的浮力驱动流淌中，参考密度的定义差不多上专门重要的。
在默认的情形下，FLUENT会通过对所有单元取平均来运算操作密度。
在某些算例中假如你明确指定操作密度而不是让解算器来运算密度，你可能会得到更好的结果。
比方说，假如你用压力边界条件解自然对流问题，明白你所指定的压力是方程３中的p_s^'是专门重要的。
。
即使你明白真实压力p_s，你依旧需要明白操作密度r_0，以便于从p_s确定p_s^'。
因此，你应该明确定义操作密度而不使用运算的平均值。
但不管如何你所指定的密度都应该是对平均值的描述。
在某些情形下，指定操作密度会提高解的收敛性而可不能改善实际的结果。
关于这种情形，使用近似bulk密度值作为操作密度，并保证你所选的值关于区域的特点温度是合适的。
注意：
假如你使用Boussinesq近似，就可不能使用操作密度了，因此你也不必指定它。
浮力驱动流淌的解策略
关于高瑞利数流淌，你需要考虑下面的解决方针。
除此之外，在解决其它热传导问题的处理过程中所介绍的指导原则也能够用于浮力驱动流淌。
然而，需要注意的是关于高瑞利数的某些层流流淌是没有定常解存在的。
解决高瑞利（Rayleigh）数流淌的方针
关于高瑞利数流淌（Ra>10^8），为了得到最好的结果你应该遵循下面所介绍的某一处理程序：
第一个程序使用定常状态方法：
1.开始解决时使用较低的瑞利数（如：
10^7），然后使用一阶格式运行直到收敛。
2.改变有效瑞利数，改变重力加速度的数值（如：
从９．８改为０．０９８来使瑞利数减少两个量级）．
3.使用上面的结果作为高瑞利数流淌的初始推测，然后用一阶格式开始高瑞利数流淌的运算。
4.用一阶格式获得解之后，你能够采纳高阶格式连续运算。
第二个程序使用时刻相关方法来猎取定常解[62]：
1.使用相同或较低瑞利数时得到的定常状态解开始运算。
2.估量时刻常数为[14]：
其中L和U分别是长度和速度。
使用时刻步长Dt：
假如使用更大的时刻步长Dt可能会导致发散。
mp
3.当频率ft=0.05--0.09的振动衰减之后，就达到了定常状态。
注意，t是方程１中估量的时刻常数，ｆ是单位为Ｈｚ的振动频率。
一样说来，要达到定常状态一样要进行５０００个时刻步。
注意：
除非使用Boussinesq近似，否则非定常方法不能用于封闭区域。
它总是用于具有入口和出口的区域。
浮力驱动流淌的后处理
浮力驱动流淌的后处理报告和其它热传导运算的报告一样。
详情请参阅热传导的报告一节
周期性流淌和热传导
周期流是指流淌和热的解具有周期性重复的特点。
周期性流淌分两种：
一种是在周期性平面内没有压降的周期流；第二种是流向周期流。
本解讨论流向周期流以及周期性热传导，关于没有压降的周期流请参阅周期性边界条件一节。
引言
FLUENT提供流向周期流的运算。
这种流淌具有广泛的应用，如热交换管道以及通过水箱的管流。
在这些流淌模式中，几何外形沿流淌方向上具有重复性的特点，从而导致了周期性完全进展的流淌。
这些周期性条件在足够的入口长度后就会形成，具体与雷诺数和几何外形有关。
流向周期流会在足够长度L之后形成，在流向的每一个重复流淌模式之间会有一个一定的压降。
下图确实是一例。
Figure1:
二维热交换几何外形的周期性流淌例子
●交换的几何外形
●当壁面是常温或者热流一定时会产生周期性热条件。
在这一问题中，温度场是周期性进展的。
至于周期性流淌，分析这类问题能够将数值模型限制为单一模块
升级会员 