Télécharger qond.eso

Retour à la liste

Numérotation des lignes :

  1. C QOND SOURCE CHAT 05/01/13 02:40:30 5004
  2. SUBROUTINE QOND
  3. C************************************************************************
  4. C
  5. C OBJET : To compute the condensed mass flux density of H2O at
  6. C the wall from a mixture of air and vapour
  7. C
  8. C SYNTAXE : M = QOND CP ALFAB ALFAT H TP PTOT XVAP XH2O <BETA>
  9. C
  10. C M : CHPOINT SCAL (condensed mass flux density of H2O (kg/m2/s))
  11. C CP : FLOTTANT (vapour specific heat (J/kg/K))
  12. C ALFAB : FLOTTANT (brownian diffusivity (m2/s))
  13. C ALFAT : FLOTTANT (thermal diffusivity (m2/s))
  14. C H : CHPOINT SCAL (heat exchange coefficient at the wall (W/m2/K))
  15. C TP : CHPOINT SCAL (wall temperature (K))
  16. C PTOT : CHPOINT SCAL (total pressure of the mixture (Pa))
  17. C XVAP : CHPOINT SCAL (vapour quality ((vapour mass)/(water mass)))
  18. C XH2O : CHPOINT SCAL ((water mass)/(mixture mass))
  19. C BETA : FLOTTANT (coefficient (option, default beta=1))
  20. C
  21. C
  22. C M = BETA * H/CP *(ALFAB/ALFAT)**.666 *Ln((PTOT-PSAT(TP))/(PTOT-PVAP))
  23. C M = 0 IF PVAP < PSAT(TP)
  24. C
  25. C IMPORTANT: All the CHPOINTs MUST! have the same SPG
  26. C --------- S.I. units MUST! be used
  27. C M will be defined on the same SPG of H,TP,....
  28. C
  29. C************************************************************************
  30. IMPLICIT INTEGER(I-N)
  31. IMPLICIT REAL*8 (A-H,O-Z)
  32. CHARACTER*8 TYPE
  33. C***
  34.  
  35. -INC PPARAM
  36. -INC CCOPTIO
  37. -INC SMCHPOI
  38. POINTEUR MH.MCHPOI, MTP.MCHPOI,MPTOT.MCHPOI,
  39. 1 MXVAP.MCHPOI,MXH2O.MCHPOI,MFVAP.MCHPOI
  40. POINTEUR LH.MSOUPO, LTP.MSOUPO,LPTOT.MSOUPO,
  41. 1 LXVAP.MSOUPO,LXH2O.MSOUPO,LFVAP.MSOUPO
  42. POINTEUR NH.MPOVAL, NTP.MPOVAL,NPTOT.MPOVAL,
  43. 1 NXVAP.MPOVAL,NXH2O.MPOVAL,NFVAP.MPOVAL
  44.  
  45. C Reading operands
  46.  
  47. CALL LIRREE(CP,1,IRET)
  48. IF(IRET.EQ.0) RETURN
  49. CALL LIRREE(ALFAB,1,IRET)
  50. IF(IRET.EQ.0) RETURN
  51. CALL LIRREE(ALFAT,1,IRET)
  52. IF(IRET.EQ.0) RETURN
  53. CALL LIROBJ('CHPOINT ',MH,1,IRET)
  54. IF(IRET.EQ.0) RETURN
  55. CALL LIROBJ('CHPOINT ',MTP,1,IRET)
  56. IF(IRET.EQ.0) RETURN
  57. CALL LIROBJ('CHPOINT ',MPTOT,1,IRET)
  58. IF(IRET.EQ.0) RETURN
  59. CALL LIROBJ('CHPOINT ',MXVAP,1,IRET)
  60. IF(IRET.EQ.0) RETURN
  61. CALL LIROBJ('CHPOINT ',MXH2O,1,IRET)
  62. IF(IRET.EQ.0) RETURN
  63. SEGACT MXH2O
  64. BETA=1.0
  65. CALL LIRREE(BETA,0,IRET)
  66.  
  67. C Verfying operands
  68.  
  69. CALL LICHT(MH,NH,TYPE,IGEO)
  70. CALL LICHT(MTP,NTP,TYPE,IGEO0)
  71. IF(IGEO.NE.IGEO0) THEN
  72. C Indice %m1:8 : L'objet %m9:16 n'a pas le bon support géométrique
  73. MOTERR(1: 8) = ' TP '
  74. CALL ERREUR(788)
  75. RETURN
  76. ENDIF
  77. CALL LICHT(MPTOT,NPTOT,TYPE,IGEO0)
  78. IF(IGEO.NE.IGEO0) THEN
  79. C Indice %m1:8 : L'objet %m9:16 n'a pas le bon support géométrique
  80. MOTERR(1: 8) = ' PT '
  81. CALL ERREUR(788)
  82. RETURN
  83. ENDIF
  84. CALL LICHT(MXVAP,NXVAP,TYPE,IGEO0)
  85. IF(IGEO.NE.IGEO0) THEN
  86. C Indice %m1:8 : L'objet %m9:16 n'a pas le bon support géométrique
  87. MOTERR(1: 8) = ' YVAP'
  88. CALL ERREUR(788)
  89. RETURN
  90. ENDIF
  91. CALL LICHT(MXH2O,NXH2O,TYPE,IGEO0)
  92. IF(IGEO.NE.IGEO0) THEN
  93. C Indice %m1:8 : L'objet %m9:16 n'a pas le bon support géométrique
  94. MOTERR(1: 8) = ' YH2O '
  95. CALL ERREUR(788)
  96. RETURN
  97. ENDIF
  98.  
  99. C Creating the chpoint result
  100.  
  101. NC=1
  102. CALL CRCHPT(TYPE,IGEO,NC,MFVAP)
  103. CALL LICHT(MFVAP,NFVAP,TYPE,IGEO)
  104. N=NFVAP.VPOCHA(/1)
  105.  
  106. C Computing the condensed mass flux
  107.  
  108. CALL XQOND(CP,ALFAB,ALFAT,N,BETA,
  109. 1 NH.VPOCHA, NTP.VPOCHA, NPTOT.VPOCHA,
  110. 2 NXVAP.VPOCHA,NXH2O.VPOCHA,NFVAP.VPOCHA,
  111. 3 IRET)
  112.  
  113.  
  114. IF(IRET.NE.0) THEN
  115. C Tache impossible. Probablement données erronées
  116. CALL ERREUR(26)
  117. RETURN
  118. ENDIF
  119.  
  120. CALL ECROBJ('CHPOINT ',MFVAP)
  121.  
  122. SEGDES NFVAP,MH,NH,NTP,NPTOT,NXVAP,NXH2O
  123.  
  124. RETURN
  125. 100 FORMAT(2X,'Les CHPOINTs n''ont pas le meme SPG')
  126. END
  127.  
  128.  
  129.  
  130.  
  131.  
  132.  

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