|
| |
|
Главная Журналы TAIR - температура воздуха Г, [см. (8.9)]; TGR - температура грунта [см. (8.10а)]; TW - температура стенки [см. (8.10а)]. 8.6.4. Листинг подпрограммы ADAPT ссссссссссссссссссссссссссссссссссссссссссссссссссссссссс SUBROUTINE ADAPT с--EXAMPLE 6 - UNSTEADY CONDUCTION AROUND AN С-UNDERGROUND BUILDING С cccccccccccccccccccccccccccccccccccccccccccccccccccccccc $INCLUDE:COMMON DIMENSION T(NI,NJ) EQUIVALENCE (F(1,1,1),T(1,1)) ENTRY GRID HEADER=UNSTEADY CONDUCTION AROUND AN UNDERGROUND 1 BUILDING PRINTF=PRINTS PI=3.14159 CALL INTA3(NZX,2,NCVX(1),9,NCVX(2),3) CALL DATA2(XZONE(1),8.,XZONE(2),2.) CALL INTA3(NZY,2,NCVY(1) , 9,NCVY(2) ,3) CALL DATA2(YZONE(1),8.,YZONE(2),2.) CALL DATA2(POWRX(1),-1.5,POWRY(1),-1.5) CALL ZGRID RETURN ENTRY BEGIN TITLE(1)= TEMPERATURE CALL INTA2(KSOLVE(1),1,LAST, 120) CALL DATA7(COND,2.,RHOCP,3.,HE,5.,P,12.,DT,0.5, 1 TGR,50.,TW,70.) DO 100 J=1,M1 DO 100 1=1,Ll T (I,J)=TGR 100 CONTINUE DO 110 1=2,L2 IF(X(I).GT.8.) T(I,M1)=TW 110 CONTINUE RETURN ENTRY OUTPUT DO 200 IUNIT=IU1,IU2 IF(ITER.EQ.0) WRITE(lUNIT,210) 210 FORMAT(IX, ITER,3X, TIME,5X, T(2,14) ,3X, 1 T(2,13),3X,T(2,12),3X,T(2,ll),3X,T(2,10)) IF(MOD(ITER,3).EQ.0) THEN WRITE(lUNIT,240) ITER,TIME,T(2,14),T(2,13),T(2,12), 1 Т(2,11),Т(2,10) 240 FORMAT(IX,13,1Р6Е10.2) IF(MOD(ITER,24).EQ.O) WRITE(lUNIT,250) 250 FORMATdX, 69(*)) ENDIF 200 CONTINUE RETURN ENTRY PHI DO 300 J=2,M2 DO 300 1 = 2, L2 ALAM(I,J)=RHOCP GAM(I,J)=COND IF(X(I).GT.8.0.AND.Y(J).GT.8.0) GAM(I,J)=BIG 300 CONTINUE COME HERE TO SPECIFY BOUNDARY CONDITIONS DO 310 J=2,M2 KBCIl(J)=2 KBCLl(J)=2 310 CONTINUE ARGU=2.*PI*(TIME+DT)/Р TAIR=50.+50.*SIN(ARGU) DO 320 1=2,L2 IF(X(I).LT.8.) THEN KBCMl(I)=2 FLXCMl(I)=HE*TAIR FLXPMl(I)=-HE ENDIF 320 CONTINUE RETURN END С CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 8.6.5. Результаты расчетов RESULTS OF CONDUCT FOR CARTESIAN COORDINATE SYSTEM •k-k-k-k-kit****k*k-kit*********-k-k-k-k****-k-k*-k-k****-k-k*-k-k**** UNSTEADY CONDUCTION AROUND AN UNDERGROUND BUILDING ITER TIME T(2,14) T(2,13) T(2,12) T(2,ll) T(2,10) 0 0.OOE+OO 5.00E+01 5.OOE+01 5.00E+01 5.OOE+01 5.OOE+01 ***************************************************** 3 1.50E+00 7.47E+01 6.68E+01 5.77E+01 5.34E+01 5.19E+01 6 3.OOE+00 8.96E+01 8.15E+01 6.88E+01 6.06E+01 5.70E+01 9 4.50E+00 8.28E+01 8.04E+01 7.29E+01 6.58E+01 6.17E+01 12 6.OOE+00 5.77E+01 6.29E+01 6.61E+01 6.46E+01 6.24E+01 15 7.50E+00 2.85E+01 3.88E+01 5.14E+01 5.70E+01 5.80E+01 18 9.OOE+00 1.24E+01 2.20E+01 3.72E+01 4.68E+01 5.06E+01 21 1.05E+01 1.85E+01 2.21E+01 3.14E+01 3.99E+01 4.43E+01 24 1.20E+01 4.34E+01 3.90E+01 3.74E+01 3.99E+01 4.26E+01 ic****-k-k***-k-k*-k-k-k******-k-k-k*-**-k***-k-k-k**-k**-k*-k-k*-k******-k ITER TIME T(2,14) T(2,13) T(2,12) T(2,ll) T(2,10) 0 O.OOE+00 5.OOE+01 5.OOE+01 5.OOE+01 5.OOE+01 5.OOE+01 27 1.35Е+01 30 1.50Е+01 33 1.65Е+01 36 1.80Е+01 39 1.95Е+01 42 2.10Е+01 45 2.25Е+01 48 2.40Е+01 7.23Е+01 8.84Е+01 8.21Е+01 5.72Е+01 2.83Е+01 1.22Е+01 1.84Е+01 4.33Е+01 6.28Е+01 7.94Е+01 7.91Е+01 6.21Е+01 3.83Е+01 2.16Е+01 2.19Е+01 3.89Е+01 5.14Е+01 6.53Е+01 7.08Е+01 6.47Е+01 5.05Е+01 3.66Е+01 3.10Е+01 3.71Е+01 4.69Е+01 5.65Е+01 6.32Е+01 6.29Е+01 5.58Е+01 4.61Е+01 3.94Е+01 3.96Е+01 4.63Е+01 5.32Е+01 5.91Е+01 6.06Е+01 5.68Е+01 4.98Е+01 4.38Е+01 4.23Е+01 51 54 57 60 63 66 69 72 2.55E+01 2. 70E + 01 2.85E+01 3.OOE + 01 3. 15E + 01 3. 30E + 01 3.45E+01 3.60E+01 7.23E+01 8.84E+01 8.21E+01 5.72E+01 2.83E+01 1 . 22E+01 1. 84E+01 4 .33E+01 6.27E+01 7.93E+01 7.91E+01 6.21E+01 3.83E+01 2.17E+01 2.19E+01 3.89E+01 13E+01 52E+01 07E+01 47E+01 05E+01 66E+01 llE+01 72E+01 4.67E+01 5. 64E + 01 6.31E+01 6.29E+01 5.58E+01 4.61E+01 3.94E+01 3.97E+01 4.61E+01 5.31E+01 5.91E+01 6.07E+01 5.68E+01 4.99E+01 4.39E+01 4.23E+01 75 78 81 84 87 90 93 96 75E+01 90E+01 05E+01 20E+01 35E+01 50E+01 65E+01 80E+01 .23E+01 ,84E+01 .21E+01 .72E+01 .83E+01 .22E+01 .84E+01 .33E+01 6.27E+01 7.93E+01 91E+01 21E+01 83E+01 17E+01 19E+01 89E+01 13E+01 52E+01 08E+01 47E+01 06E+01 67E+01 llE+01 72E+01 4.67E+01 5.65E+01 6.32E+01 6.30E+01 5.59E+01 4. 62E + 01 3.95E+01 3.97E+01 62E+01 31E+01 92E+01 07E+01 69E+01 OOE+01 39E+01 24E+01 4.95E+01 5.lOE+01 5.25E+01 5.40E+01 5.55E+01 5.70E+01 5.85E+01 6.OOE+01 .23E+01 ,84E+01 ,21E+01 ,73E+01 ,83E+01 ,22E+01 ,84E+01 4.33E+01 6.27E+01 7.94E+01 7.91E+01 6.21E+01 3.83E+01 2.17E+01 2.19E+01 3.89E+01 5.13E+01 4. 6.53E+01 5. 7.08E+01 6. 6.47E+01 6. 5.06E+01 5. 3.67E+01 4. 3.11E+01 3. 3.72E+01 3. 68E+01 65E+01 32E+01 30E+01 59E+01 62E+01 95E+01 97E+01 4.62E+01 5.32E+01 5.92E+01 6.08E+01 5.70E+01 5.OOE+01 4.40E+01 4.24E+01 8.6.6. Обсуждение результатов Из анализа результатов, видно, что поле температуры становится почти полностью периодическим к пятому циклу. В периодическом распределении температуры может быть заметен эффект тепловой инерции грунта. При ITER = 99 температура поверхности почвы Т(2, 14) достаточно велика, но в глубине можно обнаружить значения температур меньше 50. Это может быть уподоблено «памяти» о предыдущей зиме, сохраненной грунтом. Точно так же в первый зимний месяц (соответствующий ITER = 1П), несмотря на то что температура Т (2,14) низкая, в глубине можно обнаружить значения выше 50. Это показывает, что, хотя температура в каждой точке меняется во времени периодически, тепловая инерция создает сдвиг по фазе между колебаниями температуры в различных точках. 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 |