Discussion :

[Méliodas]

Bonjour,

En mécanique des fluides, j'ai calculé une fonction de courant qui est la suivante: Psy=Vsin(theta)r -(Vsin(theta)a²)/r.

Elle modélise un flux de courant d'une vitesse V sur un cylindre de rayon a en coordonnée polaire.

Je souhaiterai représenter sous matlab cette fonction ainsi que ces isovaleurs.

Pour l'instant mon code est le suivant:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
x=linspace(-10,10,200);
y=linspace(-10,10,200);
a=2.5;
V=100;
[X,Y] = meshgrid(x,y);
[theta,R] = cart2pol(X,Y)
psy= V*sin(theta)*R - a^2./R*sin(theta);
contour(X,Y,psy,30)

Cependant, ma figure me donne un résultat erroné. Pouvais-vous m’indiquer mes éventuelles erreurs s'il vous plait ?

Merci.

[le fab]

Salut

Cependant, ma figure me donne un résultat erroné.

erroné comment ?

[phryte]

Bonjour,

A la ligne 7 selon le texte il faudrait :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

-V*a.^2/R*sin(theta)

[Méliodas]

Merci, pour cette petite erreur.

Est-il possible de modéliser en plus de cet écoulement, le cylindre en question sous matlab ?

Techniquement, je ne devrais pas avoir de lignes de courants passant par le centre de mon graphique car elle ne traverse pas le cylindre dans la réalité.

[phryte]

Bonjour,
Une petite avancée, il doit y avoir encore une erreur :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
clear
x=linspace(-10,10,10);
y=linspace(-10,10,10);
a=2.5;
V=100;
[X,Y] = meshgrid(x,y);
[theta,R] = cart2pol(X,Y);
psy= V*sin(theta).*R.^2 - V*a^2./R.^2*sin(theta);
contour(X,Y,psy,10)
hold on
viscircles([0 0],a);
axis square
grid

[Méliodas]

Super ! Merci beaucoup.