#Décomposition du peroxyde de di tertiobutyle#

import numpy as np
import numpy as np
import matplotlib.pyplot as plt


# Données  cinétiques  : réaction d'ordre 1 #
#contante de vitesse : k = k0exp(-Ea/RT)#
k0=10**15             #unité : s-1#
Ea=157000             #unité : J/mol#

#Donnée thermodynamique#
DrH0=- 150000         #unité : J/mol#
R = 8.314           #unité : J/K/mol#
Cp = 2100           #capacité thermique massique moyenne du système chimique  , J/K/kg #
mu = 900            # masse volumique   , kg/m3 #

#Donnée sur le réacteur : RPAC #
VR = 5.2e-4      #volume du réacteur en m3#
#Alimentation #
CRe=6e3        #concentration molaire en entrée du réactif (peroxyde) en mol/m3#
tau = 600    #temps de passage en s#
Te1=600      #-température d'entrée  en kelvin#

#1-Variation du taux de conversion X en fonction de la température  T  : aspect cinétique"
#A compléter pour pouvoir tracer Xcin en fonction de la température




#tracé des varaiations du taux de conversion en fonction de la température #
plt.figure(0)
plt.plot(Temp,Xcin,'+b',markersize=1,label="aspect cinétique")
plt.xlabel("T(K)")
plt.ylabel("X cinétique")
plt.grid()
plt.legend()
plt.show()

#2-Aspect thermodynamique , fonctionnement adiabatique "
#A compléter pour pouvoir tracer Xadia en fonction de la température



plt.figure(1)
plt.plot(Temp,Xadia,'+r',markersize=1,label="fonctionnement adiabatique")
plt.xlabel("T(K)")
plt.ylabel("X adiabatique")
plt.grid()
plt.legend()
plt.show()

plt.figure(2)
plt.plot(Temp,Xadia,'+r',markersize=1,label="aspect thermodynamique")
plt.plot(Temp,Xcin,'+b',markersize=1,label="aspect cinétique")
plt.xlabel("T(K)")
plt.ylim(0,1)
plt.ylabel("X adiabatique")
plt.grid()
plt.legend()
plt.show()

#3-Aspect thermodynamique, avec refroidissement par l'eau #
S= 314e-4  #surface d'échange  en m2#
U = 80  # coefficient d'écange    , W/m2/K
Teau=293    #température de l'eau en K#

#A compléter pour pouvoir tracer les varaiations du taux d econcersion Xechane  en fonction de la température








plt.figure(3)
plt.plot(Temp,Xechange,'+r',markersize=1,label="aspect thermodynamique")
plt.plot(Temp,Xcin,'+b',markersize=1,label="aspect cinétique")
plt.xlabel("T(K)")
plt.ylim(0,1)
plt.ylabel(" taux de conversion")
plt.grid()
plt.legend()
plt.show()