Bonjour,
J'essaye du modéliser l'écoulement plan d'un fluide parfait irrotationnel autour d'un rectangle, en supposant que la vitesse du fluide parfait loin en amont du corps est connue et parallèle aux faces horizontales du rectangle (voir schéma ci joint).
J'ai fait un petit code résolvant une équation de Laplace mais les coins de la boîte introduisent des instabilités, j'ignore comment stabiliser ma solution (en fait c'est bêtement une méthode des différences finies de gros bourrin pour l'instant)
Une autre idée serait d'utiliser le principe de superposition et de dessiner un corps fictif (éventuellement coder une méthode, une fois que l'on sait comment approximer le contour cela doit être faisable, à l'aide de doublets, quadripôles, puits /sources,...). Mais j'avoue de je ne suis pas super à l'aise pour dessiner ma ligne de courant sur approximant le contour du rectangle (c'était surtout des configurations courbes et du "conformal mapping" à partir de cylindres en cours).
Merci d'avance
-----
