Accueil | Cast3M

   1 : $$$$ KEPSILON NOTICE  CHAT      11/09/12    21:16:43     7124           
   2 :                                              DATE     11/09/12
   3 : 
   4 :    Procedure KEPSILON                       Voir aussi :
   5 :     ------------------   
   6 :     SYNTAXE ( EQEX ) : Cf operateur EQEX
   7 :     _________________
   8 : 
   9 :     'OPER' 'KEPSILON' RO UN MU DT (RGB TN) 'INCO' 'KN' 'EN'
  10 :                                                                   
  11 : 
  12 :          
  13 :     OBJET :
  14 :     -----
  15 :          
  16 :  Cette procedure calcule la viscosite effective (tourbillonnaire +
  17 :  moleculaire) obtenue par la resolution en transitoire (un pas de temps)
  18 :  d'un modele K-Epsilon. Le resultat est place dans la table 'INCO' a
  19 :  l'indice 'MUF' pour la viscosite dynamique effective (kg/m/s).
  20 :  Pour la version standard le desequilibre (Production/Dissipation) est
  21 :  limite a 10 (Limitation proposee par Menter).
  22 :  Precautions: L'algorithme ne converge pas (sans qu'il diverge) lorsque 
  23 :  l'allongement des mailles du maillage est superieur a 20!!
  24 :          
  25 :     Commentaires
  26 :     ____________ 
  27 :                                                                                 
  28 :     RO     Densite
  29 :            FLOTTANT ou MOT
  30 :                                                                                 
  31 :     UN     Champ de vitesse
  32 :            CHPOINT VECT SOMMET ou MOT
  33 : 
  34 :     MU     Viscosite dynamique (Kg/m/s)
  35 :            FLOTTANT ou MOT
  36 : 
  37 :     DT     Pas de temps
  38 :            FLOTTANT ou MOT
  39 :                                                                                 
  40 :     RGB    Coefficient du terme de flottabilite
  41 :            VECTEUR ou MOT
  42 :                                                                                 
  43 :     TN     Champ de temperature
  44 :            CHPOINT SCAL SOMMET ou MOT
  45 : 
  46 :  Un coefficient de type MOT indique que l'operateur va chercher le
  47 :  coefficient dans la table INCO a l'indice donne.
  48 : 
  49 :  Les options (parametres) de cette procedure doivent se trouver dans un
  50 :  LISTMOTS a l'entree 'ALGO_KEPSILON' de la procedure RV.
  51 :  ex : RV.'ALGO_KEPSILON'= MOTS 'Bw' 'Cnu';
  52 :  En l'absence de cette entree les valeurs par defaut sont prises.
  53 : 
  54 :  La liste des Options/Parametres est : IMPR,RNG,Filtre,Bw,Cnu,Nut,Fi,
  55 :  M2M,CSTE,Ret,KL,KLbr,Chien,Sharma,Jones,Lam
  56 :  Les options par defaut sont : Nut et le modele k - epsilon standard.
  57 : 
  58 :  IMPR: permet d'afficher les options utilisees.
  59 : 
  60 :  RNG : modele RNG k - epsilon.
  61 : 
  62 :  Filtre : Modele K-epsilon filtre. L'echelle de longueur est
  63 :          filtree a une valeur precisee dans RV.'INCO'.'Echl'. Ce modele
  64 :          permet de mieux capter des instationarites ou des instabilites
  65 :          a grande echelle (de taille superieure a celle du filtre). On
  66 :          peut prendre comme echelle de longueur la taille des elements
  67 :          du maillage.
  68 : 
  69 :  Bw  : declenche une condition de realisabilite sur les contraintes de
  70 :        cisaillement maximum. On verifie que (u'v')/k < 0.3   (Bradshaw)
  71 :        DEFAUT = FAUX.
  72 : 
  73 :  Cnu : La 'constante' Cnu est reliee au rapport Ksi=(Nut P)/epsilon :
  74 :        Cnu = F(1/Ksi) (Voir Rodi)
  75 :        Pour Ksi = 1 Cnu=0.09.
  76 :        DEFAUT = FAUX.
  77 : 
  78 :  Nut : Les variables de resolution intermediaires sont Teta et Nut
  79 :        DEFAUT = VRAI.
  80 : 
  81 :  Fi  : Les variables de resolution intermediaires sont Teta et Fi
  82 :        DEFAUT = FAUX.
  83 : 
  84 :  M2M : Constantes de Mohammadi et Medic.
  85 :        DEFAUT = FAUX.
  86 : 
  87 :  CSTE: Les constantes du model sont lues dans la table INCO aux entrees
  88 :        'cnu' 'c2' 'sgk' 'sge' La constante c1 en est deduite
  89 :        DEFAUT = FAUX.
  90 : 
  91 :  KL  : Modele K-L. L'echelle de longueur doit etre specifiee dans la
  92 :        table 'INCO' (entree 'INCO'.'Echl'). Numeriquement ce modele est
  93 :        obtenu en remplaçant l'equation (EDP) sur epsilon par
  94 :        epsilon = k**1.5 / L. Il est necessaire de renseigner le modele
  95 :        comme si on resolvait le modele K-epsilon. La condition limite
  96 :        sur epsilon doit verifier l'equation ci-dessus
  97 : 
  98 :  KLbr: Modele K-L bas Reynolds. C'est le modele de Wolfshtein (1967) et
  99 :        Yap (1987). A la place de l'echelle de longueur il faut donner la
 100 :        distance a la paroi: indice 'dparoi' (CHPOINT SCAL SOMMET) dans
 101 :        la table INCO. Les conditions limites sont U=0,K=0 et Epsilon=0 a
 102 :        la paroi. La premiere maille doit se trouver a un y+ < 1.
 103 :        Voir le jeux de donnee  canalKLbr.dgibi
 104 : 
 105 :  Chien: Modele k-epsilon Bas-Reynolds de Chien. Ce modele
 106 :        necessite la donnee de la distance a la paroi:
 107 :        entree 'INCO'.'dparoi' (meme support que k ou epsilon) et le
 108 :        calcul de y+ (entree 'INCO' 'yplus'). Cette derniere entree
 109 :        doit etre recalculee a chaque pas de temps via une procedure.
 110 :        Les conditions limites et les contraintes de maillage pres de la
 111 :        paroi sont les memes que precedemment.
 112 :        Voir le jeux de donnee  canal-Chien.dgibi
 113 : 
 114 :  Sharma: Modele k-epsilon Bas-Reynolds de Launder Sharma. Ce modele ne
 115 :        necessite aucune donnee supplementaire ce qui presente un
 116 :        avantage certain.
 117 :        Les conditions limites et les contraintes de maillage pres de la
 118 :        paroi sont les memes que precedemment.
 119 :        Voir le jeux de donnee  canal-Sharma.dgibi
 120 : 
 121 :  Jones: Modele k-epsilon Bas-Reynolds de Jones-Launder Sharma. Ce modele 
 122 :        est quasiment identique au precedent aux valeurs des constantes
 123 :        pres.
 124 : 
 125 :   Lam:  Modele k-epsilon Lam et Bremhorst (en test)
 126 : 
 127 :  Liste des CHPOINTs (SCAL SOMMET) crees dans la table 'INCO'
 128 :        TKTE      teta=k/epsilon
 129 :        NUTI      valeur intermediaire de NUT
 130 :        FI        inconnue Fi
 131 :        PRODT     Prodution turbulente / NUT = (grad U + grad^t U)grad U
 132 :        TKTI      teta=k/epsilon intermediaire
 133 :        Ksi       facteur de desequilibre : nut PRODT / epsilon
 134 :        MUF       viscosite dynamique effective (tourbillonnaire+moleculaire)
 135 :