Afficher cette notice en
$$$$ TRANGEOL NOTICE CB215821 15/06/25 21:15:05 8554 DATE 15/06/25 Procedure TRANGEOL Voir aussi : ------------------ RES1 RES2 = TRANGEOL modarcy TRANS1 (TRANS2); APPELE par TRANSGEN - DECONSEILLE AUX UTILISATEURS NON DEVELOPPEURS resout : POROSITE * DC/DT = DIV (DIFFUSIVITE GRAD C - CONVECTION C) + SOURCE ou C est la CONCENTRATION, eventuellement a plusieurs composantes (multiespeces non couplees). La diffusivite, la convection ... doivent etre toutefois les memes pour toutes les especes pour avoir des matrices de discretisation communes. seul le terme source est donne par espece. Condition aux limites Neumann DIrichlet Mixtes flux total DIFFUSIVITE = DIFFUSIVITE entree + Dispersivite calculee Decentrement automatique (peut etre desactive en EFMH). Resolution sur 1 avancee en TEMPS. Discretisation VF ou EFMH Solveur KRES ------------------------------------------------------------------- |modarcy Objet modele (MMODEL cree par MODE) DARCY | ------------------------------------------------------------------- ------------------------------------------------------------------- | TRANS1 Table contenant les indices suivants : | |---------------------------------------------------------------- | |'DIFFUSIVITE' Donnees physiques et materielles : | | diffusivite effective - CHAMPOINT de COMPOSANTES | | K11 K21 K22 K31 K32 K33 au centre des elements | | | |'POROSITE' Valeur du coef devant D/DT (Type Champoint, Comp | | 'SCAL', ou FLOTTANT) - Defaut 1. | | | | DELTAT pas de temps | | | |'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. | | OPTIONNEL | | | |'VITELEM' Vitesse au centre des elements (Type CHPO centre, | | comp. VX VY VZ). Utilise uniquement si DECENTREMENT | | ou si dispersion. OPTIONNEL donc | | | |'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 : | |---------------------- | | | |'CONCENTRATION' concentration en debut de pas de temps | | (quantite d'element par unite de volume d'eau) | | (Type CHPO Centre, Comp libre 4 lettres au plus) | | la concentration peut avoir plusieurs composante | | la resolution etant alors faite pour chaque | | composante | | | |-------------------------------------- | |Conditions aux limites / chargements : | |-------------------------------------- | ------------------------------------------------------------------- | | | 'CLIMITES' table contenant les indices suivants : | | | |'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) | | | | | |-------------------- | |Donnees numeriques : | |-------------------- | | | | | | '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 - OBLIGATOIRE | | | | 'DECENTREMENT' 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 | | OBLIGATOIRE | | | | 'TYPDISCRETISATION' 'VF' si VF et 'EFMH' si EFMH | | | | | |'THETA_DIFFUSION' Coefficient de relaxation compris entre 0. et 1| | (theta-methode diffusion) ('FLOTTANT' - OBLIGATOIRE)| | | |'THETA_CONVECTION' Idem pour la convection | | ('FLOTTANT', OBLIGATOIRE) | ------------------------------------------------------------------- | Solveurs : - OPTIONNELS |----------- | | | 'METHINV' TABLE DE PARAMETRE du solveur KRES, cf KRES | | peut etre remplie partiellement | | deux indices importants : | | 'TYPINV' 1 pour direct 3 pour BiCGSTAB | | 'PRECOND' 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 | | METHINV est OPTIONNELLE | | | | 'SOLVEUR' type de solveur 1 pour direct | | 2 grad conjugue | | 3 BiCGSTAB | | 4 GMRES | | ecrase le choix dans METHINV si double entree | | OPTIONNEL EGALEMENT | | | | 'PRECONDITIONNEUR' type de precond. 0 rien | | 1 diag. | | 2 D-ILU | | 3 ILU0 | | 4 MILU | | 5 ILUT | | 6 ILUT2 | | OPTIONNEL - Ecrase choix dans METHINV si doub| | le entree. | | | |----------- | | recalculs - OPTIONNELS | |----------- | | 'ABANDON' Si VRAI alors les resultats sont nuls lorsque | | la concentration est inferieur a 'SEUILCALC' | | un calcul non couteux est effectue a chaque pas de| | temps. Le calcul reprend normalement si la conc. | | refranchit a nouveau le seuil (source accrue....)| | | | 'SEUILCALC' Si ABANDON est vrai il faut preciser le seuil | | d'arret des calculs | | | |----------- | | pilotage -OPTIONNELS | |----------- | | 'RECALCUL' si VRAI alors toutes les matrices seront recalculees | | | | 'INITIALI' Si vrai on considere que c'est le premier appel | | On a RECALCUL mis a VRAI et rien de TRANS2 ne sert | | ie tous les objets de preconditionnement sont recalc | | | | 'MODIFICATI' TABLE plus detaillee que 'RECALCUL' contient les | | indices : | | 'POROSITE' = VRAI si recalcul lie a variation de | | la porosite | | 'CONVECTI' = VRAI si recalcul lie a variation de | | la CONVECTIOn | | 'DELTAT' = VRAI si recalcul lie a variation de | | DELTAT | | 'DIFFUSIV' = VRAI si recalcul lie a variation de | | DIFFUSIVITE | | | | En cas d'absence de tout indice on reinitialise | | REINIT domine RECALCUL qui domine MODIFICATI | ------------------------------------------------------------------ ------------------------------------------------------------------- | TRANS2 - OPTIONNEL | ------------------------------------------------------------------ | Table optionnelle qui contient tous les preconditionnements | | comme les matrices de discretisation, le surfaces des elements| | etc... calules lors d'un precedent appel | | | | TRANS2 est donc en general le resultat RES2 d'un precedent | | calcul reinjecte en entree de TRANGEOL | | selon les options de pilotages precedentes, tout | | ou partie de TRANS2 sera utilisee. | | TRANS2 contient notamment les traces de concentration | | en EFMH | | | | Ne pas donner TRANS2 revient a faire un premier calcul | | C'est identique a REINIT = VRAI | | | | Ce preconditionneur n'a d'interet que pour des calculs a pas | | de temps constant (en EFMH) ou relativement lineaires (en EFMH | | et VF) | |_________________________________________________________________| SORTIES RES1 table identique a TRANS1 mais donc la concentration est remise a jour (nouvelle apres calcul). L'indice INITIALI est mis a FAUX (cree si absent auparavent). Les indices FLUXDIFF et FLUXCONV contenant les CHPO de flux convectifs et diffusifs integres aux faces sont ajoutes. RES2 table contenant les objets de preconditionnement pour une utilisation plus rapide lors d'appels suivants. UTILISATION premier appel : RES1 RES2 = TRANGEOL MODARCY TRANS1 ; apres : RES1 RES2 = TRANGEOL MODARCY RES1 RES2 ; point fixe : reflechir !!!
© Cast3M 2003 - Tous droits réservés.
Mentions légales