logo资料库

CFD各种压力总结.docx

第1页 / 共47页
第2页 / 共47页
第3页 / 共47页
第4页 / 共47页
第5页 / 共47页
第6页 / 共47页
第7页 / 共47页
第8页 / 共47页
资料共47页,剩余部分请下载后查看
可视化 k
表压指的是压力表测压值。根据目前压力表的工作原理很容易知 道表压是一种相对压力,为真实压力与大气压间的差值。大气压 101325Pa,若表压值为零,则意味着此处真实压力为101325Pa。静 压就是我们寻常意义上的压力,可以用压力表测量获得。其值是一个 绝对压力值,但是在软件中表现为一个相对于操作压力的相对值。在 fluent 中静压的英文名称为 static pressure。 总压是静压与动压的和。在 FLUENT 的压力入口中常要求用户 输入总压值,其实这里是同时考虑了压力与速度的。因此在压力入口 位置仅仅只是输入压力表读数是不对的,因当考虑速度的。当然若此 处速度为 0 的话,总压值与静压值相等。 CFD 软件计算的压力值都是相对值。若想得到绝对压力值,可 设置参考压力值为 0。FLUENT 默认参考压力值为一个大气压 101325Pa。 至于驻点压力,指的是流体从无穷远处等熵滞止于相对于流场运 动的钝体前驻点处时的压力,也就是(相对)总压。 驻点压力是指总压,总压等于动压加上静压。即里面的 total pressure。这里的gauge pressure 是指静压.而且total pressure 和 gauge pressure 都是是相对于“operating condition”栏里设置的压强 的差值。 可压缩流动中,入口可以设置压力也可以设置流量。若为压力入 口,则需要设置静压值,此时若出口为静压出口时,则压力出口失效, 出口的压力是通过内部迭代计算得到。可压缩流动中,若入口为流量
边界,则计算域内总压不守恒,流量守恒。求解器通过调整总压值以 满足流量要求。此时收敛会很难。因此,若流量与压力均已知的情况 下,优先使用压力入口。 可压缩流动: (1)总压和总温是给定的; (2)如果入口为亚音速流动,静压由计算域中紧邻边界的内部 流体单元值外插得到;如果入口处的流动为超音速的,静压必须是给 定的。 在 Turbulence Specification Method (湍流定义方法)下拉列表 中,可以简单地用一个常数来定义湍流参数,即通过给定湍流强度、 湍流粘度比、水力直径或湍流特征长在边界上的值来定义流场边界上 的湍流。 (1)湍流强度(Turbulence Intensity) 湍流强度 I 的定义如下: 16.0I 1- 8 (Re) 3 2 ( u 湍动能 k  I 2) , avgu 为平均流速。 avg 湍流的长度尺度: 007 L , L 取管道的水力直径。 .0 l kC 43  23 l 湍流耗散率: 值约等于 0.09。 , C 为湍流模型中的一个经验常数,其 在使用压强入口边界条件时需要输入下列参数: (1)总压。 (2)总温。
(3)流动方向。 (4)静压。 (5)用于湍流计算的湍流参数。 上述变量均在 Pressure Inlet(压强入口)面板中输入。在 Boundary Conditions(边界条件)面板中选择压强入口边界,然后点击 set(设 置)按钮就可以进入压强入口条件的设置面板。 1. 定义总压 在 Pressure Inlet(压强入口)面板中的 Gauge Total Pressure(表 总压)栏中输入总压的值,在 Total Temperature(总温)栏中输入总 温。这里总压输入仍然采用表压形式,即输入的总压值实际上是绝对 总压与操作压强之差。 不可压流中总压用伯努力方程定义,即: P 0  sP  2 1   2 (1) 在可压流中用满足等熵关系的伯努力方程定义,即: P 0  P s 1     1   2 2 M  1      (2) 在公示(1)和(2)中, 0P 为实际总压, sP 为实际静压, M 为马赫 数,为比热比。 2.定义流动方向的步骤如下: (1)选择流动方向的定义方法,即在 Direction Specification Method(方向定义方法)下拉列表中选择 Direction Vector(方向矢量) 或 Normal to Boundary(垂直于边界)。
(2)如果选择了 Normal to Boundary(垂直于边界),并且准备 计算的是轴对称旋转流问题,则要在 Tangential-Component of Flow Direction(流动方向的分量)中输入适当的数值。如果选择的是 Normal to Boundary(垂直于边界),同时采用了直角坐标,并且不是轴对称 旋转流动问题,则无需输入任何参数。 (3)如果在前面步骤中选择的是 Direction Vector(方向矢量), 并且所计算的问题是三维问题,随后就需要选择矢量分量的坐标系形 式,即在 Coordinate System(坐标系)下拉列表中选择 Cartesian(X,Y,Z) (直角坐标)、Cylindrical(Radial,Tangential,Axial)(柱坐标)或 Local Cylindrical(Radial,Tangential,Axial)(局部柱坐标)。 这里选用的直角坐标系就是网格划分时使用的直角坐标系。 3.定义静压 静压在 FLUENT 中被称为Supersonic/Initial Gauge Pressure (超音速/初始表压),如果入口流动是超音速的,或者用户准备用 压强入口边界条件进行计算的初始化工作,则必须定义静压。
在流场为亚音速流场时,FLUENT 将忽略 Supersonic/Initial Gauge Pressure(超音速/初始表压)的输入数据,而用驻点参数求出 静压。然而如果使用压强入口条件进行流场的初始化,则需要同时输 入静压和驻点参数,以便可以用(2)式或者(1)式计算出初始流场。 因为上述两个公式中使用了马赫数和速度,所以在定义边界条件时需 要给出合理的入口马赫数(可压流)或入口速度(不可压流)。 设想的气流速度为 0 时候的截面,此时的气流的参数称驻点参 数,是人为的规定一种状态,借用这种状态的气流参数来表征实际流 体的运动状态。 压强入口边界的计算流程 FLUENT 处理压强入口边界条件的方式是在假定没有流动损失 的前提下,由驻点条件 计算出入口条件。对于不可压流动,是从(1)式推导出入口边 界条件。对于可压流,则用(2)式推导出入口边界条件。 可压流压强入口边界的计算 在可压流计算中,压强入口边界的总压、静压和速度可以用理想 气体的等熵关系联系起来。用户在入口边界上输入的总压 0 'P 和相邻 单元上的静压 SP ' 用下列等熵关系联系在一起: ' ' P P 0 s   P op P op     1  1   2 2 M  1      其中马赫数定义为: M   c   sRT  (3) (4)
c 为音速, opP 为操作压强,因为输入的边界条件是以操作压强为 参照的相对压强,所以实际压强的值应该是输入压强与操作压强之 和。 压强出口边界条件 压强出口边界条件在流场出口边界上定义静压,而静压的值仅在 流场为亚音速时使用。如果在出口边界上流场达到超音速,则边界上 的压强将从流场内部通过插值得到。其他流场变量均从流场内部插值 获得。 在压强出口边界上还需要定义“回流(backflow)”条件。回流条 件是在压强出口边界上出现回流时使用的边界条件。推荐使用真实流 场中的数据做回流条件,这样计算将更容易收敛。 压强出口边界的输入参数如下: (1)静压。 (2)回流条件,其中包括: 能量计算中的总温。 回流方向定义方法。 湍流计算中的湍流参数 操作压强的设定 1.操作压强的意义 操作压强对于不可压理想气体流动和低马赫数可压流动来说是 十分重要的,因为不可压理想气体的密度是用操作压强通过状态方程
直接计算出来的(参见公式  P op R M w T ),而在低马赫数可压流动中 操作压强则起到了避免截断误差的负面影响的重要作用。 对于高马赫数可压缩流动,操作压强的意义不大。在这种计算中 使用绝对压力会更方便。因为 FLUENT 总是使用表压进行计算,所 以需要在这类问题的计算中将操作压强设置为零,而使表压和绝对压 力相等。 压缩空气的压力入口条件的设置(喷气织机为例) 图 1.主喷嘴 图 2.主喷嘴三维模型 图 3.主喷嘴流场区域名称图
分享到:
收藏