Afficher cette notice en
$$$$ TRANSGEN NOTICE CHAT 11/09/12 21:18:24 7124 DATE 11/09/12 Operateur TRANSGEN Voir aussi : DARCYTRA, TRANGEOL ------------------ TRANSGEN TABLE ; ISSUE de la procedure DARCYTRA ! La syntaxe est conservee a l'exception de quelques points : 1 - quelques nouvelles fonctionnalites supplementaires - numeriques : choix VF EFMH, decentrement, mass lumping, solveur KRES accessible - physique : dispersivite calculee, nouvelles conditions aux limites (mixtes, flux total) 2 - syntaxe modifiee pour porosite et caracteristiques, et FLUXDIFF FLUXCONV au lieu de 'FLUX', CONVECTION est maintenant une vitesse et non un flux convectif, VITELEM la vitesse au centre est rajoutee pour les calculs de dispersivite. 3 - plus d'histoire de composantes 'H' 'TH' pour la concentration et sa trace car pas lieu d'etre en VF et relativement incompatible avec une gestion multiespece. La composante de concentration est libre (ex I129) et les conditions aux limites doivent avoir le meme nom de composante, ainsi que pour toutes les variables de meme dimension que la concentration (Concentration de saturation etc ...). Voire notice detaillee en dessous pour les noms de composantes. En gros, les jeux de donnees Darcytra tournent si les 'TH' sont transformes en 'H' plus modifs -1- et -2- (au plus quelques lignes dans les jeux de donnees), voir les jeux transport*.dgibi etc .. Fonction -------- Resolution de l'equation de transport de Radio nucleides en milieu poreux par une methode d'elements finis mixtes hybrides ou VF. Les inconnues du probleme sont - en EFMH, la concentration, la trace de concentration et le debit diffusif. - en VF, la concentration Gere pas de temps, retard, diffusion, dispersion, convection, source, preicipitation dissolution, decroissance, conditions aux limites (Dirichlet, Neumann, Mixtes, flux total) numerique : solveurs directs et iteratifs, decentrement, VF et EFMH, implicte explicite krank-Nicholson, mass lump Remarque -------- TRANSGEN remplace DARCYTRA pour le transport, DARCYTRA reste pour la resolut de l'equation de DARCY. Les personnes qui tiennent a utiliser des VF pour resoudre DARCY peuvent utiliser TRANSGEN en mettant une porosite nulle (annule le terme en temps) et une convection nulle. ---------------------------------- Operandes (a mettre dans TABLE) : ---------------------------------- ___________________________________________________________________ | | | Indice Contenu | | | ------------------------------------------------------------------- | | |------------------------------------------------ | |Donnees physiques, geometriques et materielles : | |------------------------------------------------ | | | |'MODELE' Objet modele (MMODEL cree par MODE) DARCY | | | |'CARACTERISTIQUES' Donnees physiques et materielles : | | diffusivite effective - CHAMPOINT de COMPOSANTES | | K11 K21 K22 K31 K32 K33 au centre des elements | | | |'POROSITE' Valeur de la porosite (Type Champoint, Comp | | 'SCAL', ou FLOTTANT) - Defaut 1. | | | |'DECROISSANCE' Valeur du terme de decroissance (Type FLOTTANT) | | Tel que dC/dt = - Lambda * C - Defaut 0. | | | |'COEF_RETARD' Coefficient de retard lineaire dans le cas simple, | | ou Pente a l'origine de la fonction F(C) dans le | | cas d'isotherme non lineaire de Langmuir | | ou Coefficient K de l'isotherme de Freundlich | | (Type CHPO Centre 'SCAL', ou FLOTTANT) | | | |'LANGMUIR' Quantite maximale 'Fsat' adsorbee sur le solide | | rapportee a l'unite de volume du fluide et exprimee| | dans la meme unite que le solute. | | (Type CHPO Centre 'SCAL', ou FLOTTANT). | | F = (R-1) C / [1 + ((R-1) C / Fsat)] | | Si cet indice et le suivant sont absents, | | l'equilibre d'adsorption est lineaire. Cet indice a| | priorite sur l'indice FREUNDLICH. | | | |'FREUNDLICH' Exposant de la loi de Freundlich F = K (C ^ 1/n) | | (Type FLOTTANT). | | Dans ce cas (et si l'indice LANGMUIR n'existe pas),| | l'indice 'COEF_RETARD' contient le coefficient | | K ramene a une unite de volume de fluide. | | - Non disponible pour l'instant - | | | |'LIMITE_SOLUBILITE' Limite de solubilite (Type chpoin), composante | identique a la concentration | | si absente pas de precipitation dissolution | | | |'COEF_DISSOLUTION' Coef. de dissolution (Type CHPO Centre, Comp | | 'SCAL'). Tel que dC/dt = Coef * (Csat - C) | | Si absent pas de dissolution precipitation | | | |'CONVECTION' vitesse au face. C'est le debit integre aux faces | | multiplie par la normale sortante de l'element | | et divise par la longueur de la face. | | Il s'agit de la projection du vecteur vitesse sur | | direction normale a la face. (Type CHPO Face, comp.| | VX VY VZ). L'interet est que cette vitesse est | | desormais intependante de l'orientation des normale| | ce qui est utile dans certains cas. | | |'VITELEM' Vitesse au centre des elements (Type CHPO centre, | | comp. VX VY VZ). Utilise uniquement si DECENTREMENT | | ou si dispersion | | | |'ALPHAL' coefficient de dispersivite longitudinale (CHPO de | | composante SCAL) - 0 si absent | | | |'ALPHAT' coefficient de dispersivite transverse (CHPO de | | composante SCAL) - 0 si absent | | Rque : si ALPHAL ou ALPHAT est present les deux | | doivent etre renseignes. | |---------------------- | |Conditions initiales : | |---------------------- | | | |'TEMPS' TABLE contenant a l'indice 0 la valeur du temps | | initial (FLOTTANT) | | | | | |'CONCENTRATION' TABLE contenant a l'indice 0 la concentration | | (quantite d'element par unite de volume d'eau) | | (Type CHPO Centre, Comp libre 4 lettres au plus) | | | | | |'PRECIPITE' TABLE contenant a l'indice 0 la quantite initiale | | de precipite par unite de volume de milieu solide | | (Type CHPO Centre, Comp identique a concentration) | | | |'FLUXDIFF' TABLE contenant a l'indice 0 les debits diffusifs | | integres | | par face (Type CHPO Face, nom de la concentration) | | | |'FLUXCONV' TABLE contenant les debits convectifs integres | | par face (Type CHPO Face, nom de la concentration) | | a l'indice 0. | | | |'DISSOLUTION' TABLE contenant a l'indice 0 la quantite initiale | | pour estimer la dissolution au premier pas de temps| | (Type CHPO, Comp de la concentration). | | | | | |-------------------------------------- | |Conditions aux limites / chargements : | |-------------------------------------- | | | | | |'TRACE_IMPOSE' Valeurs des traces imposees (charge ou concentra- | | -tion) - nom de la concentration | | | |'FLUX_IMPOSE' Valeurs des flux imposes integres par face | | (Type CHARGEMENT de CHPO Face) - nom concentration | | | |'FLUXTOT_IMP' Valeurs des flux totaux imposes integres par face | | (Type CHARGEMENT de CHPO Face, comp. nom de la | | concentration ) | | | |'MIXTES' Table : - indice C contient les valeurs des flux | | mixtes imposes integres par face | | (Type CHARGEMENT de CHPO Face, | | comp. idem concentration defaut 0.)| | - indices A et B sont des reels | | | | la condition mixte s'ecrit | | C = A * flux diffusif + B * Concentration | | | | | |'SOURCE' Valeurs du terme source par maille et par unite de | | temps (ex : puits, filiation) | | Les valeurs a l'indice i sont les valeurs entre | | les temps i-1 et i. | | (CHARGEMENT de CHPO Centre, comp de conc. ini) | | | | | |'DISSOLUTION_IMPOSEE' Valeurs des dissolutions imposees par unite| | de temps et par maille. (Type CHARGEMENT de CHPO, | | Comp conc). Les valeurs a l'indice i sont les | | valeurs moyennes de dissolution par unite de temps | | entre les temps i-1 et i. | | Priorite de la dissolution imposee sur les | | cinetiques. | | | |-------------------- | |Donnees numeriques : | |-------------------- | | | | | |'TEMPS_CALCULES' Valeur des temps calcules (LISTREEL) | | Contient obligatoirement le temps final. | | | |'TEMPS_SAUVES' Valeur des temps sauvegardes (LISTREEL - defaut : | | on sauve tous les pas de temps) | | | | 'LUMP' FAUX SI pas de mass lumping, VRAI sinon. | | VRAI seulement sur des maillages de rectangles et | | parallelepipedes rectangles et tenseur de dissusion| | orthotrope. Permet de rendre les schemas monotone | | pour la diffusion-instationnaire | | | | 'DECENTR' VRAI si diffusion numerique pour Peclet = 2, permet| | de stabiliser (en explicite) voire rendre monotone | | le schema de convection. | | FAUX si schema sans convection, ou en implicite et | | absence d'oscillations - plus precis | | | | 'TYPDISCRETISATION' 'VF' si VF et 'EFMH' si EFMH | | | | 'METHINV' TABLE DE PARAMETRE du solveur KRES, cf KRES | | peut etre remplie partiellement | | deux indices importants : | | 'TYPINV' OBLIGATOIRE 1 pour direct 3 pour BiCGSTAB | | 'PRECOND' obligatoire 3 pour ILU0, 5 pour ILUT | | conseil TYPINV = 1 en 2D ou petits calculs, 3 sinon| | conseil PRECOND = 3 sauf si problemes mettre 5 | | DERNIER CONSEIL : si message du type convergence | | breakdown, Pivot nul ... mettre la tolerence | | BCGSTOL a la precision machine 1.D-300, cf notice | | de KRES | | | |'THETA_DIFF' Coefficient de relaxation compris entre 0. et 1. | | (theta-methode diffusion) ('FLOTTANT' - defaut 1.) | | | |'THETA_CONVECTION' Idem pour la convection | | ('FLOTTANT', Defaut = THETA_DIFF) | |'THETA_DEC' Idem mais pour la decroissance | | ('FLOTTANT' - defaut 1/2) | | | |'THETA_DISS' Idem mais pour la dissolution | | ('FLOTTANT' - defaut 1.) | | | |'EPSI_RET' Precision relative d'arret pour la resolution | | iterative (Picard) de l'adsorption non lineaire | | (Type FLOTTANT) - Defaut 1.D-4 | | | |'EPSI_COR' Petit saut de concentration pour calculer le coef. | | de retard par la methode de la corde lorsque le | | retard est non-lineaire. | | (Type FLOTTANT) - Defaut 1.D-4 | | | |'ITMAX_RET' Nombre maxi d'iterations correspondant au retard | | non lineaire avant d'abandonner. | | (Type 'ENTIER') - Defaut 20 | |_________________________________________________________________| --------------------------------- Resultats (stockes dans TABLE) : --------------------------------- ___________________________________________________________________ | | | Indice Contenu | | | ------------------------------------------------------------------- | | | | |'TEMPS' TABLE contenant les temps sauvegardes (FLOTTANT) | | | |'CONCENTRATION' TABLE contenant les concentrations | | (Type CHPO Centre, Comp libre, celle de C ini) | | | |'FLUXDIFF' TABLE contenant les debits diffusifs integres | | par face (Type CHPO Face, nom de la concentration) | | | |'FLUXCONV' TABLE contenant les debits convectifs integres | | par face (Type CHPO Face, nom de la concentration) | | | |'PRECIPITE' TABLE contenant la quantite de precipite par maille| | (Type CHPO Centre, Comp de la concentration) | | | | | |'RETARD' Si cet indice a ete prealablement defini comme une | | TABLE, alors il contient les valeurs du coefficient| | de retard (Type 'CHPO' centre, Comp 'SCAL'). Sinon,| | les valeurs du coefficient de retard ne sont pas | | sauvegardees. | |_________________________________________________________________| ___________________________________________________________________ | | | Les tables resultats sont indicees par des entiers variant de 0 | | a N . | | A l'indice 0 on stocke les valeurs initiales, aux indices | | suivants les champs correspondant au temps de sortie TEMPS.I . | | Les champs servant en cas de reprise sont ceux correpondant au | | dernier indice. | |_________________________________________________________________|
© Cast3M 2003 - Tous droits réservés.
Mentions légales