Télécharger @SATURAT.procedur

Retour à la liste

Numérotation des lignes :

  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.  

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