Télécharger updaEFMH.notice

Retour à la liste thématique

Afficher cette notice en

Numérotation des lignes :
   1 : $$$$ UPDAEFMH NOTICE  CHAT      11/09/12    21:18:27     7124           
   2 :                                              DATE     11/09/12
   3 : 
   4 :   Operateur UPDAEFMH                   Voir aussi :
   5 :         ------------------  
   6 : 
   7 :   ATTENTION La vitesse est optionnelle, L'ordre est important
   8 :   et les types d'arguments qui se suivent aussi pour tester leur
   9 :  presence
  10 :  
  11 :   APPELE PAR TRANGEOL - PAS POUT UTILISATEUR
  12 :  
  13 :    |-----------------------------------------------------------------|
  14 :    | Phrase d'appel (en GIBIANE)                                     |
  15 :    |-----------------------------------------------------------------|
  16 :    |                                                                 |
  17 :    |SMTr MatrTr TbDarTra MassEFMH Difftot = UPDAEFMH MoDARCY Porosite|
  18 :    |                     MateDiff difftot ChPSour cini tcini deltat  |
  19 :    |                     (Qface) nomespec nbespece nbsource LMLump   |
  20 :    |                     DECENTR massEFMH | mattr tbdartra TABMODI;  |
  21 :    |                                      | mattm                    |
  22 :    |-----------------------------------------------------------------|
  23 :    | Generalites : UPDAEFMH construit la matrice de discretisation   |
  24 :    |               du probleme de transport convection-diffusion pour|
  25 :    |               le premier pas de tps d'un algorithme transitoire.|
  26 :    |               Le second membre et les Conditions limites de flux|
  27 :    |               sont pris en compte.                              |
  28 :    |            RESTE TCINI, DECENTR et TERME LIN                    |
  29 :    |-----------------------------------------------------------------|
  30 :    |                                                                 |
  31 :    |-----------------------------------------------------------------|
  32 :    |                         ENTREES                                 |
  33 :    |-----------------------------------------------------------------|
  34 :    | MoDARCY  : modele Darcy.                                        |
  35 :    |                                                                 |
  36 :    | Porosite : champ par elements de composante 'CK'                |
  37 :    |                                                                 |
  38 :    | MateDiff : Tenseur de diffusion  (type iso, ..) champ par       |
  39 :    |            points de composante 'K' en isotrope, 'K11', 'K21',  |
  40 :    |            'K22' en anisotrope 2d et  'K11', 'K21', 'K22', 'K31'|
  41 :    |            'K32', 'K33' en anisotrope 3d. Type 'CARACTERISTIQUE'|
  42 :    |                                                                 |
  43 :    | Diffdisp : Tenseur de dispersion  (type iso, ..) champ par      |
  44 :    |            points de composante 'K' en isotrope, 'K11', 'K21',  |
  45 :    |            'K22' en anisotrope 2d et  'K11', 'K21', 'K22', 'K31'|
  46 :    |            'K32', 'K33' en anisotrope 3d. Type 'CARACTERISTIQUE'|
  47 :    |                                                                 |
  48 :    | ChPSour  : Champ par points des sources volumiques par unite de |
  49 :    |            temps (support maillage centre). Composante ??????   |
  50 :    |                                                                 |
  51 :    | Cini     : Concentration initiale, CHPOINT centre.              |
  52 :    |            Composante 'H'.                                      |
  53 :    |                                                                 |
  54 :    | Tcini    : Trace de concentration aux faces (eventuellement a   |
  55 :    |            plusieurs composantes (especes)                      |
  56 :    |                                                                 |
  57 :    | Deltat   : Pas de temps                                         |
  58 :    |                                                                 |
  59 :    | Qface    : vitesse aux faces, CHPO face de composantes Vx, Vy   |
  60 :    |            en 2d et Vx, Vy, Vz en 3d. Il s'agit plus exatement  |
  61 :    |            de (V.n)n, c'est a dire de la composante normale de  |
  62 :    |            la vitesse aux faces. ???????? (je pressens que      |
  63 :    |            castem va sortir des flux, cad integres sur surfaces)|
  64 :    |                                                                 |
  65 :    | nomespec : liste des noms de composante des especes dans Cini   |
  66 :    |                                                                 |
  67 :    | nbespece : nombre de composante de Cini, soit nombre d'especes  |
  68 :    |                                                                 |
  69 :    | nbsource : nombre de composantes du terme source qd X especes   |
  70 :    |                                                                 |
  71 :    | LMLump   : Logique. Si vrai on effectue une condensation de     |
  72 :    |            masse de la matrice EFMH                             |
  73 :    |                                                                 |
  74 :    | DECENTR  : Logique. Vrai veut dire schemas decentres et faux    |
  75 :    |            veut dire schema convectif centre.                   |        
  76 :    |                                                                 |
  77 :    | MatTm    : matrice globale sur les traces. MATRIK en entree     |
  78 :    |            sort MATRIK si non modifiee, RIGIDITE sinon          |
  79 :    |            Soit on rentre cet argument soit le suivant Mattr    |
  80 :    |                                                                 |
  81 :    | MatTr    : idem mais rigidite en entree on ressort cette matrice|
  82 :    |            inchangee si les options MATMODI indiquent aucune    |
  83 :    |            modif. Optionnel. On rentre Mattm si absent.         |
  84 :    |                                                                 |
  85 :    | TbDarTra : table Darcy transitoire utilisee par MHYB, SMTP ...  |
  86 :    |                                                                 |
  87 :    | TABMODI  : table contenant des logiques indiquant la necessite  |
  88 :    |            ou non de reclalculer certains termes.               |
  89 :    |            'POROSITE' : VRAI si le coefficient devant D/DT      |
  90 :    |                         (porosite) est modifie depuis le dernier|
  91 :    |                         appel                                   |
  92 :    |            'DELTAT'   : VRAI si le pas de tps a change          |
  93 :    |            'CONVECTI' : VRAI si la vitesse a change             |
  94 :    |            'COEF_LIN' : VRAI si le coef en facteur de C a change|
  95 :    |            'DIFFUSI'  : VRAI si les diffusivites ont change     |
  96 :    |                                                                 |
  97 :    | CHCLIM   : table d'indice 'NEUMANN' et 'DIRICHLET' contenant les|
  98 :    |            Chpoint a n composantes contenant les conditions aux |
  99 :    |            limites de Neumann et Dirichlet par espece.          |
 100 :    |                                                                 |
 101 :    |                                                                 |
 102 :    |-----------------------------------------------------------------|
 103 :    |                         ENTREES-SORTIES                         |
 104 :    |-----------------------------------------------------------------|
 105 :    |                                                                 |
 106 :    | MassEFMH : matrice elementaire EFMH                             |
 107 :    |                                                                 |
 108 :    | Remarque                                                        |
 109 :    | --------                                                        |
 110 :    | On a toujours interet a rentrer Mattm si on l'a et qu'il n'y a  |
 111 :    | pas de modification, afin de conserver les factorisations LU    |
 112 :    | dans solvEFMH en cas de resolution par methode directe.         |
 113 :    |                                                                 |
 114 :    |-----------------------------------------------------------------|
 115 :    |                         SORTIES                                 |
 116 :    |-----------------------------------------------------------------|
 117 :    |                                                                 |
 118 :    | SMTr     : second membre sur les traces                         |
 119 :    |                                                                 |
 120 :    | Matrtr   : Matrice du probleme, rigidite si nouvelle, egale a la|
 121 :    |            matrice entree (mattm ou mattr) sinon                |
 122 :    |                                                                 |
 123 :    |-----------------------------------------------------------------|
 124 : 
 125 : 

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