Opérateur CHITRNSP DATE 98/12/29 Voir aussiCHI1 BMTD CHI2 DMTD TRACHIT DMMU TRACHIS NOESPCHI PECHE NOCOMCHI Procédure CHITRNSP ------------------ CHITRNSP TAB1 ; Objet ----- Cette procédure permet d'effectuer un calcul couplé transport/ chimie. Le transport utilise les éléments finis mixte_hybrides. Commentaires ------------- TAB1 est une table de soustype 'GEOCHIMIE'. En entrée, TAB1 sert à définir les options et les paramètres du calcul. En sortie TAB1 contient les données d'entrées et les résultats de façon à permettre une reprise du calcul. Les indices de la table TAB1 sont des mots (à écrire en toutes lettres, et en majuscules s'ils sont mis entre cotes) dont voici la description : ------------------------------------------------ Données physiques, géométriques et matérielles : ------------------------------------------------ Indices: 'SOUSTYPE' 'MODELE' 'DIFFUSION' <'POROSITE'> 'DOMAINE' 'CONVECTION' 'CHIMI1' <'ITERC'> <'PRECISION'> <'DECROISSANCE'> 'SOUSTYPE' mot 'GEOCHIMIE' 'MODELE' Objet modèle (MMODEL créé par MODE,formulation DARCY) 'DIFFUSION' Données physiques et matérielles : conductivité hydraulique (MCHAML créé par MATE) 'POROSITE' Contient la porosité au centre de l'élément (CHPOIN de support DOMAINE.CENTRE) La valeur par défaut est 1. 'DOMAINE' Références géométriques (TABLE créée par DOMA) 'CONVECTION' Flux de la vitesse convective (CHAMPOIN de support DOMAINE.FACE) 'CHIMI1' Table issue de CHI1 'ITERC' nombre max d'itérations de couplage (défaut 100) 'PRECISION' précision critère de convergence pour le couplage valeur par défaut 1.E-3 'DECROISSANCE' table TAB2 contenant les données relatives à la décroissance/filiation. Cette table est indicée par le mot 'TETA' et des entiers de 1 à N. N étant le nombre de couples père fils. TAB2.TETA est un réel le coefficient d'implicitation. TAB2.i est une table d'indices 'PERE' 'FILS' et 'LAMBDA' TAB2.i .'PERE' est un entier identifiant du père. TAB2.i .'FILS' est un entier identifiant du fils. TAB2.i .'LAMBDA' est un réel la constante de décroissance. -------------------------------------- paramètres de calcul de chimie -------------------------------------- indices: <'EPS'> <'ITMAX'> <'ITERSOLI'> <'PRECPE'> <'IAFFICHE'> <'NITERPE'> <'DELPE'> <'MDELPE'> <'NFI'> <'TEMPE'> <'CLIM'> <'SORTIE'> <'IMPRIM'> 'EPS' un REEL, la précision du calcul. Valeur par défaut 1.E-4. 'ITMAX' un ENTIER nombre maximal d'itérations dans la résolution du système chimique. Valeur par défaut 20. 'ITERSOLI' un ENTIER nombre maximal d'itérations, pour trouver les minéraux précipités. Valeur par défaut 10. 'IAFFICHE' 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. 'PRECPE' un REEL, précision sur le calcul redox. Valeur par défaut 1.E-10 'NITERPE' un ENTIER nombre maximal d'itérations de dichotomie. Valeur par défaut 50. 'DELPE' un REEL, l'intervalle initial des itérations de dichotomie. La valeur par défaut est 1. 'MDELPE' 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 trés loin de la solution) 'NFI' 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 |--- 'TEMPE' objet de type CHPOIN contenant la température. '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. 'SORTIE' un objet de type LISTMOTS. Ces mots doivent être pris dans la liste: 'PREC' 'FION' 'TYP6' 'TYP3' 'NTY4' 'TYP5' 'SURF' 'SOLU' 'POLE' 'LOGK' Ils servent à préciser les éléments que l'on veut voir figurer dans la TABLE TAB1. 'IMPRIM' 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. -------------------------------------- Conditions aux limites / chargements : -------------------------------------- Indices : <'BLOCAGE'> <'TRACE_IMPOSE'> <'FLUX_IMPOSE'> <'SOURCE'> 'BLOCAGE' Contient les matrices de blocage (RIGIDITE) 'TRACE_IMPOSE' Valeurs des traces de concentrations imposées chargement de CHPOIN au centre des faces 'FLUX_IMPOSE' Valeurs des flux de concentrations imposées chargement de CHPOIN au centre des faces 'SOURCE' Valeurs de la source Chargement relatifs aux CHPOIN centre ( défaut CHPOIN centre 0.) --------------------- Algorithme en temps : --------------------- Indices: 'PAS_DE_TEMPS' 'TEMPS_FINAL' <'TEMPS_SAUVES'> <'THETA'> 'PAS_DE_TEMPS' Valeur du pas de temps de calcul (FLOTTANT) 'TEMPS_FINAL' Valeur du temps final (FLOTTANT) 'TEMPS_SAUVES' Valeur des temps sauvegardés (LISTREEL - défaut on sauve tous les pas de temps) Tous les temps de cette liste doivent être des multiples du pas de temps, en cas d'erreur on sauve au moins le temps initial et le temps final. 'THETA' Coefficient de relaxation des termes de diffusion_ convection. Compris entre 0. et 1. (theta-méthode) (FLOTTANT - défaut 1.) Valeurs de theta généralement utilisées : Schéma de Euler explicite : 0. Schéma de Crank-Nicholson : 1/2 Schéma de Euler implicite : 1. ---------------------------------------------- Conditions initales et Résultats ------------------------------------------- A chacun des indices suivants on trouve une table indicée par des entiers de 0 à N. L'indice i contient la valeur du champ au temps TEMPS.i. L'indice 0 contient les valeurs initiales,qui doivent être données lors d'un premier calcul. Les noms des composantes des CHPOIN sont ceux utilisés dans CHI2 (et attribués par le code). On peut obtenir les numéros des identifiants correspondants à l'aide des procédures NOCOMCHI et NOESPCHI. Indices indispensables en entrée: 'TOT' 'TAQU' 'LOGC' Indices toujours présents en sortie: 'TOT' 'TAQU' 'LOGC' 'TEMPS' 'AQUE' 'FIXE' 'FLUX' Indices optionnels en entrée: 'TEMPS' 'FION' 'NTY4' Indices optionnels en sortie: <'PREC'> <'FION'> <'TYP6'> <'TYP3'> <'NTY4'> <'TYP5'> <'SURF'> <'SOLU'> <'POLE'> <'LOGK'> 'TEMPS' TABLE contenant les temps sauvegardés (FLOTTANT) La valeur par défaut lors d'un premier calcul est TAB1. 'TEMPS'. 0 = 0. Lors d'une reprise le temps initial est le temps de plus grand indice de cette table. 'TOT' TABLE contenant les concentrations totales aux centres des éléments. CHPOIN ayant une composante par composant chimique. (La valeur initiale est indispensable). 'AQUE' TABLE contenant les concentrations des aqueux aux centres des éléments .CHPOIN ayant une composante par composant chimique ( cette table est générée par la procédure) 'FIXE' TABLE contenant les concentrations des fixées aux centres des éléments .CHPOIN ayant une composante par composant chimique ( cette table est générée par la procédure) 'TAQU' Table contenant les concentrations des aqueux au centre des faces. CHPOIN possédant une composante par composant chimique.(La valeur initiale est indispensable) 'LOGC' Table contenant pour chaque composant chimique le log de la concentration de l'espèce simple associée aux centres des éléments.CHPOIN possédant une composante par composant chimique. (La valeur initiale est indispensable) 'FLUX' Table contenant les flux totaux (convectif+ diffusif) au centre des faces. CHPOIN possédant une composante par composant chimique.( cette table est générée par la procédure) 'PREC' TABLE contenant la concentration des espèces précipitées. CHPOIN centre ayant une composante par espèce chimique précipité. 'FION' TABLE contenant la force ionique. CHPOIN centre à une composante ( nom de la composante SCAL) 'TYP6' TABLE contenant la concentration des espèces de type 6. CHPOIN centre ayant une composante par espèce chimique de type 6.( cette table est générée par la procédure si TYP6 figure dans la liste TAB1.SORTIE) 'TYP3' TABLE contenant la concentration des espèces de type 3. ( cette table est générée par la procédure si TYP3 figure dans la liste TAB1.SORTIE) 'NTY4' TABLE contenant des CHPOIN centre ayant une composante pour chaque espèce de précipité potentiel. En chaque point centre d'élément on indiquera si le minéral est précipité ( =1) ou non( =0). Sert à initialiser le choix des minéraux. 'TYP5' TABLE contenant la concentration des espèces de type 5. CHPOIN centre ayant une composante par espèce chimique de type 5.( cette table est générée par la procédure si TYP5 figure dans la liste TAB1.SORTIE) 'SURF' TABLE contenant la concentration des espèces de surface. CHPOIN centre ayant une composante par espèce chimique de surface.( cette table est générée par la procédure si SURF figure dans la liste TAB1.SORTIE) 'SOLU' TABLE contenant la concentration des espèces en solution. CHPOIN centre ayant une composante par espèce chimique en solution.( cette table est générée par la procédure si SOLU figure dans la liste TAB1.SORTIE) 'POLE' TABLE contenant les fractions molaires des poles de solutions solides.CHPOIN centre ayant une composante par pole de solution solide.( cette table est générée par la procédure si POLE figure dans la liste TAB1.SORTIE) 'LOGK' TABLE contenant les logK des solutions solides. CHPOIN centre ayant une composante par solution solide. ( cette table est générée par la procédure si LOGK figure dans la liste TAB1.SORTIE)