Accueil | Cast3M

   1 : $$$$ NS       NOTICE  CHAT      11/09/12    21:17:26     7124           
   2 :                                              DATE     11/09/12
   3 :                                                                                 
   4 :    Operateur NS                             Voir aussi :
   5 :     ------------  
   6 :     SYNTAXE ( EQEX ) :   Cf operateur EQEX
   7 :     __________________
   8 : 
   9 :   1/ Formulation non conservative
  10 : 
  11 : a/
  12 :       'OPER'  'NS'   ro un mu              'INCO' 'UN'
  13 : 
  14 : b/
  15 :       'OPER'  'NS'   nu                    'INCO' 'UN'
  16 :       'OPER'  'NS'   nu  s                 'INCO' 'UN'
  17 : 
  18 :      approximation de Boussinesq :
  19 : a/
  20 :       'OPER'  'NS'   ro un mu  gb tn tref  'INCO' 'UN'
  21 : 
  22 : b/
  23 :       'OPER'  'NS'   nu gb tn tref         'INCO' 'UN'
  24 : 
  25 :   2/ Formulation conservative (dilatable)
  26 : 
  27 :       'OPER'  'NS'   mu un    'INCO' 'GN'
  28 :       'OPER'  'NS'   mu un S  'INCO' 'GN'
  29 : 
  30 : 
  31 :                                                                                 
  32 :                                                                                 
  33 :     OBJET :
  34 :     -----
  35 : 
  36 :  Cet operateur discretise les termes de diffusion, de convection et
  37 : eventuellement le terme source de l'equation de Navier - Stokes.
  38 :  Pour une dicretisation element finis EFM1 (algorithme explicite),
  39 : il calcule l'increment.
  40 :  Pour une dicretisation element finis EF (algorithme implicite ou semi
  41 : implicite), il calcule les matrices elementaires et les second membres.
  42 : 
  43 :  Suivant l'option les equations sont traitees sous forme conservative
  44 : ou non conservative.
  45 : 
  46 : 1/ Formulation non conservative
  47 : 
  48 : a/
  49 : ro(du/dt + u Grad u) = mu Lapl u -  Grad p  < + s (=S) >
  50 :    -----                            ------  < + ro*g*beta(T-Tref) >
  51 : 
  52 : b/
  53 : du/dt + u Grad u  = nu Lapl u - 1/ro Grad p  < + s (=S/ro) >
  54 : -----                           -----------  < + g*beta(T-Tref) >
  55 : 
  56 : 2/ Formulation conservative (avec la vitesse massique comme inconnue).
  57 : 
  58 : dG/dt + Div ( u X G ) = mu (Lapl u + 1/3 (Grad Div u))
  59 : -----
  60 :                                                - Grad p  < + S >
  61 :                                                  ------
  62 : (Les termes soulignes ne sont pas discretises dans NS voir procedure EXEC)
  63 : 
  64 :     Commentaires
  65 :     ____________
  66 : 
  67 :  ro,nu,mu  densite, viscosite cinematique (resp. dynamique) moleculaire
  68 :          FLOTTANT ou CHPOINT SCAL CENTRE ou CHPOINT SCAL SOMMET ou MOT
  69 :  s,S     Source volumique de quantite de mouvement. (s=S/ro)
  70 :          POINT ou CHPOINT VECT CENTRE ou MOT
  71 : 
  72 :  approximation de Boussinesq :
  73 :  gb      Coefficient de flottabilite (g*beta ou g est l'accelleration
  74 :          de la pesanteur et beta le coefficient de dilatabilite)
  75 :          POINT ou CHPOINT VECT CENTRE ou MOT
  76 :  tn      Champ de temperature
  77 :          CHPOINT SCAL SOMMET ou MOT
  78 :  tref    temperature de reference
  79 :          FLOTTANT ou CHPOINT SCAL SOMMET ou MOT
  80 : 
  81 :  un      Champ de vitesse transportant
  82 :          CHPOINT VECT SOMMET ou MOT
  83 :  gn      Champ de vitesse massique
  84 :          CHPOINT VECT SOMMET ou MOT
  85 : 
  86 :  Un coefficient de type MOT indique que l'operateur va chercher le
  87 :  coefficient dans la table INCO a l'indice MOT.
  88 : 
  89 : 
  90 :     Options : (EQEX)
  91 :     _________
  92 : 
  93 : 
  94 : 
  95 :  1/ Discretisation EFM1 :  OPTI EFM1
  96 :  ___________________________________
  97 : 
  98 :  Algorithme Explicite
  99 : 
 100 :  La discretisation des termes de convection peut etre :
 101 : 
 102 :  centree                              CENTREE
 103 :  decentree                            SUPG
 104 :  decentree avec capture de choc       SUPGCC   Option par defaut
 105 : 
 106 :  Formulation non conservative         NOCONS   Option par defaut
 107 :  Formulation conservative             CONS
 108 : 
 109 :  2/ Discretisation EF :  OPTI EF
 110 :  _______________________________
 111 : 
 112 :  Algorithme IMPLICITE    OPTI IMPL ou SEMI omega
 113 : 
 114 :  La discretisation des termes de convection peut etre :
 115 : 
 116 :  centree                              CENTREE
 117 :  decentree                            SUPG
 118 :  decentree avec capture de choc       SUPGCC   Option par defaut
 119 :  tenseur visqueux (ordre 2 en temps)  TVISQ
 120 : 
 121 :  Formulation non conservative         NOCONS   Option par defaut
 122 :