Afficher cette notice en
$$$$ SPAL NOTICE CHAT 11/09/12 21:18:10 7124 DATE 11/09/12 Procedure SPAL Voir aussi : KEPSILON SYNTAXE (cf. operateur EQEX) ____________________________ 'ZONE' $MD 'OPER' 'SPAL' 'RHO' 'UN' 'MU' 'DT' ('PERIODIC' GEOM1 GEOM2) 'INCO' 'NU0' Objet : _______ Calcule le champ de viscosite dynamique turbulente grace au modele de Spalart-Allmaras. Commentaires : ______________ 1) LES PARAMÈTRES REQUIS sont: RHO*[MOT|FLOTTANT|CHPOINT] : Masse volumique (kg/m3) UN *[MOT|CHPOINT] : Vitesse d'advection (m/s) MU *[MOT|FLOTTANT|CHPOINT] : Viscosite moleculaire dyn. (Pa.s) DT *[MOT|FLOTTANT] : Duree du pas de temps (s) NU0*[MOT] : Nom attribue a la viscosite modifiee Un objet de type MOT indique que l'on va chercher la valeur dans la table 'INCO'. /!\ ATTENTION: Ce modele necessite aussi la donnee de la distance a la paroi dans (RV.'PAROIS'.'DIST') !! 2) Le champ de viscosite effective (moleculaire + turbulente) est renvoye dans la table 'INCO' a l'indice 'MUFN' par defaut, mais l'utilisateur peut definir ce nom lui-meme (cf. definition des parametres avances, remarque 4) /!\ ATTENTION: Les conditions aux limites de Dirichlet ainsi que les conditions initiales devront porter sur 'NU0' et non pas sur 'MUFN' !! 3) LE PARAMÈTRE OPTIONNEL 'PERIODIC' permet d'imposer des conditions de periodicite sur 'MUFN' entre les maillages GEOM1 et GEOM2. 4) LES PARAMÈTRES AVANCÉS du modele peuvent etre personnalises en ajoutant une table nommee 'SPALART_ALLMARAS' dans RV: +----------------------------------------------------------------+ | TABLE PRINCIPALE 'SPALART_ALLMARAS' | +-----------+-----------+----------------------------------------+ | Indice | Valeur | Description | +-----------+-----------+----------------------------------------+ | 'KVERS' | [MOT-cle] | Variante du modele a utiliser: | | | 'ORIG' | - Modele original de base (par defaut) | | | 'TRIP' | - Modele original avec ft1 et ft2 | <= À FAI | | 'SALSA' | - Modele modifie par Rung et al. | <= À FAI | | | | | 'NOMMUF' | MOT | Nom de l'inconnue de viscosite totale | | | | ('MUFN' par defaut) | | | | | | 'KCONST' | TABLE | Constantes du modele | | | | | | 'KTGRAD' | [MOT-cle] | Mesure scalaire du tenseur gradient: | | | 'TOROT' | - Taux de rotation (Par defaut) | | | 'TODEF' | - Taux de deformation | | | 'COMPL' | - Norme euclidienne du tenseur complet | | | 'MIXTE' | - Expression de Dacles-Mariani et al. | | | | | | 'KMUFN' | [MOT-cle] | Instant auquel est renvoye 'MUFN' | | | 'APRES' | - Fin du pas de temps (Par defaut) | | | 'AVANT' | - Debut du pas de temps | | | 'DEMI' | - Apres le demi pas de temps ('ALGO1') | | | | | | 'KSRC' | [MOT-cle] | Algo. utilise pour les termes sources | | | 'ALGO1' | - Methode de Newton (Par defaut) | | | 'ALGO2' | - Separation S+/S- | | | | | | 'NEWTON' | TABLE | Parametres de l'algorithme de Newton | | | | si le parametre 'KSRC' vaut 'ALGO1' | | | | | | 'METHINV' | TABLE | Options de la methode d'inversion. | | | | Par defaut, ce sont celles definies | | | | pour le probleme global dans RV | | | | | | 'KOPT2' | TABLE | Options de discretisation temporelle. | | | | Par defaut, les parametres passes par | | | | EQEX 'OPTI' sont appliques a la fois a | | | | l'operateur de convection/diffusion | | | | TSCA et a la derivee temporelle DFDT | | | | (sauf que DFDT est 'CENTREE') | | | | | | 'VERROU' | TABLE | Etat des verrous numeriques | | | | | | 'DUMP' | LOGIQUE | Sauver les variables internes au modele| | | | SA dans la table 'INCO' sous l'indice | | | | 'SPAL'? (Par defaut: FAUX) | +-----------+-----------+----------------------------------------+ +----------------------------------------------------------------+ | SOUS-TABLE 'KCONST' | +-----------+-----------+----------------------------------------+ | Indice | Valeur | Description | Defaut | +-----------+-----------+-------------------------------+--------+ | 'SIGMA' | FLOTTANT | Nombre de Prandtl turbulent | 2/3 | | 'CB1' | FLOTTANT | Taux de production turbulente | 0.1355 | | 'CB2' | FLOTTANT | Diffusion non conservative | 0.622 | | 'KAPPA' | FLOTTANT | Constante de Von Karman | 0.41 | | 'CV1' | FLOTTANT | Épai. ss-couche visq. (Bas-Re)| 7.1 | | 'CW1' | FLOTTANT | Equilibre Prod/Dest zone log. | 3.2391 | | 'CW2' | FLOTTANT | Controle du coef. frottement | 0.3 | | 'CW3' | FLOTTANT | Borne sup. de fw (environ) | 2. | +-----------+-----------+----------------------------------------+ +----------------------------------------------------------------+ | SOUS-TABLE 'NEWTON' | +-----------+-----------+----------------------------------------+ | Indice | Valeur | Description | Defaut | +-----------+-----------+-------------------------------+--------+ | 'CRIT' | FLOTTANT | Critere d'arret en norme inf. | 1.E-10 | | 'IMAX' | ENTIER | Nombre max. d'iterations | 10 | | 'OMEGA' | FLOTTANT | Facteur de relaxation | 1. | +-----------+-----------+----------------------------------------+ +----------------------------------------------------------------+ | SOUS-TABLE 'VERROU' | +-----------+-----------+----------------------------------------+ | Indice | Valeur | Description | Defaut | +-----------+-----------+-------------------------------+--------+ | 'POSITIF' | LOGIQUE | Maintenir NU0 et S0 positifs | VRAI | | 'DURBIN' | LOGIQUE | Realisabilite selon Durbin | FAUX | <= À FAI +-----------+-----------+----------------------------------------+ Complements : PRÉCISIONS SUR LES ALGORITHMES DISPONIBLES _____________ Les termes sources sont tous non lineaires, et peuvent devenir negatifs (du moins dans sa version originale). Deux algorithmes ont ete mis en place pour ameliorer la robustesse: 'ALGO1' => Le traitement des termes sources est effectue avant l'etape de convection/diffusion (methode de splitting). On avance jusqu'a un pas de temps intermediaire grace a une methode de Newton prenant en compte uniquement les termes sources. Le demi-pas de temps restant concerne alors la resolution de l'EDO de transport, sans termes sources. 'ALGO2' => La partie negative du terme source est linearisee afin d'etre implicitee, tandis que la partie positive est laissee au second membre (explicite). On resout alors l'equation de transport complete, avec sources.
© Cast3M 2003 - Tous droits réservés.
Mentions légales