Afficher cette notice en
1 : $$$$ FIMP NOTICE CHAT 11/09/12 21:16:12 7124 2 : DATE 11/09/12 3 : 4 : Operateur FIMP Voir aussi : NAVI 5 : -------------- 6 : I) 7 : 8 : SYNTAXE (EQEX) : Cf operateur EQEX 9 : _______________ 10 : 11 : 'ZONE' $paroi 'OPER' FIMP COEF 'INCO' 'TN' 12 : 13 : 14 : 15 : 16 : OBJET : 17 : ______ 18 : 19 : Discretise une densite de flux ou une source et calcule l'increment 20 : 21 : EN 2D 22 : elements lignes (SEG2 ou SEG3) -> Flux (en K/ms) 23 : elements massifs (TRI3 TRI7 etc) -> Source volumique (en K/m2s) 24 : EN 3D 25 : elements lignes (SEG2 ou SEG3) -> Pas de sens !! 26 : elements coques (TRI3 TRI7 etc) -> Flux (en K/m2s) 27 : elements massifs (CUB8 CU27 etc) -> Source volumique (en K/m3s) 28 : 29 : Commentaires : 30 : ______________ 31 : 32 : 33 : $paroi Objet MMODEL de type 'NAVIER_STOKES' associe a la 34 : surface sur laquelle porte le flux 35 : 36 : COEF densite de flux CHPOINT SCAL CENTRE 37 : ou CHPOINT SCAL SOMMET 38 : ou FLOTTANT ou MOT 39 : (par convention un flux entrant est compte positivement) 40 : 41 : TN Champ de temperature (en K) CHPOINT SCAL SOMMET 42 : 43 : L'operateur permet de calculer un terme source soit surfacique soit 44 : volumique suivant la nature du support geometrique, pour l'inconnue 45 : TN. Cette inconnue doit etre un CHPOINT SCAL SOMMET. Cependant une 46 : extension est possible si on veut rajouter une source volumique 47 : a l'equation Div U = 0 48 : pour cela on fait porter l'inconnue sur la pression (INCO PRES) 49 : et on precise le support de la pression avec l'option INCOD. 50 : Dans ce cas seul les supports geometriques volumiques sont autorises. 51 : ex : 52 : 'OPTI' 'INCOD' KPRES 53 : 'ZONE' $MT 'OPER' 'FIMP' COEF 'INCO' 'PRES' 54 : ou 55 : KPRES = 'CENTRE' 56 : ou 'CENTREP1' 57 : ou 'MSOMMET' 58 : 59 : 60 : Un coefficient de type MOT indique que l'operateur va chercher le 61 : coefficient dans la table INCO a l'indice MOT. 62 : 63 : 64 : OPTION : (EQEX) 65 : ________ 66 : 67 : Formulation element finis EF et EFM1 (par defaut) 68 : Inconnue duale INCOD SOMMET (option par defaut) 69 : 70 : 71 : II Discretisation des Equations d'Euler 72 : _________________________________________ 73 : 74 : 75 : IIa : gaz parfait mono-constituent polytropique 76 : ________________________________________________ 77 : 78 : 79 : Discretisation en VF "cell-centered" des equations d'Euler pour un gaz 80 : parfait mono-constituent polytropique 81 : 82 : Inconnues: 83 : 84 : densite, quantite de mouvement (qdm), energie totale par unite de volume 85 : (variables conservatives) 86 : 87 : On peut calculer: 88 : 89 : IIa.1 Calcul de la contribution de la force de gravite au residu 90 : 91 : IIa.2 Calcul de la contribution de la force de gravite a la matrice 92 : jacobienne 93 : 94 : 95 : IIa.1 Le residu 96 : ________________ 97 : 98 : 99 : RCHRES = 'FIMP' 'VF' 'GRAVMONO' 'RESI' LISTINCO CHPRN CHPGN CHPGRAV ; 100 : 101 : 102 : LISTINCO : objet de type LISTMOTS 103 : Noms de composantes du resultat (RCHRES) 104 : Il contient dans l'ordre suivant: le noms de la densite, 105 : de la qdm, de l'energie totale par unite de volume 106 : 107 : CHPRN : CHPOINT contenant la masse volumique (une 108 : composante, 'SCAL'). 109 : 110 : CHPGN : CHPOINT contenant le qdm (deux composantes en 2D, 'UX ', 111 : 'UY ', meme SPG que CHPRN). 112 : 113 : CHPGRAV : CHPOINT contenant la gravite (deux composantes en 2D, 114 : 'UX ', 'UY ', meme SPG que CHPRN). 115 : 116 : RCHRES : objet de type CHPOINT (composantes = LISTINCO, meme SPG 117 : que CHPRN) 118 : 119 : 120 : IIa.2 La matrice jacobienne 121 : ___________________________ 122 : 123 : 124 : RJAC = 'FIMP' 'VF' 'GRAVMONO' 'JACOCONS' LISTINCO CHPRN CHPGN CHPGRAV ; 125 : 126 : 127 : LISTINCO : objet de type LISTMOTS 128 : Il contient dans l'ordre suivant: le noms de la densite, 129 : de la qdm, de l'energie totale par unite de volume 130 : 131 : CHPRN : CHPOINT contenant la masse volumique (une 132 : composante, 'SCAL'). 133 : 134 : CHPGN : CHPOINT contenant le qdm (deux composantes en 2D, 'UX ', 135 : 'UY ', meme SPG que CHPRN). 136 : 137 : CHPGRAV : CHPOINT contenant la gravite (deux composantes en 2D, 138 : 'UX ', 'UY ', meme SPG que CHPRN). 139 : 140 : RJAC : objet de type MATRIK 141 : (meme SPG que CHPRN) 142 : (inconnues primales = inconnues duales = LMOT1) 143 : Il contient le jacobien du residu par rapport aux variables 144 : conservatives. 145 : 146 : 147 : 148 : III Discretisation des Equations de Navier-Stokes avec le 149 : modele turbulent k-epsilon 150 : __________________________________________________________ 151 : 152 : 153 : IIIa : gaz thermiquement parfait multi-constituent 154 : ___________________________________________________ 155 : 156 : 157 : Discretisation des Equations de Navier-Stokes multi-constituent 158 : avec le modele turbulent k-epsilon 159 : 160 : Inconnues: 161 : 162 : densite, quantite de mouvement (qdm), energie totale par unite de volume, 163 : densites des especes qui sont dans (TABGAS.'ESPEULE'), 164 : energie cinetique de la turbulence par unite de volume, 165 : taux de dissipation de l'energie turbulente par unite de volume 166 : (variables conservatives) 167 : 168 : On peut calculer: 169 : 170 : IIIa.1 la contribution de la force de la gravite, 171 : de les termes sources des equations de conservation 172 : des especes et des termes sources des equations 173 : d'energie cinetique de la turbulence et de taux de 174 : dissipation de l'energie turbulente au residu 175 : 176 : Le Residu 177 : ____________ 178 : 179 : 180 : RCHRES = 'FIMP' 'VF' 'GRAVKEPS' 'RESI' TABGAS LISTINCO 181 : CHPRN CHPGN CHPGRAV CHPRKA CHPREPS GRADR GRADV ; 182 : 183 : TABGAS : la table vec des proprietes de gas; 184 : ici on a besoin de: 185 : a) le nombre \sigma_t (vois les 186 : terms de la gravite dans 187 : les equations k-\eps) 188 : nom: TABGAS . 'SIGT' 189 : b) le nombre C_{\eps 1} (vois 190 : l'equation de \eps) 191 : nom: TABGAS . 'CEPS1' 192 : c) le nombre C_{\eps 2} (vois 193 : l'equation de \eps) 194 : nom: TABGAS . 'CEPS2' 195 : d) les taux des productions des especes 196 : ( le term source); type 'CHPOINT' 197 : (SPG = 'CENTRE', le nombre de composants 198 : egal a nombre des especes dans 199 : TABGAS . 'ESPEULE') 200 : e) TABGAS . 'ESPEULE' - les noms des 201 : especes qui sont explicitement 202 : dans les equations 203 : f) TABGAS . 'MUTURB' - la viscosite turbulente 204 : (SPG = 'CENTRE'; number de compos. = 1) 205 : g) les enthalpies des formations 206 : des especes a la temperature T=0K 207 : 208 : 209 : LISTINCO : objet de type LISTMOTS 210 : Noms de composantes du resultat (RCHRES) 211 : Il contient dans l'ordre suivant: le noms de la densite, 212 : de la qdm, de l'energie totale par unite de volume, 213 : des especes qui sont dans (TABGAS.'ESPEULE'), 214 : d'energie cinetique de la turbulence par unite de volume, 215 : de taux de dissipation de l'energie turbulente par unite de 216 : volume 217 : 218 : CHPRN : CHPOINT contenant la masse volumique (une 219 : composante, 'SCAL'). 220 : 221 : CHPGN : CHPOINT contenant le qdm (deux composantes en 2D, 'UX ', 222 : 'UY ', meme SPG que CHPRN). 223 : 224 : CHPGRAV : CHPOINT contenant la gravite (deux composantes en 2D, 225 : 'UX ', 'UY ', meme SPG que CHPRN). 226 : 227 : CHPRKA : CHPOINT contenant l'energie cinetique turbulente 228 : par unite de volume (une composante, meme SPG que 229 : CHPRN). 230 : 231 : CHPREPS : CHPOINT contenant le taux de dissipation 232 : d'energie de la turbulence par unite de volume 233 : (une composante, meme SPG que CHPRN). 234 : 235 : GRADR : gradient de la densite 236 : (SPG = 'CENTRE', composantes: 'P1DX','P1DY',('P1DZ')) 237 : 238 : GRADV : gradient de la vitesse 239 : (SPG = 'CENTRE', composantes: 'P1DX','P1DY', ('P1DZ') 240 : 'P2DX','P2DY', ('P2DZ') 241 : ('P3DX','P3DY', ('P3DZ'))) 242 : 243 : RCHRES : objet de type CHPOINT (composantes = LISTINCO, meme SPG 244 : que CHPRN) 245 : 246 :
© Cast3M 2003 - Tous droits réservés.
Mentions légales