Télécharger la plaquette Cast3M


fr en

Télécharger @satura.procedur

Retour à la liste

Numérotation des lignes : 

  1. * @SATURA   PROCEDUR  FD218221  26/02/16    21:15:08     12474          
  2. *---------------------------------------------------------------------*
  3. * PROCEDURE @SATURA                                                   *
  4. *---------------------------------------------------------------------*
  5. * NOM         : @SATURA                                               *
  6. *                                                                     *
  7. * LANGAGE     : GIBIANE-CAST3M                                        *
  8. * AUTEUR      : G. Sciumè (I2M - University of Bordeaux)              *
  9. * COURRIEL    : giuseppe.sciume@u-bordeuax.fr                         *
  10. *---------------------------------------------------------------------*
  11. * COMMENTAIRE: Calcul le degré de saturation en liquide du milieu     *
  12. *              poreux                                                 *
  13. *---------------------------------------------------------------------*
  14. * VERSION    : v1, 01/11/2014, version initiale                       *
  15. * HISTORIQUE :                                                        *
  16. *---------------------------------------------------------------------*
  17. * APPELLE PAR    : @BETONJA  @PARAMEC                                 *
  18. * DONNES D'ENTRE :                                                    *
  19. * - CHPC1 = pression capillaire                                       *
  20. * - HYDR  = degre d'hydratation                                       *
  21. * - VGA   = parametre de calcul du parametre "a", relatif a T=20°C    *
  22. * - VGB   = exposant "b"                                              *
  23. * - VGC   = nouveau parametre c                                       *
  24. * - VGHY  = nouveau parametre h0                                      *
  25. *---------------------------------------------------------------------*
  26. DEBP @SATURA CHPC1*'CHPOINT' CHTK1*'CHPOINT' HYDR*'CHPOINT' 
  27.        VGA*'CHPOINT' VGB*'CHPOINT' VGC*'CHPOINT' VGHY*'CHPOINT';
  28. *
  29. TKC = 647.15;
  30.  
  31. * Impac of evolution of microstructure
  32. HYD_EFF = EXP (VGC*(LOG((HYDR + VGHY)*((1. + VGHY)**(-1.)))));
  33.  
  34. * Impact of temperature on water surface tension      
  35. TEMP1 = 'BORNER' (1. - (CHTK1*(1./TKC))) 'SCAL' 'MINIMUM' 0.;
  36. GAW   = 0.2358 * (TEMP1**1.256) * (1. - (0.625 * TEMP1));
  37. GAW0  = 0.07273;
  38. TEM_EFF = (GAW + (0.05 * GAW0))*(1. / (1.05 * GAW0));  
  39.  
  40. * Coupling the two effects
  41. VGA_mod = VGA * HYD_EFF * TEM_EFF;
  42. *LIST 'RESUME' VGA;
  43. *LIST 'RESUME' HYD_EFF;
  44. *LIST 'RESUME' TEM_EFF;
  45.  
  46. * Van Gennuchten law
  47. temp1 = ABS ((VGA_mod**(-1.))*CHPC1);
  48. temp1 = 'BORNER' TEMP1 'SCAL' 'MINIMUM' 1.E-20;
  49. temp2 = VGB*((VGB - 1.)**(-1.));
  50. temp3 = EXP (temp2*(LOG(temp1)));
  51. temp4 = ABS (temp3 + 1.);
  52. temp5 = (-1.)*(VGB**(-1.));
  53. *SW0 = temp4 ** temp5;
  54. SW0 = EXP (temp5*(LOG(temp4)));
  55. SW0 = 'BORNER' SW0 'SCAL' 'MAXIMUM' 0.999;
  56.  
  57. 'FINP' SW0 ;
  58. ***********************************************************************
  59. * FIN PROCEDURE @SATURA
  60. ***********************************************************************
  61.  
  62.  

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