Télécharger transgen.notice

Retour à la liste

Afficher cette notice en

Numérotation des lignes :
   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