Discussion :

[Nutsek]

Bonjour,
je suis nouveau et etudiant en classe Prépa.
Je désire modéliser un champ de vecteur vitesse d'une equation issue de la mécanique des fluides a l'aide de Python(x,y) et de quiver.
J'ai les équations mais mon résultats n'est pas bon... je n'arrive pas a obtenir le bon champ...Voici le code :

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
import numpy as np
import matplotlib.pyplot as plt
 
a=1
N=50
u=np.linspace(-5,5,N)
Vx=np.zeros([N,N])
Vy=np.zeros([N,N])
for kx,x in enumerate(u):
    for ky,y in enumerate(u):
        Vx[kx,ky]=5*a*(5**(a-1))*(np.sin((1-a)*3*np.pi)/2)
        Vy[kx,ky]=5*a*(5**(a-1))*(-np.cos((a-1)*3*np.pi)/2)
 
plt.quiver(Vx,Vy)
 
plt.show()

j'ai pris 5 pour la constante Ko qui y a devant et 5 pour r apres changement d'origine (passage d'un r,teta a x,y)
Si vous pouvez m'aider merci beaucoup^!
Paul

[VinsS]

Salut,

Je ne pourrais peut-être pas t'être d'une grande aide dans ton domaine mais une chose me saute aux yeux dans ton code.

La valeur de a ne varie jamais, est-ce normal ?

[Julien N]

Salut,

La valeur de a ne varie jamais, est-ce normal ?

Tout à fait. Tous les vecteurs sont strictements identiques dans ce cas. A quoi doit ressembler ce champs de vecteur? Quelle est l'expression des vitesses. Là je vois (simplifié):

Quelque soit les coordonées (x,y). ça me semble louche.

Ju