|
| |
|
Главная Журналы ЛР (I,J)=0. ALAM(I,J)=1. GAM{I,J)=1. IBLOCK(I,J)=0 DO 60 N=1,NFMAX 60 F(I,J,H)=0. 50 CGNTINQE DO 70 1=2,NI KBCJl(I)=1 KBCMl(I)=1 FLXCJl(I)=0. FLXCMl(I)=0. FLXPJl(I)=0. FLXPMl(I)=0. 7 0 CONTINUE DO 80 J=2,NJ KBCIl(J)=1 KBCLl(J)=1 FLXCIl(J)=0. FLXCLl(J)=0. FLXPll(J)=0. FLXPLl(J)=0. 8 0 CONTINUE RETURN ENTRY READY IF(KOUT.NE.1) OPEN(UNIT=7,FILE=PRINTF) IU1 = 6 IF(K0UT.EQ.2) IU1=7 IU2 = 7 IF(KOUT.EQ.l) IU2=6 CREATE INITIAL OUTPUT DO 100 IUNIT=IU1,IU2 IF(MODE.EQ.l) WRITE(lUNIT,1) 1 FORMAT 1IX, RESULTS OF CONDUCT FOR CARTESIAN COORDINATE SYSTEM 1/lX,50(1H*) ) IF(M0DE.EQ.2) WRITE(lUNIT,2) 2 FORMAT(IX,RESULTS OF CONDUCT FOR AXISYMMETRIC COORDINATE SYSTEM 1/lX,53(IH*) ) IF(MODE.EQ.3) WRITE(lUNIT,3) 3 FORMAT(IX,RESULTS OF CONDUCT FOR POLAR COORDINATE SYSTEM 1/lX,46(IH*) ) WRITE(lUNIT,5) HEADER 5 FORMAT(IX,64(-)/lX,A64/lX,64(-) ) IF(L1.GT.NI.OR.Ml.GT.NJ.OR.Ll.LT.4.OR.Ml.LT.4) THEN WRITE(lUNIT,6) 6 FORMAT{IX,EXECUTION TERMINATED DUE TO ONE(OR MORE) OF THE FOLLOWI ING REASON(3)/2Х,I) Ll GREATER THAN NI/2X,2) Ml GREATER THAN NJ 2 /2X,3) Ll LESS THAN 4/2X,4) Ml LESS THAN 4/) KST0P=1 ENDIF 100 CONTINUE IF(KSTOP.NE.0) STOP CALCULATE GEOMETRICAL QUANTITIES L2=L1-1 L3=L2-1 M2=M1-1 M3=M2-1 X(1)=XU(2) DO 110 1=2,L2 110 X(I)=0.5*(XU(I+1)+XU(I)) X(L1)=XU(L1) Y{1)=YV(2) DO 120 J=2,M2 120 Y(J)=0.5*(YV(J+1)+YV(J)) Y{M1)=YV(M1) DO 130 1=2,L2 130 XCV(I)=X0(I+l)-XU(I) DO 14 0 J=2,M2 140 YCV(J)=YV(J+1)-YV(J) IF(MODE.EQ.l) THEN DO 150 J=1,M1 RV(J)=1 150 R(J)=1 ELSE RY1=R(1)-Y(1) DO 160 J=2,M1 160 R(J)=Y(J)+RY1 RV(2)=R(1) DO 170 J=3,M2 170 RV(J)=RV(J-1)+YCV(J-1) RV(M1)=R(M1) ENDIF IF(M0DE.EQ.3) THEN DO 180 J=1,M1 180 SX(J)=R(J) ELSE DO 190 J=1,M1 SX(J)=1. 190 CONTINUE ENDIF DO 200 J=2,M2 YCVR(J)=R(J)*YCV(J) IF(M0DE.EQ.3) THEN ARX (J) =YCV (J) ELSE ARX(J)=YCVR(J) ENDIF 200 CONTINUE С RETURN END CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC SUBROUTINE HEART Q + + + + * + + * + + * + * + + + $INCLUDE;COMMON CONSTRUCT LOOP FOR ALL EQUATIONS DO 999 N=I,NFMAX NF=N IF(KSOLVE(NF).EQ.O) GO TO 999 CALL PHI CALCULATE COEFFICIENTS IN THE DISCRETIZATION EQUATION С BETA=4./3. IF(KORD.EQ.I) BETA=I. RLX=(1.-RELAX(NF))/RELAX(NF) CONSIDER VOLUMETRIC TERMS DO 10 J=2,M2 DO 10 1 = 2, L2 VOL=YCVR(J)*XCV(I) APT=ALAM(I, J)/DT CON(I,J) = (CON(I,J)+APT*F(I, J, NF) )*VOL AP(I,J)=(APT-AP(I,J))*VOL 10 CONTINUE COEFFICIENTS FOR X-DIRECTION DIFFUSION DO 20 J=2,M2 DO 20 1=2,L3 DIFF=ARX(J)*2.*GAM(I, J)*GAM(I+1, J) /( (XCV(I)*GAM(I+ 1, J) + 1 XCV(I+l)*GAM(I,J)+SMALL)*SX(J)) AIP(I, J)=DIFF+SMALL AIM(I+1,J)=AIP(I,J) 20 CONTINUE DO 30 J=2,M2 CONSIDER 1=1 BOUNDARY DIFF=GAM(2,J)/(0.5*XCV(2)*SX(J))+SMALL AIM(2,J)=BETA*DIFF 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 |