Bonjour !
Il me faut résoudre un problème 2D pour équation du type parabolique
du/dt-div(Kgradu)=0 ;
avec les conditions limites de Neumann
Kdu/dn=J ;
en utilisant la méthode des éléments finis. Il fallait que la géométrie du domaine de calcul puisse être arbitraire et que le coefficient soit une fonction constante en morceaux dans le domaine. J'ai essayé d'écrire un code dans le milieu Matlab (je sais que dans Matlab il y a des fonctions intégrées qui peuvent le faire, mais je voulais le faire mais je voulais y arriver moi-même pour pouvoir faire le même sans Matlab), je l'ai pas réussi. Je Vous demande s'il existe des codes déjà faits qui font ça, qui peuvent résoudre mon problème ? Avec les commentaires bien détaillés, préférable que ça soit en langage de Matlab, comme j'ai peu de connaissances en programmation. Tout ce qu'on peut trouver dans les livres ce sont des trucs trop fondamentaux, là on donne peu d'attention aux questions : comment approcher les conditions limites, un code existant était plus pratique et pourrait m'aider plus qu'un livre. S'il y en a, dites-moi où je peux le trouver. Merci