fluent边界条件(一) 联系客服

发布时间 : 星期四 文章fluent边界条件(一)更新完毕开始阅读

例如,在雷诺数为50000是湍流强度为4%

湍流尺度l是和携带湍流能量的大涡的尺度有关的物理量。在完全发展的管流中,l被管道的尺寸所限制,因为大涡不能大于管道的尺寸。L和管的物理尺寸之间的计算关系如下: l?0.07L

其中L为管道的相关尺寸。因子0.07是基于完全发展湍流流动混合长度的最大值的,对于非圆形截面的管道,你可以用水力学直径取代L。

如果湍流的产生是由于管道中的障碍物等特征,你最好用该特征长度作为湍流长度L而不是用管道尺寸。

注意:公式l?0.07L并不是适用于所有的情况。它只是在大多数情况下得很好的近似。对于特定流动,选择L和l的原则如下: ? 对于完全发展的内部流动,选择强度和水力学直径指定方法,并在水力学直径流场中指

定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-Allmaras模型)或者强度和粘性比(对于k-e模型或者RSM)。

推导湍流量的关系式

要获得更方便的湍流量的输运值,如:I, L,或者m_t/m,你必须求助于经验公式,下面是FLUENT中常用的几个有用的关系式。要获得修改的湍流粘性,它和湍流强度I长度尺度l有如下关系:

~?3uIl vavg2在Spalart-Allmaras模型中,如果你要选择湍流强度和水力学直径来计算l可以从前面的公式中获得。

湍动能k和湍流强度I之间的关系为: k?3?uavgI?2 2其中u_avg为平均流动速度

除了为k和e指定具体的值之外,无论你是使用湍流强度和水力学直径,强度和长度尺度或者强度粘性比方法,你都要使用上述公式。

如果你知道湍流长度尺度l你可以使用下面的关系式:

??C?34k l32其中C?是湍流模型中指定的经验常数(近似为0.09),l的公式在前面已经讨论了。

除了为k和e制定具体的值之外,无论你是使用湍流强度和水力学直径还是强度和长度尺度,你都要使用上述公式。

E的值也可以用下式计算,它与湍流粘性比m_t/m以及k有关:

k2??t??????C? ??????其中C?是湍流模型中指定的经验常数(近似为0.09)。

除了为k和e制定具体的值之外,无论你是使用湍流强度和水力学直径还是强度和长度

尺度,你都要使用上述公式。

如果你是在模拟风洞条件,在风洞中模型被安装在网格和/或金属网格屏下游的测试段,你可以用下面的公式:

?1???kU? L?其中,?k是你希望的在穿过流场之后k的衰减(比方说k入口值的10%), U?自由流的速度

L?是流域内自由流的流向长度Equation 9是在高雷诺数各向同性湍流中观察到的幂率衰减

的线性近似。它是基于衰减湍流中k的精确方程U ?k/?x = - e.

如果你用这种方法估计e,你也要用方程7检查结果的湍流粘性比m_t/m,以保证它不是太大。

虽然这不是FLUENT内部使用的方法,但是你可以用它来推导e的常数自由流值,然后你可以用湍流指定方法下拉菜单中选择K和Epsilon直接指定。在这种情况下,你需要使用方程3从I来计算k。

当使用RSM时,如果你不在雷诺应力指定方法的下拉列表中使用雷诺应力选项,明显的制定入口处的雷诺应力值,它们就会近似的由k的指定值来决定。湍流假定为各向同性,保证

uiuj?0

以及

u?u??2k 3(下标a不求和).

如果你在雷诺应力指定方法下拉列表中选择K或者湍流强度,FLUENT就会使用这种方法。

对大涡模拟(LES)指定入口湍流

大涡模拟模型一节中所描述的LES速度入口中指定的的湍流强度值,被用于随机扰动入口处速度场的瞬时速度。它并不指定被模拟的湍流量。正如大涡模拟模型中介绍的边界条件中所描述的,通过叠加每个速度分量的随机扰动来计算流动入口边界处的随机成分.

压力入口边界条件

压力入口边界条件用于定义流动入口的压力以及其它标量属性。它即可以适用于可压流,也可以用于不可压流。压力入口边界条件可用于压力已知但是流动速度和/或速率未知的情况。这一情况可用于很多实际问题,比如浮力驱动的流动。压力入口边界条件也可用来定义外部或无约束流的自由边界。对于流动边界条件的概述,请参阅流动入口和出口一节。

压力入口边界条件的输入 综述

对于压力入口边界条件你需要输入如下信息 ? 驻点总压 ? 驻点总温 ? 流动方向 ? 静压

? 湍流参数(对于湍流计算)

? 辐射参数(对于使用P-1模型、DTRM模型或者DO模型的计算) ? 化学组分质量百分比(对于组分计算) ? 混合分数和变化(对于PDF燃烧计算) ? 程序变量(对于预混和燃烧计算) ? 离散相边界条件(对于离散相的计算) ? 次要相的体积分数(对于多相计算)

所有的值都在压力入口面板中输入(Figure 1),该面板是从边界条件打开的。

Figure 1: 压力入口面板

压力输入和静压头

压力场(p_s^')和压力输入(p_s^' or p_0^')包括静压头r_0 g x。也就是FLUENT 以下式定义的压力:

p?s??0gx?ps

或者

?p??ps??0g?s ?xx这一定义允许静压头放进体积力项(r - r_0)g中考虑,而且当密度一致时,从压力计算中排除了。因此你的压力输入不因该考虑静压的微分,压力(p^'_s)的报告也不会显示静压的任何影响。有关浮力驱动流动的内容请参阅浮力驱动流动和自然对流的信息

定义总压和总温

在压力入口面板中的Gauge Total Pressure field输入总压值。总温会在Total Temperature field中设定。记住,总压值是在操作条件面板中定义的与操作压力有关的的总压值。不可压流体的总压定义为:

p0?ps??v

2