边界条件定义边界条件概述边界条件包括流动变量和热变量在边界处的值。它是FLUENT分析得很关键的一部分,设定边界条件必须小心谨慎。边界条件的分类:进出口边界条件:压力、速度、质量进口、进风口、进气扇、压力出口、压力远场边界条件、质量出口、通风口、排气扇;壁面、repeating,andpoleboundaries:壁面,对称,周期,轴;内部单元区域:流体、固体(多孔是一种流动区域类型);内部表面边界:风扇、散热器、多孔跳跃、壁面、内部。(内部表面边界条件定义在单元表面,这意味着它们没有有限厚度,并提供了流场性质的每一步的变化。这些边界条件用来补充描述排气扇、细孔薄膜以及散热器的物理模型。内部表面区域的内部类型不需要你输入任何东西。)下面一节将详细介绍上面所叙述边界条件,并详细介绍了它们的设定方法以及设定的具体合适条件。周期性边界条件在本章中介绍,模拟完全发展的周期性流动将在周期性流动和热传导一章中介绍。使用边界条件面板边界条件(Figure1)对于特定边界允许你改变边界条件区域类型,并且打开其他的面板以设定每一区域的边界条件参数菜单:Define/BoundaryConditions...Figure1:边界条件面板改变边界区域类型设定任何边界条件之前,必须检查所有边界区域的区域类型,如有必要就作适当的修改。比方说:如果你的网格是压力入口,但是你想要使用速度入口,你就要把压力入口改为速度入口之后再设定。改变类型的步骤如下::1.在区域下拉列表中选定所要修改的区域2.在类型列表中选择正确的区域类型3.当问题提示菜单出现时,点击确认确认改变之后,区域类型将会改变,名字也将自动改变(如果初始名字时缺省的请参阅边界条件区域名字一节),设定区域边界条件的面板也将自动打开。!注意:这个方法不能用于改变周期性类型,因为该边界类型已经存在了附加限制。创建边界条件一节解释了如何创建和分开周期性区域。需要注意的是,只能在图一中每一个类别中改变边界类型(注意:双边区域表面是分离的不同单元区域.)Figure1:区域类型的分类列表设定边界条件在FLUENT中,边界条件和区域有关而与个别表面或者单元无关。如果要结合具有相同边界条件的两个或更多区域请参阅合并区域一节。设定每一特定区域的边界条件,请遵循下面的步骤:1.在边界条件区域的下拉列表中选择区域。2.点击Set...按钮。或者,1.在区域下拉列表中选择区域。2.在类型列表中点击所要选择的类型。或者在区域列表中双击所需区域.,选择边界条件区域将会打开,并且你可以指定适当的边界条件在图像显示方面选择边界区域在边界条件中不论你合适需要选择区域,你都能用鼠标在图形窗口选择适当的区域。如果你是第一次设定问题这一功能尤其有用,如果你有两个或者更多的具有相同类型的区域而且你想要确定区域的标号(也就是画出哪一区域是哪个)这一功能也很有用。要使用该功能请按下述步骤做:1.用网格显示面板显示网格。2.用鼠标指针(默认是鼠标右键——参阅控制鼠标键函数以改变鼠标键的功能)在图形窗口中点击边界区域。在图形显示中选择的区域将会自动被选入在边界条件面板中的区域列表中,它的名字和编号也会自动在控制窗口中显示改变边界条件名字每一边界的名字是它的类型加标号数(比如pressure-inlet-7)。在某些情况下你可能想要对边界区域分配更多的描述名。如果你有两个压力入口区域,比方说,你可能想重名名它们为small-inlet和large-inlet。(改变边界的名字不会改变相应的类型)重名名区域,遵循如下步骤:1.在边界条件的区域下拉列表选择所要重名名的区域。2.点击Set...打开所选区域的面板。3.在区域名字中输入新的名字4.点击OK按钮。注意:如果你指定区域的新名字然后改变它的类型,你所改的名字将会被保留,如果区域名字是类型加标号,名字将会自动改变。边界条件的非一致输入每一类型的边界区域的大多数条件定义为轮廓函数而不是常值。你可以使用外部产生的边界轮廓文件的轮廓,或者用自定义函数(UDF)来创建。具体情况清参阅相关内容流动入口和出口FLUENT有很多的边界条件允许流动进入或者流出解域。下面一节描述了每一种边界条件的类型的使用以及所需要的信息,这样就帮助你适当的选择边界条件。下面还提供了湍流参数的入口值的确定方法。使用流动边界条件下面对流动边界条件的使用作一概述对于流动的出入口,FLUENT提供了十种边界单元类型:速度入口、压力入口、质量入口、压力出口、压力远场、质量出口,进风口,进气扇,出风口以及排气扇。下面是FLUENT中的进出口边界条件选项:速度入口边界条件用于定义流动入口边界的速度和标量压力入口边界条件用来定义流动入口边界的总压和其它标量。质量流动入口边界条件用于可压流规定入口的质量流速。在不可压流中不必指定入口的质量流,因为当密度是常数时,速度入口边界条件就确定了质量流条件。压力出口边界条件用于定义流动出口的静压(在回流中还包括其它的标量)。当出现回流时,使用压力出口边界条件来代替质量出口条件常常有更好的收敛速度。压力远场条件用于模拟无穷远处的自由可压流动,该流动的自由流马赫数以及静态条件已经指定了。这一边界类型只用于可压流。质量出口边界条件用于在解决流动问题之前,所模拟的流动出口的流速和压力的详细情况还未知的情况。在流动出口是完全发展的时候这一条件是适合的,这是因为质量出口边界条件假定出了压力之外的所有流动变量正法向梯度为零。对于可压流计算,这一条件是不适合的。进风口边界条件用于模拟具有指定的损失系数,流动方向以及周围(入口)环境总压和总温的进风口。进气扇边界条件用于模拟外部进气扇,它具有指定的压力跳跃,流动方向以及周围(进口)总压和总温。通风口边界条件用于模拟通风口,它具有指定的损失系数以及周围环境(排放处)的静压和静温。排气扇边界条件用于模拟外部排气扇,它具有指定的压力跳跃以及周围环境(排放处)的静压。决定湍流参数在入口、出口或远场边界流入流域的流动,FLUENT需要指定输运标量的值。本节描述了对于特定模型需要哪些量,并且该如何指定它们。也为确定流入边界值昀为合适的方法提供了指导方针。使用轮廓指定湍流参量在入口处要准确的描述边界层和完全发展的湍流流动,你应该通过实验数据和经验公式创建边界轮廓文件来完美的设定湍流量。如果你有轮廓的分析描述而不是数据点,你也可以用这个分析描述来创建边界轮廓文件,或者创建用户自定义函数来提供入口边界的信息。一旦你创建了轮廓函数,你就可以使用如下的方法:Spalart-Allmaras模型:在湍流指定方法下拉菜单中指定湍流粘性比,并在在湍流粘性比之后的下拉菜单中选择适当的轮廓名。通过将m_t/m和密度与分子粘性的适当结合,FLUENT为修改后的湍流粘性计算边界值。k-e模型:在湍流指定方法下拉菜单中选择K和Epsilon并在湍动能(Turb.KineticEnergy)和湍流扩散速度(Turb.DissipationRate)之后的下拉菜单中选择适当的轮廓名。雷诺应力模型:在湍流指定方法下拉菜单中选择K和Epsilon并在湍动能(Turb.KineticEnergy)和湍流扩散速度(Turb.DissipationRate)之后的下拉菜单中选择适当的轮廓名。在湍流指定方法下拉菜单中选择雷诺应力部分,并在每一个单独的雷诺应力部分之后的下拉菜单中选择适当的轮廓名。湍流量的统一说明在某些情况下流动流入开始时,将边界处的所有湍流量指定为统一值是适当的。比如说,在进入管道的流体,远场边界,甚至完全发展的管流中,湍流量的精确轮廓是未知的。在大多数湍流流动中,湍流的更高层次产生于边界层而不是流动边界进入流域的地方,因此这就导致了计算结果对流入边界值相对来说不敏感。然而必须注意的是要保证边界值不是非物理边界。非物理边界会导致你的解不准确或者不收敛。对于外部流来说这一特点尤其突出,如果自由流的有效粘性系数具有非物理性的大值,边界层就会找不到了。你可以在使用轮廓指定湍流量一节中描述的湍流指定方法,来输入同一数值取代轮廓。你也可以选择用更为方便的量来指定湍流量,如湍流强度,湍流粘性比,水力直径以及湍流特征尺度,下面将会对这些内容作一详细叙述。湍流强度I定义为相对于平均速度u_avg的脉动速度u^'的均方根。小于或等于1%的湍流强度通常被认为低强度湍流,大于10%被认为是高强度湍流。从外界,测量数据的入口边界,你可以很好的估计湍流强度。例如:如果你模拟风洞试验,自由流的湍流强度通常可以从风洞指标中得到。在现代低湍流风洞中自由流湍流强度通常低到0.05%。.对于内部流动,入口的湍流强度完全依赖于上游流动的历史,如果上游流动没有完全发展或者没有被扰动,你就可以使用低湍流强度。如果流动完全发展,湍流强度可能就达到了百分之几。完全发展的管流的核心的湍流强度可以用下面的经验公式计算:81Re16.0HDavguuI例如,在雷诺数为50000是湍流强度为4%湍流尺度l是和携带湍流能量的大涡的尺度有关的物理量。在完全发展的管流中,l被管道的尺寸所限制,因为大涡不能大于管道的尺寸。L和管的物理尺寸之间的计算关系如下:Ll07.0其中L为管道的相关尺寸。因子0.07是基于完全发展湍流流动混合长度的昀大值的,对于非圆形截面的管道,你可以用水力学直径取代L。如果湍流的产生是由于管道中的障碍物等特征,你昀好用该特征长度作为湍流长度L而不是用管道尺寸。注意:公式Ll07.0并不是适用于所有的情况。它只是在大多数情况下得很好的近似。对于特定流动,选择L和l的原则如下:对于完全发展的内部流动,选择强度和水力学直径指定方法,并在水力学直径流场中指定L=D_H。对于旋转叶片的下游流动,穿孔圆盘等,选择强度和水力学直径指定方法,并在水力学直径流场中指定流动的特征长度为L对于壁面限制的流动,入口流动包含了湍流边界层。选择湍流强度和长度尺度方法并使用边界层厚度d_99来计算湍流长度尺度l,在湍流长度尺度流场中输入l=0.4d_99这个值湍流粘性比m_t/m直接与湍流雷诺数成比例(Re_t?k^2/(en))。Re_t在高湍流数的边界层,剪切层和完全发展的管流中是较大的(100到1000)。然而,在大多数外流的自由流边界层中m_t/m相当的小。湍流参数的典型设定为1m_t/m10。要根据湍流粘性比来指定量,你可以选择湍流粘性比(对于Spalart-Allmaras模型)或者强度和粘性比(对于k-e模型或者RSM)。推导湍流量的关系式要获得更方便的湍流量的输运值,如:I,L,或者m_t/m,你必须求助于经验公式,下面是FLUENT中常用的几个有用的关系式。要获得修改的湍流粘性,它和湍流强度I长度尺度l有如下关系:Iluvavg23~在Spalart-Allmaras模型中,如果你要选择湍流强度和水力学直径来计算l可以从前面的公式中获得。湍动能k和湍流强度I之间的关系为:223Iukavg其中u_avg为平均流动速度除了为k和e指定具体的值之外,无论你是使用湍流强度和水力学直径,强度和长度尺度或者强度粘性比方法,你都要使用上述公式。如果你知道湍流长度尺度l你可以使用下面的关系式:lkC2343其中C是湍流模型中指定的经验常数(近似为0.09),l的公式在前面已经讨论了。除了为k和e制定具体的值之外,无论你是使用湍流强度和水力学直径还是强度和长度尺度,你都要使用上述公式。E的值也可以用下式计算,它与湍流粘性比m_t/m以及k有关:12tkC其中C是湍流模型中指定的经验常数(近似为0.09)。除了为k和e制定具体的值之外,无论你是使用湍流强度和水力学直径还是强度和长度尺度,你都要使用上述公式。如果你是在模拟风洞条件,在风洞中模型被安装在网格和/或金属网格屏下游的测试段,你可以用下面的公式:LkU其中,k是你希望的在穿过流场之后k的衰减(比方说k入口值的10%),U自由流的速度L是流域内自由流的流向长度Equation9是在高雷诺数