Accueil | Cast3M

   1 : $$$$ SPAL     NOTICE  CHAT      11/09/12    21:18:10     7124           
   2 :                                              DATE     11/09/12
   3 : 
   4 : Procedure  SPAL                                  Voir aussi : KEPSILON
   5 : 
   6 :   SYNTAXE (cf. operateur EQEX)
   7 :   ____________________________
   8 : 
   9 :         'ZONE' $MD 'OPER' 'SPAL' 'RHO' 'UN' 'MU' 'DT'
  10 :                                  ('PERIODIC' GEOM1 GEOM2)
  11 :                                  'INCO' 'NU0'
  12 : 
  13 : 
  14 : 
  15 : 
  16 :    Objet :
  17 :    _______
  18 : 
  19 :    Calcule le champ de viscosite dynamique turbulente grace au modele
  20 :    de Spalart-Allmaras.
  21 : 
  22 : 
  23 :    Commentaires :
  24 :    ______________
  25 : 
  26 :    1) LES PARAMÈTRES REQUIS sont:
  27 : 
  28 :       RHO*[MOT|FLOTTANT|CHPOINT] : Masse volumique            (kg/m3)
  29 :       UN *[MOT|CHPOINT]          : Vitesse d'advection        (m/s)
  30 :       MU *[MOT|FLOTTANT|CHPOINT] : Viscosite moleculaire dyn. (Pa.s)
  31 :       DT *[MOT|FLOTTANT]         : Duree du pas de temps      (s)
  32 :       NU0*[MOT]                  : Nom attribue a la viscosite modifiee
  33 : 
  34 :       Un objet de type MOT indique que l'on va chercher la valeur
  35 :       dans la table 'INCO'.
  36 : 
  37 :    /!\ ATTENTION: Ce modele necessite aussi la donnee de la distance
  38 :                   a la paroi dans (RV.'PAROIS'.'DIST') !!
  39 : 
  40 : 
  41 : 
  42 :    2) Le champ de viscosite effective (moleculaire + turbulente) est
  43 :       renvoye dans la table 'INCO' a l'indice 'MUFN' par defaut, mais
  44 :       l'utilisateur peut definir ce nom lui-meme (cf. definition des
  45 :       parametres avances, remarque 4)
  46 : 
  47 :    /!\ ATTENTION: Les conditions aux limites de Dirichlet ainsi que
  48 :                   les conditions initiales devront porter sur 'NU0'
  49 :                   et non pas sur 'MUFN' !!
  50 : 
  51 : 
  52 : 
  53 :    3) LE PARAMÈTRE OPTIONNEL 'PERIODIC' permet d'imposer des
  54 :       conditions de periodicite sur 'MUFN' entre les maillages GEOM1
  55 :       et GEOM2.
  56 : 
  57 : 
  58 : 
  59 :    4) LES PARAMÈTRES AVANCÉS du modele peuvent etre personnalises en
  60 :       ajoutant une table nommee 'SPALART_ALLMARAS' dans RV:
  61 : 
  62 :    +----------------------------------------------------------------+
  63 :    | TABLE PRINCIPALE 'SPALART_ALLMARAS'                            |
  64 :    +-----------+-----------+----------------------------------------+
  65 :    | Indice    | Valeur    | Description                            |
  66 :    +-----------+-----------+----------------------------------------+
  67 :    | 'KVERS'   | [MOT-cle] | Variante du modele a utiliser:         |
  68 :    |           | 'ORIG'    | - Modele original de base (par defaut) |
  69 :    |           | 'TRIP'    | - Modele original avec ft1 et ft2      |   <= À FAI
  70 :    |           | 'SALSA'   | - Modele modifie par Rung et al.       |   <= À FAI
  71 :    |           |           |                                        |
  72 :    | 'NOMMUF'  | MOT       | Nom de l'inconnue de viscosite totale  |
  73 :    |           |           |                   ('MUFN' par defaut)  |
  74 :    |           |           |                                        |
  75 :    | 'KCONST'  | TABLE     | Constantes du modele                   |
  76 :    |           |           |                                        |
  77 :    | 'KTGRAD'  | [MOT-cle] | Mesure scalaire du tenseur gradient:   |
  78 :    |           | 'TOROT'   | - Taux de rotation (Par defaut)        |
  79 :    |           | 'TODEF'   | - Taux de deformation                  |
  80 :    |           | 'COMPL'   | - Norme euclidienne du tenseur complet |
  81 :    |           | 'MIXTE'   | - Expression de Dacles-Mariani et al.  |
  82 :    |           |           |                                        |
  83 :    | 'KMUFN'   | [MOT-cle] | Instant auquel est renvoye 'MUFN'      |
  84 :    |           | 'APRES'   | - Fin du pas de temps (Par defaut)     |
  85 :    |           | 'AVANT'   | - Debut du pas de temps                |
  86 :    |           | 'DEMI'    | - Apres le demi pas de temps ('ALGO1') |
  87 :    |           |           |                                        |
  88 :    | 'KSRC'    | [MOT-cle] | Algo. utilise pour les termes sources  |
  89 :    |           | 'ALGO1'   | - Methode de Newton (Par defaut)       |
  90 :    |           | 'ALGO2'   | - Separation S+/S-                     |
  91 :    |           |           |                                        |
  92 :    | 'NEWTON'  | TABLE     | Parametres de l'algorithme de Newton   |
  93 :    |           |           | si le parametre 'KSRC' vaut 'ALGO1'    |
  94 :    |           |           |                                        |
  95 :    | 'METHINV' | TABLE     | Options de la methode d'inversion.     |
  96 :    |           |           | Par defaut, ce sont celles definies    |
  97 :    |           |           | pour le probleme global dans RV        |
  98 :    |           |           |                                        |
  99 :    | 'KOPT2'   | TABLE     | Options de discretisation temporelle.  |
 100 :    |           |           | Par defaut, les parametres passes par  |
 101 :    |           |           | EQEX 'OPTI' sont appliques a la fois a |
 102 :    |           |           | l'operateur de convection/diffusion    |
 103 :    |           |           | TSCA et a la derivee temporelle DFDT   |
 104 :    |           |           | (sauf que DFDT est 'CENTREE')          |
 105 :    |           |           |                                        |
 106 :    | 'VERROU'  | TABLE     | Etat des verrous numeriques            |
 107 :    |           |           |                                        |
 108 :    | 'DUMP'    | LOGIQUE   | Sauver les variables internes au modele|
 109 :    |           |           | SA dans la table 'INCO' sous l'indice  |
 110 :    |           |           | 'SPAL'? (Par defaut: FAUX)             |
 111 :    +-----------+-----------+----------------------------------------+
 112 : 
 113 :    +----------------------------------------------------------------+
 114 :    | SOUS-TABLE 'KCONST'                                            |
 115 :    +-----------+-----------+----------------------------------------+
 116 :    | Indice    | Valeur    | Description                   | Defaut |
 117 :    +-----------+-----------+-------------------------------+--------+
 118 :    | 'SIGMA'   | FLOTTANT  | Nombre de Prandtl turbulent   | 2/3    |
 119 :    | 'CB1'     | FLOTTANT  | Taux de production turbulente | 0.1355 |
 120 :    | 'CB2'     | FLOTTANT  | Diffusion non conservative    | 0.622  |
 121 :    | 'KAPPA'   | FLOTTANT  | Constante de Von Karman       | 0.41   |
 122 :    | 'CV1'     | FLOTTANT  | Épai. ss-couche visq. (Bas-Re)| 7.1    |
 123 :    | 'CW1'     | FLOTTANT  | Equilibre Prod/Dest zone log. | 3.2391 |
 124 :    | 'CW2'     | FLOTTANT  | Controle du coef. frottement  | 0.3    |
 125 :    | 'CW3'     | FLOTTANT  | Borne sup. de fw (environ)    | 2.     |
 126 :    +-----------+-----------+----------------------------------------+
 127 : 
 128 :    +----------------------------------------------------------------+
 129 :    | SOUS-TABLE 'NEWTON'                                            |
 130 :    +-----------+-----------+----------------------------------------+
 131 :    | Indice    | Valeur    | Description                   | Defaut |
 132 :    +-----------+-----------+-------------------------------+--------+
 133 :    | 'CRIT'    | FLOTTANT  | Critere d'arret en norme inf. | 1.E-10 |
 134 :    | 'IMAX'    | ENTIER    | Nombre max. d'iterations      | 10     |
 135 :    | 'OMEGA'   | FLOTTANT  | Facteur de relaxation         | 1.     |
 136 :    +-----------+-----------+----------------------------------------+
 137 : 
 138 :    +----------------------------------------------------------------+
 139 :    | SOUS-TABLE 'VERROU'                                            |
 140 :    +-----------+-----------+----------------------------------------+
 141 :    | Indice    | Valeur    | Description                   | Defaut |
 142 :    +-----------+-----------+-------------------------------+--------+
 143 :    | 'POSITIF' | LOGIQUE   | Maintenir NU0 et S0 positifs  | VRAI   |
 144 :    | 'DURBIN'  | LOGIQUE   | Realisabilite selon Durbin    | FAUX   |   <= À FAI
 145 :    +-----------+-----------+----------------------------------------+
 146 : 
 147 : 
 148 : 
 149 :    Complements : PRÉCISIONS SUR LES ALGORITHMES DISPONIBLES
 150 :    _____________
 151 : 
 152 :    Les termes sources sont tous non lineaires, et peuvent devenir
 153 :    negatifs (du moins dans sa version originale). Deux algorithmes ont
 154 :    ete mis en place pour ameliorer la robustesse:
 155 : 
 156 :    'ALGO1' => Le traitement des termes sources est effectue avant
 157 :               l'etape de convection/diffusion (methode de splitting).
 158 :               On avance jusqu'a un pas de temps intermediaire grace
 159 :               a une methode de Newton prenant en compte uniquement
 160 :               les termes sources. Le demi-pas de temps restant
 161 :               concerne alors la resolution de l'EDO de transport,
 162 :               sans termes sources.
 163 : 
 164 :    'ALGO2' => La partie negative du terme source est linearisee afin
 165 :               d'etre implicitee, tandis que la partie positive est
 166 :               laissee au second membre (explicite). On resout alors
 167 :               l'equation de transport complete, avec sources.
 168 : 
 169 : 
 170 :