|
| |
|
Главная Журналы AIP(1,J)=AIM(2, J) AIM(2,J)=AIM(2,J)*ARX(J) AIM(1,J)=(BETA-1.)*AIP(2,J)/ARX(J) AIP(2, J)=AIP(2, J) +AIM ( 1, J) *ARX (J) IF(KBCI1(J).EQ.1) THEN C0N(2, J)=C0N(2, J)+AIM(2, а)Т(1, J,NF) ELSE AP(1,J)=AIP(1,J)-FLXPIl(J) C0N(1,J)=FLXCI1(J) TEMP=AIM(2,J)/AP(1,J) AP(2,J)=AP(2,J)-AIP(1,J)*TEMP AIP(2,J)=AIP(2,J)-AIM(1,J)*TEMP CON (2, J) =CON (2, J) +CON (1, J) *TEMP ENDIF AP(2,J)=AP(2,J)+AIM(2,J) AIM(2,J)=0. CONSIDER I=L1 BOUNDARY DIFF=GAM(L2,J)/(0.5*XCV(L2)*SX(J))+SMALL AIP(L2,J)=BETA*DIFF AIM(L1,J)=AIP(L2,J) AIP(L2,J)=AIP(L2,J)*ARX(J) AIP(L1,J)=(BETA-1.)*AIM(L2,J)/ARX(J) AIM(L2,J)=AIM(L2,J)+AIP(Ll,J)*ARX(J) IF(KBCL1(J).EQ.1) THEN C0N(L2,J)=CON(L2,J)+AIP(L2,J)*F(L1,J,NF) ELSE AP(L1,J)=AIM(L1,J)-FLXPLl(J) C0N(L1,J)=FLXCL1(J) TEMP=AIP(L2,J)/AP(L1,J) AP(L2,J)=AP(L2,J)-AIM(L1,J)*TEMP AIM(L2,J)=AIM(L2,J)-AIP(L1,J)*TEMP CON(L2,J)=C0N(L2, J)+CON(Ll,J)*TEMP ENDIF AP(L2,J)=AP(L2,J)+AIP(L2,J) AIP(L2,J)=0. 30 CONTINUE COEFFICIENTS FOR Y-DIRECTION DIFFUSION DO 4 0 J=2,M3 DO 40 1=2,L2 AREA=RV(J+1)*XCV(I) DIFF=AREA*2.*GAM(I,J)*GAM(I,J+1)/(YCV(J)*GAM(I,J+1) 1 YCV(J+1)*GAM(I,J)+SMALL) AJP(I,J)=DIFF+SMALL AJM(I,J+1)=AJP(I,J) 40 CONTINUE DO 50 1=2,L2 CONSIDER J=I BOUNDARY AREA=RV(2)*XCV(I) DIFF=GAM(I,2)/(О.5*YCV(2))+SMALL AJM(I,2)=BETA*DIFF AJP(I,1)=AJM(I,2) AJM(1,2)=AJM(1,2)*AREA AJM (1, 1) = (BETA- 1 . ) AJP (1, 2 ) / ( RV ( 3 ) *XCV (I) ) AJP(1,2)=AJP(1,2)+AJM(1,1)*AREA IF(KBCJ1(I).EQ.1) THEN CON(I,2)=C0N(I,2)+AJM(I, 2)*F(I,1,NF) ELSE AP(I,1)=AJP(I,1)-FLXPJl(I) CON(I,1)=FLXCJ1(I) TEMP = AJM(I,2)/АР(I, 1) AP(I,2)=AP(1,2)-AJP(1,1)*TEMP AJP(1,2)=AJP(1,2)-AJM(1,1)*TEMP CON(1,2)=CON(1,2)+CON(1,1)*TEMP ENDIF AP(I,2)=AP(1,2)+AJM(1,2) AJM(I,2)=0. CONSIDER J=M1 BOUNDARY AREA=RV(M1)*XCV(I) DIFF=GAM(I,M2)/(0.5*YCV(M2))+SMALL AJP(I,M2)=BETA*DIFF AJM(I,M1)=AJP(I,M2) AJP(I,M2)=AJP(I, M2)*AREA AJP(I,MI)=(BETA-1.)*AJM(I,M2)/(RV(M2)*XCV(I)) AJM(I,M2)=AJM(I,M2)+AJP(I,Ml)*AREA IF(KBCM1(I).EQ.l) THEN CON(I,M2)=CON(I,M2)+AJP(I,M2)*F(I,M1,NF) ELSE AP(I,M1)=AJM(I,Ml)-FLXPMl(I) C0N(I,M1)=FLXCM1(I) TEMP=AJP(I,M2)/АР(I,M1) AP(I,M2)=AP(I,M2)-AJM(I,Ml)*TEMP AJM(I,M2)=AJM(I,M2)-AJP(I,Ml)*TEMP CON(I,M2)=CON(I,M2)+CON(I,Ml)*TEMP ENDIF AP(I,M2)=AP(I,M2)+AJP(I,M2) AJP(I,M2)=0. 50 CONTINUE COME HERE TO INTRODUCE UNDERRELAXATION CONSTRUCT AP(I,J) AND CON(I,J) IN THEIR FINAL FORM DO 60 J=2,M2 DO 60 1=2,L2 ANB=AIP(I,J)+AIM(I,J)+AJP(I,J)+AJM(I,J) AINR=ANB*RLX AP(I,J)=AP(I,J)+ANB+AINR CON(I,J)=CON(I,J)+AINR*F(I,J,NF) 60 CONTINUE CALL THE SOLVE ROUTINE TO OBTAIN THE SOLUTION OF THE С DISCRETIZATION EQUATIONS CALL SOLVE 999 CONTINUE С TIME=TIME+DT ITER=ITER+1 IF(ITER.GE.LAST) KSTOP=I RETURN С cccccccccccccccccccccccccccccccccccccccccccccccccc SUBROUTINE SOLVE Q + * + + + + * + * + + + * + * + + + $INCLUDE:COMMON DIMENSION RT(6) Q + + * + + + + + * + * + + * + + + * + * + + + + * + * + + + + + BIG1=1.E+10 SMALLl=1.0E-5 LL2 = 2-*L2 LL=LL2-2 MM2=2*M2 MM=MM2-2 N=NF NTM=NTIMES(N) DO 999 NT=1,NTM NTT=NT ICON=I COME HERE TO PERFORM THE I-DIRECTION BLOCK CORRECTION С--- PTX(1)=0. QTX(1)=0. DO 10 1=2,L2 BL=SMALL BLP=0 . BLM=0. BLC=0. DO 20 J=2,M2 IF(AP(I,J).LT.BIGl) THEN BL=BL+AP(I,J) IF(AP(I,J+1).LT.BIGl) BL=BL-AJP(I,J) IF(AP(I,J-1).LT.BIGl) BL=BL-AJM(I,J) IF(AP(I+1,J).LT.BIGl) BLP=BLP+AIP(I,J) IF(AP(I-1,J).LT.BIGl) BLM=BLM+AIM(I,J) CONVERGENCE CRITERION FOR THE SOLUTION ROUTINE RT(1)=AIP(I,J)*F(I+1,J,N) RT(2)=AIM(I,J)*F(I-1, J,N) 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 |