排气扇边界条件

排气扇边界条件用于模拟具有指定压力跳跃和周围（流出）环境压力的外部排气扇

排气扇边界条件的输入

排气扇边界条件需要输入：
l       静压
l       回流条件
l       总温即驻点温度（用于能量计算）。
l       湍流参数（对于湍流计算）
l       化学组分质量百分数（对于组分计算）。
l       混合分数和变化（对于PDE燃烧计算）。
l       发展变量（对于预混和燃烧计算）。
l       二级相的体积分数(对于多相流计算)
l       辐射参数(对于P-1模型、DTRM或者DO模型的计算)
l       离散相边界条件（对于离散相计算）
l       压力跳跃

上面的所有值都由排气扇面板输入（Figure 1），它是从边界条件打开的（见设定边界条件
一节）。

前四项的指定方法和压力出口边界的方法相同。详情请参阅压力出口边界的输入一节。压力
跳跃的指定在指定压力跳跃一节中描述。

Figure 1: The Exhaust Fan Panel

指定压力跳跃
FLUENT中模拟了排气扇，排气扇被假定为无限薄，并且通过排气扇具有不连续的压力升高，
它是垂直于排气扇的当地流体速度的函数。你可以定义通过排气扇的压力跳跃为常量、多项
式、分段线性函数或者分段多项式函数。定义这些函数的面板和定义温度相关属性的面板相
同，详情请参阅使用温度相关函数定义属性一节。

模拟排气扇必须小心谨慎，要保证通过排气扇向前的流动压力有所升高。在回流算例中，排
气扇被看成具有同一损失系数的进气口。

壁面边界条件

壁面边界条件用于限制流体和固体区域。在粘性流动中，壁面处默认为非滑移边界条件，但
是你也可以根据壁面边界区域的平动或者转动来指定切向速度分量，或者通过指定剪切来模
拟滑移壁面(你也可以在FLUENT中用对称边界类型来模拟滑移壁面，但是使用对称边界就需
要在所有的方程中应用对称条件。详情请参阅对称边界条件一节)。
在当地流场的详细资料基础上可以计算出流体和壁面之间的剪应力和热传导。

壁面边界的输入
概述
壁面边界条件需要输入下列信息：
l       热边界条件（对于热传导计算）
l       速度边界条件（对于移动或旋转壁面）
l       剪切（对于滑移壁面，此项可选可不选）
l       壁面粗糙程度（对于湍流，此项可选可不选）
l       组分边界条件（对于组分计算）
l       化学反应边界条件（对于壁面反应）
l       辐射边界条件(对于P-1模型、DTRM或者DO模型的计算)
l       离散相边界条件（对于离散相计算）

在壁面处定义热边界条件
如果你在解能量方程，你就需要在壁面边界处定义热边界条件。在FLUENT中有五种类型的热
边界条件：
l       固定热流量
l       固定温度
l       对流热传导
l       外部辐射热传导
l       外部辐射热传导和对流热传导的结合
如果壁面区域是双边壁面（在两个区域之间形成界面的壁面，如共轭热传导问题中的流/固
界面）就可以得到这些热条件的子集，但是你也可以选择壁面的两边是否耦合。详情请参阅
在壁面处定义热边界条件。
下面各节介绍了每一类型的热条件的输入。如果壁面具有非零厚度，你还应该设定壁面处薄
壁面热阻和热生成的相关参数，详情请参阅在壁面处定义热边界条件。

热边界条件由壁面面板输入（Figure 1），它是从边界条件打开的（见设定边界条件一节）
。

Figure 1:壁面面板
对于固定热流量条件，在热条件选项中选择热流量。然后你就可以在热流量框中设定壁面处
热流量的适当数值。设定零热流量条件就定义了绝热壁，这是壁面的默认条件。
选择固定温度条件，在壁面面板中的热条件选项中选择温度选项。你需要指定壁面表面的温
度。壁面的热传导可以用温度边界条件一节中的方程1或3来计算。
对于对流热传导壁面，在热条件中选择对流。输入热传导系数以及自有流温度，FLUENT就会
用对流热传导边界条件中的方程1来计算壁面的热传导。
如果你所模拟的是从外界而来的辐射热传导，你可以在壁面面板中xx辐射选项，然后设定
外部发射率以及外部辐射温度。
如果选择混合选项，你就可以选择对流和辐射结合的热条件。对于这种条件，你需要设定热
传导系数、自由流温度、外部发射率以及外部辐射温度。
默认情况下壁面厚度为零。然而你可以结合任何的热条件来模拟两个区域之间材料的薄层。
例如：你可以模拟两个流体区域之间的薄金属片的影响，固体区域上的薄层或者两个固体区
域之间的接触阻力。FLUENT会解一维热传导方程来计算壁面所提供的热阻以及壁面内部的热
生成。
在热传导计算中要包括这些影响，你就需要指定材料的类型，壁面的厚度以及壁面的热生成
速度。在材料名字下拉菜单中选择材料类型，然后在壁面厚度框中指定厚度。壁面的热阻为
D x/k，其中k是壁面材料的热传导系数，D x是壁面厚度。你所设定的热边界条件将在薄壁
面的外部指定，如图2所示，其中T_b壁面处所指定的固定温度。

Figure 2: 热条件被指定在薄壁面的外侧
在热生成速度框中指定壁面内部热生成速度。这一选项是非常有用的，比方说，模拟已知电
能分布的印刷电路板。
如果壁面区域的每一边是流体或者固体区域。当你具有这类壁面区域的网格读入到FLUENT，
一个阴影区域会自动产生，以便于壁面的每一边都是清楚的壁面区域。在壁面区域面板中，
阴影区域的名字将在阴影表面区域框中显示出来。你可以选择在每一个区域指定不同的热条
件或者将两个区域耦合：
l       要耦合壁面的两个边，在热条件选项中选择耦合选项(只有壁面是双边时这一选项
才会出现在壁面面板中)。不需要输入任何附加的热边界信息，因为解算器会直接从相邻单
元的解中计算出热传导。然而你可以指定材料类型、壁面厚度以及热生成速度来计算壁面热
阻，详情请参阅壁面处热边界条件的定义一节。注意，你所设定的壁面每一边的阻抗参数会
自动分配给它的阴影壁面区域。指定壁面内的热生成速度是很有用的，比如，模拟已知电能
分布但是不知道热流量或者壁面温度的印刷电路板。
l       要解耦壁面的两个边，并为每一个边指定不同的热条件，在热条件类型中选择温度
或者热流作为热条件类型(对于双边壁面，不应用对流和热辐射)。壁面和它的阴影之间的关
系会被保留，以便于你在以后可以再次耦合它们。你需要设定所选的热条件的相关参数，前
面对这方面的内容已经叙述过了不再重复。两个非耦合壁面具有不同的厚度，并且相互之间
有效地绝缘。如果对于非耦合壁面指定非零厚度的壁面，你所设定的热边界条件就会在两个
薄壁的外边的那个边指定，如图3所示，其中T_b1和T_b2分别是两个壁面的温度或者热流量
。k_w1和k_w2时耦合薄壁面的热传导率。注意图3中两个壁面之间的缺口并不是模型的一部
分，它只是在图形中用来表明每一个非耦合壁面的热边界条件在哪里应用。


Figure 3: 热条件在非耦合薄壁的外边指定

对移动壁面定义速度条件

如果你希望在计算中包括壁面的切向运动，你就需要定义平动或者转动速度。壁面速度条件
在壁面面板的运动部分输入，在这里你可以xx面板底部的移动壁面选项来显示和编辑，此
时壁面面板会扩大显示为下图：

Figure 1: 移动壁面的壁面面板
如果邻近壁面的单元区域是移动的，（比如你使用移动参考系或者滑动网格）你可以xx相
对邻近单元区域选项来选择指定的相对移动区域的移动速度。如果指定相对速度，那么相对
速度为零意味着在相对坐标系中壁面是静止的，因此在{jd1}坐标系中以相对于邻近单元的速
度运行。如果选择{jd1}速度（xx{jd1}选项），速度为零就意味着避免在{jd1}坐标系中是静
止的，而且以相对于邻近单元的速度以动，但是在相对坐标系中方向相反。
如果你使用一个或多个移动参考系、滑动网格或者混合平面，并且你希望壁面固定在移动参
考系上。推荐你指定相对速度（默认）而不是{jd1}速度。然后，如果你修改邻近单元区域的
速度，就像你指定{jd1}速度一样，你就不需要对壁面速度做任何改变。
注意：如果邻近单元不是移动的那么它和相对选项是等同的。
对于包括线性，壁面边界是平动的问题（如以移动带作为壁面的矩形导管），你可以xx平
动选项，并指定壁面速度和方向(X,Y,Z矢量)。作为默认值，通过指定平动速度为零，壁面
移动是未被xx的。
对于包括转动壁面运动的问题，你可以xx转动选项，并对指定的旋转轴定义旋转速度。要
定义轴，请设定旋转轴方向和和旋转轴原点。这一轴和邻近单元区域所使用的旋转轴是无关
的，而且和其它的壁面旋转轴无关。对于三维问题旋转轴是通过指定坐标原点的矢量，它平
行于在旋转轴方向框中指定的从(0,0,0)到(X,Y,Z)的矢量。对于二维问题，你只需要指定旋
转轴起点，旋转轴是通过指定点的z向矢量。对于二维轴对称问题，你不必定义旋转轴：通
常是绕x轴旋转，起点为(0,0)。

需要注意的是，只有在壁面限制表面的旋转时，模拟切向旋转运动才是正确的(比如圆环或
者圆柱)。还要注意只有对静止参考系内的壁面才能指定旋转运动。

如定义壁面处热边界条件所讨论的，当你读入具有双边壁面的网格时（它在流/固区域形成
界面），会自动形成阴影区域来区分壁面区域的每一边。对于双边壁面，壁面和阴影区域可
能指定不同的运动，而不管它们耦合与否。然而需要注意的是，你不能指定邻近固体区域的
壁面（或阴影）的运动。

模拟滑移壁面

作为默认，无粘流动的壁面是非滑移条件，但是在FLUENT中，你可以指定零或非零剪切来模
拟滑移壁面。要指定剪切，在壁面面板中选择指定剪切应力项（见下图），然后你可以在剪
切应力项中输入剪切的x, y, 和z分量指定剪切应力选项不是用壁面函数。

Figure 1: 滑移壁面的壁面面板

在湍流壁面限制的流动中模拟壁面粗糙度的影响

流过粗糙表面的流体会有各种各样的情况。比如流过机翼表面、船体、涡轮机、换热器以及
管系统的流动，还有具有各种粗糙度的地面上的大气边界层。壁面粗糙度影响了壁面处的阻
力、热传导和质量输运。

如果你是在模拟具有壁面限制的湍流流动，壁面粗糙度的影响是很大的，你可以通过修改壁
面定律的粗糙度来考虑避免粗糙度影响。
粗糙管和隧道的实验表明了当用半对数规则画图时，近粗糙壁面的平均速度分布具有相同的
坡度(1/k)但是具有不同的截止点(在对数定律中附加了常数B)。对于粗糙壁面，平均速度的
壁面定律具有的形式为：

其中u^* = C_m^1/4k^1/2； 是粗糙度函数，它衡量了由于粗糙影响而导致的截止点的转移
。一般说来， 依赖于粗糙的类型(相同的沙子、铆钉、螺纹、肋、铁丝网等)和尺寸。对于
各种类型的粗糙情况没有统一而有效的公式。然而，对于沙粒粗糙情况和各种类型的统一粗
糙单元，人们发现 和无量纲高度K_s ^+ = r K_s u^*/m具有很好的相关性，其中K_s 是物
理粗糙高度u^* = C_m^1/4k^1/2。实验数据分析表明粗糙函数 并不是K_s^+的单值函数，而
是依赖于K_s^+的值有不同的形式。观察表明有三种不同的类型：
l       液体动力光滑(K_s^+ < 3 ~ 5)
l       过渡区(3 ~ 5 < K_s^+ < 70 ~ 90)
l       xx粗糙(K_s^+ > 70 ~ 90)
根据上述数据，在光滑区域内粗糙度的影响可以忽略，但是在过渡区域就越来越重要了，在
xx粗糙区域具有xx的影响。

在FLUENT中，整个粗糙区域分为三个区域。粗糙函数 的计算源于Nikuradse's数据[27]基础
上的由Cebeci和Bradshaw提出的公式：
对于液体动力光滑区域(K_s^+ < 2.25)：

对于过渡区(2.25 < K_s^+ < 90)：

其中C_K_s为粗糙常数，依赖于粗糙的类型。
在xx粗糙区域(K_s^+ > 90)：

在解算器中，给定粗糙参数之后，粗糙函数 (K_s^+)用相应的公式计算出来。方程1中的修
改之后的壁面定律被用于估计壁面处的剪应力以及其它的对于平均温度和湍流量的壁面函数
。
要模拟壁面粗糙的影响，你必须指定两个参数：粗糙高度K_s和粗糙常数C_K_s。默认的粗糙
高度为零，这符合光滑壁面。对于产生影响的粗糙度，你必须指定非零的K_s。对于同沙粒
粗糙情况，沙粒的高度可以简单的被看作K_s。然而，对于非同一沙粒平均直径(D_50)应该
是最有意义的粗糙高度。对于其它类型的粗糙情况，需要用同等意义上的沙粒粗糙高度K_s
。
适当的粗糙常数(C_K_s)主要由给定的粗糙情况决定。默认的粗糙常数(C_K_s = 0.5)是用来
满足在使用k-e湍流模型时，它可以在具有同一沙粒粗糙的充满流体的管中再现Nikuradse'
s阻力数据。当你模拟和同一沙粒粗糙不同的情况时，你就需要调解粗糙常数了。例如，有
些实验数据表明，对于非同一沙粒、肋和铁丝网，粗糙常数(C_K_s =  0.5 ~ 1.0)具有更高
的值。不幸的是，对于任意类型的粗糙情况还没有一个清楚的选择粗糙常数C_K_s的指导方
针。
需要注意的是，要求邻近壁面单元应该小于粗糙高度并不是物理意义上的问题。对于{zh0}的
结果来说，要保证从壁面到质心的距离要比K_s大。

定义壁面的组分边界条件

FLUENT默认所有的组分在壁面处具有零梯度条件（除了参加表面化学反应的组分），但是可
以指定壁面处的组分质量分数。也就是如同在入口处指定的Dirichlet边界条件，也可以用
于壁面。
如果你希望保留默认的零梯度条件，你就不必输入任何东西了。如果你希望指定壁面处的组
分质量分数，步骤如下：
1.      在壁面面板的组分边界条件中，选择组分名字右边的下拉列表指定的质量分数（而
不是零梯度），此时面板会扩展为包含组分质量分数的对话框。

Figure 1: 组分边界条件输入的壁面面板
2.      输入相应的组分质量分数。
每一组分的边界条件类型是分别指定的，所以对于不同的组分你可以采用不同的方法。
注意：如果在湍流流动中你使用组分的Dirichle条件，FLUENT就不会是用壁面函数来计算壁
面处的组分扩散流量。

定义壁面的反应边界条件

如果你在组分模型面板中xx了表面反应的模拟，你就可以表明在壁面处表面反应是否被激
活。xx或关闭表面反应，壁面面板就会相应地打开或关闭表面反应选项。
注意：组分在壁面处是假定为零梯度条件的，它不参加任何表面反应。

定义壁面的辐射边界条件

如果你打算使用P-1辐射模型、DTRM或者DO模型，你就需要设定壁面的（内部）发散率以及
（可选）黑体温度。详情请参阅设定边界条件一节(Rosseland不需要任何边界条件的输入，
因为FLUENT假定发射率为1，如果你使用DO模型你也要定义壁面为漫反射、镜面反射或者半
透明，详情请参阅设定辐射边界条件)

定义壁面的离散相边界条件

如果你是在模拟粒子的离散相，你就可以在壁面处设定粒子轨道详情请参阅离散向模型的边
界设定。

壁面边界的默认设定

默认热边界条件为固定的热流为零，壁面默认为不移动。

壁面处的剪应力计算程序

对于非滑移壁面条件，FLUENT使用邻近壁面或者流体边界的流动性质来预测壁面处流体的剪
应力。在层流流动制，这一计算简单地依赖于壁面处的速度梯度，在湍流流动中则使用壁面
限制湍流流动的近壁面处理方法。
对于指定剪切的壁面， FLUENT会在边界处计算切向速度。
如果是无粘流动，所有的壁面都使用滑移条件，所以它们是无摩擦的而且对邻近流体单元不
施加剪应力。

层流中的剪应力计算
在层流流动中壁面剪应力和法向速度梯度的关系为：

当壁面处的速度梯度很大时，你必须保证网格足够精细，这样才能解出边界层的xx结果。
层流流动中近壁面节点放置的指导方针在节点密度和节点束中介绍。

湍流中的剪应力计算

湍流流动的壁面处理，在壁面限制的湍流流动的近壁面处理一节中叙述。

壁面边界的热传导计算

温度边界条件

当在壁面处应用固定温度条件，从流体单元到壁面的热传导，由下式计算：

其中：
h_f     =流体边界当地热传导系数
T_w     =壁面表面温度
T_f     =当地流体温度
q^"=壁面处传来的对流热流量
q^"_rad=辐射热流量
注意：流体边界热传导系数是基于当地流场条件计算得来的（比如说湍流层次、温度以及速
度轮廓），请参阅流体边界热传导计算一节的方程1，以及标准壁面函数9。
从固体单元到壁面边界的热传导公式为：

其中：
k_s     =固体的热传导率
T_s     =当地固体温度
D n     =壁面表面和固体单元中心的距离。

热流边界条件

当你在壁面处定义热流边界条件时，你需要在壁面表面指定热流量。FLUENT使用温度边界条
件中的方程1，然后你就可以输入热流量来确定邻近流体单元的壁面表面温度：

其中，流体边界热传导系数已经在温度边界条件中叙述了，它是基于当地流场条件计算得到
的。当壁面和固体区域交界时，壁面表面的温度为：

上述两式的变量请参阅温度边界条件一节。

对流热传导边界条件

当你在壁面处指定对流热传导系数作为边界条件时，FLUENT使用你所输入的外部热传导系数
以及外部热沉（heat sink）温度来计算到壁面的热流量：

其中：
h_ext   =你所定义的外部热传导系数
T_ext   =你所定义的外部热沉温度
q^"_rad =辐射热流量
上述方程假定壁面零厚度。

外部辐射边界条件
当使用外部辐射条件时，流入壁面的热流量为：

其中：
e_ext=你所定义的外部壁面表面的发射率
s=Stefan-Boltzmann常数
T_w     =壁面的表面温度
T_?=区域外部的温度的辐射源或者消失（sink）处
q^"_rad=从内部去向壁面辐射的热流量
Equation 1假定壁面厚度为零。

外部对流和辐射结合的边界条件
当你选择组合的外部热传导方程条件时，到壁面的热流量为：

其中的变量已经在对流热传导边界条件和外部辐射边界条件中定义了。Equation 1假定壁面
厚度为零。

流动边界热传导系数的计算

在层流流动中，壁面处流体边界热传导是用应用于壁面的Fourier定律计算得到的，FLUENT
使用它的离散格式为：

其中n是垂直于壁面的当地坐标。
对于湍流流动，FLUENT对于从热和动量迁移中类比得到的温度使用壁面定律[93]。详细内容
请参阅标准壁面函数。

对称边界条件

对称边界条件用于所计算的物理外形以及所期望的流动/热解具有镜像对称的特征的情况中
。也可以用它们来模拟粘性流动的滑移壁面。本节描述了对称平面内流动的处理，并提供了
一些使用对称边界的例子。在对称边界条件中你不需要定义任何边界条件，但是你必须谨慎
地定义对称边界的位置。
在对称外形的中线处，你应该使用轴边界类型而不是对称边界类型，如轴边界条件一节中的
的图1，详细内容请参阅轴边界条件。

对称边界的计算程序

FLUENT假定所有量通过对称边界的流量为零。经过对称平面的对流流量为零，因此对称边界
的法向速度为零。通过对称平面没有扩散流量：因此所有流动变量的法向梯度在对称平面内
为零。因此对称边界条件可以总结如下：
l       对称平面内法向速度为零
l       对称平面内所有变量的法向梯度为零
如上所述，对称的定义要求这些条件决定流过对称平面的流量为零。因为对称边界的剪应力
为零，所以在粘性流动计算中它也可以用滑移壁面来解释。

对称边界的例子

对称边界用于减少计算模拟的范围，它只需要模拟所有物理系统的一个对称子集。下面两个
图是通过该种方法使用对称边界的例子。

Figure 1:使用对称边界模拟三维管道的四分之一

Figure 2: 使用对称边界模拟圆形截面的四分之一
下面的图则是误用对称平面的两个例子，在这两个例子中，虽然几何外形是对称的，但是流
动本身却不符合对称边界条件的要求。在{dy}个例子中浮力产生了非对称流动。在第二个例
子中，流动中的涡流产生了一个垂直于应该是对称平面的流动。。需要注意的是，这两个粒
子都要使用旋转周期性边界（请参阅周期性边界一节的图一）

Figure 3: 对称的误用

周期性边界条件

周期性边界条件用来解决，物理模型和所期待的流动的流动/热解具有周期性重复的特点。F
LUENT提供了两种类型的周期性边界条件。{dy}种类型不允许通过周期性平面具有压降（对
于FLUENT4用户来说：这一类型的周期性边界是指FLUENT4中的圆柱形边界）。第二种类型允
许通过平移周期性边界具有压降，它是你能够模拟xx发展的周期性流动（在FLUENT4中是
周期性边界）。
本节讨论了无压降的周期性边界条件。在周期性流动和热传导一节中，xx发展的周期性模
拟能力得到了详尽的描述。

周期性边界的例子

周期性边界条件用于模拟通过计算模型内的两个相反平面的流动是相同的情况。下图是周期
性边界条件的典型应用。在这些例子中，通过周期性平面进入计算模型的流动和通过相反的
周期性平面流出流场的流动是相同的。正如这些例子所示，周期性平面通常是成对使用的。

Figure 1: 在圆柱容器中使用周期性边界定义涡流

周期性边界的输入

对于没有任何压降的周期性边界，你只需要输入一个东西，那就是你的所模拟的几何外形是
旋转性周期还是平移性周期。（对于有周期性压降的周期流还要输入其它的东西，请参阅周
期性流动和热传导一节。）

旋转性周期边界是指关于旋转对称几何外形中线形成了一个包括的角度。本节中的图一就是
旋转性周期。平移性周期边界是指在直线几何外形内形成周期性边界。下面两图是平移性周
期边界：

Figure 1: 物理区域

Figure 2: 所模拟的区域
对于周期性边界，你需要在周期性面板（下图）中指定平移性边界还是旋转性边界，该面板
是从设定边界条件菜单中打开的。

Figure 3: 周期性面板

(对于耦合解算器，周期性面板中将会有附加的选项，这一选项允许你指定压力跳跃，详细
内容请参阅周期性流动和热传导一节。)
如果区域是旋转性区域，请选择旋转性区域类型。如果是平移性就选择平移性区域类型。对
于旋转性区域，解算器会自动计算通过周期性区域的旋转角度。旋转轴是为邻近单元指定的
旋转轴。
注意：对于使用旋转周期性边界来说，你不必指定邻近单元区域为移动的。例如，你能够使
用具有管的平切片的非旋转坐标系来模拟三维管流，管的切片需要具有旋转性周期。
你可以使用Grid/Check菜单选项（参阅检查网格一节）来计算和显示周期性边界所有表面的
旋转角度的{zd0}值、最小值和平均值。如果{zd0}值、最小值和平均值之间的差别可以忽略
，那么网格有一个问题：对于指定轴来说网格几何外形不是周期性的。

周期性边界的默认设定

默认为平移周期性边界条件

周期性边界的计算程序

FLUENT在周期性边界处理流动就像反向周期性平面是和前面的周期性边界直接相邻一样，因
此，当计算流过邻近流体单元的周期性边界时，就会使用与反向周期性平面相邻的流体单元
的流动条件。

轴边界的计算程序

轴边界条件

轴边界类型必须使用在对称几何外形的中线处（见下图）。它也可以用在圆柱两极的四边形
和六面体网格的中线上（比如：像FLUENT4之类的结构网格生成代码所产生的网格）。在轴
边界处，你不必定义任何边界条件。

Figure 1: 在轴对称几何外形的中线处轴边界条件的使用

轴边界的计算程序
要确定轴上特定点的适当物理值，FLUENT使用邻近单元中的单元值。