On considère le mécanisme chimique suivant :
A → X (k1)
2X+Y → 3X (kz)
B + X → Y + D (k3)
X → E (k4)
En considérant
k1=k4=1 s^-1 et k2= k3 =1 M-1.s-1
[A]0= 2 M, [B]o = 5 M, [X]0= 3 M, [Y]0= 0 M, [D]0= [E]o = 0 M,
Avec un pas de temps de 2.10-3s, calculez les variations de concentration de A, B, X, Y, D et E en fonction du temps, avec les méthodes d'Euler explicite et de Runge Kutta à l'ordre 2 jusqu'à atteindre un temps final de 10 secondes.
Dans un cas où A = constante = 1 M et B = constante = 3 M, tracez X et Y en fonction du temps, ainsi que X en fonction de Y, avec des pas de temps de 4 10-3, 4 10-2 et 4 10-1 s. Que remarquez-vous ?
Etudiez la sensibilité des résultats aux différents paramètres initiaux.
Les tracés sont à effectuer sous Excel ou un tableur équivalent, ou en programmation sous Python.
Vous configurerez votre feuille de calcul ou programme de manière à pouvoir faire varier les valeurs des concentrations initiales, constantes et pas de temps en changeant uniquement les valeurs dans une case ou des variables initiales.