整理fluent求解器资料.docx

1、整理fluent求解器资料(完整版)fluent求解器资料 编辑整理：尊敬的读者朋友们：这里是精品文档编辑中心，本文档内容是由我和我的同事精心编辑整理后发布的，发布之前我们对文中内容进行仔细校对，但是难免会有疏漏的地方，但是任然希望（(完整版)fluent求解器资料）的内容能够给您的工作和学习带来便利。同时也真诚的希望收到您的建议和反馈，这将是我们进步的源泉，前进的动力。本文可编辑可修改，如果觉得对您有帮助请收藏以便随时查阅，最后祝您生活愉快 业绩进步，以下为(完整版)fluent求解器资料的全部内容。(完整版)fluent求解器资料 编辑整理：张嬗雒老师尊敬的读者朋友们：这里是精品文档编辑中

2、心，本文档内容是由我和我的同事精心编辑整理后发布到文库，发布之前我们对文中内容进行仔细校对,但是难免会有疏漏的地方，但是我们任然希望 (完整版）fluent求解器资料 这篇文档能够给您的工作和学习带来便利.同时我们也真诚的希望收到您的建议和反馈到下面的留言区,这将是我们进步的源泉,前进的动力.本文可编辑可修改，如果觉得对您有帮助请下载收藏以便随时查阅，最后祝您生活愉快 业绩进步,以下为 (完整版)fluent求解器资料 这篇文档的全部内容。压力基求解器在压力基求解器中，控制方程是依次求解的。压力基求解器是从原来的分离式求解器发展来的，按顺序仪次求解动量方程、压力修正方程、能量方程和组分方程及其

3、他标量方程，如湍流方程等，和之前不同的是，压力基求解器还增加了耦合算法,可以自由在分离求解和耦合求解之间转换， 需要注意的是，在压力基求解器中提供的几个物理模型，在密度基求解器中是没有的。这些物理模型包括：流体体积模型（VOF)，多项混合模型，欧拉混合模型,PDF燃烧模型，预混合燃烧模型，部分预混合燃烧模型，烟灰和NOx模型，Rosseland辐射模型,熔化和凝固等相变模型，指定质量流量的周期流动模型，周期性热传导模型和壳传导模型等。 与密度基求解器的区别：区别1：压力基求解器主要用于低速不可压缩流动的求解，而密度基求解器则主要针对高速可压缩流动而设计,但是现在两种方法都已经拓展成为可以求解很

4、大流动速度范围的求解方法.两种求解方法的共同点是都使用有限容积的离散方法，但线性化和求解离散方程的方法不同. 区别2：密度基求解器从原来的耦合求解器发展来的,同时求解连续性方程、动量方程、能量方程和组分方程。然后依次再求解标量方程。（注:密度基求解器不求解压力修正方程，因为其压力是由状态方程得出的）.密度基求解器收敛速度快,需要内存和计算量比压力基求解器要大！ 特点:适用于压力基但不适用于密度基的模型: （1) 空化模型 (2） VOF模型 （3） Mixture多相流模型 （4） Eulerian多相流模型 （5） 非预混燃烧模型 （6） 预混燃烧模型 （7） 部分预混燃烧模型 （8) 组合

5、PDF传输模型密度基求解器（Coupled Sover)是同时fluent求解连续方程、动量方程、能量方程及组分输运方程的耦合方程组，然后逐一地求解湍流标量方程。由于控制方程是非线性的,且相互之间是耦合的,因此,在得到收敛解之前，要经过多轮迭代： 1)根据当前的解的结果，更新所有流动变量.如果计算刚刚开始，则用初始值来更新。 2）同时求解连续方程、动量方程、能量方程及组分输运方程的耦合方程组（后两个方程视需要进行求解） 3)根据需要，逐一地求解湍流、辐射等标量方程。注意在求解之前，方程中用到的有关变量要用到前面得到的结果更新。 4)对于包含离散相的模拟,当内部存在相间耦合时，根据离散相的轨迹计

6、算结果更新连续相的源项。 5)检查方程组是否收敛，若不收敛,回到第1)步,重新计算。 格式:密度基求解器有两种格式：隐式和显式.密度基显式与隐式求解器依次求解额外的标量方程(如湍流和辐射等).两种格式求解器的主要不同点在于对于耦合方程的线性化上.由于隐式格式具有很好的稳定性，因此使用隐式求解器能够比显式格式更快的获得收敛的稳定解。然而，隐式格式要比显式格式消耗更多的内存。具体情况可以查看Fluent理论手册。 特点:对于高速可压流动，由强体积力(如浮力或者旋转力)导致的强耦合流动，或者在非常精细的网格上求解的流动，需要考虑密度基求解器。密度基求解器耦合了流动和能量方程，常常很快便可以收敛.密度

7、基求解器所需要的内存约是压力基求解器的1.5到2倍，选择时可以根据这一情况来权衡利弊.在需要耦合隐式的时候,如果计算机内存不够，就可以采用压力基或密度基显式。密度基显式虽然也耦合了流动和能量方程，但是它还是比密度基隐式需要的内存少，当然它的收敛性也相应差一些。FLUENT求解器设置主要包括:1、压力-速度耦合方程格式选择2、对流插值 3、梯度插值 4、压力插值下面对这几种设置做详细说明.一、压力速度耦合方程求解算法FLUENT中主要有四种算法：SIMPLE，SIMPLEC,PISO，FSM（1）SIMPLE(semi-implicit method for pressurelinked equ

8、ations）半隐式连接压力方程方法,是FLUENT的默认格式.(2）SIMPLEC(SIMPLEconsistent）。对于简单的问题收敛非常快速，不对压力进行修正,所以压力松弛因子可以设置为1(3)Pressure-Implicit with Splitting of Operators (PISO).对非定常流动问题或者包含比平均网格倾斜度更高的网格适用(4）Fractional Step Method (FSM)对非定常流的分步方法.用于NITA格式，与PISO具有相同的特性.二、对流插值（动量方程）FLUENT有五种方法：一阶迎风格式、幂率格式、二阶迎风格式、MUSL三阶格式、QUI

9、CK格式（1）FLUENT默认采用一阶格式。容易收敛,但精度较差，主要用于初值计算。（2)Power Lar。幂率格式,当雷诺数低于5时，计算精度比一阶格式要高.（3）二阶迎风格式.二阶迎风格式相对于一阶格式来说,使用更小的截断误差,适用于三角形、四面体网格或流动与网格不在同一直线上；二阶格式收敛可 能比较慢。(4）MUSL(monotone upstream-centered schemes for conservation laws).当地3阶离散格式.主要用于非结构网格，在预测二次流，漩涡，力等时更精确。（5）QUICK（Quadratic upwind interpolation）格式

10、。此格式用于四边形/六面体时具有三阶精度，用于杂交网格或三角形/四面体时只具有二阶精度。三、梯度插值梯度插值主要是针对扩散项。FLUENT有三种梯度插值方案：greengauss cellbased，Green-gauss node-based，least-quares cell based。（1)格林高斯基于单元体.求解方法可能会出现伪扩散.(2）格林-高斯基于节点. 求解更精确，最小化伪扩散,推荐用于三角形网格上（3）基于单元体的最小二乘法插值。推荐用于多面体网格，与基于节点的格林高斯格式具有相同的精度和格式。四、压力插值压力基分离求解器主要有五种压力插值算法。（1）标准格式（Standa

11、rd)。为FLUENT缺省格式，对大表妹边界层附近的曲线发现压力梯度流动求解精度会降低（但不能用于流动中压力急 剧变化的地方此时应该使用PRESTO!格式代替）(2）PRESTO!主要用于高旋流，压力急剧变化流（如多孔介质、风扇模型等），或剧烈弯曲的区域。（3）Linear（线性格式）。当其他选项导致收敛困难或出现非物理解时使用此格式。(4)second order（二阶格式)。用于可压缩流动,不能用于多孔介质、阶跃、风扇、VOF/MIXTURE多相流.（5）Body Force Weighted体积力。当体积力很大时，如高雷诺数自然对流或高回旋流动中采用此格式。+=首先,所谓的steady和

12、unsteady就表述了流动状态是否随时间变化的含义，这是定性的问题，你自己建立模型一定要明确的。这在计算过程中就是体现在方程不同，unsteady流多了时间变量，那么unsteady流就要进行时间离散。其次,两种求解的结果对比而言：一种情况是你要求解的物理问题是steady的，从理论上来说那么两种求解方式收敛之后的结果都是一样或者近似的.只不过需要注意的是，在unsteady的求解中dt的选取会影响你的计算结果,有可能会计算发散，而且还必须要足够的计算步达到收敛才能和steady的结果进行比较。换句话说，如果你用unsteady的方法去求解steady流,如果计算本身就没收敛，就取结果进行比

13、较，那么肯定是不行的。一般而言,在进行unsteady求解的时候,前面一段时间的计算结果基本上是不予采用的，因为有一个数值收敛的过程。另一种是你要求解的物理问题是unsteady的，那么你用steady的求解方法得出的结果就是一堆垃圾了，没有任何价值。看你的帖子里面说有物体的移动,我不清楚具体的物理模型,但估计应该是unsteady流.generally， the default setting is choosing for solver：FLUENT provides three di erent solver formulations:segregatedcoupled implicit

14、coupled explicit(显式格式主要用于激波等波动解的捕捉问题)The segregated solver traditionally has been used for incompressible and mildly compressible flows。 The coupled approach， on the other hand, was originally designed for highspeed compressible flows。By default, FLUENT uses the segregated solver, for high-speed com

15、pressible flows （as discussed above）, highly coupled flows with strong body forces (e.g.， buoyancy or rotational forces), or flows being solved on very fine meshes, you may want to consider the coupled implicit solver instead.For cases where the use of the coupled implicit solver is desirable, but y

16、our machine does not have sufficient memory, you can use the segregated solver or the coupled explicit solver instead.(explicit save memory use,but need more iterations for converged solution.Choosing the Discretization Scheme1)first-order upwind vs secondorder upwindWhen the flow is aligned with th

17、e grid the firstorder upwind discretization may be acceptable。 For triangular and tetrahedral grids, since the flow is never aligned with the grid， you will generally obtain more accurate results by using the second-order discretization. For quad/hex grids， you will also obtain better results using

18、the second-order discretization， especially for complex flows。 For most cases, you will be able to use the secondorder scheme from the start of the calculation. In some cases, however， you may need to start with the first-order scheme and then switch to the secondorder scheme after a few iterations.

19、 For example， if you are running a highMachnumber flow calculation that has an initial solution much different than the expected final solution， Finally， if you run into convergence diffculties with the second-order scheme， you should try the firstorder scheme instead.2)Quick vs upwind（Quick适用于结构网格,

20、流动方向与网格一致，对于非结构网格推荐用2阶迎风）The QUICK discretization scheme may provide better accuracy than the secondorder scheme for rotating or swirling flows solved on quadrilateral or hexahedral meshes。 For compressible flows with shocks, using the QUICK scheme for all variables， including density， is highly rec

21、ommended for quadrilateral, hexahedral, or hybrid meshes.3）centraldifferencing scheme vs upwindThe central-differencing scheme is available only when you are using the LES turbulence model， and it should be used only when the mesh spacing（网格间距)is fine enough so that the magnitude of the local Peclet

22、 number （Equation 26。25) is less than 1。4)power law vs upwindA power law scheme is also available， but it will generally yield the same accuracy as the first-order scheme.Choosing the Pressure Interpolation Scheme(压力离散格式）a number of pressure interpolation schemes are available when the segregated so

23、lver is used in FLUENT。 For most cases the standard（default) scheme is acceptable， but some types of models may benenit from one of the other schemes:For problems involving large body forces， the body-forceweighted scheme is recommended.For flows with high swirl numbers， high-Rayleighnumber natural

24、convection, highspeed rotating flows， flows involving porous media, and flows in strongly curved domains, use the PRESTO! scheme. 对于可压流，应该使用二阶格式Use the second-order scheme for improved accuracy when one of the other schemes is not applicable。Choosing the Density Interpolation Scheme which is availab

25、le at solve a singlephase compressible flow。If you are calculating a compressible flow with shocks, the first-order upwind scheme may tend to smooth the shocks； you should use the second-orderupwind or QUICK scheme for such flows。Choosing the Pressure-Velocity Coupling Method（压力速度方程耦合方法）SIMPLE vs。 S

26、IMPLECSIMPLE is the default, but many problems will benenit from the use of SIMPLEC， For relatively uncomplicated problems （laminar ows with no additional models activated) in which convergence is limited by the pressurevelocity coupling， you can often obtain a converged solution more quickly using

27、SIMPLEC. With SIMPLEC, the pressurecorrection under-relaxation factor is generally set to 1.0, which aids in convergence speedup。 In some problems， however， increasing the pressure-correction under-relaxation to 1。0 can lead to instability due to high grid skewness. For such cases， you will need to

28、use one or more skewness correction schemes, use a slightly more conservative underrelaxation value （up to 0.7)， or use the SIMPLE algorithm。 The SIMPLEC skewness correction allows FLUENT to obtain a solution on a highly skewed mesh in approximately the same number of iterations as required for a mo

29、re orthogonal mesh。PressureImplicit with Splitting of Operators （PISO)The PISO algorithm with neighbor correction is highly recommended for all transient flow calculations, especially when you want to use a large time step。 （For problems that use the LES turbulence model, which usually requires smal

30、l time steps, using PISO may result in increased computational expense, so SIMPLE or SIMPLEC should be considered instead.） PISO can maintain a stable calculation with a larger time step and an underrelaxation factor of 1。0 for both momentum and pressure。For steady-state problems, PISO with neighbor

31、 correction does not provide any noticeable advantage over SIMPLE or SIMPLEC with optimal underrelaxation factors。When you use PISO neighbor correction， underrelaxation factors of 1.0 or near 1.0 are recommended for all equations.If you use just the PISO skewness correction for highlydistorted meshe

32、s （without neighbor correction)， set the underrelaxation factors for momentum and pressure so that they sum to 1 (e。g。, 0。3 for pressure and 0.7 for momentum)。 If you use both PISO methods, follow the under-relaxation recommendations for PISO neighbor correction, above.Fractional Step MethodThe Fractional Step method (FSM） is av