Problema 4

Concretamente, se llama cinética o cinemática a la especialidad de la física centrada en el análisis del movimiento, dejando de lado el estudio de las fuerzas que lo originan. … La cinética, en este marco, suele estudiar la trayectoria de un cuerpo en movimiento en función del tiempo.

library(tidyverse)

## -- Attaching packages --------------------------------------- tidyverse 1.3.1 --

## v ggplot2 3.3.5     v purrr   0.3.4
## v tibble  3.1.6     v dplyr   1.0.7
## v tidyr   1.1.4     v stringr 1.4.0
## v readr   2.1.1     v forcats 0.5.1

## -- Conflicts ------------------------------------------ tidyverse_conflicts() --
## x dplyr::filter() masks stats::filter()
## x dplyr::lag()    masks stats::lag()

Problema 4

A + B -> C

A=0.2
B=0.15
C=0
k1=0.1
dt=0.01
t=0

CC=c()
tt=c()
AA=c()
BB=c()
while(t<1000){
  AA=c(AA,A)
  BB=c(BB,B)
  CC=c(CC,C)
  tt=c(tt,t)
  
  dC= k1*A*B*dt
  dB= -k1*A*B*dt
  dA= -k1*A*B*dt
  
  C=C+dC
  t=t+dt
  A=A+dA
  B=B+dB
}
plot(tt,CC, type="l", lwd=2, xlim=c(0,500), ylim=c(0,0.2))
lines(tt,AA, type="l", col=2, lwd=2)
lines(tt,BB, type="l", col=6, lwd=2)

A+B-> C+S-> X

A=0.2
B=0.15
C=0
S=0.02
X=0
k1=0.1
k2=0.2
dt=0.01
t=0

CC=c()
tt=c()
AA=c()
BB=c()
SS=c()
XX=c()

while(t<1000){
  AA=c(AA,A)
  BB=c(BB,B)
  CC=c(CC,C)
  XX=c(XX,X)
  SS=c(SS,S)
  tt=c(tt,t)
  
  
  
  dC= k1*A*B*dt
  dB= -k1*A*B*dt
  dA= -k1*A*B*dt
  dS= dC - k2*C*S
  dX= -dC + k2*C*S
  
  C=C+dC
  t=t+dt
  A=A+dA
  B=B+dB
  S=S+dS
  X=X+dX
}

plot(tt,CC, type="l", lwd=2, xlim=c(0,500), ylim=c(0,0.2), col= "magenta4")
lines(tt,AA, type="l", col=2, lwd=2)
lines(tt,BB, type="l", col=6, lwd=2)
lines(tt,SS, type="l", col="tan", lwd=2)
lines(tt,XX, type="l", col=5, lwd=2)
legend(450,0.2, legend=c("[C]","[A]","[B]","[S]","[X]"), col=c("magenta4",2,6,"tan",5),lwd=3 , cex=0.5)

A+B-> C+S-> X -> S+D

A=0.2
B=0.15
C=0
S=0.02
X=0
D=0
k1=0.1
k2=0.2
k3=0.05
dt=0.01
t=0

CC=c()
tt=c()
AA=c()
BB=c()
SS=c()
XX=c()
DD=c()

while(t<500){
  AA=c(AA,A)
  BB=c(BB,B)
  CC=c(CC,C)
  XX=c(XX,X)
  SS=c(SS,S)
  DD=c(DD,D)
  tt=c(tt,t)
  
  
  
  dC= k1*A*B*dt
  dB= -k1*A*B*dt
  dA= -k1*A*B*dt
  dS= dC - k2*C*S
  dX= -dC + k2*C*S - k3*S*D*dt
  dD= k3*S*D*dt
  
  C=C+dC
  t=t+dt
  A=A+dA
  B=B+dB
  S=S+dS
  X=X+dX
  D=D+dD
}

plot(tt,CC, type="l", lwd=2, xlim=c(0,500), ylim=c(0,0.2), col= "magenta4", main="A+B-> C+S-> X -> S+D", xlab=c("Tiempo (s)"), ylab=c("Concentraciones Molares (M)"))
lines(tt,AA, type="l", col=2, lwd=2)
lines(tt,BB, type="l", col=6, lwd=2)
lines(tt,SS, type="l", col="tan", lwd=2)
lines(tt,XX, type="l", col=5, lwd=2)
lines(tt,DD, type="l", col=7, lwd=2)
legend(450,0.2, legend=c("[C]","[A]","[B]","[S]","[X]"), col=c("magenta4",2,6,"tan",5),lwd=3 , cex=0.5)