Accueil | Cast3M

   1 : $$$$ CHI2     NOTICE  CHAT      11/09/12    21:15:24     7124           
   2 :                                              DATE     11/09/12
   3 :   
   4 :    Operateur CHI2                                Voir aussi : 
   5 :    --------------  
   6 :        TAB4 = CHI2  TAB1  TAB2  <TAB3> ;
   7 : 
   8 : 
   9 :         Objet
  10 :        -------
  11 :        Le but est de calculer la speciation d'une eau, en tout point
  12 :        d'un domaine a partir de la donnee des concentrations analytiques
  13 :        de chaque composant chimique du systeme. Le calcul se fait en
  14 :        deux temps a l'aide des operateurs CHI1 et CHI2.
  15 :        CHI1 rassemble toutes les donnees relatives a un systeme chimique,
  16 :        et CHI2 effectue la speciation.
  17 :         La terminologie est celle de Mineql.
  18 : 
  19 :           Toutes les concentrations sont donnees en moles par litre.
  20 : 
  21 :        Commentaires
  22 :        ------------
  23 :        TAB1 est un objet de type TABLE et de sous type chimi1 
  24 :             (cf operateur CHI1)
  25 : 
  26 : 
  27 :        TAB2 est un objet de type TABLE de sous-type 'DONNEES_CHIMIQUES'.
  28 :             Elle est indicee par les mots :
  29 :            'LOGC' 'TOT' <'FIONI'> <'NTY4'> <'TEMPE'>  <'CLIM'>
  30 : 
  31 :        TAB2.LOGC est un objet de type CHPOIN qui possede une composante
  32 :              par composant chimique. Pour chaque composant chimique il 
  33 :              contiendra le log de l'estimation de l'espece simple
  34 :              associee. Le nom de ces composantes est un mot de 4 
  35 :              caracteres, forme par X suivi eventuellement de 0 ou 00 et
  36 :              du numero identifiant le composant chimique.
  37 : 
  38 :        TAB2.TOT est un objet de type CHPOIN qui possede une composante
  39 :              par composant chimique. Pour chaque composant chimique il
  40 :              contiendra la concentration totale  (ou analytique) ( en
  41 :              solution + mineraux). Le nom de ces composantes est un mot
  42 :              de 4 caracteres, forme par X suivi eventuellement de 0 ou 
  43 :              00 et du numero identifiant le composant chimique.
  44 : 
  45 :        TAB2.FIONI objet de type CHPOIN ayant une composante scalaire, et
  46 :              contenant une estimation de la force ionique en chaque 
  47 :              point du maillage.
  48 : 
  49 :        TAB2.NTY4 objet de type CHPOIN ayant une composante pour chaque
  50 :              espece de precipite potentiel. En chaque point du maillage
  51 :              on indiquera si le mineral est precipite ( =1) ou non( =0).
  52 :              Les nom des composantes sont ceux figurant dans la liste
  53 :              TAB1.IDEN.NOMPRECI.
  54 : 
  55 :        TAB2.TEMPE objet de type CHPOIN contenant la temperature.
  56 : 
  57 :        TAB2.CLIM valeur de l'activite imposee des especes de type 3.
  58 :              Objet de type CHPOIN ayant une composante pour chaque
  59 :              espece dont l'activite est imposee.Les noms des composantes
  60 :              sont ceux figurant dans la liste TAB1.IDEN.NOMTYP3.
  61 : 
  62 :        TAB3 est un objet de type TABLE contenant les parametres de
  63 :             calcul. Elle est indicee par les mots :
  64 :             <'EPS'> <'ITMAX'> <'ITERSOLI'> <'PRECPE'> <'IAFFICHE'>
  65 :            <'NITERPE'> <'DELPE'> <'MDELPE'>  <'NFI'>  <'SORTIE'>
  66 :            <'IMPRIM'>
  67 : 
  68 :        TAB3.EPS  est un REEL, la precision du calcul.
  69 :              Valeur par defaut 1.E-4.
  70 : 
  71 :        TAB3.ITMAX est un ENTIER nombre maximal d'iterations dans la
  72 :              resolution du systeme chimique. Valeur par defaut 20.
  73 : 
  74 :        TAB3.ITERSOLI est un ENTIER nombre maximal d'iterations, pour
  75 :              trouver les mineraux precipites. Valeur par defaut 10.
  76 : 
  77 :        TAB3.IAFFICHE est un ENTIER permettant le choix d'affichage des
  78 :              resultats pour les solutions solides.
  79 :              1  coefficients stoechiometriques des solutions solides
  80 :              2  fractions molaires des solutions solides
  81 :              Valeur par defaut 2.
  82 : 
  83 : 
  84 :        TAB3.CALCLOG est un ENTIER, indication sur le calcul en log
  85 :               de concentration ou pas dans la speciation 
  86 :              0 calcul en concentration     
  87 :              1 calcul en log de concentration
  88 :              Valeur par defaut 0
  89 :              
  90 :        TAB3.PRECPE est un REEL, precision sur le calcul redox.
  91 :              Valeur par defaut 1.E-10
  92 : 
  93 :        TAB3.NITERPE   est un ENTIER nombre maximal d'iterations de
  94 :               dichotomie. Valeur par defaut 50.
  95 : 
  96 :        TAB3.DELPE est un REEL, l'intervalle initial des iterations de
  97 :                dichotomie. La valeur par defaut est 1.
  98 : 
  99 :        TAB3.MDELPE est un ENTIER nombre maximal de pas dans la recherche
 100 :              de l'intervalle de dichotomie. Valeur par defaut 20.
 101 :             ( evite de cycler lorsque l'on est tres loin de la solution)
 102 : 
 103 :        TAB3.NFI est un ENTIER nombre de cycles de chimie.
 104 :               Valeur par defaut 4. Un cycle correspond a la sequence:
 105 :                     * calcul de la force ionique
 106 :                     * modification des logk
 107 :                                            |---
 108 :                     * boucle mineraux a    |* resolution ( iterative )
 109 :                         precipiter         |
 110 :                                            |* verification des mineraux
 111 :                                            |  precipites
 112 :                                            |---
 113 : 
 114 :        TAB3.SORTIE est un objet de type LISTMOTS. Ces mots doivent
 115 :              etre pris dans la liste:
 116 :              'PREC' 'FION' 'TYP6' 'TYP3' 'NTY4' 'TYP5' 'SURF' 'SOLU'
 117 :              'POLE' 'LOGK'
 118 :              Ils servent a preciser les elements que l'on veut voir
 119 :              figurer dans la TABLE TAB4.
 120 : 
 121 :        TAB3.IMPRIM est un objet de type LISTENTI . Dans le cas ou
 122 :              l'on demande un niveau de message superieur a 0
 123 :              ( OPTION IMPI 1 ), ceci permet de limiter les impressions
 124 :              aux seuls noeuds du maillage dont le numero figure dans 
 125 :              la liste.
 126 : 
 127 : 
 128 : 
 129 : 
 130 :        TAB4 objet de type TABLE de sous type CHIMI2, contenant la
 131 :             speciation du systeme. Elle est indicee par les mots:
 132 :             'AQUE' 'FIXE' 'LOGC' 'ERR' <'PREC'>  <'FION'> <'TYP6'>
 133 :               <'TYP3'> <'NTY4'> <'TYP5'> <'SURF'> <'SOLU'> <'POLE'>
 134 :               <'LOGK'> 
 135 : 
 136 :        TAB4.AQUE est un objet de type CHPOIN qui a une composante par
 137 :              composant chimique. Le nom de ces composantes est un mot 
 138 :              de 4 caracteres, forme par, X suivi eventuellement de 0
 139 :              ou 00 et du numero identifiant le composant chimique.
 140 :              Pour chaque composant il contiendra la concentration  
 141 :              en solution .
 142 : 
 143 :        TAB4.FIXE est un objet de type CHPOIN qui a une composante par
 144 :              composant chimique. Pour chaque composant il contiendra la
 145 :              concentration du mineral precipite+surface.
 146 : 
 147 :        TAB4.LOGC est un objet de type CHPOIN qui possede une composante
 148 :              par composant chimique. Pour chaque composant chimique il 
 149 :              contiendra le log de la concentration de l'espece simple
 150 :              associee.
 151 : 
 152 :        TAB4.PREC est un objet de type CHPOIN qui a une composante par
 153 :               espece chimique precipite. Le nom de ces composantes sont
 154 :               ceux figurant dans la liste TAB1.IDEN.NOMPRECI.
 155 :               Il contiendra la concentration des especes precipitees.
 156 : 
 157 :        TAB4.ERR est un objet de type CHPOIN contenant les erreurs
 158 :                 eventuellement survenue lors du calcul
 159 :                 ( 0 : resultats ok,
 160 :                   1 : nb d'iteration de Newton superieur a ITMAX,
 161 :                   7 : violation de la regle des phases,
 162 :                   8 : pb dans la resolutionn du systeme lineaire
 163 :                      dans le newton)
 164 :        
 165 :        TAB4.FION est un objet de type CHPOIN contenant la force ionique
 166 : 
 167 :        TAB4.TYP6 est un objet de type CHPOIN qui a une composante par
 168 :              espece chimique de type 6. Le nom de ces composantes sont
 169 :              ceux figurant dans la liste TAB1.IDEN.NOMTYP6.
 170 :              Il contiendra la concentration des especes de type 6,
 171 :              c'est a dire non prises en compte dans le probleme.
 172 : 
 173 :        TAB4.TYP3 est un objet de type CHPOIN qui a une composante par
 174 :              espece chimique de type 3. Le nom de ces composantes sont
 175 :              ceux figurant dans la liste TAB1.IDEN.NOMTYP3.
 176 :               Il contiendra la concentration des especes de type 3.
 177 :              C'est a dire la quantite totale necessaire de chaque espece
 178 :              pour que la contrainte soit verifiee.
 179 :              Par exemple si on veut imposer un pH de 8 ,la concentration
 180 :              trouvee en type 3 donne l'oppose de la quantite d H+
 181 :              necessaire pour avoir un pH de 8.
 182 : 
 183 :        TAB4.NTY4 objet de type CHPOIN ayant une composante pour chaque
 184 :              espece de precipite potentiel. En chaque point du maillage
 185 :              on indiquera si le mineral est precipite ( =1) ou non( =0).
 186 :              Les nom des composantes sont ceux figurant dans la liste
 187 :              TAB1.IDEN.NOMPRECI.
 188 : 
 189 :        TAB4.TYP5 est un objet de type CHPOIN qui a une composante par
 190 :              espece chimique de type 5.
 191 :               Il contiendra la concentration des especes de type 5.
 192 :               Les especes de type 5 sont les precipites potentiels qui
 193 :               ne sont pas precipites.
 194 : 
 195 :        TAB4.SURF est un objet de type CHPOIN qui a une composante par
 196 :              espece chimique de surface.
 197 :              Il contiendra la concentration des especes de surface.
 198 : 
 199 :        TAB4.SOLU est un objet de type CHPOIN qui a une composante par
 200 :              espece chimique en solution.
 201 :              Il contiendra la concentration des especes en solution.
 202 : 
 203 :        TAB4.POLE est un objet de type CHPOIN qui a une composante par
 204 :              pole de solution solide. Le nom de ces composantes sont 
 205 :              ceux figurant dans la liste TAB1.IDEN.NOMPOLE. Il contiendra 
 206 :              les fractions molaires des poles de solutions solides.
 207 : 
 208 :        TAB4.LOGK est un objet de type CHPOIN qui a une composante par
 209 :              solution solide. Le nom de ces composantes sont ceux 
 210 :              figurant dans la liste TAB1.IDEN.NOMSOSO.
 211 :              Il contiendra les logK des solutions solides.
 212 : 
 213 : 
 214 :  
 215 : 
 216 :  
 217 : 
 218 :  
 219 :  
 220 :  
 221 :  
 222 :  
 223 :  
 224 :  
 225 :