Afficher cette notice en
1 : $$$$ TRANSGEN NOTICE CHAT 11/09/12 21:18:24 7124 2 : DATE 11/09/12 3 : 4 : 5 : Operateur TRANSGEN Voir aussi : DARCYTRA, TRANGEOL 6 : ------------------ 7 : TRANSGEN TABLE ; 8 : 9 : 10 : 11 : 12 : ISSUE de la procedure DARCYTRA ! 13 : La syntaxe est conservee a l'exception de quelques points : 14 : 15 : 1 - quelques nouvelles fonctionnalites supplementaires 16 : - numeriques : choix VF EFMH, decentrement, mass lumping, 17 : solveur KRES accessible 18 : - physique : dispersivite calculee, nouvelles conditions 19 : aux limites (mixtes, flux total) 20 : 21 : 2 - syntaxe modifiee pour porosite et caracteristiques, et FLUXDIFF 22 : FLUXCONV au lieu de 'FLUX', CONVECTION est maintenant une vitesse 23 : et non un flux convectif, VITELEM la vitesse au centre est rajoutee 24 : pour les calculs de dispersivite. 25 : 26 : 3 - plus d'histoire de composantes 'H' 'TH' pour la concentration 27 : et sa trace car pas lieu d'etre en VF et relativement incompatible 28 : avec une gestion multiespece. La composante de concentration est 29 : libre (ex I129) et les conditions aux limites doivent avoir 30 : le meme nom de composante, ainsi que pour toutes les variables de meme 31 : dimension que la concentration (Concentration de saturation etc ...). 32 : Voire notice detaillee en dessous pour les noms de composantes. 33 : En gros, les jeux de donnees Darcytra tournent si les 'TH' sont 34 : transformes en 'H' plus modifs -1- et -2- (au plus quelques lignes 35 : dans les jeux de donnees), voir les jeux transport*.dgibi etc .. 36 : 37 : Fonction 38 : -------- 39 : 40 : Resolution de l'equation de transport de Radio nucleides en milieu 41 : poreux par une methode d'elements finis mixtes hybrides ou VF. 42 : Les inconnues du probleme sont 43 : - en EFMH, la concentration, la trace de 44 : concentration et le debit diffusif. 45 : - en VF, la concentration 46 : Gere pas de temps, retard, diffusion, dispersion, convection, 47 : source, preicipitation dissolution, decroissance, conditions aux 48 : limites (Dirichlet, Neumann, Mixtes, flux total) 49 : numerique : solveurs directs et iteratifs, decentrement, VF et 50 : EFMH, implicte explicite krank-Nicholson, mass lump 51 : 52 : 53 : Remarque 54 : -------- 55 : 56 : TRANSGEN remplace DARCYTRA pour le transport, DARCYTRA reste pour la resolut 57 : de l'equation de DARCY. Les personnes qui tiennent a utiliser des VF 58 : pour resoudre DARCY peuvent utiliser TRANSGEN en mettant une porosite 59 : nulle (annule le terme en temps) et une convection nulle. 60 : 61 : 62 : 63 : ---------------------------------- 64 : Operandes (a mettre dans TABLE) : 65 : ---------------------------------- 66 : 67 : ___________________________________________________________________ 68 : | | 69 : | Indice Contenu | 70 : | | 71 : ------------------------------------------------------------------- 72 : | | 73 : |------------------------------------------------ | 74 : |Donnees physiques, geometriques et materielles : | 75 : |------------------------------------------------ | 76 : | | 77 : |'MODELE' Objet modele (MMODEL cree par MODE) DARCY | 78 : | | 79 : |'CARACTERISTIQUES' Donnees physiques et materielles : | 80 : | diffusivite effective - CHAMPOINT de COMPOSANTES | 81 : | K11 K21 K22 K31 K32 K33 au centre des elements | 82 : | | 83 : |'POROSITE' Valeur de la porosite (Type Champoint, Comp | 84 : | 'SCAL', ou FLOTTANT) - Defaut 1. | 85 : | | 86 : |'DECROISSANCE' Valeur du terme de decroissance (Type FLOTTANT) | 87 : | Tel que dC/dt = - Lambda * C - Defaut 0. | 88 : | | 89 : |'COEF_RETARD' Coefficient de retard lineaire dans le cas simple, | 90 : | ou Pente a l'origine de la fonction F(C) dans le | 91 : | cas d'isotherme non lineaire de Langmuir | 92 : | ou Coefficient K de l'isotherme de Freundlich | 93 : | (Type CHPO Centre 'SCAL', ou FLOTTANT) | 94 : | | 95 : |'LANGMUIR' Quantite maximale 'Fsat' adsorbee sur le solide | 96 : | rapportee a l'unite de volume du fluide et exprimee| 97 : | dans la meme unite que le solute. | 98 : | (Type CHPO Centre 'SCAL', ou FLOTTANT). | 99 : | F = (R-1) C / [1 + ((R-1) C / Fsat)] | 100 : | Si cet indice et le suivant sont absents, | 101 : | l'equilibre d'adsorption est lineaire. Cet indice a| 102 : | priorite sur l'indice FREUNDLICH. | 103 : | | 104 : |'FREUNDLICH' Exposant de la loi de Freundlich F = K (C ^ 1/n) | 105 : | (Type FLOTTANT). | 106 : | Dans ce cas (et si l'indice LANGMUIR n'existe pas),| 107 : | l'indice 'COEF_RETARD' contient le coefficient | 108 : | K ramene a une unite de volume de fluide. | 109 : | - Non disponible pour l'instant - | 110 : | | 111 : |'LIMITE_SOLUBILITE' Limite de solubilite (Type chpoin), composante 112 : | identique a la concentration | 113 : | si absente pas de precipitation dissolution | 114 : | | 115 : |'COEF_DISSOLUTION' Coef. de dissolution (Type CHPO Centre, Comp | 116 : | 'SCAL'). Tel que dC/dt = Coef * (Csat - C) | 117 : | Si absent pas de dissolution precipitation | 118 : | | 119 : |'CONVECTION' vitesse au face. C'est le debit integre aux faces | 120 : | multiplie par la normale sortante de l'element | 121 : | et divise par la longueur de la face. | 122 : | Il s'agit de la projection du vecteur vitesse sur | 123 : | direction normale a la face. (Type CHPO Face, comp.| 124 : | VX VY VZ). L'interet est que cette vitesse est | 125 : | desormais intependante de l'orientation des normale| 126 : | ce qui est utile dans certains cas. 127 : | | 128 : |'VITELEM' Vitesse au centre des elements (Type CHPO centre, | 129 : | comp. VX VY VZ). Utilise uniquement si DECENTREMENT | 130 : | ou si dispersion | 131 : | | 132 : |'ALPHAL' coefficient de dispersivite longitudinale (CHPO de | 133 : | composante SCAL) - 0 si absent | 134 : | | 135 : |'ALPHAT' coefficient de dispersivite transverse (CHPO de | 136 : | composante SCAL) - 0 si absent | 137 : | Rque : si ALPHAL ou ALPHAT est present les deux | 138 : | doivent etre renseignes. | 139 : |---------------------- | 140 : |Conditions initiales : | 141 : |---------------------- | 142 : | | 143 : |'TEMPS' TABLE contenant a l'indice 0 la valeur du temps | 144 : | initial (FLOTTANT) | 145 : | | 146 : | | 147 : |'CONCENTRATION' TABLE contenant a l'indice 0 la concentration | 148 : | (quantite d'element par unite de volume d'eau) | 149 : | (Type CHPO Centre, Comp libre 4 lettres au plus) | 150 : | | 151 : | | 152 : |'PRECIPITE' TABLE contenant a l'indice 0 la quantite initiale | 153 : | de precipite par unite de volume de milieu solide | 154 : | (Type CHPO Centre, Comp identique a concentration) | 155 : | | 156 : |'FLUXDIFF' TABLE contenant a l'indice 0 les debits diffusifs | 157 : | integres | 158 : | par face (Type CHPO Face, nom de la concentration) | 159 : | | 160 : |'FLUXCONV' TABLE contenant les debits convectifs integres | 161 : | par face (Type CHPO Face, nom de la concentration) | 162 : | a l'indice 0. | 163 : | | 164 : |'DISSOLUTION' TABLE contenant a l'indice 0 la quantite initiale | 165 : | pour estimer la dissolution au premier pas de temps| 166 : | (Type CHPO, Comp de la concentration). | 167 : | | 168 : | | 169 : |-------------------------------------- | 170 : |Conditions aux limites / chargements : | 171 : |-------------------------------------- | 172 : | | 173 : | | 174 : |'TRACE_IMPOSE' Valeurs des traces imposees (charge ou concentra- | 175 : | -tion) - nom de la concentration | 176 : | | 177 : |'FLUX_IMPOSE' Valeurs des flux imposes integres par face | 178 : | (Type CHARGEMENT de CHPO Face) - nom concentration | 179 : | | 180 : |'FLUXTOT_IMP' Valeurs des flux totaux imposes integres par face | 181 : | (Type CHARGEMENT de CHPO Face, comp. nom de la | 182 : | concentration ) | 183 : | | 184 : |'MIXTES' Table : - indice C contient les valeurs des flux | 185 : | mixtes imposes integres par face | 186 : | (Type CHARGEMENT de CHPO Face, | 187 : | comp. idem concentration defaut 0.)| 188 : | - indices A et B sont des reels | 189 : | | 190 : | la condition mixte s'ecrit | 191 : | C = A * flux diffusif + B * Concentration | 192 : | | 193 : | | 194 : |'SOURCE' Valeurs du terme source par maille et par unite de | 195 : | temps (ex : puits, filiation) | 196 : | Les valeurs a l'indice i sont les valeurs entre | 197 : | les temps i-1 et i. | 198 : | (CHARGEMENT de CHPO Centre, comp de conc. ini) | 199 : | | 200 : | | 201 : |'DISSOLUTION_IMPOSEE' Valeurs des dissolutions imposees par unite| 202 : | de temps et par maille. (Type CHARGEMENT de CHPO, | 203 : | Comp conc). Les valeurs a l'indice i sont les | 204 : | valeurs moyennes de dissolution par unite de temps | 205 : | entre les temps i-1 et i. | 206 : | Priorite de la dissolution imposee sur les | 207 : | cinetiques. | 208 : | | 209 : |-------------------- | 210 : |Donnees numeriques : | 211 : |-------------------- | 212 : | | 213 : | | 214 : |'TEMPS_CALCULES' Valeur des temps calcules (LISTREEL) | 215 : | Contient obligatoirement le temps final. | 216 : | | 217 : |'TEMPS_SAUVES' Valeur des temps sauvegardes (LISTREEL - defaut : | 218 : | on sauve tous les pas de temps) | 219 : | | 220 : | 'LUMP' FAUX SI pas de mass lumping, VRAI sinon. | 221 : | VRAI seulement sur des maillages de rectangles et | 222 : | parallelepipedes rectangles et tenseur de dissusion| 223 : | orthotrope. Permet de rendre les schemas monotone | 224 : | pour la diffusion-instationnaire | 225 : | | 226 : | 'DECENTR' VRAI si diffusion numerique pour Peclet = 2, permet| 227 : | de stabiliser (en explicite) voire rendre monotone | 228 : | le schema de convection. | 229 : | FAUX si schema sans convection, ou en implicite et | 230 : | absence d'oscillations - plus precis | 231 : | | 232 : | 'TYPDISCRETISATION' 'VF' si VF et 'EFMH' si EFMH | 233 : | | 234 : | 'METHINV' TABLE DE PARAMETRE du solveur KRES, cf KRES | 235 : | peut etre remplie partiellement | 236 : | deux indices importants : | 237 : | 'TYPINV' OBLIGATOIRE 1 pour direct 3 pour BiCGSTAB | 238 : | 'PRECOND' obligatoire 3 pour ILU0, 5 pour ILUT | 239 : | conseil TYPINV = 1 en 2D ou petits calculs, 3 sinon| 240 : | conseil PRECOND = 3 sauf si problemes mettre 5 | 241 : | DERNIER CONSEIL : si message du type convergence | 242 : | breakdown, Pivot nul ... mettre la tolerence | 243 : | BCGSTOL a la precision machine 1.D-300, cf notice | 244 : | de KRES | 245 : | | 246 : |'THETA_DIFF' Coefficient de relaxation compris entre 0. et 1. | 247 : | (theta-methode diffusion) ('FLOTTANT' - defaut 1.) | 248 : | | 249 : |'THETA_CONVECTION' Idem pour la convection | 250 : | ('FLOTTANT', Defaut = THETA_DIFF) | 251 : |'THETA_DEC' Idem mais pour la decroissance | 252 : | ('FLOTTANT' - defaut 1/2) | 253 : | | 254 : |'THETA_DISS' Idem mais pour la dissolution | 255 : | ('FLOTTANT' - defaut 1.) | 256 : | | 257 : |'EPSI_RET' Precision relative d'arret pour la resolution | 258 : | iterative (Picard) de l'adsorption non lineaire | 259 : | (Type FLOTTANT) - Defaut 1.D-4 | 260 : | | 261 : |'EPSI_COR' Petit saut de concentration pour calculer le coef. | 262 : | de retard par la methode de la corde lorsque le | 263 : | retard est non-lineaire. | 264 : | (Type FLOTTANT) - Defaut 1.D-4 | 265 : | | 266 : |'ITMAX_RET' Nombre maxi d'iterations correspondant au retard | 267 : | non lineaire avant d'abandonner. | 268 : | (Type 'ENTIER') - Defaut 20 | 269 : |_________________________________________________________________| 270 : 271 : 272 : 273 : --------------------------------- 274 : Resultats (stockes dans TABLE) : 275 : --------------------------------- 276 : 277 : ___________________________________________________________________ 278 : | | 279 : | Indice Contenu | 280 : | | 281 : ------------------------------------------------------------------- 282 : | | 283 : | | 284 : |'TEMPS' TABLE contenant les temps sauvegardes (FLOTTANT) | 285 : | | 286 : |'CONCENTRATION' TABLE contenant les concentrations | 287 : | (Type CHPO Centre, Comp libre, celle de C ini) | 288 : | | 289 : |'FLUXDIFF' TABLE contenant les debits diffusifs integres | 290 : | par face (Type CHPO Face, nom de la concentration) | 291 : | | 292 : |'FLUXCONV' TABLE contenant les debits convectifs integres | 293 : | par face (Type CHPO Face, nom de la concentration) | 294 : | | 295 : |'PRECIPITE' TABLE contenant la quantite de precipite par maille| 296 : | (Type CHPO Centre, Comp de la concentration) | 297 : | | 298 : | | 299 : |'RETARD' Si cet indice a ete prealablement defini comme une | 300 : | TABLE, alors il contient les valeurs du coefficient| 301 : | de retard (Type 'CHPO' centre, Comp 'SCAL'). Sinon,| 302 : | les valeurs du coefficient de retard ne sont pas | 303 : | sauvegardees. | 304 : |_________________________________________________________________| 305 : 306 : 307 : ___________________________________________________________________ 308 : | | 309 : | Les tables resultats sont indicees par des entiers variant de 0 | 310 : | a N . | 311 : | A l'indice 0 on stocke les valeurs initiales, aux indices | 312 : | suivants les champs correspondant au temps de sortie TEMPS.I . | 313 : | Les champs servant en cas de reprise sont ceux correpondant au | 314 : | dernier indice. | 315 : |_________________________________________________________________| 316 : 317 :
© Cast3M 2003 - Tous droits réservés.
Mentions légales