Probleme capteur solaire sphérique

Fermé
HSANAAKAHLIL Messages postés 1 Date d'inscription vendredi 7 juin 2013 Statut Membre Dernière intervention 7 juin 2013 - 7 juin 2013 à 17:59
bonjour
j'ai besoin un programme sous Matlab pour résoudre le bilan thermique d'un capteur solaire sphérique , j'ai trouvé dans la littérature que la méthode qui peut être résoudre ce problème là est la méthode itératif , j'ai aucun idée a cette méthode , si quelque un me aider ça me fais un grand plaisir . mrc
voila les équations que j'ai trouvé sous forme matricielle :
KT(1,1)=-(sigma*4*pi*(Ray2)^2/((1/Emmis2)+((1-Emmis1)*(Ray2)^2)/(Emmis1*(Ray1-Epaiss1)^2)))*T(1,1)^3-0.697*4*pi*(Ray2)^2*((T(2,1)-T(1,1))^(1/4))-hcv1*4*pi*(Ray1)^2-sigma*Emmis1*4*pi*(Ray1)^2*T(1,1)^3;
KT(1,2)=(sigma*4*pi*(Ray2)^2/((1/Emmis2)+((1-Emmis1)*(Ray2)^2)/(Emmis1*(Ray1-Epaiss1)^2)))*T(2,1)^3+0.697*4*pi*(Ray2)^2*((T(2,1)-T(1,1))^(1/4));
KT(1,3)=0;
KT(2,1)=0.697*4*pi*(Ray2)^2*((T(2,1)-T(1,1))^(1/4))+(sigma*4*pi*(Ray2)^2/((1/Emmis2)+((1-Emmis1)*(Ray2)^2)/(Emmis1*(Ray1-Epaiss1)^2)))*T(1,1)^3;
KT(2,2)=-hcv3*4*pi*(Ray2-Epaiss2)^2-0.697*4*pi*(Ray2)^2*((T(2,1)-T(1,1))^(1/4))-(sigma*4*pi*(Ray2)^2/((1/Emmis2)+((1-Emmis1)*(Ray2)^2)/(Emmis1*(Ray1-Epaiss1)^2)))*T(2,1)^3;
KT(2,3)=hcv3*4*pi*(Ray2-Epaiss2)^2;
KT(3,1)=0;
KT(3,2)=hcv3*4*pi*(Ray2-Epaiss2)^2;
KT(3,3)=-hcv3*4*pi*(Ray2-Epaiss2)^2-2*debit*cp4;
%
F(1,1)= -abs1*2*pi*(Ray1)^2*Gin*-hcv1*pi*(Ray1)^2*Tamb-sigma*Emmis1*4*pi*(Ray1)^2* Tciel^4;
F(2,1)= -2*pi*(Ray2)^2*trans1*abs2*Gin;
F(3,1)= -2*debit*cp4*Tfe;