Разделы сайта

Читаемое

Обновления Mar-2024

Промышленность Ижоры -->  Динамика жидкости: уравнения 

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 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182

1С ,

2 С ALGEK APPLIES А HODXFtED MULTI(4)-SURFACE TECHNIQUE TO

3 С THE GENERATION OF A GRID ABOUT A NACA-00*T AEROFOIL

4 С AT ZERO INCIDENCE. THE UPPER HALF GRID IS GENERATED

6 DIMENSION X(51,51> Д(51,51) ДВ(б) ,YB(6) *XS(4r51) ,YS(4,51)

7 DIMENSION RAC<51),RFD(51),SAF<51),SCD{51),SH(4),XA(51),RA(51)

8 COMMON tA,XA

10 OPEN(1,FILE- ALGEM.DAT *)

11 OPEN(6,FILES ALGEH.OUT)

12 READ(1Д)аНАХ,KMAX,IPR,IRFLД,S2,S3,AV

13 1 FORMAT(4I5,4E10.3)

14 READ(1,2)PAC, QAC, PFD,QFD,PAF,QAF,PCD, QCD

15 2 FORMAT(8E10.3)

16 WRITE(6,3)

17 VRITE(6,4)JMAX,KMAX,IPR*T*S2rS3,AV

18 3 FORMAT(* MULTISURFACE GRID GENERATION*) 19, 4 FORMATC JMAX KHAX 42I3, IPR=tM2,

20 1 T *,El0.3,5Xr S2,S3*\2E10.3, AW=\F6.3f, >

21 . VRITE(6,5)PAC,QAC.PFD,QFD

22 VRITEt6,6)PAF,QAF,PCD,QCD

23 5 FORMATC PAC \E10,3,* QAC ,E10.3, PFD*%E10.3, QFD SE10.3)

24 6 FORMATC PAF::s\E10.3/ QAF**,E10.3r PCD \B10.3/ 0CD ,E10.3)

25 С *

26 С DEFINE CORNER POINTS OF THE BOUNDARY

27 С

28 DATA XB/2.00,3.00,5.00,5.00,2.25,0.00/

29 DATA YB/0.00,0.00,0.00,2.25,2.25,0.00/

30 PI 3.14159265

31 С

32 С GENERATE STRETCHING FUNCTIONS

33 с

34 CALL STRECH(JMAX,PAC,QAC,RAC)

35 с

36 CALL STRECH(JMAX,PFD,OFD,RFD)

37 С

38 CALL STRECH(KMAX,PAF,QAF,SAF)

39 С

40 CALL STRECH(KMAX,PCD,OCD,SCD)

41 С

42 WRITE(6,7)(RAC(J),J l,JMAX)

43 WRITE(6,8)(RFD(J),J 1,JMAX)

44 VRITE(6,9)(SAF(J),J-1,JMAX)

45 VRITE(6,10)(SCD(J),J 1,JMAX)

46 7 FORMATC RAC M8F7.4)

47 8 FORMATC RFD M8F7.4)

48 9 FORMATC SAF M8F7.4)

49 10 FORMATC SCD=M8F7.4,/)

50 С

51 С OBTAIN SURFACE COORDINATES OF BODYr AB

52 С

63 CALL FOIL(0,T,RAB,XD,YD) 54 С



106 С

Рис. 13.26 (продолжение).

55 RACMX = RAB + ХВ(3) - ХВ(2)

56 RFE = 0.5*PI*(YB(5) - YB(1))

57 RFDMX = RFE + ХВ(4) - ХВ(5)

58 IFdPR .EQ. OGOTO 13

59 VRITE(6,11)(xa(L),L=1.51)

60 VRITE(6,12)(RA(L),L=1,51)

61 11 FORMATC xa=M8F7.4)

62 12 FORMATC RA=M8F7.4)

63 С

64 С GENERATE BOUNDING SURFACES 1 AND 4

65 С

66 13 do 17 J = 1,JMAX

67 RACD = RAC{J)*RACMX

68 IF(RACD .LT. RAB)GOTO 14

69 XS(1,j) = XB{2) + (RACD - RAB)*{XB(3) - XB(2))/(RACMX-RAB)

70 ys{1,j) = 0.

71 g6t0 15

72 С

73 14 CALL F0IL(1Л,RACD,XD,УD)

74 С

75 xs(1,j) = XD*{XB(2)-XB(1)> + XB(1) , 76 ys{1,j) = yd*(XB{2)-XB(1)) +YB(1)

77 15 RFDD = RFD(j)*RFDMX

78 IF(RFDD .LT. RFE)G0T0 16

79 xs(4,j) = XB(5) + {RFDD-RFE)*{XB{4)-XB(5))/(RFDMX-RFE>

80 ys(4,j) = YB(5)

81 GOTO 17

82 16 RR = YB(5) - YB(1)

83 THE = RFDD/RR

84 xs(4,j) XB(5) - RR*COS(THE)

85 ys(4,j) = RR*SIN(THE)

86 17 CONTINUE

87 С

88 с SURCH GENERATES SURFACES 2 AND 3 SO THAT THE GRID

89 с ADJACENT TO SURFACES 1 AND 4 IS ORTHOGONAL

90 с

91 CALL SURCH(JMAX.s2,s3,XS,ys)

92 С

93 DO 21 L = 1.4

94 VRITE(6,18)L

95 18 FORMATdOH SURFACE Д1)

96 DO 21 J = 1,JMAX,18

97 JA = J

98 JB = JA + 17

99 VRITE(6,19) (XS{L JC),JC=JA,JB)

100 VRITE(6,20)(YS(L,JC),aC aA,JB)

101 19 FORMATC XS=M8F7.4)

102 20 FORMATC YS=M8F7.4,/)

103 21 CONTINUE

104 С

105 С GENERATE INTERIOR GRID



107 . А1 2./(3.*AW-l.)

108 А2= 2./(3.*(l.-AW) - 1.)

109 AJH JMAX - 1

110 DZI l./AJM

111 DO 24 К = 1,KMAX

112 DO 23 J 1,JMAX

113 AJ J - 1

114 ZI AJ*DZI

115 S SAF(K) + ZI*(SCD(K)-SAF(K))

116 SH{1> = (l.-S)**2*(l.-Al*S)

117 SH(2) (l.-S)**2*S*(Al+2.)

118 SH(3) (l.-S)*S*S*(A2+2.)

119 SH(4). S*S*(1.-A2*(1.-S))

120 X(J,K) 0.

121 Y(J,K) 0.

122 DO 22 L * 1,4

123 X(J,K) X(J,K) + SH(L)*XS(L,J)

124 Y(J,K> = Y(J,K) + SH(L)*YS(L,J)

125 22 COHTINUE

126 23 CONTINUE

127 24 CONTINUE

128 с

129 с REFLECT GRID ABOUT X-AXIS

130 С

131 IFdRFL .EQ. OIGOTO 28

132 JMAP = JMAX - 1

133 DO 27 К * 1,KHAX

134 DO 25 J 1,JMAX

135 JA = 2*JMAX - J

136 JB = JA - JMAP

137 X(JA,K) X(JB,K)

138 Y(JA,K) = y(JB,K)

139 25 CONTINUE

140 DO 26 J = 1,JMAP

141 JA 2*JMAX * J

142 X(J,K) X{JA,K)

143 Y(J,K) -Y(JA,K)

144 26 CONTINUE

145 27 CONTINUE

146 JMAX JMAX + JMAP

147 с

148 28 DO 33 К 1,KMAX

149 WRITE(6,29)K,SAF(K),SCD(K)

150 29 FORMATC K M3,5X,* SAF \E10.3/ SCD ,E10.3>

151 DO 32 J 1,JMAX,18

152 JA J

153 JB = JA + 17

154 WRITE(6,30)(X(JC,K),JCxJA,JB)

155 WRITE(6,31)(Y(JC,K),JC JA,JB>

156 30 FORMATC XxM8F7.4)

157 31 FORMATC YxM8F7.4,/)

158 32 CONTINUE

159 33 CONTINUE

160 STOP

161 END

Рис. 13.26 (окончание).



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 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182

© 2003 - 2024 Prom Izhora
При копировании текстов приветствуется обратная ссылка