# Exemple simple de calcul de l'évolution de la température dans un réacteur
# chimique avec une cinétique d'ordre 1 et une capacité calorifique fixe

import numpy as np
import matplotlib.pyplot as plt

drH=-122e3
C=1e5

duree=600       # sur 10 minutes
N=1000
dt=duree/N

t=np.linspace(0,duree,N)
T=np.zeros(N)
n=np.zeros(N)    # moles de réactif
T[0]=298
n[0]=500


# loi d'Arrhenius pour la constante k de vitesse
def k(T):
    return np.exp(-(20e3/(8.31*T)))

# méthode d'Euler
for i in range(N-1):
    #-------------------------à compléter-------------------------------------






    #-------------------------------------------------------------------------


# affichage du graphe T(t)
plt.figure(figsize=(10,4))
plt.title("Température dans le réacteur")
plt.xlabel("Temps en s")
plt.ylabel("Température en K")
plt.plot(t,T)
plt.grid()
plt.show()

# affichage du graphe n(t)
plt.figure(figsize=(10,4))
plt.title("nombre de moles de réactifs restant")
plt.xlabel("Temps en s")
plt.ylabel("mol")
plt.plot(t,n)
plt.grid()
plt.show()