Iterate simultaneous equations till it reach convergence
2 visualizzazioni (ultimi 30 giorni)
Mostra commenti meno recenti
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;
Expected Output Required:
Here the output is generated for 20 time steps. But we need to iterate till the difference between T_f upto 0.01 or 0.02.
Risposta accettata
Più risposte (0)
Vedere anche
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!