燃烧过程数值模拟.pptxVIP

燃烧过程的数值模拟;§1 燃烧过程数值模拟的研究内容和方法;燃烧过程数值模拟的发展阶段;SIMPLE算法问世(1972)： 求解不可压流动时，如果对包含速度与压力的代数方程直接求解，则可同时得到速度场与压力场 ? 对计算机要求高.SIMPLE算法 ? 分离式的求解方法的基本思想是在流场迭代的求解的任何层次上，速度场必须满足质量守恒方程，从而保证流场迭代计算的收敛。 美国学者Thompson等人提出采用微分方程生成适体坐标的方法(TTM方法)： 为有限差分法与有限容积法处理不规则边界问题提供了一条新路：通过交换把物理平面上的不规则区域(二维问题)变换到计算平面上的规则区域，从而在计算平面上完成计算，再将结果传递到物理平面上。;2.走向应用期 (1975 – 1984);PHOENICS正式投放市场(1981) Rhie与Chou提出同位网格方法(1982)：吸取了交错网格成功的经验，又把所有变量的求解置于同一网格上。目前在非正交曲线坐标系中应用广泛。 80年代初，一批改革开放之初出国进修的访问学者相继学成回国，并在国内开展了CFD/NHT的教学与科研工作。;3. 蓬勃发展期 (1985 ? );一批有关CFD/NHT的新教材与参考书及期刊出版或创刊: 国外如： Anderson D A et al., Computational fluid mechanics and heat transfer, Washington: Hemisphere. 2nd ed. 1997 Minkowycz W J, Sparrow E M eds. Advances in numerical heat transfer. New York: Taylor Francis, Vol.1, 1997 国内如： 岑可法，樊建人，工程气固多相流动的理论及计算，杭州：浙江大学出版社，1990 范维澄，万跃鹏，流动及燃烧的模型与计算，合肥：中国科学技术大学出版社，1992 周力行，湍流气固两相流动和燃烧的理论与数值模拟，陈文芳，林文漪译，北京：科学出版社，1994 ;多个大型商业通用软件投放市场，比较著名的有: PHOENICS (1981);FLUENT (1983);FIDAP (1983);STAR-CD (1987);COMPACT (1989);FLOW-3D (1991,现改名为CFD)。 其中FIDAP采用有限元法（其余为有限容积法），后来并FLUENT，在上述软件中，目前在我国设立代理商的有PHEONICS,STAR-CD,CFX及FLUENT等 数值计算方法不断发展： 在网格生成技术方面，同位网方法得到进一步发展，非结???化网格的研究蓬勃展开；对流项格式的精度不断提高； 压力与速度耦合关系的处理中，提出了算子分裂算法PISO，SIMPLE算法系列化并推广到了可压流。 ;二、燃烧过程模拟的主要步骤和研究方法;在所建立的物理模型的基础上建立数学模型，构造基本守恒方程：连续方程、动量方程、能量方程、组分方程等。 上述基本方程通常不封闭。由物理概念或某些假设出发，提出模拟理论。 需要模化的分过程：湍流流动、湍流燃烧、辐射换热、多相流动和燃烧;2. 选择坐标系;3. 建立网格;4. 建立离散方程;有限差分法(Finite difference method, FDM): (1)这种方法将求解区域用节点所组成的点的集合来代替。每个节点所描述的流动与传热问题的偏微分方程中的导数项用相应的差分表达式来代替，从而在每个节点上形成一个代数方程，其中包含了本节点及其邻点上所求量的未知值; (2)在规则区域的结构化网格上，有限差分法十分简便而有效，并且很容易引入对流项的高阶格式； (3)难以保证离散方程的守恒特性；对不规则区域的适应性很差。 有限容积法(Finite volume method, FVM) (1)这种方法从描述流动与传热问题的守恒型控制方程出发，对它在控制容积上积分，在积分过程中需要对界面上被求函数本身(对流通量)及其一阶导数(扩散通量)的构成方式作出假设，从而形成不同的离散格式。其中扩散项一般采用相当于二阶精度的线性插值，故离散格式的区别主要体现在对流项上。 (2)有限容积法导出的离散方程可以保证具有守恒性；对区域形状的适应性较有限差分法好。;有限元法（Finite element method, FEM）