湍流-多相流-数值理论

1、CFD 模拟CAD实验验证优化和放大设计实验+半经验计算.模型原型nmPrReANu (Q-Q3-Q4)/Q=燃烧效率; Q1/Q=热效率P Pe e T Te e Y Ye e燃料燃料+ +空气空气P P0 0 T T0 0 Y Y0 0Q QQ Q1 1水水Q Q3 3+Q+Q4 4Q Q5 5蒸汽蒸汽Q Q2 2烟气烟气只能预报各变量沿轴线方向的变化燃料燃料 + +空气空气烟气烟气Y YF0F0T T0 0u u0 0T TF FY YFeFeT Te eu ue e(2D) 21 Mar1999 VECTOR DAT FILE00.05 0.1 0.15 0.2Z00.050.10.1

2、50.20.250.30.350.40.450.50.55X386386472472472558558558644644644730730730816816902902902902988988988107410741160116012471247133313331419141915051505T(K) Y=0.04m(2D) 21 Mar1999 VECTOR DAT FILE(2D) 22 Mar 1999 VECTOR DAT FILE00.05 0.1 0.15 0.2Z00.10.20.30.40.5X888888888888881761761765291762652652652653

3、53353353529618Left Wall(Kw/m )2(2D) 22 Mar 1999 VECTOR DAT FILE建立或选择数学物理模型发展或选择数值方法建立或选择计算程序用实验数据检验模拟结果弥散型多相流泡状流塞状流分层流波状流弹状流环状流泡状流出现于比较大的液体速度和不太大的气体速度的情况泡状流出现于高液体速度和低气体速度的情况时间平均周期 要比湍流脉动的时间尺度大得多，同时又比宏观非定常流动，例如波浪运动的周期小得多00/; 可压缩流动的Favre平均（密度加权平均）通用变量可以是速度分量v、温度T或组分浓度YPDF-p(f) 的定义是0到1之间变化的函数f在f 到 f+df

4、 之间出现的概率是 p(f)df )(limlim)x, a(pTTT10a2101022221010)df)f (fp(df)f (pf)f (ffdf)f (fpfdf)f (p1210102210)df)f (p)f (df)f (p)f (df)f (p)f ()f ()f ()1()f ()f (p2/ )ff (f222)ff ()ff (fgff2)ff ()2(1222/12exp)f ()1 ()f ()f (p101b1a1b1adf)f1 (f)f1 (f)f (p)x()x ;t ()x , t (2)x,(R0d)x,(R)x(T)x(b)x(a)x , t (b)x

5、; t (aab22)x; t (R)x()hx; t ()x , t (2)x, h(R11hdh)x;h(R)x(Lm1同一变量在同一位置不同时刻的关联系数积分时间尺度同一位置同一时刻两个变量的关联系数同一变量同一时刻两个位置的关联系数积分长度尺度 /fL v /fL v /fL v2/(18 )rppmd2/311(1 Re/6)rrp /Tl uk /Tl uk211/() ()pppp ppluc n ru1/rfSt/rTHt1/rpSl 1/rpSlSt1 冻结流动或强滑移流动Ht1 扩散冻结流动Sl1 稠密流动 )/(exp)(nkkdddR尺寸大于dk 的颗粒的百分数平均粒径

6、Sauter直径的应用最广泛 3(/6)mpkkkpnd pppfppp/11;/)/1 (pp)/(0000uupp)/(0000uupp各种密度容积分数稀疏流动质量负荷 pdcRe/24 Re/24)6/Re1 (3/2pdccd = 0.44 3/23/gp BBccdd/ )1ln(0)/()1ln(NuDdmBp0.50.33p0.50.33pNu20.6RePrSh20.6ReSc rpdtduuup/ )()/exp()(0rpptuuuu)/exp(0rpptvv)1 ()1 ()(/0/0rrtrpptrppevyeuutux当 t up u vp 0 y=vpr 0)/1

7、(gpfvdtdvrpp)e1 ()/1 (gvr/ tpfrp流体是液体或气体颗粒也可以是|3pppMvvdF 2/122/1|)(6 . 1yvppsvvdF qEdFdTFppExTppTjj32)6/()2()/(5 . 4 331/223()0.5(/6)()(/6)1.5()()()pippidvpdivmipiBiMisiTiEidtdpipipipidttdvdppipidtdMisiTiEimFFFFFFFFdvvdvvddvvt dFFFFVjxtVjxtSntdtdV)v(V)v(dSvjjVV0S)v(0VSjxtVdtdjjxjtdtdv00Sdtddtdxjtxjx

8、vjxjxtjjjjjjvv)v(0)v(0jjxvdtdssYVsSsVVwdSjVSsxYxjsxstw)D()vY()Y(jsjjsxYxdtdYw)D(jsjsivVVsisxpVvidtdV)F(VSVvjijijxvxvijVijVxitidtdp)(pV)vv()v(VvijjijsisxvxvxxpijxitF)()vv()v(ijjijj温度形式的能量方程ijjpipipipprpixvxvxiijji vvx-vnv-vnm+g+)(+xp-=vvxvtjiijj hvx-QW+xTx=hvx+htjjssjjjjjsjSjSjSjjsvYx-W-xYDx=Yvx+Yt0=

9、vxtjjkkdVkdVdVdVkkkkkdVkdVkkVdVdVdVdVVkkkk/11 Svjxtj)( Svjxtj)( kkjkxtSvjk)( ()()()jijijkrkpij iitxxxkiiikvv vgvvv S kikiiikkikjkxkiktSvvvgvvvrkkj)()()(体积平均后的层流多相流动的连续和的动量方程 TScQnqQwTcvTcpkkrssxTxtppjxptjjj)()()( kprkkhkkkkjkxkkktTScQQQnTcvTcj)()()( SwDYvYssxYxsjxstjsjj)()()(体积平均后的层流多相流动的能量和组分方程ijjpipipipprpixvxvxiijji vvx-vnv-vnm+ g+)(+xp-=vvxvtjiijj hvx-QW+Qn+q-xTx=hvx+htjjssppprjjjjjsjSjSjSjjsvYx-W-xYDx=Yvx+Yt0=vxtjj0=vnx+tnpjpjp pjppipippjpipjpjpipkiiprippipjpjpipvnv+vnv+vvnx- vn-vvn1+gn=vvnx+vnt pjpppppjppjpjrpphppppjpjpp v nh+ h nv+ h vnx- Q-Q-Qmn=hvnx+hnt rpdtduuuu