% LEP 14-3 : Combustion time for a single particle
%

rhoc=2226000;
DAB=0.0001;
dp0=100;
phic=1;
CAinf=linspace(1,20,100);

for i=1:100
rhocm = rhoc/12;
Ks(i) = 8*DAB*CAinf(i)/(rhocm*phic)*10^9;
tc(i) = dp0^2/Ks(i);
end
tiledlayout(1,2)
set (gcf,'Position',[0,0.1,600,300])
nexttile
plot(CAinf,tc(1:100),'LineWidth',1.5);
xlabel('C_{A\infty} (mol/m^3)')
ylabel('t_c (ms)')
axis([0 20 0 1000])
title('Combustion time')
grid on
nexttile
plot(CAinf,Ks(1:100),'LineWidth',1.5);
xlabel('C_{A\infty} (mol/m^3)')
ylabel('K_s (\mum^2/ms)')
axis([0 20 0 100])
title('Burning Rate constant')
grid on