Télécharger pent.notice

Retour à la liste thématique

Afficher cette notice en

Numérotation des lignes :
   1 : $$$$ PENT     NOTICE  CHAT      11/09/12    21:17:35     7124           
   2 :                                              DATE     11/09/12
   3 :                                              
   4 :    Operateur PENT                           Voir aussi :
   5 :     --------------  
   6 : 
   7 :     Objet :
   8 :     _______
   9 : 
  10 :     Évaluation du gradient d'un champ dans le cadre d'une 
  11 :     discretisation de type volumes finis (variables aux 
  12 :     centres)
  13 : 
  14 :     -----------------------------------------------------------
  15 :     | 1ere possibilite : creation d'un gradient aux CENTRES   |
  16 :     -----------------------------------------------------------
  17 : 
  18 :    RCHPO1 RCHPO2 RCHELEM1 = 'PENT' MOD1
  19 :            'CENTRE' MCLE1 MCLE2 LMOT1 CHPO1 ('CLIM' CHPO2) ;
  20 :  
  21 :    ou
  22 :  
  23 :    RCHPO1 RCHPO2  = 'PENT' MOD1
  24 :            'CENTRE' MCLE1 MCLE2 LMOT1 CHPO1 ('CLIM' CHPO2) 
  25 :            'GRADGEO' RCHELEM1 ;
  26 : 
  27 : 
  28 :     Commentaire :
  29 :     _____________
  30 : 
  31 :     MOD1    : Objet MODELE.
  32 : 
  33 :     MCLE1    : MOT; indique la façon de considerer la frontiere;
  34 :                4 choix possibles:
  35 :                * 'BORDNULL': reconstruction lineaire exacte; 
  36 :                   le gradient du CHPOINT est nul sur les elements 
  37 :                   de frontiere;
  38 :                * 'LINEXACT': reconstruction lineaire exacte; le 
  39 :                   gradient est calcule sur les elements de frontiere
  40 :                   par interpolation lineaire exacte.
  41 :                * 'EULESCAL': reconstruction lineaire exacte; le 
  42 :                   gradient est calcule en utilisant des conditions 
  43 :                   aux limites de type mur pour un champ scalaire 
  44 :                   (etat miroir a l'element de bord). Ceci dans le 
  45 :                   cadre des equations d'Euler.
  46 :                * 'EULEVECT': reconstruction lineaire exacte; le 
  47 :                   gradient est calcule en utilisant des conditions 
  48 :                   aux limites de type mur pour un champ vectoriel 
  49 :                   (etat miroir a l'element de bord). Ceci dans le 
  50 :                   cadre des equations d'Euler.
  51 : 
  52 :     MCLE2    : MOT; indique le type de limiteur de gradient a 
  53 :                calculer. 
  54 :                * 'LIMITEUR', on calcule le limiteur de Barth-Jespersen; 
  55 :                * 'NOLIMITE', les coefficients de limiteur sont egal a
  56 :                  1.0
  57 : 
  58 :     LMOT1    : LISTMOTS, composantes de CHPO1 et CHPO2
  59 : 
  60 :     CHPO1    : CHPOINT 'CENTRE' (i composantes, 1 <= i <= 9) dont on
  61 :                souhaite calculer le gradient.
  62 : 
  63 :     CHPO2    : CHPOINT (meme composantes que CHPO1): champoint
  64 :                qui specifie les conditions limites de type Dirichlet 
  65 :                sur certains points de type 'FACE'
  66 :               
  67 :     RCHELEM1 : Champ par element des coefficients geometriques pour le
  68 :                calcul du gradient
  69 : 
  70 :     RCHPO1   : CHPOINT 'CENTRE' (NDIM * i composantes); contient le
  71 :                gradient du CHPO1; le gradient associe a la i-eme
  72 :                composante a pour noms de composantes 'PiDX', 'PiDY'
  73 :                ('PiDZ').
  74 : 
  75 :     RCHPO2  :  CHPOINT 'CENTRE'  (i composantes); contient les
  76 :                coefficients multiplicateurs compris entre 0 et 1 par
  77 :                lesquels il faut multiplier le gradient si on souhaite
  78 :                que ce dernier soit limite. Le nom des composantes est 
  79 :                'Pi', avec la meme convention que pour RCHPO1.
  80 : 
  81 : 
  82 :    Remarques :
  83 :    __________
  84 : 
  85 :    1) Le gradient calcule est exact a l'interieur du domaine si la 
  86 :       fonction est lineaire. Cette propriete est vraie egalement sur 
  87 :       le bord avec l'option 'LINEXACT'.
  88 : 
  89 :    2) Les options 'EULESCAL' et 'EULEVECT' traite la frontiere du 
  90 :       domaine comme un mur.
  91 : 
  92 :    3) Si on utilise l'option 'EULEVECT', CHPO1 (et CHPO2) doit avoir 
  93 :       deux composantes en 2D ('UX','UY') et trois composantes en 3D 
  94 :       ('UX','UY','UZ')
  95 : 
  96 : 
  97 :     -----------------------------------------------------------
  98 :     | 2eme possibilite : creation d'un gradient aux FACEs     |
  99 :     -----------------------------------------------------------
 100 : 
 101 : 
 102 :  
 103 :    RCHPO1  RCHELEM1 = 'PENT' MOD1 'FACE' 'DIAMAN2'  LMOT1 LMOT2
 104 :                        CHPO1 CHPO2 CHPO3 ;
 105 : 
 106 :    ou
 107 : 
 108 :    RCHPO1  = 'PENT' MOD1 'FACE' 'DIAMAN2' LMOT1 LMOT2
 109 :               CHPO1 CHPO2 CHPO3 'GRADGEO' RCHELEM1 ;
 110 : 
 111 : 
 112 :    Commentaire :
 113 :    _____________
 114 : 
 115 : 
 116 : 
 117 :    MOD1     :  Objet MODELE.
 118 : 
 119 :    LMOT1    :  LISTMOTS, composantes de CHPO1 et CHPO2
 120 : 
 121 :    LMOT2    :  LISTMOTS, composantes de CHPO3 et RCHPO1
 122 : 
 123 :    CHPO1    :  CHPOINT 'CENTRE' dont on
 124 :                souhaite calculer le gradient.
 125 : 
 126 :    CHPO2    :  CHPOINT qui specifie les conditions limites de type
 127 :                Dirichlet sur certains points de type 'FACE'
 128 :              
 129 :    CHPO3    :  CHPOINT qui specifie les conditions limites de type
 130 :                von Neumann sur certains points de type 'FACE'
 131 :              
 132 :    RCHELEM1 :  Champ par element des coefficients geometriques pour le
 133 :                calcul du gradient.
 134 : 
 135 :    RCHPO1   :  CHPOINT 'FACE' (NDIM * i composantes); contient le
 136 :                gradient du CHPO1
 137 : 
 138 :    Remarques :
 139 :    __________
 140 : 
 141 :    1) La condition limite de type von Neumann prise en compte est donne 
 142 :       par le produit scalaire de CHPO3 et des normales aux faces
 143 :    
 144 : 
 145 : 
 146 :     --------------------------------------------------------------------
 147 :     | 3eme possibilite : creation d'un gradient aux FACEs  
 148 :     | en 2 dimensions avec tenseur symetrique                                   
 149 :       -------------------------------------------------------------------
 150 : 
 151 : 
 152 :  
 153 :    RCHPO1  RCHELEM1 = 'PENT' 'FACE' MCLE1 MOD1 CHPO1 
 154 :            ('DISPDIF CHPO3) ('CLIM' CHPO2) 
 155 :            ('NEUM' CHPO4) ('MIXT' CHPO5) ;
 156 : 
 157 :    ou
 158 : 
 159 :    RCHPO1  = 'PENT' 'FACE'  MCLE1 MOD1 CHPO1 ('DISPDIF CHPO3) 
 160 :              ('CLIM' CHPO2) ('NEUM' CHPO4) ('MIXT' CHPO5) 
 161 :               'GRADGEO' RCHELEM1 ;
 162 : 
 163 : 
 164 :    Commentaire :
 165 :    _____________
 166 : 
 167 : 
 168 : 
 169 :    MOD1     :  Objet MODELE.
 170 : 
 171 :    MCLE1    :  Methode pour le calcul du gradient. Options possibles : 
 172 :                'MPFA'
 173 : 
 174 :    CHPO1    :  CHPOINT 'CENTRE'  dont on souhaite calculer le gradient.
 175 : 
 176 :    CHPO2    :  CHPOINT 'FACE' (meme composantes que CHPO1): 
 177 :                champoint qui specifie les conditions aux limites de type
 178 :                Dirichlet sur certains points de type 'FACE'
 179 : 
 180 :    CHPO4    :  CHPOINT 'FACE' (1 composante)  
 181 :                champoint qui specifie les conditions aux limites de type
 182 :                flux sur certains points de type 'FACE'
 183 : 
 184 :    CHPO5    :  CHPOINT 'FACE' (4 composantes, lambda1,lambda2,qlimx,qlimy)  
 185 :                champoint qui specifie les conditions aux limites de type
 186 :                mixtes sur certains points de type 'FACE'
 187 :                lambda1 (d grad T . n) + lambda2 T = (qlimx*nx) + (qlim*ny)
 188 : 
 189 :    CHPO3    :  CHPOINT 'CENTRE' (3 composantes 'K11','K22','K21') 
 190 :                du tenseur 
 191 :              
 192 :    RCHELEM1 :  Champ par element des coefficients geometriques pour le
 193 :                calcul du gradient.
 194 : 
 195 :    RCHPO1   :  CHPOINT 'FACE' : contient le produit scalaire du gradient 
 196 :                de CHPO1 avec les normales aux faces du maillage.  
 197 :                Il a pour nom de composante 'FLUX'.
 198 : 
 199 : 
 200 :     --------------------------------------------------------------------
 201 :     | 4eme possibilite : creation d'un gradient aux FACEs  
 202 :     |  avec un tenseur symetrique, en suivant la methode decrite
 203 :     | par Christophe Le Potier dans "Schema volumes finis pour des operateurs 
 204 :     |  de diffusion fortement anisotropes sur des maillages non structures,
 205 :     |  C. R. Acad. Sci. Ser. I \textbf{340}, 2005, pp. 921--926." 
 206 :     |  Le schema obtenu est symetrique.
 207 :     |  -------------------------------------------------------------------
 208 : 
 209 : 
 210 :    RCHPO1  RCHELEM1 = 'PENT' 'FACE' MCLE1 MOD1 CHPO1 
 211 :            ('DISPDIF CHPO3) ('CLIM' CHPO2) 
 212 :            ('NEUM' CHPO4) ;
 213 : 
 214 :    Commentaire :
 215 :    _____________
 216 : 
 217 : 
 218 : 
 219 :    MOD1     :  Objet MODELE.
 220 : 
 221 :    MCLE1    :  Methode pour le calcul du gradient. Options possibles : 
 222 :                'VFSYM'
 223 : 
 224 :    CHPO1    :  CHPOINT 'CENTRE'  dont on souhaite calculer le gradient.
 225 : 
 226 :    CHPO2    :  CHPOINT 'SOMMET' (meme composantes que CHPO1): 
 227 :                champoint qui specifie les conditions aux limites de type
 228 :                Dirichlet sur certains points de type 'SOMMET'
 229 : 
 230 :    CHPO4    :  CHPOINT 'FACE' (1 composante)  
 231 :                champoint qui specifie les conditions aux limites de type
 232 :                flux sur certains points de type 'FACE'
 233 : 
 234 :    RCHELEM1 :  Champ par element des coefficients geometriques pour le
 235 :                calcul du gradient (en 2 dimensions).
 236 : 
 237 :    RCHPO1   :  CHPOINT 'FACE' : contient le produit scalaire du gradient 
 238 :                de CHPO1 avec les normales aux faces du maillage.  
 239 :                Il a pour nom de composante 'FLUX'.
 240 : 

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