FLUENT模拟时将圆柱个性字体转换器可粘贴为二维的，进口流量要怎么重新设置
点击联系发帖人
时间：2018-04-10 08:25
您还可以使用以下方式登录
当前位置：&>&&>& > 中科大Fluent讲稿(完整版,完美整理).
中科大Fluent讲稿(完整版,完美整理).
导读：就爱阅读网友为您分享以下“中科大Fluent讲稿(完整版,完美整理).”的资讯，希望对您有所帮助，感谢您对92to.com的支持!2.1.2 动量守恒方程 惯性坐标系下，i方向的动量守恒方程为：
??t(?ui)???xj(?uiuj)???p?xi???ij?cj??gi?Fi
???u?uji???????x?xi?j????2?ul?????ij ，?gi，Fi?3?xl???式中，p是静压；?ij是应力张量，定义为：?ij是重力体积力和其它体积力（如源于两相之间的作用），Fi还可以包括其它模型源项或者用户自定义源项。 对于二维轴对称几何条件，轴向和轴向的动量守恒方程分别为： ??t(?u)?1?r?x(r?uu)?1?r?r(r?vu)???p?x ????1????u2r?2?(??v)?? ??r?x???x3??1????u?v??r??r?r???r?x?????Fx
2－4 ???和 ??t(?v)?1?r?x(r?uv)?1?r?r(r?vv)???p?r ?1????v?u????? ?r??r?x???x?r?????1????v2r?2?(??v)?? ??r?r???x3????2?vr2?w?(??v)???Fr
2－5 3rr2?2w是旋流速度。
2.1.3 能量方程
FLUENT可以计算流体和（或者）固体区域之间的传热问题。如果是周期性换热流动，则流动边界要给定周期边界条件。如果计算计算模型包括两个流动区域，中间被固体或者墙壁隔开的换热问题，则要特别注意：1，两个流体都不能用流出边界条件（outflow）；2，两个区域的流动介质可以不同，但要分别定义流体性质（如果计算组分，只能给一个混合组分）。
ix 流体1 流体2
FLUENT求解的能量方程形式如下： ??t(?E)???xi(ui(?E?p))???xi(keff?T?xi??hj?Jj??uj(?ij)eff?Sh
2－6 j?式中，keff?kt?k，为有效导热系数(湍流导热系数根据湍流模型来定义)。Jj?是组分j?的扩散通量。方程右边前三项分别为导热项，组分扩散项和粘性耗散项。Sh是包括化学反应热和其它体积热源的源项。其中， 2pui
2－7 E?h???2对于理想气体，焓定义为：h??mj?hj?；对于不可压缩气体，焓定义为：j?h??mj?hj??j?pT?。mj?是组分j?的质量分数，组分j?的焓定义为：hj???cp,j?dT，其中TrefTref?298.15K。
2.1.4 PDF模型的能量方程
如果在非绝热PDF燃烧模型模式下，FLUENT求解的总焓方程为： ??t(?H)???xi(?uiH)??(kt?H)??ik?ui?xk?Sh
2－8 ?xicp?xi假定刘易斯数为1，方程右边第一项为组分扩散和导热项的合并项；第二项为粘性耗散，为非守恒形式。总焓H定义为: H??mj?Hj?j? 组分j?的总焓定义为： THj???cp,j?dT?hj?(TrefTref,j?0,j?)
2－9 其中hj?(Tref,j?)是组分j?基于参考温度Tref,j?的生成焓。 虽然能量的标准形式里包括了压力做功和动能项，但在采用segregated solver求解不可压问题时候都可以忽略掉。当然，如果想不忽略它们的作用，可以在define/models/energy x 0中设置。对于可压缩流动问题，在用coupled solvers求解时总是考虑压力做功和动能项。 粘性耗散项是考虑流体中的粘性剪切作用产生的热量。如果用segregated solver求解，默认设置并没有考虑。如果Brinkman数（Br??Ue2k?T性加热一定不能忽略。这时候一定要设置Viscous Heating选项。对于可压缩流动，一般Br>1，如果还用segregated solver求解，一定要考虑粘性加热。如果是coupled solver求解，粘性加热会自动考虑。 Fluent求解焓方程时，组分扩散项都已经包括。用segregated solver求解，如果想不考，?T是系统温度差）大于1时，粘虑该项，可以在组分模型面板（Species Model Panel）中关闭能量扩散项。如果采用了非绝热的PDF燃烧模型，方程中并不明确出现该项，应为导热和组分扩散项合并为一项了。当用coupled solver求解时，能量方程总会考虑该项。
2.1.5 化学反应源项
化学反应源项如下： Sh,reaction0?h0j?????j?M?j????cp,j?dT?Rj?
2－10 Tref,j???Tref其中，hj?是组分j?的生成焓；Rj?是组分j?生成的体积率。对于非绝热PDF燃烧模型，（2－9），生成热定义在总焓中，所以化学反应热不包含在源项中。
2.1.6 固体区域的能量方程
在固体区域，FLUENT采用的能量方程为如下形式： ??t?h???xi(ui?h)???xi(k?T?xi????
2－11 )?q????体积热源。方程左边第二项表示式中，?是密度；h是显焓；k是导热系数；T是温度；q由于固体旋转或者平移运动热传输。方程右边两相分别为固体导热和体积热源。
2.1.7 固体内部导热各向异性的影响
当用segregated solver求解时，FLUENT允许你指定材料的各向导热系数。固体导热各向异性方程形式如下： ??xi(kij?T?xi) 其中，kij是导热系数矩阵。
2.1.8 进口热扩散
进口的净能量输运包括对流和扩散两部分。指定进口温度就可以确定对流部分，但扩散项取决于计算出来的温度场梯度。因此我们不能给定扩散分量或者净能量输运。但在一些问题中，我们更希望能给定净能量输运，而不是给定进口温度。如果用segregated solver求解时，可以在dfine/models/energy中去掉进口能量扩散，从而达到给定净进口能量输运。但是 xi 我们用coupled solver时，不能去掉能量扩散部分。
第二节 计算传热过程中用户输入
如果用FLUENT计算有传热的问题时候，必须击活相关模型和提供热边界条件，并且给出材料物性。这一系列过程如下： 1， 击活能量面板。Define-Models-Energy
2， （对于segregated solver）如果模拟粘性流动过程，而且要考虑粘性加热，击活Viscous Heating；Define-Models-Viscous Heating
3， 定义热边界条件（包括流体进口，出口和壁面）Define-Boundary Conditions。在流动进口和出口要给定温度，但壁面可以有如下边界条件选择： （1） 指定热流量 （2） 指定温度 （3） 对流换热 （4） 外部辐射 （5） 对流换热＋辐射换热
4， 定义材料热物性。Define-Materials. 比热和导热系数都要给出，并且可以用温度函数的形式给出。
2.2.1 温度限制
为了计算的稳定性，FLUENT对计算出来的温度给了范围限制。给定温度限制，一方面是为了计算稳定的需要，同时，真实温度也有其相应的范围。由于给定材料物性不好，或者其它原因，计算出的中间超过了物理应该达到的温度。FLUENT中，给定的最高温度5000K，最小温度1K，如果计算过程中的温度超过这个范围，那么就在这最高温度或最低温度值处锁定。如果你觉得这个限制不合理，你可以自己调节。Solve-control-limits
2.2.2 传热问题求解过程
对于一些简单的传热过程FLUENT的默认设置可以成功进行模拟，但如果要加快你的问题的收敛速度或者提高计算过程的稳定性，下面的一些过程就比较重要了。
2.2.3 松弛因子确定
如果用segregated solver求解能量方程，在solve-controls-solution处定义松弛因子。如果你采用非绝热PDF模型，也必须和通常一样设置包括温度项在内地松弛因子。在求解温度和焓时候，FLUENT默认设置能量方程松弛因子为1。在一些问题里，能量场影响流动场（物性随温度变化，或者有浮力），这时候松弛因子要小些，比如在0.8到1之间。如果流动场和温度场不是耦合的（没有随温度变化的热物性或者浮力影响），松弛因子就可以采用1。 如果我们求解的是焓方程（非绝热PDF燃烧模型），温度需要设置松弛因子。焓的变化中不是所有的都用来计算温度的变化。这对于一些问题，你需要流动场焓变化快，而温度不能变化太快（影响流体热物性太快）的解决很有好处。
xii 2.2.4 组分扩散项
如果用segregated solver求解组分输运方程，如果考虑组分扩散，计算收敛会比较困难。为了提高收敛性，可以在define-models-species处取消对组分扩散的考虑。这时候组分扩散对能量的影响就被忽略了。如果我们选择coupled solver求解，那么组分扩散一定是存在的。
2.2.5耦合和非耦合流动场与温度场计算
如果流动和传热不是耦合的（没有温度变化的热物性或者浮力影响），那么我们可以先求解绝热流动场，然后加进能量方程。这时候可以暂时先关闭动量或者能量方程中的一个，先求解另外的一个。Solve-controls-solution. 如果流动和温度场是耦合的，你可以先求解流动方程，收敛后再击活能量方程，一起求解。需要注意的是，Coupled solver 总是同时求解流动与能量方程。
2.2.6 传热计算结果输出
FLUENT提供了几种传热结果的输出形式。可以以图形的形式输出，也可以用下列参数或者函数形式输出。可以输出的参数包括：静温、总温、静焓、相对总温、壁面温度（内或外表面）、总焓、总焓方差、熵、总能量、内能、表面热流量、表面换热系数、表面Nusselt数和表面stanton数。定义焓或者能量输出的参数跟求解的是可压或不可压问题有关。 可以通过Report-Fluxes给出控制体的每个边界或者通过壁面的总的换热量。必须注意，要检测一下是否能量平衡，这可以检查是否求解已经收敛。 也可以通过Report-surface Integrals-Enthalpy给出某个边界或壁面的总的换热量。焓流??率定义为Q???HV?dA。
2－12 Report-Surface integrals-surface Heat transfer Codf. (wall heat flux)可以给出某个表面的平?hdA均换热系数(h?) 。
2－13 A 第三节 浮力驱动的流动和自然对流
对于混合对流问题，我们定义一个浮力参数为??GrRe2???gh?v2，当?接近或者超过1时，浮力对流动有很强的作用。相反，如果?很小，浮力的作用将可以忽略。而纯粹自然对流问题中，浮力诱发流动强弱我们用另外一个参数来衡量，即Rayleigh数。定义为： Ra?g??TL?/??
2－14 3其中，???1????T8为热膨胀系数；??10k?cp热扩散系数。如果Ra?10，自然对流处于8层流状态，在10?Ra?10为层流到湍流的过渡区域。
2.3.1 Boussinesq模型
对于许多的自然对流问题，采用Boussinesq模型比定义密度是温度的函数有更好的收敛性。该模型在所有求解方程中，认为密度是常数。但是，在动量方程中的浮力项中，密度 xiii
FLUENT商用程序可以模拟许多的工程实际问题，包括可压缩、不可压缩流动，牛顿流体、非牛顿流体，单相、多相流动，有旋、无旋流动，惯性坐标系、非惯性坐标系下的流动，有化学反应、无化学反应的流动问题等。其生成无结构网格的程序把计算复杂几何条件下的流动及传热传质问题变的简单。同时，软件还提供了许多的湍流模型、壁面处理及燃烧、传热模型供针对特定问题选择。用户自定义函数也为改进和完善模型，处理个性化问题和给出更合理的边界条件提供了可能。 本讲义以FLUENT5说明书为主要参考资料，介绍了该软件的基本功能、基本物理模型、湍流模型、湍流模拟的近壁处理及边界条件，并且对燃烧过程的模拟和用户自定义函数做了描述。通过本课程学习，可以掌握和利用FLUENT程序在流体及传热传质等领域进行数值研究。 燃烧模型部分由董刚副教授编译，方海生同学编译了用户自定义函数。基本物理模型，湍流模型及湍流模拟近壁处理及边界条件由刘明侯副教授编译。由于时间非常仓促（一个暑假时间），只能用不完全的内容作为计算流体和传热传质课程的内容。还有些内容来不及加入讲义内，希望以后逐步完善。文字没有很好地校对，一定会由错误、疏漏或不妥的地方，请同学们校正。
第一章， 概述…………………………………………………………………..….…...(1)
第二章， 基本物理模型………………………………………………………..…...….(5) 第一节， 连续和动量方程…………………………………………………….…(5) 第二节， 计算传热过程用户输入………………………………………….……(9) 第三节， 浮力驱动的流动和自然对流…………………………………………(10) 第四节， 有旋和旋转流动问题…………………………………………….……(15) 第五节， 可压流动…………………………………………………..……..……(18) 第六节， 无粘流动……………………………………………………..…..……(19) 第七节， 用户自定义标量输运模型…………………………………..….….…(22)
第三章， 湍流模型……………………………………………………….……..…….(24) 第一节， 前言……………………………………………………………………(24) 第二节， 平均量输运方程…………………………………………….. ……….(25) 第三节， 湍流模型……………………………………………………………….(26) 第四节， 湍流模型算例及其设置……………………………………………….(39) 第四章， 湍流流动地近壁处理…………………………………………………...….(44) 第五章， 边界条件…………………………………………………………………….(52) 第六章， FLUENT中地燃烧模拟…………………………………………..…….….(62) 第一节， 燃烧模拟的重要性……………………………………………..….….(62) 第二节， FLUENT燃烧模拟方法概要…………………………………..…….(62) 第三节， 气相燃烧模型…………………………………………………......…..(63) 第四节， 污染物模型…………………………………………………….………(75) 第五节， FLUENT中燃烧模拟计算的步骤和原则……………………………(77) 第七章， 自定义函数………………………………………………………………….(79) 第一节， 概要…………………………………………………………………….(79) 第二节， 书写UDFs …………………………………………………………….(80) 第三节， 编译连接UDFs………………………………………………………..(98) 第四节， 在FLUENT模型中使用UDFs………………………………………(106) 第五节， UDFs实例…………………………………………………………….(112) 第六节， UDFs的应用………………………………………………………….(128) 附录A，udf. H中的宏解释…………………………………………………………….(146) 附录B，常用C库函数…………………………………………………………………(147)
FLUENT简介
FLUENT是用于计算复杂几何条件下流动和传热问题的程序。它提供的无结构网格生成程序，把计算相对复杂的几何结构问题变得容易和轻松。可以生成的网格包括二维的三角 ii 形和四边形网格；三维的四面体、六面体及混合网格。并且，可以根据计算结果调整网格。这种网格的自适应能力对于精确求解有较大梯度的流场如自由剪切流和边界层问题有很实际的作用。同时，网格自适应和调整只是在需要加密的流动区域里实施，而非整个流动场，因此可以节约计算时间。
程序的结构
FLUENT程序软件包应该包括以下几个部分： 1， FLUENT解法器 2， prePDF，用于模拟PDF燃烧过程 3， GAMBIT，网格生成 4， TGrid，额外的处理器，用于从现有的边界网格生成体网格。 5， Filters(Translators)，转换其它程序生成的网格，用于FLUENT计算。可以接口的程序包括：ANSYS, I-DEAS, NASTRAN,PATRAN等。
基本程序结构示意图。
PDF程序 几何形状或GAMBIT 设置几何形状 生成2D或3D网格 网格 其它软件包，如CAD，CAE等 prePDF
2D或3D网格 FLUENT
网格输入及调整 物理模型 边界条件 流体物性确定 计算 结果后处理 边界网格边界和（或）体网格网格
TGrid 2D三角网格 3D四面体网格 2D和3D混合网格
FLUENT程序的用途
1， 采用三角形、四边形、四面体、六面体及其混合网格计算二维和三位流动问题； 计算过程中，网格可以自适应。 2， 可压缩与不可压缩流动问题；
iii 3， 稳态和瞬态流动问题； 4， 无粘流，层流及湍流问题； 5， 牛顿流体及非牛顿流体； 6， 对流换热问题（包括自然对流和混合对流）； 7， 导热与对流换热耦合问题； 8， 辐射换热； 9， 惯性坐标系和非惯性坐标系下的流动问题模拟； 10，多运动坐标系下的流动问题； 11，化学组分混合与反应； 12，可以处理热量、质量、动量和化学组分的源项； 13，用Lagrangian 轨道模型模拟稀疏相（颗粒，水滴，气泡等）； 14，多孔介质流动； 15，一维风扇、热交换器性能计算； 16，两相流问题； 17，复杂表面形状下的自由面流动；
二维网格：
triangle quadrilateral
tetrahedron 三维网格：
hexahedron prism or wedge pyramid
图1－1，FLUENT的基本控制体形状
用FLUENT程序求解问题的步骤
1， 确定几何形状，生成计算网格（用GAMBIT，也可以读入其它指定程序生成的网格）； 2， 选择2D或3D来模拟计算； 3， 输入网格； 4， 检查网格； 5， 选择解法器； 6， 选择求解的方程：层流或湍流（或无粘流），化学组分或化学反应，传热模型等。确定其它需要的模型如：风扇、热交换器、多孔介质等模型。 7， 确定流体物性；
iv 8， 指定边界条件； 9， 条件计算控制参数； 10，流场初始化； 11，计算； 12，检查结果； 13，保存结果，后处理等。
关于FLUENT解法器的说明
1， FLUENT 2D，
二维单精度解法器 2， FLUENT 3D，
三维单精度解法器 3， FLUENT 2ddp， 二维双精度解法器 4， FLUENT 3ddp， 三维双精度解法器
FLUENT程序启动方法
1， WINDOEWS NT下，点击FLUENT5。 2， 在MS－DOS下，键入命令。
FLUENT 2D/3D/2ddp/3ddp。平行计算命令为：FLUENT 2D/3D/2ddp/3ddp
-t x 。x 是处理器编号。如，我们用3号处理器计算三维双精度问题，命令为： FLUENT 3ddp –t3
FLUENT命令的一般形式为： FLUENT [version] [-help] [options]
FLUENT求解方法的选择
1．非耦合求解 ( Segregated ) 2．耦合隐式求解 ( Coupled Implicit ) 3．耦合显式求解 ( Coupled Explicit ) 非耦合求解方法主要用于不可压缩或压缩性不强的流体流动。耦合求解则可以用在高速可压缩流动。FLUENT默认设置是非耦合求解，但对于高速可压流动，有强的体积力（浮力或离心力）的流动，求解问题时网格要比较密，建议采用耦合隐式求解方法，可以耦合求解能量和动量方程，能比较快地得到收敛解。缺点是需要的内存比较大（是非耦合求解迭代时间的1.5-2倍）。如果必须要耦合求解，但是你的机器内存不够，这时候可以考虑用耦合显式解法器求解问题。该解法器也耦合了动量，能量及组分方程，但内存却比隐式求解方法小。缺点是收敛时间比较长。 这里需要指出的是非耦合求解的一些模型在耦合求解解法器里并不都有。耦合解法器没有的模型包括：多相流模型，混合分数/PDF燃烧模型，预混燃烧模型，污染物生成模型，相变模型，Rosseland辐射模型，确定质量流率的周期性流动模型及周期性换热模型等。 隐式( Implicit )：对于给定变量，单元内的未知值用邻近单元的已知和未知值计算得出。因此，每一个未知值会在不止一个方程中出现，这些方程必须同时解来给出未知量。 显式( Explicit )：对于给定变量，每一个单元内的未知量用只包含已知量的关系式计算得到。因此未知量只在一个方程中出现，而且每一个单元内的未知量的方程只需解一次就可以给出未知量的值。
v 2.2.4 组分扩散项
如果用segregated solver求解组分输运方程，如果考虑组分扩散，计算收敛会比较困难。为了提高收敛性，可以在define-models-species处取消对组分扩散的考虑。这时候组分扩散对能量的影响就被忽略了。如果我们选择coupled solver求解，那么组分扩散一定是存在的。
2.2.5耦合和非耦合流动场与温度场计算
如果流动和传热不是耦合的（没有温度变化的热物性或者浮力影响），那么我们可以先求解绝热流动场，然后加进能量方程。这时候可以暂时先关闭动量或者能量方程中的一个，先求解另外的一个。Solve-controls-solution. 如果流动和温度场是耦合的，你可以先求解流动方程，收敛后再击活能量方程，一起求解。需要注意的是，Coupled solver 总是同时求解流动与能量方程。
2.2.6 传热计算结果输出
FLUENT提供了几种传热结果的输出形式。可以以图形的形式输出，也可以用下列参数或者函数形式输出。可以输出的参数包括：静温、总温、静焓、相对总温、壁面温度（内或外表面）、总焓、总焓方差、熵、总能量、内能、表面热流量、表面换热系数、表面Nusselt数和表面stanton数。定义焓或者能量输出的参数跟求解的是可压或不可压问题有关。 可以通过Report-Fluxes给出控制体的每个边界或者通过壁面的总的换热量。必须注意，要检测一下是否能量平衡，这可以检查是否求解已经收敛。 也可以通过Report-surface Integrals-Enthalpy给出某个边界或壁面的总的换热量。焓流??率定义为Q???HV?dA。
2－12 Report-Surface integrals-surface Heat transfer Codf. (wall heat flux)可以给出某个表面的平?hdA均换热系数(h?) 。
2－13 A 第三节 浮力驱动的流动和自然对流
对于混合对流问题，我们定义一个浮力参数为??GrRe2???gh?v2，当?接近或者超过1时，浮力对流动有很强的作用。相反，如果?很小，浮力的作用将可以忽略。而纯粹自然对流问题中，浮力诱发流动强弱我们用另外一个参数来衡量，即Rayleigh数。定义为： Ra?g??TL?/??
2－14 3其中，???1????T8为热膨胀系数；??10k?cp热扩散系数。如果Ra?10，自然对流处于8层流状态，在10?Ra?10为层流到湍流的过渡区域。
2.3.1 Boussinesq模型
对于许多的自然对流问题，采用Boussinesq模型比定义密度是温度的函数有更好的收敛性。该模型在所有求解方程中，认为密度是常数。但是，在动量方程中的浮力项中，密度 xiii 才随温度变化。(???0)g???0?(T?T0)g，因而用???0(1???T)计算浮力项。这样的近似对密度变化很小的流动问题有较好计算结果。 该模型对封闭区域里的自然对流问题适合，如果模拟温度变化很小的流动场也同样适用。但是，如果计算组分，燃烧或者有化学反应的问题时，该方法不适合。
2.3.2浮力驱动流计算的用户输入
在计算混合对流和自然对流问题时，需要输入如下参数： 1， 求解能量方程；define-models-Energy
2， 击活重量选项面板define-operating conditions；在笛卡儿坐标系的每个方向上设置重力加速度。默认的值是零。 3， 如果采用不可压理想气体法则，注意检查define-operating conditions 中的压力设置到一个合适值（非零）。 4， 取决于你是否采用了Boussinesq模型，还有以下参数需要合适确定： （1） 如果不用Boussinesq模型，则要给定合适的密度；或者在define-materials 里给出密度与温度的函数关系。 （2） 如果采用了Boussinesq模型，给出参考温度（方程中的T0），在define-materials 里给定热膨胀系数。 （3） 如果计算问题包含多个流体物质，可以为单个物质考虑是否选用Boussinesq近似。 ??0gx?ps，条5， 在压力进口与出口边界条件下，应该输入的压力为等效压力p?s件是进口和出口没有外部压力梯度。Define-boundary conditions
6， solve-controls-solution，压力离散化方法选择。如果你采用的是四边形或六面体网格，并且Segregated solver 求解，建议选择PRESTO离散方法。
2.3.3 流体参考密度
如果不用Boussinesq近似假设，流体在动量方程中体积力项中出现参考密度?0，其形?ps??0gx，根据流体静力学平衡，我式为：(???0)g，根据FLUENT定义的压力p?s们有： ?ps?x??g
2－15 上式可以变成：
?p?s?x?(???0)g
2－16 所以，定义流体参考密度是十分重要的。 FLUENT默认设置中，计算所有网格的平均密度为参考密度。在一些计算问题里，也许你给出一个参考密度要比自动计算一个参考密度要好。例如，你计算的是自然对流问题，ps，而且有个压力边界条件，了解上面方程中p?s的定义是十分重要的。虽然你知道确切的但你还需要知道参考密度?0，而根据ps的值来计算p?s。因此，你必须明确给出参考密度，而不能等计算机自动计算。而有时候给定参考密度虽然对最后结果没有太多影响，但对收敛速度有明显作用。这时候通常用体平均密度作为参考密度。这里需要指出的是，如果你采用的是Boussinesq近似方法，不需要用参考密度，因此不需要特别给出。
对于高Rayleigh 数流动问题的求解，需要注意以下一些问题。需要指出的是对一些层流问题而言，高Rayleigh 数流动有可能没有稳态解存在。 当求解高Rayleigh 数（Ra?108）流动问题时，根据下列步骤将能得到最好结果。 第一步是求稳态近似结果： ? 选用First-order scheme，在小Rayleigh数下求得稳态解。（可以通过变化重力加速度的方法减少Ra数（比如从9.8降低到0.098，Ra数就降低了两个数量级） ? ? 用小Ra数的收敛解为初始值，求解高Ra数下的解。 得到收敛解后，可以换higher-order scheme 继续求解。 第二步是求与时间相关的稳定解： ? 用前面的稳态解为初始条件，在相同或略小Ra数下求解。 ? 估计时间常数 ??LU~L2?(PrRa)?1/2?Lg??TL。其中，L和U是长度和速度? 尺度，采用的时间步长为：?t??/4，如果时间步长比?t大，有可能不收敛。 求解过程中会有频率为f??0.05~0.09振荡，衰减后就达到稳态解。?是上面的求出的时间常数，f是振荡频率（Hz）。通常需要超过5000步才能得到稳定解。
需要进一步指出的是除非我们采用了Boussinesq近似，上面方法不能用于封闭区域的流动问题，只能用于有进口和出口的流动问题。
2.3.4周期性流动与换热
如果我们计算的流动或者热场有周期性重复，或者几何边界条件周期性重复，就形成了周期性流动。FLUENT可以模拟两类周期性流动问题。第一，无压降的周期性平板问题（循环边界）；第二，有压降的周期性边界导致的完全发展或周期性流向流动问题（周期性边界）。 流向周期性流动模拟的条件： ? 流动是不可压的； ? 几何形状必须是周期性平移； ? 如果用coupled solver求解，则只能给定压力阶跃；如果是Segregated solver，可以给定质量流率或者压力阶跃。 ? ? ? 周期性流动中不能考虑进口和出口有质量差，也不考虑过程中的额外源项或者稀疏相源项。 只能计算进口出口没有质量流率变化的组分问题。但不能考虑化学反应。 不能计算稀疏相或者多相流动问题。 如果在这过程中计算有换热问题，则还必须满足以下条件： ? 必须用segregated solver 求解 ? 热边界条件必须是给定热流率或者给定壁面温度。对于一个具体的问题，热边界条件只能选择一个，而不能是多热边界条件问题。对于给定温度热边界条件，所有壁面的温度必须相同（不能有变化）。对于给定热流率边界条件，不同壁可以用不同值或曲线来模拟。 ? ? 对于有固体区域的问题，固体区域不能跨越周期性平板。 热力学和输运特性（热容，热导系数，粘性系数，密度等）不能是温度的函数（所以不能模拟有化学反应流动问题）。但输运特性（有效导热系数，有效粘性系数）可以随空间有周期性变化，因此可以对有周期性湍流输运特性不同的流动问题有模拟能力。
2.3.5 计算流向周期性流动问题的步骤：
通常，可以先计算周期性流动到收敛，这时候不考虑温度场。下一步，冻结速度场而计算温度场。步骤如下： 1， 建立周期性边界条件网格； 2， 输入热力学和分子输运特性参数； 3， 指定周期性压力梯度或者确定通过周期性边界的质量流量； 4， 计算周期性流动场。求解连续，动量（湍流量）方程。 5， 指定热边界条件（等温或者给定热流密度）； 6， 给定进口体平均温度； 7， 求解能量方程（其它方程不求解，只求解能量方程），得到周期性温度场。
2.3.5.1流向周期性流动理论 周期性速度定义 ?对于位置矢量r，周期性速度定义为： ?????
u(r)?u(r?L)?u(r?2L)?...
2－17 ?????
v(r)?v(r?L)?v(r?2L)?...
2－18 ?????w(r)?w(r?L)?w(r?2L)?...
2－19 ?其中，L是计算区域内周期性长度矢量。
2.3.5.2 周期性流向周期压力 上面方程中压力不是周期性的，而压力降是周期性的。即： ????????p?p(r)?p(r?L)?p(r?L)?p(r?2L)?...
2－20 如果选择coupled solver，压降?p是常数；但对Segregated solver方法求解，计算区域?L?内的压力梯度可以分解为两部分：梯度的周期性分量，?~p(r)，梯度的线性变化量，??。L即： ???L?p(r)?????~p(r)
L??周期性压力（~是但压力减去线性变化压力得到的值。压力的线性变化分量?r是p(r)）导致一个力作用于流体的动量方程。由于?并不能事先知道，只有通过在给定区域积分求出的流量与给定的的质量流量一致，才能确定。在压力修正的SIMPLE, SIMPLEC和PISO运算法则中都出现?，根据计算得到的质量流量与实际流量之差得到。
2.3.5.3 用Segregaed solver求解流向周期性流动的用户设置 要计算给定质量流量或压降的空间周期性流动问题，首先要建立计算网格。网格要求有平移周期性边界条件。在define-periodic conditions 面板上第一如下参数： 1， 决定是给定质量流量还是压力梯度。（Specify mass flow, Specify pressure Gradient） 2， 给出质量流量或压力梯度。如果你选择了给定质量流量，你也可以给出压力梯度（假定值），但不是必须的。如果给定值合理，可以加速流动场收敛速度。对计算结果也没有影响。对于轴对称问题，质量流量是2?弧度的质量流量。 3， 通过定义流向的X，Y，Z点来（或者二维时候的X，Y）定义流动方向。流动方向 xvi 是从初始点到指定点的方向上，该方向必须和周期性平移边界平行。 4， 如果确定了质量流量，FLUENT需要计算出压力梯度?。可以通过改变Relaxation Factor 和Number of Iterations（压力修正方程迭代次数）或者给出估计的?用以控制计算过程，这些都在define-periodicty conditions 面板上做。
2.3.5.4用Coupled solver方法求解周期性流动场 网格处理和前面类似，参数设置方面，需要： 1， define-boundary condition；
2， 给定周期性压降；
2.3.6对周期性流动问题计算过程的监视 如果是给定质量流量问题，可以在计算过程中观察压力降?（Statistic Monitors panel ）
2.3.7 给定温度边界条件下的周期性传热计算
要计算周期性流动中的传热问题，必须选用Segregated solver 方法。对于给定壁面温度条件，现在我们假设壁面温度是常数，则温度为：[L119] ???T(r)?Twall,inletTbulk?Twall
2－22 ???T?V?dA,inlet进口体温定义为：Tbulk?A， 这里的??，积分在周期性边界的进口上（A）???V?dAA在计算区域的L长度上满足周期性变化条件。
2.3.8 壁面温度确定周期性传热设置
1， 求解能量方程（define-models-Energy ） 2， 给出壁面温度。所有固体边界的温度相同，所有边界（除了周期性边界）必须用该温度的壁面封闭。Define-boundary conditions
3， 定义固体区域和固体温度。固体区域内不能有热源。Define-boundary conditions
4， 设定流体物性（密度，热容，粘性系数，导热系数等），但不能随温度变化。Define-Materials
5， 给定周期性进口上游的进口体平均温度（不能和壁面温度相同）define-periodic conditions
有了以上设置后就可以迭代求解了。最好的方法是先求解流场，等流场收敛后，再冻结流场只求解温度场。再给定处温度场时，温度值在壁温与进口体平均温度之间。 可以通过监视进出口平均体温度之比??monitors-per/bulk-temp-ratio
Twall?TbulkTwall?Tbulk,inlet,exit。Define-statistic 第四节 有旋与旋转流动问题
xvii 是从初始点到指定点的方向上，该方向必须和周期性平移边界平行。 4， 如果确定了质量流量，FLUENT需要计算出压力梯度?。可以通过改变Relaxation Factor 和Number of Iterations（压力修正方程迭代次数）或者给出估计的?用以控制计算过程，这些都在define-periodicty conditions 面板上做。
2.3.5.4用Coupled solver方法求解周期性流动场 网格处理和前面类似，参数设置方面，需要： 1， define-boundary condition；
2， 给定周期性压降；
2.3.6对周期性流动问题计算过程的监视 如果是给定质量流量问题，可以在计算过程中观察压力降?（Statistic Monitors panel ）
2.3.7 给定温度边界条件下的周期性传热计算
要计算周期性流动中的传热问题，必须选用Segregated solver 方法。对于给定壁面温度条件，现在我们假设壁面温度是常数，则温度为：[L119] ???T(r)?Twall,inletTbulk?Twall
2－22 ???T?V?dA,inlet进口体温定义为：Tbulk?A， 这里的??，积分在周期性边界的进口上（A）???V?dAA在计算区域的L长度上满足周期性变化条件。
2.3.8 壁面温度确定周期性传热设置
1， 求解能量方程（define-models-Energy ） 2， 给出壁面温度。所有固体边界的温度相同，所有边界（除了周期性边界）必须用该温度的壁面封闭。Define-boundary conditions
3， 定义固体区域和固体温度。固体区域内不能有热源。Define-boundary conditions
4， 设定流体物性（密度，热容，粘性系数，导热系数等），但不能随温度变化。Define-Materials
5， 给定周期性进口上游的进口体平均温度（不能和壁面温度相同）define-periodic conditions
有了以上设置后就可以迭代求解了。最好的方法是先求解流场，等流场收敛后，再冻结流场只求解温度场。再给定处温度场时，温度值在壁温与进口体平均温度之间。 可以通过监视进出口平均体温度之比??monitors-per/bulk-temp-ratio
Twall?TbulkTwall?Tbulk,inlet,exit。Define-statistic 第四节 有旋与旋转流动问题
xvii 有旋流动在一些燃烧器里用以增加流体滞驻时间和驻涡，从而稳定燃烧。旋转流动在透平机等实际过程中也应用广泛。有旋与旋转流动问题可以分为一下5类：1，轴对称有旋或旋转流动；2，完全三维有旋或旋转流动；3，一个坐标系旋转的流动；4，多坐标系旋转流动；5，流动需要滑移网格。 1， 轴对称有旋或旋转流动 流动没有周向梯度（但可以有周向速度）。周向动量方程为： ??t(?w)?1?r?x(r?uw)?1?r?r(r?vw)?1???w? r??r?x??x??1??3??w??r?2?r?rr?r?vw??
2－23 ????r??
?其中，x 和r分别是轴向和径向坐标；u, v, w分别是轴向，径向和周向速度分量。 2， 三维有旋流动 如果几何形状发生变化或者周向有流动梯度，则需要用三维模拟。对于三位问题，和解二维问题类似，没有特别的输入或者求解步骤。但是，定义速度进口条件时候需要用柱坐标系，并且，需要在求解时渐渐增加旋转速度。 3， 有旋转坐标系的流动 如果有旋转边界（螺旋桨等）必须用旋转坐标系来求解该类问题。如果有多个旋转边界，还需要多个旋转坐标系来求解。
2.4.1有旋或旋转流动的物理概念 对于有旋流动，轴向动量（rw 或r2?）往往会导致自由涡流动。周向速度w 随半径增加而增加，在粘性力起主导作用的r ?0的区域，周向速度近似为零。龙卷风就是一个很好的例子。 对于自由涡流动，流体周向运动的离心力与径向压力梯度相等：
??rr壁面旋转驱动的流动中，w/r或?是常数
有旋流动的湍流模拟
许多流动明显具有旋流（龙卷风，旋转射流等），必须考虑选用FLUENT提供的比较高级的模型，如RNG k-?模型、可实现k-?模型或者雷诺应力模型。具体选择哪个模型，取决于流动的旋流强度（旋流数）。旋流数定义为轴向与周向动量比： ???rwv?dA?
2－24 S??R?uv?dAR是水力学半径。 ?p?w2 对于弱旋和中等旋度流动问题（S0.5），必须选用雷诺应力模型。只有雷诺应力模型才能模拟该流动中的强的各向异性影响。 有旋流动问题的第二关键问的是边界条件问题。因为场的模拟好坏，主要取决于采用的模型，而壁面也参与了涡旋（涡量）的产生（由于压力梯度产生的二次流或涡流），采用非平衡壁面函数可以得到比较好的模拟结果，因为该壁面法则的平均流动速度对压力梯度敏 xviii 感。
2.4.2有旋与旋转流动问题的网格设置
坐标系限制：对于轴对称问题，旋转轴必须是x 轴，网格必须在y=0线以上。 除了上面的注意事项，对于有旋和旋转问题，网格划分还必须保证问题求解有足够分别率。特别是旋转流动中，边界层很薄，FLUENT需要在靠近旋转边界的地方网格比较细。除此之外，对于有旋流动问题，周向速度梯度比较陡（近中心线区域为典型的自由涡流动），因此要保证有很好的求解效果，网格要求比较密。
2.4.3轴对称有旋流动设置 1， 求解周向动量。Define-models-solver-Axisymmetric swirl
2， 给出进口或者壁面的旋转或者有旋速度分量，r?；define-boundary conditions (对于旋转轴，选用axis boundary边界条件)。
求解有旋或者旋转问题时的困难在于动量方程之间的耦合。如果是旋转很强的流动，会导致比较大的径向压力梯度，并驱动流体在轴向和径向的流动；动量之间的强耦合作用会导致求解过程中的不稳定性。如果要得到一个好的收敛解，需要有一定技巧。 求解步骤为： 1，（segregated solver only）如果用的是四边形和六面体网格，选用PRESTO模型（solution controls panel ） 2，为了求解大压力梯度和轴向速度梯度，网格必须足够精细。 3，(segregated solver only）改变松弛因子。径向和轴向速度为0.3-0.5，周向速度为0.8-1。 4，（segregated solver only）一步一步求解： （1） 如果包含inflow/outflow 的问题，那么先求解无旋流动。即用Axisymmetric ,而不选用Axisymmetric swirl option。并且不设置任何有旋边界条件。计算的结果作为有旋流动的初始值。 （2） 击活Axisymmetric Swirl option，设置有旋和旋转边界条件 （3） 只求解周向速度动量方程。让旋转的边界条件扩散到整个流场。如果是求解的湍流场，这适合湍流方程也应该同时求解。 （4） 关闭（冻结）求解周向动量方程，再求解连续和其它动量方程。如果是求解的湍流场，湍流量方程也同时求解。 （5） 松弛因子给定合适的值，同时求解所有方程。
除了上面描述的以为，如果求解的是有换热问题，可以先求解绝热流动场。如果求解的是湍流问题，也可以先计算层流，然后假如湍流模型继续计算。该方法对segregated or coulpled 求解都适合。 5，如果可能的话，先用小的旋转速度或者进口旋流速度进行计算。然后再增加到要计算的值。 （1） 再进口旋转速度或者旋转边界条件上给小的旋转速度，例如给需要计算值的10％。 （2） 在上面给定的条件下求解（可以用上面一步一步的求解方法） （3） 存储上面的初步结果。 （4） 更改初始和边界条件，增加旋转速度（也许是第一次的2倍）。 （5） 用上面的计算结果做初始值，开始新一轮计算。
xix （6） 进一步增加旋流速度，重复上面过程4－5，直到需要计算的旋流速度值。
第五节 可压流动
当气流速度很大，或者流场压力变化很大时，流体就收到了压速性的影响。 马赫数定义为：M?u/c，这里，c 是气体的音速，为：c? ?RT，??cp/cv比热比。当马赫数小于1时，流动为亚音速流动；当马赫数远远小于1（如M<0.1）时，流体的可压速性及压力脉动对密度变化影响都可以忽略。当马赫数接近1时候（跨音速），可压速性影响就显得十分重要了。如果马赫数大于1，流体就变为超音速流动。FLUENT对于亚音速，跨音速以及超音速等可压流动都有模拟能力。
2.5.1可压速流动的基本概念
描述可压速流动我们都离不开两个参数，总压p0和总温T0。对于理想气体，它们与静压与静温之间的关系为： ?p0ps??12???1???1?M?
2－25 2???1?T0Ts??12M
2－26 2上述关系成立的条件是等熵过程。
2.5.2 可压速流动的基本方程
FLUENT提供的标准连续和动量方程就可以描述可压速流动问题，除了下面会介绍的可压速流动处理以外，不需要其它特别的物理模型。需要指出的是，求解可压速问题，一定要求解能量方程。如果用segregated solver 求解，一定要考虑粘性耗散项（粘性加热）。 对于可压速流动，理想气体方程为： ??(pop?p)/RTs
2－27 这里，pop是运行条件里确定的运行压力，p是当地静压，R是通用气体常数，根据输入的气体分子量计算；Ts用能量方程求得。
2.5.3 可压速流动模拟的参数输入：
1，设定运行压力pop。Define-operating conditions ； 2，求解能量方程； 3，（Segregated solver only ）如果模拟的是湍流流动问题，考虑粘性耗散。Define- models- viscous.；(耦合求解，不需要，因为耦合求解自动考虑粘性耗散) 4，材料面板设置；Define-materials ； a) 选择理想气体； b) 定义物性（比如，分子量，导热系数等）。 5，设定边界条件； （1）Flow inlets
(a )压力进口：进口总温，总压；对于超音速流动，静压；
（b）质量进口：进口质量流率和总温。
xx （2）Flow exits
压力出口：出口静压（如果出口是超音速，可以忽略）； 需要特别指出的是，输入的边界条件中，无论是静压还是总压，都必须是表压（与前面给定的pop的差）。进口温度给定必须是总温（滞止温度），不是静温。
2.5.4 可压速流动求解注意事项
求解可压速问题的困难在于流体速度，密度，压力和能量的高度耦合。这样的耦合会导致求解过程中的不稳定性。因此要得到收敛解必须采取一定的手段。另外，超音速流动的激波也会导致求解的不稳定性。 下面介绍较好求解步骤： 1，（Segregated solver only ）速度的松弛因子调低（0.2-0.3）； 2，（Segregated solver only ）压力松弛因子用0.1，采用SIMPLE算法。（可压速流动不能采用SIMPLEC或者PISO） 3，对压力温度设置合适的限制条件。Solution limits。特别注意的是给定的压力和温度初始值要合理，如果给定的限制条件得到的收敛结果不好，可以更改继续计算，直至取得满意结果。 4， 有时候可以考虑用无粘收敛结果做初始场有好的收敛效果。
第六节 无粘流动
无粘流动是忽略了流体粘性作用，特别在大雷诺数流动中，惯性力起主导作用。在高速气体动力学里有比较多的应用。在这类流动中，压力作用在固体上的力比粘性力大很多，我们可以做无粘分析，快速得到作用在物体上的主要力的大小。然后，我们也可以假如粘性（包括湍流粘性）来评估对物体阻力或升力的影响。 另外一方面，也许我们求解的问题力粘性力不能忽略，但我们也可以先做无粘分析，用无粘结果作为有粘计算的初始值，这一方法，特别对一些复杂的流动，往往有好的收敛效果。
2.6.1 方程描述
无粘流动，我们求解的是Euler方程。连续方程与层流连续方程相同，而动量与能量方程由于没有考虑分子扩散而得到了简化。
连续方程：
???t???xi(?ui)?Sm
该方程是通用形式，对可压与不可压问题都适合。源项是稀疏相对连续相的贡献，也可以是用户自定义源项。 对于二维轴对称问题，连续方程可以写成： ?????v?(?u)?(?v)??Sm
2－29 ?t?x?rr动量守恒方程：
xxi （2）Flow exits
压力出口：出口静压（如果出口是超音速，可以忽略）； 需要特别指出的是，输入的边界条件中，无论是静压还是总压，都必须是表压（与前面给定的pop的差）。进口温度给定必须是总温（滞止温度），不是静温。
2.5.4 可压速流动求解注意事项
求解可压速问题的困难在于流体速度，密度，压力和能量的高度耦合。这样的耦合会导致求解过程中的不稳定性。因此要得到收敛解必须采取一定的手段。另外，超音速流动的激波也会导致求解的不稳定性。 下面介绍较好求解步骤： 1，（Segregated solver only ）速度的松弛因子调低（0.2-0.3）； 2，（Segregated solver only ）压力松弛因子用0.1，采用SIMPLE算法。（可压速流动不能采用SIMPLEC或者PISO） 3，对压力温度设置合适的限制条件。Solution limits。特别注意的是给定的压力和温度初始值要合理，如果给定的限制条件得到的收敛结果不好，可以更改继续计算，直至取得满意结果。 4， 有时候可以考虑用无粘收敛结果做初始场有好的收敛效果。
第六节 无粘流动
无粘流动是忽略了流体粘性作用，特别在大雷诺数流动中，惯性力起主导作用。在高速气体动力学里有比较多的应用。在这类流动中，压力作用在固体上的力比粘性力大很多，我们可以做无粘分析，快速得到作用在物体上的主要力的大小。然后，我们也可以假如粘性（包括湍流粘性）来评估对物体阻力或升力的影响。 另外一方面，也许我们求解的问题力粘性力不能忽略，但我们也可以先做无粘分析，用无粘结果作为有粘计算的初始值，这一方法，特别对一些复杂的流动，往往有好的收敛效果。
2.6.1 方程描述
无粘流动，我们求解的是Euler方程。连续方程与层流连续方程相同，而动量与能量方程由于没有考虑分子扩散而得到了简化。
连续方程：
???t???xi(?ui)?Sm
该方程是通用形式，对可压与不可压问题都适合。源项是稀疏相对连续相的贡献，也可以是用户自定义源项。 对于二维轴对称问题，连续方程可以写成： ?????v?(?u)?(?v)??Sm
2－29 ?t?x?rr动量守恒方程：
xxi ??t(?ui)???xj(?uiuj)???p?xi??gi?Fi
2－30 对于二维轴对称流动，上面方程可以写成： ??t??t(?u)?(?v)?1?r?x1?r?x(r?uu)?(r?uv)?1?r?r1?r?r(r?vu)??(r?vv)???p?x?p?r?Fx
2－32 能量守恒方程： ??t(?E)???xi(ui(?E?p))????xi?hj?j?Jj??Sh
2.6.2 求解无粘流动设置：
1， 设置无粘模型；Define-models-viscous ； 2， 设置边界条件和流体物；define-boundary condition /Materials ； 3， 求解，检查结果。
由于无粘问题往往是伴随高速流动，需要给出小的动量方程松弛因子。
有时间相关项的问题数值模拟
FLUENT可以解决的非定常问题包括：1，涡旋脱落等其它周期性现象；2，压缩填充与抽空问题；3，瞬时热导问题；4，瞬时化学混合与反应等。 求解时间相关项对一些定常问题求解时的收敛稳定性有时候有帮助，比如，求解Ra数比（层流向湍流过渡区域）较大的自然对流问题。对于一些问题，计算时间相关方程，是可以得到定常解的。
时间项离散
在FLUENT中，与时间相关项必须在时间和空间上离散。在空间的离散过程与求解稳态问题一样。 变量?随时间变化的过程通常可以写成如下形式： ???t?F(?)
2－34 其中，F是包含空间离散的函数。 如果时间导数项用向后差分，则一阶精度形式为： ?n?1?F(?)
2－35 ?t而二阶精度的离散形式为： n?1??n3??4?n??n?12?t时间项确定后，下面就该确定F(?)了。
xxii 一阶精度的隐式格式：二阶精度的隐式格式：?n?1??n?t3?n?1?F(?nn?1)
2－37 ?F(?n?1?4???n?12?t隐式格式的好处是对时间步无条件稳定。 )
2－38 一阶精度的显式格式：二阶精度的显式格式：?n?1??n?t3?n?1?F(?)
2－39 nn?4???n?12?t显式格式的时间步选择选择比较关键。显式格式都用来捕捉移动波，该格式比较精确，而且节约计算时间。但是下列情况不能用显式格式：1，不能用于segregated or coupled implicit solver.而只能用于coupled explicit solver .2，不能用于不可压速流动。不可压速流动必须在每个时间步长里积分收敛。
2－40 n 2.6.3 求解非定常问题的步骤：
1，击活求解非定常问题面板：define-models-solver。一阶精度的隐式格式对很多问题都适合。显式格式只能用于coupled explicit solver，用于捕捉移动波。 2，定义相关的模型和边界条件。自定义函数边界条件可以是跟时间相关。 3，如果用segregated solver, 则选用PISO来处理压力速度耦合。Solve-controls-solution。 4，监视计算结果。Solve-monitors-Statistic。 5，设置初始条件。Solve-initialize-initizlize。也可以读入稳态计算结果做初始值。File-read-data。
6，存贮中间结果。File-write-autosave。指定文件名字和存贮频率。 7，设定时间相关项处理方法 （1）如果采用一阶或者二阶隐式，需要设置如下参数： ? ? 每个时间步最多积分次数。 时间步长?t。虽然隐式格式对时间步长没有太多限制，但是对于我们模拟的一些瞬态问题，还应该把这个时间步定义为流动最小时间常数小一个数量级。最好的选择时间步的方法是看需要多少积分步才能收敛。理想的积分步是每个时间步长上为10－20次。如果需要迭代的次数很多，证明时间步太大了，需要调小点。如果迭代一两次就收敛了，证明时间步长可以调大点。对于周期下性的流动问题，时间步可以根据周期来确定。比如，一个周期可以分成20个积分步长。 （2）如果选择显式非稳定形式，则： ? ? 定义求解变量solve-controls-solution (必须没有多重网格，没有残差光滑，而且必须courant number of 1。 选择积分次数，求解。 8，存储结果，以便处理和进一步计算。File-write-data
第七节 用户自定义标量输运模型
Fluent可以在求解质量组分分数输运方程的同时，还可以求解用户自定义的任意标量输运方程。这样，我们可以求解磁流体动力学问题（MHD）。在这个过程中，流动的导电流体 xxiii 产生的磁场可以当做用户自定义标量来求解。流体产生的磁场又会对流体产生阻尼作用，这可以用用户自定义源项来处理。 对于任意的标量，?k，根据计算对流通量的不同方法，FLUENT提供了三个选择： 1，如果对流通量不需要求解，求解的方程为： ??k??(?k)?S ?
k=1, 2, ……..N
2－41 ?xi?xik 其中，?k 是扩散系数，S ?是第N个标量的源项。 k 2，如果对流通量需要计算，则求解的方程为： ??k?(?u i ?k ??k)?S ?，k=1, 2, ……..N
2－42 ?xi?xik 3，如果用户用自定义函数来确定对流通量，则方程为： ??k?(F i ?k ??k)?S ?，k=1, 2, ……..N
其中，F i 为自定义函数。 在FLUENT中，用户自定义标量只能在流体控制体中求解，而不能在固体控制体中求解。
用户自定义标量求解步骤： 1， 指定用户自定义标量数。Define-models-user-defined scalars
2， 确定对流通量确定形式。None (不需要计算)，mass flow rate (计算对流通量)，user defined（用户定义对流通量）。只能用一个确定对流通量的方法。 3， 定义自定义函数的边界条件，进口和出口边界条件。可以给固定值，也可以给出通量。 4， 如果需要用用户自定义源项，则：define-boundary conditions,给出自定义源项设置。 5， 给定初始值，求解 6， 结果检查。主要检查该标量的值scalar-n 和扩散系数diffusion coef.
Of scalar-n。
第三章 湍流模型 第一节 前言
湍流流动模型很多，但大致可以归纳为以下三类： 第一类是湍流输运系数模型，是Boussinesq于1877年针对二维流动提出的，将速度脉动的二阶关联量表示成平均速度梯度与湍流粘性系数的乘积。即： ?u1?u2???t??u1
3－1 ?x2推广到三维问题，若用笛卡儿张量表示，即有： ??ui?uj?????uiuj??t???x?xij??2???k?
3－2 ij?3?模型的任务就是给出计算湍流粘性系数?t的方法。根据建立模型所需要的微分方程的数目，可以分为零方程模型（代数方程模型），单方程模型和双方程模型。 第二类是抛弃了湍流输运系数的概念，直接建立湍流应力和其它二阶关联量的输运方程。
xxiv 第三类是大涡模拟。前两类是以湍流的统计结构为基础，对所有涡旋进行统计平均。大涡模拟把湍流分成大尺度湍流和小尺度湍流，通过求解三维经过修正的Navier-Stokes方程，得到大涡旋的运动特性，而对小涡旋运动还采用上述的模型。 实际求解中，选用什么模型要根据具体问题的特点来决定。选择的一般原则是精度要高，应用简单，节省计算时间，同时也具有通用性。 FLUENT提供的湍流模型包括：单方程（Spalart-Allmaras）模型、双方程模型（标准κ-ε模型、重整化群κ-ε模型、可实现(Realizable)κ-ε模型）及雷诺应力模型和大涡模拟。
包含更多 物理机理
Zero-Equation Models RANS-based
models One-Equation Models
Spalart-Allmaras Two-Equation Models
Standard k-? 每次迭代
RNG k-??计算量增加 ?????Realizable k-??
FLUENT提Reynolds-Stress Model 供的模型选 择 Large-Eddy Simulation
Direct Numerical Simulation
湍流模型种类示意图 xxv 第三类是大涡模拟。前两类是以湍流的统计结构为基础，对所有涡旋进行统计平均。大涡模拟把湍流分成大尺度湍流和小尺度湍流，通过求解三维经过修正的Navier-Stokes方程，得到大涡旋的运动特性，而对小涡旋运动还采用上述的模型。 实际求解中，选用什么模型要根据具体问题的特点来决定。选择的一般原则是精度要高，应用简单，节省计算时间，同时也具有通用性。 FLUENT提供的湍流模型包括：单方程（Spalart-Allmaras）模型、双方程模型（标准κ-ε模型、重整化群κ-ε模型、可实现(Realizable)κ-ε模型）及雷诺应力模型和大涡模拟。
包含更多 物理机理
Zero-Equation Models RANS-based
models One-Equation Models
Spalart-Allmaras Two-Equation Models
Standard k-? 每次迭代
RNG k-??计算量增加 ?????Realizable k-??
FLUENT提Reynolds-Stress Model 供的模型选 择 Large-Eddy Simulation
Direct Numerical Simulation
湍流模型种类示意图 xxv
第二节 平均量输运方程
雷诺平均就是把Navier-Stokes方程中的瞬时变量分解成平均量和脉动量两部分。对于速度，有： ui?ui?ui?
3－3 其中，ui和ui?分别是平均速度和脉动速度（i=1,2,3） 类似地，对于压力等其它标量，我们也有： ??????
3－4 其中，?表示标量，如压力、能量、组分浓度等。 把上面的表达式代入瞬时的连续与动量方程，并取平均（去掉平均速度ui上的横线），我们可以把连续与动量方程写成如下的笛卡儿坐标系下的张量形式： ????(?ui)?0
3－5 ?t?xi?DuDti?uj?ul??????ui2???????????ij??ui?u?j
????xi?xj???x?x3?x?xjil??j????p??上面两个方程称为雷诺平均的Navier-Stokes（RANS）方程。他们和瞬时Navier-Stokes方程有相同的形式，只是速度或其它求解变量变成了时间平均量。额外多出来的项??ui?u?j是雷诺应力，表示湍流的影响。如果要求解该方程，必须模拟该项以封闭方程。 如果密度是变化的流动过程如燃烧问题，我们可以用法夫雷（Favre）平均。这样才可以求解有密度变化的流动问题。法夫雷平均就是出了压力和密度本身以外，所有变量都用密度加权平均。变量的密度加权平均定义为： ~????/?
3－7 符号～表示密度加权平均；对应于密度加权平均值的脉动值用???表示，即有：~???????。很显然，这种脉动值的简单平均值不为零，但它的密度加权平均值等于零，即： ????0， ?????0
Boussinesq近似与雷诺应力输运模型
为了封闭方程，必须对额外项雷诺应力??ui?u?j进行模拟。一个通常的方法是应用Boussinesq假设，认为雷诺应力与平均速度梯度成正比，即：
??ui?u?j??t??2?ui??(?k??)?ij
3－8 t?3?xi?Boussinesq假设被用于Spalart-Allmaras单方程模型和k??双方程模型。Boussinesq近??ui?uj???x?xij?似的好处是与求解湍流粘性系数有关的计算时间比较少，例如在Spalart-Allmaras单方程模型中，只多求解一个表示湍流粘性的输运方程；在k??双方程模型中，只需多求解湍动能k和耗散率ε两个方程，湍流粘性系数用湍动能k和耗散率ε的函数。Boussinesq假设的缺点是认为湍流粘性系数?t是各向同性标量，对一些复杂流动该条件并不是严格成立，所以具有其应用限制性。 另外的方法是求解雷诺应力各分量的输运方程。这也需要额外再求解一个标量方程，通常是耗散率ε方程。这就意味着对于二维湍流流动问题，需要多求解4个输运方程，而三维湍流问题需要多求解7个方程，需要比较多的计算时间，对计算机内存也有更高要求。
xxvi 在许多问题中，Boussinesq近似方法可以得到比较好的结果，并不一定需要花费很多时间来求解雷诺应力各分量的输运方程。但是，如果湍流场各向异性很明显，如强旋流动以及应力驱动的二次流等流动中，求解雷诺应力分量输运方程无疑可以得到更好的结果。
第三节 湍流模型
3.3.1 单方程（Spalart-Allmaras）模型 ~Spalart-Allmaras模型的求解变量是?，表征出了近壁（粘性影响）区域以外的湍流运~动粘性系数。?的输运方程为： 1????G???Dt??~??xj?~D?~????~?(????)?xj??~???????Cb2????x?j???????Y
3－9 ?????其中，G?是湍流粘性产生项；Y?是由于壁面阻挡与粘性阻尼引起的湍流粘性的减少；??~和Cb2是常数；ν是分子运动粘性系数。 湍流粘性系数用如下公式计算： ~f ?t????1其中，f?1是粘性阻尼函数，定义为：f?1??333??C?1，并且??~??。 湍流粘性产生项，G?用如下公式模拟： ~~
3－10 G??Cb1?S?~其中，S?S?~?2kd2f?2，而f?2?1??1??f?1。其中，Cb1和k是常数，d是计算点到壁面的距离；S??ij2?ij?ij。?ij定义为： ??
3－11 ???ui1??uj???2??xj??xi由于平均应变率对湍流产生也起到很大作用，FLUENT处理过程中，定义S为： S??ij?Cprodmin(0,Sij??ij)
3－12 其中，Cprod?2.0，?ij?Sij?ui1??uj???2??xj??xi?ij?ij，Sij?2SijSij，平均应变率Sij定义为： ??
3－13 ??在涡量超过应变率的计算区域计算出来的涡旋粘性系数变小。这适合涡流靠近涡旋中心的区域，那里只有“单纯”的旋转，湍流受到抑止。包含应变张量的影响更能体现旋转对湍流的影响。忽略了平均应变，估计的涡旋粘性系数产生项偏高。 湍流粘性系数减少项Y?为： ~2???Y??Cw1?fw??
3－14 ?d?其中，fw6?1?Cw?3?g?66??g?Cw3?61/6
g?r?Cw2(r?r)
3－16 ~?r?~22
3－17 Skd xxvii ~其中，Cw1，Cw2，Cw3是常数，S?S?~?kd22f?2。在上式中，包括了平均应变率~对S的影响，因而也影响用S计算出来的r。 上面的模型常数在FLUENT中默认值为：Cb1?0.1335，Cb2?0.622，??~?2/3，C?1?7.1，Cw1?Cb1/k2?(1?Cb2)/??~，Cw2?0.3，Cw3?2.0，k?0.41。
~在壁面，湍流运动粘性?设置为零。当计算网格足够细，可以计算层流底层时，壁面切应力用层流应力－应变关系求解，即： ?u?yu?
3－18 u??如果网格粗错不能用来求解层流底层，则假设与壁面近邻的网格质心落在边界层的对数区，则根据壁面法则： uu????u?y??lnE???
3－19 k???1其中，k=0.419，E=9.793。
对流传热传质模型
在FLUENT中，用雷诺相似湍流输运的概念来模拟热输运过程。给出的能量方程为： ??t(?E)???xi[ui(?E?p)]?cp?t????k????xi?Prt????T???x?uj(?ij)eff?i???Sh
3－20 ??式中，E是总能量，(?ij)eff是偏应力张量，定义为： (?ij)eff??eff(?uj?xi??ui?xj)?23?eff?ui?xi?ij
3－21 其中，(?ij)eff表示粘性加热，耦合求解。如果默认为分开求解，FLUENT不求解处(?ij)eff。但是可以通过变化“粘性模型”面板上的湍流普朗特数（Prt），其默认值为0.85。 湍流质量输运与热输运类似，默认的Schmidt数是0.7，该值同样也可以在“粘性模型”面板上调节。 标量的壁面处理与动量壁面处理类似，分别选用合适的壁面法则。 综上所述，Spalart-Allmaras模型是相对简单的单方程模型，只需求解湍流粘性的输运方程，并不需要求解当地剪切层厚度的长度尺度。该模型对于求解有壁面影响流动及有逆压力梯度的边界层问题有很好模拟效果，在透平机械湍流模拟方面也有较好结果。 Spalart-Allmaras模型的初始形式属于对低雷诺数湍流模型，这必须很好解决边界层的粘性影响区求解问题。在FLUENT中，当网格不是很细时，采用壁面函数来解决这一问题。当网格比较粗糙时，网格不满足精确的湍流计算要求，用壁面函数也许是最好的解决方案。另外，该模型中的输运变量在近壁处的梯度要比k??中的小，这使得该模型对网格粗糙带来数值误差不太敏感。 但是，Spalart-Allmaras模型不能预测均匀各向同性湍流的耗散。并且，单方程模型没有考虑长度尺度的变化，这对一些流动尺度变换比较大的流动问题不太适合。比如，平板射流问题，从有壁面影响流动突然变化到自由剪切流，流场尺度变化明显。
xxviii 3.3.2 标准k??模型
标准k??模型需要求解湍动能及其耗散率方程。湍动能输运方程是通过精确的方程推导得到，但耗散率方程是通过物理推理，数学上模拟相似原形方程得到的。该模型假设流动为完全湍流，分子粘性的影响可以忽略。因此，标准k??模型只适合完全湍流的流动过程模拟。 标准k??模型的湍动能k和耗散率ε方程为如下形式： ?DkDtD???t????????xi???k???k??
3－22 ??x??Gk?Gb????YM?i?2???????
3－23 ??x??C1?k(Gk?C3?Gb)?C2??k?i??t?????????Dt?xi???k?在上述方程中，Gk表示由于平均速度梯度引起的湍动能产生，Gb是用于浮力影响引起的湍动能产生；YM可压速湍流脉动膨胀对总的耗散率的影响。湍流粘性系数k2?t??C??在FLUENT中，作为默认值常数，C1?＝1.44，C2?=1.92，C??0.09，湍动能k与。 耗散率ε的湍流普朗特数分别为?k＝1.0，??＝1.3。可以通过调节“粘性模型”面板来调节这些常数值。
3.3.3 重整化群κ-ε模型
重整化群κ-ε模型是对瞬时的Navier-Stokes方程用重整化群的数学方法推导出来的模型。模型中的常数与标准κ-ε模型不同，而且方程中也出现了新的函数或者项。其湍动能与耗散率方程与标准κ-ε模型有相似的形式： ???k??????k?eff???Gk?Gb????YM
3－24 Dt?xi??xi?Dk2??????????R
?????eff???C1?(Gk?C3?Gb)?C2??Dt?xi??xi?kkGk表示由于平均速度梯度引起的湍动能产生，Gb是用于浮力影响引起的湍动能产生；YMD?可压速湍流脉动膨胀对总的耗散率的影响，这些参数与标准κ-ε模型中相同。?k和??分别是湍动能k和耗散率?的有效湍流普朗特数的倒数。
湍流粘性系数计算公式为： ~??2k??~??1.72d?d?
????3~??1?C???~其中，???eff/?，C??100 对上面方程积分，可以精确得到有效雷诺数（涡旋尺度）对湍流输运的影响，这有助于处理低雷诺数和近壁流动问题的模拟。 对于高雷诺数，上面方程可以给出：?t??C?k2?，C??0.0845。这个结果非常有意思，和标准κ-ε模型的半经验推导给出的常数C??0.09非常近似。 在FLUENT中，如果是默认设置，用重整化群κ-ε模型时候是针对的高雷诺数流动 xxix 3.3.2 标准k??模型
标准k??模型需要求解湍动能及其耗散率方程。湍动能输运方程是通过精确的方程推导得到，但耗散率方程是通过物理推理，数学上模拟相似原形方程得到的。该模型假设流动为完全湍流，分子粘性的影响可以忽略。因此，标准k??模型只适合完全湍流的流动过程模拟。 标准k??模型的湍动能k和耗散率ε方程为如下形式： ?DkDtD???t????????xi???k???k??
3－22 ??x??Gk?Gb????YM?i?2???????
3－23 ??x??C1?k(Gk?C3?Gb)?C2??k?i??t?????????Dt?xi???k?在上述方程中，Gk表示由于平均速度梯度引起的湍动能产生，Gb是用于浮力影响引起的湍动能产生；YM可压速湍流脉动膨胀对总的耗散率的影响。湍流粘性系数k2?t??C??在FLUENT中，作为默认值常数，C1?＝1.44，C2?=1.92，C??0.09，湍动能k与。 耗散率ε的湍流普朗特数分别为?k＝1.0，??＝1.3。可以通过调节“粘性模型”面板来调节这些常数值。
3.3.3 重整化群κ-ε模型
重整化群κ-ε模型是对瞬时的Navier-Stokes方程用重整化群的数学方法推导出来的模型。模型中的常数与标准κ-ε模型不同，而且方程中也出现了新的函数或者项。其湍动能与耗散率方程与标准κ-ε模型有相似的形式： ???k??????k?eff???Gk?Gb????YM
3－24 Dt?xi??xi?Dk2??????????R
?????eff???C1?(Gk?C3?Gb)?C2??Dt?xi??xi?kkGk表示由于平均速度梯度引起的湍动能产生，Gb是用于浮力影响引起的湍动能产生；YMD?可压速湍流脉动膨胀对总的耗散率的影响，这些参数与标准κ-ε模型中相同。?k和??分别是湍动能k和耗散率?的有效湍流普朗特数的倒数。
湍流粘性系数计算公式为： ~??2k??~??1.72d?d?
????3~??1?C???~其中，???eff/?，C??100 对上面方程积分，可以精确得到有效雷诺数（涡旋尺度）对湍流输运的影响，这有助于处理低雷诺数和近壁流动问题的模拟。 对于高雷诺数，上面方程可以给出：?t??C?k2?，C??0.0845。这个结果非常有意思，和标准κ-ε模型的半经验推导给出的常数C??0.09非常近似。 在FLUENT中，如果是默认设置，用重整化群κ-ε模型时候是针对的高雷诺数流动 xxix 问题。如果对低雷诺数问题进行数值模拟，必须进行相应的设置。
重整化群κ-ε模型有旋修正
通常，平均运动有旋时候对湍流有重要影响。FLUENT中重整化群κ-ε模型通过修正湍流粘性系数来考虑了这类影响。 湍流粘性的修正形式为： ?t??t0f(?s,?,k?)
3－27 其中，?t0是不考虑有旋计算出来的湍流粘性系数；Ω是FLUENT计算出来的特征旋流数；?s是旋流常数，不同值表示有旋流动的强度不同。流动可以是强旋或者中等旋度的。FLUENT默认设置?s＝0.05,针对中等旋度的流动问题，对于强旋流动，可以选择较大的值。 湍动能及其耗散率的有效湍流普朗特数倒数的计算公式为： ??1..??2..??mol?eff
3－28 式中，?0＝1，在高雷诺数流动问题中，?mol/?eff??1，?k????1.393。 湍流耗散率方程右边的R为： R?C???(1??/?0)?1???332k
3－29 其中，??Sk/?，?0?4.38，??0.012。 为了更清楚体现R对耗散率的影响，我们把耗散率输运方程重写为： 22????????*???C2?? 3－30 ?????eff???C1?(Gk?C3?Gb)?C2??Dt?xi??xi?kkkD?*则：C2??C2??C???(1??/?0)1???33
3－31 **在???0的区域，R的贡献为正；C2?大于C2?。以对数区为例，??3，C2??2.0，这和标准κ-ε模型中给出的C2?＝1.92接近。因此，对于弱旋和中等旋度的流动问题，重整化群κ-ε模型给出的结果比标准κ-ε模型的结果要大。 重整化群模型中，C1??1.42，C2??1.68。
3.3.4 可实现κ-ε模型
可实现κ-ε模型的湍动能及其耗散率输运方程为： ?t??????????Dt?xj??k??DkD???k???xj????Gk?Gb????YM
??2?t?????????????????CS???C?CC3?Gb
???121??Dt?xj??t??xj?kk??????其中，C1?max?0.43,????，??Sk/? ??5?? xxx 在上述方程中，Gk表示由于平均速度梯度引起的湍动能产生，Gb是用于浮力影响引起的湍动能产生；?k，??YM可压速湍流脉动膨胀对总的耗散率的影响。C2和C1?是常数；分别是湍动能及其耗散率的湍流普朗特数。在FLUENT中，作为默认值常数，C1?＝1.44，C2=1.9，?k＝1.0，??＝1.2。 可实现κ-ε模型的湍动能的输运方程与标准κ-ε模型和重整化群κ-ε模型有相同的形式，只是模型参数不同。但耗散率方程有较大不同。首先耗散率产生项（方程右边第二项）不包含湍动能产生项Gk，现在的形式更能体现能量在谱空间的传输。另外的特色在于耗散率减少项中，不具有奇异性。并不象标准κ-ε模型模型那样把K放在分母上。 该模型适合的流动类型比较广泛，包括有旋均匀剪切流，自由流（射流和混合层），腔道流动和边界层流动。对以上流动过程模拟结果都比标准κ-ε模型的结果好，特别是可实现κ-ε模型对圆口射流和平板射流模拟中，能给出较好的射流扩张角。 ?κ-ε模型中，C?不再是个常数，而是通过如下公式计算： C??1A0?AsUK*湍流粘性系数公式为?t??C?k2，这和标准κ-ε模型相同。不同的是，在可实现
3－34 其中，U*??~~~?ij??ij??ijk?k，?ij??ij?2?ijk?k，the mean ?ij是 is SijSij??ij?ij，6cos?，而： SijSij，Sij?ui1?uj?(?) 2?xi?xjrate-of –rotation tensor viewed in a rotating reference frame with the angular velocity ?k。模型常数A0?4.04，As???13arccos(6W)，式中W＝SijS~SjkSkj~，S?我们可以发现，C?是平均应变率与旋度的函数。在平衡边界层惯性底层，可以得到C?＝0.09,与标准κ-ε模型中采用底常数一样。
双方程模型中，无论是标准κ-ε模型、重整化群κ-ε模型还是可实现κ-ε模型，三个模型有类似的形式，即都有κ和ε的输运方程，它们的区别在于：1，计算湍流粘性的方法不同；2，控制湍流扩散的湍流Prandtl数不同；3，ε方程中的产生项和Gk关系不同。但都包含了相同的表示由于平均速度梯度引起的湍动能产生Gk，用于浮力影响引起的湍动能产生Gb；可压速湍流脉动膨胀对总的耗散率的影响YM。 湍动能产生项 Gk???ui?u?j?uj?xi
Gb??gi?t?TPrt?xi式中，Prt是能量的湍流普特朗数，对于可实现κ-ε模型，默认设置值为0.85。对于重整化群κ-ε模型，Prt?1/?，??1/Pr?k/?Cp。热膨胀系数???1??????，对???T?p于理想气体，浮力引起的湍动能产生项变为： ?t??Gb??gi
3－37 ?Prt?xi xxxi 在FLUENT程序中，如果有重力作用，并且流场里有密度或者温度的梯度，浮力对湍动能的影响都是存在的。浮力对耗散率的影响不是很清楚，因此，默认设置中，耗散率方程中的浮力影响不被考虑。如果要考虑浮力对耗散率的影响，用“粘性模型”面板来控制。浮力对耗散率影响是用C3?来体现。但C3?并不是常数，而是如下的函数形式： C3??tanhvu
v是平行于重力方向的速度分量；u是垂直于重力方向的速度分量。如果流动速度与重力方向相同的剪切流动，C3?＝1，对于流动方向与重力方向垂直的剪切流，C3?＝0。
对于高马赫数的流动问题，可压速性对湍流影响在YM中体现。 YM???2M2t k其中，Mt是马赫数，定义为：Mt?（a??RT是声速）。 2a默认设置中，只要选择可压速理想气体，可压速效应都是考虑的。 在上述的双方程模型中，对流传热传质模型都是通过雷诺相似湍流动量输运方程得到的。能量方程形式为： ??t(?E)???xi?ui??E?p???式中，E是总的能量，keff(?ij)eff??eff????T?keff??Sh
?u(?)jijeff???xi??xi?是有效导热系数；(?ij)eff是偏应力张量，定义为： ?2?ui????
3－40 ?3eff?xiji???uj?ui????x?xji?(?ij)eff表示的是粘性加热，耦合求解时总是计算。如果不是耦合求解时候，作为默认设置，并不求解该量。如果有需要，需在“粘性模型”面板中设置。 对于重整化群κ-ε模型，有效导热系数为： keff??cp?eff
3－41 α用（3－28）计算，式中，?0?1/Pr?k/?Cp。事实上，?随着?mol/?eff的变化而变化，这是重整化群κ-ε模型的一个优点，因为实验中证明，湍流普朗特数随分子普朗特数及湍流而变化。 湍流质量输运处理过程与能量输运过程类似。对于标准κ-ε模型和可实现的κ-ε模型，默认的Schmidt数是0.7，重整化群模型中，是通过方程3－28来计算的，其中，?0?1/Sc，Sc是分子Schimidt数。
3.3.5 雷诺应力模型（RSM）
雷诺应力模型是求解雷诺应力张量的各个分量的输运方程。具体形式为： ??t(?uiuj)???xk(?Ukuiuj)?
对流项Cij ???xk??uuiu?p(?kjui??ikuj)?jkT??????uuij? ??xk??xk?L湍流扩撒项Dij
分子扩散Dij
xxxii ?Uj??Ui?????giuj??gjui? ???uu?uujk?ik?x?xk?k??应力产生项Pij
浮力产生项目Gij ????ui?uj??ui?uj?? ?p??2???x??xi??xk?xkj?压力应变项?ij
?2??kujum?ikm?uium??jkm?
3－42 系统旋转产生项Fij 上面方程中，Cij，Dij，Pij，Fij不需要模拟，而Dij，Gij，?ij，?ij需要模拟以封闭方程。下面简单对几个需要模拟项的模拟。 Dij可以用Delay and Harlow [L38]的梯度扩散模型来模拟，即： DTijLTT??kukul?uiu?Cs???xk???xl??uiu????t??xk??k?xk?j??
3－43 ??但这个模型会导致数值不稳定，因此FLUENT程序中采用标量湍流扩散模型： DTijj??
???k??0.82，这和标准κ-ε模型中选取1.0有所不同。 式中，湍流粘性系数用?t??C?k2来计算，根据Lien and Leschziner [L98]，根据Gibson and Launder [L58], Fu [L55], Launder [L88,L89], 压力应变项?ij可以分解为三项，即： ?ij??ij,1??ij,2??ij
3－45 ?ij,1，?ij,2和?wij分别是慢速项，快速项和壁面反射项。 w?ij,1??C1???2?，常数C1?1.8。 uu??kijij?k?3??12???ij,2??C2??Pij?Fij?Gij?Cij???ij?P?G?C??，C2?0.60，P?Pkk，23??11G?Gkk，C?Ckk。 22壁面反射项用于重新分布近壁的雷诺正应力分布，主要是减少垂直于壁面的雷诺正应力，增加平行于壁面的雷诺正应力。该项模拟为： ?wij?C1????ukumnknm?ijk??k?uiuknjnk?ujuknink? 22C?d?l?knn
3－46 jk,2ik??Cl?d3/????km,2nknm?ij??ik,2njnk???C222???0.3，nk是xk在垂直于壁面方向上的单位分量，d是到壁面的式中，C1?＝0.5，C2距离；Cl?C?3/4/k，C??0.09，k=0.41。 wij默认设置时候，FLUENT不计算? 。如果需要计算时候，在“粘性模型”面板中设xxxiii ?Uj??Ui?????giuj??gjui? ???uu?uujk?ik?x?xk?k??应力产生项Pij
浮力产生项目Gij ????ui?uj??ui?uj?? ?p??2???x??xi??xk?xkj?压力应变项?ij
?2??kujum?ikm?uium??jkm?
3－42 系统旋转产生项Fij 上面方程中，Cij，Dij，Pij，Fij不需要模拟，而Dij，Gij，?ij，?ij需要模拟以封闭方程。下面简单对几个需要模拟项的模拟。 Dij可以用Delay and Harlow [L38]的梯度扩散模型来模拟，即： DTijLTT??kukul?uiu?Cs???xk???xl??uiu????t??xk??k?xk?j??
3－43 ??但这个模型会导致数值不稳定，因此FLUENT程序中采用标量湍流扩散模型： DTijj??
???k??0.82，这和标准κ-ε模型中选取1.0有所不同。 式中，湍流粘性系数用?t??C?k2来计算，根据Lien and Leschziner [L98]，根据Gibson and Launder [L58], Fu [L55], Launder [L88,L89], 压力应变项?ij可以分解为三项，即： ?ij??ij,1??ij,2??ij
3－45 ?ij,1，?ij,2和?wij分别是慢速项，快速项和壁面反射项。 w?ij,1??C1???2?，常数C1?1.8。 uu??kijij?k?3??12???ij,2??C2??Pij?Fij?Gij?Cij???ij?P?G?C??，C2?0.60，P?Pkk，23??11G?Gkk，C?Ckk。 22壁面反射项用于重新分布近壁的雷诺正应力分布，主要是减少垂直于壁面的雷诺正应力，增加平行于壁面的雷诺正应力。该项模拟为： ?wij?C1????ukumnknm?ijk??k?uiuknjnk?ujuknink? 22C?d?l?knn
3－46 jk,2ik??Cl?d3/????km,2nknm?ij??ik,2njnk???C222???0.3，nk是xk在垂直于壁面方向上的单位分量，d是到壁面的式中，C1?＝0.5，C2距离；Cl?C?3/4/k，C??0.09，k=0.41。 wij默认设置时候，FLUENT不计算? 。如果需要计算时候，在“粘性模型”面板中设xxxiii 置。
线性压力应变模型
对于小雷诺数流动，特别是用双层模型求解近壁流动问题时，FLUENT中通过改进模?来改进压力应变项Launder [L91]。这一过程只有在选择双层流型常数C1，C2，C1?和C2模型时候，在“粘性模型”面板上调节。 C1?1?2.58AC2?0.75A ?(0.?e?Ret)2?? ?C1???23C1?1.67 1??2?C2??6? ??max?3C2,0C2??????其中，Ret??k2/(??)，参数A和张量不变量A2, A3定义为： 9??A??1?(A2?A3)? 8??A2??ik?ki A3??ik?kj?ji 式中，?ij是雷诺应力张量各向异性部分，定义为： 2????uiuj??k?ij3????k??????
3－47 ????ij 二阶压力应变模型
二阶压力应变模型由Spezible {L157}等人提出。 ?ij??(C1???C1P)bij?C2??(bikbkj??C4?k(bikS?bjkSik?23*13bmnbmn?ij)?(C3?C3*bijbij)?kSij jkbmnSmn?ij)?C5?k(bik?jk?bjk?ik)
3－48 式中，bij是雷诺应力各向异性张量，定义为： 2???uu??k?ij?ij3
bij???2?k??????
?????ui1??uj?；????ij?2??xj???xi??；模型常数??平均应变率Sij定义为: Sij*?ui1??uj???2??xj??xi*C1?3.4，C1?1.8，C2?4.2，C3?0.8，C3?1.3，C4?1.25，C5?0.4。 二阶压力应变模型不需要考虑壁面反射影响去模拟对数区湍流边界层过程。
浮力对湍流的影响
浮力引起的产生项模拟为： Gij?T?T?g???gijPrt??x?xij??t???
3－50 ??其中，Prt是能量的湍流普朗特数，默认设置值为0.85。 对于理想气体，把热膨胀系数的定义代入上式，得： Gij??????g???gij?Prt??x?xij??t??
?? 耗散项?ij的模拟 耗散张量?ij模拟为： ?ij?23?ij(???YM)
3－52 2式中，YM???2Mt，Mt是马赫数；标量耗散率?用标准k-?模型中的采用的耗散率输运方程求解。
雷诺应力模型的边界条件
在流场进口，雷诺应力模型需要各个雷诺应力分量和湍动能耗散率的值。这些值可以直接输入，也可以湍流强度和特征长度来计算。 在壁面，雷诺应力模型通过壁面函数，给出各个雷诺应力分量和耗散率的值。
雷诺应力模型的能量与质量输运方程
在雷诺应力模型中，对流传热传质模型都是通过雷诺相似湍流动量输运方程得到的。能量方程形式为： ??t(?E)???xi?ui??E?p???cp?t?T???(k?)?uj(?ij)eff??xi?Prt?xi???Sh
3－53 ??式中，E是总的能量；(?ij)eff是偏应力张量，定义为： ??uj?ui????x?xji??2?ui????
3－54 ?3eff?xiji?(?ij)eff??eff(?ij)eff表示的是粘性加热，耦合求解时总是计算。如果不是耦合求解时候，作为默认设置，并不求解该量，并且Prt=0.85。如果有需要，需在“粘性模型”面板中设置。
3.3.6 大涡模拟（LES）
湍流中包含了不同时间与长度尺度的涡旋。最大长度尺度通常为平均流动的特征长度尺度。最小尺度为Komogrov尺度。 LES的基本假设：1，动量、能量、质量及其它标量主要由大涡输运；2，流动的几何和边界条件决定了大涡的特性，而流动特性主要在大涡中体现；3，小尺度涡旋受几何和边界条件影响较小，并且各向同性；大涡模拟过程中，直接求解大涡，小尺度涡旋模拟，从而使得网格要求比DNS低。
3.3.6.1大涡模拟的控制方程
LES的控制方程是对Navier-Stokes方程在波数空间或者物理空间进行过滤得到的。过滤的过程是去掉比过滤宽度或者给定物理宽度小的涡旋，从而得到大涡旋的控制方程。 过滤变量（上横线）定义为： ?(x)??D?(x?)G(x,x?)dx?
其中，D表示流体区域；G是决定涡旋大小的过滤函数。 在FLUENT中，有限控制体离散本身暗中包括了过滤运算， ?(x)?1V?V?(x?)dx?， x??V
3－56 其中V是计算控制体体积，过滤函数为： x??V?1/VG(x,x?)??
3－57 x??V?0目前，大涡模拟对不可压流动问题得到较多应用，但在可压缩问题中的应用还很少，因此这里涉及的理论都是针对不可压流动的大涡模拟方法。在FLUENT中，大涡模拟只能针对不可压流体（当然并非说是密度是常数）的流动。 过滤不可压的Navier-Stokes方程后，可以得到LES控制方程： ???t???ui?xi?0
??t(?ui)???xj(?uiuj)???xj(??ui?xj)??p?xi???ij?xj
其中，?ij为亚网格应力，定义为： ?ij??uiuj??ui?uj
3－60 很明显，上述方程与雷诺平均方程很相似，只不过大涡模拟中的变量是过滤过的量，而非时间平均量，并且湍流应力也不同。
3.3.6.2 亚网格模型
由于LES中亚网格应力项是未知的，并且需要模拟以封闭方程。目前，采用比较多的亚网格模型为涡旋粘性模型，形式为： ?ij?13?kk?ij??2?tSij
3－61 式中，?t是亚网格湍流粘性系数；Sij是求解尺度下的应变率张量，定义为： Sij?uj?ui1????2??xj?xi???
3－62 ??求解亚网格湍流粘性系数?t时，FLUENT提供了两种方法。第一，Smagorinsky-Lilly模型；第二，基于重整化群的亚网格模型。 最基本的亚网格模型是Smagorinsky [L145]最早提出的，Lilly [L99]把它进行了改善，这就是今天的Smagorinsky-Lilly模型。该模型的涡粘性计算方程为： ?t??LsS
2式中，Ls是亚网格的混合长度；S?度Ls可以用下式计算。 Ls?min(kd,CsV1/32SijSij。Cs是Smagorinsky常数，则亚网格混合长)
3－64 其中，k=0.42，d是到最近壁面的距离，V是计算控制体体积。 Lilly通过对均匀各向同性湍流惯性子区湍流分析，得到了Cs＝0.23。但是研究中发现，对于有平均剪切或者过渡流动中，该系数过高估计了大尺度涡旋的阻尼作用。因此，对于比较多的流动问题，Cs＝0.1有比较好的模拟结果，该值是FLUENT的默认设置值。 我们再来看看基于重整化群思想的亚网格模型。人们用重整化群理论推导出了亚网格涡旋粘性系数[L182]，该方法得到的是亚网格有效粘性系数，?eff????t，而 ?eff???s2?eff??????1?H?C????3???????1/31/3
3－65 式中，?s?(CrngV)22SijSij，H(x)是Heaviside函数， x?0?xH(x)??
3－66 0x?0?V是计算控制体体积；重整化群常数Crng?0.157,而常数C=100。 对于高雷诺数流动（?t??? ），?eff??t，基于重整化群理论的亚网格模型就与Smagorinsky-Lilly模型相同，只是模型常数有区别。在流动场的低雷诺数区域，上面的函数就小于零，从而只有分子粘性起作用。所以，基于重整化群理论的亚网格模型对流动转捩和近壁流动问题有较好模拟效果。
3.3.6.3 大涡模拟的边界条件
对于给定进口速度边界条件，速度等于各个方向分量与随机脉动量的和，即： ui??ui??I?u
3.3.6.2 亚网格模型
由于LES中亚网格应力项是未知的，并且需要模拟以封闭方程。目前，采用比较多的亚网格模型为涡旋粘性模型，形式为： ?ij?13?kk?ij??2?tSij
3－61 式中，?t是亚网格湍流粘性系数；Sij是求解尺度下的应变率张量，定义为： Sij?uj?ui1????2??xj?xi???
3－62 ??求解亚网格湍流粘性系数?t时，FLUENT提供了两种方法。第一，Smagorinsky-Lilly模型；第二，基于重整化群的亚网格模型。 最基本的亚网格模型是Smagorinsky [L145]最早提出的，Lilly [L99]把它进行了改善，这就是今天的Smagorinsky-Lilly模型。该模型的涡粘性计算方程为： ?t??LsS
2式中，Ls是亚网格的混合长度；S?度Ls可以用下式计算。 Ls?min(kd,CsV1/32SijSij。Cs是Smagorinsky常数，则亚网格混合长)
3－64 其中，k=0.42，d是到最近壁面的距离，V是计算控制体体积。 Lilly通过对均匀各向同性湍流惯性子区湍流分析，得到了Cs＝0.23。但是研究中发现，对于有平均剪切或者过渡流动中，该系数过高估计了大尺度涡旋的阻尼作用。因此，对于比较多的流动问题，Cs＝0.1有比较好的模拟结果，该值是FLUENT的默认设置值。 我们再来看看基于重整化群思想的亚网格模型。人们用重整化群理论推导出了亚网格涡旋粘性系数[L182]，该方法得到的是亚网格有效粘性系数，?eff????t，而 ?eff???s2?eff??????1?H?C????3???????1/31/3
3－65 式中，?s?(CrngV)22SijSij，H(x)是Heaviside函数， x?0?xH(x)??
3－66 0x?0?V是计算控制体体积；重整化群常数Crng?0.157,而常数C=100。 对于高雷诺数流动（?t??? ），?eff??t，基于重整化群理论的亚网格模型就与Smagorinsky-Lilly模型相同，只是模型常数有区别。在流动场的低雷诺数区域，上面的函数就小于零，从而只有分子粘性起作用。所以，基于重整化群理论的亚网格模型对流动转捩和近壁流动问题有较好模拟效果。
3.3.6.3 大涡模拟的边界条件
对于给定进口速度边界条件，速度等于各个方向分量与随机脉动量的和，即： ui??ui??I?u
xxxvii 其中，I是脉动强度，?是高斯随机数，满足??0，???1。 如果网格足够密并可以求解层流底层的流动的话，壁面切应力采用线性应力应变关系，即： uu?uu???u?y?1k
3－67 如果网格不够细，则假定与壁面邻近网格质心落在边界层对数区内，则： ?lnE(?u?y?)
其中，k=0.418，E=9.793。
雷诺平均模型的比较
模型名字 Spalart-Allmaras 标准k?? RNG k?? Realizable k?? 雷诺应力模型
优点 计算量小，对一定复杂程度的 边界层问题有较好效果 应用多，计算量合适，有较多数据积累和相当精度 能模拟射流撞击，分离流，二次流，旋流等中等复杂流动 和RNG模型差不多，还可以模拟圆口射流问题 考虑的物理机理更仔细，包括了湍流各向异性影响 缺点 计算结果没有被广泛测试，缺少子模型，如考虑燃烧或浮力问题 对于流向有曲率变化，较强压力梯度，有旋问题等复杂流动模拟效果欠缺 受到涡旋粘性各向同性假设限制 受到涡旋粘性各向同性假设限制 CPU时间长（2～3倍），动量和湍流量高度耦合。
第四节 湍流模型算例及其设置
湍流模型设置命令：Define-model-viscous
无粘，层流和湍流
湍流模型选项
近壁处理方法选择
附加湍流选项
???? ???????? ???????? ?????? ??xxxix
算例一：有换热的腔道流动问题 adiabatic wall cold air V = 50 fpm T = 0°F 1 ft insulation constant temperature wall T = 100°F 10 ft
步骤： 1， 检查是否湍流：Re?5980 P 1 ft Dh2， 选择低雷诺数湍流模型－RNG k??模型；壁面处理用非平衡壁面函数，考虑压力梯度影响 3， 网格划分：（1）四边形网格；（2）由于在靠近水平板处，垂直方向梯度较大，则近壁网格加密，并保证第一个控制体在对数区内；（3）变化流向网格间距，用于捕捉边界层发展；（4），根据计算结果，自动调节网格，用于进一步计算温度梯度。
计算结果：
Velocity contours BLs on upper & lower surfaces accelerate the core flow
Temperature contours Important that thermal BL was accurately resolved as well P
算例二：圆柱绕流
air V = 4 fps 1 ft 2 ft wall 5 ft 14.5 ft Compute drag coefficient of the cylinder
步骤： 1， 确定雷诺数，ReD?24600 2， 钝体绕流，后面有不稳定的涡旋脱落。采用RNG k??模型，壁面处理是双层区模型； 3， 网格处理：近壁网格加密，由于是双层区模型，需要网格划分到y??1
计算圆柱绕流的涡旋脱落过程
xli 百度搜索“就爱阅读”,专业资料,生活学习,尽在就爱阅读网92to.com,您的在线图书馆百度搜索“就爱阅读”,专业资料、生活学习,尽在就爱阅读网92to.com,您的在线图书馆!
欢迎转载：
相关推荐：}
当前位置：&>&&>& > 中科大Fluent讲稿(完整版,完美整理).
中科大Fluent讲稿(完整版,完美整理).
导读：就爱阅读网友为您分享以下“中科大Fluent讲稿(完整版,完美整理).”的资讯，希望对您有所帮助，感谢您对92to.com的支持!2.1.2 动量守恒方程 惯性坐标系下，i方向的动量守恒方程为：
??t(?ui)???xj(?uiuj)???p?xi???ij?cj??gi?Fi
???u?uji???????x?xi?j????2?ul?????ij ，?gi，Fi?3?xl???式中，p是静压；?ij是应力张量，定义为：?ij是重力体积力和其它体积力（如源于两相之间的作用），Fi还可以包括其它模型源项或者用户自定义源项。 对于二维轴对称几何条件，轴向和轴向的动量守恒方程分别为： ??t(?u)?1?r?x(r?uu)?1?r?r(r?vu)???p?x ????1????u2r?2?(??v)?? ??r?x???x3??1????u?v??r??r?r???r?x?????Fx
2－4 ???和 ??t(?v)?1?r?x(r?uv)?1?r?r(r?vv)???p?r ?1????v?u????? ?r??r?x???x?r?????1????v2r?2?(??v)?? ??r?r???x3????2?vr2?w?(??v)???Fr
2－5 3rr2?2w是旋流速度。
2.1.3 能量方程
FLUENT可以计算流体和（或者）固体区域之间的传热问题。如果是周期性换热流动，则流动边界要给定周期边界条件。如果计算计算模型包括两个流动区域，中间被固体或者墙壁隔开的换热问题，则要特别注意：1，两个流体都不能用流出边界条件（outflow）；2，两个区域的流动介质可以不同，但要分别定义流体性质（如果计算组分，只能给一个混合组分）。
ix 流体1 流体2
FLUENT求解的能量方程形式如下： ??t(?E)???xi(ui(?E?p))???xi(keff?T?xi??hj?Jj??uj(?ij)eff?Sh
2－6 j?式中，keff?kt?k，为有效导热系数(湍流导热系数根据湍流模型来定义)。Jj?是组分j?的扩散通量。方程右边前三项分别为导热项，组分扩散项和粘性耗散项。Sh是包括化学反应热和其它体积热源的源项。其中， 2pui
2－7 E?h???2对于理想气体，焓定义为：h??mj?hj?；对于不可压缩气体，焓定义为：j?h??mj?hj??j?pT?。mj?是组分j?的质量分数，组分j?的焓定义为：hj???cp,j?dT，其中TrefTref?298.15K。
2.1.4 PDF模型的能量方程
如果在非绝热PDF燃烧模型模式下，FLUENT求解的总焓方程为： ??t(?H)???xi(?uiH)??(kt?H)??ik?ui?xk?Sh
2－8 ?xicp?xi假定刘易斯数为1，方程右边第一项为组分扩散和导热项的合并项；第二项为粘性耗散，为非守恒形式。总焓H定义为: H??mj?Hj?j? 组分j?的总焓定义为： THj???cp,j?dT?hj?(TrefTref,j?0,j?)
2－9 其中hj?(Tref,j?)是组分j?基于参考温度Tref,j?的生成焓。 虽然能量的标准形式里包括了压力做功和动能项，但在采用segregated solver求解不可压问题时候都可以忽略掉。当然，如果想不忽略它们的作用，可以在define/models/energy x 0中设置。对于可压缩流动问题，在用coupled solvers求解时总是考虑压力做功和动能项。 粘性耗散项是考虑流体中的粘性剪切作用产生的热量。如果用segregated solver求解，默认设置并没有考虑。如果Brinkman数（Br??Ue2k?T性加热一定不能忽略。这时候一定要设置Viscous Heating选项。对于可压缩流动，一般Br>1，如果还用segregated solver求解，一定要考虑粘性加热。如果是coupled solver求解，粘性加热会自动考虑。 Fluent求解焓方程时，组分扩散项都已经包括。用segregated solver求解，如果想不考，?T是系统温度差）大于1时，粘虑该项，可以在组分模型面板（Species Model Panel）中关闭能量扩散项。如果采用了非绝热的PDF燃烧模型，方程中并不明确出现该项，应为导热和组分扩散项合并为一项了。当用coupled solver求解时，能量方程总会考虑该项。
2.1.5 化学反应源项
化学反应源项如下： Sh,reaction0?h0j?????j?M?j????cp,j?dT?Rj?
2－10 Tref,j???Tref其中，hj?是组分j?的生成焓；Rj?是组分j?生成的体积率。对于非绝热PDF燃烧模型，（2－9），生成热定义在总焓中，所以化学反应热不包含在源项中。
2.1.6 固体区域的能量方程
在固体区域，FLUENT采用的能量方程为如下形式： ??t?h???xi(ui?h)???xi(k?T?xi????
2－11 )?q????体积热源。方程左边第二项表示式中，?是密度；h是显焓；k是导热系数；T是温度；q由于固体旋转或者平移运动热传输。方程右边两相分别为固体导热和体积热源。
2.1.7 固体内部导热各向异性的影响
当用segregated solver求解时，FLUENT允许你指定材料的各向导热系数。固体导热各向异性方程形式如下： ??xi(kij?T?xi) 其中，kij是导热系数矩阵。
2.1.8 进口热扩散
进口的净能量输运包括对流和扩散两部分。指定进口温度就可以确定对流部分，但扩散项取决于计算出来的温度场梯度。因此我们不能给定扩散分量或者净能量输运。但在一些问题中，我们更希望能给定净能量输运，而不是给定进口温度。如果用segregated solver求解时，可以在dfine/models/energy中去掉进口能量扩散，从而达到给定净进口能量输运。但是 xi 我们用coupled solver时，不能去掉能量扩散部分。
第二节 计算传热过程中用户输入
如果用FLUENT计算有传热的问题时候，必须击活相关模型和提供热边界条件，并且给出材料物性。这一系列过程如下： 1， 击活能量面板。Define-Models-Energy
2， （对于segregated solver）如果模拟粘性流动过程，而且要考虑粘性加热，击活Viscous Heating；Define-Models-Viscous Heating
3， 定义热边界条件（包括流体进口，出口和壁面）Define-Boundary Conditions。在流动进口和出口要给定温度，但壁面可以有如下边界条件选择： （1） 指定热流量 （2） 指定温度 （3） 对流换热 （4） 外部辐射 （5） 对流换热＋辐射换热
4， 定义材料热物性。Define-Materials. 比热和导热系数都要给出，并且可以用温度函数的形式给出。
2.2.1 温度限制
为了计算的稳定性，FLUENT对计算出来的温度给了范围限制。给定温度限制，一方面是为了计算稳定的需要，同时，真实温度也有其相应的范围。由于给定材料物性不好，或者其它原因，计算出的中间超过了物理应该达到的温度。FLUENT中，给定的最高温度5000K，最小温度1K，如果计算过程中的温度超过这个范围，那么就在这最高温度或最低温度值处锁定。如果你觉得这个限制不合理，你可以自己调节。Solve-control-limits
2.2.2 传热问题求解过程
对于一些简单的传热过程FLUENT的默认设置可以成功进行模拟，但如果要加快你的问题的收敛速度或者提高计算过程的稳定性，下面的一些过程就比较重要了。
2.2.3 松弛因子确定
如果用segregated solver求解能量方程，在solve-controls-solution处定义松弛因子。如果你采用非绝热PDF模型，也必须和通常一样设置包括温度项在内地松弛因子。在求解温度和焓时候，FLUENT默认设置能量方程松弛因子为1。在一些问题里，能量场影响流动场（物性随温度变化，或者有浮力），这时候松弛因子要小些，比如在0.8到1之间。如果流动场和温度场不是耦合的（没有随温度变化的热物性或者浮力影响），松弛因子就可以采用1。 如果我们求解的是焓方程（非绝热PDF燃烧模型），温度需要设置松弛因子。焓的变化中不是所有的都用来计算温度的变化。这对于一些问题，你需要流动场焓变化快，而温度不能变化太快（影响流体热物性太快）的解决很有好处。
xii 2.2.4 组分扩散项
如果用segregated solver求解组分输运方程，如果考虑组分扩散，计算收敛会比较困难。为了提高收敛性，可以在define-models-species处取消对组分扩散的考虑。这时候组分扩散对能量的影响就被忽略了。如果我们选择coupled solver求解，那么组分扩散一定是存在的。
2.2.5耦合和非耦合流动场与温度场计算
如果流动和传热不是耦合的（没有温度变化的热物性或者浮力影响），那么我们可以先求解绝热流动场，然后加进能量方程。这时候可以暂时先关闭动量或者能量方程中的一个，先求解另外的一个。Solve-controls-solution. 如果流动和温度场是耦合的，你可以先求解流动方程，收敛后再击活能量方程，一起求解。需要注意的是，Coupled solver 总是同时求解流动与能量方程。
2.2.6 传热计算结果输出
FLUENT提供了几种传热结果的输出形式。可以以图形的形式输出，也可以用下列参数或者函数形式输出。可以输出的参数包括：静温、总温、静焓、相对总温、壁面温度（内或外表面）、总焓、总焓方差、熵、总能量、内能、表面热流量、表面换热系数、表面Nusselt数和表面stanton数。定义焓或者能量输出的参数跟求解的是可压或不可压问题有关。 可以通过Report-Fluxes给出控制体的每个边界或者通过壁面的总的换热量。必须注意，要检测一下是否能量平衡，这可以检查是否求解已经收敛。 也可以通过Report-surface Integrals-Enthalpy给出某个边界或壁面的总的换热量。焓流??率定义为Q???HV?dA。
2－12 Report-Surface integrals-surface Heat transfer Codf. (wall heat flux)可以给出某个表面的平?hdA均换热系数(h?) 。
2－13 A 第三节 浮力驱动的流动和自然对流
对于混合对流问题，我们定义一个浮力参数为??GrRe2???gh?v2，当?接近或者超过1时，浮力对流动有很强的作用。相反，如果?很小，浮力的作用将可以忽略。而纯粹自然对流问题中，浮力诱发流动强弱我们用另外一个参数来衡量，即Rayleigh数。定义为： Ra?g??TL?/??
2－14 3其中，???1????T8为热膨胀系数；??10k?cp热扩散系数。如果Ra?10，自然对流处于8层流状态，在10?Ra?10为层流到湍流的过渡区域。
2.3.1 Boussinesq模型
对于许多的自然对流问题，采用Boussinesq模型比定义密度是温度的函数有更好的收敛性。该模型在所有求解方程中，认为密度是常数。但是，在动量方程中的浮力项中，密度 xiii
FLUENT商用程序可以模拟许多的工程实际问题，包括可压缩、不可压缩流动，牛顿流体、非牛顿流体，单相、多相流动，有旋、无旋流动，惯性坐标系、非惯性坐标系下的流动，有化学反应、无化学反应的流动问题等。其生成无结构网格的程序把计算复杂几何条件下的流动及传热传质问题变的简单。同时，软件还提供了许多的湍流模型、壁面处理及燃烧、传热模型供针对特定问题选择。用户自定义函数也为改进和完善模型，处理个性化问题和给出更合理的边界条件提供了可能。 本讲义以FLUENT5说明书为主要参考资料，介绍了该软件的基本功能、基本物理模型、湍流模型、湍流模拟的近壁处理及边界条件，并且对燃烧过程的模拟和用户自定义函数做了描述。通过本课程学习，可以掌握和利用FLUENT程序在流体及传热传质等领域进行数值研究。 燃烧模型部分由董刚副教授编译，方海生同学编译了用户自定义函数。基本物理模型，湍流模型及湍流模拟近壁处理及边界条件由刘明侯副教授编译。由于时间非常仓促（一个暑假时间），只能用不完全的内容作为计算流体和传热传质课程的内容。还有些内容来不及加入讲义内，希望以后逐步完善。文字没有很好地校对，一定会由错误、疏漏或不妥的地方，请同学们校正。
第一章， 概述…………………………………………………………………..….…...(1)
第二章， 基本物理模型………………………………………………………..…...….(5) 第一节， 连续和动量方程…………………………………………………….…(5) 第二节， 计算传热过程用户输入………………………………………….……(9) 第三节， 浮力驱动的流动和自然对流…………………………………………(10) 第四节， 有旋和旋转流动问题…………………………………………….……(15) 第五节， 可压流动…………………………………………………..……..……(18) 第六节， 无粘流动……………………………………………………..…..……(19) 第七节， 用户自定义标量输运模型…………………………………..….….…(22)
第三章， 湍流模型……………………………………………………….……..…….(24) 第一节， 前言……………………………………………………………………(24) 第二节， 平均量输运方程…………………………………………….. ……….(25) 第三节， 湍流模型……………………………………………………………….(26) 第四节， 湍流模型算例及其设置……………………………………………….(39) 第四章， 湍流流动地近壁处理…………………………………………………...….(44) 第五章， 边界条件…………………………………………………………………….(52) 第六章， FLUENT中地燃烧模拟…………………………………………..…….….(62) 第一节， 燃烧模拟的重要性……………………………………………..….….(62) 第二节， FLUENT燃烧模拟方法概要…………………………………..…….(62) 第三节， 气相燃烧模型…………………………………………………......…..(63) 第四节， 污染物模型…………………………………………………….………(75) 第五节， FLUENT中燃烧模拟计算的步骤和原则……………………………(77) 第七章， 自定义函数………………………………………………………………….(79) 第一节， 概要…………………………………………………………………….(79) 第二节， 书写UDFs …………………………………………………………….(80) 第三节， 编译连接UDFs………………………………………………………..(98) 第四节， 在FLUENT模型中使用UDFs………………………………………(106) 第五节， UDFs实例…………………………………………………………….(112) 第六节， UDFs的应用………………………………………………………….(128) 附录A，udf. H中的宏解释…………………………………………………………….(146) 附录B，常用C库函数…………………………………………………………………(147)
FLUENT简介
FLUENT是用于计算复杂几何条件下流动和传热问题的程序。它提供的无结构网格生成程序，把计算相对复杂的几何结构问题变得容易和轻松。可以生成的网格包括二维的三角 ii 形和四边形网格；三维的四面体、六面体及混合网格。并且，可以根据计算结果调整网格。这种网格的自适应能力对于精确求解有较大梯度的流场如自由剪切流和边界层问题有很实际的作用。同时，网格自适应和调整只是在需要加密的流动区域里实施，而非整个流动场，因此可以节约计算时间。
程序的结构
FLUENT程序软件包应该包括以下几个部分： 1， FLUENT解法器 2， prePDF，用于模拟PDF燃烧过程 3， GAMBIT，网格生成 4， TGrid，额外的处理器，用于从现有的边界网格生成体网格。 5， Filters(Translators)，转换其它程序生成的网格，用于FLUENT计算。可以接口的程序包括：ANSYS, I-DEAS, NASTRAN,PATRAN等。
基本程序结构示意图。
PDF程序 几何形状或GAMBIT 设置几何形状 生成2D或3D网格 网格 其它软件包，如CAD，CAE等 prePDF
2D或3D网格 FLUENT
网格输入及调整 物理模型 边界条件 流体物性确定 计算 结果后处理 边界网格边界和（或）体网格网格
TGrid 2D三角网格 3D四面体网格 2D和3D混合网格
FLUENT程序的用途
1， 采用三角形、四边形、四面体、六面体及其混合网格计算二维和三位流动问题； 计算过程中，网格可以自适应。 2， 可压缩与不可压缩流动问题；
iii 3， 稳态和瞬态流动问题； 4， 无粘流，层流及湍流问题； 5， 牛顿流体及非牛顿流体； 6， 对流换热问题（包括自然对流和混合对流）； 7， 导热与对流换热耦合问题； 8， 辐射换热； 9， 惯性坐标系和非惯性坐标系下的流动问题模拟； 10，多运动坐标系下的流动问题； 11，化学组分混合与反应； 12，可以处理热量、质量、动量和化学组分的源项； 13，用Lagrangian 轨道模型模拟稀疏相（颗粒，水滴，气泡等）； 14，多孔介质流动； 15，一维风扇、热交换器性能计算； 16，两相流问题； 17，复杂表面形状下的自由面流动；
二维网格：
triangle quadrilateral
tetrahedron 三维网格：
hexahedron prism or wedge pyramid
图1－1，FLUENT的基本控制体形状
用FLUENT程序求解问题的步骤
1， 确定几何形状，生成计算网格（用GAMBIT，也可以读入其它指定程序生成的网格）； 2， 选择2D或3D来模拟计算； 3， 输入网格； 4， 检查网格； 5， 选择解法器； 6， 选择求解的方程：层流或湍流（或无粘流），化学组分或化学反应，传热模型等。确定其它需要的模型如：风扇、热交换器、多孔介质等模型。 7， 确定流体物性；
iv 8， 指定边界条件； 9， 条件计算控制参数； 10，流场初始化； 11，计算； 12，检查结果； 13，保存结果，后处理等。
关于FLUENT解法器的说明
1， FLUENT 2D，
二维单精度解法器 2， FLUENT 3D，
三维单精度解法器 3， FLUENT 2ddp， 二维双精度解法器 4， FLUENT 3ddp， 三维双精度解法器
FLUENT程序启动方法
1， WINDOEWS NT下，点击FLUENT5。 2， 在MS－DOS下，键入命令。
FLUENT 2D/3D/2ddp/3ddp。平行计算命令为：FLUENT 2D/3D/2ddp/3ddp
-t x 。x 是处理器编号。如，我们用3号处理器计算三维双精度问题，命令为： FLUENT 3ddp –t3
FLUENT命令的一般形式为： FLUENT [version] [-help] [options]
FLUENT求解方法的选择
1．非耦合求解 ( Segregated ) 2．耦合隐式求解 ( Coupled Implicit ) 3．耦合显式求解 ( Coupled Explicit ) 非耦合求解方法主要用于不可压缩或压缩性不强的流体流动。耦合求解则可以用在高速可压缩流动。FLUENT默认设置是非耦合求解，但对于高速可压流动，有强的体积力（浮力或离心力）的流动，求解问题时网格要比较密，建议采用耦合隐式求解方法，可以耦合求解能量和动量方程，能比较快地得到收敛解。缺点是需要的内存比较大（是非耦合求解迭代时间的1.5-2倍）。如果必须要耦合求解，但是你的机器内存不够，这时候可以考虑用耦合显式解法器求解问题。该解法器也耦合了动量，能量及组分方程，但内存却比隐式求解方法小。缺点是收敛时间比较长。 这里需要指出的是非耦合求解的一些模型在耦合求解解法器里并不都有。耦合解法器没有的模型包括：多相流模型，混合分数/PDF燃烧模型，预混燃烧模型，污染物生成模型，相变模型，Rosseland辐射模型，确定质量流率的周期性流动模型及周期性换热模型等。 隐式( Implicit )：对于给定变量，单元内的未知值用邻近单元的已知和未知值计算得出。因此，每一个未知值会在不止一个方程中出现，这些方程必须同时解来给出未知量。 显式( Explicit )：对于给定变量，每一个单元内的未知量用只包含已知量的关系式计算得到。因此未知量只在一个方程中出现，而且每一个单元内的未知量的方程只需解一次就可以给出未知量的值。
v 2.2.4 组分扩散项
如果用segregated solver求解组分输运方程，如果考虑组分扩散，计算收敛会比较困难。为了提高收敛性，可以在define-models-species处取消对组分扩散的考虑。这时候组分扩散对能量的影响就被忽略了。如果我们选择coupled solver求解，那么组分扩散一定是存在的。
2.2.5耦合和非耦合流动场与温度场计算
如果流动和传热不是耦合的（没有温度变化的热物性或者浮力影响），那么我们可以先求解绝热流动场，然后加进能量方程。这时候可以暂时先关闭动量或者能量方程中的一个，先求解另外的一个。Solve-controls-solution. 如果流动和温度场是耦合的，你可以先求解流动方程，收敛后再击活能量方程，一起求解。需要注意的是，Coupled solver 总是同时求解流动与能量方程。
2.2.6 传热计算结果输出
FLUENT提供了几种传热结果的输出形式。可以以图形的形式输出，也可以用下列参数或者函数形式输出。可以输出的参数包括：静温、总温、静焓、相对总温、壁面温度（内或外表面）、总焓、总焓方差、熵、总能量、内能、表面热流量、表面换热系数、表面Nusselt数和表面stanton数。定义焓或者能量输出的参数跟求解的是可压或不可压问题有关。 可以通过Report-Fluxes给出控制体的每个边界或者通过壁面的总的换热量。必须注意，要检测一下是否能量平衡，这可以检查是否求解已经收敛。 也可以通过Report-surface Integrals-Enthalpy给出某个边界或壁面的总的换热量。焓流??率定义为Q???HV?dA。
2－12 Report-Surface integrals-surface Heat transfer Codf. (wall heat flux)可以给出某个表面的平?hdA均换热系数(h?) 。
2－13 A 第三节 浮力驱动的流动和自然对流
对于混合对流问题，我们定义一个浮力参数为??GrRe2???gh?v2，当?接近或者超过1时，浮力对流动有很强的作用。相反，如果?很小，浮力的作用将可以忽略。而纯粹自然对流问题中，浮力诱发流动强弱我们用另外一个参数来衡量，即Rayleigh数。定义为： Ra?g??TL?/??
2－14 3其中，???1????T8为热膨胀系数；??10k?cp热扩散系数。如果Ra?10，自然对流处于8层流状态，在10?Ra?10为层流到湍流的过渡区域。
2.3.1 Boussinesq模型
对于许多的自然对流问题，采用Boussinesq模型比定义密度是温度的函数有更好的收敛性。该模型在所有求解方程中，认为密度是常数。但是，在动量方程中的浮力项中，密度 xiii 才随温度变化。(???0)g???0?(T?T0)g，因而用???0(1???T)计算浮力项。这样的近似对密度变化很小的流动问题有较好计算结果。 该模型对封闭区域里的自然对流问题适合，如果模拟温度变化很小的流动场也同样适用。但是，如果计算组分，燃烧或者有化学反应的问题时，该方法不适合。
2.3.2浮力驱动流计算的用户输入
在计算混合对流和自然对流问题时，需要输入如下参数： 1， 求解能量方程；define-models-Energy
2， 击活重量选项面板define-operating conditions；在笛卡儿坐标系的每个方向上设置重力加速度。默认的值是零。 3， 如果采用不可压理想气体法则，注意检查define-operating conditions 中的压力设置到一个合适值（非零）。 4， 取决于你是否采用了Boussinesq模型，还有以下参数需要合适确定： （1） 如果不用Boussinesq模型，则要给定合适的密度；或者在define-materials 里给出密度与温度的函数关系。 （2） 如果采用了Boussinesq模型，给出参考温度（方程中的T0），在define-materials 里给定热膨胀系数。 （3） 如果计算问题包含多个流体物质，可以为单个物质考虑是否选用Boussinesq近似。 ??0gx?ps，条5， 在压力进口与出口边界条件下，应该输入的压力为等效压力p?s件是进口和出口没有外部压力梯度。Define-boundary conditions
6， solve-controls-solution，压力离散化方法选择。如果你采用的是四边形或六面体网格，并且Segregated solver 求解，建议选择PRESTO离散方法。
2.3.3 流体参考密度
如果不用Boussinesq近似假设，流体在动量方程中体积力项中出现参考密度?0，其形?ps??0gx，根据流体静力学平衡，我式为：(???0)g，根据FLUENT定义的压力p?s们有： ?ps?x??g
2－15 上式可以变成：
?p?s?x?(???0)g
2－16 所以，定义流体参考密度是十分重要的。 FLUENT默认设置中，计算所有网格的平均密度为参考密度。在一些计算问题里，也许你给出一个参考密度要比自动计算一个参考密度要好。例如，你计算的是自然对流问题，ps，而且有个压力边界条件，了解上面方程中p?s的定义是十分重要的。虽然你知道确切的但你还需要知道参考密度?0，而根据ps的值来计算p?s。因此，你必须明确给出参考密度，而不能等计算机自动计算。而有时候给定参考密度虽然对最后结果没有太多影响，但对收敛速度有明显作用。这时候通常用体平均密度作为参考密度。这里需要指出的是，如果你采用的是Boussinesq近似方法，不需要用参考密度，因此不需要特别给出。
对于高Rayleigh 数流动问题的求解，需要注意以下一些问题。需要指出的是对一些层流问题而言，高Rayleigh 数流动有可能没有稳态解存在。 当求解高Rayleigh 数（Ra?108）流动问题时，根据下列步骤将能得到最好结果。 第一步是求稳态近似结果： ? 选用First-order scheme，在小Rayleigh数下求得稳态解。（可以通过变化重力加速度的方法减少Ra数（比如从9.8降低到0.098，Ra数就降低了两个数量级） ? ? 用小Ra数的收敛解为初始值，求解高Ra数下的解。 得到收敛解后，可以换higher-order scheme 继续求解。 第二步是求与时间相关的稳定解： ? 用前面的稳态解为初始条件，在相同或略小Ra数下求解。 ? 估计时间常数 ??LU~L2?(PrRa)?1/2?Lg??TL。其中，L和U是长度和速度? 尺度，采用的时间步长为：?t??/4，如果时间步长比?t大，有可能不收敛。 求解过程中会有频率为f??0.05~0.09振荡，衰减后就达到稳态解。?是上面的求出的时间常数，f是振荡频率（Hz）。通常需要超过5000步才能得到稳定解。
需要进一步指出的是除非我们采用了Boussinesq近似，上面方法不能用于封闭区域的流动问题，只能用于有进口和出口的流动问题。
2.3.4周期性流动与换热
如果我们计算的流动或者热场有周期性重复，或者几何边界条件周期性重复，就形成了周期性流动。FLUENT可以模拟两类周期性流动问题。第一，无压降的周期性平板问题（循环边界）；第二，有压降的周期性边界导致的完全发展或周期性流向流动问题（周期性边界）。 流向周期性流动模拟的条件： ? 流动是不可压的； ? 几何形状必须是周期性平移； ? 如果用coupled solver求解，则只能给定压力阶跃；如果是Segregated solver，可以给定质量流率或者压力阶跃。 ? ? ? 周期性流动中不能考虑进口和出口有质量差，也不考虑过程中的额外源项或者稀疏相源项。 只能计算进口出口没有质量流率变化的组分问题。但不能考虑化学反应。 不能计算稀疏相或者多相流动问题。 如果在这过程中计算有换热问题，则还必须满足以下条件： ? 必须用segregated solver 求解 ? 热边界条件必须是给定热流率或者给定壁面温度。对于一个具体的问题，热边界条件只能选择一个，而不能是多热边界条件问题。对于给定温度热边界条件，所有壁面的温度必须相同（不能有变化）。对于给定热流率边界条件，不同壁可以用不同值或曲线来模拟。 ? ? 对于有固体区域的问题，固体区域不能跨越周期性平板。 热力学和输运特性（热容，热导系数，粘性系数，密度等）不能是温度的函数（所以不能模拟有化学反应流动问题）。但输运特性（有效导热系数，有效粘性系数）可以随空间有周期性变化，因此可以对有周期性湍流输运特性不同的流动问题有模拟能力。
2.3.5 计算流向周期性流动问题的步骤：
通常，可以先计算周期性流动到收敛，这时候不考虑温度场。下一步，冻结速度场而计算温度场。步骤如下： 1， 建立周期性边界条件网格； 2， 输入热力学和分子输运特性参数； 3， 指定周期性压力梯度或者确定通过周期性边界的质量流量； 4， 计算周期性流动场。求解连续，动量（湍流量）方程。 5， 指定热边界条件（等温或者给定热流密度）； 6， 给定进口体平均温度； 7， 求解能量方程（其它方程不求解，只求解能量方程），得到周期性温度场。
2.3.5.1流向周期性流动理论 周期性速度定义 ?对于位置矢量r，周期性速度定义为： ?????
u(r)?u(r?L)?u(r?2L)?...
2－17 ?????
v(r)?v(r?L)?v(r?2L)?...
2－18 ?????w(r)?w(r?L)?w(r?2L)?...
2－19 ?其中，L是计算区域内周期性长度矢量。
2.3.5.2 周期性流向周期压力 上面方程中压力不是周期性的，而压力降是周期性的。即： ????????p?p(r)?p(r?L)?p(r?L)?p(r?2L)?...
2－20 如果选择coupled solver，压降?p是常数；但对Segregated solver方法求解，计算区域?L?内的压力梯度可以分解为两部分：梯度的周期性分量，?~p(r)，梯度的线性变化量，??。L即： ???L?p(r)?????~p(r)
L??周期性压力（~是但压力减去线性变化压力得到的值。压力的线性变化分量?r是p(r)）导致一个力作用于流体的动量方程。由于?并不能事先知道，只有通过在给定区域积分求出的流量与给定的的质量流量一致，才能确定。在压力修正的SIMPLE, SIMPLEC和PISO运算法则中都出现?，根据计算得到的质量流量与实际流量之差得到。
2.3.5.3 用Segregaed solver求解流向周期性流动的用户设置 要计算给定质量流量或压降的空间周期性流动问题，首先要建立计算网格。网格要求有平移周期性边界条件。在define-periodic conditions 面板上第一如下参数： 1， 决定是给定质量流量还是压力梯度。（Specify mass flow, Specify pressure Gradient） 2， 给出质量流量或压力梯度。如果你选择了给定质量流量，你也可以给出压力梯度（假定值），但不是必须的。如果给定值合理，可以加速流动场收敛速度。对计算结果也没有影响。对于轴对称问题，质量流量是2?弧度的质量流量。 3， 通过定义流向的X，Y，Z点来（或者二维时候的X，Y）定义流动方向。流动方向 xvi 是从初始点到指定点的方向上，该方向必须和周期性平移边界平行。 4， 如果确定了质量流量，FLUENT需要计算出压力梯度?。可以通过改变Relaxation Factor 和Number of Iterations（压力修正方程迭代次数）或者给出估计的?用以控制计算过程，这些都在define-periodicty conditions 面板上做。
2.3.5.4用Coupled solver方法求解周期性流动场 网格处理和前面类似，参数设置方面，需要： 1， define-boundary condition；
2， 给定周期性压降；
2.3.6对周期性流动问题计算过程的监视 如果是给定质量流量问题，可以在计算过程中观察压力降?（Statistic Monitors panel ）
2.3.7 给定温度边界条件下的周期性传热计算
要计算周期性流动中的传热问题，必须选用Segregated solver 方法。对于给定壁面温度条件，现在我们假设壁面温度是常数，则温度为：[L119] ???T(r)?Twall,inletTbulk?Twall
2－22 ???T?V?dA,inlet进口体温定义为：Tbulk?A， 这里的??，积分在周期性边界的进口上（A）???V?dAA在计算区域的L长度上满足周期性变化条件。
2.3.8 壁面温度确定周期性传热设置
1， 求解能量方程（define-models-Energy ） 2， 给出壁面温度。所有固体边界的温度相同，所有边界（除了周期性边界）必须用该温度的壁面封闭。Define-boundary conditions
3， 定义固体区域和固体温度。固体区域内不能有热源。Define-boundary conditions
4， 设定流体物性（密度，热容，粘性系数，导热系数等），但不能随温度变化。Define-Materials
5， 给定周期性进口上游的进口体平均温度（不能和壁面温度相同）define-periodic conditions
有了以上设置后就可以迭代求解了。最好的方法是先求解流场，等流场收敛后，再冻结流场只求解温度场。再给定处温度场时，温度值在壁温与进口体平均温度之间。 可以通过监视进出口平均体温度之比??monitors-per/bulk-temp-ratio
Twall?TbulkTwall?Tbulk,inlet,exit。Define-statistic 第四节 有旋与旋转流动问题
xvii 是从初始点到指定点的方向上，该方向必须和周期性平移边界平行。 4， 如果确定了质量流量，FLUENT需要计算出压力梯度?。可以通过改变Relaxation Factor 和Number of Iterations（压力修正方程迭代次数）或者给出估计的?用以控制计算过程，这些都在define-periodicty conditions 面板上做。
2.3.5.4用Coupled solver方法求解周期性流动场 网格处理和前面类似，参数设置方面，需要： 1， define-boundary condition；
2， 给定周期性压降；
2.3.6对周期性流动问题计算过程的监视 如果是给定质量流量问题，可以在计算过程中观察压力降?（Statistic Monitors panel ）
2.3.7 给定温度边界条件下的周期性传热计算
要计算周期性流动中的传热问题，必须选用Segregated solver 方法。对于给定壁面温度条件，现在我们假设壁面温度是常数，则温度为：[L119] ???T(r)?Twall,inletTbulk?Twall
2－22 ???T?V?dA,inlet进口体温定义为：Tbulk?A， 这里的??，积分在周期性边界的进口上（A）???V?dAA在计算区域的L长度上满足周期性变化条件。
2.3.8 壁面温度确定周期性传热设置
1， 求解能量方程（define-models-Energy ） 2， 给出壁面温度。所有固体边界的温度相同，所有边界（除了周期性边界）必须用该温度的壁面封闭。Define-boundary conditions
3， 定义固体区域和固体温度。固体区域内不能有热源。Define-boundary conditions
4， 设定流体物性（密度，热容，粘性系数，导热系数等），但不能随温度变化。Define-Materials
5， 给定周期性进口上游的进口体平均温度（不能和壁面温度相同）define-periodic conditions
有了以上设置后就可以迭代求解了。最好的方法是先求解流场，等流场收敛后，再冻结流场只求解温度场。再给定处温度场时，温度值在壁温与进口体平均温度之间。 可以通过监视进出口平均体温度之比??monitors-per/bulk-temp-ratio
Twall?TbulkTwall?Tbulk,inlet,exit。Define-statistic 第四节 有旋与旋转流动问题
xvii 有旋流动在一些燃烧器里用以增加流体滞驻时间和驻涡，从而稳定燃烧。旋转流动在透平机等实际过程中也应用广泛。有旋与旋转流动问题可以分为一下5类：1，轴对称有旋或旋转流动；2，完全三维有旋或旋转流动；3，一个坐标系旋转的流动；4，多坐标系旋转流动；5，流动需要滑移网格。 1， 轴对称有旋或旋转流动 流动没有周向梯度（但可以有周向速度）。周向动量方程为： ??t(?w)?1?r?x(r?uw)?1?r?r(r?vw)?1???w? r??r?x??x??1??3??w??r?2?r?rr?r?vw??
2－23 ????r??
?其中，x 和r分别是轴向和径向坐标；u, v, w分别是轴向，径向和周向速度分量。 2， 三维有旋流动 如果几何形状发生变化或者周向有流动梯度，则需要用三维模拟。对于三位问题，和解二维问题类似，没有特别的输入或者求解步骤。但是，定义速度进口条件时候需要用柱坐标系，并且，需要在求解时渐渐增加旋转速度。 3， 有旋转坐标系的流动 如果有旋转边界（螺旋桨等）必须用旋转坐标系来求解该类问题。如果有多个旋转边界，还需要多个旋转坐标系来求解。
2.4.1有旋或旋转流动的物理概念 对于有旋流动，轴向动量（rw 或r2?）往往会导致自由涡流动。周向速度w 随半径增加而增加，在粘性力起主导作用的r ?0的区域，周向速度近似为零。龙卷风就是一个很好的例子。 对于自由涡流动，流体周向运动的离心力与径向压力梯度相等：
??rr壁面旋转驱动的流动中，w/r或?是常数
有旋流动的湍流模拟
许多流动明显具有旋流（龙卷风，旋转射流等），必须考虑选用FLUENT提供的比较高级的模型，如RNG k-?模型、可实现k-?模型或者雷诺应力模型。具体选择哪个模型，取决于流动的旋流强度（旋流数）。旋流数定义为轴向与周向动量比： ???rwv?dA?
2－24 S??R?uv?dAR是水力学半径。 ?p?w2 对于弱旋和中等旋度流动问题（S0.5），必须选用雷诺应力模型。只有雷诺应力模型才能模拟该流动中的强的各向异性影响。 有旋流动问题的第二关键问的是边界条件问题。因为场的模拟好坏，主要取决于采用的模型，而壁面也参与了涡旋（涡量）的产生（由于压力梯度产生的二次流或涡流），采用非平衡壁面函数可以得到比较好的模拟结果，因为该壁面法则的平均流动速度对压力梯度敏 xviii 感。
2.4.2有旋与旋转流动问题的网格设置
坐标系限制：对于轴对称问题，旋转轴必须是x 轴，网格必须在y=0线以上。 除了上面的注意事项，对于有旋和旋转问题，网格划分还必须保证问题求解有足够分别率。特别是旋转流动中，边界层很薄，FLUENT需要在靠近旋转边界的地方网格比较细。除此之外，对于有旋流动问题，周向速度梯度比较陡（近中心线区域为典型的自由涡流动），因此要保证有很好的求解效果，网格要求比较密。
2.4.3轴对称有旋流动设置 1， 求解周向动量。Define-models-solver-Axisymmetric swirl
2， 给出进口或者壁面的旋转或者有旋速度分量，r?；define-boundary conditions (对于旋转轴，选用axis boundary边界条件)。
求解有旋或者旋转问题时的困难在于动量方程之间的耦合。如果是旋转很强的流动，会导致比较大的径向压力梯度，并驱动流体在轴向和径向的流动；动量之间的强耦合作用会导致求解过程中的不稳定性。如果要得到一个好的收敛解，需要有一定技巧。 求解步骤为： 1，（segregated solver only）如果用的是四边形和六面体网格，选用PRESTO模型（solution controls panel ） 2，为了求解大压力梯度和轴向速度梯度，网格必须足够精细。 3，(segregated solver only）改变松弛因子。径向和轴向速度为0.3-0.5，周向速度为0.8-1。 4，（segregated solver only）一步一步求解： （1） 如果包含inflow/outflow 的问题，那么先求解无旋流动。即用Axisymmetric ,而不选用Axisymmetric swirl option。并且不设置任何有旋边界条件。计算的结果作为有旋流动的初始值。 （2） 击活Axisymmetric Swirl option，设置有旋和旋转边界条件 （3） 只求解周向速度动量方程。让旋转的边界条件扩散到整个流场。如果是求解的湍流场，这适合湍流方程也应该同时求解。 （4） 关闭（冻结）求解周向动量方程，再求解连续和其它动量方程。如果是求解的湍流场，湍流量方程也同时求解。 （5） 松弛因子给定合适的值，同时求解所有方程。
除了上面描述的以为，如果求解的是有换热问题，可以先求解绝热流动场。如果求解的是湍流问题，也可以先计算层流，然后假如湍流模型继续计算。该方法对segregated or coulpled 求解都适合。 5，如果可能的话，先用小的旋转速度或者进口旋流速度进行计算。然后再增加到要计算的值。 （1） 再进口旋转速度或者旋转边界条件上给小的旋转速度，例如给需要计算值的10％。 （2） 在上面给定的条件下求解（可以用上面一步一步的求解方法） （3） 存储上面的初步结果。 （4） 更改初始和边界条件，增加旋转速度（也许是第一次的2倍）。 （5） 用上面的计算结果做初始值，开始新一轮计算。
xix （6） 进一步增加旋流速度，重复上面过程4－5，直到需要计算的旋流速度值。
第五节 可压流动
当气流速度很大，或者流场压力变化很大时，流体就收到了压速性的影响。 马赫数定义为：M?u/c，这里，c 是气体的音速，为：c? ?RT，??cp/cv比热比。当马赫数小于1时，流动为亚音速流动；当马赫数远远小于1（如M<0.1）时，流体的可压速性及压力脉动对密度变化影响都可以忽略。当马赫数接近1时候（跨音速），可压速性影响就显得十分重要了。如果马赫数大于1，流体就变为超音速流动。FLUENT对于亚音速，跨音速以及超音速等可压流动都有模拟能力。
2.5.1可压速流动的基本概念
描述可压速流动我们都离不开两个参数，总压p0和总温T0。对于理想气体，它们与静压与静温之间的关系为： ?p0ps??12???1???1?M?
2－25 2???1?T0Ts??12M
2－26 2上述关系成立的条件是等熵过程。
2.5.2 可压速流动的基本方程
FLUENT提供的标准连续和动量方程就可以描述可压速流动问题，除了下面会介绍的可压速流动处理以外，不需要其它特别的物理模型。需要指出的是，求解可压速问题，一定要求解能量方程。如果用segregated solver 求解，一定要考虑粘性耗散项（粘性加热）。 对于可压速流动，理想气体方程为： ??(pop?p)/RTs
2－27 这里，pop是运行条件里确定的运行压力，p是当地静压，R是通用气体常数，根据输入的气体分子量计算；Ts用能量方程求得。
2.5.3 可压速流动模拟的参数输入：
1，设定运行压力pop。Define-operating conditions ； 2，求解能量方程； 3，（Segregated solver only ）如果模拟的是湍流流动问题，考虑粘性耗散。Define- models- viscous.；(耦合求解，不需要，因为耦合求解自动考虑粘性耗散) 4，材料面板设置；Define-materials ； a) 选择理想气体； b) 定义物性（比如，分子量，导热系数等）。 5，设定边界条件； （1）Flow inlets
(a )压力进口：进口总温，总压；对于超音速流动，静压；
（b）质量进口：进口质量流率和总温。
xx （2）Flow exits
压力出口：出口静压（如果出口是超音速，可以忽略）； 需要特别指出的是，输入的边界条件中，无论是静压还是总压，都必须是表压（与前面给定的pop的差）。进口温度给定必须是总温（滞止温度），不是静温。
2.5.4 可压速流动求解注意事项
求解可压速问题的困难在于流体速度，密度，压力和能量的高度耦合。这样的耦合会导致求解过程中的不稳定性。因此要得到收敛解必须采取一定的手段。另外，超音速流动的激波也会导致求解的不稳定性。 下面介绍较好求解步骤： 1，（Segregated solver only ）速度的松弛因子调低（0.2-0.3）； 2，（Segregated solver only ）压力松弛因子用0.1，采用SIMPLE算法。（可压速流动不能采用SIMPLEC或者PISO） 3，对压力温度设置合适的限制条件。Solution limits。特别注意的是给定的压力和温度初始值要合理，如果给定的限制条件得到的收敛结果不好，可以更改继续计算，直至取得满意结果。 4， 有时候可以考虑用无粘收敛结果做初始场有好的收敛效果。
第六节 无粘流动
无粘流动是忽略了流体粘性作用，特别在大雷诺数流动中，惯性力起主导作用。在高速气体动力学里有比较多的应用。在这类流动中，压力作用在固体上的力比粘性力大很多，我们可以做无粘分析，快速得到作用在物体上的主要力的大小。然后，我们也可以假如粘性（包括湍流粘性）来评估对物体阻力或升力的影响。 另外一方面，也许我们求解的问题力粘性力不能忽略，但我们也可以先做无粘分析，用无粘结果作为有粘计算的初始值，这一方法，特别对一些复杂的流动，往往有好的收敛效果。
2.6.1 方程描述
无粘流动，我们求解的是Euler方程。连续方程与层流连续方程相同，而动量与能量方程由于没有考虑分子扩散而得到了简化。
连续方程：
???t???xi(?ui)?Sm
该方程是通用形式，对可压与不可压问题都适合。源项是稀疏相对连续相的贡献，也可以是用户自定义源项。 对于二维轴对称问题，连续方程可以写成： ?????v?(?u)?(?v)??Sm
2－29 ?t?x?rr动量守恒方程：
xxi （2）Flow exits
压力出口：出口静压（如果出口是超音速，可以忽略）； 需要特别指出的是，输入的边界条件中，无论是静压还是总压，都必须是表压（与前面给定的pop的差）。进口温度给定必须是总温（滞止温度），不是静温。
2.5.4 可压速流动求解注意事项
求解可压速问题的困难在于流体速度，密度，压力和能量的高度耦合。这样的耦合会导致求解过程中的不稳定性。因此要得到收敛解必须采取一定的手段。另外，超音速流动的激波也会导致求解的不稳定性。 下面介绍较好求解步骤： 1，（Segregated solver only ）速度的松弛因子调低（0.2-0.3）； 2，（Segregated solver only ）压力松弛因子用0.1，采用SIMPLE算法。（可压速流动不能采用SIMPLEC或者PISO） 3，对压力温度设置合适的限制条件。Solution limits。特别注意的是给定的压力和温度初始值要合理，如果给定的限制条件得到的收敛结果不好，可以更改继续计算，直至取得满意结果。 4， 有时候可以考虑用无粘收敛结果做初始场有好的收敛效果。
第六节 无粘流动
无粘流动是忽略了流体粘性作用，特别在大雷诺数流动中，惯性力起主导作用。在高速气体动力学里有比较多的应用。在这类流动中，压力作用在固体上的力比粘性力大很多，我们可以做无粘分析，快速得到作用在物体上的主要力的大小。然后，我们也可以假如粘性（包括湍流粘性）来评估对物体阻力或升力的影响。 另外一方面，也许我们求解的问题力粘性力不能忽略，但我们也可以先做无粘分析，用无粘结果作为有粘计算的初始值，这一方法，特别对一些复杂的流动，往往有好的收敛效果。
2.6.1 方程描述
无粘流动，我们求解的是Euler方程。连续方程与层流连续方程相同，而动量与能量方程由于没有考虑分子扩散而得到了简化。
连续方程：
???t???xi(?ui)?Sm
该方程是通用形式，对可压与不可压问题都适合。源项是稀疏相对连续相的贡献，也可以是用户自定义源项。 对于二维轴对称问题，连续方程可以写成： ?????v?(?u)?(?v)??Sm
2－29 ?t?x?rr动量守恒方程：
xxi ??t(?ui)???xj(?uiuj)???p?xi??gi?Fi
2－30 对于二维轴对称流动，上面方程可以写成： ??t??t(?u)?(?v)?1?r?x1?r?x(r?uu)?(r?uv)?1?r?r1?r?r(r?vu)??(r?vv)???p?x?p?r?Fx
2－32 能量守恒方程： ??t(?E)???xi(ui(?E?p))????xi?hj?j?Jj??Sh
2.6.2 求解无粘流动设置：
1， 设置无粘模型；Define-models-viscous ； 2， 设置边界条件和流体物；define-boundary condition /Materials ； 3， 求解，检查结果。
由于无粘问题往往是伴随高速流动，需要给出小的动量方程松弛因子。
有时间相关项的问题数值模拟
FLUENT可以解决的非定常问题包括：1，涡旋脱落等其它周期性现象；2，压缩填充与抽空问题；3，瞬时热导问题；4，瞬时化学混合与反应等。 求解时间相关项对一些定常问题求解时的收敛稳定性有时候有帮助，比如，求解Ra数比（层流向湍流过渡区域）较大的自然对流问题。对于一些问题，计算时间相关方程，是可以得到定常解的。
时间项离散
在FLUENT中，与时间相关项必须在时间和空间上离散。在空间的离散过程与求解稳态问题一样。 变量?随时间变化的过程通常可以写成如下形式： ???t?F(?)
2－34 其中，F是包含空间离散的函数。 如果时间导数项用向后差分，则一阶精度形式为： ?n?1?F(?)
2－35 ?t而二阶精度的离散形式为： n?1??n3??4?n??n?12?t时间项确定后，下面就该确定F(?)了。
xxii 一阶精度的隐式格式：二阶精度的隐式格式：?n?1??n?t3?n?1?F(?nn?1)
2－37 ?F(?n?1?4???n?12?t隐式格式的好处是对时间步无条件稳定。 )
2－38 一阶精度的显式格式：二阶精度的显式格式：?n?1??n?t3?n?1?F(?)
2－39 nn?4???n?12?t显式格式的时间步选择选择比较关键。显式格式都用来捕捉移动波，该格式比较精确，而且节约计算时间。但是下列情况不能用显式格式：1，不能用于segregated or coupled implicit solver.而只能用于coupled explicit solver .2，不能用于不可压速流动。不可压速流动必须在每个时间步长里积分收敛。
2－40 n 2.6.3 求解非定常问题的步骤：
1，击活求解非定常问题面板：define-models-solver。一阶精度的隐式格式对很多问题都适合。显式格式只能用于coupled explicit solver，用于捕捉移动波。 2，定义相关的模型和边界条件。自定义函数边界条件可以是跟时间相关。 3，如果用segregated solver, 则选用PISO来处理压力速度耦合。Solve-controls-solution。 4，监视计算结果。Solve-monitors-Statistic。 5，设置初始条件。Solve-initialize-initizlize。也可以读入稳态计算结果做初始值。File-read-data。
6，存贮中间结果。File-write-autosave。指定文件名字和存贮频率。 7，设定时间相关项处理方法 （1）如果采用一阶或者二阶隐式，需要设置如下参数： ? ? 每个时间步最多积分次数。 时间步长?t。虽然隐式格式对时间步长没有太多限制，但是对于我们模拟的一些瞬态问题，还应该把这个时间步定义为流动最小时间常数小一个数量级。最好的选择时间步的方法是看需要多少积分步才能收敛。理想的积分步是每个时间步长上为10－20次。如果需要迭代的次数很多，证明时间步太大了，需要调小点。如果迭代一两次就收敛了，证明时间步长可以调大点。对于周期下性的流动问题，时间步可以根据周期来确定。比如，一个周期可以分成20个积分步长。 （2）如果选择显式非稳定形式，则： ? ? 定义求解变量solve-controls-solution (必须没有多重网格，没有残差光滑，而且必须courant number of 1。 选择积分次数，求解。 8，存储结果，以便处理和进一步计算。File-write-data
第七节 用户自定义标量输运模型
Fluent可以在求解质量组分分数输运方程的同时，还可以求解用户自定义的任意标量输运方程。这样，我们可以求解磁流体动力学问题（MHD）。在这个过程中，流动的导电流体 xxiii 产生的磁场可以当做用户自定义标量来求解。流体产生的磁场又会对流体产生阻尼作用，这可以用用户自定义源项来处理。 对于任意的标量，?k，根据计算对流通量的不同方法，FLUENT提供了三个选择： 1，如果对流通量不需要求解，求解的方程为： ??k??(?k)?S ?
k=1, 2, ……..N
2－41 ?xi?xik 其中，?k 是扩散系数，S ?是第N个标量的源项。 k 2，如果对流通量需要计算，则求解的方程为： ??k?(?u i ?k ??k)?S ?，k=1, 2, ……..N
2－42 ?xi?xik 3，如果用户用自定义函数来确定对流通量，则方程为： ??k?(F i ?k ??k)?S ?，k=1, 2, ……..N
其中，F i 为自定义函数。 在FLUENT中，用户自定义标量只能在流体控制体中求解，而不能在固体控制体中求解。
用户自定义标量求解步骤： 1， 指定用户自定义标量数。Define-models-user-defined scalars
2， 确定对流通量确定形式。None (不需要计算)，mass flow rate (计算对流通量)，user defined（用户定义对流通量）。只能用一个确定对流通量的方法。 3， 定义自定义函数的边界条件，进口和出口边界条件。可以给固定值，也可以给出通量。 4， 如果需要用用户自定义源项，则：define-boundary conditions,给出自定义源项设置。 5， 给定初始值，求解 6， 结果检查。主要检查该标量的值scalar-n 和扩散系数diffusion coef.
Of scalar-n。
第三章 湍流模型 第一节 前言
湍流流动模型很多，但大致可以归纳为以下三类： 第一类是湍流输运系数模型，是Boussinesq于1877年针对二维流动提出的，将速度脉动的二阶关联量表示成平均速度梯度与湍流粘性系数的乘积。即： ?u1?u2???t??u1
3－1 ?x2推广到三维问题，若用笛卡儿张量表示，即有： ??ui?uj?????uiuj??t???x?xij??2???k?
3－2 ij?3?模型的任务就是给出计算湍流粘性系数?t的方法。根据建立模型所需要的微分方程的数目，可以分为零方程模型（代数方程模型），单方程模型和双方程模型。 第二类是抛弃了湍流输运系数的概念，直接建立湍流应力和其它二阶关联量的输运方程。
xxiv 第三类是大涡模拟。前两类是以湍流的统计结构为基础，对所有涡旋进行统计平均。大涡模拟把湍流分成大尺度湍流和小尺度湍流，通过求解三维经过修正的Navier-Stokes方程，得到大涡旋的运动特性，而对小涡旋运动还采用上述的模型。 实际求解中，选用什么模型要根据具体问题的特点来决定。选择的一般原则是精度要高，应用简单，节省计算时间，同时也具有通用性。 FLUENT提供的湍流模型包括：单方程（Spalart-Allmaras）模型、双方程模型（标准κ-ε模型、重整化群κ-ε模型、可实现(Realizable)κ-ε模型）及雷诺应力模型和大涡模拟。
包含更多 物理机理
Zero-Equation Models RANS-based
models One-Equation Models
Spalart-Allmaras Two-Equation Models
Standard k-? 每次迭代
RNG k-??计算量增加 ?????Realizable k-??
FLUENT提Reynolds-Stress Model 供的模型选 择 Large-Eddy Simulation
Direct Numerical Simulation
湍流模型种类示意图 xxv 第三类是大涡模拟。前两类是以湍流的统计结构为基础，对所有涡旋进行统计平均。大涡模拟把湍流分成大尺度湍流和小尺度湍流，通过求解三维经过修正的Navier-Stokes方程，得到大涡旋的运动特性，而对小涡旋运动还采用上述的模型。 实际求解中，选用什么模型要根据具体问题的特点来决定。选择的一般原则是精度要高，应用简单，节省计算时间，同时也具有通用性。 FLUENT提供的湍流模型包括：单方程（Spalart-Allmaras）模型、双方程模型（标准κ-ε模型、重整化群κ-ε模型、可实现(Realizable)κ-ε模型）及雷诺应力模型和大涡模拟。
包含更多 物理机理
Zero-Equation Models RANS-based
models One-Equation Models
Spalart-Allmaras Two-Equation Models
Standard k-? 每次迭代
RNG k-??计算量增加 ?????Realizable k-??
FLUENT提Reynolds-Stress Model 供的模型选 择 Large-Eddy Simulation
Direct Numerical Simulation
湍流模型种类示意图 xxv
第二节 平均量输运方程
雷诺平均就是把Navier-Stokes方程中的瞬时变量分解成平均量和脉动量两部分。对于速度，有： ui?ui?ui?
3－3 其中，ui和ui?分别是平均速度和脉动速度（i=1,2,3） 类似地，对于压力等其它标量，我们也有： ??????
3－4 其中，?表示标量，如压力、能量、组分浓度等。 把上面的表达式代入瞬时的连续与动量方程，并取平均（去掉平均速度ui上的横线），我们可以把连续与动量方程写成如下的笛卡儿坐标系下的张量形式： ????(?ui)?0
3－5 ?t?xi?DuDti?uj?ul??????ui2???????????ij??ui?u?j
????xi?xj???x?x3?x?xjil??j????p??上面两个方程称为雷诺平均的Navier-Stokes（RANS）方程。他们和瞬时Navier-Stokes方程有相同的形式，只是速度或其它求解变量变成了时间平均量。额外多出来的项??ui?u?j是雷诺应力，表示湍流的影响。如果要求解该方程，必须模拟该项以封闭方程。 如果密度是变化的流动过程如燃烧问题，我们可以用法夫雷（Favre）平均。这样才可以求解有密度变化的流动问题。法夫雷平均就是出了压力和密度本身以外，所有变量都用密度加权平均。变量的密度加权平均定义为： ~????/?
3－7 符号～表示密度加权平均；对应于密度加权平均值的脉动值用???表示，即有：~???????。很显然，这种脉动值的简单平均值不为零，但它的密度加权平均值等于零，即： ????0， ?????0
Boussinesq近似与雷诺应力输运模型
为了封闭方程，必须对额外项雷诺应力??ui?u?j进行模拟。一个通常的方法是应用Boussinesq假设，认为雷诺应力与平均速度梯度成正比，即：
??ui?u?j??t??2?ui??(?k??)?ij
3－8 t?3?xi?Boussinesq假设被用于Spalart-Allmaras单方程模型和k??双方程模型。Boussinesq近??ui?uj???x?xij?似的好处是与求解湍流粘性系数有关的计算时间比较少，例如在Spalart-Allmaras单方程模型中，只多求解一个表示湍流粘性的输运方程；在k??双方程模型中，只需多求解湍动能k和耗散率ε两个方程，湍流粘性系数用湍动能k和耗散率ε的函数。Boussinesq假设的缺点是认为湍流粘性系数?t是各向同性标量，对一些复杂流动该条件并不是严格成立，所以具有其应用限制性。 另外的方法是求解雷诺应力各分量的输运方程。这也需要额外再求解一个标量方程，通常是耗散率ε方程。这就意味着对于二维湍流流动问题，需要多求解4个输运方程，而三维湍流问题需要多求解7个方程，需要比较多的计算时间，对计算机内存也有更高要求。
xxvi 在许多问题中，Boussinesq近似方法可以得到比较好的结果，并不一定需要花费很多时间来求解雷诺应力各分量的输运方程。但是，如果湍流场各向异性很明显，如强旋流动以及应力驱动的二次流等流动中，求解雷诺应力分量输运方程无疑可以得到更好的结果。
第三节 湍流模型
3.3.1 单方程（Spalart-Allmaras）模型 ~Spalart-Allmaras模型的求解变量是?，表征出了近壁（粘性影响）区域以外的湍流运~动粘性系数。?的输运方程为： 1????G???Dt??~??xj?~D?~????~?(????)?xj??~???????Cb2????x?j???????Y
3－9 ?????其中，G?是湍流粘性产生项；Y?是由于壁面阻挡与粘性阻尼引起的湍流粘性的减少；??~和Cb2是常数；ν是分子运动粘性系数。 湍流粘性系数用如下公式计算： ~f ?t????1其中，f?1是粘性阻尼函数，定义为：f?1??333??C?1，并且??~??。 湍流粘性产生项，G?用如下公式模拟： ~~
3－10 G??Cb1?S?~其中，S?S?~?2kd2f?2，而f?2?1??1??f?1。其中，Cb1和k是常数，d是计算点到壁面的距离；S??ij2?ij?ij。?ij定义为： ??
3－11 ???ui1??uj???2??xj??xi由于平均应变率对湍流产生也起到很大作用，FLUENT处理过程中，定义S为： S??ij?Cprodmin(0,Sij??ij)
3－12 其中，Cprod?2.0，?ij?Sij?ui1??uj???2??xj??xi?ij?ij，Sij?2SijSij，平均应变率Sij定义为： ??
3－13 ??在涡量超过应变率的计算区域计算出来的涡旋粘性系数变小。这适合涡流靠近涡旋中心的区域，那里只有“单纯”的旋转，湍流受到抑止。包含应变张量的影响更能体现旋转对湍流的影响。忽略了平均应变，估计的涡旋粘性系数产生项偏高。 湍流粘性系数减少项Y?为： ~2???Y??Cw1?fw??
3－14 ?d?其中，fw6?1?Cw?3?g?66??g?Cw3?61/6
g?r?Cw2(r?r)
3－16 ~?r?~22
3－17 Skd xxvii ~其中，Cw1，Cw2，Cw3是常数，S?S?~?kd22f?2。在上式中，包括了平均应变率~对S的影响，因而也影响用S计算出来的r。 上面的模型常数在FLUENT中默认值为：Cb1?0.1335，Cb2?0.622，??~?2/3，C?1?7.1，Cw1?Cb1/k2?(1?Cb2)/??~，Cw2?0.3，Cw3?2.0，k?0.41。
~在壁面，湍流运动粘性?设置为零。当计算网格足够细，可以计算层流底层时，壁面切应力用层流应力－应变关系求解，即： ?u?yu?
3－18 u??如果网格粗错不能用来求解层流底层，则假设与壁面近邻的网格质心落在边界层的对数区，则根据壁面法则： uu????u?y??lnE???
3－19 k???1其中，k=0.419，E=9.793。
对流传热传质模型
在FLUENT中，用雷诺相似湍流输运的概念来模拟热输运过程。给出的能量方程为： ??t(?E)???xi[ui(?E?p)]?cp?t????k????xi?Prt????T???x?uj(?ij)eff?i???Sh
3－20 ??式中，E是总能量，(?ij)eff是偏应力张量，定义为： (?ij)eff??eff(?uj?xi??ui?xj)?23?eff?ui?xi?ij
3－21 其中，(?ij)eff表示粘性加热，耦合求解。如果默认为分开求解，FLUENT不求解处(?ij)eff。但是可以通过变化“粘性模型”面板上的湍流普朗特数（Prt），其默认值为0.85。 湍流质量输运与热输运类似，默认的Schmidt数是0.7，该值同样也可以在“粘性模型”面板上调节。 标量的壁面处理与动量壁面处理类似，分别选用合适的壁面法则。 综上所述，Spalart-Allmaras模型是相对简单的单方程模型，只需求解湍流粘性的输运方程，并不需要求解当地剪切层厚度的长度尺度。该模型对于求解有壁面影响流动及有逆压力梯度的边界层问题有很好模拟效果，在透平机械湍流模拟方面也有较好结果。 Spalart-Allmaras模型的初始形式属于对低雷诺数湍流模型，这必须很好解决边界层的粘性影响区求解问题。在FLUENT中，当网格不是很细时，采用壁面函数来解决这一问题。当网格比较粗糙时，网格不满足精确的湍流计算要求，用壁面函数也许是最好的解决方案。另外，该模型中的输运变量在近壁处的梯度要比k??中的小，这使得该模型对网格粗糙带来数值误差不太敏感。 但是，Spalart-Allmaras模型不能预测均匀各向同性湍流的耗散。并且，单方程模型没有考虑长度尺度的变化，这对一些流动尺度变换比较大的流动问题不太适合。比如，平板射流问题，从有壁面影响流动突然变化到自由剪切流，流场尺度变化明显。
xxviii 3.3.2 标准k??模型
标准k??模型需要求解湍动能及其耗散率方程。湍动能输运方程是通过精确的方程推导得到，但耗散率方程是通过物理推理，数学上模拟相似原形方程得到的。该模型假设流动为完全湍流，分子粘性的影响可以忽略。因此，标准k??模型只适合完全湍流的流动过程模拟。 标准k??模型的湍动能k和耗散率ε方程为如下形式： ?DkDtD???t????????xi???k???k??
3－22 ??x??Gk?Gb????YM?i?2???????
3－23 ??x??C1?k(Gk?C3?Gb)?C2??k?i??t?????????Dt?xi???k?在上述方程中，Gk表示由于平均速度梯度引起的湍动能产生，Gb是用于浮力影响引起的湍动能产生；YM可压速湍流脉动膨胀对总的耗散率的影响。湍流粘性系数k2?t??C??在FLUENT中，作为默认值常数，C1?＝1.44，C2?=1.92，C??0.09，湍动能k与。 耗散率ε的湍流普朗特数分别为?k＝1.0，??＝1.3。可以通过调节“粘性模型”面板来调节这些常数值。
3.3.3 重整化群κ-ε模型
重整化群κ-ε模型是对瞬时的Navier-Stokes方程用重整化群的数学方法推导出来的模型。模型中的常数与标准κ-ε模型不同，而且方程中也出现了新的函数或者项。其湍动能与耗散率方程与标准κ-ε模型有相似的形式： ???k??????k?eff???Gk?Gb????YM
3－24 Dt?xi??xi?Dk2??????????R
?????eff???C1?(Gk?C3?Gb)?C2??Dt?xi??xi?kkGk表示由于平均速度梯度引起的湍动能产生，Gb是用于浮力影响引起的湍动能产生；YMD?可压速湍流脉动膨胀对总的耗散率的影响，这些参数与标准κ-ε模型中相同。?k和??分别是湍动能k和耗散率?的有效湍流普朗特数的倒数。
湍流粘性系数计算公式为： ~??2k??~??1.72d?d?
????3~??1?C???~其中，???eff/?，C??100 对上面方程积分，可以精确得到有效雷诺数（涡旋尺度）对湍流输运的影响，这有助于处理低雷诺数和近壁流动问题的模拟。 对于高雷诺数，上面方程可以给出：?t??C?k2?，C??0.0845。这个结果非常有意思，和标准κ-ε模型的半经验推导给出的常数C??0.09非常近似。 在FLUENT中，如果是默认设置，用重整化群κ-ε模型时候是针对的高雷诺数流动 xxix 3.3.2 标准k??模型
标准k??模型需要求解湍动能及其耗散率方程。湍动能输运方程是通过精确的方程推导得到，但耗散率方程是通过物理推理，数学上模拟相似原形方程得到的。该模型假设流动为完全湍流，分子粘性的影响可以忽略。因此，标准k??模型只适合完全湍流的流动过程模拟。 标准k??模型的湍动能k和耗散率ε方程为如下形式： ?DkDtD???t????????xi???k???k??
3－22 ??x??Gk?Gb????YM?i?2???????
3－23 ??x??C1?k(Gk?C3?Gb)?C2??k?i??t?????????Dt?xi???k?在上述方程中，Gk表示由于平均速度梯度引起的湍动能产生，Gb是用于浮力影响引起的湍动能产生；YM可压速湍流脉动膨胀对总的耗散率的影响。湍流粘性系数k2?t??C??在FLUENT中，作为默认值常数，C1?＝1.44，C2?=1.92，C??0.09，湍动能k与。 耗散率ε的湍流普朗特数分别为?k＝1.0，??＝1.3。可以通过调节“粘性模型”面板来调节这些常数值。
3.3.3 重整化群κ-ε模型
重整化群κ-ε模型是对瞬时的Navier-Stokes方程用重整化群的数学方法推导出来的模型。模型中的常数与标准κ-ε模型不同，而且方程中也出现了新的函数或者项。其湍动能与耗散率方程与标准κ-ε模型有相似的形式： ???k??????k?eff???Gk?Gb????YM
3－24 Dt?xi??xi?Dk2??????????R
?????eff???C1?(Gk?C3?Gb)?C2??Dt?xi??xi?kkGk表示由于平均速度梯度引起的湍动能产生，Gb是用于浮力影响引起的湍动能产生；YMD?可压速湍流脉动膨胀对总的耗散率的影响，这些参数与标准κ-ε模型中相同。?k和??分别是湍动能k和耗散率?的有效湍流普朗特数的倒数。
湍流粘性系数计算公式为： ~??2k??~??1.72d?d?
????3~??1?C???~其中，???eff/?，C??100 对上面方程积分，可以精确得到有效雷诺数（涡旋尺度）对湍流输运的影响，这有助于处理低雷诺数和近壁流动问题的模拟。 对于高雷诺数，上面方程可以给出：?t??C?k2?，C??0.0845。这个结果非常有意思，和标准κ-ε模型的半经验推导给出的常数C??0.09非常近似。 在FLUENT中，如果是默认设置，用重整化群κ-ε模型时候是针对的高雷诺数流动 xxix 问题。如果对低雷诺数问题进行数值模拟，必须进行相应的设置。
重整化群κ-ε模型有旋修正
通常，平均运动有旋时候对湍流有重要影响。FLUENT中重整化群κ-ε模型通过修正湍流粘性系数来考虑了这类影响。 湍流粘性的修正形式为： ?t??t0f(?s,?,k?)
3－27 其中，?t0是不考虑有旋计算出来的湍流粘性系数；Ω是FLUENT计算出来的特征旋流数；?s是旋流常数，不同值表示有旋流动的强度不同。流动可以是强旋或者中等旋度的。FLUENT默认设置?s＝0.05,针对中等旋度的流动问题，对于强旋流动，可以选择较大的值。 湍动能及其耗散率的有效湍流普朗特数倒数的计算公式为： ??1..??2..??mol?eff
3－28 式中，?0＝1，在高雷诺数流动问题中，?mol/?eff??1，?k????1.393。 湍流耗散率方程右边的R为： R?C???(1??/?0)?1???332k
3－29 其中，??Sk/?，?0?4.38，??0.012。 为了更清楚体现R对耗散率的影响，我们把耗散率输运方程重写为： 22????????*???C2?? 3－30 ?????eff???C1?(Gk?C3?Gb)?C2??Dt?xi??xi?kkkD?*则：C2??C2??C???(1??/?0)1???33
3－31 **在???0的区域，R的贡献为正；C2?大于C2?。以对数区为例，??3，C2??2.0，这和标准κ-ε模型中给出的C2?＝1.92接近。因此，对于弱旋和中等旋度的流动问题，重整化群κ-ε模型给出的结果比标准κ-ε模型的结果要大。 重整化群模型中，C1??1.42，C2??1.68。
3.3.4 可实现κ-ε模型
可实现κ-ε模型的湍动能及其耗散率输运方程为： ?t??????????Dt?xj??k??DkD???k???xj????Gk?Gb????YM
??2?t?????????????????CS???C?CC3?Gb
???121??Dt?xj??t??xj?kk??????其中，C1?max?0.43,????，??Sk/? ??5?? xxx 在上述方程中，Gk表示由于平均速度梯度引起的湍动能产生，Gb是用于浮力影响引起的湍动能产生；?k，??YM可压速湍流脉动膨胀对总的耗散率的影响。C2和C1?是常数；分别是湍动能及其耗散率的湍流普朗特数。在FLUENT中，作为默认值常数，C1?＝1.44，C2=1.9，?k＝1.0，??＝1.2。 可实现κ-ε模型的湍动能的输运方程与标准κ-ε模型和重整化群κ-ε模型有相同的形式，只是模型参数不同。但耗散率方程有较大不同。首先耗散率产生项（方程右边第二项）不包含湍动能产生项Gk，现在的形式更能体现能量在谱空间的传输。另外的特色在于耗散率减少项中，不具有奇异性。并不象标准κ-ε模型模型那样把K放在分母上。 该模型适合的流动类型比较广泛，包括有旋均匀剪切流，自由流（射流和混合层），腔道流动和边界层流动。对以上流动过程模拟结果都比标准κ-ε模型的结果好，特别是可实现κ-ε模型对圆口射流和平板射流模拟中，能给出较好的射流扩张角。 ?κ-ε模型中，C?不再是个常数，而是通过如下公式计算： C??1A0?AsUK*湍流粘性系数公式为?t??C?k2，这和标准κ-ε模型相同。不同的是，在可实现
3－34 其中，U*??~~~?ij??ij??ijk?k，?ij??ij?2?ijk?k，the mean ?ij是 is SijSij??ij?ij，6cos?，而： SijSij，Sij?ui1?uj?(?) 2?xi?xjrate-of –rotation tensor viewed in a rotating reference frame with the angular velocity ?k。模型常数A0?4.04，As???13arccos(6W)，式中W＝SijS~SjkSkj~，S?我们可以发现，C?是平均应变率与旋度的函数。在平衡边界层惯性底层，可以得到C?＝0.09,与标准κ-ε模型中采用底常数一样。
双方程模型中，无论是标准κ-ε模型、重整化群κ-ε模型还是可实现κ-ε模型，三个模型有类似的形式，即都有κ和ε的输运方程，它们的区别在于：1，计算湍流粘性的方法不同；2，控制湍流扩散的湍流Prandtl数不同；3，ε方程中的产生项和Gk关系不同。但都包含了相同的表示由于平均速度梯度引起的湍动能产生Gk，用于浮力影响引起的湍动能产生Gb；可压速湍流脉动膨胀对总的耗散率的影响YM。 湍动能产生项 Gk???ui?u?j?uj?xi
Gb??gi?t?TPrt?xi式中，Prt是能量的湍流普特朗数，对于可实现κ-ε模型，默认设置值为0.85。对于重整化群κ-ε模型，Prt?1/?，??1/Pr?k/?Cp。热膨胀系数???1??????，对???T?p于理想气体，浮力引起的湍动能产生项变为： ?t??Gb??gi
3－37 ?Prt?xi xxxi 在FLUENT程序中，如果有重力作用，并且流场里有密度或者温度的梯度，浮力对湍动能的影响都是存在的。浮力对耗散率的影响不是很清楚，因此，默认设置中，耗散率方程中的浮力影响不被考虑。如果要考虑浮力对耗散率的影响，用“粘性模型”面板来控制。浮力对耗散率影响是用C3?来体现。但C3?并不是常数，而是如下的函数形式： C3??tanhvu
v是平行于重力方向的速度分量；u是垂直于重力方向的速度分量。如果流动速度与重力方向相同的剪切流动，C3?＝1，对于流动方向与重力方向垂直的剪切流，C3?＝0。
对于高马赫数的流动问题，可压速性对湍流影响在YM中体现。 YM???2M2t k其中，Mt是马赫数，定义为：Mt?（a??RT是声速）。 2a默认设置中，只要选择可压速理想气体，可压速效应都是考虑的。 在上述的双方程模型中，对流传热传质模型都是通过雷诺相似湍流动量输运方程得到的。能量方程形式为： ??t(?E)???xi?ui??E?p???式中，E是总的能量，keff(?ij)eff??eff????T?keff??Sh
?u(?)jijeff???xi??xi?是有效导热系数；(?ij)eff是偏应力张量，定义为： ?2?ui????
3－40 ?3eff?xiji???uj?ui????x?xji?(?ij)eff表示的是粘性加热，耦合求解时总是计算。如果不是耦合求解时候，作为默认设置，并不求解该量。如果有需要，需在“粘性模型”面板中设置。 对于重整化群κ-ε模型，有效导热系数为： keff??cp?eff
3－41 α用（3－28）计算，式中，?0?1/Pr?k/?Cp。事实上，?随着?mol/?eff的变化而变化，这是重整化群κ-ε模型的一个优点，因为实验中证明，湍流普朗特数随分子普朗特数及湍流而变化。 湍流质量输运处理过程与能量输运过程类似。对于标准κ-ε模型和可实现的κ-ε模型，默认的Schmidt数是0.7，重整化群模型中，是通过方程3－28来计算的，其中，?0?1/Sc，Sc是分子Schimidt数。
3.3.5 雷诺应力模型（RSM）
雷诺应力模型是求解雷诺应力张量的各个分量的输运方程。具体形式为： ??t(?uiuj)???xk(?Ukuiuj)?
对流项Cij ???xk??uuiu?p(?kjui??ikuj)?jkT??????uuij? ??xk??xk?L湍流扩撒项Dij
分子扩散Dij
xxxii ?Uj??Ui?????giuj??gjui? ???uu?uujk?ik?x?xk?k??应力产生项Pij
浮力产生项目Gij ????ui?uj??ui?uj?? ?p??2???x??xi??xk?xkj?压力应变项?ij
?2??kujum?ikm?uium??jkm?
3－42 系统旋转产生项Fij 上面方程中，Cij，Dij，Pij，Fij不需要模拟，而Dij，Gij，?ij，?ij需要模拟以封闭方程。下面简单对几个需要模拟项的模拟。 Dij可以用Delay and Harlow [L38]的梯度扩散模型来模拟，即： DTijLTT??kukul?uiu?Cs???xk???xl??uiu????t??xk??k?xk?j??
3－43 ??但这个模型会导致数值不稳定，因此FLUENT程序中采用标量湍流扩散模型： DTijj??
???k??0.82，这和标准κ-ε模型中选取1.0有所不同。 式中，湍流粘性系数用?t??C?k2来计算，根据Lien and Leschziner [L98]，根据Gibson and Launder [L58], Fu [L55], Launder [L88,L89], 压力应变项?ij可以分解为三项，即： ?ij??ij,1??ij,2??ij
3－45 ?ij,1，?ij,2和?wij分别是慢速项，快速项和壁面反射项。 w?ij,1??C1???2?，常数C1?1.8。 uu??kijij?k?3??12???ij,2??C2??Pij?Fij?Gij?Cij???ij?P?G?C??，C2?0.60，P?Pkk，23??11G?Gkk，C?Ckk。 22壁面反射项用于重新分布近壁的雷诺正应力分布，主要是减少垂直于壁面的雷诺正应力，增加平行于壁面的雷诺正应力。该项模拟为： ?wij?C1????ukumnknm?ijk??k?uiuknjnk?ujuknink? 22C?d?l?knn
3－46 jk,2ik??Cl?d3/????km,2nknm?ij??ik,2njnk???C222???0.3，nk是xk在垂直于壁面方向上的单位分量，d是到壁面的式中，C1?＝0.5，C2距离；Cl?C?3/4/k，C??0.09，k=0.41。 wij默认设置时候，FLUENT不计算? 。如果需要计算时候，在“粘性模型”面板中设xxxiii ?Uj??Ui?????giuj??gjui? ???uu?uujk?ik?x?xk?k??应力产生项Pij
浮力产生项目Gij ????ui?uj??ui?uj?? ?p??2???x??xi??xk?xkj?压力应变项?ij
?2??kujum?ikm?uium??jkm?
3－42 系统旋转产生项Fij 上面方程中，Cij，Dij，Pij，Fij不需要模拟，而Dij，Gij，?ij，?ij需要模拟以封闭方程。下面简单对几个需要模拟项的模拟。 Dij可以用Delay and Harlow [L38]的梯度扩散模型来模拟，即： DTijLTT??kukul?uiu?Cs???xk???xl??uiu????t??xk??k?xk?j??
3－43 ??但这个模型会导致数值不稳定，因此FLUENT程序中采用标量湍流扩散模型： DTijj??
???k??0.82，这和标准κ-ε模型中选取1.0有所不同。 式中，湍流粘性系数用?t??C?k2来计算，根据Lien and Leschziner [L98]，根据Gibson and Launder [L58], Fu [L55], Launder [L88,L89], 压力应变项?ij可以分解为三项，即： ?ij??ij,1??ij,2??ij
3－45 ?ij,1，?ij,2和?wij分别是慢速项，快速项和壁面反射项。 w?ij,1??C1???2?，常数C1?1.8。 uu??kijij?k?3??12???ij,2??C2??Pij?Fij?Gij?Cij???ij?P?G?C??，C2?0.60，P?Pkk，23??11G?Gkk，C?Ckk。 22壁面反射项用于重新分布近壁的雷诺正应力分布，主要是减少垂直于壁面的雷诺正应力，增加平行于壁面的雷诺正应力。该项模拟为： ?wij?C1????ukumnknm?ijk??k?uiuknjnk?ujuknink? 22C?d?l?knn
3－46 jk,2ik??Cl?d3/????km,2nknm?ij??ik,2njnk???C222???0.3，nk是xk在垂直于壁面方向上的单位分量，d是到壁面的式中，C1?＝0.5，C2距离；Cl?C?3/4/k，C??0.09，k=0.41。 wij默认设置时候，FLUENT不计算? 。如果需要计算时候，在“粘性模型”面板中设xxxiii 置。
线性压力应变模型
对于小雷诺数流动，特别是用双层模型求解近壁流动问题时，FLUENT中通过改进模?来改进压力应变项Launder [L91]。这一过程只有在选择双层流型常数C1，C2，C1?和C2模型时候，在“粘性模型”面板上调节。 C1?1?2.58AC2?0.75A ?(0.?e?Ret)2?? ?C1???23C1?1.67 1??2?C2??6? ??max?3C2,0C2??????其中，Ret??k2/(??)，参数A和张量不变量A2, A3定义为： 9??A??1?(A2?A3)? 8??A2??ik?ki A3??ik?kj?ji 式中，?ij是雷诺应力张量各向异性部分，定义为： 2????uiuj??k?ij3????k??????
3－47 ????ij 二阶压力应变模型
二阶压力应变模型由Spezible {L157}等人提出。 ?ij??(C1???C1P)bij?C2??(bikbkj??C4?k(bikS?bjkSik?23*13bmnbmn?ij)?(C3?C3*bijbij)?kSij jkbmnSmn?ij)?C5?k(bik?jk?bjk?ik)
3－48 式中，bij是雷诺应力各向异性张量，定义为： 2???uu??k?ij?ij3
bij???2?k??????
?????ui1??uj?；????ij?2??xj???xi??；模型常数??平均应变率Sij定义为: Sij*?ui1??uj???2??xj??xi*C1?3.4，C1?1.8，C2?4.2，C3?0.8，C3?1.3，C4?1.25，C5?0.4。 二阶压力应变模型不需要考虑壁面反射影响去模拟对数区湍流边界层过程。
浮力对湍流的影响
浮力引起的产生项模拟为： Gij?T?T?g???gijPrt??x?xij??t???
3－50 ??其中，Prt是能量的湍流普朗特数，默认设置值为0.85。 对于理想气体，把热膨胀系数的定义代入上式，得： Gij??????g???gij?Prt??x?xij??t??
?? 耗散项?ij的模拟 耗散张量?ij模拟为： ?ij?23?ij(???YM)
3－52 2式中，YM???2Mt，Mt是马赫数；标量耗散率?用标准k-?模型中的采用的耗散率输运方程求解。
雷诺应力模型的边界条件
在流场进口，雷诺应力模型需要各个雷诺应力分量和湍动能耗散率的值。这些值可以直接输入，也可以湍流强度和特征长度来计算。 在壁面，雷诺应力模型通过壁面函数，给出各个雷诺应力分量和耗散率的值。
雷诺应力模型的能量与质量输运方程
在雷诺应力模型中，对流传热传质模型都是通过雷诺相似湍流动量输运方程得到的。能量方程形式为： ??t(?E)???xi?ui??E?p???cp?t?T???(k?)?uj(?ij)eff??xi?Prt?xi???Sh
3－53 ??式中，E是总的能量；(?ij)eff是偏应力张量，定义为： ??uj?ui????x?xji??2?ui????
3－54 ?3eff?xiji?(?ij)eff??eff(?ij)eff表示的是粘性加热，耦合求解时总是计算。如果不是耦合求解时候，作为默认设置，并不求解该量，并且Prt=0.85。如果有需要，需在“粘性模型”面板中设置。
3.3.6 大涡模拟（LES）
湍流中包含了不同时间与长度尺度的涡旋。最大长度尺度通常为平均流动的特征长度尺度。最小尺度为Komogrov尺度。 LES的基本假设：1，动量、能量、质量及其它标量主要由大涡输运；2，流动的几何和边界条件决定了大涡的特性，而流动特性主要在大涡中体现；3，小尺度涡旋受几何和边界条件影响较小，并且各向同性；大涡模拟过程中，直接求解大涡，小尺度涡旋模拟，从而使得网格要求比DNS低。
3.3.6.1大涡模拟的控制方程
LES的控制方程是对Navier-Stokes方程在波数空间或者物理空间进行过滤得到的。过滤的过程是去掉比过滤宽度或者给定物理宽度小的涡旋，从而得到大涡旋的控制方程。 过滤变量（上横线）定义为： ?(x)??D?(x?)G(x,x?)dx?
其中，D表示流体区域；G是决定涡旋大小的过滤函数。 在FLUENT中，有限控制体离散本身暗中包括了过滤运算， ?(x)?1V?V?(x?)dx?， x??V
3－56 其中V是计算控制体体积，过滤函数为： x??V?1/VG(x,x?)??
3－57 x??V?0目前，大涡模拟对不可压流动问题得到较多应用，但在可压缩问题中的应用还很少，因此这里涉及的理论都是针对不可压流动的大涡模拟方法。在FLUENT中，大涡模拟只能针对不可压流体（当然并非说是密度是常数）的流动。 过滤不可压的Navier-Stokes方程后，可以得到LES控制方程： ???t???ui?xi?0
??t(?ui)???xj(?uiuj)???xj(??ui?xj)??p?xi???ij?xj
其中，?ij为亚网格应力，定义为： ?ij??uiuj??ui?uj
3－60 很明显，上述方程与雷诺平均方程很相似，只不过大涡模拟中的变量是过滤过的量，而非时间平均量，并且湍流应力也不同。
3.3.6.2 亚网格模型
由于LES中亚网格应力项是未知的，并且需要模拟以封闭方程。目前，采用比较多的亚网格模型为涡旋粘性模型，形式为： ?ij?13?kk?ij??2?tSij
3－61 式中，?t是亚网格湍流粘性系数；Sij是求解尺度下的应变率张量，定义为： Sij?uj?ui1????2??xj?xi???
3－62 ??求解亚网格湍流粘性系数?t时，FLUENT提供了两种方法。第一，Smagorinsky-Lilly模型；第二，基于重整化群的亚网格模型。 最基本的亚网格模型是Smagorinsky [L145]最早提出的，Lilly [L99]把它进行了改善，这就是今天的Smagorinsky-Lilly模型。该模型的涡粘性计算方程为： ?t??LsS
2式中，Ls是亚网格的混合长度；S?度Ls可以用下式计算。 Ls?min(kd,CsV1/32SijSij。Cs是Smagorinsky常数，则亚网格混合长)
3－64 其中，k=0.42，d是到最近壁面的距离，V是计算控制体体积。 Lilly通过对均匀各向同性湍流惯性子区湍流分析，得到了Cs＝0.23。但是研究中发现，对于有平均剪切或者过渡流动中，该系数过高估计了大尺度涡旋的阻尼作用。因此，对于比较多的流动问题，Cs＝0.1有比较好的模拟结果，该值是FLUENT的默认设置值。 我们再来看看基于重整化群思想的亚网格模型。人们用重整化群理论推导出了亚网格涡旋粘性系数[L182]，该方法得到的是亚网格有效粘性系数，?eff????t，而 ?eff???s2?eff??????1?H?C????3???????1/31/3
3－65 式中，?s?(CrngV)22SijSij，H(x)是Heaviside函数， x?0?xH(x)??
3－66 0x?0?V是计算控制体体积；重整化群常数Crng?0.157,而常数C=100。 对于高雷诺数流动（?t??? ），?eff??t，基于重整化群理论的亚网格模型就与Smagorinsky-Lilly模型相同，只是模型常数有区别。在流动场的低雷诺数区域，上面的函数就小于零，从而只有分子粘性起作用。所以，基于重整化群理论的亚网格模型对流动转捩和近壁流动问题有较好模拟效果。
3.3.6.3 大涡模拟的边界条件
对于给定进口速度边界条件，速度等于各个方向分量与随机脉动量的和，即： ui??ui??I?u
3.3.6.2 亚网格模型
由于LES中亚网格应力项是未知的，并且需要模拟以封闭方程。目前，采用比较多的亚网格模型为涡旋粘性模型，形式为： ?ij?13?kk?ij??2?tSij
3－61 式中，?t是亚网格湍流粘性系数；Sij是求解尺度下的应变率张量，定义为： Sij?uj?ui1????2??xj?xi???
3－62 ??求解亚网格湍流粘性系数?t时，FLUENT提供了两种方法。第一，Smagorinsky-Lilly模型；第二，基于重整化群的亚网格模型。 最基本的亚网格模型是Smagorinsky [L145]最早提出的，Lilly [L99]把它进行了改善，这就是今天的Smagorinsky-Lilly模型。该模型的涡粘性计算方程为： ?t??LsS
2式中，Ls是亚网格的混合长度；S?度Ls可以用下式计算。 Ls?min(kd,CsV1/32SijSij。Cs是Smagorinsky常数，则亚网格混合长)
3－64 其中，k=0.42，d是到最近壁面的距离，V是计算控制体体积。 Lilly通过对均匀各向同性湍流惯性子区湍流分析，得到了Cs＝0.23。但是研究中发现，对于有平均剪切或者过渡流动中，该系数过高估计了大尺度涡旋的阻尼作用。因此，对于比较多的流动问题，Cs＝0.1有比较好的模拟结果，该值是FLUENT的默认设置值。 我们再来看看基于重整化群思想的亚网格模型。人们用重整化群理论推导出了亚网格涡旋粘性系数[L182]，该方法得到的是亚网格有效粘性系数，?eff????t，而 ?eff???s2?eff??????1?H?C????3???????1/31/3
3－65 式中，?s?(CrngV)22SijSij，H(x)是Heaviside函数， x?0?xH(x)??
3－66 0x?0?V是计算控制体体积；重整化群常数Crng?0.157,而常数C=100。 对于高雷诺数流动（?t??? ），?eff??t，基于重整化群理论的亚网格模型就与Smagorinsky-Lilly模型相同，只是模型常数有区别。在流动场的低雷诺数区域，上面的函数就小于零，从而只有分子粘性起作用。所以，基于重整化群理论的亚网格模型对流动转捩和近壁流动问题有较好模拟效果。
3.3.6.3 大涡模拟的边界条件
对于给定进口速度边界条件，速度等于各个方向分量与随机脉动量的和，即： ui??ui??I?u
xxxvii 其中，I是脉动强度，?是高斯随机数，满足??0，???1。 如果网格足够密并可以求解层流底层的流动的话，壁面切应力采用线性应力应变关系，即： uu?uu???u?y?1k
3－67 如果网格不够细，则假定与壁面邻近网格质心落在边界层对数区内，则： ?lnE(?u?y?)
其中，k=0.418，E=9.793。
雷诺平均模型的比较
模型名字 Spalart-Allmaras 标准k?? RNG k?? Realizable k?? 雷诺应力模型
优点 计算量小，对一定复杂程度的 边界层问题有较好效果 应用多，计算量合适，有较多数据积累和相当精度 能模拟射流撞击，分离流，二次流，旋流等中等复杂流动 和RNG模型差不多，还可以模拟圆口射流问题 考虑的物理机理更仔细，包括了湍流各向异性影响 缺点 计算结果没有被广泛测试，缺少子模型，如考虑燃烧或浮力问题 对于流向有曲率变化，较强压力梯度，有旋问题等复杂流动模拟效果欠缺 受到涡旋粘性各向同性假设限制 受到涡旋粘性各向同性假设限制 CPU时间长（2～3倍），动量和湍流量高度耦合。
第四节 湍流模型算例及其设置
湍流模型设置命令：Define-model-viscous
无粘，层流和湍流
湍流模型选项
近壁处理方法选择
附加湍流选项
???? ???????? ???????? ?????? ??xxxix
算例一：有换热的腔道流动问题 adiabatic wall cold air V = 50 fpm T = 0°F 1 ft insulation constant temperature wall T = 100°F 10 ft
步骤： 1， 检查是否湍流：Re?5980 P 1 ft Dh2， 选择低雷诺数湍流模型－RNG k??模型；壁面处理用非平衡壁面函数，考虑压力梯度影响 3， 网格划分：（1）四边形网格；（2）由于在靠近水平板处，垂直方向梯度较大，则近壁网格加密，并保证第一个控制体在对数区内；（3）变化流向网格间距，用于捕捉边界层发展；（4），根据计算结果，自动调节网格，用于进一步计算温度梯度。
计算结果：
Velocity contours BLs on upper & lower surfaces accelerate the core flow
Temperature contours Important that thermal BL was accurately resolved as well P
算例二：圆柱绕流
air V = 4 fps 1 ft 2 ft wall 5 ft 14.5 ft Compute drag coefficient of the cylinder
步骤： 1， 确定雷诺数，ReD?24600 2， 钝体绕流，后面有不稳定的涡旋脱落。采用RNG k??模型，壁面处理是双层区模型； 3， 网格处理：近壁网格加密，由于是双层区模型，需要网格划分到y??1
计算圆柱绕流的涡旋脱落过程
xli 百度搜索“就爱阅读”,专业资料,生活学习,尽在就爱阅读网92to.com,您的在线图书馆百度搜索“就爱阅读”,专业资料、生活学习,尽在就爱阅读网92to.com,您的在线图书馆!
欢迎转载：
相关推荐：}
我要回帖
更多推荐
- ·笔记本电脑能抖音直播吗怎么在电脑上直播
- ·谁用过的比较好用滋润的口碑最好的十大身体乳乳推荐下呗？
- ·ted演讲在哪个软件看不同国家看到的内容一样吗？
- ·腾讯健康系统怎样解除申述是换实名吗？
- ·MicrobiomeLabs益葆来的王牌益生好国产品牌哪里有卖是什么？
- ·哪位大神请指教可以指教一下，为什么我用笔刷减淡模式下擦出来的光，跟原图相比差的那么多，而且没有光感，
- ·哪位知道魅族换屏多少钱
- ·电脑卡顿求解答
- ·用心呵护自然之美中心论点美是什么？什么时候注册的？
- ·我在微信买了一天烟 韵达快递3951008477288 为什么查不到
- ·求助，我新买的手机打不开超级看看密码错误怎样处理省电和免打扰
- ·因而前8XBXB通俗的意思登的，然则最进wwW8XBXBCoM偏偏却不管了
- ·iPhone x苹果x运行内存怎么清理升级要多少钱
- ·oracl 多条数据oracle列合并成一列显示
- ·超级空投模式已经结束了吗
- ·求此图如何找一个图片的原图啊！
- ·为什么我的苹果里QQ
- ·苹果软件买错了如何申请退款买voso还可以退款吗
- ·热血传奇手游不悔石机版不悔石怎么获得 不悔石多少元宝
- ·根本不能wwWeee881cOM显示内蓉,该是不可以eee881结着使用了
- ·FLUENT模拟时将圆柱个性字体转换器可粘贴为二维的，进口流量要怎么重新设置
- ·多键盘手机哪款好
- ·曾经的中国手机品牌百家争鸣，为什么现在却越来越少了
- ·路遇街头气功表演，大家一起看看是不是唬人的把戏
- ·连云港映象西班牙怎么样样，映象西班牙最新评论，映象西
- ·长跑后膝盖疼会膝盖疼的原因有哪些
- ·做瑜伽腿部拉伸动作侧身拉伸为何另一边胸腔要往上往后翻
- ·健身房私人教练靠谱吗的私教竟是这些人，怎样找到靠谱的私教
- ·美利达勇士600多少钱付利威300多少钱
- ·2.0运动助手版入手哪个比较好
- ·怎样配合减脂期力量训练练减脂
- ·广西充气游泳池供应商白殿风夏天可以去游泳吗？有影响吗？
- ·中国男乒夺得冠军队内是哪些球员掳虑÷宸颍匪髀
- ·求助，关于职业拳击拳王金腰带保留的规则
- ·哪里是网球时时能引起高血压吗
