Lab 4 design - end, x(i-1),v(i-1), plot(x, t) p v0=0.002,

Info iconThis preview shows page 1. Sign up to view the full content.

View Full Document Right Arrow Icon
v0=0.002, k0=0.05, x(1)=0,v(1)=0,dx=0.01 for i=2:100, x(i)=x(i-1)+dx, dv=dx/(k0/v0*(1-x(i))), v(i)=v(i-1)+dv, if v(i)>0.005, break, end, end, x(i-1), v(i-1), plot (x, v); p v0=0.002, k0=0.05, x(1)=0,v(1)=0,dx=0.01,t(1)=1035, fa0=0.0376, v0=0.002, v=5, ua=16500, ha0=-216.67, hb0=-61.09, hc0=-74.81, cpa=163, cpb=83, cpc=80, t0=1035, ta=1150, deltah=hb0+hc0-ha0, deltacp=cpb+cpc-cpa, ca0=fa0/v0, for i=2:101, x(i)=x(i-1)+dx, dv=dx/(k0/v0*(1-x(i))), v(i)=v(i-1)+dv, dt=(dv*ua*(ta-t(i))+k0*ca0*(1-x(i)))*((-(deltah)+deltacp*(t(i)-298)))/ (fa0*(cpa+x(i)*(cpb+cpc-cpa))) t(i)=t(i-1)+dt, if v(i)>0.005, break,end,
Background image of page 1
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: end, x(i-1),v(i-1), plot(x, t) p v0=0.002, k0=6*exp(-5000/t(i)), x(1)=0,v(1)=0,dx=0.01,t(1)=1035, fa0=0.0376, v0=0.002, v=5, ua=16500, ha0=-216.67, hb0=-61.09, hc0=-74.81, cpa=163, cpb=83, cpc=80, t0=1035, ta=1150, deltah=hb0+hc0-ha0, deltacp=cpb+cpc-cpa, ca0=fa0/v0, for i=2:101, x(i)=x(i-1)+dx, dv=dx/(k0/v0*(1-x(i))), v(i)=v(i-1)+dv, dt=(dv*ua*(ta-t(i))+k0*ca0*(1-x(i)))*((-(deltah)+deltacp*(t(i)-298)))/ (fa0*(cpa+x(i)*(cpb+cpc-cpa))) t(i)=t(i-1)+dt, if v(i)>0.005, break,end, end, x(i-1),v(i-1), plot(x, t) p...
View Full Document

This note was uploaded on 10/19/2010 for the course SOC 4523 taught by Professor Jacobs during the Spring '10 term at Wisc Eau Claire.

Ask a homework question - tutors are online