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

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

HEADER=NON-NEWTONIAN FLOW IN A SEMICIRCULAR DUCT PRINTF=PRINTll PL0TF=PL0T11

CALL DATA2(R(1),0.,PI,3.14159)

CALL DATA2(XL,0.5*PI,YL,1.)

CALL INTA3(MODE,3,NCVLX,10,NCVLY, 16)

CALL EZGRID

RETURN

ENTRY BEGIN

TITLE(1)= W/WBAR

TITLE(2) = (T-TWAV) / (TB-TWAV)

CALL INTA6(KSOLVE(1),1,KPRINT(1),1,KPRINT(2),1, 1 KPLOT(1),1,KPL0T(2),1,LAST, 18)

CALL DATA9(DEN,1.,REGAM,0.8,RK, 1. , POWER, 0.5, 1 AMU,1.,COND,1.,CP,1.,DPDZ,-l.,QW,l.)

RHOCP=DEN*CP

С- SINCE THE ZERO DEFAULT VALUES OF W(I,J) AND T(I,J) ARE C- SATISFACTORY,

С THESE ARRAY ARE NOT FILLED HERE. С

RETURN

Q * * * * * * *- * * * -A- -Jr * * * *- *- * * * *- *- * * * -A- * *

ENTRY OUTPUT IFdTER.EQ. 15) THEN

KSOLVE(1)=0

KSOLVE(2)=1 ENDIF ASUM=0. WSUM=0. TSUM=0. DO 200 J=2,M2 DO 200 1=2,L2

AR=XCV(I)*YCVR(J)

ASUM=ASUM+AR

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

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

WBAR=WSUM/ASUM TB=TSUM/(WSUM+SMALL) WP=R(M1)+0.5*PI*R(M1) HP=0.5*PI*R(M1) DH=4.*ASUM/WP

AMUR=RK*(WBAR/DH+SMALL)*♦(POWER-1.) RE=DEN*WBAR*DH/AMUR

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

DTDZ=QW*HP/(RHOCP*WSUM+SMALL)

TWAV=0.

DO 220 1=2,L2

TWAV=TWAV+T(I,Ml)*XCV(I) 220 CONTINUE

TWAV=TWAV/XL

ANU=QW*DH/(COND*(TWAV-TB)+SMALL)



DO 240 IUNIT=IU1,IU2

IF(ITER.EQ.0) WRITE(lUNIT,250) 250 FORMAT(IX,ITER,2X,W(3,4),5X,F.RE,5X, 1 •T(5,4) ЧбХ, NU)

WRITE(lUNIT,2 60) ITER,W(3,4),FRE,T(5,4),ANU 260 FORMAT(2X,I2,1P4E10.2) 240 CONTINUE

IF(ITER.EQ.LAST) THEN DO 270 J=1,M1 DO 270 1=1,LI

W(I,J)=W(I,J)/WBAR T(I,J)=(T(I,J)-TWAV)/(TB-TWAV) 270 CONTINUE

CALL PRINT CALL PLOT ENDIF RETURN

ENTRY PHI IF(NF.EQ.l) THEN REL=1.-REGAM DO 300 J=2,M2 DO 300 1=2,L2

IF(ITER.EQ.0) THEN

GAM(I,J)=AMU ELSE

WP=0.5*(W(I+1,J)+W(I,J)) WM=0.5*(W(I,J)+W(I-1,J)) IF(I.EQ.2) WM=W(I-1,J) IF(I.EQ.L2) WP=0. DWDX=(WP-WM)/(XCV(I)*SX(J)) WP=0.5*(W(I,J+1)+W(I,J)) WM=0.5*(W(I,J)+W(I,J-1)) IF(J.EQ.2) WM=0. IF(J.EQ.M2) WP=0. DWDY=(WP-WM)/YCV(J)

GAMT=RK*(DWDX**2+DWDY**2)**((POWER-1.)*0.5) IF(ITER.EQ.1) GAM(I,J)=GAMT GAM(I,J)=REGAM*GAMT+REL*GAM(I, J) ENDIF 300 CONTINUE

COME HERE TO ADJUST THE PRESSURE GRADIENT IF(ITER.EQ.1) THEN GAMAV=0. DO 310 J=2,M2 DO 310 1=2,L2

GAMAV=GAMAV+GAM(I,J)*XCV(I)*YCVR(J) 310 CONTINUE

GAMAV=GAMAV/ASUM DPDZ=DPDZ*GAMAV/AMU ENDIF

DO 320 J=2,M2 DO 320 1=2,L2 SC (I,J)=-DPDZ



320 CONTINUE ENDIF

IF(NF.EQ.2) THEN DO 330 J=2,M2 DO 330 1=2,L2 GAM(I,J)=COND SC(I,J)=-RHOCP*W(I, J)*DTDZ 330 CONTINUE ENDIF

COME HERE TO SPECIFY BOUNDARY CONDITIONS DO 34 0 J=2,M2 KBCIl(J)=2 340 CONTINUE

DO 350 1=2,L2 KBCJl(I)=2 350 CONTINUE

IF(NF.EQ.2) THEN DO 3 60 J=2,M2 KBCLl(J)=2 360 CONTINUE

DO 370 1=2,L2 KBCMl(I)=2 FLXCMl(I)=QW 370 CONTINUE ENDIF RETURN END

С CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC

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

RESULTS OF CONDUCT FOR POLAR COORDINATE SYSTEM **********************************************

NON-NEWTONIAN FLOW IN A SEMICIRCULAR DUCT

ITER

W(3, 4)

F.RE

T(5,4)

OOE+OO

.OOE+OO

.OOE+OO

.22E+20

29E-02

.24E+01

.OOE+OO

.22E+20

23E-02

.76E+01

.OOE+OO

.22E+20

,09E-02

.64E+01

. OOE+OO

.22E+20

,65E-02

.57E+01

.OOE+OO

.22E+20

OOE-02

.52E+01

.OOE+OO

.22E+20

.21E-02

.50E+01

.OOE+OO

.22E+20

, 34E-02

.48E+01

,OOE+OO

.22E+20

,42E-02

.47E+01

.OOE+OO

.22E+20

47E-02

.47E+01

.OOE+OO

.22E+20

. 49E-02

.46E+01

.OOE+00

.22E+20

.51E-02

.46E+01

.OOE+OO

.22E+20

,52E-02

.46E+01

.OOE+OO

.22E+20

,53E-02

.46E+01

.OOE+OO

.22E+20

.53E-02

.46E+01

.OOE+OO

.22E+20

,53E-02

.46E+01

.OOE+OO

.22E+20

,53E-02

.46E+01

.97E-01

. 33E+00

,53E-02

.46E+01

.97E-01

. 33E+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