% Example 15-1: Finding the Effective Diffusivity
%

pelletporosity=0.4;
CF=0.8;
DAB=0.000001;
tortuosity=linspace(0.1,20,200);
for i=1:200
 EffDiff(i)=pelletporosity*CF*DAB/tortuosity(i);
end
plot(tortuosity,EffDiff(1:200),'LineWidth',2)
ylabel('Effective Diffusivity')
xlabel('Tortuosity')
leg=legend('D_e');
leg.ItemTokenSize=[10,10];
axis([0 20 0 1e-6])
title('Effective Diffusivity Profile')
grid on