湍流燃烧模型-PDF
PDF模型
概率密度函数PDF方法以随机的观点来对待湍流问题，对解决湍流化学反应流的问题具有很强的优势。在湍流燃烧中存在一些非输运量( 如反应速率, 密度, 温度及气相体积分数等) 的湍流封闭问题。尽管这些量没有输运方程, 但它们常常是输运变量的已知函数。平均或者过滤高度非线性的化学反应源项会引起方程的封闭问题。因此, 用PDF 的方法来解决这些非输运量的湍流封闭问题显然是一个既简单又直接的途径。
PDF方法是一种较为流行的湍流燃烧模型, 能够较为精确的模拟任何详细的化学动力学过程, 适用于预混、非预混和部分预混的任何燃烧问题。目前, 确定输运变量脉动概率密度函数的方法有输运方程和简化假定两种, 分别称之为输运方程的PDF和简化的PDF。前者建立输运变量脉动的概率密度输运方程, 通过求解该方程来获得输运变量脉动的概率分布。后者假定输运变量脉动的概率密度函数的具体形式, 通过确定其中的一些待定参数来获得输运变量脉动的概率分布。湍流燃烧中, 后者应用最为普遍和广泛。在简化的PDF 中, 输运变量脉动的概率密度函数常常采用双D 分布、截尾高斯分布和B 函数分布等形式。
PDF在理论上可以精确考虑任意详细的化学反应机理，但是其具体求解时需借助其它的模型和算法，而且计算量相对较大。PDF的方程是由N-S方程推导而来，其中的化学反应源项是封闭的，但压力脉动梯度项以及分子粘性和分子扩散引起的PDF的分子输运项是不封闭的, 需要引入模型加以封闭。例如，在速度-标量-湍流频率PDF中，必须采用小尺度混合模型、随机速度模型和湍流频率模型加以封闭。模化后的输运方程难以用有限容积、有限差分和有限元等方法来求解, 比较可行的一种方法是蒙特卡洛(MonteCarlo)方法, 在该方法中输运方程被转化为拉格朗日( Lagrangian)方程, 流体由大量遵循Lagrang ian方程的随机粒子的系统来描述, 最后对粒子作统计平均得到流场物理量和各阶统计矩。另有与有限容积法相结合的蒙特卡洛法。
PDF模型的发展
1969年Lungdren首先推导、计算了速度的联合PDF运输方程，避免了对梯度扩散模型进行模拟，对很简单的流动过程得到了简析解[1]。
随后的Pope等人推导计算了标量的联合PDF运输方程，标量的联合PDF运输方程可以对任意的化学反应机理进行精确的计算[2]。Pope在1984年推出了求解标量PDF方程的蒙特卡洛方法，该方法可以处理多变量联合PDF的能力，包括速度-标量联合PDF运输方程[3]。在1990年Pope和陈义良推出了速度-标量-频率三者联合的PDF方程[4][5]。在1997年，Pope等人提出了速度-湍流频率-压力-内能联合PDF，将其应用于可压缩流体；随后又将其应用于计算超声速湍流反应流[6]。
目前国内外的学者通过合理简化化学反应机理，减少了计算机的计算负荷，使得PDF运输方程在国内外都有很好的应用。
目前国内外的研究
目前PDF模型中取得的重要进展有速度-组分-湍流频率PDF方法，基于PDF 模型的用于大涡模拟的亚滤波尺度模型（滤波密度函数），基于PDF模型的辐射与湍流的相互作用，基于PDF解决煤粒子和液体喷射问题的模型，求解PDF运输方程的欧拉场方法等[7]。
国内的研究
国内对PDF模型的研究主要中国科学技术大学的陈义良等人，王海峰、陈义良等人将PDF方法和火焰面模型在湍流扩散燃烧中的应用性能作了比较[8]；王海峰、陈义良等人将有限容积法和MonteCarlo法结合起来求解速度-标量-湍流频率的联合PDF方程的混合算法发展到无结构网格中,并以湍流射流火焰为算例，对该混合算法数值误差进行了分析[9]；曹洪军，张会强等人基于大涡模型对湍流非预混火焰的PDF进行了研究，其预测的瞬态温度场显示了合理的湍流火焰形态[10]；吴辉霞,刘玉英等人采用κ-ε湍流模型、概率密度函数(PDF)法及详细化学反应机理模拟湍流燃烧的基础上,研究了辐射换热对湍流射流燃烧过程的影响[11]。
国外的研究
目前国外研究对PDF模型研究主要有康奈尔大学的Pope, S.B.和University
of Texas的Raman, V等人。Pope等人运用LES/PDF方法研究了CH4/H2/N2湍流喷射火焰[12]；Raman, V等人将LES/PDF方法研究天然气喷射扩散湍流火焰中碳粒子与湍流的相互影响[13]；Haworth, D.C等人将PDF模型结合蒙特卡洛法对湍流与辐射相互作用进行研究，模拟与实验相符，并讨论了辐射换热的影响[14]。
PDF模型存在的问题
PDF目前亟需解决的问题有二个方面：小尺度混合模型的问题，超声速流动的问题。小尺度混合模型是求解湍流与化学反应相互作用的关键，目前现有的小尺度混合模型没有考虑化学反应对分子混合的影响。PDF模型在超声速流动中的应用需把PDF模型发展成速度-标量-焓-湍流频率-联合PDF，但是该该模型目前存在以下问题：没有模型可以准确描述压力膨胀项对湍流动能的影响；现有模型不能反应膨胀耗散项在超声速反应流中的重要作用；目前没有合适准确的随机压力模型和随机内能模型[15]。
PDF模型计算成本高的问题往往取决于化学源项。因此，如果在机理简化和化学加速算法方面取得进展会明显降低PDF计算成本。为将PDF模型直接用于实际燃烧和更广的能源问题，有必要集成以下几个因素：适用于多种燃烧方式（包括在接近熄火时的非传统燃烧结构，至今还未研究清楚的化学热力学机理等）的多尺度、多物理场模型，湍流、气相化学、液体喷雾、煤烟、热辐射之间复杂的相互影响以及多组分燃料等。可以预见，为解决这些问题，基于PDF的方法在21世纪将会得到更广阔的应用。
参考文献
[1]Lundgren.T.S. Model Equation for Nonhomogeneous Turbulence [J]. Phys.Fluids1969(12): 485
[2]Pope S.B.The probability approach to the modelling of turbulent reacting flows[J].Combustion and Flame, V olume 27, August–December 1976:299-312
[3]Nguyen, Thong Van, and Pope S.B.Monte Carlo calculations of turbulent diffusion flames[J].Combustion Science and Technology,42, 1984,1:13-45.
[4]陈义良,斯蒂夫·波柏. 湍流速度和耗散率联合几率密度函数的输运方程[J]. 应用数学和力学,1989,01:23-32
[5]PopeS.B., and Y.L.Chen.The velocity‐dissipation probability density function model for turbulent flows[J].Physics of Fluids A: Fluid Dynamics2 (1990): 1437.
[6]Delarue, B.J., and PopeS.B.Calculations of subsonic and supersonicturbulent reacting mixing layers using probability density function methods[J].Physics of Fluids,1998,10: 487.
[7]Haworth, D.C.Progress in probability density function methods for turbulent reacting flows[J]. Progress in Energy and Combustion Science, 2010,36(2): p. 168-259.
[8]王海峰, 陈义良与刘明侯, 湍流扩散燃烧的数值研究-PDF方法和火焰面模型的性能比较[J]. 工程热物理学报, 2005,26: 241-244.
[9]朱旻明, 陈义良, 吕柳艳,王海峰.湍流射流火焰的PDF模拟及误差分析[J].中国科学技术大学学报,2004,34(2):176-182.
[10]曹红军, 张会强,林文漪.基于大涡模拟的湍流非预混燃烧混合分数概率密度函数[J]. 清华大学学报(自然科学版), 2012(07):1018-1022
[11]吴辉霞,刘玉英,薛然然,张荣春,樊未军. 湍流射流火焰热辐射影响的数值研究[J]. 工程热物理学报,2010,11:1961-1964.
[12]WangH.,https://www.360docs.net/doc/f414916743.html,rge eddy simulation/probability density function modeling of a turbulent CH4/H2/N2 jet flame[J].Proceedings of the Combustion Institute, 2011,33 (1):1319-1330
[13]Donde P., RamanV., MuellerM.E., PitschH.LES/PDF based modeling of soot-turbulence interactions in turbulent flames[J].Proceedings of the Combustion Institute, 2013,34 (1):1183-1192.
[14]WangA., ModestM.F., HaworthD.C.,Wang L.Monte Carlo simulation of radiative heat transfer and turbulence interactions in methane/air jet flames[J]. Journal of Quantitative Spectroscopy and Radiative Transfer, 2008,109(2):269-279
[15]范周琴, 孙明波,刘卫东.湍流燃烧的概率密度函数输运方程模型研究[J]. 飞航导弹, 2010(05): 90-95
湍流模型概述
大多数飞行器都是在高Re数下飞行，表面的流态是湍流。为了准确地确定湍流流态下的摩阻、热流，湍流成为一个重要而困难的研究课题。 （一）DNS 目前处理湍流数值计算问题有三种方法，第一种方法即所谓直接数值模拟方法（DNS方法），直接求解湍流运动的N-S方程，得到湍流的瞬时流场，即各种尺度的随机运动，可以获得湍流的全部信息。随着现代计算机的发展和先进的数值方法的研究，DNS方法已经成为解决湍流的一种实际的方法。但由于计算机条件的约束，目前只能限于一些低Re数的简单流动，不能用于工程应用。目前国际上正在做的湍流直接数值模拟还只限于较低的需诺数(Re～200)和非常简单的流动外形，如平板边界层、完全发展的槽道流，以及后台阶流动等。用直接数值模拟方法处理工程中的复杂流动问题，即使是当前最先进的计算机也还差三个量级。 （二）LES 另一种方法称做大涡模拟方法（LES方法）。这是一种折衷的方法，即对湍流脉动部分直接地模拟，将N-S方程在一个小空间域内进行平均（或称之为滤波），以使从流场中去掉小尺度涡，导出大涡所满足的方程。小涡对大涡的影响会出现在大涡方程中，再通过建立模型（亚格子尺度模型）来模拟小涡的影响。由于湍流的大涡结构强烈地依赖于流场的边界形状和边界条件，难以找出普遍的湍流模型来描述具有不同的边界特征的大涡结构，宜做直接模拟。相反地，小尺度涡对边界条件不存在直接依赖关系，而且一般具有各向同性性质。所以亚格子模型具有更大的普适性，比较容易构造，这是它比雷诺平均方法要优越的地方。自从1970年Deardorff第一次给出具有工程意义的LES计算以来，LES方法已经成为计算湍流的最强有力的工具之一，应用的方向也在逐步扩展，但是仍然受计算机条件等的限制，使之成为解决大量工程问题的成熟方法仍有很长的路要走。 （三）RANS 目前能够用于工程计算的方法就是模式理论。所谓湍流模式理论，就是依据湍流的理论知识、实验数据或直接数值模拟结果，对Reynolds应力做出各种假设，即假设各种经验的和半经验的本构关系，从而使湍流的平均Reynolds方程封闭。随着计算流体力学的发展，湍流模式理论也有了很大的进步，有了非常丰硕的成果。从对模式处理的出发点不同，可以将湍流模式理论分类成两大类：一类称为二阶矩封闭模式，另一类称涡粘性封闭模式。 （1）雷诺应力模式 所谓二阶矩封闭模式，是从Reynolds应力满足的方程出发，将方程右端未知的项（生成项，扩散项，耗散项等）用平均流动的物理量和湍流的特征尺度表示出来。典型的平均流动的变量是平均速度和平均温度的空间导数。这种模式理论，由于保留了Reynolds应力所满足的方程，如果模拟的好，可以较好地反映Reynolds应力随空间和时间的变化规律，因而可以较好地反映湍流运动规律。因此，二阶矩模式是一种较高级的模式，但是，由于保留了Reynolds应力的方程，加上平均运动的方程整个方程组总计15个方程，是一个庞大的方程组，应用这样一个庞大的方程组来解决实际工程问题，计算量很大，这就极大地限制了二阶矩模式在工程问题中的应用。 （2）涡粘性模式
煤粉燃烧模拟——湍流破碎模型.
煤粉燃烧——湍流破碎模型（EBU） 简介 该帮助文件主要介绍煤粉燃烧模型的设置和求解，采用湍流破碎模型（EBU）。 EBU燃烧模型，也称涡团破碎模型，假设化学反应的平均速度与化学动力学无关，而只取决于低温的反应物和高温的燃烧产物之间的湍流混合作用。 主要包括： 1）煤粉燃烧模型的建立和求解 2）湍流破碎模型（EBU）的应用 3）选择合适的求解参数 4）计算结果的后处理 问题描述 3D模型的剖面图如图1所示。左侧为两个环形入口，右侧为一个圆形出口。由于模型的对称性，取系统的1/4进行建模。煤粉和携带空气（一次风）从内环进入燃烧室，二次风从外环进入燃烧室，发生燃烧反应，产物从压力出口流出。模型建立和求解 Step1：网格 1、读取mesh文件 2、检查网格：grid check 3、显示网格：display grid a）从列表框中选取所有surfaces b）点击display，并关闭grid display 面板 Step2：模型 1、选择k-ε湍流模型 Define models viscous 2、启动能量方程 Define models energy 3、启动物质输运方程 Define models species transport & reaction
a）选择model列表中的species transport b）选择reaction列表中的volumetric c）选择mixture material 列表中的coal-hv-volatiles-air d）选择turbulence-chemistry interaction（湍流与化学反应的作用）列表中的eddy dissipation（涡流耗散） e）点击ok，关闭species model面板 4、启动Discrete ordinates模型（DO离散坐标系） a）从model列表中选择discrete ordinates b）设置flow iterations per radiation iteration（流动和辐射迭代次数）为1 c）设置angular discretization （角离散化）中的Theta Divisions and Phi Divisions为4 d）设置angular discretization （角离散化）中的Theta Pixels and Phi Pixels 为3 e）点击ok，关闭radiation model 5、启动discrete phase（离散相）模型 a）M ax number of steps（最大步数）40000 b）启动specify length scale（步长），设为0.0025 c）点击ok，关闭discrete phase面板 Step3：injections（喷口） 1、v-1入口截面设9个喷口：define injections a）点击creat，新建喷口 c）9个喷口通用性质，见表 c）点击turbulent dispersion，并启动discrete random walk model
计算流体力学
1、数值的耗散与频散： 在数值解中出现的振幅衰减波长加宽的现象叫数值耗散，与高阶偶次空间偏导数有关；在数值解中出现解得主波后有一系列频及传播速度不等的尾波的现象叫数值频散，与高阶奇次偏导数有关。 2、湍流模型理论：湍流模式理论或简称湍流模型，就是以雷诺平均运动方程与脉动运动方程为基础，依靠理论与经验的结合，引进一系列模型假设，而建立起得一组 描写湍流平均量的封闭方程组。 3、修正的偏微分方程：与差分方程相等价的微分方程称之为修正的微分方程。 4、自适应网格:为了计算具有高雷诺数的流场，必须将流场内的网格加密，但是实际计算中并不需要对全流场的网格所有部分同样加密，只需在某些部分，如物面附近、尾流区等得网格加密即可。因此需要事先估计一些变化较快的区域，但这种估计又是是正确的。有时则不正确。特别是不定常流动，流动过程本身就是变化的，所以需要不断的调整网格的位置和疏密，这样就产生了自适应网格。 5、CFL 条件：定义t C x μ ?=? ，不等式1C ≤ 称为CFL 条件，此条件一般应用于双曲线偏微分方程的显式格式。物理意义：即在时间步长内，波的位移应小于空间步长。 数学意义：差分方程解的依赖区域包含微分方程解得依赖区域。 1、简答CFD 方法求解流动问题的基本步骤 答：①确定流动模型；②计算区域离散化；③用离散节点变量代替场；④将控制方程中偏导数进行离散，得到线性方程组；⑤边界条件和初值条件离散化；⑥离散的线性方程组求解，得到离散值；⑦计算结果数据处理。 2、简述离散偏微分方程的三个原则及LAX 定理 三原则；相容性、稳定性、收敛性。 LAX 定理：对于一个选定的线性偏微分方程的初值问题，对应的差分方法是相容的，则差分方程解得收敛性和稳定性事等价的或者说稳定性是收敛性的充要条件。 3、简述差分格构造的基本规律，并应用规律方程 0t x μμ λ??+=?? 利用网格点（） （）（）构造方程的差分格式，并验证其离散格式的精度等级。 答：构造的基本规律 ：①为保证均匀流场，差分的分子各项系数之和为零 ②分母向量级与微分的阶数一致 ③构造差分级指明针对哪点构造 ④差分格式的精度 由网格点（）（）（）规律方程（）构造得 1 11 1 0n n n n j j j j x x μμμμλ +++---+=?? 令112j j x k k x μμμ-+= ? 用泰勒公 式展开的23 126j j x xx x x x x μμμμμ-??=-?+- 所以12101k k k +=?? -=? 得12 11k k =-??=? 所以1j j x x μμμ--+=? 所以具有一阶精度 4、简要概括流动的数值计算对网格的基本要求 答：①计算域边界上的网格节点都应在边界上 ②物理域上的特点与计算域上的节点要求一一对应 ③网格应尽量尺寸匀称，相邻网格长度比应小于2 ④物理域网格夹角不宜太小（≥45°） ⑤流动参数梯度大的地方网格要加密，否则稀疏。 5、简述人工压缩方法（时间相关法）的基本思想 答：用非定常流动方程来求解定常流动问题，用其稳态求解定常流动的解，将不可压缩的粘性流动的连续方程，添加到可压缩项。则与动量方程构成定常粘性流动时间相关方程，可把非定常流动的稳态解作为非定常流动的解。
湍流模型发展综述
湍流模型发展综述 摘要：在概述了湍流问题的基础上，本文简要介绍了湍流的四种模型，对湍流模型在不同情况下的模拟能力进行了对比，最后简述了湍流模型的发展方向。 关键词：湍流模型；Navier-Stokes方程组；J-K模型 Abstract：On the basis of introducing the problems of turbulence, this paper briefly analyzed four kinds of turbulence models and compared their ability of simulation in different situations. At last, the paper expounded the development direction of the turbulence model. Key words：Turbulence model; Navier-Stokes equations; J-K model 一、引言 湍流又称紊流,是自然界中常见的一种很不规则的流动现象。当粘性阻尼无法消除惯性的影响时,自然界中的绝大部分流动都是湍流。 湍流运动的实验研究表明，虽然湍流结构十分复杂，但它仍然遵循连续介质的一般动力学规律，湍流流动的各物理量的瞬时值也应该服从一般的N-S方程。对粘性流体服从的N-S方程进行时均化，就可以得到雷诺平均方程。与定常的N-S方程相比，不同之处是在该式右边多了九项与脉动量有关的项，这脉动量的乘积的平均值与密度的乘积是湍流流动中的一种应力，称为湍流应力或雷诺应力。其中，法向雷诺应力和切向雷诺应力各有三个。 湍流问题就是在给定的边界条件下解雷诺方程。由于雷诺平均方程中未知数个数远多于方程个数而出现了方程不封闭的问题，这就需要依据各种半经验理论提出相应的补充方程式，即各种湍流模型。一般按照所用湍流量偏微分方程的物理含义或者数量进行区分，分别称为梅罗尔—赫林方法和雷诺方法。而后者又将湍流模型分成四类。（1）零方程模型；（2）一方程模型；（3）二方程模型；（4）应力方程模型。下面就对这些模型进行简单的描述。 二、湍流模型简介 1、零方程模型 最初的湍流模型只考虑了一阶湍流计算统计量的动力学微分方程，即平均方程，没有引进高阶统计量的微分方程，因而称之为一阶封闭模式或零方程模型。零方程模型又称为代数模型，代数模型又可以分成以下几种模型：（1）Cebeci —Smith 模型，（2）Baldwin—Lomax 模型，（3）Johnson—King 模型。 其中，B-L与C-S模型的不同之处在于外层湍流粘性系数取法不同。后者适用于湍流边界层，而前者则可用于 N-S方程的计算。此两模型已在工程计算中
第15章 预混燃烧模拟
第十五章预混燃烧模拟FLUENT有一个预混湍流燃烧模型，基于反应过程参数方法。有关这一模型的内容按以下节次给出： ●15.1 概述和限制 ●15.2 预混燃烧模型 ●15.3 使用预混燃烧模型 15.1 概述和限制 15.1.1 概述 在预混燃烧中，燃料和氧化剂在点火之前进行分子级别的混合。火焰前锋传入未燃烧的反应物产生燃烧。预混燃烧的例子有吸气式内燃机，稀薄燃气轮机的燃烧器，气体泄露爆炸。 预混燃烧比非预混燃烧更难以模拟。原因在于（亚音速）预混燃烧通常做为薄层火焰产生，并被湍流拉伸和扭曲。火焰传播的整体速率受层流火焰速度和湍流涡旋控制。层流火焰速度由物质和热量逆流扩散到反应物并燃烧的速率决定。为得到层流火焰速度，需要确定内部火焰结构以及详细的化学动力学和分子扩散过程。由于实际的层流火焰厚度只有微米量级或更小，求解所需要的开销是不可承受的。 湍流的影响是使传播中的层流火焰层皱折、拉伸，增加了薄层的面积，并因此提高了火焰速度。大的湍流涡使火焰层皱折，而小的湍流涡，如果它们比层流火焰的厚度还小，将会穿过火焰层并改变层流火焰结构。 与之相比，非预混燃烧可以极大地简化为一个混合问题（例如，14.1节中介绍的混合物组分方法）。预混燃烧模拟的要点在于捕获湍流火焰速度，它受层流火焰速度和湍流的影响。 在预混火焰中，燃料和氧化剂在进入燃烧设备之前已经紧密混合。反应在燃烧区发生，这一区域将未燃烧的反应物和燃烧产物隔开。部分预混火焰具有预混和扩散火焰两方面的性质。它们发生在有额外的氧化剂或燃料气流进入预混系统，或是当扩散火焰离开燃烧器以在燃烧前产生某些预混的情况。 预混和部分预混火焰FLUENT的有限速率公式(见13章)模拟。还可以参阅16章了解更多有关FLUENT部分预混燃烧模型方面的信息。如果火焰是完全预混合的，则只有一股具有单一混合比的气流进入燃烧器，可以使用预混燃烧模型。 15.1.2 限制 在使用预混燃烧模型时有以下限制： ●必须使用非耦合求解器。预混燃烧模型在两种耦合求解器中都不能得到。 ●预混燃烧模型只对湍流、亚音速模型有效。这一类型的火焰成为爆燃。在爆炸中， 可燃混合物被冲击波后面的热量点燃，这一类型的燃烧可以使用非耦合和耦合求解 器用有限速率模型模拟。有关限速率模型见13章。 ●预混燃烧模型不能和污染物（如碳烟和NOx）模型一起使用。但完全预混系统可以 用部分预混模型（见16章）模拟。 ●不能用预混燃烧模型模拟反应的离散相粒子。只有惰性粒子可以使用预混燃烧模 型。 15.2 预混燃烧理论 湍流预混燃烧模型基于Zimont等人的工作[275，276，278]，涉及求解一个关于反应过
湍流模型
我们知道，描述流体运动（层流）的流体力学基本方程组是封闭的，而描述湍流运动的方程组由于采用了某种平均（时间平均或网格平均等）而不封闭，须对方程组中出现的新未知量采用模型而使其封闭，这就是CF D中的湍流模型。湍流模型的主要作用是将新未知量和平均速度梯度联系起来。目前，工程应用中湍流的数值模拟主要分三大类：直接数值模拟（D NS）；基于雷诺平均N-S方程组（RANS）的模型和大涡模拟（LES）。DNS是直接数值求解N-S方程组，不需要任何湍流模型，是目前最精确的方法。其优点在于可以得出流场内任何物理量（如速度和压力）的时间和空间演变过程，旋涡的运动学和动力学问题等。由于直接求解N-S方程，其应用也受到诸多方面的限制。第一：计算域形状比较简单，边界条件比较单一；第二：计算量大。影响计算量的因素有三个：网格数量、流场的时间积分长度（与计算时间长度有关）和最小旋涡的时间积分长度（与时间步长有关），其中网格数量是重要因素。为了得到湍流问题足够精确的解，要求能够数值求解所有旋涡的运动，因此要求网格的尺度和最小旋涡的尺度相当，即使采用子域技术，其网格规模也是巨大的。为了求解各个尺度旋涡的运动，要求每个方向上网格节点的数量与Re3/4成比例，考虑一个三维问题，网格节点的数量与Re9/4成比例。目前，DNS能够求解Re(104)的范围。 基于RANS的湍流模型采用雷诺平均的概念，将物理量区分为平均量和脉动量，将脉动量对平均量的影响用模型表示出来。目前，基于RANS方程已经发展了许多模型，几乎能对所有雷诺数范围的工程问题求解，并得出一些有用的结果。其缺点在于：第一：不同的模型解决不同类型的问题，
湍流燃烧模型-PDF
PDF 模型 概率密度函数PDF方法以随机的观点来对待湍流问题，对解决湍流化学反应流的问题具有很强的优势。在湍流燃烧中存在一些非输运量( 如反应速率, 密度, 温度及气相体积分数等) 的湍流封闭问题。尽管这些量没有输运方程, 但它们常常是输运变量的已知函数。平均或者过滤高度非线性的化学反应源项会引起方程的封闭问题。因此，用PDF的方法来解决这些非输运量的湍流封闭问题显然是一个既简单又直接的途径。 PDF方法是一种较为流行的湍流燃烧模型，能够较为精确的模拟任何详细的化学动力学过程, 适用于预混、非预混和部分预混的任何燃烧问题。目前, 确定输运变量脉动概率密度函数的方法有输运方程和简化假定两种, 分别称之为输运方程的PDF和简化的PDF。前者建立输运变量脉动的概率密度输运方程，通过求解该方程来获得输运变量脉动的概率分布。后者假定输运变量脉动的概率密度函数的具体形式, 通过确定其中的一些待定参数来获得输运变量脉动的概率分布。湍流燃烧中, 后者应用最为普遍和广泛。在简化的PDF 中, 输运变量脉动的概率密度函数常常采用双 D 分布、截尾高斯分布和B 函数分布等形式。 PDF在理论上可以精确考虑任意详细的化学反应机理，但是其具体求解时需借助其它的模型和算法，而且计算量相对较大。PDF的方程是由N-S方程推导而来，其中的化学反应源项是封闭的，但压力脉动梯度项以及分子粘性和分子扩散引起的PDF的分子输运项是不封闭的，需要引入模型加以封闭。例如，在速度- 标量-湍流频率PDF中，必须采用小尺度混合模型、随机速度模型和湍流频率模型加以封闭。 模化后的输运方程难以用有限容积、有限差分和有限元等方法来求解，比较可行的一种方法是蒙特卡洛(MonteCarlo)方法,在该方法中输运方程被转化为拉格朗日(Lagrangian)方程,流体由大量遵循Lagrang ian方程的随机粒子的系统来描述, 最后对粒子作统计平均得到流场物理量和各阶统计矩。另有与有限容积法相结合的蒙特卡洛法。 PDF 模型的发展 1969年Lungdren首先推导、计算了速度的联合PDF运输方程，避免了对梯度扩散模型进行模拟，对很简单的流动过程得到了简析解[1]。
计算流体力学实例
汽车外部气体流动模拟 振动和噪声控制研究所 1.模型概述 在汽车外部建立一个较大的长方体几何空间，长度约为30m，宽度和高度约为5m，在空间内部挖出汽车形状的空腔，汽车尺寸参照本田CRV为4550mm*1820mm*1685mm。由于汽车向前开进，气体从车头流向车尾，因此将汽车前方空间设为气体入口，后方空间设为气体出口，模拟气体在车外的流动。另外为了节省计算成本将整个模型按1:100的比例缩小，考虑到模型和流体均是对称的，因此仅画出几何模型的一半区域，建立对称面以考虑生成包含理想气体的流体域。在Catia中建立的模型如图1.1所示。 图1.1几何模型 2.利用ICEM CFD进行网格划分 a)导入有Catia生成的stp格式的模型； b)模型修复，删除多余的点、线、面，允许公差设为0.1； c)生成体，由于本模型仅为流体区域，因此将全部区域划分为一个体，选取方法可以 使用整体模型选取； d)为了后面的设置边界方便，因此将具有相同特性的面设为一个part,共设置了in， out，FreeWalls，Symmetry和Body； e)网格划分，设置Max element=2，共划分了1333817个单元，有225390个节点； f)网格输出，设置求解器为ANSYS CFX，输出cfx5文件。 3.利用ANSYS CFX求解 a)生成域，物质选定Air Ideal Gas，参考压强设为1atm，浮力选项为无浮力模型，
域运动选项为静止，网格变形为无；流体模型设定中的热量传输设定为Isothermal，流体温度设定为288k，湍流模型设定为Shear Stress Transport模型，壁面函数 选择Automatic。 b)入口边界设定，类型为Inlet，位置选定在in，质量与栋梁选定Normal Speed，设 定为15m/s，湍流模型设定类型为Intensity and Length Scale=0.05，Eddy Len.Scale=0.1m。 c)出口边界设定，边界类型为Outlet，位置选out。质量与动量选项为Static Pressure，相对压强为0pa。 d)壁面边界设定，边界类型为Wall，位置选在FreeWalls。壁面边界详细信息中指定 WallInfluence On Flow为Free Slip。 e)对称边界设定，边界类型为Symmetry，位置选在Symmetry。 f)汽车外壁面设定，边界类型为Wall，位置设在Body，壁面详细信息选项中指定Wall Influence On Flow为No Slip，即汽车壁面为无滑移壁面。 g)初始条件设定，初始速度分量设为U方向为15m/s，其他两个方向的速度为零。 h)求解设置，残差类型选为RMS，残差目标设定为1e-5，当求解达到此目标时，求解 自动终止。求解之前的模型如图3.1所示。 图3.1求解之前的模型 4.结果后处理 从图4.1中可以看出计算收敛。
fluent湍流模型
第十章湍流模型 本章主要介绍Fluent所使用的各种湍流模型及使用方法。 各小节的具体内容是： 10．1 简介 10．2 选择湍流模型 10．3 Spalart-Allmaras 模型 10．4 标准、RNG和k-e相关模型 10．5 标准和SST k-ω模型 10．6 雷诺兹压力模型 10．7 大型艾迪仿真模型 10．8 边界层湍流的近壁处理 10．9 湍流仿真模型的网格划分 10．10 湍流模型的问题提出 10．11 湍流模型问题的解决方法 10．12 湍流模型的后处理 10．1 简介 湍流出现在速度变动的地方。这种波动使得流体介质之间相互交换动量、能量和浓度变化，而且引起了数量的波动。由于这种波动是小尺度且是高频率的，所以在实际工程计算中直接模拟的话对计算机的要求会很高。实际上瞬时控制方程可能在时间上、空间上是均匀的，或者可以人为的改变尺度，这样修改后的方程耗费较少的计算机。但是，修改后的方程可能包含有我们所不知的变量，湍流模型需要用已知变量来确定这些变量。 FLUENT 提供了以下湍流模型： ·Spalart-Allmaras 模型 ·k-e 模型 －标准k-e 模型 －Renormalization-group (RNG) k-e模型 －带旋流修正k-e模型 ·k-ω模型 －标准k-ω模型 －压力修正k-ω模型 －雷诺兹压力模型 －大漩涡模拟模型 10．2 选择一个湍流模型 不幸的是没有一个湍流模型对于所有的问题是通用的。选择模型时主要依靠以下几点：流体是否可压、建立特殊的可行的问题、精度的要求、计算机的能力、时间的限制。为了选择最好的模型，你需要了解不同条件的适用范围和限制 这一章的目的是给出在FLUENT中湍流模型的总的情况。我们将讨论单个模型对cpu 和内存的要求。同时陈述一下一种模型对那些特定问题最适用，给出一般的指导方针以便对于你需要的给出湍流模型。 10．2．1 雷诺平均逼近vs LES 在复杂形体的高雷诺数湍流中要求得精确的N-S方程的有关时间的解在近期内不太可能实现。两种可选择的方法用于把N-S方程不直接用于小尺度的模拟：雷诺平均和过滤。
湍流燃烧数值模拟研究
湍流燃烧及其数值模拟研究 1. 湍流燃烧 1.1 湍流燃烧基本概念 当流动雷诺数数较小时，由于流体粘性的作用，流体呈层流流态。当流动的特征雷诺数超过相应的临界值，流动从层流转捩到湍流。湍流燃烧是指湍流流动 中可燃气的燃烧，在能源、动力、航空和航天等工程领域,经常遇到的实际燃烧过程几乎全部都是湍流燃烧过程。湍流燃烧实质是湍流，化学反应和传热传质等过程相耦合的结果。湍流对燃烧的影响与湍流强度和湍流涡旋尺度有关。小尺度湍流通过湍流扩散使火焰区内的输运效应增加，从而使化学反应速率增加。但气流脉动不会火焰面产生皱褶，只能把火焰变成波纹状。大尺度湍流对火焰内部结构没有影响，但使火焰阵面出现皱褶，增加其燃烧面积，造成火焰表现传播速度增加。当湍流强度及湍流尺度均较大时，火焰前沿不再连续而分裂成四分五裂。 燃烧对湍流的影响主要表现在燃烧释放的热流流团膨胀，影响气体的密度和运动速度，从而影响当地的涡旋，湍流强度和湍流结构。 1.2 湍流燃烧分类湍流燃烧按其燃料和氧化剂的初始混合状态可以分类为：湍流非预混燃烧、预混燃烧和部分预混燃烧。在湍流非预混燃烧燃料和氧化剂事先是分离的，燃料和氧化剂一边混合一边燃烧，燃烧速率主要受湍流混合过程控制，而在湍流预混燃烧中，燃料和氧化剂在进入核心燃烧区以前已经充分混合，化学反应的速率由火焰前缘从炽热的燃烧区向冷态无反应区的传播所控制。上面两种燃烧方式是湍流燃烧的两个极限情形，很多情况下两种燃烧模式是并存的，称为部分预混燃烧。部分预混燃烧可出现在下列情形中叫：(1) 在一个完全以非预混燃烧为配置的燃烧装置发牛了局部熄火；(2) 当预混火焰前缘穿过非均匀的混气时；(3) 射流非预混火焰发生抬举，其根部是一。个典型的部分预混火焰。这三种部分预混燃烧情形涉及了经常受到关注的燃烧研究话题如局部熄火、火焰稳定等，它们对研究湍流燃烧过程的机理有很大意义。 在湍流燃烧中,湍流流动过程和化学反应过程有强烈的相互关联和相互影响. 湍流通过强化混合而影响着时平均化学反应速率,同时化学反应放热过程又影响着湍流,如何定量地来描述和确定这种相互作用是湍流燃烧研究的一个重要内容. 湍流是非常复杂的，它包括湍流问题，湍流与燃烧的相互作用，流动参数与
第二章计算流体力学的基本知识
第二章计算流体力学的基本知识 流体流动现象大量存在于自然界及多种工程领域中，所有这些工程都受质量守恒、动量守恒和能量守恒等基本物理定律的支配。这章将首先介绍流体动力学的发展和流体力学中几个重要守恒定律及其数学表达式，最后介绍几种常用的商业软件。 2.1 计算流体力学简介 2.1.1计算流体力学的发展 流体力学的基本方程组非常复杂，在考虑粘性作用时更是如此，如果不靠计算机，就只能对比较简单的情形或简化后的欧拉方程或N-S方程进行计算。20世纪30～40年代，对于复杂而又特别重要的流体力学问题，曾组织过人力用几个月甚至几年的时间做数值计算，比如圆锥做超声速飞行时周围的无粘流场就从1943年一直算到1947年。 数学的发展，计算机的不断进步，以及流体力学各种计算方法的发明，使许多原来无法用理论分析求解的复杂流体力学问题有了求得数值解的可能性，这又促进了流体力学计算方法的发展，并形成了"计算流体力学"。 从20世纪60年代起，在飞行器和其他涉及流体运动的课题中，经常采用电子计算机做数值模拟，这可以和物理实验相辅相成。数值模拟和实验模拟相互配合，使科学技术的研究和工程设计的速度加快，并节省开支。数值计算方法最近发展很快，其重要性与日俱增。 自然界存在着大量复杂的流动现象，随着人类认识的深入，人们开始利用流动规律来改造自然界。最典型的例子是人类利用空气对运动中的机翼产生升力的机理发明了飞机。航空技术的发展强烈推动了流体力学的迅速发展。 流体运动的规律由一组控制方程描述。计算机没有发明前，流体力学家们在对方程经过大量简化后能够得到一些线形问题解析解。但实际的流动问题大都是复杂的强非线形问题，无法求得精确的解析解。计算机的出现以及计算技术的迅速发展使人们直接求解控制方程组的梦想逐步得到实现，从而催生了计算流体力
第章预混燃烧模拟
第十五章预混燃烧模拟 FLUENT 有一个预混湍流燃烧模型，基于反应过程参数方法。有关这一模型的内容按以下节次给出： 15.1概述和限制 15.2 预混燃烧模型 15.3 使用预混燃烧模型 15.1 概述和限制 15.1.1概述在预混燃烧中，燃料和氧化剂在点火之前进行分子级别的混合。火焰前锋传入未燃烧的反应物产生燃烧。预混燃烧的例子有吸气式内燃机，稀薄燃气轮机的燃烧器，气体泄露爆炸。 预混燃烧比非预混燃烧更难以模拟。原因在于（亚音速）预混燃烧通常做为薄层火焰产生，并被湍流拉伸和扭曲。火焰传播的整体速率受层流火焰速度和湍流涡旋控制。层流火焰速度由物质和热量逆流扩散到反应物并燃烧的速率决定。为得到层流火焰速度，需要确定内部火焰结构以及详细的化学动力学和分子扩散过程。由于实际的层流火焰厚度只有微M 量级或更小，求解所需要的开销是不可承受的。 湍流的影响是使传播中的层流火焰层皱折、拉伸，增加了薄层的面积，并因此提高了火焰速度。大的湍流涡使火焰层皱折，而小的湍流涡，如果它们比层流火焰的厚度还小，将会穿过火焰层并改变层流火焰结构。 与之相比，非预混燃烧可以极大地简化为一个混合问题（例如，14.1节中介绍的混 合物组分方法）。预混燃烧模拟的要点在于捕获湍流火焰速度，它受层流火焰速度和湍流的影响。 在预混火焰中，燃料和氧化剂在进入燃烧设备之前已经紧密混合。反应在燃烧区发生，这一区域将未燃烧的反应物和燃烧产物隔开。部分预混火焰具有预混和扩散火焰两方面的性质。它们发生在有额外的氧化剂或燃料气流进入预混系统，或是当扩散火焰离开燃烧器以在燃烧前产生某些预混的情况。 预混和部分预混火焰FLUENT的有限速率公式（见13章濮拟。还可以参阅16章了解更多有关FLUENT部分预混燃烧模型方面的信息。如果火焰是完全预混合的，则只有一股具有单一混合比的气流进入燃烧器，可以使用预混燃烧模型。 15.1.2限制 在使用预混燃烧模型时有以下限制：必须使用非耦合求解器。预混燃烧模型在两种耦合求解器中都不能得到。预混燃烧模型只对湍流、亚音速模型有效。这一类型的火焰成为爆燃。在爆炸中，可燃混合物被冲击波后面的热量点燃，这一类型的燃烧可以使用非耦合和耦合求解器用有限速率模型模拟。有关限速率模型见13章。预混燃烧模型不能和污染物（如碳烟和 NOx ）模型一起使用。但完全预混系统可以用部分预混模型（见16 章）模拟。 不能用预混燃烧模型模拟反应的离散相粒子。只有惰性粒子可以使用预混燃烧模型。 15.2预混燃烧理论 湍流预混燃烧模型基于Zimont 等人的工作[275，276，278]，涉及求解一个关于反应过程变量的输运方程。这一方程的封闭基于湍流火焰速度的定义。 15.2.1 火焰前锋的传播 在许多工业预混系统中，燃烧发生在一个非常薄的火焰层中。当火焰前锋移动时，未燃的反应物燃烧，变为燃烧产物。因此预混燃烧模型用火焰层将反应的流场分为已燃物区和未燃物区。反应的传播等同于火焰前锋的传播。 火焰前锋传播的模拟通过求借一个关于标量c的输送方程，c为(Favre平均)反应进 程变量。
四种湍流模型介绍
由于航发燃烧室中的流动特性极其复杂，要想提高数值计算的预测能力,必须要慎重选择湍流模型。用四种不同的湍流模型对带双径向旋流杯的下游流场进行数值模拟,将计算结果与实验结果作对比，比较各湍流模型的原理和物理基础，优劣，并分析流场速度分布和回流区特性。 涉及的湍流模型： 标准k-ε湍流模型（SＫE) １标准k-ε湍流模型有较高的稳定性,经济性和计算精度，应用广泛,适合高雷诺数湍流，但不适合旋流等各向异性较强的流动。 2简单的湍流模型是两个方程的模型，需要解两个变量，即速度和长度。在fｌｕｅｎt中，标准 k-ε湍流模型自从被Lauｎdeｒand Spalｄing 提出之后，就变成流场计算中的主要工具。其在工业上被普遍应用,其计算收敛性和准确性都非常符合工程计算的要求。 3但其也有某些限制,如ε方程包含不能在壁面计算的项，因此必须使用壁面函数。另外，其预测强分离流,包含大曲率的流动和强压力梯度流动的结果较弱。 它是个半经验的公式,是从实验现象中总结出来的。 动能输运方程是通过精确的方程推导得到,耗散率方程是通过物理推理，数学上模拟相似原型方程得到的。 应用范围：该模型假设流动为完全湍流，分子粘性的影响可以忽略,此标准κ－ε模型只适合完全湍流的流动过程模拟。 可实现的ｋ-ε模型是才出现的,比起标准k-ε模型来有两个主要的不同点：·可实现的k－ε模型为湍流粘性增加了一个公式。 ·为耗散率增加了新的传输方程,这个方程来源于一个为层流速度波动而作的精确方程。 术语“rｅalizablｅ”,意味着模型要确保在雷诺压力中要有数学约束，湍流的连续性。 应用范围： 可实现的k-ε模型直接的好处是对于平板和圆柱射流的发散比率的更精确的预测。而且它对于旋转流动、强逆压梯度的边界层流动、流动分离和二次流有很好的表现。 可实现的k－ε模型和RＮG k－ε模型都显现出比标准k－ε模型在强流线弯曲、漩涡和旋转有更好的表现。由于带旋流修正的k-ε模型是新出现的模型，所以还没有确凿的证据表明它比ＲＮＧｋ－ε模型有更好的表现。但是最初的研究表明可实现的k-ε模型在所有k-ε模型中流动分离和复杂二次流有很好的作用。 该模型适合的流动类型比较广泛,包括有旋均匀剪切流,自由流(射流和混合层）,腔道流动和边界层流动。对以上流动过程模拟结果都比标准k-ε模型的结果好，特别是可再现k－ε模型对圆口射流和平板射流模拟中,能给出较好的射流扩张。
fluent中燃烧模型分类
FLUENT燃烧模型 化学反应 模拟方法 方法描述计算反应的选择 有限速率模型需要求解组分质量分数的 输运方程，化学反应机理 由用户自己定义。反应速 率在组分输运方程中作为 源项，并由阿累尼乌斯公 式计算。应用范围最广泛。 应用：模拟化学组分混合、 输运和反应的问题；壁面 或粒子表面反应问题 层流有限 速率模型 使用Arrhenius公式计算化学源项，忽略湍流脉动的影响。对于化学动力学控制的燃烧（如层流燃 烧），或化学反应相对缓慢的湍流燃烧是准确的。但对一般湍流火焰中Arrhenius化学动力学的高度 非线性一般不精确；对于化学反应相对缓慢、湍流脉动较小的燃烧（如超音速火焰）可能可以接受。 漩涡破碎模型 Eddy Dissipation 大部分燃料快速燃烧，整体反应速率由湍流混合控制。复杂且常是未知的化学反应动力学速率可以 完全的被忽略掉。化学反应速率由大尺度涡混合时间尺度k/ε控制。只要k/ε（湍流）出现，燃烧 即可进行，不需要点火源来启动燃烧。（缺点：未能考虑分子输运和化学动力学因素的影响） 适用条件：高雷诺数湍流预混燃烧过程。 EBU-Arrehenius 模型 EDC模型 假定化学反应都发生在小涡中（精细涡），反应时间由小涡生存时间和化学反应本身需要的时间共 同控制。EDC模型能够在湍流反应中考虑详细的化学反应机理。但是他们的数值积分计算开销很大。 使用条件：只有在快速化学反应假定无效的情况下才能使用这一模型（如快速熄灭火焰中缓慢的 CO烧尽、选择性非催化还原中的NO转化问题）。 非预混燃烧模型不求解每个组分的质量分数输运方程，求解混合分数输运方程和一个或两个守恒标量的方程，然后从预测的混合分数公布推导出每一个组分的浓度。通过概率密度函数或PDF来考虑湍流的影响。 应用：主要用于模拟湍流扩散火焰的反应系统。这个系统要求接近化学平衡，氧化物和燃料以两个或者三个进口进入计算域。 预混燃烧模型主要用于单一、完全预先混合好的燃烧系统。反应物和燃烧产物被火焰前沿分开。求解出反应发展变量来预测前沿的位置。湍流的影响通过湍流火焰速度计算。 部分预混燃烧模型描述非预混燃烧完全预混燃烧相结合的系统。结合混合分数方程和反应物发展变量来分别确定组分浓度和火焰前沿位置。适用于计算域内具有变化等值比率的预混火焰情况。通过求解混合分数方程和反应过程参数来确定火焰峰面的位置。 PDF输运方程模型结合CHEMKIN可以考虑详细的化学反应机理，高度的非线性化学反应项是精确模拟，无须封闭模型，可以合理的模拟湍流和详细化学反应动力学之间的相互作用，是模拟湍流燃烧的精确模拟方法。但计算量特别大。 优点：可以计算中间组分；考虑分裂影响；考虑湍流-化学反应之间的作用；无需求解组分输运方程 缺点：系统须满足（靠近）局部平衡；不能用于可压缩或非湍流流动；不能用于预混燃烧。
湍流模型介绍
湍流模型介绍 因为湍流现象是高度复杂的，所以至今还没有一种方法能够全面、准确地对所有流动问题中的湍流现象进行模拟。在涉及湍流的计算中，都要对湍流模型的模拟能力以及计算所需系统资源进行综合考虑后，再选择合适的湍流模型进行模拟。FLUENT 中采用的湍流模拟方法 包括Spalart-Allmaras模型、standard（标准）k ?ε模型、RNG（重整化群）k ?ε模型、Realizable（现实）k ?ε模型、v2 ?f 模型、RSM（Reynolds Stress Model，雷诺应力模型）模型和LES（Large Eddy Simulation，大涡模拟）方法。 7.2.1 雷诺平均与大涡模拟的对比 因为直接求解NS 方程非常困难，所以通常用两种办法对湍流进行模拟，即对NS 方程进行雷诺平均和滤波处理。这两种方法都会增加新的未知量，因此需要相应增加控制方程的数量，以便保证未知数的数量与方程数量相同，达到封闭方程组的目的。雷诺平均NS 方程是流场平均变量的控制方程，其相关的模拟理论被称为湍流模式理论。湍流模式理论假定湍流中的流场变量由一个时均量和一个脉动量组成，以此观点处理NS 方程可以得出雷诺平均NS 方程（简称RNS 方程）。在引入Boussinesq 假设，即认为湍流雷诺应力与应变成正比之后，湍流计算就归结为对雷诺应力与应变之间的比例系数（即湍流粘性系数）的计算。根据计算中使用的变量数目和方程数目的不同，湍流模式理论中所包含的湍流模型又被分为二方程模型、一方程模型和零方程模型（代数模型）等大类。 FLUENT 中使用的三种k ?ε模型、Spalart-Allmaras 模型、k ?ω模型及雷诺应力模型RSM）等都属于湍流模式理论。 大涡模拟（LES）方法是通过滤波处理计算湍流的，其主要思想是大涡结构（又称拟 序结构）受流场影响较大，小涡则可以认为是各向同性的，因而可以将大涡计算与小涡计算分开处理，并用统一的模型计算小涡。在这个思想下，大涡模拟通过滤波处理，首先将小于某个尺度的旋涡从流场中过滤掉，只计算大涡，然后通过求解附加方程得到小涡的解。过滤尺度一般就取为网格尺度。显然这种方法比直接求解NS 方程的DNS 方程效率更高，消耗系统资源更少，但却比湍流模式方法更精确。尤其应该注意的是，湍流模式理论无法准确模拟大涡结构，因此在需要模拟大涡结构时，只能采用LES 方法1。 尽管大涡模拟理论比湍流模式理论更精确，但是因为大涡模拟需要使用高精度的网格，对计算机资源的要求比较高，所以还不能在工程计算中被广泛使用。在绝大多数情况下，湍流计算还要采用湍流模式理论，大涡模拟则可以在计算资源足够丰富的时候尝试使用。 7.2.2 Spalart-Allmaras 模型 Spalart-Allmaras 模型是一方程模型里面最成功的一个模型，最早被用于有壁面限制情 况的流动计算中，特别在存在逆压梯度的流动区域内，对边界层的计算效果较好，因此经常被用于流动分离区附近的计算，后来在涡轮机械的计算中也得到广泛应用。 最早的Spalart-Allmaras 模型是用于低雷诺数流计算的，特别是在需要准确计算边界层 粘性影响的问题中效果较好。FLUENT 对Spalart-Allmaras 进行了改进，主要改进是可以在网格精度不高时使用壁面函数。在湍流对流场影响不大，同时网格较粗糙时，可以选用这个模型。 Spalart-Allmaras 模型是一种新出现的湍流模型，在工程应用问题中还没有出现多少成
Fluent软件的燃烧模型介绍(精)
Fluent软件的燃烧模型介绍 Fluent软件中包含多种燃烧模型、辐射模型及与燃烧相关的湍流模型,适用于各种复杂情况下的燃烧问题,包括固体火箭发动机和液体火箭发动机中的燃烧过程、燃气轮机中的燃烧室、民用锅炉、工业熔炉及加热器等。燃烧模型是FLUENT软件优于其它CFD软件的最主要的特征之一。下面对Fluent软件的燃烧模型作一简单介绍: 一、气相燃烧模型 ·有限速率模型 这种模型求解反应物和生成物输运组分方程,并由用户来定义化学反应机理。反应率作为源项在组分输运方程中通过阿累纽斯方程或涡耗散模型。有限速率模型适用于预混燃烧、局部预混燃烧和非预混燃烧。 应用领域:该模型可以模拟大多数气相燃烧问题,在航空航天领域的燃烧计算中有广泛的应用。 PDF模型 该模型不求解单个组分输运方程,但求解混合组分分布的输运方程。各组分浓度由混合组分分布求得。PDF模型尤其适合于湍流扩散火焰的模拟和类似的反应过程。在该模型中,用概率密度函数PDF来考虑湍流效应。该模型不要求用户显式地定义反应机理,而是通过火焰面方法(即混即燃模型或化学平衡计算来处理,因此比有限速率模型有更多的优势。 应用领域:该模型应用于非预混燃烧(湍流扩散火焰,可以用来计算航空发动机的环形燃烧室中的燃烧问题及液体/固体火箭发动机中的复杂燃烧问题。 非平衡反应模型
层流火焰模型是混合组分/PDF模型的进一步发展,从而用来模拟非平衡火焰燃烧。在模拟富油一侧的火焰时,典型的平衡火焰假设失效。该模型可以模拟形成Nox的中间产物。 应用领域:该模型可以模拟火箭发动机的燃烧问题和RAMJET及SCRAMJET 的燃烧问题。 预混燃烧模型 该模型专用于燃烧系统或纯预混的反应系统。在此类问题中,充分 混合的反应物和反应产物被火焰面隔开。通过求解反应过程变量来预测火焰面的位置。湍流效应可以通过层流和湍流火焰速度的关系来考虑。 应用领域:该模型可以用来模拟飞机加力燃烧室中的复杂流场模拟、气轮机、天然气燃炉等。 二、分散相燃烧模型 除了可以模拟各种气相燃烧问题以外,FLUENT5还提供了模拟分散相燃烧问题(液体燃料燃烧、喷射燃烧、固体颗粒燃烧等的燃烧模型:在拉格朗日坐标下,模拟分散相(包括固体颗粒/油滴/气泡等在瞬态和稳态下的运动轨迹 多种球形和非球形粒子的曳力规律 线性分布或Rosin-Rammler方程的粒子大小分布 连续相的湍流效应对粒子传播的影响 分散相的加热/冷却 液滴的汽化和蒸发 燃烧粒子,包括油滴的挥发过程和焦碳的燃烧
煤粉燃烧模拟湍流破碎模型
煤粉燃烧一一湍流破碎模型(EBU) 简介 该帮助文件主要介绍煤粉燃烧模型的设置和求解，采用湍流破碎模型 (EBU)。 EBU燃烧模型，也称涡团破碎模型，假设化学反应的平均速度与化学动力学无关，而只取决于低温的反应物和高温的燃烧产物之间的湍流混合作用。 主要包括： 1) 煤粉燃烧模型的建立和求解 2) 湍流破碎模型(EBU)的应用 3) 选择合适的求解参数 4) 计算结果的后处理 问题描述 3D模型的剖面图如图1所示。左侧为两个环形入口，右侧为一个圆形出口由于模型的对称性，取系统的1/4进行建模。煤粉和携带空气(一次风)从内环进入燃烧室，二次风从外环进入燃烧室，发生燃烧反应，产物从压力出口流出。模型建立和求解Stepl:网格 1、读取mesh文件 2、检查网格：grid -—? check 3、显示网格：display——> grid a) 从列表框中选取所有surfaces b) 点击display，并关闭grid display 面板 Step2模型 1、选择k- &湍流模型 Define —f models ——+ viscous 2、启动能量方程 Define ——? models ——? en ergy 3、启动物质输运方程
Define + models 一—k species-―tran sport & react ion a）选择model 列表中的species transport b）选择reaction 列表中的volumetric c）选择mixture material 歹U表中的coal-hv-volatiles-air d）选择turbulence-chemistry interaction （湍流与化学反应的作用）列表中的eddy dissipation （涡流耗散） e）点击ok, 关闭species modeI面板 4、启动Discrete ordinates模型（DO离散坐标系） a）从model 歹U表中选择discrete ordinates b）设置flow iterations per radiation iteration （流动和辐射迭代次数）为1 c）设置angular discretization （角离散化）中的Theta Divisions and Phi Divisions 为4 d）设置angular discretization （角离散化）中的Theta Pixels and Phi Pixels 为3 e）点击ok, 关闭radiati on model 5、启动discrete phase（离散相）模型 a）M ax number of steps （最大步数）40000 b）启动specify length scale （步长），设为0.0025 c）点击ok, 关闭discrete phase面板 Step3: injections （喷口） 1、v-1 入口截面设9 个喷口：define k injections a）点击creat,新建喷口 c）9个喷口通用性质，见表