密相液固两相流动的数值研究.pdfVIP

1=) A辑第 12卷第 3期 水 动 力 学 研 究 与 进 展 Ser．A，Vo1．12，No．3 1997年 9月 JOURNALOFH-YDRODYNAMICS Sep．，1997 密相液固两相流动的数值研究 魏进家 姜培正 、／王长安 。 j 大学 王 安710049) 弓 摘 要 本文推导了适用于考虑浓度变化影响的密相藏固两相流动数值计算的DIPSAR算 法 ，对竖直上 流动进行了数值计算，计算值与实测值较为符合。 关键词 DIPSAR算法 0 引 言 液固两相流动广泛存在于 自然界和许多工业部门，在诸如环保、水利、化工、海洋工程等领 域均可遇到，因此，研究液固两相流动规律具有重要的意义。对液固两相流动进行数值研究，是 两相流研究的一个热门课题。SpaldingEI]在计算不可压单相流体的SIMPLE算法的基础上，提 出了适用于两相流动的IPSA(Inter—phase—Slip—Algorithm)算法。王长安 又借用了SIMPLER 算法的思想，对IPSA进行了改进，提出了IPSAR(IPSA—Revised)算法。IPSA和IPSAR算法仅 考虑了速度修正值而得到了浓度加权压力修正方程，没有考虑浓度修正值的影响，然而当颗粒 相浓度较高时，浓度变化也较大 ]，因此该算法只适用于浓度变化不大的流动。本文在 IPSAR 算法的基础上，考虑浓度变化的影响，推导了适用于密相液固两相流动的DIPSAR(Dense—IP— SAR)算法。 简单地说，DIPSAR算法与 IPSAR算法的主要区别是在浓度加权压力修正方程中引入了 由于浓度修正值所导致的附加项 ，因而通过求解浓度加权压力修正方程能够反映出浓度变化 的影响。经过这样的改进后的DIPSAR算法可适用于密相液固两相流动的计算。 1 控制方程 任意曲线坐标系下，稳态的液固两相流动方程可表示为0： J(Ul 一 g )]一J (1) 其中． 为Jacob行列式；矿 为逆变度量张量；U1为k相速度的逆变分量； 为 的扩散系 数f 分别为 1、U (一1，2，3)、K，、及 U 为相速度的协变分量 ；K1为流体相紊动能； 为流体相紊动能耗散率；T 为颗粒相温度； 为方程源项。 将方程(1)按体积积分法对控制体积分并离散得： [ J (： 一 g,- )]△xz△x；+ [ (i 一v*ig )]：△ △群 本文于 1996年 1o月31日收到。国家 自然科学基金资助项 目。 魏进家等 ：密相液固两相流动的数值研究 351 + [ J(ul 一 g )1nx~nxl= △ △ ；△ ； (2) 式中，e、w…tb 、分别表示该控制体的东、西、前、后、北、南界面，△ 、nxl、△ 表示东 西、前后、南北方向控制体长度 。 在(2)式中，令 =1，并由Ui一 得： [ Jg”U ]axIaxl一 [ Jg“un]△ ；△ ； + [ Jg。 ]I△ ；△ ；一 [ Jg”Hl△ ；△ ； + 以。]△ }△ ；一 ”以 △ △ ；一J △ ；△ △ ；一6一 (3) 式中～b 为连续方程的曲率源项。 6一 一 [ J(g” z+g”U的) △ ；△ ；：4-[ J(guH+g仉。)△ ；△ ； + [ (guH+gu此)卫△ △ ； (4) 对于 (3)式左边第一项，将U 一u ．+U ¨， ．= + 代入并