Accueil | Cast3M

   1 : $$$$ ELNO     NOTICE  CHAT      11/09/12    21:16:01     7124           
   2 :                                              DATE     11/09/12
   3 : 
   4 :   Operateur ELNO                           Voir aussi : MODE DOMA
   5 :     --------------                                        KCHT   PENT  
   6 : 
   7 :     CHP2 = ELNO MODE1 CHP1 | ('CENTRE')                ;
   8 :                            | ('CENTREP1')
   9 :                            | ('MSOMMET')
  10 :                            | ('VOLUMF' GRCHP1 LIMCHP1)                      
  11 : 
  12 : 
  13 : 
  14 : 
  15 : 
  16 :     Objet :
  17 :     _______
  18 : 
  19 :  L'operateur ELNO transforme un CHPOINT definit sur les points CENTRE
  20 : ou CENTREP1 en un CHPOINT definit sur les points SOMMET.
  21 : Cet operateur concerne uniquement les CHPOINTs s'appuyant sur une
  22 : formulation 'NAVIER_STOKES'.
  23 : 
  24 : 
  25 :     Commentaires :
  26 :     ______________
  27 : 
  28 :     MODE1   : Objet modele (type MMODEL). On attend une formulation
  29 :               'NAVIER_STOKES'
  30 : 
  31 :     CHP1    : Objet de type CHPOINT (points CENTRE)
  32 :               L'objet ne doit pas etre partitionne et ne doit avoir
  33 :               qu'une composante
  34 : 
  35 :     CHP2    : Objet de type CHPOINT (points SOMMET)
  36 : 
  37 :     'CENTRE'           (option par defaut).
  38 :       ou
  39 :     'CENTREP1' : Mot-cle du type de CHPOINT de CHP1.
  40 :       ou
  41 :     'MSOMMET'
  42 : 
  43 :     'VOLUMF'   : Mot-cle pour une reconstruction de type Volumes finis
  44 :      Dans ce cas il faut donner les informations suivantes :
  45 : 
  46 :      GRCHP1 : Objet de type CHPOINT (points CENTRE) contenant les 
  47 :               valeurs  du gradient de CHP1.
  48 : 
  49 :      LIMCHP1 : Objet de type CHPOINT (points CENTRE) contenant les 
  50 :                valeurs du limiteur du gradient : Methode MUSCL
  51 :                (voir operateur PENT).
  52 :                
  53 :     Complements d'information :
  54 :     ___________________________
  55 : 
  56 :       On cherche CHP2 minimisant au sens des distributions l'ecart
  57 :     (CHP2 - CHP1) avec la contrainte supplementaire que l'integrale
  58 :     du champ sur le domaine soit conservee.
  59 :       On ne peut appliquer cette transformation qu'aux grandeurs
  60 :     physiques intensives (temperature vitesse pression etc).
  61 :     Cette derniere propriete n'est pas controlee par l'operateur.
  62 :       Les matrices masses sont condensees sur la diagonale.
  63 : 
  64 :     Exemple :
  65 :       KPRES = 'CENTREP1' ;
  66 :       pe=exco rv.inco.'PRESSION' 'PRES' ;
  67 :       pe= kcht $bell scal KPRES pe ;
  68 :       pn= elno $bell pe KPRES ;
  69 : 
  70 :     Cas CHP1 MSOMMET :
  71 :     ------------------
  72 : 
  73 :     Dans ce cas CHP1, connu aux sommets de l'element, est interpole
  74 :     sur les autres noeuds selon l'element lineaire correspondant.
  75 : 
  76 : 
  77 :      Remarques pour les volumes finis
  78 :      ________________________________
  79 : 
  80 :        La valeur en un point sommet est calcule en effectuant la 
  81 :      moyenne des valeurs reconstruites a partir des Centres des 
  82 :      elements auquel il appartient.
  83 : 
  84 :        Si on utilise le MOTCLE 'VOLUMF', la donnee du gradient et 
  85 :      du limiteur est obligatoire.
  86 : 
  87 :      Si le gradient est calcule avec l'operateur PENT et le limiteur
  88 :     impose a 1, la reconstruction au sommet d'un CHPOINT CENTRE 
  89 :     lineaire est exact a l'interieur du domaine. 
  90 : 
  91 :      Si la valeur du limiteur en tout point est 0, la reconstruction 
  92 :     consiste simplement a effectuer la moyenne des valeurs
  93 :     de CHP1 au centre voisins d'un sommet.
  94 : 
  95 :