Ecoulement autour d'un cylindre en rotation
Nom de fichier: /home/buffat/COURS/FEMLAB/CHAP1/cylrot.fl
Modes d'application et modules utilisés dans ce modèle:
| Propriété | Valeur |
| Nom du modèle | Ecoulement autour d'un cylindre en rotation |
| Auteur | M. BUFFAT |
| Entreprise | UFR de mecanique, UCB Lyon |
| Service | |
| Référence | |
| Date Enregistrée | 7 mars 2005 10:18:48 |
| Date de création | 6 mars 2005 11:21:46 |
| version FEMLAB | FEMLAB 3.1.0.157 |
Calcul de la fonction de courant pour un ecoulement potentiel autour d'un cylindre en rotation
Calcul de la force de portance
| Nom | Expression | Valeur |
| U0 | 1 | 1 |
| R | 0.1 | 0.1 |
| Psi1 | 0.4 | 0.4 |
Nombres de géométries: 1
Dimensions d’espace: 2D
| Nom | Expression |
| Ue | U0*y*(1-R*R/(x*x+y*y)) |
| Pr | -1/2*(ux*ux+uy*uy) |
| Nombre de degrés de liberté | 2240 |
| Nombre d’éléments de frontière | 72 |
| Nombre d’éléments | 1084 |
| Qualité minimale des éléments | 0.5959 |
Type de mode d’application: Équation de Laplace
Nom du mode d’application: lp
| Propriété | Valeur |
| Type d’élément par défaut | Lagrange - Quadratique |
| Contraintes faibles | Off |
Variables dépendantes: u
Variables indépendantes: x, y, z
Fonctions de forme: shlag(2,'u')
IFrontières internes non actives
| Point | 1-8 |
| terme faible (weak) | 0 |
| Terme dweak (dweak) | 0 |
| terme contraint (constr) | '0' |
| weakconstr | 1 |
| Fonctions de forme (wcshape) | [] |
| Condition initiale (wcinit) | 0 |
| Frontière | 1-2, 5, 8 | 3-4, 6-7 |
| Type | Conditions aux limites de Dirichlet | Conditions aux limites de Dirichlet |
| terme faible (weak) | 0 | 0 |
| Terme dweak (dweak) | 0 | 0 |
| terme contraint (constr) | 0 | 0 |
| (q) | 0 | 0 |
| (h) | 1 | 1 |
| (g) | nx | 0 |
| (r) | U0*y | Psi1 |
| weakconstr | 1 | 1 |
| Fonctions de forme (wcshape) | [] | [] |
| Ordre d’integration (wcgporder) | 2 | 2 |
| Condition initiale (wcinit) | 0 | 0 |
| Sous-Domaine | 1 |
| Fonctions de forme (shape) | shlag(2,'u') |
| Ordre d’integration (gporder) | 4 |
| Ordre de contrainte (cporder) | 2 |
| terme faible (weak) | 0 |
| Terme dweak (dweak) | 0 |
| terme contraint (constr) | 0 |
| weakconstr | 1 |
| Condition initiale de sous-domaine | 1 |
| u | 0 |
Résoudre à l’aide d’un script: off
| Solveur selectionné automatiquement | on |
| Solveur | Stationnaire linéaire |
| Forme de Résolution | coefficient |
| Symmetric | off |
| Adaption | off |
Type de solveur: Solveur de systèmes linéaires
| Paramètre | Valeur |
| Valeur seuil du pivot | 0.1 |
| Facteur d’allocation de la memoire | 0.7 |
| Paramètre | Valeur |
| Méthode de prise en compte des contraintes | Eliminate |
| Null-space function | Auto |
| Taille des blocs d’assemblage | 5000 |
| Utiliser l’Hermitien transposé | On |
| Utiliser les fonctions complexes à arguments réels | Off |
| Type de scaling | Auto |
| Scaling manuel | |
| Equilibrage des colonnes | On |
| Contrôle manuel du réassemblage | Off |
| Charge constante | On |
| Contraintes constantes | On |
| Masse constante | On |
| Jacobien constant | On |
| Jacobien des contraintes constant | On |
| Nom | Description | Expression |
| absux_lp | |grad(u)| | sqrt(ux^2+uy^2) |
| abscu1x_lp | |c*grad(u)| | sqrt(cu1x^2+cu1y^2) |