Manual pages | Cast3M

$$$$ K_PRO    NOTICE  CHAT      11/09/12    21:16:47     7124           
                                             DATE     11/09/12

     Procedure K_PRO                         Voir aussi : DARCYSAT 
     ---------------                                      HT_PRO   
        K1 = K_PRO TAB1 TAB2 ;                   


    Description :
    _____________

    This procedure enables the user to calculated effective permeability
    in a non saturated porous media. 


     K1  (type CHPOINT) permeability field

     TAB1 : (type TABLE) table structure containing at the following indexes

         .'PNS_PROV' (type CHPOINT) : pressure head field (negative or nul)  
         .'TH2O_PROV' (type CHPOINT): water content field 
         .'SATURATION_PROV' (type CHPOINT) : saturation field


     TAB2 : (type TABLE) table structure containing permeability parameters.
            standard permeability law is : K = K_sat * (S **  alpha)  
            with K_sat, permeability for saturated media
            and S, saturation field

           TAB2 indexes  are :

      .'ALPHA' (type FLOTTANT or CHPOINT) addimentional coefficient 
      .'PERMSAT' (type FLOTTANT or CHPOINT) permeability for saturated media


     Example :                                                             
     _______                                                             

   This user procedure K_PRO have differents permeability law for each 
   part of the modeled media

********************************************
'DEBPROC' KR_PRO TAB1*'TABLE' PRECED*'TABLE' ;
*---- part 'SITE'
  si (EGA PRECED.'NOMZONE' 'SITE') ;  
    K1 = 'KOPS' TAB1.'SATURATION_PROV' '**' PRECED.'ALPHA' ;
    K1 = 'KOPS' K1 '*' PRECED.'PERMSAT' ;
  finsi ;
*---- part 'BO'
  si (EGA PRECED.'NOMZONE' 'BO') ;  
    K1 = (kops TAB1.'PNS_PROV' '*' PRECED.'ALPHA') EXP  ; 
    K1 = kops K1 '*' PRECED.'PERMSAT' ;
  finsi ;
'FINPROC' K1 ;
********************************************

 it can be used by DARCYSAT procedure if parameters of each  
 permeability law are entered in the table SATUR in following way :

*---- permeability definition
SATUR.'LOI_PERMEABILITE' = TABLE 'MULTIZONE' ;
*-  site zone 
SATUR.'LOI_PERMEABILITE'. 'SITE' = TABLE 'STANDARD' ;
SATUR.'LOI_PERMEABILITE'. 'SITE'. 'ALPHA' = ... ;
SATUR.'LOI_PERMEABILITE'. 'SITE'. 'PERMSAT'= ... ;
SATUR.'LOI_PERMEABILITE'. 'SITE'. 'MODELE' = ... ;
*-  BO zone
SATUR.'LOI_PERMEABILITE'. 'BO' = TABLE 'PERSONNELLE' ;
SATUR.'LOI_PERMEABILITE'. 'BO'. 'ALPHA' = ... ;
SATUR.'LOI_PERMEABILITE'. 'BO'. 'PERMSAT'= ... ;
SATUR.'LOI_PERMEABILITE'. 'BO'. 'MODELE' = ... ;
... 
...
*----calculation
DARCYSAT SATUR