Bonjour à tous,
J'essaie de résoudre l'équation de la chaleur en 1D instationnaire avec convection.
En fait, on retombe sur l'équation de l'ailette :
kSd²T/dx²=hP(T-Text)
Avec k la conductivité du matériau de la barre 1D, S la section, h le coefficient de convection et Text=90 la température extérieure.
Je trouve :
T(x)=A*exp(ax)+B*exp(-ax)+Text
Avec a = (hP/kS)^0.5
Jusque là tout va bien. Ensuite je détermine A et B en fonction des conditions initiales, et là...
Lorsque je veux poser T(0)=160 par exemple et dT(0)/dx=0 (barre isolé au début), ça déconne : la température démarre bien à 160 mais elle... augmente! Sachant que le milieu est à Text=90, soit un milieu plus froid que la température imposée au début, ça devrait refroidir. Et pourtant ça augmente. Je comprends pas.
Du coup j'ai dégainé Matlab pour voir et... ça fait la même chose.
Je comprends pas. Est-ce à dire que je peux pas imposer la température et sa dérivée au même point?
Je pense pas, parce que dans une équa diff d'ordre 2 comme la trajectoire d'un point par exemple, on peut très bien spécifier position et vitesse (fonction et dérivée, donc) au même point. Bon, certes, c'est pas la même équation...
Parce que quand j'impose les températures aux deux extrémités, ou la température d'un côté et la dérivée de l'autre, ça se comporte bien. Il y a quelque chose que je capte pas.
Par exemple, j'ai a=5. J'obtiens A=B=35. Quand je trace ma fonction, je démarre bien de 160 avec une pente nulle, mais ça augmente. J'ai une courbure positive.
Où est-ce que je me plante?
Je vois bien pourquoi ça explose : le terme en exp(+x) devient prépondérant, je devrais donc trouver un terme B beaucoup plus faible...
Mais j'ai A=B car je regarde la dérivée en 0. Si je la regarde un poil plus loin, le terme B va être plus faible grâce à l'expo en -x...
Merci!
-----