Afficher cette notice en
$$$$ CHITRNSP NOTICE CHAT 11/09/12 21:15:24 7124 DATE 11/09/12 Procedure CHITRNSP Voir aussi : NOCOMCHI BMTD ------------------ CHI2 DMTD TRACHIT DMMU TRACHIS CHI1 NOESPCHI PECHE CHITRNSP TAB1 ; Objet ----- Cette procedure permet d'effectuer un calcul couple transport/ chimie. Le transport utilise les elements finis mixte_hybrides. Commentaires ------------- TAB1 est une table de soustype 'GEOCHIMIE'. En entree, TAB1 sert a definir les options et les parametres du calcul. En sortie TAB1 contient les donnees d'entrees et les resultats de façon a permettre une reprise du calcul. Les indices de la table TAB1 sont des mots (a ecrire en toutes lettres, et en majuscules s'ils sont mis entre cotes) dont voici la description : ------------------------------------------------ Donnees physiques, geometriques et materielles : ------------------------------------------------ Indices: 'SOUSTYPE' 'MODELE' 'DIFFUSION' <'POROSITE'> 'DOMAINE' 'CONVECTION' 'CHIMI1' <'ITERC'> <'PRECISION'> <'DECROISSANCE'> 'SOUSTYPE' mot 'GEOCHIMIE' 'MODELE' Objet modele (MMODEL cree par MODE,formulation DARCY) 'DIFFUSION' Donnees physiques et materielles : conductivite hydraulique (MCHAML cree par MATE) 'POROSITE' Contient la porosite au centre de l'element (CHPOIN de support DOMAINE.CENTRE) La valeur par defaut est 1. 'DOMAINE' References geometriques (TABLE creee par DOMA) 'CONVECTION' Flux de la vitesse convective (CHAMPOIN de support DOMAINE.FACE) 'CHIMI1' Table issue de CHI1 'ITERC' nombre max d'iterations de couplage (defaut 100) 'PRECISION' precision critere de convergence pour le couplage valeur par defaut 1.E-3 'DECROISSANCE' table TAB2 contenant les donnees relatives a la decroissance/filiation. Cette table est indicee par le mot 'TETA' et des entiers de 1 a N. N etant le nombre de couples pere fils. TAB2.TETA est un reel le coefficient d'implicitation. TAB2.i est une table d'indices 'PERE' 'FILS' et 'LAMBDA' TAB2.i .'PERE' est un entier identifiant du pere. TAB2.i .'FILS' est un entier identifiant du fils. TAB2.i .'LAMBDA' est un reel la constante de decroissance. -------------------------------------- parametres de calcul de chimie -------------------------------------- indices: <'EPS'> <'ITMAX'> <'ITERSOLI'> <'PRECPE'> <'IAFFICHE'> <'NITERPE'> <'DELPE'> <'MDELPE'> <'NFI'> <'TEMPE'> <'CLIM'> <'SORTIE'> <'IMPRIM'> 'EPS' un REEL, la precision du calcul. Valeur par defaut 1.E-4. 'ITMAX' un ENTIER nombre maximal d'iterations dans la resolution du systeme chimique. Valeur par defaut 20. 'ITERSOLI' un ENTIER nombre maximal d'iterations, pour trouver les mineraux precipites. Valeur par defaut 10. 'IAFFICHE' un ENTIER permettant le choix d'affichage des resultats pour les solutions solides. 1 coefficients stoechiometriques des solutions solides 2 fractions molaires des solutions solides Valeur par defaut 2. 'PRECPE' un REEL, precision sur le calcul redox. Valeur par defaut 1.E-10 'NITERPE' un ENTIER nombre maximal d'iterations de dichotomie. Valeur par defaut 50. 'DELPE' un REEL, l'intervalle initial des iterations de dichotomie. La valeur par defaut est 1. 'MDELPE' un ENTIER nombre maximal de pas dans la recherche de l'intervalle de dichotomie. Valeur par defaut 20. ( evite de cycler lorsque l'on est tres loin de la solution) 'NFI' un ENTIER nombre de cycles de chimie. Valeur par defaut 4. Un cycle correspond a la sequence: * calcul de la force ionique * modification des logk |--- * boucle mineraux a |* resolution ( iterative ) precipiter | |* verification des mineraux | precipites |--- 'TEMPE' objet de type CHPOIN contenant la temperature. 'CLIM' valeur de l'activite imposee des especes de type 3. Objet de type CHPOIN ayant une composante pour chaque espece dont l'activite est imposee. 'SORTIE' 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 a preciser les elements que l'on veut voir figurer dans la TABLE TAB1. 'IMPRIM' un objet de type LISTENTI . Dans le cas ou l'on demande un niveau de message superieur a 0 ( OPTION IMPI 1 ), ceci permet de limiter les impressions aux seuls noeuds du maillage dont le numero 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 imposees chargement de CHPOIN au centre des faces 'FLUX_IMPOSE' Valeurs des flux de concentrations imposees chargement de CHPOIN au centre des faces 'SOURCE' Valeurs de la source Chargement relatifs aux CHPOIN centre ( defaut 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 sauvegardes (LISTREEL - defaut on sauve tous les pas de temps) Tous les temps de cette liste doivent etre 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-methode) (FLOTTANT - defaut 1.) Valeurs de theta generalement utilisees : Schema de Euler explicite : 0. Schema de Crank-Nicholson : 1/2 Schema de Euler implicite : 1. ---------------------------------------------- Conditions initales et Resultats ------------------------------------------- A chacun des indices suivants on trouve une table indicee par des entiers de 0 a N. L'indice i contient la valeur du champ au temps TEMPS.i. L'indice 0 contient les valeurs initiales,qui doivent etre donnees lors d'un premier calcul. Les noms des composantes des CHPOIN sont ceux utilises dans CHI2 (et attribues par le code). On peut obtenir les numeros des identifiants correspondants a l'aide des procedures NOCOMCHI et NOESPCHI. Indices indispensables en entree: 'TOT' 'TAQU' 'LOGC' Indices toujours presents en sortie: 'TOT' 'TAQU' 'LOGC' 'TEMPS' 'AQUE' 'FIXE' 'FLUX' Indices optionnels en entree: 'TEMPS' 'FION' 'NTY4' Indices optionnels en sortie: <'PREC'> <'FION'> <'TYP6'> <'TYP3'> <'NTY4'> <'TYP5'> <'SURF'> <'SOLU'> <'POLE'> <'LOGK'> 'TEMPS' TABLE contenant les temps sauvegardes (FLOTTANT) La valeur par defaut 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 elements. CHPOIN ayant une composante par composant chimique. (La valeur initiale est indispensable). 'AQUE' TABLE contenant les concentrations des aqueux aux centres des elements .CHPOIN ayant une composante par composant chimique ( cette table est generee par la procedure) 'FIXE' TABLE contenant les concentrations des fixees aux centres des elements .CHPOIN ayant une composante par composant chimique ( cette table est generee par la procedure) 'TAQU' Table contenant les concentrations des aqueux au centre des faces. CHPOIN possedant une composante par composant chimique.(La valeur initiale est indispensable) 'LOGC' Table contenant pour chaque composant chimique le log de la concentration de l'espece simple associee aux centres des elements.CHPOIN possedant une composante par composant chimique. (La valeur initiale est indispensable) 'FLUX' Table contenant les flux totaux (convectif+ diffusif) au centre des faces. CHPOIN possedant une composante par composant chimique.( cette table est generee par la procedure) 'PREC' TABLE contenant la concentration des especes precipitees. CHPOIN centre ayant une composante par espece chimique precipite. 'FION' TABLE contenant la force ionique. CHPOIN centre a une composante ( nom de la composante SCAL) 'TYP6' TABLE contenant la concentration des especes de type 6. CHPOIN centre ayant une composante par espece chimique de type 6.( cette table est generee par la procedure si TYP6 figure dans la liste TAB1.SORTIE) 'TYP3' TABLE contenant la concentration des especes de type 3. ( cette table est generee par la procedure si TYP3 figure dans la liste TAB1.SORTIE) 'NTY4' TABLE contenant des CHPOIN centre ayant une composante pour chaque espece de precipite potentiel. En chaque point centre d'element on indiquera si le mineral est precipite ( =1) ou non( =0). Sert a initialiser le choix des mineraux. 'TYP5' TABLE contenant la concentration des especes de type 5. CHPOIN centre ayant une composante par espece chimique de type 5.( cette table est generee par la procedure si TYP5 figure dans la liste TAB1.SORTIE) 'SURF' TABLE contenant la concentration des especes de surface. CHPOIN centre ayant une composante par espece chimique de surface.( cette table est generee par la procedure si SURF figure dans la liste TAB1.SORTIE) 'SOLU' TABLE contenant la concentration des especes en solution. CHPOIN centre ayant une composante par espece chimique en solution.( cette table est generee par la procedure 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 generee par la procedure 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 generee par la procedure si LOGK figure dans la liste TAB1.SORTIE)
© Cast3M 2003 - Tous droits réservés.
Mentions légales