Notices | Cast3M

$$$$ KR_PRO   NOTICE  CHAT      11/09/12    21:16:48     7124           
                                             DATE     11/09/12

                                         Voir aussi : DARCYSAT
                                                      HT_PRO    
    ------------------                  

    CHPO2 = KR_PRO  TAB1 CHPO1  ;



    Objet :
    _______

    Cette procedure calcule la permeabilite a l'eau
    en fonction de la saturation reduite S. 
    Cette procedure est utilisee a partir de la procedure DARCYSAT.
    La loi est donnee par le sous type de TAB1 :
    - PUISSANCE
       k = ks S^B
    - MUALEM
       k = ks * s^0.5 * (1 - (1 - (s^(1/M)))^M)^2
    - BURDINE
       k = ks * s^2 * (1 - (1 - (s^(1/M)))^M)
    - MUALEM_BURDINE
       k = ks * s^A * (1 - (1 - (s^(1/M)))^M)^B
    - BROOKS_COREY
       k = ks * s^A * S^(2/lambda + B)
    - EXPONENTIELLE
       k = ks * C / (C - 1 + exp (- alpha Pw N))
    - LOGARITHMIQUE
       k = ks * C / (C + (log (- alpha Pw))^N)
       ou -alpha Pw tronque a 1 quand devient inferieur a 1.

    
    Commentaires
    -------------

    TAB1  : table contenant les caracteristiques physiques, ayant
            pour indices 

            SI loi 'PUISSANCE'
            'ALPHA'   : coef. B (s.d.)

            SI loi 'MUALEM' 'BURDINE' et 'MUALEM_BURDINE'
            'COEF_M' : coef M

            SI loi 'MUALEM_BURDINE' et 'BROOKS_COREY'
            'COEF_A' : coef A
            'COEF_B' : coef B        
    
            SI 'BROOKS_COREY'
            'LAMBDA' : coef lambda

            SI 'EXPONENTIELLE' ou 'LOGARITHMIQUE'
            'ALPHA'  : alpha
            'COEF_C' : C
            'COEF_N' : N


            'PERMSAT' : coef. Ks, permeabilite a saturation (m/s)
                      ('FLOTTANT', ou 'CHPOINT'
                      dont les composantes sont celles demandees par la
                      formulation, voir 'MODE' 'DARCY')
                      sous-formulations autorisees : ISOTROPE, 
                      ORTHOTROPE, ANISOTROPE 2D ou 3D.   
            'MODELE'  : objet modele correspondant au domaine concerne

    CHPO1 : saturation reduite ('CHPO' centre ou 'FACE') pour les 
            lois PUISSANCE, MUALEM, BURDINE, MUALEM_BURDINE, BROOKS_COREY
            pression en eau pour les lois EXPONENTIELLE et LOGARITHMIQUE


    CHPO2 : permeabilite totale en eau (m/s)

    
    Remarques :
    -----------

    1. Le champ en sortie ont les meme points d'appui que le champ en
       entree.

    2. Il a autant de composantes que l'indice PERMSAT.

    3. Les directions d'anisotropie sont (0 0 1), (0 1 0) et (0 0 1)
       dans l'ordre. Charge a l'utilisateur de definir le tenseur
       PERMSAT dans ce repere-la.