Télécharger fimp.notice

Retour à la liste thématique

Afficher cette notice en

Numérotation des lignes :
   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