Барилович В.А. Основы термогазодинамики двухфазных потоков и их численное решение
Подождите немного. Документ загружается.
360
TVS=5.677*P0**0.25+273.
CPLT1=CPIT(TVS)
DHR=CPLT1*((TL0-TVS)-TVS*ALOG(TL0/TVS))
GTR=TG*BL/RBL
PTI=GTR*DHR
CALL SPLINE (5,BAR,VPA,BVPA,CVPA,DVPA)
CALL SPLINE (5,DAVL,ENTN,BENS,CENS,DENS)
CALL SPLINE (5,BAR,TPA,BTPA,CTPA,DTPA)
CALL SPLINE (5,TEMP,ENT,BE,CE,DE)
CALL SPLINE (5,TEMP,ENTS,BES,CES,DES)
CALL SPLINE (5,TEMP,CP,BCP,CCP,DCP)
CALL SPLINE (5,BAR,VKAP,BVKAP,CVKAP,DVKAP)
CALL SEVAL
(5,P0,BAR,VPA,BVPA,CVPA,DVPA,VP,VPP,VPPP)
CALL SEVAL
(5,TK0,TEMP,VKAP,BVKAP,CVKAP,DVKAP,VK,VKP,VKPP)
CALL SEVAL
(5,TK0,TEMP,ENTS,BES,CES,DES,EPTK,EPTKP,EPTKPP)
CALL SEVAL
(5,TK0,TEMP,ENT,BE,CE,DE,ENTK,ENTKP,ENTKPP)
CALL SEVAL (5,P0,BAR,TPA,BTPA,CTPA,DTPA,TP,TPP,TPPP)
CALL SEVAL
(5,TK0,TEMP,CP,BCP,CCP,DCP,CPK,CPKP,CPKPP)
CALL SEVAL
(5,P0,DAVL,ENTN,BENS,CENS,DENS,ENTP,ENTPP,ENTPPP)
CALL SEVAL
(5,TK0,TEMP,VKAP,BVKAP,CVKAP,DVKAP,VPL,VPLP,VPLPP)
CALL SEVAL
(5,TK0,TEMP,CP,BCP,CCP,DCP,CPPL,CPPLP,CPPLPP)
CALL SEVAL (5,TK0,TEMP,ENT,BE,CE,DE,EPL,EPLP,EPLPP)
CALL SEVAL
(5,TK0,TEMP,ENTS,BES,CES,DES,EPPL,EPPLP,EPPLPP)
FID=FIP/DELF
GK=TG*BL/(RBL*PROP)
DGK=GK*(1.-X)
VGK=GK*X
PNV=6.*VK*DGK/(DK0**3*PI)
DRS=(R2-R1)/NK
VOL=DRS*FIP
WKIN=CK*COS(ALF)-U
WPIN=CV*COS(ALF)-U
DPSI=PSI(CK)
V
VI
361
VPSI=PSI(CV)
WRITE(4,*) ' Угол атаки капель =',DPSI*PF/PI
WRITE(4,*) ' Угол атаки пapa =',VPSI*PF/PI
WRITE(4,*) ' GK=',GK,'DGK=',DGK,'VGK=',VGK,'PNV=',PNV
NR=NK
PNVA=6.*(GK-
SPEED(CV)*COS(VPSI)*DRS*NK*HB/VP)/(PI*DK0**3*
& (COS(DPSI)/VK-SPEED(CV)*COS(VPSI)/SPEED(CK)/VP))
SK=PNVA/(NK*HB)
C
VGKA=SPEED(CV)*COS(VPSI)/VP*DRS*NK*HB/(1.+SPEED(CV)*
COS(VPSI)/VP
C & *(1.-X)*VK/X/SPEED(CK)/COS(DPSI))*BL/RBL
C DGKA=VGKA*(1.-X)/X
C GKA=VGKA+DGKA
C PNVA=6.*VK*DGKA/(DK0**3*PI)
WRITE(4,*) '
GK*=',GKA,'DGK*=',DGKA,'VGK*=',VGKA,'PNV*=',PNVA
C DELTAV=ABS(VGK-VGKA)/VGK*100.
C IF (DELTAV.GE.EPS) PRINT *, 'Необходима балансировка
сопла'
C IF (DELTAV.GE.EPS) GO TO 116
A0=SPEED(CV)*COS(VPSI)*(R2-R1)/ALOG(R2/R1)
DO 1 I=1,NK
R(I)=R1+(R2-R1)/NK/2.+(R2-R1)/NK*(I-1)
IF(I.EQ.NK)PRINT *,'R(NK)=',R(NK)
RWK(I)=SPEED(CK)*SIN(DPSI)
FWK(I)=SPEED(CK)*COS(DPSI)
WP(I)=A0/R(I)
TK(I)=TK0
DK(I)=DK0
CALL SEVAL
(5,P0,BAR,VPA,BVPA,CVPA,DVPA,VP,VPP,VPPP)
P(I)=P0-A0**2/VP/2.*(1./R(I)**2-1./((R2-R1)/ALOG(R2/R1))**2)
CALL SEVAL
(5,P(I),BAR,VPA,BVPA,CVPA,DVPA,VP,VPP,VPPP)
P(I)=P0-A0**2/VP/2.*(1./R(I)**2-1./((R2-R1)/ALOG(R2/R1))**2)
1 CONTINUE
WPL=0.
DELPL=0.
TPL=0.
VI
362
L=0
GPL=0
WRITE(4,100) L
WRITE(4,104)
DO 500 I=1,NK
WRITE(4,103) I,RWK(I),FWK(I),R(I)*1.E3,WP(I),P(I)*1.E-5,TK(I),
& DK(I)*1.E6
500 CONTINUE
WRITE(4,101)
NFI=INT((FI1+FI2)/FIS)
IF(FLOAT(NFI).LT.(FI1+FI2)/FIS) NFI=NFI+1
IX=1
DO 26 J=1,NFI
L=J*INT(FIS)
WRITE(4,100)L
100 FORMAT(79(' ')/T2'I',T4,I3,T8,'I',T12,'WRK',T18,'I',T22,'WKF',
& T29,'I',T35,'R',T40,'I',T45,'WP',T50,'I',T55,'P',T59,'I',
& T63,'TK',T69,'I',T74,'DK',T79,'I')
WRITE(4,104)
101
FORMAT(T2,'I_____I_________I__________I__________I_________I
',
&'________I_________I_________I')
104
FORMAT(T2,'I_ГРАДI___M/C___I___M/C____I____MM____I___M/
C___I',
&'___ATM__I____K____I___MKM___I')
NJ=J
RV1=0.
RV2=0.
RV3=0.
NZ=0
IY=0
DO 23 I=1,NK
NI=I
PRINT *,NJ,NI
PK=WP(I)
DD=DK(I)
RS=R1+I*DRS
IF(R(I).NE.0.) GO TO 4
NEQN=2
T=0.
VII
VIII
363
RELERR=0.6E-6
ABSERR=0.
TFINAL=FIP
TPRINT=FID
IFLAG=1
TOUT=T
K=1
Y(1)=P(I)
Y(2)=WP(I)
CALL SEVAL
(5,P(I),BAR,VPA,BVPA,CVPA,DVPA,VP,VPP,VPPP)
2 CALL
RKF45(EMPTY,NEQN,Y,T,TOUT,RELERR,ABSERR,IFLAG,WORK,
IWORK)
CALL SEVAL
(5,Y(1),BAR,VPA,BVPA,CVPA,DVPA,VP,VPP,VPPP)
GOTO (22,3,17,18,19,20,21,22),IFLAG
3 TOUT=T+TPRINT
IF(T.LT.TFINAL) GOTO 2
WRITE(4,102) I,Y(2),Y(1)*1E-5
102 FORMAT(T2,'I',T5,I2,T8,'I',T18,'I',T29,'I',T40,'I',T42,F7.3,T50
& 'I',T52,F6.4,T59,'I',T69,'I',T79,'I')
DP=(Y(1)-P(I))/FIP
WP(I)=Y(2)
P(I)=Y(1)
IF(I.EQ.NK) WRITE(4,101)
IF(RS.GE.R2-DELPL) GOTO 14
GOTO 23
4 CV2=RV1
CV3=RV2
CV4=RV3
NEQN=4
T=0.
Y(1)=RWK(I)
Y(2)=FWK(I)
Y(3)=R(I)
Y(4)=TK(I)
DD=DK(I)
RELERR=0.6E-6
ABSERR=0.
TFINAL=FIP
TPRINT=FID
IX
IX
364
IFLAG=1
TOUT=T
K=2
CALL SEVAL
(5,P0,BAR,VPA,BVPA,CVPA,DVPA,VP,VPP,VPPP)
CALL SEVAL (5,P0,BAR,TPA,BTPA,CTPA,DTPA,TP,TPP,TPPP)
5 CALL
RKF45(DROP,NEQN,Y,T,TOUT,RELERR,ABSERR,IFLAG,WORK,I
WORK)
TD=Y(4)
CALL SEVAL
(5,P0,BAR,VPA,BVPA,CVPA,DVPA,VP,VPP,VPPP)
CALL SEVAL
(5,TD,TEMP,VKAP,BVKAP,CVKAP,DVKAP,VK,VKP,VKPP)
CALL SEVAL
(5,TD,TEMP,ENTS,BES,CES,DES,EPTK,EPTKP,EPTKPP)
CALL SEVAL
(5,TD,TEMP,ENT,BE,CE,DE,ENTK,ENTKP,ENTKPP)
CALL SEVAL (5,P0,BAR,TPA,BTPA,CTPA,DTPA,TP,TPP,TPPP)
CALL SEVAL (5,TD,TEMP,CP,BCP,CCP,DCP,CPK,CPKP,CPKPP)
CALL SEVAL
(5,P0,DAVL,ENTN,BENS,CENS,DENS,ENTP,ENTPP,ENTPPP)
CALL SEVAL
(5,TD,TEMP,ENTS,BES,CES,DES,EPTK,EPTKP,EPTKPP)
CALL SEVAL
(5,TD,TEMP,ENT,BE,CE,DE,ENTK,ENTKP,ENTKPP)
CALL SEVAL
(5,Y(4),BAR,TPA,BTPA,CTPA,DTPA,TP,TPP,TPPP)
CALL SEVAL (5,TD,TEMP,CP,BCP,CCP,DCP,CPK,CPKP,CPKPP)
CALL SEVAL
(5,Y(4),DAVL,ENTN,BENS,CENS,DENS,ENTP,ENTPP,ENTPPP)
GOTO (22,6,17,18,19,20,21,22),IFLAG
6 TOUT=T+TPRINT
IF(T.LT.TFINAL) GOTO 5
RT=Y(3)
IF(I.EQ.NK)PRINT *,' DROP: R(NK)=',R(NK)
TCKR=Y(1)
TCKF=Y(2)
NEQN=7
Y(1)=RWK(I)
Y(2)=FWK(I)
Y(3)=R(I)
365
Y(4)=P(I)
Y(5)=WP(I)
Y(6)=TK(I)
Y(7)=DK(I)
T=0.
TFINAL=FIP
TPRINT=FID
IFLAG=1
TOUT=T
K=3
CALL SEVAL
(5,P(I),BAR,VPA,BVPA,CVPA,DVPA,VP,VPP,VPPP)
CALL SEVAL (5,P(I),BAR,TPA,BTPA,CTPA,DTPA,TP,TPP,TPPP)
7 CALL
RKF45(VAP,NEQN,Y,T,TOUT,RELERR,ABSERR,IFLAG,WORK,IW
ORK)
TD=Y(6)
CALL SEVAL
(5,Y(4),BAR,VPA,BVPA,CVPA,DVPA,VP,VPP,VPPP)
CALL SEVAL
(5,TD,TEMP,VKAP,BVKAP,CVKAP,DVKAP,VK,VKP,VKPP)
CALL SEVAL
(5,TD,TEMP,ENTS,BES,CES,DES,EPTK,EPTKP,EPTKPP)
CALL SEVAL
(5,TD,TEMP,ENT,BE,CE,DE,ENTK,ENTKP,ENTKPP)
CALL SEVAL
(5,Y(4),BAR,TPA,BTPA,CTPA,DTPA,TP,TPP,TPPP)
CALL SEVAL (5,TD,TEMP,CP,BCP,CCP,DCP,CPK,CPKP,CPKPP)
CALL SEVAL
(5,Y(4),DAVL,ENTN,BENS,CENS,DENS,ENTP,ENTPP,ENTPPP)
GOTO (22,8,17,18,19,20,21,22),IFLAG
8 TOUT=T+TPRINT
IF(T.LT.TFINAL) GOTO 7
RT=Y(3)
TCKR=Y(1)
TCKF=Y(2)
IF(I.EQ.NK)PRINT *,'VAP1: R(NK)=',R(NK)
CV2=RV1
CV3=RV2
CV4=RV3
NEQN=7
Y(1)=RWK(I)
X
366
Y(2)=FWK(I)
Y(3)=R(I)
Y(4)=P(I)
Y(5)=WP(I)
Y(6)=TK(I)
Y(7)=DK(I)
T=0.
TFINAL=FIP
TPRINT=FID
IFLAG=1
TOUT=T
K=4
CALL SEVAL
(5,P(I),BAR,VPA,BVPA,CVPA,DVPA,VP,VPP,VPPP)
CALL SEVAL (5,P(I),BAR,TPA,BTPA,CTPA,DTPA,TP,TPP,TPPP)
9 CALL
RKF45(VAP,NEQN,Y,T,TOUT,RELERR,ABSERR,IFLAG,WORK,IW
ORK)
TD=Y(6)
CALL SEVAL
(5,Y(4),BAR,VPA,BVPA,CVPA,DVPA,VP,VPP,VPPP)
CALL SEVAL
(5,TD,TEMP,VKAP,BVKAP,CVKAP,DVKAP,VK,VKP,VKPP)
CALL SEVAL
(5,TD,TEMP,ENTS,BES,CES,DES,EPTK,EPTKP,EPTKPP)
CALL SEVAL
(5,TD,TEMP,ENT,BE,CE,DE,ENTK,ENTKP,ENTKPP)
CALL SEVAL
(5,Y(4),BAR,TPA,BTPA,CTPA,DTPA,TP,TPP,TPPP)
CALL SEVAL (5,TD,TEMP,CP,BCP,CCP,DCP,CPK,CPKP,CPKPP)
CALL SEVAL
(5,Y(4),DAVL,ENTN,BENS,CENS,DENS,ENTP,ENTPP,ENTPPP)
GOTO (22,10,17,18,19,20,21,22),IFLAG
10 TOUT=T+TPRINT
IF(T.LT.TFINAL) GOTO 9
WRITE(4,103)I,Y(1),Y(2),Y(3)*1E3,Y(5),Y(4)*1E-5,Y(6),Y(7)*1E6
103 FORMAT(T2,'I',T5,I2,T8,'I',T10,F7.4,T18,'I',T20,F8.4,T29,'I',
& T31,F8.4,T40,'I',T42,F7.3,T50,'I',T52,F6.4,T59,'I',T61,F7.3,
& T69,'I',T71,F7.4,T79,'I')
DP=(Y(4)-P(I))/FIP
CM=Y(2)
RWK(I)=Y(1)
XI
367
FWK(I)=Y(2)
R(I)=Y(3)
IF(I.EQ.NK)PRINT *,'VAP2: R(NK)=',R(NK)
P(I)=Y(4)
WP(I)=Y(5)
TK(I)=Y(6)
DK(I)=Y(7)
A3=TK(I)
RV1=CV2
RV2=CV3
RV3=CV4
IF(I.EQ.NK) WRITE(4,101)
IF(RS.LT.R2-DELPL) GOTO 23
DO 11 IS=1,NK
IF(R(IS).GE.R2-DELPL) NZ=NZ+1
11 CONTINUE
PRINT *,'NZ=',NZ
IF(NZ.EQ.0) GOTO 14
NR=NR-NZ
DO 12 IT=1,NR
R(NK-IT+1)=R(NK-NZ-IT+1)
RWK(NK-IT+1)=RWK(NK-NZ-IT+1)
FWK(NK-IT+1)=FWK(NK-NZ-IT+1)
TK(NK-IT+1)=TK(NK-NZ-IT+1)
12 CONTINUE
DO 13 IU=1,NK-NR
R(IU)=0.
13 CONTINUE
14 DD=DK(I)
DG=NZ*SK*PI*DD**3/6./FIP
PRINT*,'DG=',DG,'DELPL=',DELPL
IF(DG.EQ.0..AND.DELPL.EQ.0.) GOTO 23
NEQN=3
T=0.
TFINAL=FIP
TPRINT=FID
IFLAG=1
TOUT=T
K=5
PRINT *,' K=5'
Y(1)=WPL
IF(Y(1).EQ.0.) Y(1)=A1*CM
XI
XII
XIII
368
Y(2)=DELPL
IF(Y(2).EQ.0.) Y(2)=DG/Y(1)*FIP*A2
PK=WP(I)
DD=DK(I)
Y(3)=TPL
IF(Y(3).EQ.0.) Y(3)=A3
P1=P(I)
CALL SEVAL
(5,P1,BAR,VPA,BVPA,CVPA,DVPA,VP,VPP,VPPP)
CALL SEVAL
(5,TD,TEMP,VKAP,BVKAP,CVKAP,DVKAP,VK,VKP,VKPP)
CALL SEVAL
(5,TD,TEMP,ENTS,BES,CES,DES,EPTK,EPTKP,EPTKPP)
CALL SEVAL
(5,TD,TEMP,ENT,BE,CE,DE,ENTK,ENTKP,ENTKPP)
CALL SEVAL (5,P1,BAR,TPA,BTPA,CTPA,DTPA,TP,TPP,TPPP)
CALL SEVAL
(5,TPL,TEMP,VKAP,BVKAP,CVKAP,DVKAP,VPL,VPLP,VPLPP)
CALL SEVAL
(5,TPL,TEMP,CP,BCP,CCP,DCP,CPPL,CPPLP,CPPLPP)
CALL SEVAL (5,TD,TPL,ENT,BE,CE,DE,EPL,EPLP,EPLPP)
CALL SEVAL
(5,TD,TPL,ENTS,BES,CES,DES,EPPL,EPPLP,EPPLPP)
CALL SEVAL
(5,P1,DAVL,ENTN,BENS,CENS,DENS,ENTP,ENTPP,ENTPPP)
15 CALL
RKF45(FILM,NEQN,Y,T,TOUT,RELERR,ABSERR,IFLAG,WORK,IW
ORK)
TPL=Y(3)
CALL SEVAL
(5,P1,BAR,VPA,BVPA,CVPA,DVPA,VP,VPP,VPPP)
CALL SEVAL
(5,TD,TEMP,VKAP,BVKAP,CVKAP,DVKAP,VK,VKP,VKPP)
CALL SEVAL
(5,TD,TEMP,ENTS,BES,CES,DES,EPTK,EPTKP,EPTKPP)
CALL SEVAL
(5,TD,TEMP,ENT,BE,CE,DE,ENTK,ENTKP,ENTKPP)
CALL SEVAL (5,P1,BAR,TPA,BTPA,CTPA,DTPA,TP,TPP,TPPP)
CALL SEVAL
(5,TPL,TEMP,VKAP,BVKAP,CVKAP,DVKAP,VPL,VPLP,VPLPP)
CALL SEVAL
(5,TPL,TEMP,CP,BCP,CCP,DCP,CPPL,CPPLP,CPPLPP)
XIII
369
CALL SEVAL (5,TPL,TEMP,ENT,BE,CE,DE,EPL,EPLP,EPLPP)
CALL SEVAL
(5,TPL,TEMP,ENTS,BES,CES,DES,EPPL,EPPLP,EPPLPP)
CALL SEVAL
(5,P1,DAVL,ENTN,BENS,CENS,DENS,ENTP,ENTPP,ENTPPP)
GOTO (22,16,17,18,19,20,21,22),IFLAG
16 TOUT=T+TPRINT
IF(T.LT.TFINAL) GOTO 15
WPL=Y(1)
DELPL=Y(2)
TPL=Y(3)
GPL=WPL*DELPL/VPL*HB
WRITE(4,*) ' WPL=',Y(1),'DEL=',Y(2),'GPL=',GPL,'TPL=',Y(3)
PRINT*, ' WPL=',Y(1),'DEL=',Y(2),'GPL=',GPL,'TPL=',Y(3)
GOTO 24
17 PRINT *,RELERR,ABSERR
GOTO (2,5,7,9,15),K
18 PRINT *,' МНОГО ШАГОВ'
GOTO (2,5,7,9,15),K
19 ABSERR=1E-6
PRINT *,RELERR,ABSERR
GOTO (2,5,7,9,15),K
20 RELERR=10*RELERR
PRINT *,RELERR,ABSERR
IFLAG=2
GOTO (2,5,7,9,15),K
21 PRINT *, ' МНОГО ДАННЫХ'
IFLAG=2
GOTO (2,5,7,9,15),K
22 PRINT *, ' НЕВЕРНЫЙ ВЫЗОВ'
STOP
23 CONTINUE
PRINT *,'R(NK)=',R(NK),' R2=',R2
24 IF(DELPL.GE.IX*DRS) IY=1
IF(IY.EQ.1) IX=IX+1
NK=NK-IY
IF(FIS*J.NE.FI1) GOTO 26
R5=R1-R3
R1=R3
R2=R4
DO 25 I=1,NK
IF(R(I).NE.0.) R(I)=R(I)-R5
XIII