Discussion :

[flotty]

Bonjour,

Je cherche à établir les courbes des fonctions de courant et potentielle à partir de la vitesse complexe W(z) avec z un nombre complexe comme donnée en pièce jointe (écoulement autour d'un cylindre).

Je connais le résultat à obtenir (troisième cas en partant de la gauche) mais n'arrive pas à y parvenir.

Voici le code que j'ai écrit pour tenter d'obtenir les courbes:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
 
function [y,z,r] = cylinder(Vinf,Alpha,M)
hold on;
Vinf=3;
Alpha=2;
M=5;
Rz=sqrt(M/(2*pi*Vinf));
Gamma=-2*Rz;
z=complex(2,2);
W=Vinf*(z*(cos(Alpha)-i*sin(Alpha))+(Rz.^2)*(cos(Alpha)+i*sin(Alpha))/z)-((i*Gamma)/(2*pi))*log10(z);
for Alpha=0:pi/10:pi/2;
    Psi=imag(W);
    Phi=real(W);
    y=plot(-5:5/10:5,Psi,'r');
    z=plot(-5:5/10:5,Phi,'g');
    r=plot(-5:5/10:5,Rz,'b');
end

Mes erreurs sont certainement très grossières mais je n'arrive pas à m'en sortir.
Pourriez-vous m'expliquer où sont ces erreurs svp,
Merci d'avance