Accueil | Cast3M

   1 : $$$$ FPU      NOTICE  CHAT      11/09/12    21:16:19     7124           
   2 :                                              DATE     11/09/12
   3 :                                                                                 
   4 :    Operateur FPU                           Voir aussi :
   5 :      -------------  
   6 :      Syntaxe (EQEX) : Cf Operateur EQEX
   7 :      ________________
   8 : 
   9 : 
  10 :  1ere syntaxe
  11 :  
  12 :      Formulation EFM1 :
  13 :      
  14 :         'OPER' 'FPU' NU UET YP           'INCO' 'UN' 'KN' 'EN'
  15 :  
  16 :  2eme syntaxe
  17 : 
  18 :      Formulation EF :
  19 :      
  20 :         'OPER' 'FPU' RO UN MU UET YP     'INCO' 'UN' <'KN' 'EN'>
  21 : 
  22 :  3eme syntaxe
  23 : 
  24 :      Formulation EF :
  25 :      
  26 :         'OPER' 'FPU' RO UN MU UET $mt NUEFF 'INCO' 'UN'
  27 : 
  28 : 
  29 : 
  30 :                                                                                 
  31 :                                                                                 
  32 :      Objet :
  33 :      _______
  34 :                                                                                 
  35 :  I/ 1ere syntaxe - Formulation EFM1
  36 : 
  37 :      Discretise une condition de tension a la paroi suivant un modele de
  38 :      fonction de paroi (modele de longueur de melange de Van Driest [1])
  39 :      La solution U+ de l'equation est tabulee pour Y+ = 1 a 100. Loi Log
  40 :      standard au dela.
  41 :      Pour l'equation de QDM la condition est une condition de Neumann
  42 :      sur la frontiere maillee du domaine. La fonction de paroi modelise
  43 :      une partie de l'ecoulement du fluide (zone d'epaisseur Yp) qui se
  44 :      fait en dehors du domaine maille. Il convient d'en tenir compte
  45 :      lorsqu'on fait des bilans.
  46 : 
  47 :      Les conditions limites correspondantes sur K et Epsilon sont
  48 :      calculees si presence des inconnues 'KN' et 'EN'.
  49 : 
  50 :      les valeurs de K et Epsilon sont imposees comme des conditions de
  51 :      Dirichlet.
  52 :      Dans le cas EFM1 les inconnues K et Epsilon sont obligatoires.
  53 : 
  54 :  II/ 2eme syntaxe - Formulation EF
  55 : 
  56 :      Les fonctionnalites de cette option different sur deux points de la
  57 :      precedente. La fonction de paroi est implicitee (methode iterative)
  58 :      et la loi U+ = F(Y+) est donnee par la loi Reichardt valable jusqu'a
  59 :      Y+ = 0. Pour le reste cette option est identique a la precedente.
  60 :      Les conditions limites correspondantes sur K et Epsilon sont
  61 :      calculees si presence des inconnues 'KN' et 'EN'. La loi de paroi est
  62 :      exterieure au domaine maille.
  63 : 
  64 :  III/ 3eme syntaxe - Formulation EF
  65 : 
  66 :      Les fonctionnalites de cette option different sur deux points de la
  67 :      precedente. La fonction de paroi est integree a la premiere rangee
  68 :      d'element du maillage. La valeur de Yp (distance a la frontiere) est
  69 :      calculee automatiquement. La condition limite sur la QDM est une
  70 :      condition d'adherence (u=v=w=0). C'est FPU qui l'impose. La viscosite
  71 :      effective est modifiee sur la premiere rangee d'elements pour assurer
  72 :      la continuite de la contrainte avec le reste de l'ecoulement.
  73 :      Pour l'instant les valeurs de K et Epsilon ne sont pas calculees.
  74 :      Pour le reste cette option est identique a la precedente.
  75 : 
  76 : 
  77 :  [1] On Turbulent Flow Near a Wall. R.H. Van Driest.
  78 :      Journal of the Aeronautical Sciences (Nov 1956)
  79 : 
  80 :      Commentaires :
  81 :      ______________
  82 : 
  83 :  I/ 1ere syntaxe - Formulation EFM1 - 
  84 : 
  85 :   NU     Viscosite cinematique moleculaire (m**2/s) FLOTTANT
  86 :   YP     distance a la paroi   (m)                  FLOTTANT
  87 :          Cette distance doit etre telle que le premier point du maillage
  88 :          se situe dans la couche limite. Un Y+ compris entre 30 et 300
  89 :          est l'ideal. (Verification a posteriori par l'utilisateur)
  90 :   UET    Vitesse de frottement (m/s):
  91 :          en formulation EFM1                       CHPOINT SCAL CENTRE
  92 :                                                    MOT
  93 :   UN     Champ de vitesse   (m/s)                  CHPOINT VECT SOMMET
  94 :  <KN>    Energie turbulente                        CHPOINT SCAL SOMMET
  95 :  <EN>    Taux de dissipation de K                  CHPOINT SCAL SOMMET
  96 : 
  97 :      La vitesse de frottement UET, doit etre initialisee en formulation
  98 :      EFM1 a une valeur physiquement admissible sous peine de divergence
  99 :      de la solution. On peut s'aider de la formule ci-dessous.
 100 :      Choisir de preference une valeur d'epaisseur de couche limite YP
 101 :      telle que sa valeur adimentionnee Y+ soit comprise entre 30 et 300
 102 :      lors des calculs. (Y+=YP*UET/NU)
 103 : 
 104 :  II/ 2eme syntaxe - Formulation EF
 105 : 
 106 :   RO     Densite (Kg/m**3)                         FLOTTANT
 107 :                                                    CHPOINT SCAL SOMMET
 108 :                                                    MOT
 109 :   UN     Champ de vitesse   (m/s)                  CHPOINT VECT SOMMET
 110 :                                                    MOT
 111 :   MU     Viscosite dynamique moleculaire (Kg/m/s)  FLOTTANT
 112 :                                                    CHPOINT SCAL SOMMET
 113 :                                                    MOT
 114 :   UET    Vitesse de frottement (m/s):
 115 :          en formulation EF                         CHPOINT SCAL SOMMET
 116 :                                                    MOT
 117 :   YP     distance a la paroi   (m)                 FLOTTANT
 118 : 
 119 :  <KN>    Energie turbulente                        CHPOINT SCAL SOMMET
 120 :  <EN>    Taux de dissipation de K                  CHPOINT SCAL SOMMET
 121 : 
 122 :  III/ 3eme syntaxe - Formulation EF (Fonction de paroi integree au maillage)
 123 : 
 124 :   RO     Densite (Kg/m**3)                         FLOTTANT
 125 :                                                    CHPOINT SCAL SOMMET
 126 :                                                    MOT
 127 :   UN     Champ de vitesse   (m/s)                  CHPOINT VECT SOMMET
 128 :                                                    MOT
 129 :   MU     Viscosite dynamique moleculaire (Kg/m/s)  FLOTTANT
 130 :                                                    CHPOINT SCAL SOMMET
 131 :                                                    MOT
 132 :   UET    Vitesse de frottement (m/s):
 133 :          en formulation EF                         CHPOINT SCAL SOMMET
 134 :                                                    MOT
 135 :   $MD    Modele Navier-Stokes de la frontiere      MMODEL
 136 :   NUEFF  Viscosite effective qui sera modifiee     CHPOINT SCAL SOMMET
 137 :          par l'operateur.                          MOT
 138 : 
 139 : 
 140 : 
 141 : 
 142 :  Un coefficient de type MOT indique que l'operateur va chercher le
 143 :  coefficient dans la table INCO a l'indice MOT.
 144 : 
 145 : 
 146 :      OPTION
 147 :      ______
 148 : 
 149 :      Formulation Element Fini                   OPTION EFM1 (default)
 150 :                                                 OPTION EF
 151 :