Главная Журналы TITLE(2)= (T-T2)/(T1-T2) CALL INTA6(KSOLVE(2),1,KPRINT(1),1,KPRINT(2) ,1, 1 KPLOT (1) , 1,KPLOT (2) , 1, LAST, 6) CALL DATA8(AMU,1.,BMU,0.04,COND,1.,CP,1. , DEN, 1., 1 DPDZ,- l.,Tl,100.,-l T2,0.) RHOCP=DEN*CP DO 100 J=1,M1 DO 100 1=1,Ll W(I,J)=0 . T(I, J)=T1 T(I,M1)=T2 100 CONTINUE RETURN q - - - - - ENTRY OUTPUT IF(ITER.EQ.3) THEN KSOLVE(1)=1 KSOLVE(2)=0 ENDIF ASUM=0. WSUM=0. DO 200 J=2,M2 DO 200 1=2,L2 AR=XCV(I)*YCVR(J) ASUM=ASUM+AR WSUM=WSUM+W(I,J)*AR 200 CONTINUE WBAR=WSUM/ASUM WP=XL* (ROUT+RIN)+YL DH=4.*ASUM/WP RE=DH*WBAR*DEN/AMU FRE=-2.*DPDZ*DH/(DEN*WBAR**2+SMALL)*RE DO 210 IUNIT=IU1,IU2 IFdTER.EQ. 0) WRITE (lUNIT, 220) 22 0 FORMAT(/,IX, •ITER,2X, •W(6,8) ,5X, W(4,ll) , 1 4X, T (6,8) •,5X, T (4,11) ,6X, FRE) WRITE (lUNIT, 230) ITER,W(6,8),W(4,11),T(6,8), 1 T( 4,11),FRE 230 F0RMAT(2X,I2,1P6E11.3) 210 CONTINUE IF(ITER.EQ.LAST) THEN DO 24 0 J=1,M1 DO 240 1=1,Ll W(I,J)=W(I,J)/WBAR T(I,J)=(T(I,J)-T2)/(T1-T2) 24 0 CONTINUE CALL PRINT CALL PLOT ENDIF RETURN ENTRY PHI IF(NF.EQ.l) THEN DO 300 J=2,M2 DO 300 1=2,L2 GAM(I,J)=AMU+BMU*T(I,J) SC(I,J)=-DPDZ 300 CONTINUE ENDIF IF(NF.EQ.2) THEN DO 310 J=2,M2 DO 310 1=2,L2 GAM(I,J)=COND 310 CONTINUE ENDIF COME HERE TO SPECIFY BOUNDARY CONDITIONS DO 32 0 J=2,M2 KBCLl(J)=2 320 CONTINUE RETURN С CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 11.2.5. Результаты расчетов RESULTS OF CONDUCT FOR POLAR COORDINATE SYSTEM ********************************************** DUCT FLOW WITH TEMPERATURE-DEPENDENT VISCOSITY ITER 0 1 2 3 4 5 6 1 = TH= W(6,8) W(4,ll) O.OOOE+00 O.OOOE+00 1. O.OOOE+00 O.OOOE+00 5. O.OOOE+00 0.OOOE+00 5. 0.OOOE+00 1.106E-02 I.106E-02 1.I06E-02 I O.OOE+00 5 T(6, 8) OOOE+02 513E+01 513E+01 T(4, 11) ,000E+02 , 108E+01 FRE OOOE+00 ,00OE+00 OOOE+00 5. 614E-03 5. 6I4E-03 5. 614E-03 5. 2 3 108E+01 0.OOOE+00 513E+01 1.108E+01 O.OOOE+00 108E+01 1.850E+02 108E+01 1.850E+02 5I3E+01 513E+01 513E+01 108E+01 1.850E+02 24E-02 1.57E-01 2.62E-0I 3.67E-01 4.71E-01 5.24E-01
11.2.6. Обсуждение результатов Решения линейных уравнений для Г и w сходятся за одну итерацию каждое. Вычисленное значение JRe ~ 185, что намного больше значения /Re = 58 (см. пример 9). Основная причина такого отличия заключается в том, что при определении числа Рейнольдса мы использовали значение вязкости, равное единице, хотя она меняется в расчетной области от 1 до 5. Если же при вычислении числа Рейнольдса использовать а = 3, то /Re будет равно 61,7, что гораздо ближе к результату, найденному в примере 9. Основное влияние переменной вязкости проявляется в распределении безразмерной скорости w/w. Скорости намного выше у внешней стенки, где вязкость меньше. Положение максимума скорости сместилось в сторону внешней стенки. Скорости вблизи внутренней стенки, напротив, уменьшились. 11.2.7. Заключительные замечания Этот пример имеет много особенностей. Здесь мы впервые столкнулись с ситуацией, когда поле скорости зависит от распределения температуры, а поле температуры не зависит от скоростей. В итоге в уравнении для температуры отсутствует источниковый член. По результатам этого примера можно изучить влияние переменной вязкости на поле скорости. Переменная вязкость другого рода возникает при турбулентном течении, которое мы рассмотрим в следующем примере. 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 [ 79 ] 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 |