Bonjour, je cherche à trouver une approximation de la résolution de l'équation de la chaleur avec pdsolve (en coordonnées cylindriques), et j'ai un problème avec les conditions initiales .

Pour vous expliquer mon problème, je vais d'abord vous montrer celles avec lesquelles ça marche :
- T(r,0) = 280
- (dT/dr)(r=0,t) = 0
- (dT/dr)(r=3,t) = 0

Dans ce cas tout se passe sans soucis, j'ai une belle courbe en sortie et tout va bien .

Seulement je voudrais modifier la dernière condition (en l'occurence pour modéliser les pertes par rayonnement ) et mettre :
- (dT/dr)(r=3) = -T(r=3,t)^4
(les valeurs numériques sont bidons, c'est juste pour simplifier ici)

Et là, à nouveau ça marche, aucun message d'erreur, seulement j'obtiens exactement la même courbe, ce qui n'est pas normal, le -T(r=3,t)^4 est considéré nul (ce qui est loin d'être le cas ) .

Essayant de comprendre le problème, j'essaie différentes manipulations, en voilà les résultats :
- en changeant cette condition en : (dT/dr)(r=3,t) = -T(r=3,t), cette condition est bien prise en compte, et j'ai une pente assez raide en fin de courbe
- dès que l'on met une puissance sur T(r=3,t), Maple l'ignore