Di,MiVi,M1-Vi,m2 = k = constui+1M1Vi,M2Vi,M = Vi,M2(1+ k)= f Vi,M2f is determined according to total mass conservation :L2L2ZPi,miVimAx,=)VViM2△x, = FLOWINi=2i=2FLOWINVi,M1 =.L2Vi.M2ZpimiVi,m2Ax,It is taken as the boundary condition for next iteration11/35
11/35 , 1 , 2 , 2 i M i m i M v v k const v − = = , 1 , 2 2 , (1 ) i M i M i M v = + = v k f v f is determined according to total mass conservation : 2 2 , 1 , 1 , 1 , 2 2 2 L L i M i M i i M i M i i i v x v x FL f OWIN = = = = 2 , 1 , 2 2 L i M i M i FLOWIN f v x = , 1 * i M i M, 2 v = • f v It is taken as the boundary condition for next iteration
t.M(2) The 1st derivatives at outlet =const.Wi,Mui+IMIVi,M1 - Vi,M2 = k = constAyVi,MI = Vi,M2 + kAy = Vi,M2 + CC is determined according to total mass conservationZp,m(v.M2 +C)Ax, = FLOWINi=2FLOWIN-EPi.M'i.M2AxCZpi,mAx,'rMi = vi.M2 + C is taking as boundary condition for nextiteration.When flow is fully developed at outlet,: f=1,C-0;Otherwisethereissomedifferencesbetweenthetwotreatments, but such difference is not important.12/35
12/35 (2) The 1st derivatives at outlet =const. i M i M , 1 , 2 v v k const y − = = i M i M i M , 1 , 2 , 2 v v k y v C = + = + C is determined according to total mass conservation. 2 , 1 , 2 2 ( ) L i M i M i i v C x FLOWIN = + = , 1 , 2 , 1 i M i M i i M i FLOWIN v x C x − = When flow is fully developed at outlet,：f =1,C=0； Otherwise there is some differences between the two treatments, but such difference is not important. * i M, 1 i M, 2 v = + v C is taking as boundary condition for next iteration