Afficher cette notice en
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