Accueil | Cast3M

   1 : $$$$ KR_PRO   NOTICE  CHAT      11/09/12    21:16:48     7124           
   2 :                                              DATE     11/09/12
   3 : 
   4 :   Voir aussi : DARCYSAT
   5 :                                                       HT_PRO    
   6 :     ------------------                  
   7 : 
   8 :     CHPO2 = KR_PRO  TAB1 CHPO1  ;
   9 : 
  10 : 
  11 : 
  12 :     Objet :
  13 :     _______
  14 : 
  15 :     Cette procedure calcule la permeabilite a l'eau
  16 :     en fonction de la saturation reduite S. 
  17 :     Cette procedure est utilisee a partir de la procedure DARCYSAT.
  18 :     La loi est donnee par le sous type de TAB1 :
  19 :     - PUISSANCE
  20 :        k = ks S^B
  21 :     - MUALEM
  22 :        k = ks * s^0.5 * (1 - (1 - (s^(1/M)))^M)^2
  23 :     - BURDINE
  24 :        k = ks * s^2 * (1 - (1 - (s^(1/M)))^M)
  25 :     - MUALEM_BURDINE
  26 :        k = ks * s^A * (1 - (1 - (s^(1/M)))^M)^B
  27 :     - BROOKS_COREY
  28 :        k = ks * s^A * S^(2/lambda + B)
  29 :     - EXPONENTIELLE
  30 :        k = ks * C / (C - 1 + exp (- alpha Pw N))
  31 :     - LOGARITHMIQUE
  32 :        k = ks * C / (C + (log (- alpha Pw))^N)
  33 :        ou -alpha Pw tronque a 1 quand devient inferieur a 1.
  34 : 
  35 :     
  36 :     Commentaires
  37 :     -------------
  38 : 
  39 :     TAB1  : table contenant les caracteristiques physiques, ayant
  40 :             pour indices 
  41 : 
  42 :             SI loi 'PUISSANCE'
  43 :             'ALPHA'   : coef. B (s.d.)
  44 : 
  45 :             SI loi 'MUALEM' 'BURDINE' et 'MUALEM_BURDINE'
  46 :             'COEF_M' : coef M
  47 : 
  48 :             SI loi 'MUALEM_BURDINE' et 'BROOKS_COREY'
  49 :             'COEF_A' : coef A
  50 :             'COEF_B' : coef B        
  51 :     
  52 :             SI 'BROOKS_COREY'
  53 :             'LAMBDA' : coef lambda
  54 : 
  55 :             SI 'EXPONENTIELLE' ou 'LOGARITHMIQUE'
  56 :             'ALPHA'  : alpha
  57 :             'COEF_C' : C
  58 :             'COEF_N' : N
  59 : 
  60 : 
  61 :             'PERMSAT' : coef. Ks, permeabilite a saturation (m/s)
  62 :                       ('FLOTTANT', ou 'CHPOINT'
  63 :                       dont les composantes sont celles demandees par la
  64 :                       formulation, voir 'MODE' 'DARCY')
  65 :                       sous-formulations autorisees : ISOTROPE, 
  66 :                       ORTHOTROPE, ANISOTROPE 2D ou 3D.   
  67 :             'MODELE'  : objet modele correspondant au domaine concerne
  68 : 
  69 :     CHPO1 : saturation reduite ('CHPO' centre ou 'FACE') pour les 
  70 :             lois PUISSANCE, MUALEM, BURDINE, MUALEM_BURDINE, BROOKS_COREY
  71 :             pression en eau pour les lois EXPONENTIELLE et LOGARITHMIQUE
  72 : 
  73 : 
  74 :     CHPO2 : permeabilite totale en eau (m/s)
  75 : 
  76 :     
  77 :     Remarques :
  78 :     -----------
  79 : 
  80 :     1. Le champ en sortie ont les meme points d'appui que le champ en
  81 :        entree.
  82 : 
  83 :     2. Il a autant de composantes que l'indice PERMSAT.
  84 : 
  85 :     3. Les directions d'anisotropie sont (0 0 1), (0 1 0) et (0 0 1)
  86 :        dans l'ordre. Charge a l'utilisateur de definir le tenseur
  87 :        PERMSAT dans ce repere-la. 
  88 :         
  89 :