Télécharger trangeol.notice

Retour à la liste thématique

Afficher cette notice en

Numérotation des lignes :
   1 : $$$$ TRANGEOL NOTICE  CB215821  15/06/25    21:15:05     8554           
   2 :                                              DATE     15/06/25
   3 : 
   4 :   Procedure TRANGEOL                     Voir aussi :
   5 :       ------------------  
   6 :       RES1 RES2 = TRANGEOL modarcy TRANS1 (TRANS2);
   7 : 
   8 : 
   9 : 
  10 : 
  11 :    APPELE par TRANSGEN - DECONSEILLE AUX UTILISATEURS NON DEVELOPPEURS
  12 : 
  13 :    resout : POROSITE * DC/DT = DIV (DIFFUSIVITE GRAD C - CONVECTION C)
  14 :                                + SOURCE
  15 : 
  16 :              ou C est la CONCENTRATION, eventuellement a plusieurs
  17 :                   composantes (multiespeces non couplees).
  18 :                   La diffusivite, la convection ... doivent etre
  19 :                   toutefois les memes pour toutes les especes
  20 :                   pour avoir des matrices de discretisation communes.
  21 :                   seul le terme source est donne par espece.
  22 : 
  23 :             Condition aux limites Neumann DIrichlet Mixtes flux total
  24 : 
  25 :             DIFFUSIVITE = DIFFUSIVITE entree + Dispersivite calculee
  26 : 
  27 :             Decentrement automatique (peut etre desactive en EFMH).
  28 : 
  29 :             Resolution sur 1 avancee en TEMPS.
  30 : 
  31 :             Discretisation VF ou EFMH
  32 :  
  33 :             Solveur KRES
  34 : 
  35 :   -------------------------------------------------------------------
  36 :   |modarcy       Objet modele (MMODEL cree par MODE) DARCY          |
  37 :   -------------------------------------------------------------------
  38 : 
  39 :   -------------------------------------------------------------------
  40 :   | TRANS1       Table contenant les indices suivants :             |
  41 :   |---------------------------------------------------------------- |
  42 :   |'DIFFUSIVITE' Donnees physiques et materielles :                 |
  43 :   |              diffusivite effective - CHAMPOINT de COMPOSANTES   |
  44 :   |              K11 K21 K22 K31 K32 K33 au centre des elements     |
  45 :   |                                                                 |
  46 :   |'POROSITE'    Valeur du coef devant D/DT (Type Champoint, Comp   |
  47 :   |              'SCAL', ou FLOTTANT) - Defaut 1.                   |
  48 :   |                                                                 |
  49 :   | DELTAT       pas de temps                                       |
  50 :   |                                                                 |
  51 :   |'CONVECTION'  vitesse au face. C'est le debit integre aux faces  |
  52 :   |              multiplie par la normale sortante de l'element     |
  53 :   |              et divise par la longueur de la face.              |
  54 :   |              Il s'agit de la projection du vecteur vitesse sur  |
  55 :   |              direction normale a la face. (Type CHPO Face, comp.|
  56 :   |              VX VY VZ). L'interet est que cette vitesse est     |
  57 :   |              desormais intependante de l'orientation des normale|
  58 :   |              ce qui est utile dans certains cas.                |
  59 :   |              OPTIONNEL                                          |
  60 :   |                                                                 |
  61 :   |'VITELEM'    Vitesse au centre des elements (Type CHPO centre,   |
  62 :   |             comp. VX VY VZ). Utilise uniquement si DECENTREMENT |
  63 :   |             ou si dispersion. OPTIONNEL donc                    |
  64 :   |                                                                 |
  65 :   |'ALPHAL'      coefficient de dispersivite longitudinale (CHPO de |
  66 :   |              composante SCAL) - 0 si absent                     |
  67 :   |                                                                 |
  68 :   |'ALPHAT'      coefficient de dispersivite transverse (CHPO de    |
  69 :   |              composante SCAL) - 0 si absent                     |
  70 :   |              Rque : si ALPHAL ou ALPHAT est present les deux    |
  71 :   |              doivent etre renseignes.                           | 
  72 :   |                                                                 |
  73 :   |----------------------                                           |
  74 :   |Conditions initiales :                                           |
  75 :   |----------------------                                           |
  76 :   |                                                                 |
  77 :   |'CONCENTRATION' concentration en debut de pas de temps           |
  78 :   |              (quantite d'element par unite de volume d'eau)     |
  79 :   |              (Type CHPO Centre, Comp libre 4 lettres au plus)   |
  80 :   |               la concentration peut avoir plusieurs composante  |
  81 :   |               la resolution etant alors faite pour chaque       |
  82 :   |               composante                                        |
  83 :   |                                                                 |
  84 :   |--------------------------------------                           |
  85 :   |Conditions aux limites / chargements :                           |
  86 :   |--------------------------------------                           |
  87 :   -------------------------------------------------------------------
  88 :   |                                                                 |
  89 :   | 'CLIMITES'   table contenant les indices suivants :             |
  90 :   |                                                                 |
  91 :   |'TRACE_IMPOSE' Valeurs des traces imposees (charge ou concentra- |
  92 :   |              -tion) - nom de la concentration                   |
  93 :   |                                                                 |
  94 :   |'FLUX_IMPOSE' Valeurs des flux imposes integres par face         |
  95 :   |              (Type CHARGEMENT de CHPO Face) - nom concentration |
  96 :   |                                                                 |
  97 :   |'FLUXTOT_IMP' Valeurs des flux totaux imposes integres par face  |
  98 :   |              (Type CHARGEMENT de CHPO Face, comp. nom de la     |
  99 :   |               concentration )                                   |
 100 :   |                                                                 |
 101 :   |'MIXTES'      Table : - indice C contient les valeurs des flux   |
 102 :   |                               mixtes imposes integres par face  |
 103 :   |                             (Type CHARGEMENT de CHPO Face,      |
 104 :   |                              comp. idem concentration defaut 0.)|
 105 :   |                      - indices A et B sont des reels            |
 106 :   |                                                                 |
 107 :   |              la condition mixte s'ecrit                         |
 108 :   |                   C =  A * flux diffusif +  B * Concentration   |
 109 :   |                                                                 |
 110 :   ------------------------------------------------------------------
 111 :   |                                                                 |
 112 :   |'SOURCE'      Valeurs du terme source par maille et par unite de |
 113 :   |              temps (ex : puits, filiation)                      |
 114 :   |              Les valeurs a l'indice i sont les valeurs entre    |
 115 :   |              les temps i-1 et i.                                |
 116 :   |              (CHARGEMENT de CHPO Centre, comp de conc. ini)     |
 117 :   |                                                                 |
 118 :   |                                                                 |
 119 :   |--------------------                                             |
 120 :   |Donnees numeriques :                                             |
 121 :   |--------------------                                             |
 122 :   |                                                                 |
 123 :   |                                                                 |
 124 :   | 'LUMP'       FAUX SI pas de mass lumping, VRAI sinon.           |
 125 :   |              VRAI seulement sur des maillages de rectangles et  |
 126 :   |              parallelepipedes rectangles et tenseur de dissusion|
 127 :   |              orthotrope. Permet de rendre les schemas monotone  |
 128 :   |              pour la diffusion-instationnaire - OBLIGATOIRE     |
 129 :   |                                                                 |
 130 :   | 'DECENTREMENT' VRAI si diffusion numerique pour Peclet = 2,     |
 131 :   |                permet                                           |
 132 :   |              de stabiliser (en explicite) voire rendre monotone |
 133 :   |              le schema de convection.                           |
 134 :   |              FAUX si schema sans convection, ou en implicite et |
 135 :   |              absence d'oscillations - plus precis               |
 136 :   |              OBLIGATOIRE                                        |
 137 :   |                                                                 |
 138 :   | 'TYPDISCRETISATION'   'VF' si VF et 'EFMH' si EFMH              |
 139 :   |                                                                 |
 140 :   |                                                                 |
 141 :   |'THETA_DIFFUSION' Coefficient de relaxation compris entre 0. et 1|
 142 :   |              (theta-methode diffusion) ('FLOTTANT' - OBLIGATOIRE)|
 143 :   |                                                                 |
 144 :   |'THETA_CONVECTION' Idem pour la convection                       |
 145 :   |              ('FLOTTANT', OBLIGATOIRE)                          |
 146 :   -------------------------------------------------------------------
 147 :   | Solveurs : - OPTIONNELS
 148 :   |-----------
 149 :   |                                                                 |
 150 :   | 'METHINV'    TABLE DE PARAMETRE du solveur KRES, cf KRES        |
 151 :   |              peut etre remplie partiellement                    |
 152 :   |              deux indices importants :                          |
 153 :   |              'TYPINV'  1 pour direct 3 pour BiCGSTAB            |
 154 :   |              'PRECOND' 3 pour ILU0, 5 pour ILUT                 |
 155 :   |              conseil TYPINV = 1 en 2D ou petits calculs, 3 sinon|
 156 :   |              conseil PRECOND = 3 sauf si problemes mettre 5     |
 157 :   |              DERNIER CONSEIL : si message du type convergence   |
 158 :   |              breakdown, Pivot nul ... mettre la tolerence       |
 159 :   |              BCGSTOL a la precision machine 1.D-300, cf notice  |
 160 :   |              de KRES                                            |
 161 :   |              METHINV est OPTIONNELLE                            |
 162 :   |                                                                 |
 163 :   | 'SOLVEUR'    type de solveur 1 pour direct                      |
 164 :   |                              2 grad conjugue                    |
 165 :   |                              3 BiCGSTAB                         |
 166 :   |                              4 GMRES                            |
 167 :   |              ecrase le choix dans METHINV si double entree      |
 168 :   |              OPTIONNEL EGALEMENT                                |
 169 :   |                                                                 |
 170 :   | 'PRECONDITIONNEUR' type de precond. 0 rien                      |
 171 :   |                                     1 diag.                     |
 172 :   |                                     2 D-ILU                     |
 173 :   |                                     3 ILU0                      |
 174 :   |                                     4 MILU                      |
 175 :   |                                     5 ILUT                      |
 176 :   |                                     6 ILUT2                     |
 177 :   |                    OPTIONNEL - Ecrase choix dans METHINV si doub|
 178 :   |                                le entree.                       |
 179 :   |                                                                 |
 180 :   |-----------                                                      |
 181 :   | recalculs  - OPTIONNELS                                         |
 182 :   |-----------                                                      |
 183 :   | 'ABANDON'      Si VRAI alors les resultats sont nuls lorsque    |
 184 :   |                la concentration est inferieur a 'SEUILCALC'     |
 185 :   |                un calcul non couteux est effectue a chaque pas de|
 186 :   |                temps. Le calcul reprend normalement si la conc. |
 187 :   |                refranchit a nouveau le seuil (source accrue....)|
 188 :   |                                                                 |
 189 :   | 'SEUILCALC'    Si ABANDON est vrai il faut preciser le seuil    |
 190 :   |                d'arret des calculs                              |
 191 :   |                                                                 |
 192 :   |-----------                                                      |
 193 :   | pilotage   -OPTIONNELS                                          |
 194 :   |-----------                                                      |
 195 :   | 'RECALCUL' si VRAI alors toutes les matrices seront recalculees |
 196 :   |                                                                 |
 197 :   | 'INITIALI' Si vrai on considere que c'est le premier appel      |
 198 :   |            On a RECALCUL mis a VRAI et rien de TRANS2 ne sert   |
 199 :   |            ie tous les objets de preconditionnement sont recalc |
 200 :   |                                                                 |
 201 :   | 'MODIFICATI' TABLE plus detaillee que 'RECALCUL' contient les   |
 202 :   |               indices :                                         |
 203 :   |             'POROSITE' = VRAI si recalcul lie a variation de    |
 204 :   |                               la porosite                       |
 205 :   |             'CONVECTI' = VRAI si recalcul lie a variation de    |
 206 :   |                               la CONVECTIOn                     |
 207 :   |             'DELTAT'   = VRAI si recalcul lie a variation de    |
 208 :   |                               DELTAT                            |
 209 :   |             'DIFFUSIV' = VRAI si recalcul lie a variation de    |
 210 :   |                               DIFFUSIVITE                       |
 211 :   |                                                                 |
 212 :   | En cas d'absence de tout indice on reinitialise                 |
 213 :   | REINIT domine RECALCUL qui domine MODIFICATI                    |
 214 :   ------------------------------------------------------------------
 215 : 
 216 :   -------------------------------------------------------------------
 217 :   |                        TRANS2 - OPTIONNEL                       |
 218 :   ------------------------------------------------------------------
 219 :   |   Table optionnelle qui contient tous les preconditionnements   |
 220 :   |   comme les matrices de discretisation, le surfaces des elements|
 221 :   |   etc... calules lors d'un precedent appel                      |
 222 :   |                                                                 |
 223 :   |   TRANS2 est donc en general le resultat RES2 d'un precedent    |
 224 :   |          calcul reinjecte en entree de TRANGEOL                 |
 225 :   |          selon les options de pilotages precedentes, tout       |
 226 :   |          ou partie de TRANS2 sera utilisee.                     |
 227 :   |          TRANS2 contient notamment les traces de concentration  |
 228 :   |          en EFMH                                                |
 229 :   |                                                                 |
 230 :   |  Ne pas donner TRANS2 revient a faire un premier calcul         |
 231 :   |  C'est identique a REINIT = VRAI                                |
 232 :   |                                                                 |
 233 :   |  Ce preconditionneur n'a d'interet que pour des calculs a pas   |
 234 :   |  de temps constant (en EFMH) ou relativement lineaires (en EFMH |
 235 :   |  et VF)                                                         |
 236 :   |_________________________________________________________________|
 237 : 
 238 : 
 239 :                         SORTIES
 240 : 
 241 :   RES1 table identique a TRANS1 mais donc la concentration est remise
 242 :        a jour (nouvelle apres calcul).
 243 :        L'indice INITIALI est mis a FAUX (cree si absent auparavent).
 244 :        Les indices FLUXDIFF et FLUXCONV contenant les CHPO de flux
 245 :        convectifs et diffusifs integres aux faces sont ajoutes.
 246 : 
 247 : 
 248 :   RES2 table contenant les objets de preconditionnement pour une 
 249 :        utilisation plus rapide lors d'appels suivants.
 250 : 
 251 : 
 252 : 
 253 :                       UTILISATION
 254 : 
 255 :       premier appel :   RES1 RES2 = TRANGEOL MODARCY TRANS1 ;
 256 :       apres         :   RES1 RES2 = TRANGEOL MODARCY RES1 RES2 ;
 257 :       point fixe : reflechir !!!
 258 : 
 259 : 

© Cast3M 2003 - Tous droits réservés.
Mentions légales