Iterate simultaneous equations till it reach convergence

Question

Senthil 2019-5-17

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/462560-iterate-simultaneous-equations-till-it-reach-convergence

回答： Senthil 2019-5-17

采纳的回答： Senthil

在 MATLAB Online 中打开

Question: How to iterate T_R, T_f and T_a for multiple times till T_f reaches convergence.

Code below:

Q_R=935.9432834;
h_f=2500;
A_f=0.017671459;
A_a=0.017671459;
h_cs=13.3;
L_c=0.15;
h_rs=8.8;
R_c=0.075;
T_amb=35;
E_bm=0.05;
SIGMA=5.67E-08;
T_sky=29;
h_a=10;
h_c=10;
h_ct=10;
m_f=2;
c_f=4186;
A_s=0.070685835;
A_top=0.017671459;
hct=8.8;
t_lid=0.01;
k_lid=237;
m_a=0.000797131;
c_a=718;
DELTAt=1;
m_R=0.5;
c_R=910;
T_R0=25; %initial condition"	
T_f0=25; %initial condition"
T_a0=25; %initial condition"
T_R=(Q_R-(h_f*A_f)*(T_R0-T_f0)-((h_cs*2*L_c)+(h_rs*R_c))*3.14*R_c*(T_R0-T_amb)-(2*L_c+R_c)*(3.14*R_c*E_bm*SIGMA*(T_R0^4-T_sky^4))/(T_R0-T_amb))*(DELTAt)/(m_R*c_R)+T_R0;
T_f=(((h_f*A_f*(T_R0-T_f0))-(h_c*A_a*(T_f0-T_a0)))*DELTAt/(m_f*c_f))+T_f0;
T_a=(((h_a*A_a)*(T_f0-T_a0))+((h_a*A_s)*(T_R0-T_a0))-(A_top*(T_a0-T_amb))/(1/(h_ct+((E_bm*SIGMA*(T_R0^4-T_sky^4))/(T_a0-T_amb)))+(t_lid/k_lid)))*(DELTAt)/(m_a*c_a)+T_a0;