Opérateur CHI2 DATE 97/12/30 Opérateur CHI2 -------------- TAB4 = CHI2 TAB1 TAB2 <TAB3> ; Objet ------- Le but est de calculer la spéciation d'une eau, en tout point d'un domaine à partir de la donnée des concentrations analytiques de chaque composant chimique du système. Le calcul se fait en deux temps à l'aide des opérateurs CHI1 et CHI2. CHI1 rassemble toutes les données relatives à un système chimique, et CHI2 effectue la spéciation. La terminologie est celle de Mineql. Toutes les concentrations sont données en moles par litre. Commentaires ------------ TAB1 est un objet de type TABLE et de sous type chimi1 (cf opérateur CHI1) TAB2 est un objet de type TABLE de sous-type 'DONNEES_CHIMIQUES'. Elle est indicée par les mots : 'LOGC' 'TOT' <'FIONI'> <'NTY4'> <'TEMPE'> <'CLIM'> TAB2.LOGC est un objet de type CHPOIN qui possède une composante par composant chimique. Pour chaque composant chimique il contiendra le log de l'estimation de l'espèce simple associée. Le nom de ces composantes est un mot de 4 caractères, formé par X suivi éventuellement de 0 ou 00 et du numéro identifiant le composant chimique. TAB2.TOT est un objet de type CHPOIN qui possède une composante par composant chimique. Pour chaque composant chimique il contiendra la concentration totale (ou analytique) ( en solution + minéraux). Le nom de ces composantes est un mot de 4 caractères, formé par X suivi éventuellement de 0 ou 00 et du numéro identifiant le composant chimique. TAB2.FIONI objet de type CHPOIN ayant une composante scalaire, et contenant une estimation de la force ionique en chaque point du maillage. TAB2.NTY4 objet de type CHPOIN ayant une composante pour chaque espèce de précipité potentiel. En chaque point du maillage on indiquera si le minéral est précipité ( =1) ou non( =0). Les nom des composantes sont ceux figurant dans la liste TAB1.IDEN.NOMPRECI. TAB2.TEMPE objet de type CHPOIN contenant la température. TAB2.CLIM valeur de l'activité imposée des espèces de type 3. Objet de type CHPOIN ayant une composante pour chaque espèce dont l'activité est imposée.Les noms des composantes sont ceux figurant dans la liste TAB1.IDEN.NOMTYP3. TAB3 est un objet de type TABLE contenant les paramètres de calcul. Elle est indicée par les mots : <'EPS'> <'ITMAX'> <'ITERSOLI'> <'PRECPE'> <'IAFFICHE'> <'NITERPE'> <'DELPE'> <'MDELPE'> <'NFI'> <'SORTIE'> <'IMPRIM'> TAB3.EPS est un REEL, la précision du calcul. Valeur par défaut 1.E-4. TAB3.ITMAX est un ENTIER nombre maximal d'itérations dans la résolution du système chimique. Valeur par défaut 20. TAB3.ITERSOLI est un ENTIER nombre maximal d'itérations, pour trouver les minéraux précipités. Valeur par défaut 10. TAB3.IAFFICHE est un ENTIER permettant le choix d'affichage des résultats pour les solutions solides. 1 coefficients stoechiométriques des solutions solides 2 fractions molaires des solutions solides Valeur par défaut 2. TAB3.PRECPE est un REEL, précision sur le calcul redox. Valeur par défaut 1.E-10 TAB3.NITERPE est un ENTIER nombre maximal d'itérations de dichotomie. Valeur par défaut 50. TAB3.DELPE est un REEL, l'intervalle initial des itérations de dichotomie. La valeur par défaut est 1. TAB3.MDELPE est un ENTIER nombre maximal de pas dans la recherche de l'intervalle de dichotomie. Valeur par défaut 20. ( évite de cycler lorsque l'on est tres loin de la solution) TAB3.NFI est un ENTIER nombre de cycles de chimie. Valeur par défaut 4. Un cycle correspond à la séquence: * calcul de la force ionique * modification des logk |--- * boucle minéraux à |* résolution ( itérative ) précipiter | |* vérification des minéraux | précipités |--- TAB3.SORTIE est un objet de type LISTMOTS. Ces mots doivent etre pris dans la liste: 'PREC' 'FION' 'TYP6' 'TYP3' 'NTY4' 'TYP5' 'SURF' 'SOLU' 'POLE' 'LOGK' Ils servent à préciser les élements que l'on veut voir figurer dans la TABLE TAB4. TAB3.IMPRIM est un objet de type LISTENTI . Dans le cas où l'on demande un niveau de message supérieur à 0 ( OPTION IMPI 1 ), ceci permet de limiter les impressions aux seuls noeuds du maillage dont le numéro figure dans la liste. TAB4 objet de type TABLE de sous type CHIMI2, contenant la spéciation du système. Elle est indicée par les mots: 'AQUE' 'FIXE' 'LOGC' <'PREC'> <'FION'> <'TYP6'> <'TYP3'> <'NTY4'> <'TYP5'> <'SURF'> <'SOLU'> <'POLE'> <'LOGK'> TAB4.AQUE est un objet de type CHPOIN qui a une composante par composant chimique. Le nom de ces composantes est un mot de 4 caractères, formé par, X suivi éventuellement de 0 ou 00 et du numéro identifiant le composant chimique. Pour chaque composant il contiendra la concentration en solution . TAB4.FIXE est un objet de type CHPOIN qui a une composante par composant chimique. Pour chaque composant il contiendra la concentration du minéral précipité+surface. TAB4.LOGC est un objet de type CHPOIN qui possède une composante par composant chimique. Pour chaque composant chimique il contiendra le log de la concentration de l'espèce simple associée. TAB4.PREC est un objet de type CHPOIN qui a une composante par espèce chimique précipité. Le nom de ces composantes sont ceux figurant dans la liste TAB1.IDEN.NOMPRECI. Il contiendra la concentration des espèces précipitées. TAB4.FION est un objet de type CHPOIN contenant la force ionique TAB4.TYP6 est un objet de type CHPOIN qui a une composante par espèce chimique de type 6. Le nom de ces composantes sont ceux figurant dans la liste TAB1.IDEN.NOMTYP6. Il contiendra la concentration des espèces de type 6, c'est à dire non prises en compte dans le problème. TAB4.TYP3 est un objet de type CHPOIN qui a une composante par espèce chimique de type 3. Le nom de ces composantes sont ceux figurant dans la liste TAB1.IDEN.NOMTYP3. Il contiendra la concentration des espèces de type 3. C'est à dire la quantité totale nécéssaire de chaque espèce pour que la contrainte soit vérifiée. Par exemple si on veut imposer un pH de 8 ,la concentration trouvée en type 3 donne l'opposé de la quantite d H+ nécéssaire pour avoir un pH de 8. TAB4.NTY4 objet de type CHPOIN ayant une composante pour chaque espèce de précipité potentiel. En chaque point du maillage on indiquera si le minéral est précipité ( =1) ou non( =0). Les nom des composantes sont ceux figurant dans la liste TAB1.IDEN.NOMPRECI. TAB4.TYP5 est un objet de type CHPOIN qui a une composante par espèce chimique de type 5. Il contiendra la concentration des espèces de type 5. Les espèces de type 5 sont les précipités potentiels qui ne sont pas précipités. TAB4.SURF est un objet de type CHPOIN qui a une composante par espèce chimique de surface. Il contiendra la concentration des espèces de surface. TAB4.SOLU est un objet de type CHPOIN qui a une composante par espèce chimique en solution. Il contiendra la concentration des espèces en solution. TAB4.POLE est un objet de type CHPOIN qui a une composante par pole de solution solide. Le nom de ces composantes sont ceux figurant dans la liste TAB1.IDEN.NOMPOLE. Il contiendra les fractions molaires des poles de solutions solides. TAB4.LOGK est un objet de type CHPOIN qui a une composante par solution solide. Le nom de ces composantes sont ceux figurant dans la liste TAB1.IDEN.NOMSOSO. Il contiendra les logK des solutions solides.