Главная  Журналы 

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

ASUM=0. WSUM=0. TSOM=0. DO 200 J=2,M2 DO 200 1=2,L2

AR=XCV(I)*YCV(J)

ASUM=ASUM+AR

WSUM=WSUM+W(I,J)*AR

TSUM=TSOM+W(I,J)*T(I,J)*AR 200 CONTINUE

WBAR=WSUM/ASUM TB=TSUM/(WSUM+SMALL) WP=2.*X(L1)+Y(M1) DH=4.*ASUM/WP RE=DH*WBAR*DEN/AMU

FRE=-2.*DPDZ*DH/(DEN*WBAR**2+SMALL)*RE

TWAV=0.

DO 210 1=2,L2

TWAV=TWAV+XCV(I)*T(I,1) 210 CONTINUE

TWAV=TWAV/X(LI)

ANU=QW*DH/(COND*(TWAV-TB)+SMALL) DO 220 IUNIT=IU1,IU2

IF(ITER.EQ.O) WRITE(IUNIT,230) 2 30 FORMAT(2X,ITER,4X,W(5,8),4X,W(3,7),

1 4X,T(5,8),4X,T(3,7),5X,•FRE,8X,-NU)

WRITE(lUNIT,24 0) ITER,W(5,8),W(3,7),T(5,8), 1 T (3,7) ,FRE,ANU

240 FORMAT(3X,I2,2X,1P6E10.2) 220 CONTINUE

IF(ITER.EQ.LAST) THEN COME HERE TO CALCULATE LOCAL NU ON BOTTOM WALL DO 250 1=2,L2

ANULOC=QW*DH/(COND*(T(I,1)-ТВ)+SMALL) DO 260 IUNIT=IU1,IU2

IF(I.EQ.2) WRITE(IUNIT,270) 270 FORMAT( , I ,8X, X(I) , 7X,

1 -LOCAL NU (BOTTOM WALL)•)

WRITE(lUNIT,280) I,X(I),ANULOC 280 FORMAT(IX,12,5X,1PE9.2, IIX, 1PE9.2)

260 CONTINUE

250 CONTINUE

DO 290 J=1,M1 DO 290 1=1,LI

F(I,J,3)=W(I,J)/WBAR F(I,J,4)=(T(I,J)-TWAV)/(TB-TWAV) 290 CONTINUE

CALL PRINT CALL PLOT ENDIF RETURN

С CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC

ENTRY PHI IF(NF.EQ.l) THEN

DO 300 J=2,M2

DO 300 1=2,L2



GAM(I,J)=AMO SC(I,J)=-DPDZ 300 CONTINUE ENDIF

IF(NF.EQ.2) THEN

DTDZ=QW*X(L1)/(WSUM*RHOCP) DO 320 J=2,M2 DO 320 1=2,L2 GAM(I,J)=COND SC(I,J)=-RHOCP*DTDZ*W(I,J) 320 CONTINUE ENDIF

COME HERE TO SPECIFY BOUNDARY CONDITIONS DO 310 J=2,M2 KBCLl(J)=2 310 CONTINUE

IF(NF.EQ.2) THEN DO 330 J=2,M2 KBCIl(J)=2 3 30 CONTINUE

DO 340 1=2,L2 KBCJl(I)=2 FLXCJl(I)=QW KBCMl(I)=2 34 0 CONTINUE ENDIF RETURN END

С CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC

10.1.5. Результаты расчетов

RESULTS OF CONDUCT FOR CARTESIAN COORDINATE SYSTEM

RECTANGULAR DUCT WITH ONE HEATED WALL

ITER

W(5,8)

W(3,7)

T(5,8) T(3,7)

. OOE+00

. OOE+00

0.OOE+00 O.OOE+00

OOE+00

.33E+20

.46E-02

,40E-02

O.OOE+00 O.OOE+00

llE+01

.33E+20

.46E-02

,40E-02

0.OOE+00 O.OOE+00

llE+01

.33E+20

,46E-02

,40E-02

O.OOE+00 O.OOE+00

llE+01

.33E+20

,46E-02

.40E-02

4.41E-02 1.43E-01

llE+01

.49E+00

. 46E-02

.40E-02

4.41E-02 1.43E-01

llE+01

.49E+00

.4 6E-02

, 40E-02

4.41E-02 1.43E-01

llE+01

.49E+00

X(I)

LOCAL NU (BOTTOM WALL)

.25E-02

3.14E+00

,39E-01

3.24E+00

.37E-01

3.45E+00

,53E-01

3.68E+00

,83E-01

3.84E+00

3 4

X= 0

.OOE+00 7

.25E-02 2

.39E-01 4.37E-01 6

.53E-01 8

3E-01 1



J= 1 2 3 4 5 6 7

Y= O.OOE+00 3.62E-02 1.19E-01 2.19E-01 3.27E-01 4.41E-01 5.59E-01 J= a 9 10 11 12

Y= 6.73E-01 7.81E-01 8.81E-01 9.64E-01 1.OOE+OO

****** AXIAL VELOCITY W ******

J 12

OOE+OO

OOE+OO

OOE+00

OOE+00

OOE+OO

OOE+OO

OOE+00

OOE+00

77E-03

03E-02

36E-02

53E-02

61E-02

62E-02

OOE+OO

25E-02

99E-02

05E-02

61E-02

85E-02

88E-02

OOE+OO

84E-02

65E-02

44E-02

41E-02

82E-02

88E-02

OOE+00

22E-02

81E-02

17E-02

46E-02

OOE-01

OlE-01

OOE+OO

41E-02

40E-02

09E-02

06E-01

12E-01

13E-01

OOE+OO

41E-02

40E-02

09E-02

06E-01

12E-01

13E-01

OOE+00

22E-02

81E-02

17E-02

46E-02

OOE-01

OlE-Ol

OOE+00

84E-02

65E-02

44E-02

41E-02

82E-02

88E-02

OOE+OO

25E-02

99E-02

05E-02

61E-02

85E-02

88E-02

OOE+OO

77E-03

1.03E-02

36E-02

53E-02

61E-02

62E-02

OOE+00

OOE+OO

OOE+00

OOE+00

OOE+00

OOE+00

OOE+00

TEMPERATURE

OOE+00

18E-02

6.89E-02

47E-

OlE-

02

59E-03

OOE+OO

37E-02

23E-02

6.92E-02

49E-

OlE-

58E-03

59E-03

76E-02

60E-02

7.17E-02

60E-

02

04E-

02

50E-03

46E-C3

99E-02

81E-02

8.15E-02

35E-

02

62E-

43E-03

29E-03

24E-01

22E-01

1.04E-01

32E-

02

41E-

65E-02

42E-02

65E-01

63E-01

1.43E-01

llE-

01

06E-

24E-02

OOE-02

24E-01

22E-01

2.02E-01

70E-

01

39E-

21E-01

19E-01

98E-01

96E-01

2.77E-01

47E-

17E-

OOE-01

97E-01

81E-01

79E-01

3.63E-01

35E-

01

07E-

91E-01

89E-01

68E-01

67E-01

4.52E-01

27E-

01

OlE-

85E-01

83E-01

47E-01

46E-01

5.33E-01

09E-

84E-

68E-01

66E-01

OOE+00

82E-01

5.69E-01

45E-

01

20E-

05E-01

OOE+OO

****** w/WBAR ******

1 = J

OOE+00

OOE+00

OOE+00

OOE+00

0.OOE+00

OOE+OO

OOE+OO

OOE+00

20E-02

78E-01

34E-01

2.63E-01

76E-01

78E-01

OOE+00

16E-01

15E-01

96E-01

7.93E-01

34E-01

39E-01

OOE+OO

16E-01

OOE-01

llE+00

1.27E+00

35E+00

35E+00

OOE+00

82E-01

98E-01

40E+00

1.63E+00

72E+00

74E+00

OOE+00

15E-01

lOE+00

56E+00

1.82E+00

93E+00

94E+00

OOE+OO

15E-01

lOE+00

56E+00

1.82E+00

93E+00

94E+00

OOE+00

82E-01

98E-01

40E+00

1.63E+00

72E+00

74E+00

OOE+OO

16E-01

OOE-01

llE+00

1.27E+00

35E+00

35E+00

OOE+00

16E-01

15E-01

96E-01

7.93E-01

34E-01

39E-01

OOE+00

20E-02

78E-01

34E-01

2.63E-0I

76E-01

78E-01

OOE+00

OOE+OO

OOE+OO

OOE+00

O.OOE+00

OOE+00

OOE+00





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