Manual pages | Cast3M

$$$$ HT_PRO   NOTICE  FANDEUR   13/03/15    21:15:00     7712           
                                             DATE     13/03/15

                                         Voir aussi : DARCYSAT
                                                      KR_PRO   
    ----------------

    CHPO2 CHPO3 CHPO4 = HT_PRO  TAB1 (CHPO1) ;


    Description :
    _____________

    This procedure is used with DARCYSAT.

    It calculates the saturation, the water yield and the capacity as 
    functions of water pressure P, using the law given by TAB1 subtype :
    - VAN_GENUCHTEN
       Reduced Saturation in [0,1] : S = [1 + (-beta.Pw)**n]**(-m)
    - EXPONENTIELLE
       S = C / ([exp(-beta.Pw)]**N + C - 1)
    - LOGARITHMIQUE
       (-beta.Pw) set to 1 when lower than 1.
       S = C / ([ln(-beta.Pw)]**N + C)

    with S = (TH2O - teta_r) / (poros - teta_r)


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

    TAB1  : table containing the physical properties, which indexes are
            BHETA, NEXP, MEXP : law coefficients
            BHETA, COEF_N, COEF_C : coef for log or exp. laws
            PORO              : porosity (s. d.), poros
                                (FLOTTANT or CHAMP-POINT centers,
                                 comp 'SCAL')
            TERESIDU          : residual water yield (n. d.), teta_r
                                (FLOTTANT or CHAMP-POINT centers,
                                 comp 'SCAL')

    CHPO1 : Water pressure (negative if non saturated)
            = succion if Pgaz = 0 (1 atm for convention)
            (CHAMP-POINT centers or faces according to MOT1,
             component 'SCAL')

    CHPO2 : Reduced saturation , S
            (CHAMP-POINT centers or faces according to MOT1, comp 'SCAL',
             component 'SCAL', s.d. between 0 and 1)

    CHPO3 : water yield, TH2O
            (CHAMP-POINT centers, component 'SCAL', comp 'SCAL',
             s.d. between teta_r and poros)

    CHPO4 : capacity, analytically calculated, CAPP
            (CHAMP-POINT centers, component 'SCAL', comp 'SCAL',

    Remarks :
    ---------

    If the CHPO1 objet is not given, only testing procedure's arguments 
    withoud making calculations.