Résolution numérique d'une EDP shéma implicite.
Bonjour,
Je bloque. Je dois continuer un schéma implicite de résolution d'une EDP par itération. C'est pour modéliser un champ de vitesse dans une artère avec l'équation de Navier Stokes (puis après autour d'une bosse dans l'artère). J'ai appliqué des conditions de vitesse initiales mais:
- comment faire pour que à chaque point atteint on puisse appliquer une modification sur le point précèdent, le point suivant, mais aussi le point du dessus et du haut? Car dans les exemples avec l'équation de la chaleur on applique seulement une modification aux points suivants et non précédents ou au dessus:
ex là ça ne va pas:
Code:
1
2
3
4
5
6
7
8
9
10
11
12
| for i in range(nx + 1):
for j in range(ny + 1):
if i == j:
A[i, j] = A5
if i + 1 == j:
A[i, j] = A6
if i - 1 == j:
A[i, j] = A4
if j + 1 == i:
A[i, j]= A8
if j - 1 == i:
A[i, j] = A2 |
Merci beaucoup, je n'arrive vraiment pas à trouver...
