Esope source files

   1 : C XQOND     SOURCE    CHAT      06/06/01    21:22:17     5450           
   2 :       SUBROUTINE XQOND(CP,ALFAB,ALFAT,N,BETA,
   3 :      1                H,TP,PTOT,XVAP,XH2O,FVAP,IRET)
   4 : 
   5 : C*********************************************************************
   6 : C
   7 : C OBJET   : To compute the condensed mass flux density of H2O at
   8 : C           the wall from a mixture of air and vapour.
   9 : C
  10 : C SYNTAXE : M = QOND CP ALFAB ALFAT H TP PTOT XVAP XH2O 
  11 : C
  12 : C M     : CHPOINT SCAL (condensed mass flux density of H2O (kg/m2/s))
  13 : C alias FVAP
  14 : C CP    : FLOTTANT     (vapour specific heat (J/kg/K))
  15 : C ALFAB : FLOTTANT     (brownian diffusivity (m2/s))
  16 : C ALFAT : FLOTTANT     (thermal diffusivity (m2/s))
  17 : C H     : CHPOINT SCAL (heat exchange coefficient at the wall (W/m2/K))
  18 : C TP    : CHPOINT SCAL (wall temperature (K))
  19 : C PTOT  : CHPOINT SCAL (total pressure of the mixture (Pa))
  20 : C XVAP  : CHPOINT SCAL (vapour quality ((vapour mass)/(water mass)))
  21 : C XH2O  : CHPOINT SCAL                 ((water mass)/(mixture mass))
  22 : C BETA  : FLOTTANT     (coefficient (option, default beta=1))
  23 : C IRET  : ENTIER       (return test ( 0 OK ) (ne 0 error) )
  24 : C
  25 : C
  26 : C M = BETA * H/CP *(ALFAB/ALFAT)**.666 *Ln((PTOT-PSAT(TP))/(PTOT-PVAP))
  27 : C M = 0   IF PVAP < PSAT(TP)
  28 : C
  29 : C IMPORTANT:   All the CHPOINTs MUST! have the same SPG
  30 : C ---------    S.I. units MUST! be used
  31 : C              M will be defined on the same SPG of H,TP,....
  32 : C
  33 : C
  34 : C*********************************************************************
  35 : 
  36 :       IMPLICIT INTEGER(I-N)
  37 :       IMPLICIT REAL*8 (A-H,O-Z)                                                 
  38 : -INC CCREEL
  39 :       DIMENSION H(*),TP(*),PTOT(*),XVAP(*),XH2O(*),FVAP(*)
  40 : *      SAVE RAIR,RVAP
  41 : *      DATA RAIR,RVAP/287.0D0,462.0D0/
  42 : 
  43 :       IRET=0
  44 : 
  45 :       IF(CP.LE.0.D0.OR.ALFAB.LE.0.D0.OR.ALFAT.LE.0.D0.OR.BETA.LE.0.D0)
  46 :      & THEN
  47 :       IRET=6
  48 :       RETURN
  49 :       ENDIF
  50 : 
  51 : 
  52 :       XLEWIS = (ALFAB/ALFAT)**(2.0D0/3.0D0)
  53 :       XLEWIS = XLEWIS*BETA/CP
  54 : *      RRATIO = RAIR/RVAP
  55 :       RRATIO= 287.d0/462.d0
  56 : 
  57 :       DO I=1,N
  58 : 
  59 :         TWALL  = TP(I)
  60 :         PX     = PTOT(I)
  61 :         XSTEAM = XVAP(I)
  62 :         XWATER = XH2O(I)
  63 :         XAIR   = 1.0D0 - XWATER
  64 :         XH     = H(I)
  65 : 
  66 :         IF(XWATER.LT.0.D0.OR.XWATER.GE.(1.D0-XPETIT))THEN
  67 :          IRET=1
  68 :          RETURN
  69 :         ENDIF
  70 :         IF(XSTEAM.LT.0.D0.OR.XSTEAM.GT.(1.D0+1.D-15))THEN
  71 :          IRET=2
  72 :          RETURN
  73 :         ENDIF
  74 :         IF(XH.LE.0.D0)THEN
  75 :          IRET=3
  76 :          RETURN
  77 :         ENDIF
  78 :         IF(TWALL.LE.0.D0)THEN
  79 :          IRET=4
  80 :          RETURN
  81 :         ENDIF
  82 :         IF(PX.LE.0.D0)THEN
  83 :          IRET=5
  84 :          RETURN
  85 :         ENDIF
  86 : 
  87 :         PSAT = 1.055D26 / (TWALL*TWALL*TWALL*TWALL*TWALL)
  88 :         PSAT = 1.013D0  * PSAT*EXP (-7.0D3/TWALL)
  89 : 
  90 :         IF(PX.LE.PSAT)THEN
  91 :          FVAP(I) = 0.D0
  92 :          ELSE
  93 : 
  94 :          PVAP = PX * XSTEAM * XWATER
  95 :          PVAP = PVAP / ((XSTEAM * XWATER) + (XAIR*RRATIO))
  96 :          XFV  = (PX - PSAT) / (PX - PVAP)
  97 :          XFV  = XLEWIS * XH * LOG(XFV)
  98 :          XFV  = MAX(XFV,0.0D0)
  99 :          FVAP(I) = XFV
 100 :         ENDIF
 101 : C       write(6,*)'i,ptot,pvap,psat,j :',i,px,pvap,psat,xfv
 102 : 
 103 :       ENDDO
 104 : 
 105 :       RETURN
 106 :       END
 107 :  
 108 :  
 109 :  
 110 :