专利名称：一种测量气固两相流动中颗粒运动速度的方法
技术领域：
本发明涉及气固两相流动领域，具体涉及一种测量气固两相流动中颗粒运动速度的方法。
背景技术：
气固两相流动广泛存在于能源、化工、环保以及气象等领域。对气固两相流动的研究和应用，除了采用理论和实验方法外，随着计算机计算速度的快速提高，气固两相流动的计算机数值模拟在科研及工程实际中正发挥出越来越大的作用。在气固两相流动的计算机数值模拟过程中，其关键环节是建立描述气固两相流动的数学模型，目前在稀疏气固两相流动(固相容积份额在5%以下)的数值模拟中，由于固相之间的作用相对较弱，数学模型的构建已比较成熟，但对于稠密气固两相流动(固相容积份额大于5%)而言，固相间相互作用的数学描述仍是该领域的一个难点。目前，描述固相间相互作用的数学模型较多，但应用最为广泛和具有代表性的可归纳为以下三种模型双流体模型、软球模型和硬球模型。双流体模型基于欧拉方法，参照描述气相(牛顿流体)的数学模型来描述固相，其主要不足是随固相浓度增加，固相流动特征逐渐偏离牛顿流体，以至于该模型不适用于稠密气固两相流动的数值模拟；软球模型通过运动轨迹跟踪流场内的每一颗粒，准确、细腻地描述颗粒碰撞接触过程的作用力，可以模拟稠密气固两相流动，但关键不足是计算量巨大，目前常用PC 机只能对数万个颗粒所构成的气固两相流动进行数值模拟；硬球模型往往结合蒙特卡诺方法，基于概率来判别颗粒是否发生碰撞，可进行大规模气固两相流动的数值计算，但由于根据概率判别颗粒间是否发生碰撞，使得颗粒碰撞过程模拟欠缺真实性，并且硬球模型的基本假设是两体瞬时碰撞，而当固相浓度增大到一定程度时，例如发生堆积情况下，该模型已不适宜用于数值模拟。综上所述，目前对稠密气固两相流动缺乏能够进行高效、大规模数值模拟的数学模型。
发明内容
发明目的针对现有模拟方法存在的不足，本发明的目的是提供利用一种描述气固两相流动中固相之间相互作用对颗粒运动影响的数学模型，对稠密气固两相流动进行高效、大规模数值计算的方法。技术方案一种测量气固两相流动中颗粒运动速度的方法，该方法包括以下步骤
1)将气固两相计算区域划分成网格；
2)设定气相和固相的入口、出口、壁面边界条件并采用欧拉法模拟得到计算区域中包括流场速度u在内的气相流场的信息；
3)在单个时间步长中，根据颗粒受到的除固相间相互作用力之外的力，计算出计算区域内每个网格中所含颗粒的在该时间步长中的加速度a及速度从而得出颗粒新的位置；
4)根据步骤3)中计算得到的各颗粒的位置，扫描计算区域内的所有网格，得出各个网格的固相容积份额/，并建立数组储存；将步骤3)中得到的计算区域内每个网格中所含颗粒在单个时间步长中的速度V,建立数组储存；
5)将步骤4)得到的各网格固相容积份额/以及各网格中所含颗粒的速度V，作为颗粒因受固相之间相互作用而影响运动的主要因素，由
权利要求
1.一种测量气固两相流动中颗粒运动速度的方法，其特征在于该方法包括以下步骤1)将气固两相计算区域划分成网格；2)设定气相和固相的入口、出口、壁面边界条件并采用欧拉法模拟得到计算区域中包括流场速度u在内的气相流场的信息；3)在单个时间步长中，根据颗粒受到的除固相间相互作用力之外的力，计算出计算区域内每个网格中所含颗粒在该时间步长中的加速度a及速度K/)，从而得出颗粒新的位置；4)根据步骤3)中计算得到的各颗粒新的位置，扫描计算区域内的所有网格，得出各个网格的固相容积份额八并建立数组储存；将步骤3)中得到的计算区域内每个网格中所含颗粒在单个时间步长中的速度V,建立数组储存；5)将步骤4)得到的各网格固相容积份额/以及各网格中所含颗粒的速度V，作为颗粒因受固相之间相互作用而影响运动的主要因素，由 Sv :aux ^,=1 - WifnJxW - ι Σ τ/ (0) ι / Σ ]}求得Fr^r^3-1K -1K“W“ W，进而得到各颗粒新的速度Vj ；式中 分别代表n ζ方向；为本方法中新定义的H ^方向的运动约束系数；为各网格中所含颗粒受固相间相互作用前H^方向的速度；Vim为各网格中所含颗粒受固相间相互作用后X、_7、ζ方向的速度；i代表颗粒的编号； /为当地网格的固相容积份额；为当地网格中颗粒达到堆积饱和状态时的固相容积份额； Ini为当地网格中各个颗粒的质量，根据所计算颗粒的密度及平均粒径来确定；6)首先判断已执行步骤3—步骤5循环的次数是否达到设定值nl，设定值nl=计算时间/时间步长，若达到设定值则执行步骤7)；若未达到设定值，则判断已执行步骤3—步骤 5循环的次数是否为所需更新气相场次数n2的整数倍，n2<nl，若是，则回到步骤2，若否，则回到步骤3 ；7)完成数值计算，输出结果。
2.根据权利要求1所述的测量气固两相流动中颗粒运动速度的方法，其特征在于，所述的步骤3)中，所述的除固相间相互作用力之外的力至少包括流体力及、重力dXg，颗粒在该时间步长中的加速度及速度KirwIKitoitVaXi^，式中Fsum每个颗粒所受的合力； Fswn=FS+G, Dt为时间步长； ν/init)为颗粒上一个时间步长计算得到的速度。
3.根据权利要求2所述的测量气固两相流动中颗粒运动速度的方法，其特征在于，流体力F=fi ~ 一 W，式中u是流场速度，κ是颗粒当前计算中的速度； β=0. 75 C0 \ ν - u \ P (I - ε) ε ^1/ Dp,G是流体的阻力系数，当彻<=1000时，Cd =24(1+0. 15 RH讼，当7fe>1000时，Cd =0. 43 -,Re= Dp \ ur \ / v, Re为流动雷诺数，&为每个颗粒的直径，u=u 一 κ ； P为流体密度，ε为空隙率1 一八
4.根据权利要求1或2所述的测量气固两相流动中颗粒运动速度的方法，其特征在于， 所述的时间步长为KT9-KT1秒。
全文摘要
本发明提供了利用一种描述气固两相流动中固相之间相互作用对颗粒运动影响的数学模型，对稠密气固两相流动进行高效、大规模数值计算的方法。本发明利用的数学模型构建了气固两相流动中颗粒间发生的碰撞、接触及摩擦对颗粒运动所产生的影响与颗粒相占据的容积份额、包括速度大小、方向在内的颗粒群运动特征之间的数学关系，利用该模型能够对稠密气固两相流动进行较为准确地数值模拟计算，即使在颗粒发生堆积的情况下也能使数值模拟结果与实验结果较好地吻合。此外，该模型所涉及的计算量与现有模型相比大大减少，利用目前常用计算硬件设备也可进行大规模数值模拟计算。
文档编号G01P7/00GK102230943SQ20111008765
公开日2011年11月2日 申请日期2011年4月8日 优先权日2011年4月8日
发明者彭小敏, 朱立平, 袁竹林 申请人:东南大学