第十三讲煤燃烧过程数学模化幻灯片.pptVIP

控制容积公式的最吸引人的特征是：所得到的结果(解)将意味着任何一组的控制容积内，当然也就是在整个计算域内，诸如质量、动量以及能量这样一些物理量的积分守恒都可以精确地得到满足。对于任意数目的网格结点，这一特征都存在(不只是限于网格结点数变得很大时的极限意义上)。因而，即便是粗网格的解也照样显示准确的积分平衡. 交错网格系统 在这种交错网格系统中，除速度以外的各个变量，例如压力、密度、温度、组分浓度、比热、交换系数等均采用同一个网格系统；三个速度分量分别采用它独用的网格系统。如果压力和温度等变量的网格点为P，那么速度分量的网格点将位于网格P的与x方向垂直的界面的中心，v分量的网格点位于网格P的与y方向垂直的界面的中心，w分量的网格点位于网格P的与z方向垂直的界面的中心。因此，在同一个三维积分区域内就有四组不同的网格交错排列。 差分方程的求解 三对角线矩阵法 SIMPLE方法 . 估计速度分布； . 用式(9-214)计算准速度分布； . 计算压力方程(9-216)的系数，并解出压力分布； . 把解出的压力分布作为p*，解动量方程得到u*和v*； . 建立并求解压力校正方程，得到p?； . 求出速度校正值u?和v?，校正速度分布，但不用p?校正压力分布； . 如果需要，求解其它变量?的方程； . 重复上述第2至第7步，直至收敛。 煤粉火焰综合求解及示例 求解程序 在给定一个合理的，收敛的流场，调用颗粒两相流的计算，首先只考虑动量藕合，而不考虑燃烧和反应，这个藕合通常是很微弱的，但是若在入口处煤与气之间存在着很大的速度差，则在流动与两相流动之间迭代几次，再加上颗粒反应，这时再调用气相湍流燃烧、修正气体温度、化学组成，混合物分子量及流体质量密度。在开始使用气相湍流燃烧程序时，可以根据进气是预混的还是未混合的来选择两种不同的方法。迭代直接开始去计算需要的平衡态解或化学动力学解。这样迭代计算，如果收敛且计算结果合理，则可将辐射换热加入到循环之中，最后，要利用各个模型的最精确的公式，得到满足全部方程及边界条件的最终稳定解。 模型计算结果示例 杨伯翰大学一维模型计算结果 二维模型计算及结果 三维计算结果实例 颗粒运动的轨迹 水分蒸发规律 挥发份析出规律和着火时间、距离 焦炭燃烬规律 一次风反切时炉内温度场的数值模拟 煤粉颗粒扩散及两相流模型 描述两相流动的基本方法 一类是通过拉格朗日方法，一类是欧拉方法 稀相流动的基本分析 容积比V颗粒/V汽?1%； 颗粒碰撞不频繁； 单个颗粒的传热、传质、阻力系数不直接受邻近颗粒影响。 最早进行气-固两相问题的完整尝试可以追溯到Migdal和Agoata[852]1967年给出的方法，拉格朗日方法提出的考虑两相耦合问题的方法主要是PSIC法，PSIC的主要思想是： 首先解无颗粒存在的气相场，在此场中求解各组颗粒的轨迹，沿轨迹求颗粒、尺寸、速度、温度等经历变化，在气相单元边界上记录质量，动量和能量、组分、提供气相求解的源项，再解气相场，反复直到收敛。 颗粒在气流中的受力分析 即是颗粒受力的总和，一般包括：稳态气动阻力，浮力，虚假质量效应，Basset力，Magnus效应，Saffman升力和热泳力等，在颗粒反应过程中由于表面质量流动不均匀引起的推力作用 颗粒的湍流扩散 离散相本身的湍流输运(湍流扩散)； 由于离散相存在的输运特性对连续相湍流特性的修正； 湍流脉动引起的相间传输率特性的修正。 两相模型 数值求解方法 各项相应称为时间导数项、对流项、扩散项和源项 . 形式相同。这一点使我们有可能发展建立求解的通用程序。 . 非线性。即存在因变量或它们的导数的非一次项。它主要表现在对流项和化学反应源项。其它源项和扩散项也可能包含非线性，例如湍流脉动动能输运方程的源项，考虑交换系数为组分浓度和温度函数时的扩散项。在变密度的过程中，动量、能量和组分方程中的时间导数项也具有非线性。这表明基本方程一般来说无法用解析方法而必须用数值方法求解。 . 耦合。构成方程组的各方程不互相独立，因变量交错地存在于各个方程中。因此不能直接求解，而必须选代求解。 离散化方法 一个微分方程的数值解系一组?的空间和时间分布的离散的数的集合。这与在实验室中进行的实验相类似。数值分析方法和实验方法两者所获得的数据都只能是一些有限数量的数值；就实际需要而言，这个数量的数值是够实用就可以了。 所谓的数值方法就是把计算域内有限数量位置(叫做网格结点)上的因变量当作为基本的未知量来处理。该方法的任务是提供一组关于这些未知量的代数方程并规定求解这组方程的算法。 离散化的概念 用离散的值取代了包含在微分方程精确解中的连续信息，就将注意力集中在网格结点处的值。这实际上已经离散了?的分布，这一类数值方法叫做离散化方法。 这样，就从另一个