Разделы сайта
Читаемое
Обновления Mar-2024
|
Промышленность Ижоры --> Динамика жидкости: уравнения 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 (окончание).
|
© 2003 - 2024 Prom Izhora
При копировании текстов приветствуется обратная ссылка |