Bonjour!

Dans le cas d'un problème couplé, je cherche à identifier le terme source dans l'équation de la chaleur. J'ai utilisé un schéma de discrétisation de Crank-Nicholson. Voilà pour l'intro!
Je rappelle la forme de l'équation : rho.Cp.(dT/dt)=div(k.grad(T))+Q

l'expression de Q dans mon cas est : Q= rho*Energie*(d(alpha)/dt) avec alpha qui représente un avancement.

Le problème, c'est que l'expression de alpha est assez compliquée, et il n'y a pas de solution "simple" pour exprimer alpha à l'instant t+dt. Donc, je souhaite utiliser la fonction fminsearch, (ou une autre, je ne suis pas bornée) sachant que alpha est tjrs compris entre 0 et 1.

Voila l'expression en code matlab de alpha à l'instant (k+1) :

Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
 
alpha1=3*dt*K*(alpha2-1)(-ln(1-alpha2))^2/3;
alpha1 est alpha à l'instant k et alpha2 est alpha à l'instant k+1
dt est mon pas de temps
K est une fonction qui dépend de la température

Je connais l'expression de K(T), et je connais la température, du coup, tout devrait rouler, mais je n'arrive pas à écrire bien le bout de code... Pourriez vous m'aider svp?

Je vous remercie par avance!