Afficher cette notice en
1 : $$$$ INITVF NOTICE CHAT 11/09/12 21:16:36 7124 2 : DATE 11/09/12 3 : 4 : Operateur INITVF Voir aussi : 5 : ---------------- 6 : 7 : 8 : ATTENTION La vitesse est optionnelle, L'ordre est important 9 : et les types d'arguments qui se suivent aussi pour tester leur 10 : presence 11 : 12 : Appele par TRANGEOL. Pas pour utilisateurs. 13 : 14 : |-----------------------------------------------------------------| 15 : | Phrase d'appel (en GIBIANE) | 16 : |-----------------------------------------------------------------| 17 : | | 18 : | RESI Matot jaco Mpor2 Mchamt mchamt1 difftot nomespec nbespece | 19 : | nbsource TABRES TABMODI NOUVMAT = INITVF MoDARCY Porosite | 20 : | MateDiff ChPSour DeltaT Cini TetaDiff | 21 : | TetaConv TetaLin (QFACE) QELEM | 22 : | DISPL DISPT CHCLIM optresol ; | 23 : | | 24 : |-----------------------------------------------------------------| 25 : | Generalites : INITVF 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 : champ par elements de composante 'CK' | 39 : | | 40 : | MateDiff : Tenseur de diffusion (type iso, ..) champ par | 41 : | elements de composante 'K' en isotrope, 'K11', 'K21',| 42 : | 'K22' en anisotrope 2d et 'K11', 'K21', 'K22', 'K31'| 43 : | 'K32', 'K33' en anisotrope 3d. Type 'CARACTERISTIQ 44 : | | 45 : | ChPSour : Champ par points des sources volumiques par unite de | 46 : | temps (support maillage centre). Composante ?????? | 47 : | | 48 : | DeltaT : Pas de temps. | 49 : | | 50 : | Cini : Concentration initiale, CHPOINT centre. | 51 : | | 52 : | TetaDiff : Valeur de theta pour theta-schema en temps, operateur| 53 : | de diffusion. Entre 0 et 1. 0 = explicite, 1 = euler | 54 : | implicite. | 55 : | | 56 : | TetaConv : Valeur de theta pour theta-schema en temps, operateur| 57 : | de convection. Entre 0 et 1. 0 = explicite, 1 = Euler| 58 : | implicite. | 59 : | | 60 : | TetaLin : valeur de theta pour theta-schema en temps, operateur| 61 : | lineaire du type coef * C, ou C est l'inconnue. | 62 : | Entre 0 et 1. 0 = explicite, 1 = euler implicite. | 63 : | ??????????? A voir car peut etre identique a Tetadiff| 64 : | | 65 : | | 66 : | Qface : vitesse aux faces, CHPO face de composantes Vx, Vy | 67 : | en 2d et Vx, Vy, Vz en 3d. Il s'agit plus exatement | 68 : | de (V.n)n, c'est a dire de la composante normale de | 69 : | la vitesse aux faces. ???????? (je pressens que | 70 : | castem va sortir des flux, cad integres sur surfaces)| 71 : | | 72 : | CHCLIM : table d'indice 'NEUMANN' et 'DIRICHLET' contenant les| 73 : | Chpoint a n composantes contenant les conditions aux | 74 : | limites de Neumann et Dirichlet par espece. | 75 : | | 76 : | OPTRESOL : Table dont l'entree est optionnelle definissant | 77 : | les options de resolution pour 'KRES'. | 78 : | | 79 : |-----------------------------------------------------------------| 80 : | SORTIES | 81 : |-----------------------------------------------------------------| 82 : | | 83 : | | 84 : | RESI : second membre | 85 : | | 86 : | Matot : matrice globale de discretisation en VF | 87 : | | 88 : | Jaco : matrice globale de discretisation en VF pour le probleme 89 : | stationnaire | 90 : | | 91 : | Mpor : matrice globale de discretisation en VF pour le probleme 92 : | stationnaire | 93 : | | 94 : | Mchamt : Coef permettant de calculer le flux total | 95 : | | 96 : | Mchamt1 : Coef permettant de calculer le flux diffusif | 97 : | | 98 : | Difftot : Coefficient de diffusion totale, integre decentrement| 99 : | | 100 : | Diffdisp : Coefficient de dispersivite | 101 : | | 102 : | nomespc : liste des noms de composante des especes dans Cini | 103 : | | 104 : | nbespece : nombre de composante de Cini, soit nombre d'especes | 105 : | | 106 : | nbsource : nombre de composantes du terme source qd X especes | 107 : | | 108 : | TABRES : Table complete definissant les options de resolution | 109 : | pour 'KRES'. | 110 : | | 111 : | TABMODI : table contenant des logiques indiquant la necessite | 112 : | ou non de reclalculer certains termes. | 113 : | 'POROSITE' : VRAI si le coefficient devant D/DT | 114 : | (porosite) est modifie depuis le dernier| 115 : | appel | 116 : | 'DELTAT' : VRAI si le pas de tps a change | 117 : | 'CONVECTI' : VRAI si la vitesse a change | 118 : | 'COEF_LIN' : VRAI si le coef en facteur de C a change| 119 : | 'DIFFUSI' : VRAI si les diffusivites ont change | 120 : | | 121 : | NOUVMAT : Logique affecte a VRAI lorsque que Matot vient 122 : | d'etre calculee 123 : |-----------------------------------------------------------------| 124 : 125 :
© Cast3M 2003 - Tous droits réservés.
Mentions légales