Télécharger initEFMH.notice

Retour à la liste thématique

Afficher cette notice en

Numérotation des lignes :
   1 : $$$$ INITEFMH NOTICE  CHAT      11/09/12    21:16:36     7124           
   2 :                                              DATE     11/09/12
   3 : 
   4 :   Operateur INITEFMH                   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 :     Cette procedure est appelee par TRANGEOL.
  12 :     Pas de notice car pas pour utilisateur lambda.
  13 : 
  14 :   |-----------------------------------------------------------------| 
  15 :   | Phrase d'appel (en GIBIANE)                                     |
  16 :   |-----------------------------------------------------------------|
  17 :   |                                                                 |
  18 :   | SMTr MatrTr CoefDt TbDarTra MassEFMH nomespec                   |
  19 :   | nbespece Difftot Tcini TABRES TABMODI= INITEFMH MoDARCY Porosite|
  20 :   |                      MateDiff ChPSour DeltaT Cini TetaDiff      |
  21 :   |                      TetaConv TetaLin fluimp dircli (QFACE)     |
  22 :   |                      LMLump DECENTR CHCLIM optresol ;           |
  23 :   |                                                                 |
  24 :   |-----------------------------------------------------------------|
  25 :   | Generalites : MATTEFMH construit la matrice de discretisation   |
  26 :   |               du probleme de transport convection-diffusion pour|
  27 :   |               le premier pas de tps d'un algorithme transitoire.|
  28 :   |               Le second membre et les Conditions limites de flux|
  29 :   |               sont pris en compte.                              |
  30 :   |            RESTE TCINI, DECENTR et TERME LIN                    |
  31 :   |-----------------------------------------------------------------|
  32 :   |                                                                 |
  33 :   |-----------------------------------------------------------------|
  34 :   |                         ENTREES                                 |
  35 :   |-----------------------------------------------------------------|
  36 :   | MoDARCY  : modele Darcy.                                        |
  37 :   |                                                                 |
  38 :   | Porosite : champoint de composante 'CK'                         | 
  39 :   |                                                                 |
  40 :   | MateDiff : Tenseur de diffusion  (type iso, ..) champoint       |
  41 :   |             de composante 'K' en isotrope, 'K11', 'K21',        |
  42 :   |            'K22' en anisotrope 2d et  'K11', 'K21', 'K22', 'K31'|
  43 :   |            'K32', 'K33' en anisotrope 3d. Type 'CARACTERISTIQUE'|
  44 :   |                                                                 |
  45 :   | DeltaT   : Pas de temps.                                        |
  46 :   |                                                                 |
  47 :   | ChPSour  : Champ par points des sources volumiques par unite de |
  48 :   |            temps (support maillage centre). Composante ??????   |
  49 :   |                                                                 |
  50 :   | Cini     : Concentration initiale, CHPOINT centre.              |
  51 :   |            Composante 'H'.                                      |
  52 :   |                                                                 |
  53 :   | Qface    : vitesse aux faces, CHPO face de composantes Vx, Vy   |
  54 :   |            en 2d et Vx, Vy, Vz en 3d. Il s'agit plus exatement  |
  55 :   |            de (V.n)n, c'est a dire de la composante normale de  |
  56 :   |            la vitesse aux faces. ???????? (je pressens que      |
  57 :   |            castem va sortir des flux, cad integres sur surfaces)|
  58 :   |                                                                 |
  59 :   | TetaDiff : Valeur de theta pour theta-schema en temps, operateur|
  60 :   |            de diffusion. Entre 0 et 1. 0 = explicite, 1 = euler |
  61 :   |            implicite.                                           |
  62 :   |                                                                 |
  63 :   | TetaConv : Valeur de theta pour theta-schema en temps, operateur|
  64 :   |            de convection. Entre 0 et 1. 0 = explicite, 1 = Euler|
  65 :   |            implicite.                                           |
  66 :   |                                                                 |
  67 :   | TetaLin  : valeur de theta pour theta-schema en temps, operateur|
  68 :   |            lineaire du type coef * C, ou C est l'inconnue.      |
  69 :   |            Entre 0 et 1. 0 = explicite, 1 = euler implicite.    |
  70 :   |            ??????????? A voir car peut etre identique a Tetadiff|
  71 :   |                                                                 | 
  72 :   | LMLump   : Logique. Si vrai on effectue une condensation de     |
  73 :   |            masse de la matrice EFMH                             |
  74 :   |                                                                 |
  75 :   | DECENTR  : Logique. Vrai veut dire schemas decentres et faux    |
  76 :   |            veut dire schema convectif centre.                   |        
  77 :   |                                                                 |          
  78 :   | CHCLIM   : table d'indice 'NEUMANN' et 'DIRICHLET' contenant les|
  79 :   |            Chpoint a n composantes contenant les conditions aux |
  80 :   |            limites de Neumann et Dirichlet par espece.          |
  81 :   |            L'indice 'FLUXTOT' contient les conditions limites   |
  82 :   |            de flux total et 'FLUMIXTE' concerne une condition   |
  83 :   |            de flux mixte : 'FLUMIXTE' . 'VAL' contient le champ |
  84 :   |            a n composantes indiquant le flux, 'FLUMIXTE' . 'A'  |
  85 :   |            et 'FLUMIXTE' . 'B' les coef (champoints SCAL) tels  |
  86 :   |            que A D grad (C) + B (C) = VAL                       | 
  87 :   |                                                                 |
  88 :   | OPTRESOL : Table dont l'entree est optionnelle definissant      |     
  89 :   |               les options de resolution pour 'KRES'.            |        
  90 :   |                                                                 |
  91 :   |-----------------------------------------------------------------|
  92 :   |                         SORTIES                                 |
  93 :   |-----------------------------------------------------------------|
  94 :   |                                                                 |
  95 :   |                                                                 |
  96 :   | MassEFMH : matrice elementaire EFMH                             |
  97 :   |                                                                 |
  98 :   | MatrTr   : matrice globale sur les traces                       |
  99 :   |                                                                 |
 100 :   | SMTr     : second membre sur les traces                         |
 101 :   |                                                                 |
 102 :   | TbDarTra : table Darcy transitoire utilisee par MHYB, SMTP ...  |
 103 :   |                                                                 |
 104 :   | nomespec : liste des noms de composante des especes dans Cini   |
 105 :   |                                                                 |
 106 :   | nbespece : nombre de composante de Cini, soit nombre d'especes  |
 107 :   |                                                                 |
 108 :   | nbsource : nombre de composantes du terme source qd X especes   |
 109 :   |                                                                 |
 110 :   | Diffdisp : Dipersivite, tenseur chpoint K11 K22 K33 K21 K31 K32 |
 111 :   |                                                                 |
 112 :   | TABRES   : Table complete definissant les options de resolution |    
 113 :   |            pour 'KRES'.                                         |
 114 :   |                                                                 |
 115 :   | Tcini    : Trace de concentration aux faces (eventuellement a   |
 116 :   |            plusieurs composantes (especes)                      |
 117 :   |                                                                 |
 118 :   | TABMODI  : table contenant des logiques indiquant la necessite  |
 119 :   |            ou non de reclalculer certains termes.               |
 120 :   |            'POROSITE' : VRAI si le coefficient devant D/DT      |
 121 :   |                         (porosite) est modifie depuis le dernier|
 122 :   |                         appel                                   |
 123 :   |            'DELTAT'   : VRAI si le pas de tps a change          |
 124 :   |            'CONVECTI' : VRAI si la vitesse a change             |
 125 :   |            'COEF_LIN' : VRAI si le coef en facteur de C a change|
 126 :   |            'DIFFUSI'  : VRAI si les diffusivites ont change     |
 127 :   |                                                                 |
 128 :   |-----------------------------------------------------------------|
 129 :   |                  VARIABLES INTERNES                             |
 130 :   |-----------------------------------------------------------------|
 131 :   |                                                                 |
 132 :   | CoefDt   : coeff devant dC/dt integre sur les elements          |
 133 :   |                                                                 |
 134 :   | PCONV    : Logique indiquant VRAI si presence de convection     |
 135 :   |                                                                 |
 136 :   | toltheta : 1.D-4  seuil en dessous duquel on considere que la   |
 137 :   |            valeur de theta du theta-schema est nulle (schema    |
 138 :   |            explicite) ou au contraire euler-implicite si        |
 139 :   |            theta > 0.9999                                       |
 140 :   |                                                                 |
 141 :   | Tccini   : Trace de concentration aux faces (une composante)    |
 142 :   |                                                                 |
 143 :   | SSource  : Source aux centre (une composante)                   |
 144 :   |                                                                 |
 145 :   | CCini    : concentration aux centres (une composante)           |
 146 :   |                                                                 |
 147 :   | lstcps   : liste des noms de composante des especes dans Chpsour|
 148 :   |                                                                 |
 149 :   | SSMTr     : second membre sur les traces pour une espece        |
 150 :   |                                                                 |
 151 :   | MatConv   : matrice globale sur les traces pour la convection   |
 152 :   |                                                                 |
 153 :   | Numdiff  : diffusivite numerique due au decentrement            |
 154 :   |                                                                 |
 155 :   | FLUNEU   : LOGIQUE valant VRAI si conditions de Neumann         |
 156 :   |                                                                 |
 157 :   | CLFLUX   : Chpoint a n composantes contenant les flux imposes   |
 158 :   |            pour chaque espece chimique. nul si pas de flux      |
 159 :   |            OPTIONNEL                                            |
 160 :   |                                                                 |
 161 :   *******************************************************************
 162 : 

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