Accueil | Cast3M

   1 : $$$$ CHITRNSP NOTICE  CHAT      11/09/12    21:15:24     7124           
   2 :                                              DATE     11/09/12
   3 :                                              
   4 :   Procedure CHITRNSP                      Voir aussi : NOCOMCHI BMTD
   5 :      ------------------                                   CHI2      DMTD  
   6 :                                                           TRACHIT   DMMU  
   7 :                                                           TRACHIS   CHI1  
   8 :                                                           NOESPCHI  PECHE        
   9 :             CHITRNSP TAB1  ;
  10 : 
  11 : 
  12 :         Objet
  13 :         -----
  14 :         Cette procedure permet d'effectuer un calcul couple transport/
  15 :        chimie. Le transport utilise les elements finis mixte_hybrides.
  16 : 
  17 :         Commentaires
  18 :         -------------
  19 :     TAB1    est une table de soustype 'GEOCHIMIE'.
  20 :     En entree, TAB1 sert a definir les options et les parametres du
  21 :     calcul.
  22 :     En sortie TAB1 contient les donnees d'entrees et les resultats de
  23 :     façon a permettre une reprise du calcul.
  24 :     Les indices de la table TAB1 sont des mots (a ecrire en
  25 :     toutes lettres, et en majuscules s'ils sont mis entre cotes)
  26 :     dont voici la description :
  27 :                                                                       
  28 :     ------------------------------------------------                  
  29 :     Donnees physiques, geometriques et materielles :
  30 :     ------------------------------------------------                  
  31 : 
  32 :     Indices: 'SOUSTYPE' 'MODELE' 'DIFFUSION' <'POROSITE'> 'DOMAINE'
  33 :          'CONVECTION' 'CHIMI1' <'ITERC'> <'PRECISION'> <'DECROISSANCE'>
  34 : 
  35 :     'SOUSTYPE'     mot 'GEOCHIMIE'
  36 : 
  37 :     'MODELE'    Objet modele (MMODEL cree par MODE,formulation DARCY)
  38 :                                                                       
  39 :     'DIFFUSION' Donnees physiques et materielles :
  40 :               conductivite hydraulique (MCHAML cree par MATE)
  41 :                                                                       
  42 :     'POROSITE'    Contient la porosite au centre de l'element
  43 :               (CHPOIN de support DOMAINE.CENTRE)
  44 :                La valeur par defaut est 1.
  45 :                                                                       
  46 :     'DOMAINE'     References geometriques (TABLE creee par DOMA)      
  47 :                                                                       
  48 :     'CONVECTION'  Flux de la vitesse convective (CHAMPOIN de support  
  49 :                   DOMAINE.FACE)                                       
  50 :                                                                       
  51 :                                                                       
  52 :     'CHIMI1'  Table issue de CHI1
  53 :                                                                       
  54 :                                                                       
  55 :     'ITERC' nombre max d'iterations de couplage (defaut 100)
  56 :                                                                       
  57 :                                                                       
  58 :     'PRECISION' precision critere de convergence pour le couplage
  59 :                  valeur par defaut  1.E-3                             
  60 :                                                                       
  61 :     'DECROISSANCE' table TAB2 contenant les donnees relatives a la
  62 :           decroissance/filiation. Cette table est indicee par le mot
  63 :           'TETA' et des entiers de 1 a N. N etant le nombre de couples
  64 :            pere fils.
  65 :            TAB2.TETA est un reel le coefficient d'implicitation.
  66 :            TAB2.i    est une table d'indices 'PERE' 'FILS' et 'LAMBDA'
  67 :            TAB2.i .'PERE' est un entier identifiant du pere.
  68 :            TAB2.i .'FILS' est un entier identifiant du fils.
  69 :            TAB2.i .'LAMBDA' est un reel la constante de decroissance.
  70 : 
  71 :     --------------------------------------                            
  72 :        parametres de calcul de chimie
  73 :     --------------------------------------                            
  74 :                                                                       
  75 :                                                                       
  76 :      indices: <'EPS'> <'ITMAX'> <'ITERSOLI'> <'PRECPE'> <'IAFFICHE'>
  77 :            <'NITERPE'> <'DELPE'> <'MDELPE'> <'NFI'>  <'TEMPE'>
  78 :            <'CLIM'> <'SORTIE'> <'IMPRIM'>
  79 : 
  80 :     'EPS'  un REEL, la precision du calcul.
  81 :              Valeur par defaut 1.E-4.
  82 : 
  83 :     'ITMAX'  un ENTIER nombre maximal d'iterations dans la
  84 :              resolution du systeme chimique. Valeur par defaut 20.
  85 : 
  86 :     'ITERSOLI'  un ENTIER nombre maximal d'iterations, pour
  87 :              trouver les mineraux precipites. Valeur par defaut 10.
  88 : 
  89 :     'IAFFICHE'  un ENTIER permettant le choix d'affichage des
  90 :              resultats pour les solutions solides.
  91 :              1  coefficients stoechiometriques des solutions solides
  92 :              2  fractions molaires des solutions solides
  93 :              Valeur par defaut 2.
  94 : 
  95 :     'PRECPE'  un REEL, precision sur le calcul redox.
  96 :              Valeur par defaut 1.E-10
  97 : 
  98 :     'NITERPE'   un ENTIER nombre maximal d'iterations de
  99 :               dichotomie. Valeur par defaut 50.
 100 : 
 101 :     'DELPE' un REEL, l'intervalle initial des iterations de
 102 :                dichotomie. La valeur par defaut est 1.
 103 : 
 104 :     'MDELPE' un ENTIER nombre maximal de pas dans la recherche
 105 :              de l'intervalle de dichotomie. Valeur par defaut 20.
 106 :            ( evite de cycler lorsque l'on est tres loin de la solution)
 107 : 
 108 :     'NFI' un ENTIER nombre de cycles de chimie.
 109 :               Valeur par defaut 4. Un cycle correspond a la sequence:
 110 :                     * calcul de la force ionique
 111 :                     * modification des logk
 112 :                                            |---
 113 :                     * boucle mineraux a    |* resolution ( iterative )
 114 :                         precipiter         |
 115 :                                            |* verification des mineraux
 116 :                                            |  precipites
 117 :                                            |---
 118 : 
 119 :     'TEMPE'  objet de type CHPOIN contenant la temperature.
 120 : 
 121 :     'CLIM'  valeur de l'activite imposee des especes de type 3.
 122 :              Objet de type CHPOIN ayant une composante pour chaque
 123 :              espece dont l'activite est imposee.
 124 : 
 125 :     'SORTIE'  un objet de type LISTMOTS. Ces mots doivent
 126 :              etre pris dans la liste:
 127 :              'PREC' 'FION' 'TYP6' 'TYP3' 'NTY4' 'TYP5' 'SURF' 'SOLU'
 128 :              'POLE' 'LOGK'
 129 :              Ils servent a preciser les elements que l'on veut voir
 130 :              figurer dans la TABLE TAB1.
 131 : 
 132 :     'IMPRIM' un objet de type LISTENTI . Dans le cas ou l'on demande
 133 :              un niveau de message superieur a 0 ( OPTION IMPI 1 ),
 134 :              ceci permet de limiter les impressions aux seuls noeuds
 135 :              du maillage dont le numero figure dans la liste.
 136 : 
 137 :     --------------------------------------                            
 138 :     Conditions aux limites / chargements :                            
 139 :     --------------------------------------                            
 140 : 
 141 :     Indices : <'BLOCAGE'> <'TRACE_IMPOSE'> <'FLUX_IMPOSE'> <'SOURCE'>
 142 :                                                                       
 143 :     'BLOCAGE'     Contient les matrices de blocage (RIGIDITE)         
 144 :                                                                       
 145 :     'TRACE_IMPOSE' Valeurs des traces de concentrations imposees      
 146 :                  chargement de CHPOIN au centre des faces             
 147 :                                                                       
 148 :     'FLUX_IMPOSE' Valeurs des flux de concentrations imposees         
 149 :                  chargement de CHPOIN au centre des faces             
 150 :                                                                       
 151 :     'SOURCE'      Valeurs de la source                                
 152 :                   Chargement relatifs aux CHPOIN centre               
 153 :                   ( defaut  CHPOIN centre 0.)                         
 154 :                                                                       
 155 :     ---------------------
 156 :     Algorithme en temps :
 157 :     ---------------------
 158 :                                                                       
 159 :     Indices: 'PAS_DE_TEMPS' 'TEMPS_FINAL' <'TEMPS_SAUVES'> <'THETA'>
 160 : 
 161 :     'PAS_DE_TEMPS' Valeur du pas de temps de calcul (FLOTTANT)        
 162 :                                                                       
 163 :     'TEMPS_FINAL' Valeur du temps final (FLOTTANT)
 164 :                                                                       
 165 :     'TEMPS_SAUVES' Valeur des temps sauvegardes (LISTREEL - defaut    
 166 :                   on sauve tous les pas de temps)                     
 167 :                  Tous les temps de cette liste doivent etre des
 168 :                  multiples du pas de temps, en cas d'erreur on
 169 :                  sauve au moins le temps initial et le temps final.
 170 :                                                                       
 171 :     'THETA'       Coefficient de relaxation des termes de diffusion_
 172 :                   convection. Compris entre 0. et 1.
 173 :                   (theta-methode) (FLOTTANT - defaut 1.)              
 174 :                   Valeurs de theta generalement utilisees :           
 175 :                        Schema de Euler explicite : 0.                 
 176 :                        Schema de Crank-Nicholson : 1/2                
 177 :                        Schema de Euler implicite : 1.                 
 178 :                                                                       
 179 :   ----------------------------------------------
 180 :    Conditions initales et Resultats
 181 :   -------------------------------------------
 182 : 
 183 :         A chacun des indices suivants on trouve une table indicee par
 184 :         des entiers de 0 a N. L'indice i contient la valeur du champ au
 185 :         temps TEMPS.i.
 186 :         L'indice 0 contient les valeurs initiales,qui doivent etre
 187 :         donnees lors d'un premier calcul.
 188 :         Les noms des composantes des CHPOIN sont ceux utilises
 189 :         dans CHI2 (et attribues par le code). On peut obtenir les
 190 :         numeros des identifiants correspondants a l'aide des procedures
 191 :         NOCOMCHI et NOESPCHI.
 192 : 
 193 :         Indices indispensables en entree:  'TOT' 'TAQU' 'LOGC'
 194 : 
 195 :         Indices toujours presents en sortie: 'TOT' 'TAQU' 'LOGC'
 196 :                                       'TEMPS' 'AQUE' 'FIXE'  'FLUX'
 197 : 
 198 :         Indices optionnels en entree: 'TEMPS' 'FION' 'NTY4'
 199 : 
 200 :         Indices optionnels en sortie: <'PREC'>  <'FION'> <'TYP6'> 
 201 :         <'TYP3'> <'NTY4'> <'TYP5'> <'SURF'> <'SOLU'> <'POLE'> <'LOGK'>
 202 : 
 203 : 
 204 :     'TEMPS'   TABLE contenant les temps sauvegardes (FLOTTANT)
 205 :          La valeur par defaut lors d'un premier calcul est
 206 :          TAB1. 'TEMPS'. 0 = 0.
 207 :          Lors d'une reprise le temps initial est le temps de plus
 208 :          grand indice de cette table.
 209 : 
 210 :     'TOT' TABLE contenant les concentrations totales aux centres
 211 :          des elements. CHPOIN  ayant une composante par composant
 212 :          chimique. (La valeur initiale est indispensable).
 213 : 
 214 :     'AQUE' TABLE contenant les concentrations des aqueux aux centres
 215 :           des elements .CHPOIN ayant une composante par composant
 216 :           chimique ( cette table est generee par la procedure)
 217 : 
 218 :     'FIXE' TABLE contenant les concentrations des fixees aux centres
 219 :           des elements .CHPOIN ayant une composante par composant
 220 :           chimique ( cette table est generee par la procedure)
 221 : 
 222 :     'TAQU' Table contenant les concentrations des aqueux
 223 :           au centre des faces. CHPOIN possedant une composante
 224 :           par composant chimique.(La valeur initiale est
 225 :           indispensable)
 226 : 
 227 :     'LOGC' Table contenant pour  chaque composant chimique le log
 228 :           de la concentration de l'espece simple associee aux centres
 229 :           des elements.CHPOIN possedant une composante par composant
 230 :           chimique. (La valeur initiale est  indispensable)
 231 : 
 232 :     'FLUX' Table contenant les flux totaux (convectif+ diffusif)
 233 :           au centre des faces. CHPOIN possedant une composante
 234 :           par composant chimique.( cette table est generee par la
 235 :           procedure)
 236 : 
 237 : 
 238 :     'PREC' TABLE contenant la concentration des especes precipitees.
 239 :           CHPOIN centre ayant une composante par espece chimique
 240 :           precipite.
 241 : 
 242 :     'FION' TABLE  contenant la force ionique. CHPOIN centre a une
 243 :           composante ( nom de la composante SCAL)
 244 : 
 245 :     'TYP6' TABLE contenant la concentration des especes de type 6.
 246 :            CHPOIN centre ayant une composante par espece chimique
 247 :            de type 6.( cette table est generee par la procedure si TYP6
 248 :            figure dans la liste TAB1.SORTIE)
 249 : 
 250 :     'TYP3' TABLE contenant la concentration des especes de type 3.
 251 :            ( cette table est generee par la procedure si TYP3 figure
 252 :            dans la liste TAB1.SORTIE)
 253 : 
 254 :     'NTY4' TABLE contenant des CHPOIN centre ayant une composante
 255 :           pour chaque espece de precipite potentiel.
 256 :           En chaque point centre d'element on indiquera si le mineral
 257 :           est precipite ( =1) ou non( =0). Sert a initialiser le choix
 258 :           des mineraux.
 259 : 
 260 :     'TYP5' TABLE contenant la concentration des especes de type 5.
 261 :            CHPOIN centre ayant une composante par espece chimique
 262 :            de type 5.( cette table est generee par la procedure si TYP5
 263 :            figure dans la liste TAB1.SORTIE)
 264 : 
 265 :     'SURF' TABLE contenant la concentration des especes de surface.
 266 :            CHPOIN centre ayant une composante par espece chimique
 267 :            de surface.( cette table est generee par la procedure si
 268 :            SURF figure dans la liste TAB1.SORTIE)
 269 : 
 270 :     'SOLU' TABLE contenant la concentration des especes en solution.
 271 :            CHPOIN centre ayant une composante par espece chimique
 272 :            en solution.( cette table est generee par la procedure si
 273 :            SOLU figure dans la liste TAB1.SORTIE)
 274 : 
 275 :     'POLE' TABLE contenant les fractions molaires des poles de
 276 :            solutions solides.CHPOIN centre ayant une composante par
 277 :            pole de solution solide.( cette table est generee par la
 278 :            procedure si POLE figure dans la liste TAB1.SORTIE)
 279 : 
 280 :     'LOGK' TABLE contenant les logK des solutions solides.
 281 :            CHPOIN centre ayant une composante par solution solide.
 282 :            ( cette table est generee par la procedure si LOGK
 283 :            figure dans la liste TAB1.SORTIE)
 284 : 
 285 : 
 286 : 
 287 : 
 288 : 
 289 :  
 290 : 
 291 :  
 292 : 
 293 :  
 294 :  
 295 :  
 296 :  
 297 :  
 298 : 
 299 : 
 300 :  
 301 :  
 302 :  
 303 :