Журналы. GSS системы.

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

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

ENTRY GRID

HEADER=POTENTIAL FLOW OVER A BLOCK PRINTF=PRINT14 PL0TF=PL0T14

CALL INTA2(NCVLX,12,NCVLY,10) CALL DATA2(XL,9.,YL,4.) CALL EZGRID RETURN

ENTRY BEGIN TITLE(1)= VEL. POT. TITLE(2)= и VELOCITY TITLE(3)= V VELOCITY

CALL INTA8(KSOLVE(1),1,KPRINT(1),1,KPLOT (1),1, 1 KPRINT (2) ,1, KPLOT(2),1,KPRINT (3) ,1,KPLOT(3) ,1,LAST,3) DO 100 J=1,M1 DO 100 1=1,Ll POT(I,J)=-X(I) 100 CONTINUE RETURN

ENTRY OUTPUT

CALCULATE U AND V VELOCITIES AT THEIR STAGGERED LOCATIONS DO 200 J=3,M2 DO 200 1=2,L2

V(I,J)=-(POT(I,J)-POT(I,J-1))/(Y(J)-Y(J-l)) 200 CONTINUE

DO 210 J=2,M2 DO 220 1=3,L2

и(I,J)=-(POT(I,J)-POT(I-l, J))/(X(I)-X(I-l)) 220 CONTINUE

U{2,J)=FLUXI1(J,1) и(Ll,J)=-FLUXL1(J,1) 210 CONTINUE

DO 230 J=2,M2 DO 230 1=2,L2

IF(X(I).GT.3..AND.X(I).LT.6..AND.Y(J).LT.2.) THEN U(I,J)=0. U(I+1,J)=0. V(I,J)=0. V(I,J+1)=0. ENDIF 230 CONTINUE

DO 240 IUNIT=IU1,IU2

IF(ITER.EQ.O) WRITe(IUNIT,250) 250 FORMAT(IX, ITER,2X, PHI(3,5) ,3X, U ( 3,5) ,4X,

1 •V(3,5) ,3X, •PHI(5,5) •,3X, •U(5,5) ,4X, V(5,5) )

WRITE(lUNIT,2 60) ITER,POT(3,5),U(3,5),V(3,5), 1 P0T(5,5),U(5,5),V(5,5)

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

IF(ITER.EQ.LAST) THEN CALL PRINT

COME HERE TO FILL IBLOCK(I,J) BEFORE CALLING PLOT DO 270 J=2,M2

DO 270 1=2,L2

IF(X(I).GT.3..AND.X(I).LT.6..AND.Y(J).LT.2. 1 IBLOCK(I,J)=0.

270 CONTINUE С-

CALL PLOT ENDIF

RETURN

ENTRY PHI DO 300 J=2,M2 DO 300 1=2,L2 GAM(I,J)=l.

IF(X(I) .GT.3. .AND.X(I) .LT.6. .AND.Y (J) .LT.2.)

1 GAM(I,J)=0.

300 CONTINUE

COME HERE TO SPECIFY BOUNDARY CONDITIONS DO 310 1=2,L2 KBCJl(I)=2 KBCMl (I)=2 310 CONTINUE RETURN END

С CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC

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

RESULTS OF CONDUCT FOR CARTESIAN COORDINATE SYSTEM **************************************************

POTENTIAL FLOW OVER A BLOCK

ITER PHI(3,5) U(3,5)

V(3,5) PHI (5, 5) U(5,5) V(5,5)

0 -1.13E+00 1.OOE+00 O.OOE+00 -2.63E+00 1.OOE+00 O.OOE+00

1 -6.23E-01 5.51E-01 I.22E-01 -1.28E+00 3.71E-01 4.OlE-01

2 -6.23E-01 5.51E-01 1.22E-01 -1.28E+00 3.71E-01 4.OlE-01

3 -6.23E-01 5.51E-01 I.22E-01 -1.28E+00 3.71E-01 4.OlE-01 1= I 2 3 4 5 6 7

X= 0.OOE+00 3.75E-01 1.13E+00 1.88E+00 2.63E+00 3.38E+00 4.13E+00 1= 8 9 10 11 12 13 14

X= 4.88E+00 5.63E+00 6.38E+00 7.13E+00 7.88E+00 8.63E+00 9.OOE+00

J= 1 2 3 4 5 6 7

Y= O.OOE+00 2.OOE-01 6.OOE-01 1.OOE+00 1.40E+00 1.80E+00 2.20E+00

J= 8 9 10 11 12

Y= 2.60E+00 3.OOE+00 3.40E+00 3.80E+00 4.OOE+00

****** vEL. POT. ******

1= 1 2 3 4 5

12 0.OOE+OO -2.93E-01 -8.93E-01 -1.55E+00 -2.29E+00

11 O.OOE+00 -2.93E-01 -8.91E-01 -1.54E+00 -2.28E+00

10 0.OOE+OO -2.88E-01 -8.76E-01 -1.52E+00 -2.25E+00

9 O.OOE+00 -2.78E-01 -8.46E-01 -1.46E+00 -2.19E+00

8 O.OOE+00 -2.64E-0I -8.OlE-01 -1.38E+00 -2.07E+00

7 0.OOE+OO -2.47E-01 -7.45E-01 -1.26E+00 -I.87E+00

-3.14E+00 -4.O4E+00 -3.13E+00 -4.04E+00 -3.12E+00 -4.04E+00 -3.09E+00 -4.03E+00 -3.05E+00 -4.02E+00 -3.OlE+00 -4.0IE+00

OOE+00 -

-2.28E-01

82E-01

13E+00

52E+00

77E+00

93E+00

O.OOE+00 -

-2.10E-01

23E-01

0ОЕЮ0

28E+00

61E+00

88E+00

OOE+00 -

-1.95E-01

74E-01

02E-01

12E+00

50E+00

84E+00

OOE+OO -

-1.85E-01

40E-01

37E-01

02E+00

44E+00

82E+00

OOE+OO -

-1.79E-01

23E-01

04E-01

71E-01

41E+00

81E+00

OOE+OO -

-1.78E-01

20E-01

OOE-01

65E-01

40E+00

80E+00

J 12

.96E+00

-5.86E+00

.71E+00

.45E+00

.llE+00

71E+00

OOE+00

.96E+00

-5.87E+00

72E+00

.4 6E+00

.llE+00

71E+00

OOE+00

.96E+00

-5.88E+00

75E+00

.48E+00

.12E+00

71E+00

OOE+00

.97E+00

-5.91E+00

.81E+00

.54E+00

.15E+00

72E+00

OOE+00

.98E+00

-5.95E+00

.93E+00

.62E+00

.20E+00

74E+00

OOE+00

.99E+00

-5.99E+00

.13E+00

.74E+00

.26E+00

75E+00

OOE+00

.07E+00

-6.23E+00

.48E+00

.87E+00

.32E+00

77E+00

OOE+00

.12E+00

-6.39E+00

.72E+00

.OOE+00

.38E+00

.79E+00

OOE+00

. 16E+00

-6.50E+00

.88E+00

.lOE+00

.43E+00

80E+00

OOE+00

. 18E+00

-6.56E+00

.98E+00

.16E+00

.4 6E+00

.82E+00

OOE+00

. 19E+00

-6.59E+00

.03E+00

.20E+00

.48E+00

.82E+00

OOE+00

.20E+00

-6.60E+00

.03E+00

.20E+00

.48E+00

.82E+00

OOE+00

****** и VELOCITY ******

J 12

OOE+OO

OOE+00

O.OOE+00

OOE+00

75E-01

98E-01

8.69E-01

90E-01

13E+00

21E+00

OOE+00

62E-01

84E-01

8.54E-01

83E-01

15E+00

22E+00

OOE+00

37E-01

57E-01

8.24E-01

65E-01

20E+00

25E+00

OOE+00

OlE-01

16E-01

7.71E-01

20E-01

30E+00

29E+00

OOE+00

57E-01

64E-01

6.93E-01

07E-01

52E+00

34E+00

OOE+00

lOE-01

06E-01

5.90E-01

24E-01

OOE+00

64E-01

51E-01

5.02E-01

71E-01

OOE+00

26E-01

06E-01

4.37E-01

87E-01

OOE+00

98E-01

74E-01

3.96E-01

42E-01

OOE+00

84E-01

58E-01

3.76E-01

23E-01

OOE+00

OOE+OO

OOE+00

0.OOE+00

OOE+OO

OOE+00

O.OOE+00

OOE+00

23E+00

21E+00

13E+00

9.90E-01

69E-01

98E-01

75E-01

24E+00

22E+00

15E+00

9.83E-01

54E-01

84E-01

62E-01

26E+00

25E+00

20E+00

9.65E-01

24E-01

57E-01

37E-01

28E+00

29E+00

30E+00

9.20E-01

71E-01

16E-01

OlE-01

30E+00

34E+00

52E+00

8.07E-01

93E-01

64E-01

57E-01

OOE+00

5.24E-01

90E-01

06E-01

lOE-01

OOE+00

OOE+OO

OOE+00

3.71E-01

02E-01

51E-01

64E-01

OOE+OO

OOE+00

2.87E-01

37E-01

06E-01

26E-01

OOE+00

2.42E-01

96E-01

74E-01

98E-01

OOE+00

2.23E-01

76E-01

58E-01

84E-01

OOE+00

O.OOE+00

OOE+00

♦***♦♦ V VELOCITY ******

0.OOE+00

OOE+00

O.OOE+00

OOE+00

OOE+OO

0.OOE+00

O.OOE+00

23E-02

79E-02

6.46E-02

67E-02

92E-02

07E-02

O.OOE+00

41E-02

54E-02

1.33E-01

67E-01

58E-02

95E-02

OOE+00

48E-02

llE-01

2.09E-01

94E-01

02E-01

33E-02

OOE+00

30E-02

41E-01

2.88E-01

99E-01

71E-02

78E-02

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