多相流动的基本理论教程.pptVIP

不过随机轨道模型计算时需要跟踪大量的颗粒轨道，因而造成计算机的存储量和计算量都很大，从而使其在工程应用上受到一定程度的限制。 从已有的研究来看，在湍流气固两相流动的数值模拟方法中，颗粒轨道模型的应用最为广泛。它的优点在于计算工作量小，能够模拟有蒸发、挥发、两相化学反应和在不同阶段有不同质量损失率的颗粒相的复杂经历，而且颗粒相采用拉格朗日坐标系处理可以避免伪扩散。 Crowe等[183]和先后对气固两相湍流流动的数值模拟方法进行过概括总结。在Mashayek等[184]的综述中，他们对最新的气固、气液两相流动的数值模拟方法进行了详细的介绍，包括了拉格朗日描述的直接数值模拟、大涡模拟和统计模型，以及欧拉方法描述的RANS模型和PDF模型等， 颗粒轨道法 对稀疏两相流来说，颗粒的存在对气相影响很小，可不予考虑，这种情况被称为单向耦合（One-way Coupling），即只认为气相运动特性单方面影响着颗粒的运动情况。 而对于浓度较高的气固两相流动，不仅气相影响着颗粒的运动，而且颗粒对气相运动也有明显的影响，不应被忽略。这种同时考虑颗粒和流体间相互作用的情况被称为双向耦合（Two-way Coupling）。 如果再进一步考虑颗粒间的相互碰撞，则被称为四向耦合（Four-way Coupling）。 流体相被看作为连续介质，而颗粒相被看作与流体有滑移的，沿自身轨道运动的分散群; 颗粒相自身无湍流扩散; 颗粒群按初始尺寸分组，各组颗粒群沿各自轨道运动，互不干扰; 颗粒群对流体的质量、动量和能量相互影响当作是某种等价的连续分布于多相流空间中的物质源、动量源和能量源。 拉格郎日轨道法流体相方程 直角坐标系中三维流动微分方程式各项的意义 直角坐标系中三维流动微分方程式各项的意义 方程名称 ? 连续性方程 1 0 0 X方向动量方程 u Y方向动量方程 v Z方向动量方程 w 湍流动能 k 湍流动能耗散率 ? 表中，u、v、w为x、y、z方向的速度分量，为湍流脉动能的产生项： 有效粘性系数 ，其中湍流粘性系数 ，C?=0.09， C1=1.47， C2=1.92， ??=1.3， ?k=1.0. 气相流体控制微分方程组的数值解法 对气相流体控制控制微分方程组的求解采用SIMPLEST方法。主要步骤如下： 1. 估计整个积分区域的压力分布P*； 2. 用雅克比逐点校正法解动量方程，得到速度场u*，v*，w*； 3. 建立和求解压力校正方程，得到 ； 4. 求速度校正值 ， 和 ，得到校正后的速度分布u=u*+ 等； 5. 校正压力分布，p=p*+α ，其中α为松驰因子； 6. 把求出的p作为下次迭代的估计值，重复(1)到(5)，直到收敛。 计算中采用低松驰，即α1。 拉格郎日轨道法 颗粒相方程 颗粒在湍流脉动中的扩散 湍流脉动对颗粒运动的影响.doc 脉动频谱随机轨道模型 脉动频谱随机轨道模型.doc 1.进口条件 在进口处须给出颗粒位置的计算站j，计算站越多，则最后求解出的颗粒的浓度场和速度场就越接近于实际情况，但这样做的缺点是计算时间增大。同时，我们一般以几档离散的颗粒直径来表示颗粒尺寸的连续分布，一般取i=3~5个尺寸数，以充分地描述分散颗粒群的运动规律， 颗粒的进口速度有四种设定方法： (1) 设颗粒进口速度为零，即 ，这相当于颗粒由静止状态被气流曳引加速。 (2) 设颗粒气口速度和气流进口速度一样， 。相当于颗粒在管道内已被气流充分加速。 (3) 大于或小于 ，视颗粒在管道中的加速或减速情况而定。 (4) 设颗粒和气流的速度相差颗粒终端沉降速度， 即这是达到稳定的一种假设。 一般情况下，我们都认为在进口处颗粒的速度及温度都均匀分布，当然也可以设颗粒速度按一定规律分布，并且不同尺寸组的颗粒具有不同的初始速度，这样更接近于实际的情况，但要消耗更多的计算时间。 颗粒相进口条件 ， 颗粒相边界条件 Grant.G等经验公式 Sommerfeld冲量法 无滑移 有滑移 考虑粗糙度时对壁面的处理方式 随机数模拟虚拟壁面倾角 虚拟壁面与不规则反弹 正弦表面法 颗粒相边界条件 Reflect Escape Trap Interior 颗粒发生弹性或非弹性碰撞反射 穿过壁面而逃逸（颗粒的轨道计算在此处终止） 在壁面处被捕集，非挥发性颗粒在此处终止计算，颗粒或液滴中的挥发性物质在此处被释放到气相中 穿过内部的诸如辐射或多孔介质