多相流动的基本理论.ppt

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,*,*,多相流体动力学,第二章 多相流动基础理论,2.1.4,多相流动理论模型和数值方法,特征时间,流动时间（停留时间）：,扩散驰豫时间：,平均运动驰豫时间：,流体脉动时间：,颗粒间碰撞时间：,无滑移流（平衡流）,强滑移流（冻结流）,扩散,冻结流,扩散,平衡流,稀疏悬浮流,稠密悬浮流,主要内容（气固多相流）,长期以来，气固两相流动的研究中按照对颗粒的处理方式不同，主要有两大类模型,确定轨道模型,随机轨道模型,单颗粒动力学模型（,SPD,模型）,颗粒轨道模型（,PT,模型）,小滑移模型（,SS,模型）,无滑移模型（,NS,模型）,拟流体（多流体）模型（,MF,模型）,离散介质模型,连续介质模型,本章要义,各种颗粒模型的一些基本观点,颗粒相模型,基本观点,颗粒对流体的影响,相间滑移,坐标系,颗粒相输运性质,单颗粒动力学模型,离散体系,不考虑,有,拉格朗日,无，扩散冻结,颗粒轨道模型,离散体系,考虑,有,拉格朗日,无（确定轨道）；有（随机轨道模型）,小滑移模型,连续介质,不考虑,有,（滑移,=,扩散）,欧拉,有,（扩散,=,滑移）,无滑移模型,连续介质,部分考虑,无（动力学平衡，热力学平衡或冻结）,欧拉,有,（扩散平衡）,拟流体（多流体）模型,连续介质,全部考虑,有,欧拉,有,按各种模型提出的时间大致顺序,无滑移模型,小滑移连续介质模型,滑移-扩散的颗粒群模型,分散颗粒群模型,双流体模型,颗粒轨道模型,拟流体模型（连续,-,连续介质模型）,前提：,在流体中弥散的颗粒相也是一种连续的流体；,气相和颗粒相是两种相互渗透的连续相，各自满足连续性方程、动量方程和能量守恒方程。,无滑移模型,（,No-slip Model),颗粒群看作连续介质，颗粒群只有尺寸差别，不同尺寸代表不同相；,颗粒与流体相间无相对速度；,各颗粒相的湍流扩散系数取流体相扩散系数相等；,相间相互作用等同于流体混合物间各成分相互作用，相间阻力不计。,基本假设：,小滑移连续介质模型,（,Soo,-drew Slip Model),颗粒群看作连续介质，不同尺寸组代表不同相；,各组尺寸颗粒群速度不等于当地的流体相速度，各颗粒相之间的速度亦不相等，即各颗粒相间、与流体相间有相对速度；,相间的相互作用类似于流体混合物中各种组分之间的相互作用，颗粒相和流体相间的阻力忽略不计；,颗粒的运动是由流体的运动而引起的，颗粒相的滑移是由于颗粒相对于多相流整体的湍流扩散所致，故这种小滑移也称为湍流飘移;,多相混合物整体与各相之间的关系，仍类似于多组分流体混合物和各流体组分间的关系.,基本假设：,滑移,-,扩散的颗粒群模型,(Slip-diffusion Model),各相时均速度差异造成滑移的主要部分，由于各相的初始动量不同引起；,扩散漂移造成滑移的小部分；,空间各点各尺寸组的速度、尺寸、温度等物理参数均不相同。,基本假设：,拟流体模型小结,无,滑移模型,：颗粒相的宏观运动而引起的质量迁移是由流体运动引起的；,小滑移模型,：混合物运动引起的,滑移,-,扩散模型,：颗粒相自身的宏观运动引起了质量迁移,拟流体模型数值方法,常用数值模拟方法,传统模式理论,直接模拟,大涡模拟,离散涡方法,格子气,湍流流场数值模拟方法简介,双方程模型,非线性 模型,多尺度 模型,RNG,模型,Reynolds,应力模型（,RSM,）,代数应力模型（,ASM,）,FLT,模型,SSG,模型,湍流模式理论以,Reynolds,时均运动方程和脉动运动方程为基础，依靠理论与经验的接合，引进一系列模型假设，从而建立一组描写湍流平均量的方程组。,湍流模式理论简介,对经验数据的依赖性;,将脉动运动的全部细节一律抹平从而丢失大量重要信息;,目前各种模型，都只能适用于解决一种或者几种特定的湍流运动。,湍流模式理论局限性,计算机发展,数值算法发展,直接模拟（,DNS),技术的应用,出现大型并行计算机,Petaflops,（,10,15,）,级,有限差分,有限元,谱方法,小波变换,自适应网格,并行计算技术,方程本身是精确的，不含任何认为假设和经验常数，仅有的误差只是由数值方法引入的误差。,计算包括脉动运动在内的湍流所有瞬时流动量在三维流场中的时间演变;,不用任何湍流模型，直接数值求解完整的三维非定常的,N-S,方程组;,湍流,直接模拟（,DNS）,简介,湍流旋涡结构包括大尺度涡和小尺度涡,大尺度涡,小尺度涡,湍 流 流 场 涡 结 构 图,直接模拟计算量太大，很难计算工程实际高雷诺数湍流流场。,为 什 么 要大 涡 模 拟？,湍 流 大 涡 模 拟 简 介,流,场,大尺度涡,小尺度涡,决定湍流流场的基本形态和性质；,流场质量、能量的主要携带者；,高度各向异性,无法建立统一模型,。,由大涡非线性作用产生；,流场能量的主要耗散者；,近似各向同性,可以考虑建立统一模型。,小尺度涡对大涡的影响用模型进行模拟,大涡模拟思想,对大尺度涡进行直接模拟,拟流体模型现状,为了能更完整地考虑颗粒相各种湍流输运特性以及相间的滑移和耦合，,Spalding,等,1,首先提出了双流体模型。,周力行教授对双流体模型进行了深入的研究。他们针对各向同性流动，提出了颗粒湍动能输运方程的模型,2,。针对各向异性流动，则将单相湍流流动的,RSM,模型推广至气固两相流中，提出了统一二阶矩模型（,USM）,3,。,拟流体模型现状,概率密度函数（,PDF）,方法被引用于构造双流体模型的两相湍流模型。,Reeks,4,从稳态流场中的颗粒运动方程出发，得到了颗粒相的,PDF,输运方程，同时还用,PDF,方法研究了近壁区颗粒的运动和自然边界条件的处理，克服了一般双流体模型难以描述的颗粒在壁面沉降、反弹过程的缺陷。,Zaichik,等,5,用,Rurutsu,-,Novikov,定理和泛函分析的方法，实现了采用,PDF,方法对流体湍流和颗粒相的模拟。,拟流体模型现状,Simonin,6,则运用流体涡团的,Lagrangian,模型来构造颗粒轨道上的流体涡团,Lagrangian,方程，从而得到了颗粒相的连续、动量和,Reynolds,应力方程。,周力行教授等采用了二阶矩封闭的思路来封闭,PDF,输运方程中湍流与颗粒的相间作用项，将颗粒相的,PDF,模型与流体运动的各类模型相结合，提出了,k,PDF,模型,7,和,FDSMPDF,模型,8,。,双流体模型,双流体模型,固相压力,固相的剪切粘度,固相的体积粘度,固相的应力张量,气相牛顿粘性应力方程,离散颗粒模型,确定轨道模型,随机轨道模型,单颗粒动力学模型（,SPD,模型）,颗粒轨道模型（,PT,模型）,.,共同特点,气相,颗粒相,欧拉系,拉格郎日系,不同之处,单颗粒动力学模型,颗粒轨道模型,考虑已知流场中颗粒平均运动或对流运动的轨道，忽略颗粒对流场的影响。,充分考虑气相和颗粒相间的相互作用。,分散颗粒群模型,基本假设：,在欧拉坐标系中考察流体相的运动情况，而在拉格朗日坐标系中研究颗粒群的运动情况。,即把颗粒群按初始尺寸分组，各组颗粒沿其自身轨道运动。,由于颗粒的蒸发、挥发及燃烧、流体的阻力作用和传热等原因，颗粒群沿轨道会发生速度、质量、温度、密度和尺寸的变化，同时对流体造成了分布于整个体积中的物质源、动量源和能量源。,该方法能研究颗粒群和流体相之间的较大滑移，并把复杂的颗粒变化情况耦合进来。,按照是否考虑,颗粒群的湍流扩散,，又可把颗粒轨道模型分为两类：,一类是不考虑颗粒群湍流扩散的颗粒确定轨道模型，,一类是考虑颗粒群湍流扩散的颗粒随机轨道模型。,颗粒确定轨道模型,处理颗粒群的方法较简单，能够考虑相间速度与温度的滑移，,并可以追踪比较复杂的颗粒经历，,数值计算不会产生伪扩散。,但其存在一个缺点，就是对颗粒的湍流扩散缺乏较好的处理。,考虑到,湍流脉动对颗粒轨迹,造成的影响，,Yuu,等,142,首先提出了涡作用模型。,在经过,Gosman,等,143,和,Berlemont,等,144,改进以后，得到了广泛的应用。,Sommerfeld,145,和,Shuen,146,等采用此模型进行数值求解，得到了比较满意的结果。,浙江大学热能工程研究所的岑可法院士和樊建人教授,147,提出的随机频谱颗粒轨道（,FSRT）,模型，,颗粒随机轨道模型。,模型小结,各种不同的气固两相流动模型，从不同的角度对真实的气固两相流动过程做了近似和简化，因而具有不同的适用范围。,对稀疏多相流动中固体颗粒，液体颗粒以及气泡运动的计算方法，,Loth,159,做过较为详细的介绍和分类。一般情况下可通过,判断颗粒相对浓度和相间滑移量的大小,来选择合适的模型。,不过随机轨道模型计算时需要跟踪大量的颗粒轨道，因而造成计算机的存储量和计算量都很大，从而使其在工程应用上受到一定程度的限制。,从已有的研究来看，在湍流气固两相流动的数值模拟方法中，,颗粒轨道模型,的应用最为广泛。它的优点在于计算工作量小，能够模拟有蒸发、挥发、两相化学反应和在不同阶段有不同质量损失率的颗粒相的复杂经历，而且颗粒相采用拉格朗日坐标系处理可以避免伪扩散。,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=35,个尺寸数，以充分地描述分散颗粒群的运动规律，,颗粒的进口速度有四种设定方法：,(1),设颗粒进口速度为零，即 ，这相当于颗粒由静止状态被气流曳引加速。,(2),设颗粒气口速度和气流进口速度一样，。相当于颗粒在管道内已被气流充分加速。,(3),大于或小于 ，视颗粒在管道中的加速或减速情况而定。,(4),设颗粒和气流的速度相差颗粒终端沉降速度，即这是达到稳定的一种假设。,一般情况下，我们都认为在进口处颗粒的速度及温度都均匀分布，当然也可以设颗粒速度按一定规律分布，并且不同尺寸组的颗粒具有不同的初始速度，这样更接近于实际的情况，但要消耗更多的计算时间。,颗粒相进口条件,，,颗粒相边界条件,Grant.G,等经验公式,Sommerfeld,冲量法,无滑移,有滑移,考虑粗糙度时对壁面的处理方式,随机数模拟虚拟壁面倾角,虚拟壁面与不规则反弹,正弦表面法,颗粒相边界条件,Reflect,Escape,Trap,Interior,颗粒发生弹性或非弹性碰撞反射,穿过壁面而逃逸（颗粒的轨道计算在此处终止）,在壁面处被捕集，非挥发性颗粒在此处终止计算，颗粒或液滴中的挥发性物质在此处被释放到气相中,穿过内部的诸如辐射或多孔介质间断面区域,颗粒轨道法（浓相）,硬球模型,软球模型,软球模型,基于软颗粒模型的,DEM,方法最早是由,Cundall,和,Strack,提出的（,Cundall,and,Strack,1979）,用于计算土壤力学的，它认为颗粒在碰撞时会产生变形，是,“,软,”,的，颗粒间的碰撞力由颗粒的变形和颗粒的弹性模量决定，随着颗粒变形的增加颗粒间的相互作用力也相应增加，颗粒在碰撞过程中是一个变加速度的过程。,两个颗粒在碰撞过程中可以有第三个颗粒再次碰撞过来，可以计算多体碰撞的情况。,softsphear,.doc,硬球模型,硬颗粒模型主要由颗粒的动量守恒方程，结合颗粒的恢复系数，牛顿第二定律计算颗粒碰撞前、后的速度，它认为颗粒间的碰撞是完全弹性的碰撞，碰撞前后能量、动量守恒。,同时这个模型还假设颗粒间的碰撞只有两两碰撞，因此不能处理多体碰撞的情况。,hardsphear,.doc,拉格郎日轨道法,颗粒源项计算,在用拉格朗日方法求解颗粒运动时，每一条轨迹代表一组粒径相同的颗粒的运动情况。设共有,n,个轨道穿过网格,k,，,其中第,j,个轨道含,n,j,个质量为,m,pj,的颗粒。则该组颗粒通过网格,k,时的动量源项为：,所有轨道留下的源项为：,按上式计算出所有网格的动量源项，并代入气相控制微分方程组即可进行动量耦合。,颗粒的浓度场,颗粒相的浓度场和速度场的确定.,doc,颗粒的速度场,颗粒相的浓度场和速度场的确定.,doc,FRST,模型迭代计算步骤是,(1),用,SIMPLEST,方法求解气相场,(,初次计算时可不考虑颗粒相的影响,),，采用,k-,双方程湍流模型考虑气相湍流运动，得出 、，,k,，,值。,(2),由,k,和,值计算气相脉动速度的幅值,u,i,，,v,i,，,w,i,。,(3),用随机的,Fourier,级数模拟气流的脉动速度，,并求解颗粒的速度,u,p,，,v,p,和,w,p,及颗粒的轨迹,x,p,，,y,p,，,z,p,。,(4),由体积平均法求解颗粒的速度及浓度场。,(5),颗粒相对气相的耦合作用以源项表示，将源项 、代入气相方程进行求解。,(6),重复上述迭代过程，直至气相场和颗粒相速度场及浓度场收敛为止。,拉格郎日轨道法,程序开始,（设初始条件和进口条件）,计算单相气相场到收敛,根据所求出的气相场,计算颗粒的速度、轨迹、温度等,计算颗粒源项,将,源项代入气相场，并再次迭代气相场到收敛,收敛,计算颗粒速度和浓度场,停机,参考文献,Spalding D B.A general purpose computer program for multidimensional one-and two-phase flows.J.Math,Comput,.,Simul,.,1981,23:267-276,周力行，黄晓晴三维湍流气粒两相流的,k,kp,模型工程热物理学报，,1991,，,12(4),：,428-433,周力行湍流两相流动及燃烧的统一关联矩封闭模型工程热物理学报，,1991,12(2),：,203-209,Reeks M W.PDF modeling of gas-particle flows.In:2,nd,Int.,Symp,.On Multi-phase Fluid,Non-Newtonian Fluid and Physicochemical Fluid Flows97,1997,Beijing,Zaichik,L I,Alipchenkov,V M.Simulation of transport of colliding particles suspended in turbulent shear flows.In:2,nd,Int.,Symp,.On Turbulence,Heat and Mass Transfer,1997,Delft University,Simonin,O.Continuum modeling of dispersed turbulent two-phase flows.VKI Lectures“Combustion in Two Phase Flow”,1996,李勇，周力行,kPDF,两相湍流模型和台阶后方气粒两相流动的模拟工程热物理学报，,1996,17(2),：,234-238,周力行，李勇旋流两相流动的,DSMPDF,两相湍流模型工程热物理学报，,1999,20(2),：,252-257,周力行，湍流气粒两相流动和燃烧的理论与数值模拟,【,专著,】,，科学出版社，,1994,，北京。,岑可法，樊建人，工程气固多相流动的理论与计算,【,专著,】,，浙江大学出版社，,1990,，杭州。,