Télécharger bsur.eso

Retour à la liste

Numérotation des lignes :

  1. C BSUR SOURCE CB215821 16/04/21 21:15:24 8920
  2. SUBROUTINE BSUR (X,DX,XL,RUG,XW,XN,TN,EN,BN,KIMP,PSLIM,REL,RINDEX,
  3. & P1,T1,QAE,QEE1,PHI1,P2,T2,U2,QEE2,PHI2,QW2,RE,H,PSQ,
  4. & NPP,ITP,PF,PP,DPF,DPP,RECU,XKUL,XKUT1,XKUT2,XKUT3,XKUT4)
  5.  
  6. IMPLICIT INTEGER(I-N)
  7. IMPLICIT REAL*8 (A-H,O-Z)
  8. DIMENSION XN(NPP),TN(NPP),EN(NPP),BN(NPP)
  9.  
  10. C operateur FUITE
  11. C>>> superheated vapour
  12. C version with wall condensation
  13. C
  14. C QAE: QA.E constant air flowrate
  15. C QEEi: QEi.E water flowrate
  16. C indice 1,2 : entree et sortie de troncon
  17. C RECU,XKUL,XKUT1,XKUT2,XKUT3,XKUT4 : coef lois de frot utilisateur
  18. C
  19. C H correspond au flux d'energie total cede a la paroi
  20. C HM (convection) est utilise pour le calcul de T2
  21. C
  22. HW=0.D0
  23. RINDEX = 1.0
  24. DQEE = 0.D0
  25.  
  26. CALL BPHYS(T0,P0,RA,RV,CA,CV,CL,XLAT0,ROL,XKL,XKT,XREL)
  27.  
  28. C interpolation de la temperature de paroi
  29. CALL BTPAR(XN,TN,X,TP,NPP,ITP)
  30. TP=TP+T0
  31. C epaisseur a l'entree
  32. CALL BTPAR(XN,EN,X,E1,NPP,ITP)
  33. C etendue a l'entree
  34. CALL BTPAR(XN,BN,X,B1,NPP,ITP)
  35.  
  36. QA1=QAE/E1/B1
  37. QE1=QEE1/E1/B1
  38. Q = QA1 + QE1
  39. C calcul de R et CP
  40. CALL BCAR(RA,RV,PHI1,CV,CA,Q,RQ,R,QA0,QV,CP)
  41. IF(KIMP.GE.1) THEN
  42. write(6,*) 'entree bsur : E1,B1,PHI1= ',E1,B1,PHI1
  43. ENDIF
  44.  
  45. RO=P1/R/T1
  46. U1=Q/RO
  47.  
  48. U=U1
  49. T=T1
  50. P=P1
  51. QE=QE1
  52. QA=QA1
  53. QX=QA+QE
  54. PHI=PHI1
  55.  
  56. IF ((X+DX).GT.1.) DX=1-X
  57. X=X+DX
  58.  
  59. DU=0
  60. T2=0
  61. P2=0
  62. c ouverture en sortie de maille
  63. CALL BTPAR(XN,EN,X,E2,NPP,ITP)
  64. c etendue en sortie de maille
  65. CALL BTPAR(XN,BN,X,B2,NPP,ITP)
  66. QA2=QAE/E2/B2
  67. E=E1
  68. B=B1
  69.  
  70. NITER=0
  71.  
  72. 10 CONTINUE
  73.  
  74. NITER=NITER+1
  75.  
  76. c write(6,*) ' bsur X,NITER,T1,P1,RE',X,NITER,T1,P1,RE
  77. TT2=T2
  78. PP2=P2
  79.  
  80. XMU=BMUG(T,PHI,T0)
  81. RE=QX*2*E/XMU
  82. PR=XMU*CP/BLA(T,T0)
  83. DPP = XL*QX*QX/4/E
  84. BK = BKFRO(RE,REL,XKL,XKT,2*E,RUG,RECU,XKUL,XKUT1,XKUT2,XKUT3,
  85. *XKUT4)
  86. DP=DX*BK*DPP/RO
  87.  
  88. P2=P1-DP
  89. IF (P2.LT.PSLIM) THEN
  90. PSQ = -1.
  91. RETURN
  92. ENDIF
  93.  
  94. S=0.D0
  95. C*
  96. XLA=BLA(T,T0)
  97.  
  98. HM = BHECH(RE,REL,PR,XLA,E)
  99.  
  100. C if vapour at the inlet
  101.  
  102. IF(PHI1.GT.1.D-2) THEN
  103. C
  104. C wall condensation
  105. D=2.5D-5
  106. A=2.3D-5
  107. PP=BPSAT(TP)
  108. PV=PHI1*P1
  109. RAPP=(P1-PP)/(P1-PV)
  110.  
  111. IF (RAPP.GT.(1.0001)) THEN
  112. XJ=XW*(HM/CV)*(D/A)**(2./3.)*LOG(RAPP)
  113. ELSE
  114. XJ=0.D0
  115. ENDIF
  116.  
  117. DQEE=-2.*XJ*DX*XL*B
  118. XLAT=XLAT0+(CV-CL)*T
  119. SW=(XLAT/QX/CP/E)*DQEE
  120. C**
  121. DELTAT = T-TP
  122. IF (ABS(DELTAT).GT.1D-3) THEN
  123. HW=XJ*XLAT/(T-TP)
  124. C write(6,*) ' T TP HW ',T,TP,HW
  125. ENDIF
  126.  
  127. QEE2=QEE1+DQEE
  128. QE2=QEE2/E2/B2
  129.  
  130. IF ((QEE2/QEE1).LT.(-1.D-5)) THEN
  131. write(6,*) 'bsur : relative water flowrate ',QE2/QE1
  132. RINDEX = 0.5
  133. RETURN
  134. ELSE
  135. A=RA*QA2/RV/QE2
  136. PHI2=1./(1.+A)
  137. ENDIF
  138. * no vapour at the inlet
  139. ELSE
  140. XLAT=XLAT0+(CV-CL)*T
  141. HW=0.D0
  142. PHI2=0.D0
  143. QEE2=0.D0
  144. QE2=0.D0
  145. ENDIF
  146. IF(KIMP.GE.1) THEN
  147. write(6,*) 'bsur: NITER,PHI2,QA2,QE2= ',NITER,PHI2,QA2,QE2
  148. ENDIF
  149. QEEX=(QEE1+QEE2)/2.
  150.  
  151. T2 = BT2(T1,TP,S,HM,PR,RE,XLA,XL,DX)
  152. IF(KIMP.GE.1) WRITE(*,*) 'bsur : T2= ',T2
  153.  
  154. U2 = (QA2+QE2)*R*T2/P2
  155. RO=(P1+P2)/R/(T1+T2)
  156. T=(T1+T2)/2
  157. P=(P1+P2)/2
  158. PHI=(PHI1+PHI2)/2
  159.  
  160. E=(E1+E2)/2
  161. B=(B1+B2)/2
  162. QA=QAE/E/B
  163. QE=QEEX/E/B
  164. QX=QA+QE
  165.  
  166. QW2=(QEE2-QEE1)/DX/XL/B
  167.  
  168. ERT=ABS((T2-TT2)/T2)
  169. ERP=ABS((P2-PP2)/P2)
  170.  
  171. IF(((ERT.GT.1E-4).OR.(ERP.GT.1E-4))
  172. $ .AND.(NITER.LE.10)) GOTO 10
  173.  
  174. IF(KIMP.GE.2) THEN
  175. write(6,*) ' HM HW HW/HM ',HM,HW,HW/HM
  176. ENDIF
  177.  
  178. C calcul des puissances echangees
  179.  
  180. RAPH=(CV*T+XLAT0)/XLAT
  181. H=HM+HW*RAPH
  182.  
  183. C WRITE(6,*) ' bsur ',H
  184.  
  185. HH1=(QAE*CA*T1)+(QEE1*(CV*T1+XLAT0))
  186. HH2=(QAE*CA*T2)+(QEE2*(CV*T2+XLAT0))
  187. C ! PF = HH2-HH1
  188. PF = (HH2-HH1) - DQEE*(CV*T1+XLAT0)
  189. PP = H*(T-TP)*2*DX*XL*B
  190. IF(KIMP.EQ.-3) THEN
  191. write(6,*) ' DQEE CV T1 XLAT0 ',DQEE,CV,T1,XLAT0
  192. write(6,1000) X,RAPH,PF,PP,HH1,HH2
  193. 1000 FORMAT(1X,' bsur X RH DH EC HH1 HH2 ',6E12.5)
  194. ENDIF
  195. IF(KIMP.NE.0.AND.X.EQ.1.) THEN
  196. write(6,2110) X,RE,BK
  197. 2110 format(1X,'bsur X RE BKRO ',3E12.5)
  198. ENDIF
  199. DPF=PF/DX/XL/B
  200. DPP=PP/DX/XL/B
  201.  
  202. RETURN
  203. END
  204.  
  205.  
  206.  
  207.  
  208.  
  209.  
  210.  
  211.  
  212.  
  213.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales