《工程中的概率方法》Section3Article1

人 4人∞。S P灯 w o.ffeenhalE t2. co 以或02+Dm=2 人x +古.7千=k72< x、4aa0 docit fiel0 Ul,y,ti+v(xy. =80( L oCOe T alsace 木ac!小 +可VT=○

个个 ailton( apf 儿mA DE4a边Qok Dh从以对.平n+w4从 co、Ua0m2a c心Q时成以人,,头 L-2 J-+ 下E&x、sy 久&L )<0+点Q以 △X k一Damd —人 ⊥

从x qS2、S ap奶m+b2bw Q. abortional△t +G(△x 士△X △ △ +△ + 人 +○(△

Oh stand T ⑦(△x 丁 △ + △ + △ △ Da com Wus cor e tvt arr 飞a(m:力人m0(a PDE 过9+tca OOE Mothed to untaa tho cape& is Souse utm e fo x0e、m,mAE1 △t +,=人区

Q Court a以大ms 来 he dhoo he △dL以又 yy人 。ae Q 三XQ2:e 大 s+1八 sHlet o 文。 Tpl=5a°K 0大C丈)大 viQ waa a paulos 人C人义 (mona that aA△ at tee owet

Mon feb2510:27:412002 1 8 Set-up Dirichlet conditions at bottom wa kx0=11·(j-1}·(Nx+1 r(k01=74i,j}-Ta11 A(kk0.kk0}=1.0 for ii = 1: Nx+1 for 33= 1: Ny+l 8 Set-up Dirichlet conditions at top boundary (11,j1} 1)*dx; ii 1: Nx-1 kkO= ii +(11-11 r(kkD T(ii, 13)-Tinlet 8 Set-up matrix to solve for convection-diffusion equation g Set-up parabolized boundary conditions at outlet kO (j-1}(Nx+1 6 Pirst, discretize the interior points kk=ii-1+(j-1}* kkN= i1 +(31 )*(Nx+ 2}*(Mx+1 (KkD ulii,]3)/dx"(T(i1 T(1-1,j 11,11)/d·(T41 k/dy^2·{T(i,5+1)-2*1,)·T(i1,j-1 (-2}*(Nx+1} A(kkD, kkOl u(i1, 33)/dx +2*k/dy2 k0)=0.5+u11,j/dx*T(11+1,1)·T(11-1,j A(Kk0,kNN)=0.5v(11,1}/dy 5*v(11,j)/dy*t 5v(i, k/dy^2·{T ) T(i Solve for dm * kdx k/dx"2 为 Update T for j]=l: Ny (-1}·(Mx+1); Set-up Dirichlet conditions at inlet 5,j}=T,j》+d(kk0) kk0= ii +(11-15(NX+1) g Plot result r(kkD ) Tlli, jj)-Tinlet:

2x2C∠ k=1 k=0.1 0.8 0.6 0.2 0.2 0.4 0.6 0.8 0.2 0.4 0.60.8 k=001 k=0.001 0.8 08 0.6 06 04 0.2 0.2 0.4 0.6 020.40.60.8

k=0001,20X20grd 0.6 0.4 0.2 05 k=0001,20X60grd 08 0.6 0.4 05

20by20grd,k=0.001 02r 0.1 600 700 800 9001000110012001300140015001600 20by60grd,k=0.001 .2 >0.1 600 08009001000110012001300140015001600