-
边界条件
定义边界条件概述
边界条件包括流动
变量和热变量在边界处的值。
它是
FLUENT
分析得很关键的一部分,
设定边界条件必须小心谨慎。
边界条件的分类:进出口边界条件:压力、速度、质量进口、进风口、进气扇、压力出<
/p>
口、压力远场边界条件、质量出口、通风口、排气扇;壁面、
re
peating, and pole boundaries:
壁面,对称,周期,轴
;内部单元区域:流体、固体
(
多孔是一种流动区域类型
)
;内部表
面边界:风扇、散热器、多孔跳
跃、壁面、内部。
(
内部表面边界条件定义在单元表面,这
p>
意味着它们没有有限厚度,
并提供了流场性质的每一步的变化。
p>
这些边界条件用来补充描述
排气扇、细孔薄膜以及散热器的物理模型
。内部表面区域的内部类型不需要你输入任何东
西。
)
下面一节将详细介绍上面所叙述边界条件,
< br>并详细介绍了它们的设定方法以及设定的具
体合适条件。
周期性边界条件在本章中介绍,
模拟完全发展的周期性流动将在周期性流动和
热传导一章中介绍。
使用边界条件面板
边界条件
(Figure 1)
对于特
定边界允许你改变边界条件区域类型,
并且打开其他的面板以
设
定每一区域的边界条件参数
菜单:
Define/Boundary
Conditions...
Figure 1:
边界条件面板
改变边界区域类型
设定任何边界条件
之前,
必须检查所有边界区域的区域类型,
如有必要就作适当的
修改。
比方说:
如果你的网格是压力入口,
但是你想要使用速度入口,
你就要把压力入口改为速度
入口
之后再设定。
改变类型的步骤如下:
:
1.
在区域下拉列表中选定所要修改的区域
2.
在类型列表中选择正确的区域类型
3.
当问题提示菜单出现时,点击确认
确认改变之后,
区域类型将会改变,
名字也将自动改变
(
如果初始名字时缺省的请参阅
边界条件区域名字一节
)
,
设定区域边界条件的面板也将自动打开。
< br>!注意:这个方法不能用于改变周期性类型,
因为该边界类型已经存在了附加限制
。创
建边界条件一节解释了如何创建和分开周期性区域。
需要注
意的是,
只能在图一中每一个类
别中改变边界类型
(
注意:双边区域表面是分离的不同单元区域
.)
Figure 1:
区域类型的分类列表
设定边界条件
在
FLUENT
中,边界条件和区域有关而与个别表面或者单元无关。如果要结合具有
相
同边界条件的两个或更多区域请参阅合并区域一节。
设定每一特定区域的边界条件,请遵循下面的步骤:
1.
在边界条件区域的下拉列表中选
择区域。
2.
点击
Set...
p>
按钮。
或者,
1.
在区域下拉列表
中选择区域。
2.<
/p>
在类型列表中点击所要选择的类型。
或者在区域列表中双击所需区
域
.
,
选择边界条件
< br>区域将会打开,并且你可以指定适当的边界条件
在图像显示方面选择边界区域
在边界
条件中不论你合适需要选择区域,
你都能用鼠标在图形窗口选择适当的区域。
如
果你是第一次设定问题这一功能尤其有用,
如果
你有两个或者更多的具有相同类型的区域而
且你想要确定区域的标号
(也就是画出哪一区域是哪个)
这一功能也很有用。
要使用
该功能
请按下述步骤做:
1.
用网格显示面板显示网格。
2.
用鼠标指针(
默认是鼠标右键——参阅控制鼠标键函数
以改变鼠标键的功能)
在图形窗口中点击边界区域。
在图形显示中选择的区域将会自动被选
入在边界条件面板中的区域列表中,它的名字和编号也会自动在控制窗口中显示
改变边界条件名字
每一边界的名字是
它的类型加标号数
(比如
pressure-inlet-7<
/p>
)
。
在某些情况下你可能想要
对边界区域分配更多的描述名。
如果你有两个压力入口区域,
比方说,
你可能想重名名它们
为
small-inlet
和
large-inlet
。
(
改变边界的名字不会改变相应的类型
p>
)
重名名区域,遵循如下步骤:
1.
在边界条件的区域下拉列表选择所要重名名的区域。
2.
点击
S
et...
打开所选区域的面板。
3.
在区域名字中输入新的名字
4.
点击
O
K
按钮。
注意:
如果你指定区域的新名字然后改变它的类型,
你所改的名字将会被保留,
如果区域名
字是类型加标号,名字将会自动改变。
边界条件的非一致输入
每一类型的边
界区域的大多数条件定义为轮廓函数而不是常值。
你可以使用外部产生的
边界轮廓文件的轮廓,或者用自定义函数
(UDF)
来
创建。具体情况清参阅相关内容
流动入口和出口
< br>FLUENT
有很多的边界条件允许流动进入或者流出解域。下面一节描述了每一
种边界
条件的类型的使用以及所需要的信息,
这样就帮助你适当
的选择边界条件。
下面还提供了湍
流参数的入口值的确定方法。
使用流动边界条件
下面对流动边界条件的使用作一概述
对于流动的出入口,
FLUENT
提供
了十种边界单元类型:速度入口、压力入口、质量
入口、压力出口、压力远场、质量出口
,进风口,进气扇,出风口以及排气扇。
< br>下面是
FLUENT
中的进出口边界条件选项:
?
速度入口边界条件用于定义流动入口边界的速度和标量
?
压力入口边界条件用来定义流动入口边界的总压和其它标量。
?
质量流动入口边界条件用于可压流
规定入口的质量流速。
在不可压流中不必指定入口的
质量流,因
为当密度是常数时,速度入口边界条件就确定了质量流条件。
?
压力出口边界条件用于定义流动出
口的静压(在回流中还包括其它的标量)
。当出现回
流时,使用
压力出口边界条件来代替质量出口条件常常有更好的收敛速度。
?
压力远场条件用于模拟无穷远处的
自由可压流动,
该流动的自由流马赫数以及静态条件
已经指定了
。这一边界类型只用于可压流。
?
质量出口边界条件用于在解决流动问题之前,
所模拟的流动出口
的流速和压力的详细情
况还未知的情况。
在流动出口是完全发展
的时候这一条件是适合的,
这是因为质量出口
边界条件假定出了
压力之外的所有流动变量正法向梯度为零。
对于可压流计算,
这
一条
件是不适合的。
?
进风口边界条件用于模拟具有指定
的损失系数,
流动方向以及周围
(入口)
环境总压和
总温的进风口。
?
进气扇边界条件用于模拟外部进气
扇,
它具有指定的压力跳跃,
流动方向以及周围
(进
口)总压和总温。
?
通风口边界条件用于模拟通风口,
它具有指定的损失系数以及周围环境
(排放处)
的静
压和静温。
?
排气扇边界条件用于模拟外部排气
扇,它具有指定的压力跳跃以及周围环境(排放处)
的静压。
决定湍流参数
在入口、出口或远场边界流入流域的流动,
< br>FLUENT
需要指定输运标量的值。本节描
述了对于特
定模型需要哪些量,
并且该如何指定它们。
也为确定流入边界值
最为合适的方法
提供了指导方针。
使用轮廓指定湍流参量
在入口处要准确的描述边界层和完全发展的湍流流动,
你应该通过实验数据
和经验公式
创建边界轮廓文件来完美的设定湍流量。
如果你有轮
廓的分析描述而不是数据点,
你也可以
用这个分析描述来创建边
界轮廓文件,
或者创建用户自定义函数来提供入口边界的信息。
一
旦你创建了轮廓函数,你就可以使用如下的方法:
?
Spalart-Allmara
s
模型:在湍流指定方法下拉菜单中指定湍流粘性比,并在在湍流粘性
< br>比之后的下拉菜单中选择适当的轮廓名。通过将
m_t/m
和密度与分子粘性的适当结合,
FLUENT
为修改后的湍流粘性计算边界值。
?
k-e
模
型:在湍流指定方法下拉菜单中选择
K
和
Epsilon
并在湍动能(
Turb.
Kinetic
Energy
)
p>
和湍流扩散速度
(
Turb.
Dissipation Rate
)
之后的下拉菜单中选择适
当的轮廓名。
?
< br>雷诺应力模型:
在湍流指定方法下拉菜单中选择
K
和
Epsilon
并在湍动能
(
Turb. Kinetic
Energy
)
和湍流扩散速度
(
Turb. Dissipation Rate
)
之后的下拉
菜单中选择适当的轮廓名。
在湍流指定方法下拉菜单中选择雷诺应力部分,
并在每一个单独的雷诺应力部分之后的
下拉菜单中选择适当的轮廓名。
湍流量的统一说明
在某些情况下流动流入开始时,
将边界处的所有湍流量指定为统一值是适当的。
比如说,
在进入管道的流体,远场边界,甚至完全发展的管流中
,湍流量的精确轮廓是未知的。
在大多数湍流流动中,湍流的
更高层次产生于边界层而不是流动边界进入流域的地方,
因此这就导致了计算结果对流入
边界值相对来说不敏感。
然而必须注意的是要保证边界值不
是非
物理边界。
非物理边界会导致你的解不准确或者不收敛。
对于外
部流来说这一特点尤其
突出,如果自由流的有效粘性系数具有非物理性的大值,边界层就
会找不到了。
你可以在使用轮廓指定湍流量一节中描述的湍流
指定方法,来输入同一数值取代轮廓。
你也可以选择用更为方便的量来指定湍流量,
p>
如湍流强度,
湍流粘性比,
水力直径以及湍
流
特征尺度,下面将会对这些内容作一详细叙述。
湍流强度
I
定义为相对于平均速度
u_avg
的脉动速度
< br>u^'
的均方根。
小于或等于
1%
的湍流强度通常被认为低强度湍
流,大于
10%
被认为是高强度湍流。从
外界,测量数据的入口边界,你可以很好的估计湍流强度。例如:如果你模拟风洞试验,自
由流的湍流强度通常可以从风洞指标中得到。
在现代低湍流风洞中自由流湍流强度通
常低到
0.05%
。
.
对于内部流动,
入口的湍流强度完全
依赖于上游流动的历史,
如果上游流动没有完全发
展或者没有被
扰动,
你就可以使用低湍流强度。
如果流动完全发展,
湍流强度可能就达到了
百分之几。完全发展的管流的核心的湍流强度可以
用下面的经验公式计算:
I
?
u
?
?
0
.
16
Re
D
H
u
avg
?
?
?
1
8
例如,在雷诺数为
50000
是湍流强度为
4%
湍流尺度
l
p>
是和携带湍流能量的大涡的尺度有关的物理量。在完全发展的管流中,
l
被
管道的尺寸所限制,
因为大涡不
能大于管道的尺寸。
L
和管的物理尺寸之间的计算关系如下:<
/p>
l
?
0
.
07
L
其中
L
为管道的相关
尺寸。因子
0.07
是基于完全发展湍流流动混合长度的最大值
的,对于
非圆形截面的管道,你可以用水力学直径取代
L
。
如果湍流的产生是由于管道中的障碍物等
特征,你最好用该特征长度作为湍流长度
L
而不是用管道尺寸。
注意:
公
式
l
?
0
.<
/p>
07
L
并不是适用于所有的情况。
它只是在大多数情况下得很好的近似。
对于特定流动,选择
L
和
l
的原则如下:
p>
?
对于完全发
展的内部流动,
选择强度和水力学直径指定方法,
并在水力学直
径流场中指
定
L=D_H
。
?
对于旋转叶片的下
游流动,
穿孔圆盘等,
选择强度和水力学直径指定方法,
并在水力学
直径流场中指定流动的特征长度为
L
?
对
于壁面限制的流动,
入口流动包含了湍流边界层。
选择湍流强度
和长度尺度方法并使
用边界层厚度
d_99
来计算湍流长度尺度
l
,
在湍流长
度尺度流场中输入
l=0.4
d_99
这个
值
湍流粘性比
m_t/m
直接与湍流雷诺数成比例
(Re_t
?k^2/(e
n))
。
Re_t
在高湍流数的边界
层,剪切层和完全发展的管流中是较大的
(100
到
1000)
。然而,在大多数外流的自由流边界
层中
m_t/m
相当的小。湍流参数的典型设定为
1 < m_t/m <10
。
要根据湍流粘性比来指定量,你可以选择湍流粘性比(对于
Spalart-A
llmaras
模型)或
者强度和粘性比(对于
k-e
模型或者
RSM
)
p>
。
推导湍流量的关系式
要获得更方便的湍流量的输运值,如:
I, L,
或者
m_t/m
,你必须求助于经验公式,下面
p>
是
FLUENT
中常用的几个有用的关系式
。要获得修改的湍流粘性,它和湍流强度
I
长度尺
度
l
有如下关系:
~
?
v
3
u
avg
Il
2
在
Spalart-Allmara
s
模型中,如果你要选择湍流强度和水力学直径来计算
l
可以从前面
的公式中获得。
湍动能
k
和湍流强度
I
之间的关系为:
k
?
3
?
u<
/p>
avg
I
?
2<
/p>
2
其中
u_a
vg
为平均流动速度
除了为
k
和
e
指定具体的值之
外,
无论你是使用湍流强度和水力学直径,
强度和长度尺
度或者强度粘性比方法,你都要使用上述公式。
如果你知道湍流长度尺度
l
你可以使用下面的关系式:
p>
?
?
C
?
3
4
k
l
3
2
< br>其中
C
?
是湍流模型中指定的经
验常数(近似为
0.09
)
,
l
的公式在前面已经讨论了。
除了为
k
和
e
制定具体的值之外,
无论你是使用湍流强度和水力学直径还是强度和长度
尺度,你都要使用上述公式。
E
的值也可以用下式计算,它与湍流粘性比
m_t/m
以及
k
有关:
k
2
?
?
t
?
?
?
?<
/p>
C
?
?
?
?
?
?
?
?
?
?
1
其中
C
?
< br>是湍流模型中指定的经验常数(近似为
0.09
)
。
除了为
k
和
e
制定具体的值之外,
无论你是使用湍流强度和水力学直径还是强度和长度
尺度,你都要使用上述公式。<
/p>
如果你是在模拟风洞条件,
在风洞中模型被安装在网格和
/
或金属网格屏下
游的测试段,
你可以用下面的公式:
?
?
?
kU
?
L
?
其中,
?
k
是你希望的在穿过流场之后
k
的衰减
(
比方说<
/p>
k
入口值的
10%),
U
?
自由流的速度
L
?
是流域内自由流的流向长度
Equation
9
是在高雷诺数各向同性湍流中观察到的幂率衰减
的线性近似。它是基于衰减湍流中
k
的精确方程
p>
U ?k/?x = - e.
如果你
用这种方法估计
e
,你也要用方程
7<
/p>
检查结果的湍流粘性比
m_t/m
,以保
证它不
是太大。
虽然这不是
FLUENT
内部使用的方法,但是你可以用它来推导
< br>e
的常数自由流值,然
后你可以用湍流指定方法下拉菜单
中选择
K
和
Epsilon
直接指定。在这种情况下,你需要使
用方程
3
p>
从
I
来计算
k
p>
。
当使用
RSM
时,如果你不在雷诺应力指定方法的下拉列表中使用雷诺应力选项,明显
的制定入口处的雷诺应力值,它们就会近似的由
k
的指
定值来决定。湍流假定为各向同性,
保证
u
i
u
j
?
0
以及
u
?
u
?
p>
?
2
k
3
(
下标
a
不求和
).
如果你在
雷诺应力指定方法下拉列表中选择
K
或者湍流强度,
FLUENT
就会使用这种方法。
对大涡模拟(
LES
)指定入口湍流
大涡模拟
模型一节中所描述的
LES
速度入口中指定的的湍流强度值,被
用于随机扰动入口
处速度场的瞬时速度。
它并不指定被模拟的湍
流量。
正如大涡模拟模型中介绍的边界条件中
所描述的,通过叠
加每个速度分量的随机扰动来计算流动入口边界处的随机成分
.
压力入口边界条件
压力入口边界条件用于定义流动入口的压力以及其它标量属性。它即可以适用于可压
流,也可以用于不可压流。
压力入口边界条件可用于压力已知但是流动速度和
/
或速率未知
的情况。
这一情况可用于很多实际问题,
比如浮力驱动的流动。
压力入
口边界条件也可用来
定义外部或无约束流的自由边界。对于流动边界条件的概述,请参阅
流动入口和出口一节。
压力入口边界条件的输入
综述
对于压力入口边界条件你需要输入如下信息
?
驻点总压
?
驻点总温
?
流动方向
?
静压
?
湍流参数(对于湍流计算)
?
辐射参数
(
对于使用
P-1
模型、
DTRM
模型或者
DO
模型
的计算
)
?
化学组分质量百分比
(
对于组分计算
)
?
混合分数和变化
< br>(
对于
PDF
燃烧计算
)
?
程序变量
p>
(
对于预混和燃烧计算
)
?
离散相边界条件
< br>(
对于离散相的计算
)
?
次要相的体积分数
(
对于多相计算
)
所有的值都在压力入口面板中输入
(Figure
1)
,该面板是从边界条件打开的。
Figure 1:
压力入口面板
压力输入和静压头
压力场
(p_s^')
和压力输入
(p_s^' or
p_0^')
包括静压头
r_0 g
x
。也就是
FLUENT
以下式定
义的压力:
p
?
s
?
?
0
gx
?
p
s
或者
?
p
?
?
p>
p
s
?
?
0
g
?
s
?
x
x
< br>这一定义允许静压头放进体积力项
(r - r_0)g
中考虑,
而且当密度一致时,
从压力计算中
排除了。因此你的压力输入不因该考虑静压的微分,压力
(p^'_s)
的报告也不会显示静压的任
何影响。有关浮力驱动流动的内容请参阅浮力驱
动流动和自然对流的信息
定义总压和总温
在压力入口面板中的
Gauge Total
Pressure field
输入总压值。
总温会在
Total Temperature
field
中设定。记住,总压值是在操作条件面板中定义的与操作压力有关的的总压值。不可
压流体的总压定义为:
p
0
?
p
s
?
?
v
2
对于可压流体为:
?
?
?
1
2
?
p
0
?
p
s
?
1<
/p>
?
M
?
2
?
?
?
?
?
?
1
?
其中:
p_0
=
总压
p_s
=
静压
M
=
马赫数
c
=
比热比
(c_p/c_v)
p>
如果模拟轴对称涡流,方程
1
中的
v
包括了旋转分量。如果相邻区域是移动的(即:如
果使用旋转参考坐标系,多重参考坐标系,混合平面或者滑移网格)
,而且你是使用
分离解
算器。那么方程
1
中的速度(或
者方程
3
中的马赫数)将是绝对的,或者相对与网格速度。
p>
这依赖于解算器面板中绝对速度公式是否激活。
对于耦合解算器,<
/p>
方程
1
中的速度
(或者方
程
3
中的马赫数)通常是在绝
对坐标系下的速度。
定义流动方向
你可以在压力入口明确的定义流动的方向,
或者定义流动垂直于边界。
如果你选择指定
方向矢量,你既可以设定笛卡尔坐标
x, y,
和
z
的分量,也可以设
p>
(
圆柱坐标的
)
半
径,切线和
轴向分量。
对于使用分离解算器计算移动区域问题,
流动方向将是绝对速度或者相对于网格
相对速度,
这取决于解算器面板中的绝对速度公式是否被激活。
对于耦合解算器,
流动方向
通常是绝对坐标系中的。
定义流动方向的步骤如下,总结请参考
Figure
1
。
1.
在方向指定下拉菜单中选择指定流动方向的方法,或者是方向矢量或者是垂直于边界。<
/p>
2.
如果你
在第一步中选择垂直于边界,
并且是在模拟轴对称涡流,
请输入
流动适当的切向
速度,如果不是模拟涡流就不需要其它的附加输入了。
< br>
3.
如果第一步中你选择指
定方向矢量,
并且你的几何外形是
3
维
的,
你就需要选择定义矢
量分量的坐标系统。在坐标系下拉菜单
中选择笛卡尔
(X, Y, Z)
坐标,柱坐标(半径,切
p>
线和轴)
,或者局部柱坐标。
?
笛卡尔坐标系是基于几何图形所使用的笛卡尔坐标系。
?
柱坐标在下面的坐标系统的基础上
使用轴、角度和切线三个分量。
?
对于包含一个单独的单元区域时,坐标系由旋转轴和在流体面板中原来的指定来定义。<
/p>
?
对于包含多重区域的问题(比如多重参考坐标或滑动网格)
,坐标系由
流体(固体)面
板中为临近入口的流体(固体)区域的旋转轴来定义。
< br>
对于上述所有柱坐标的定义,
正径向速度指向旋转轴的
外向。
正轴向速度和旋转轴矢量
的方向相同,正切向方向用右手
定则来判断。参阅下图一目了然。
Figure 1:
在二维、三维和轴对称区域的柱坐标速度分量
当地柱坐标系统允许你对特定的入口定义坐标系,
在压力入口面板中你就可以定
义该坐
标系统。如果你对于不同的旋转轴有几个入口,那么当地坐标系会很有用的。
p>
4.
如果你在第一步中指定方向矢量,用如下的方法定义矢量分量:
?
如果是二维非对称图形或者你在第
三步中选择矢量分量,
请输入适当的
X,
Y
,
和
(in 3D)
Z
分量。
?
如果是二维轴对称图形或者第三部
分选择了柱坐标,
请输入适当的半径,
角度以及切线
方向的分量。
?
p>
如果使用当地柱坐标系,
请输入适当的半径,
角度以及切线方向的分量,
并指定轴向的
X, Y
,
和
Z
向分量,以及坐标
起点的坐标。
图一就是各个坐标系统的矢量分量。
定义静压
如果入口流动是超声速的,
或者你打算用压力入口边界条件来对解进行初始化,
那么你
必须指定静压
(termed the
Supersonic/Initial Gauge
Pressure)
。
需要记住的是
这个静压和你在操作条件面板中的操作压力是相关的。
请参阅有关于压力
输入和静压头相关输入的解释。
只要流动是压声速的
,
FLUENT
会忽略
Superso
nic/Initial Gauge Pressure
,
它
是由指定的
驻点值来计算的。
如果你打算使用压力入口边界条件
来初始化解域,
Supersonic/Initial Gauge
Pressure
是与计算初始值的指定驻点压力相联系的,
< br>计算初始值的方法有各向同性关系式
(对
于可压流)或者
贝努力方程(对于不可压流)
。因此,对于压声速入口,它是在关于入口马
赫数(可压流)或者入口速度(不可压流)合理的估计之上设定的。
定义湍流参数
对于湍流计算,
有几种方法来定义湍
流参数。
至于哪种方法合适请参阅决定湍流参数一
节。湍流模型
是在“湍流模型”一章中介绍
定义辐射参数
如果你打算使用
P-1
辐射模型、
DTRM
或者
DO
模型,
你就需要设定内部发散率以及
(可
选)黑体温度。详
情请参阅设定边界条件一节
(Rosseland
不需要任何边
界条件的输入
)
。
定义组分质量百分比
如果你是用有限速度模型来模拟组分输运,
< br>你就需要设定组分质量百分比。
详情请参阅
组分边界条件
的定义。
定义
PDF/
混合分数参数
如果你用
PDF
模型模拟燃烧,你就需要设定
平均混合分数以及混合分数变化(如果你
是用两个混合分数就还包括二级平均混合分数和
二级混合分数变化)
。具体情况如第三步定
义边界条件所述。<
/p>
定义预混和燃烧边界条件
如果使用与混合燃烧模型,你就需要设定发展变量。请见发展变量的边界条件设定。
定义离散相边界条件
如果你是在模拟粒子的离散相,
你就
可以在压力入口设定粒子轨道详情请参阅离散向模
型的边界设定。
定义多相边界条件
对于多相流如果使用
VOF
,
cavitation
或者代数滑移混合模型
,你就需要指定所有二级
相的体积分数。详情请参阅
VOF
p>
模型、
cavitation
模型或者代数
滑移混合模型的边界设定。
压力入口边界条件的默认设定
p>
压力入口边界条件的默认设定如下(国际标准单位)
:
Gauge Total Pressure
0
Supersonic/Initial Gauge
Pressure
0
Total
Temperature
300
X-Component of Flow Direction
1
Y-Component of
Flow Direction
0
Z-Component of Flow Direction
0
Turb. Kinetic
Energy
1
Turb.
Dissipation Rate
1
压力入口边界处的计算程序
FLUENT
压力入口边界条件的处理可以描述为从驻点条件到入口条
件的非自由化的过
渡。
对于不可压流是通过入口边界贝努力方程
的应用来完成的。
对于可压流,
使用的是理想
< br>气体的各向同性流动关系式。
压力入口边界处的不可压流动计算
流动进入压力入口边界时,
FLUENT
使用边界条件压力,该压力是作为入口平面
p_0
的总压输入
的。
在不可压流动中,
入口总压,
静压
和速度之间有如下关系:
p
0
?
p
s
?
1
2
?
v
。
2
通过你在出口分配的速度大小和流动方向可以计算出速度的各个分量。
入口质量流速以及动
量、能量和组分的流量可以作为计算程序在速度入
口边界的大纲用来计算流动
对于不可压流,
< br>入口平面的速度既可以是常数也可以是温度或者质量分数的函数。
其中
质量分数是你输入作为入口条件的值。
在通过压力出口流出的流动,
用指定的总压作为静压
来使用。对于不可压流动来说,总温和静温相等
。
压力入口边界的可压流动计算
p>
对于可压流,
应用理想气体的各向同性关系可以在压力入口将总压,
静压和速度联系起
来。在入口处输入总压,在临近流体单元中输
入静压,有关系式如下:
?
?
p
0
p
?
?
?
1
2
?
?
p
0
?
?
1
?
M
?
p
?
?<
/p>
p
2
?
?
s
0
p
其中马赫数定
义为:
?
?
?
1
?
M<
/p>
?
v
v
?
c
?
RT
s
马赫数的定义就不详述了。需要注意的是上面的方程中出现了操作压力
p_op
这是因为边界
条件的输入是和
操作压力有关的压力。
给定
p_0^'
和
p_s^'
上面的方程就可以用于计算入口平
面流体的速度范围。
入口处的各个速度分量用方向矢量来计算。
对于可压流,
入口平面的密
度由理想气体定律来计算:
?
?
p
?
p>
s
?
p
0
p
?
?
RT
s
。
R
由压力入口边界条件定义的组分质量百分比来计算。
入口静温和总温的关系由下
式计
算:
T
0
?
?
1
2<
/p>
?
1
?
M
。
T
s
2
速度入口边界条件
速度入口边界条件用于定义流动速度以及流动入口的流动属性相关标量。
在这个边界条
件中,流动总的(驻点)的属性不是固定的,所以无论什么时候提
供流动速度描述,它们都
会增加。
这一边界条件适用于不可压流,
如果用于可压流它会导致非物理
结果,
这是因为它允许
驻点条件浮动。
你也应该小心不要让速度入口靠近固体妨碍物,
因为这会导致流动入口驻点
属性具有太高的非一致性。
对于特定的例子,
p>
FLUENT
可能会使用速度入口在流动出口处定义流动速度(在这
种
情况下不使用标量输入)
。在这种情况下,必须保证区域内的
所有流动性。对于流动的概述
请参阅流动入口和出口。
速度入口边界条件的输入
概述
速度入口边界条件需要输入下列信息
?
速度大小与方向或者速度分量。
?
旋转速度(对于具有二维轴对称问
题的涡流)
。
?
温度(用于能量计算)
。
?
Outflow gauge
pressure (for calculations with the coupled
solvers)
?
湍流参数(对于湍流计算)
?
辐射参数
(
对于
P-1
模型、
< br>DTRM
或者
DO
模型的计算<
/p>
)
?
化学组分质量百分数(对于组分计
算)
。
?
混合分数和变化(对于
PDE
燃烧计算
)
。
?
发展变量(对于预混和燃烧计算)
。
?
离散相边界条件(对于离散相计算)
?
二级相的体积分数
(
对于多相流计算
)
p>
上面的所有值都有速度面板输入,它是从边界条件打开的(见设定边界条件一节)
。
Figure
1:
速度入口面板
定义速度
你可以通过定义来确定入口速度。
如果临近速度入口的单元区域是移动的
(也就是说你
使用旋转参考坐标系,多重坐标系或者滑动网格)
,你也可以指定相对速度和绝对速度。对
于
FLUEN
T
中的涡流轴对称问题,你还要指定涡流速度。
定义流入速度的程序如下:
1. <
/p>
选择指定流动方向的方法:
在速度指定方法下拉菜单中选择速度大
小和方向、
速度分量
或者垂直于边界的速度大小。
2.
如果临近速度入口的单元区域是移动的,<
/p>
你可以指定相对或绝对速度。
相对于临近单元
区域或者参考坐标系下拉列表的绝对速度。
如果临近单元区域是固定的,
相对速度和绝对速
度是相等的,这个时候不用察看下拉列表。
3.
如果你想要设定速度的大小和方向或者速度
分量,
而且你的几何图形是三维的,
下一步
你就要选择定义矢量和速度分量的坐标系。坐标系就是前面所述的三种。
4.
设定适当的速度参数,下面将会介绍每一个指定方法。
如果第一步中选择的是速度的大小和方向,
< br>你需要在流入边界条件中输入速度矢量的大
小以及方向。
?
如果是二维非轴对称问题,
或者你在第三步中选择笛卡尔坐标系,
你需要定义流动
X, Y
,
和
(
在三维问题中
)
Z
三个分量的大小。
?
如果是二维轴对称问题
,
,或者第三步中使用柱坐标系,请输入流动方向的径向、轴向
和切向的三个分量值。
?
如果你在第三步中选择当地柱坐标
系,
请输入流动方向的径向、
轴向和切向的三个分量
值。并指定轴向的
X, Y
,
< br>和
Z-
分量以及坐标轴起点的
X
, Y
,
和
Z-
坐标的值。
定义流动方向的
Figure
1
表明这些不同坐标系矢量分量。
如
果你在定义速度的第一步中选择速度大小以及垂直的边界,
你需要在流入边界处输入
p>
速度矢量的大小。
如果你模拟二维轴对称涡流,
你也要输入流向的切向分量。
如果你在定义
速度的第一步中
选择速度分量,你需要在流入边界中输入速度矢量的分量。
?
如果是二维非轴对称问题,
或者你在第三步中选择笛卡尔坐标系,
你需要定义流动
X, Y
,
和
(
在三维问题中
)
Z
三个分量的大小。
?
如果是模拟涡流的二维轴对称问题
,
你需要在速度设定中设定轴向、
径向和旋转速度,
。
?
如果是第三步中使用柱坐标系,
请输入流动方向的径向、
轴向和切向的三个分量值,
以
及(可选)旋转角速度。
?
如果你在第三
步中选择当地柱坐标系,
请输入流动方向的径向、
轴向和切向的
三个分量
值。并指定轴向的
X, Y
,
和
Z-
分量以及坐标轴起点的
X, Y
,
和
Z-<
/p>
坐标的值。
记住速度的正负分量和坐标
方向的正负是相同的。柱坐标系下的速度的正负也是一样。
如
果你在第一步中定义的是速度分量,
并在模拟轴对称涡流,
你可
以指定除了涡流速度之外
的入口涡流角速度
W
< br>。相似地,如果你在第三步中使用柱坐标或者当地柱坐标系,你可以
指定除切向速
度之外的入口角速度
W
。
如果你指定
W,
v_q
作为每个单元的
W r
,
其中
r
从起点到单元的距离。
如果你指定涡流
速度和涡流角速度或者切向速度和角速度,
FLUENT
会将
v_q
和
W
r
加起来获取每个单元
< br>的旋转速度或者切向速度。
定义温度
在解能量方程时,你需要在温度场中的速度入口边界设定流动的静温。
定义流出标准压力
如果你是用一种耦合解算器,
你可以
为速度入口边界指定流出标准压力。
如果在流动要
在任何表面边
界处流出区域,
表面会被处理为压力出口,
该压力出口为流出标
准压力场中规
定的压力。
(
注意:这一
影响和
RAMPANT
中得到的速度远场边界相似。
定义湍流参数
对于湍流计算,
有几种定义湍流参数
的方法。
至于选取哪种方法以及相关的输入值请参阅确
定湍流参
数一节。湍流模型的相关内容请参阅湍流模型一章。
定义辐射参数
如果你打算使用
P-1
辐射模型、
DTRM
或者
DO
模型,
你就需要设定内部发散率以及
(可
选)黑体温度。详
情请参阅设定边界条件一节
(Rosseland
不需要任何边
界条件的输入
)
。
定义组分质量百分比
如果你是用有限速度模型来模拟组分输运,
< br>你就需要设定组分质量百分比。
详情请参阅
组分边界条件
的定义。
定义
PDF/
混合分数参数
如果你用
PDF
模型模拟燃烧,你就需要设定
平均混合分数以及混合分数变化(如果你
是用两个混合分数就还包括二级平均混合分数和
二级混合分数变化)
。具体情况如第三步定
义边界条件所述。<
/p>
定义预混和燃烧边界条件
如果使用与混合燃烧模型,你就需要设定发展变量。请见发展变量的边界条件设定。
定义离散相边界条件
如果你是在模拟粒子的离散相,
你就
可以在速度入口设定粒子轨道详情请参阅离散向模
型的边界设定。
定义多相边界条件
对于多相流如果使用
VOF
,
cavitation
或者代数滑移混合模型
,你就需要指定所有二级
相的体积分数。详情请参阅
VOF
p>
模型、
cavitation
模型或者代数
滑移混合模型的边界设定。
速度入口边界条件的默认设定
速度入
口边界条件的默认设定(国际单位)
:
Temperature
300
Velocity
Magnitude
0
X-Component of Flow Direction
1
Y-Component of
Flow Direction
0
Z-Component of Flow Direction
0
X-Velocity
0
Y-Velocity
0
Z-Velocity
0
Turb. Kinetic
Energy
1
Turb.
Dissipation Rate
1
Outflow Gauge Pressure 0
速度入口边界的计算程序
FLUEN
T
使用速度入口的边界条件输入计算流入流场的质量流以及入口的动量、能量和组
分流量。
本节介绍了通过速度入口边界条件流入流场的算例,
以及通过速度入口边界条件流
出流场的算例。
流动入口的速度入口条件处理
p>
使用速度入口边界条件定义流入物理区域的模型,
FLUENT
p>
既使用速度分量也使用标量。
这些标量定义为边界条件来计算入口质
量流速,动量流量以及能量和化学组分的流量。
邻近速度入口边界流体单元的质量流速由下式计算:
?
?
?
?
v
?
dA
m
注意只有垂直于控制体表面的流动分量才对流入质量流速有贡献。
流动出口的速度入口条件处理
p>
有时速度入口边界条件用于流出物理区域的流动。比如通过某一流域出口的流速已知,
或者被强加在模型上,就需要用这一方法。
注意:这种方法在使用之前必须保证流域内的全部连续性。
<
/p>
在分离解算器中,当流动通过速度入口边界条件流出流场时,
FL
UENT
在边界条件中
使用速度垂直于出口区域的速度分量。<
/p>
它不使用任何你所输入的其它的边界条件。
除了垂直
速度分量之外的所有流动条件,都被假定为逆流的单元。
< br>在耦合解算器中,
如果流动流出边界处的任何表面的区域,
那一表面就会被看成压力出
口,这一压力为
Outflow
Gauge Pressure field
中所规定的压力。
密度计算
入口平面的密度既可以是常数也可以是温度、压力和
/
或组分质量百分数(你在入口条
件中输入的)的函数。
< br>
质量入口边界条件
该边界条件用于规定入口的质量流量。
为了实现规定的质量流量中需要的速度,
就要调
节当地入口总
压。
这和压力入口边界条件是不同的,
在压力入口边界条件中,
规定的是流入
驻点的属性,质量流量的变化依赖于内部解。
p>
当匹配规定的质量和能量流速而不是匹
配流入的总压时,通常就会使用质量入口边界条件。
比如:
一个
小的冷却喷流流入主流场并和主流场混合,
此时,
主流的流速主
要的由
(不同的)
压力入口
/
出口边界条件对控制。
调节入口总压可能会导致节的收敛,所以如果压力入口边界条件和质量入口条件都可以接
受,你应该选择压力入口边界条件。
在不可压流中不必使
用质量入口边界条件,
因为密度是常数,
速度入口边界条件就已
经确定
了质量流。关于流动边界条件的概述请参阅流动入口和出口一节。
质量入口边界条件的输入
概述
质量入口边界条件需要输入:
?
质量流速和质量流量
?
总温(驻点温度)
?
静压
?
流动方向
?
湍流参数(对于湍流计算)
?
辐射参数
(
对于
P-1
模型、
< br>DTRM
或者
DO
模型的计算<
/p>
)
?
化学组分质量百分数(对于组分计
算)
。
?
混合分数和变化(对于
PDE
燃烧计算
)
。
?
发展变量(对于预混和燃烧计算)
。
?
离散相边界条件(对于离散相计算)
上面的所有值都由质量入口面板输入,它是从边界条件打开的(见设定边界条件一节)
。
Figure 1:
质量流动入口面板
定义质量流速度和流量
你可以输入通过质量入口的质量流速,然后
< br>FLUENT
将这个值转换为质量流量,或者直接
指定质
量流量。
如果你设定规定的质量流速,
它将在内部转换为区域上
的规定的统一质量流
量,
这一区域由流速划分。
你也可以使用边界轮廓或者自定义函数来定义质量流量
(不是质
量流速)
。
质量流速或者流量的输入如下:
1.
选择质量流速的方法:质量流速或者质量流量
2.
如果是质量流速(默认)
,在质量流速框中输入规定的质量流速。
注
意:对于轴对称问题,这一质量流速是通过完整区域
(2p-radian)
而不是
1-radian
部分的流
速。
如果选择质量流量。请在
Mass
Flux
框中输入质量流量。
注意:
对于轴对称问题,这一质量流量是通过完整区域
(2p-radian)
而不是
1-radian
部分的流
量。
定义总温
在质量流入口面板中的流入
流体的总温框中输入总温(驻点温度)值。
定义静压
如果入口流动是超声速的,
或者你打算用压力入口边界条件来对解进行初始化,
那么你
p>
必须指定静压
(termed the
Supersonic/Initial Gauge
Pressure)
。
只要流动是压
声速的,
FLUENT
会忽略
Supe
rsonic/Initial Gauge Pressure
,
它是由指定的
驻点值来计算的。
如果你打算使用压力入口边
界条件来初始化解域,
Supersonic/Initial Gauge
Pressure
是与计算初始值的指定驻点压力相联系的,
计算初始值的方法有各向同性关系式
(对
于可压流
)或者贝努力方程(对于不可压流)
。因此,对于压声速入口,它是在关于入口马
赫数(可压流)或者入口速度(不可压流)合理的估计之上设定的。
需要记住的是这个静压和你在操作条件面板中的操作压力是相关的。
请参阅有关于压力
输入和静压头相关输入的解释。
定义流动方向
你可以在压力入口明确的定义流动的方向,
< br>或者定义流动垂直于边界。
对于使用分离解
算器计算移动
区域问题,
流动方向将是绝对速度或者相对于网格相对速度,
这
取决于解算器
面板中的绝对速度公式是否被激活。对于耦合解算器,流动方向通常是绝对
坐标系中的。
定义流动方向的步骤如下,总结请参考概述中的
Figure
1
。
1.
在方向指定下拉菜单中选择指定流动方向的方法,或者是方向矢量或者是垂直于边界。<
/p>
2.
如果你
在第一步中选择垂直于边界,
并且是在模拟轴对称涡流,
请输入
流动适当的切向
速度,
如果你选择垂直于边界并且你的流动是二
维或者三维轴对称涡流,
那就不需要流
动方向上的其它的附加输
入了。
3.
如果第一步中你选择指定方向矢量,
并且你的几何外形是
3<
/p>
维的,
你就需要选择定义矢
量分量的坐标
系统。在坐标系下拉菜单中选择笛卡尔
(X, Y, Z)
坐标
,柱坐标(半径,切
线和轴)
,或者局部柱坐标。
?
如果是二维非轴对称问题或者三维问题,你需要定义流动
X,
Y,
和
(
在三维问题中
) Z
三
个分量的大小。
?
如果是二维轴对称问题
,
,请输入流动方向的径向、轴向和切向的三个分量值。
< br>
定义湍流参数
对于湍流计算,
有几种定义湍流参数
的方法。
至于选取哪种方法以及相关的输入值请参阅确
定湍流参
数一节。湍流模型的相关内容请参阅湍流模型一章。
定义辐射参数
如果你打算使用
P-1
辐射模型、
DTRM
或者
DO
模型,
你就需要设定内部发散率以及
(可
选)黑体温度。详
情请参阅设定边界条件一节
(Rosseland
不需要任何边
界条件的输入
)
。
定义组分质量百分比
如果你是用有限速度模型来模拟组分输运,
< br>你就需要设定组分质量百分比。
详情请参阅
组分边界条件
的定义。
定义
PDF/
混合分数参数
如果你用
PDF
模型模拟燃烧,你就需要设定
平均混合分数以及混合分数变化(如果你
是用两个混合分数就还包括二级平均混合分数和
二级混合分数变化)
。具体情况如第三步定
义边界条件所述。<
/p>
定义预混和燃烧边界条件
如果使用与混合燃烧模型,你就需要设定发展变量。请见发展变量的边界条件设定。
定义离散相边界条件
如果你是在模拟粒子的离散相,
你就
可以在速度入口设定粒子轨道详情请参阅离散向模
型的边界设定。
质量流入口边界的默认设定
质量入口边界条件的默认设定(国际标准单位)为:
Mass Flow-Rate
1
Total Temperature
300
Supersonic/Initial Gauge Pressure
0
X-Component of
Flow Direction
1
Y-Component of Flow Direction
0
Z-Component of
Flow Direction
0
Turb. Kinetic Energy
1
Turb. Dissipation Rate
1
质量流入口边界的计算程序
对入口区域使用质量入口边界条件,
该区域的每一个表面的速度被计算
出来,
并且这一
速度用于计算流入区域的相关解变量的流量。<
/p>
对于每一步迭代,
调节计算速度以便于保证正
确的质量流的数值。
你需要使
用质量流速、流动方向、静压以及总温来计算这个速度。
有两
种指定质量流速的方法。
第一种方法是指定入口的总质量流速
m
(dot)
。
第二种方法是指
定质量流
量
r v (
每个单位面积的质量流速
)
。如果指定总质量流速,
FLUENT
会在内部通过
将总流量除以垂直于流向区域的总入口面积得到统一质量流量:
?
v
?
?
m
A
< br>如果使用直接质量流量指定选项,
可以使用轮廓文件或者自定义函数来指定边界处
的各种质
量流量。
一旦在给定表面的
r v
值确定了,<
/p>
就必须确定表面的密度值
r
,
以找到垂直速度
v
。
密度
获取
的方法依赖于所模拟的是不是理想气体。下面检查了各种情况:
理想气体的质量流边界的流动计算
如果是理想气体,要用下式计算密度:
p
?
?
RT
如果入口是超音速,
所使用的静压是设为边界条件静压值。
p>
如果是亚音速静压是从入口表面
单元内部推导出来的。
入口的静温是从总焓推出的,总焓是从边界条件所设的总温推出的。
入口的密度是从理想气体定律,使用静压和静温推导出来的。
不可压流动的质量流边界的流动计算
如果是模拟非理想气体或者液体,
静温和总温相同。
入口处的密度很容易从温度函数和
(可
选)组分质
量百分比计算出来的。速度用质量入口边界的计算程序中的方程计算出。
质量流边界的流量计算
要计算所有变量在入口处的流量,
流速
v
和方程中变量的入口值一起使用。
例如,
质量流
量
为
r
v
,湍流动能的流量为
r k
v
。这些流量用于边界条件来计算解过程的守恒方程。
进气口边界条件
进气口边界条件用于模拟具有指定损失系数、
流动方向以及环境
(入口)
压力和温度的进气
口。
进气口边界的输入
进气口边界需要输入:
?
总压即驻点压力
?
总温即驻点温度。
?
流动方向
?
静压
?
湍流参数(对于湍流计算)
?
辐射参数
(
对于
P-1
模型、
< br>DTRM
或者
DO
模型的计算<
/p>
)
?
化学组分质量百分数(对于组分计
算)
。
?
混合分数和变化(对于
PDE
燃烧计算
)
。
?
发展变量(对于预混和燃烧计算)
。
?
离散相边界条件(对于离散相计算)
?
二级相的体积分数
(
对于多相流计算
)
?
损失系数
上面的所有值都由进气口面板输入,它是从边界条件打开的(见设定边界条件一节)
p>
。
上面的前十
一项的设定和压力入口边界的设定一样。下面介绍一下损失系数的设定:
Figure 1:
进气口面板
指定损失系数
FLUENT
中的进气口模型,进气口假定为无限薄,通过进气口的压降假定和流体的动压成
< br>比例,
并以经验公式确定你所应用的损失系数。
也就是说
压降
D p
和通过进气口速度的垂直
分
量的关系为:
?
p
< br>?
k
L
1
2
?
v
2
其中
r<
/p>
是流体密度,
k_L
为无量纲的损失系数
。
注意:
D
p
是流向压降,因此即使是在回流中,进气口都会出现阻力。
你可以定义通过进气口的损失系数为常量、
多项式、
分段线性函数或者垂向速度的分段多项
式函数。
定
义这些函数的面板和定义温度相关属性的面板相同,
详情请参阅使用温度相关函
数定义属性一节。
进气扇边界条件
< br>进气扇边界条件用于定义具有特定压力跳跃、
流动方向以及环境
< br>(进气口)
压力和温度
的外部进气扇流动。
进气扇边界的输入
进气扇边界需要输入:
?
总压即驻点压力
?
总温即驻点温度。
?
流动方向
?
静压
?
湍流参数(对于湍流计算)
?
辐射参数
(
对于
P-1
模型、
< br>DTRM
或者
DO
模型的计算<
/p>
)
?
化学组分质量百分数(对于组分计
算)
。
?
混合分数和变化(对于
PDE
燃烧计算
)
。
?
发展变量(对于预混和燃烧计算)
。
?
离散相边界条件(对于离散相计算)
?
二级相的体积分数
(
对于多相流计算
)
?
压力跳跃
上面的所有值都由进气扇面板输入,它是从边界条件打开的(见设定边界条件一节)
p>
。
上面的前十一项的设定和压力入口边界
的设定一样。下面介绍一下压力跳跃的设定:
-
-
-
-
-
-
-
-
-
上一篇:户口本翻译模板及填写注意事项讲解
下一篇:simulink常用模块