管内湍流的数值模拟

1、管内湍流的数值模拟摘 要：当Rey nolds数大于临界值时，平滑流动会出现一系列复杂的变化，最终 会导致流动特征的本质变化， 流动呈无序的混乱状态， 这种状态称为湍流。 计算 流体力学是通过计算机数值计算和图像显示， 对包含有流体流动和热传导等相关 物理现象的系统所做的分析。本文以湍流管流模型为例，借助 Flue nt软件进行空 气动力学分析，对该管内湍流流动进行模拟。关键词：计算流体力学；Flue nt;管内湍流；数值模拟1 引言流体试验表明，当Rey no Ids数大于临界值时，平滑流动会出现一系列复杂的 变化，最终会导致流动特征的本质变化，流动呈无序的混乱状态。这时，即使是 边界条件保

2、持不变， 流动也是不稳定的， 速度等流动特性都随机变化， 这种状态 称为湍流。随着高速电子计算机的出现，数值模拟越来越多地应用于流场的模拟。计算 流体力学（Computational Fluid Dynamics ,简称为CFD）就是其中一种有效的研究 流体动力学的数值模拟方法， 它是通过计算机数值计算和图像显示， 对包含有流 体流动和热传导等相关物理现象的系统所做的分析； 是基于计算机技术的一种数 值计算工具 , 用于求解流体的流动和传热问题。 它能够描述几何体边界的复杂的 流动现象， 能够在设计的初期快速地评价设计并做出修改； 在设计的中期， 用来 研究设计变化对流动的影响，减少未预料到的

3、负面影响；设计完成后，CFD提供各种数据和图像，证实设计目的。CFD大大减少了费用、时间以及新设计带来的 风险。近年来，CFD越来越多地应用于翼型设计和流场的分析中，成为一种重要 的设计和计算方法。Flue nt软件是用于模拟和分析在复杂几何区域内的流体流动与热交换问题的 专用CFD软件。它用于计算计算流体流动和传热问题的软件，其应用的范围有一 般流体的流场、自由表面的问题、紊流、非牛顿流流场、化学反应等。Flue nt提供了灵活的网格特性， 用户可以方便的使用结构网格和非结构网格对各种复杂区 域进行网格划分。本文以湍流管流模型为例，借助 Flue nt软件进行空气动力学分 析，对该管内湍流流

4、动进行模拟， 并分析了模型内的中心速度分布、 表面摩擦系 数和流速剖面。数学及物理模型的建立2.1 数学模型Flue nt软件的理论基础是计算流体力学，数值计算的控制方程是二维连续性 方程、二维N-S方程和二维能量方程。连续性方程：cP + 必 PUx ） + & Puy）;:tN-S方程:(PUx ) . u;:t.(-Ux)y；xI uxux2矽丿_ x2:x.(My )Uy一(My ).7x能量方程：匸(PE )+ 匕盯一无 hjjj +(“ff v) +ShFLUENT的能量方程式如下：jt这里keff是有效的传导率，Jj为组分j的扩散流量，等式右边分别表示了由 于传导，组分扩散和粘性

5、扩散导致的能量转移。 Sh包含了化学反应热和以定义 的其他体积热。2式中E=h-P+，其中理想气体的焓h如下定义：h=v Yjhj，如为不可压 p 2 j缩流体则定义为h二Yjhj+P ,以上两式中Yj为组分j的质量分数，且 jPThj 二Jcp,jdT Tref = 298.15。Tef2.2物理模型的建立及网格划分该湍流模型如图1所示，根据管道直径和平均流速可以计算得出:PVDRed10000d .1因此，流动属于完全湍流。湍流显示的是瞬时的小规模波动，这在CFD中常常不易解决计算，因此，本模型采用速度、压力的时均值计算。我们将现有的模 型微调一下就可以建立本次模型，这次我们将利用原来的k

6、- 模型。通常，标准的k- 莫型是湍流应用的首选模型，它是在关于湍动能k的方程的基础上，又引入 了一个关于湍动耗散率&的方程形成的。它通常能够提供流动的真实情况，尤其 是适合于计算管道和通道中的湍流流动。由于本文所研究的空间是管内，所以采用此模型来进行数学模拟。通过Flue nt的前置处理软件Gambit进行模型的创建（即划分网格），采取由 低阶元素到高阶元素（点一线一面）的方式。该模型的网格划分见图2。图1湍流几何模型图2模型网格划分2.3边界条件及初始条件边界条件如表1所示，左侧为流速进口，进口速度为1 m/s,右侧设置为压力 出口，出口的压力值设置为1 atm.。上方为壁面，下部为轴对称

7、的对称中心线。表1定义的边界条件EdgePositi onNameTypeLeftin letVELOCITY INLETRightoutletPRESSURE OUTLETTopwallWALLBottomcen terl ineAXIS初始的运行条件设置为抽对称空间的粘性流体流动， 工作压力采用系统的默 认值。开始运行时将轴向速度初始值设置为1 m/s，径向速度设置为0 m/s，然后 就可以开始计算了。如图3所示，计算至第229次时结果便收敛。(J FLUEBT 0 FliieM IncJ回区ReSi duals： cointlfiuliyvt：l 白皿丫 k usil onle-MX -

8、5 le-HB -=IterationsJibDSIDFllfENTGlfnl.i.p&ir.rtiJ图3计算结果数值模拟结果3.1中心速度分析沿中心对称轴的速度分布如图4所示:图4中心线处的速度分布我们可以看出，流速充分发展地区是从 x=5m处开始的，之后的中心处的流 速约是个常数值，大小为1.195 m/so这显然要低于层流状态下的2m/s的流速。图6表面摩擦系数分布In lx1Zk*CD1J06+ED1.*0)1 I444DPosition (iti)FLUEIT 0 Fluent IncJnQ.aiDT比ENTiiJ輸中,海顷rh)图5 在100X60的网格下的中心线速度分布此图是在划

9、分网格为100X60的状态下的中心流速，由图可以看出，二者的 区别不大。3.2表面摩擦系数分析表面摩擦系数的计算公式如下:TCfw一仁2ref Vref2其中TW是壁面切应力，变化如图6所示：Pef为参考密度，Vyf为参考速度。表面摩擦系数的FLUE1T E01 Fluent Inc2Oe-O2 -iSkinFriction应-2 me-02 -1fDe-D2 -畑血-1Zte-CQ -D11I123i5PostHbn (m)STSFLWENT6J l.pRpbhfefl)Co evident可以看出，在充分发展区的表面摩擦系数为0.0085,并且保持稳定3.3流速剖面 此处描述的是轴向流速在

10、剖面的分布情况,-:JKi*-02 -Posmciin6皿在-$ -IJW-0C =4剛憾畑鼻F LMBIMT 6J 削血pii Hi)2J&-CC lIEt心-0曲如!110业二0.G吋OBAxial Velocity图7轴向流速剖面图由图7可知，轴向速度在中心线处为最大值， 在近壁面处为最小值0。此速度 分布规律满足无滑边界条件的粘性流状况。 近壁面处的速度梯度要比层流状态的 更大些。同样，在100X60的网格中得出的轴向速度分布如图8所示，可以看出二 者区别不大。PasHion*曲曲PositionHiGIFWEhTTM 阿華 ptaH.睛，图8在100X60的网格中的轴向流速剖面图湍流

11、会被壁面状态显著影响，k- 模型要求在近壁面为其主要有效距离。近 壁面模型对网格划分特别敏感，图9既是表现y+的近壁面处对网格划分的敏感程 度。Atll YfMEJl hEKl.SiaFLUENTS m 中.|耕.图9网格划分对y+处的影响从图9中可以看出，当采用100X60的网格时，计算的精确性显著增加。结论本文采用具有模拟复杂外形的流体流动及热传导的CFD软件FLUENT，建立了一个管内湍流模型；选用k- &湍流模型对二维空间进行数值模拟。通过计算 结果分析得出管内湍流模型流场的速度分布和压力分布，管内中心线处的速度在进入充分发展区之后便是以稳定值，表面摩擦系数在进入充分发展区域也是一稳 定值，轴向速度在中心线处为最大值，在近壁面处为最小值0。采用