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