Manual pages | 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  ;



    Description :
    _____________

    This procedure is used with DARCYSAT.
    It calculates water permeability, as function of
    reduced saturation S, using a law given by TAB1 subtype :
    - 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)
       where -alpha Pw is set to 1 when lower than 1.


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

    TAB1  : table containing the physical properties, which indexes are 

            if law 'PUISSANCE'
            'ALPHA'   : coef. B (s.d.)

            if law 'MUALEM' 'BURDINE' et 'MUALEM_BURDINE'
            'COEF_M' : coef M

            if law 'MUALEM_BURDINE' et 'BROOKS_COREY'
            'COEF_A' : coef A
            'COEF_B' : coef B        
    
            if law 'BROOKS_COREY'
            'LAMBDA' : coef lambda

            if law 'EXPONENTIELLE' ou 'LOGARITHMIQUE'
            'ALPHA'  : alpha
            'COEF_C' : C
            'COEF_N' : N


            ALPHA   : coef. B (s.d.)
            PERMSAT : coef. Ks, saturated permeability (m/s)
                      ('FLOTTANT', or 'CHPOINT'
                      which components are those requiered by the
                      formulation, see 'MODE' 'DARCY')
                      authorized under-formulations : ISOTROPE,
                      ORTHOTROPE, ANISOTROPE.
            MODELE  : model objet corresponding to the current geometry

    CHPO1 : Reduced saturation ('CHPO' centre or 'FACE') for the laws
            PUISSANCE, MUALEM, BURDINE, MUALEM_BURDINE, BROOKS_COREY
            water pressure for  'EXPONENTIELLE' ou 'LOGARITHMIQUE' laws

    CHPO2 : total water permeability

    Remarks :
    ---------

    1. Output fields meshes are those of the input one.
        
    2. the output field has as many components as the PERMSAT index

    3. Anisotropy directions are (0 0 1), (0 1 0) and (0 0 1)
       respectively. It's up to the user to define PERMSAT tensor in
       that basis.