Procedures | Cast3M

   1 : * @SATURAT  PROCEDUR  PASCAL    08/10/06    21:15:05     6176           
   2 : * PROCEDURE @SATURAT
   3 : *---------------------------------------------------------------------*
   4 : * NOM         : @SATURAT                                              *
   5 : *                                                                     *
   6 : * DESCRIPTION :                                                       *
   7 : *                                                                     *
   8 : * SYNTAXE     :                                                       *
   9 : *                                                                     *
  10 : *                                                                     *
  11 : * LANGAGE     : GIBIANE-CAST3M                                        *
  12 : * AUTEUR      : S. PASCAL (CEA/DEN/DM2S/SEMT/LM2S)                    *
  13 : * COURRIEL    : serge.pascal@cea.fr                                   *
  14 : *---------------------------------------------------------------------*
  15 : * VERSION    : v1, xx/xx/xx08, version initiale                       *
  16 : * HISTORIQUE : v1, xx/xx/xx08, creation                               *
  17 : * HISTORIQUE :                                                        *
  18 : * HISTORIQUE :                                                        *
  19 : *---------------------------------------------------------------------*
  20 : * Priere de PRENDRE LE TEMPS de completer les commentaires            *
  21 : * en cas de modification de ce sous-programme afin de faciliter       *
  22 : * la maintenance !                                                    *
  23 : *---------------------------------------------------------------------*
  24 : 'DEBP' @SATURAT ;
  25 : 'ARGU' CHTK1*'CHPOINT' CHPC1*'CHPOINT' ;
  26 : * Ok !
  27 : * ENTREES                                                             *
  28 : *  - Temperature (K)                                                  *
  29 : *  - Pression capillaire (Pa)                                         *
  30 : * SORTIE                                                              *
  31 : *  - Saturation (-) : SW0                                             *
  32 : * PARAMETRES                                                          *
  33 : * |-------------|---------|---------|---------|-----------|
  34 : * |coefficient  | BO      |   BHP   |   MO    |     MH    |
  35 : * ||------------|---------|---------|---------|-----------|
  36 : * |a[MPa]       |18.6237  |46.9364  |37.5479  |96.2837    |
  37 : * |-------------|---------|---------|---------|-----------|
  38 : * |b[-]         |2.2748   |2.0601   |2.1684   |1.9540     |
  39 : * |-------------|---------|---------|---------|-----------|
  40 : * Q2       : parametre de calcul du parametre "a", relatif a T>Tcrit  *
  41 : * Q3       : parametre de calcul du parametre "a", relatif a T=20Â°C   *
  42 : * PB1      : exposant "b"                                             *
  43 : * GN1      : exposant "N" de E(theta)                                 *
  44 : * PZ1      : parametre "z" de E(theta) qd T > TCRI1                   *
  45 : * E0       : valeur inconnue => 1. !                                  *
  46 : * TREF1    : Temperature (K) "ambiante" (initiale ?)                  *
  47 : * TSAT1    : Temperature (K) au-dessus de laquelle le param. "a" vari *
  48 : * TCRI1    : Temperature (K) critique                                 *
  49 : Q2       = 7.E6 ;
  50 : Q3       = 18.62E6 ;
  51 : PB1      = 2.27 ;
  52 : GN1      = 1.2 ;
  53 : PZ1      = 0.5 ;
  54 : E0       = 1. ;
  55 : TREF1    = 293.15 ;
  56 : TSAT1    = 373.15 ; 
  57 : TCRI1    = 647.15 ;
  58 : *                                                                     *
  59 : * Calcul du parametre "a" : PA0                                       *
  60 : Q0       = (CHTK1 - TSAT1) / (TCRI1 - TSAT1) ;
  61 : Q02      = -3. * (Q0 ** 2) ;
  62 : Q03      =  2. * (Q0 ** 3) ;
  63 : Q0       = (Q3 - Q2) * (1. + Q02 + Q03) ;
  64 : MASQI1   = CHTK1 'MASQ' 'EGIN' TSAT1 ;
  65 : MASQS1   = CHTK1 'MASQ' 'SUPE' TSAT1 ;
  66 : PA0      = (Q3 * MASQI1) + ((Q0 + Q2) * MASQS1) ;
  67 : *                                                                     *
  68 : * Calcul de E(theta) : ETHETA0                                        *
  69 : MASQI2   = CHTK1 'MASQ' 'INFE' TCRI1 ;
  70 : MASQS2   = CHTK1 'MASQ' 'EGSU' TCRI1 ;
  71 : CHTKSC1  = MASQI2 * CHTK1 ;
  72 : ETHETA1  = (((TCRI1 - CHTKSC1) ** -1) * (TCRI1 - TREF1)) ** GN1 ;
  73 : ETHETA2  = GN1 / PZ1 * (CHTK1 - TCRI1 + PZ1) ;
  74 : ETHETA2  = E0 * (1. + ETHETA2) ;
  75 : * ATTENTION : jeux de donnees SdP ne distinguent pas T>TSAT1...       *
  76 : *ETHETA0  = MASQI1 + (ETHETA1 * MASQI2 * MASQS1) + (ETHETA2 * MASQS2) ;
  77 : ETHETA0  = ETHETA1 ;
  78 : *                                                                     *
  79 : * La saturation :                                                     *
  80 : M1SPB1   = -1. / PB1 ;
  81 : PA0      = PA0 ** -1 ;
  82 : SW0      = (ETHETA0 * PA0 * CHPC1) ** (1. / (1. + M1SPB1)) + 1. ;
  83 : SW0      = SW0 ** M1SPB1 ;
  84 : *                                                                     *
  85 : 'RESP' SW0 ;
  86 : *                                                                     *
  87 : 'FINP' ;
  88 : *---------------------------------------------------------------------*
  89 : * FIN PROCEDURE @SATURAT
  90 :