Cinética Química

Cinética de Primer Orden

\[ \frac{d[A]}{dt}=-k_1*[A] \]

A=1.0
P=0
k1=0.01
dt=0.1
t=0

AA=c()
tt=c()
PP=c()
while(t<500){
  AA=c(AA,A)
  PP=c(PP,P)
  tt=c(tt,t)
  
  dA=-k1*A*dt
  dP=k1*A*dt
  
  A=A+dA
  P=P+dP
  t=t+dt
}

plot(tt,AA,type="l",col="blue",ylim=c(0,1.1),xlab="t(min)",ylab="[A]",main="Cinética de Primer Orden")
lines(tt,PP,type="l",col="red")

Cinética de Segundo Orden

\[ \frac{d[A]}{dt}=-k_1*[A]^2 \]

A=1.0
P=0
k1=0.1
dt=0.1
t=0

AA=c()
tt=c()
PP=c()
while(t<50){
  AA=c(AA,A)
  PP=c(PP,P)
  tt=c(tt,t)
  
  dA=-k1*A*A*dt
  dP=k1*A*A*dt
  
  A=A+dA
  P=P+dP
  t=t+dt
}

plot(tt,AA,type="l",col="blue",ylim=c(0,1.1),xlab="t(min)",ylab="[A]",main="Cinética de Primer Orden")
lines(tt,PP,type="l",col="red")

Cinetica de Reacción consecutiva A –> B –> P

A=1.0
B=0
P=0
k1=0.1
k2=0.05
dt=0.1
t=0

AA=c()
tt=c()
PP=c()
BB=c()

while(t<50){
  AA=c(AA,A)
  PP=c(PP,P)
  BB=c(BB,B)
  tt=c(tt,t)
  
  dA=-k1*A*dt
  dB=k1*A*dt-k2*B*dt
  dP=k2*B*dt
  
  A=A+dA
  B=B+dB
  P=P+dP
  t=t+dt
}

plot(tt,AA,type="l",col="blue",ylim=c(0,1.1),xlab="t(min)",ylab="[A]",main="Cinética de Primer Orden")
lines(tt,BB,type="l",col="green")
lines(tt,PP,type="l",col="red")