TP5 Rouen

Université de Rouen Calcul Scientifique
L3-Mathématiques I. DANAILA
2011–2012 ionut.danaila@univ-rouen.fr
À la fin de cette séance vous fichiers seront copiés sur la clé USB du professeur. Toute
connexion à la messagerie électronique est strictement interdite pendant l’examen.
TP 5 : Contrôle continu
Exercice 1 Calcul d’intégrales 2D (programme CC_ex1.edp)
Calculer numériquement avec FreeFem++ l’aire d’une ellipse de demi-axes a = 2.5 et b = 4.
Exercice 2 Écoulement potentiel autour d’un cylindre (programme CC_ex2.edp)
Considérons l’écoulement plan, stationnaire d’un fluide parfait sans viscosité, autour d’un cylindre.
Γ3
La vitesse d’une particule située en un point (x, y) est dérivée à H
partir de la fonction potentiel ϕ
 
∂ϕ Γ4 r Γ2
→
−  ∂x 
v = ∇ϕ =  − ∆ϕ = 0.
 ∂ϕ  , avec ϕ solution de : 0 L

∂y
Le domaine de calcul sera le carré Ω = [0, L] × [−H, H]. Ω

-H
Le cylindre de rayon r est placé au milieu du domaine. Γ1
Nous supposons que les frontières du domaines ont suffisamment éloignées du cylindre pour retrouver
un courant uniforme, de vitesse constante. Les conditions aux limites seront donc :
sur Γ1 ∪ Γ2 ∪ Γ3 ∪ Γ4 : ~v = v0 ~ex , avec ~ex le vecteur unitaire suivant x.
sur la frontière du cylindre Γ5 : ~v · ~n = 0, la vitesse normale est nulle.
• Exprimer les conditions aux limites sur chaque frontière Γ1 , Γ2 , Γ3 , Γ4 par rapport à ∂φ/∂n.
Attention, la normale ~n est dirigée toujours vers l’extérieur du domaine !
Rappelons la définition de la dérivée suivant la normale : ∂ϕ/∂n = ∇ϕ · ~n.
• Écrire la formulation variationnelle du problème −∆ϕ = 0, en tenant compte des conditions aux
limites établies précédemment.
• Résoudre le problème avec FreeFem++. Tracer la solution ϕ (graphique à envoyer).
Données numériques : L = 8, H = 4, r = 1, v0 = 1.
• Calculer les deux composantes de la vitesse et tracer les vecteurs vitesse (graphique à envoyer).
Indication : on peut utiliser la syntaxe FreeFem++ pour calculer les dérivées u=dx(phi) ; v =dy(phi);
Pour représenter les vecteurs, on utilise plot([u,v]);
1
Exercice 3 Modélisation d’une ailette de refroidissement
Considérons une ailette de refroidissement de forme rectangulaire Ω = [0, L] × [−H, H]. La troisième
dimension W est considérée beaucoup plus grande pour que le problème soit bidimensionnel.
À l’extrémité Γ4 l’ailette est en contact avec une source de cha- Γ3

H
leur de température Tc constante. La température à l’équilibre est
la solution de l’EDP : Γ4 Γ2
−∆θ = 0, pour (x, y) ∈ Ω. L
La variable θ = (T − Ta )/Ta est la température adimensionnée par -H Ω

rapport à l’ambiance Ta (attention, températures en kelvins). Γ1
Les conditions aux limites sont imposées comme suit :

(
sur Γ4 θ = θc température imposée (cond. de Dirichlet)
sur Γ1 ∪ Γ2 ∪ Γ3 ∂θ/∂n + αθ = 0 pertes par convection (cond. de Fourier)
1. Ecrire la formulation variationnelle correspondante.

2. Résoudre le problème avec FreeFem++ (programme CC_ex3A.edp). Tracer la distribution de
température T (en o C) (iso-contours et en 3D). Commenter le résultat.
Données numériques : L = 1, H = 0.25, Ta = 20o C, Tc = 46o C, α = 10.
Question bonus (programme CC_ex3B.edp) Γ3

On suppose maintenant qu’à l’intérieur du l’ailette se trouve H
une inclusion (domaine circulaire Ω1 dans la figure) consti-
Γ4 r Γ2
tuée d’un matériau de conductivité k différente. L’équation à
résoudre devient : Ω1 L
Ω
(
k1 = 1 sur Ω\Ω1 -H
−∇.(k∇u) = 0, aveck(x, y) = Γ1
k2 = 10 sur Ω1
On positionnera le domaine Ω1 (disque de rayon r = H ∗ 0.8) au milieu de l’ailette.

R R R
• Écrire la nouvelle formulation variationnelle. Indication : utiliser le fait que Ω = Ω\Ω1 + Ω1 et que
la fonction k est constante sur chaque sous-domaine.
• Résoudre ce problème et tracer la distribution de température T (en o C).
Programmes à préparer à la fin de l’examen :

CC_ex1.edp
CC_ex2.edp et deux images (fonction potentiel, vecteurs vitesse)
CC_ex3A.edp et la visualisation de la température en o C
CC_ex3B.edp et la visualisation de la température en o C

TP5 Rouen

Transféré par

Droits d'auteur :

Formats disponibles

TP5 Rouen

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

TP5 Rouen

Transféré par

Droits d'auteur :

Formats disponibles

Université de Rouen Calcul Scientifique

Exercice 1 Calcul d’intégrales 2D (programme CC_ex1.edp)

Calculer numériquement avec FreeFem++ l’aire d’une ellipse de demi-axes a = 2.5 et b = 4.

Exercice 2 Écoulement potentiel autour d’un cylindre (programme CC_ex2.edp)

Le domaine de calcul sera le carré Ω = [0, L] × [−H, H]. Ω

sur Γ1 ∪ Γ2 ∪ Γ3 ∪ Γ4 : ~v = v0 ~ex , avec ~ex le vecteur unitaire suivant x.

sur la frontière du cylindre Γ5 : ~v · ~n = 0, la vitesse normale est nulle.

• Résoudre le problème avec FreeFem++. Tracer la solution ϕ (graphique à envoyer).

À l’extrémité Γ4 l’ailette est en contact avec une source de cha- Γ3

−∆θ = 0, pour (x, y) ∈ Ω. L

La variable θ = (T − Ta )/Ta est la température adimensionnée par -H Ω

Les conditions aux limites sont imposées comme suit :

1. Ecrire la formulation variationnelle correspondante.

Question bonus (programme CC_ex3B.edp) Γ3

On positionnera le domaine Ω1 (disque de rayon r = H ∗ 0.8) au milieu de l’ailette.

Programmes à préparer à la fin de l’examen :

Vous aimerez peut-être aussi

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.