Télécharger fimp.notice

Retour à la liste

Afficher cette notice en

Numérotation des lignes :
$$$$ FIMP     NOTICE  CHAT      11/09/12    21:16:12     7124           
                                             DATE     11/09/12
                                                                                
    Operateur FIMP                           Voir aussi : NAVI
    --------------  
    I)                                                                          

    SYNTAXE (EQEX) : Cf operateur EQEX
    _______________

  'ZONE' $paroi 'OPER' FIMP  COEF 'INCO' 'TN'

                                                                                

                                                                                
    OBJET :
    ______

    Discretise une densite de flux ou une source et calcule l'increment

        EN 2D
        elements lignes (SEG2 ou SEG3)    -> Flux (en K/ms)
        elements massifs (TRI3 TRI7 etc)  -> Source volumique (en K/m2s)
        EN 3D
        elements lignes (SEG2 ou SEG3)    -> Pas de sens !!
        elements coques (TRI3 TRI7 etc)   -> Flux (en K/m2s)
        elements massifs (CUB8 CU27 etc)  -> Source volumique (en K/m3s)

     Commentaires :
     ______________
                                                                                

     $paroi   Objet MMODEL de type 'NAVIER_STOKES' associe a la
              surface sur laquelle porte le flux

     COEF     densite de flux                CHPOINT SCAL CENTRE
                                          ou CHPOINT SCAL SOMMET
                                          ou FLOTTANT ou MOT
              (par convention un flux entrant est compte positivement)

     TN       Champ de temperature (en K)      CHPOINT SCAL SOMMET

 L'operateur permet de calculer un terme source soit surfacique soit
 volumique suivant la nature du support geometrique, pour l'inconnue
 TN. Cette inconnue doit etre un CHPOINT SCAL SOMMET. Cependant une
 extension est possible si on veut rajouter une source volumique
 a l'equation Div U = 0
  pour cela on fait porter l'inconnue sur la pression (INCO PRES)
 et on precise le support de la pression avec l'option INCOD.
 Dans ce cas seul les supports geometriques volumiques sont autorises.
 ex :
 'OPTI' 'INCOD' KPRES
 'ZONE' $MT 'OPER' 'FIMP' COEF 'INCO' 'PRES'
 ou
  KPRES = 'CENTRE'
  ou      'CENTREP1'
  ou      'MSOMMET'


 Un coefficient de type MOT indique que l'operateur va chercher le
 coefficient dans la table INCO a l'indice MOT.


    OPTION : (EQEX)
    ________

    Formulation element finis EF et EFM1 (par defaut)
    Inconnue duale INCOD SOMMET  (option par defaut)

   
 II Discretisation des Equations d'Euler 
 _________________________________________


 IIa : gaz parfait mono-constituent polytropique
 ________________________________________________


 Discretisation en VF "cell-centered" des equations d'Euler pour un gaz 
 parfait mono-constituent polytropique

 Inconnues:

 densite, quantite de mouvement (qdm), energie totale par unite de volume 
 (variables conservatives)

 On peut calculer:

 IIa.1 Calcul de la contribution de la force de gravite au residu

 IIa.2 Calcul de la contribution de la force de gravite a la matrice 
       jacobienne


 IIa.1 Le residu
 ________________


 RCHRES = 'FIMP' 'VF' 'GRAVMONO' 'RESI' LISTINCO CHPRN CHPGN CHPGRAV ;


 LISTINCO : objet de type LISTMOTS
            Noms de composantes du resultat (RCHRES)
            Il contient dans l'ordre suivant: le noms de la densite,
            de la qdm, de l'energie totale par unite de volume

 CHPRN    : CHPOINT  contenant la masse volumique (une 
            composante, 'SCAL').

 CHPGN    : CHPOINT  contenant le qdm (deux composantes en 2D, 'UX  ',
            'UY  ', meme SPG que CHPRN).

 CHPGRAV  : CHPOINT  contenant la gravite (deux composantes en 2D, 
            'UX  ', 'UY  ', meme SPG que CHPRN).

 RCHRES   : objet de type CHPOINT (composantes =  LISTINCO, meme SPG 
            que CHPRN)


 IIa.2 La matrice jacobienne
 ___________________________


 RJAC  = 'FIMP' 'VF' 'GRAVMONO' 'JACOCONS' LISTINCO CHPRN CHPGN CHPGRAV ;


 LISTINCO : objet de type LISTMOTS
            Il contient dans l'ordre suivant: le noms de la densite,
            de la qdm, de l'energie totale par unite de volume

 CHPRN    : CHPOINT  contenant la masse volumique (une 
            composante, 'SCAL').

 CHPGN    : CHPOINT  contenant le qdm (deux composantes en 2D, 'UX  ',
            'UY  ', meme SPG que CHPRN).

 CHPGRAV  : CHPOINT  contenant la gravite (deux composantes en 2D, 
            'UX  ', 'UY  ', meme SPG que CHPRN).

 RJAC     : objet de type MATRIK
           (meme SPG que CHPRN)
           (inconnues primales = inconnues duales = LMOT1)
           Il contient le jacobien du residu par rapport aux variables 
           conservatives.
       


 III Discretisation des Equations de Navier-Stokes avec le
      modele turbulent k-epsilon          
 __________________________________________________________


 IIIa : gaz thermiquement parfait multi-constituent 
 ___________________________________________________


 Discretisation des Equations de Navier-Stokes multi-constituent
  avec le modele turbulent k-epsilon  

 Inconnues:

 densite, quantite de mouvement (qdm), energie totale par unite de volume,
 densites des especes qui sont dans (TABGAS.'ESPEULE'),
 energie cinetique de la turbulence par unite de volume,
 taux de dissipation de l'energie turbulente par unite de volume
 (variables conservatives)

 On peut calculer:

 IIIa.1 la contribution de la force de la gravite,
        de les termes sources des equations de conservation
        des especes et des termes sources des equations 
        d'energie cinetique de la turbulence et de taux de 
        dissipation de l'energie turbulente au residu

 Le Residu
 ____________


 RCHRES = 'FIMP' 'VF' 'GRAVKEPS' 'RESI' TABGAS LISTINCO
          CHPRN CHPGN CHPGRAV CHPRKA CHPREPS GRADR GRADV ;

 TABGAS  :  la table vec des proprietes de gas;
            ici on a besoin de:
                        a) le nombre \sigma_t (vois les
                           terms de la gravite dans
                           les equations k-\eps)
                           nom: TABGAS . 'SIGT'
                        b) le nombre C_{\eps 1} (vois 
                           l'equation  de \eps)
                           nom: TABGAS . 'CEPS1'
                        c) le nombre C_{\eps 2} (vois
                           l'equation de  \eps)
                           nom: TABGAS . 'CEPS2'
                        d) les taux des  productions des especes
                           ( le term source); type 'CHPOINT'
                           (SPG = 'CENTRE', le nombre de composants
                            egal a nombre des especes dans
                            TABGAS . 'ESPEULE')
                        e) TABGAS . 'ESPEULE' - les noms des
                           especes qui sont explicitement
                           dans les equations
                        f) TABGAS . 'MUTURB' -  la viscosite turbulente
                           (SPG = 'CENTRE'; number de compos. = 1)
                        g) les enthalpies des formations
                           des especes a la temperature T=0K
                        

 LISTINCO : objet de type LISTMOTS
            Noms de composantes du resultat (RCHRES)
            Il contient dans l'ordre suivant: le noms de la densite,
            de la qdm, de l'energie totale par unite de volume,
            des especes qui sont dans (TABGAS.'ESPEULE'),
            d'energie cinetique de la turbulence par unite de volume,
            de taux de dissipation de l'energie turbulente par unite de
            volume

 CHPRN    : CHPOINT  contenant la masse volumique (une 
            composante, 'SCAL').

 CHPGN    : CHPOINT  contenant le qdm (deux composantes en 2D, 'UX  ',
            'UY  ', meme SPG que CHPRN).

 CHPGRAV  : CHPOINT  contenant la gravite (deux composantes en 2D, 
            'UX  ', 'UY  ', meme SPG que CHPRN).

 CHPRKA   : CHPOINT contenant l'energie cinetique turbulente 
            par unite de volume (une composante, meme SPG que
            CHPRN).

 CHPREPS  : CHPOINT contenant le taux de dissipation
            d'energie de la turbulence par unite de volume
            (une composante, meme SPG que CHPRN). 
 
 GRADR    : gradient de la densite
            (SPG = 'CENTRE', composantes: 'P1DX','P1DY',('P1DZ'))

 GRADV    : gradient de la vitesse
            (SPG = 'CENTRE', composantes: 'P1DX','P1DY', ('P1DZ')
                                          'P2DX','P2DY', ('P2DZ')
                                         ('P3DX','P3DY', ('P3DZ')))   
           
 RCHRES   : objet de type CHPOINT (composantes =  LISTINCO, meme SPG 
            que CHPRN)


© Cast3M 2003 - Tous droits réservés.
Mentions légales