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