#!/usr/bin/env python3
# -*- coding: utf-8 -*-
"""
Created on Tue Nov 22 16:13:40 2022

@author: yvestalbourdel&maximeperrier
"""
import numpy as np
import matplotlib.pyplot as plt

## valeurs numériques des grandeurs liées au problème étudié

Cperox=2.7 #mol/l
Cdichro=0.003 #mol/l
V=0.5#l
rho=1000 #g/l
m=rho*V #g
n0 = Cperox*V #mol
gamma=0.085 # rapport sans dimension
cp=4.18 # J/g/K
ctot=(gamma+1)*m*cp # J/K
G=0.27 #W/K
Tp=291

Delta_rHo = -95 * 10**3 # J/mol


def K(T) :
    return ##########

def k(T) :
    return ###########

def d_ksi(T,ksi,dt):      # fonction qui calcule d_ksi pour une durée dt à partir d'un état donné
    return ##############

def d_T(T,dksi,dt):     # fonction qui calcule d_T pour une durée dt et un déplacement dksi à partir d'un état donné
    return ##################

T_i = 291                    # température initale (K)
t_i = 0                      # instant initial (s)
t_f = 3600                     # instant final (min)
dt = (t_f-t_i)/10000

Instant=[t_i]
Temperature=[T_i]
Avancement=[0]
t=Instant[-1]
ksi=Avancement[-1]
T=Temperature[-1]

while t < t_f : 
    t = t + dt
    Instant.append(####)
    dksi = d_ksi(########)
    ksi = ksi + dksi
    Avancement.append(######)
    dT = d_T(############)
    T = T + dT
    Temperature.append(######)
       
plt.figure(1,figsize=(10,8))
plt.plot(Instant,Temperature,label=f"[H2O2]={Cperox}mol/L [Cr2O7]={Cdichro}mol/L")
plt.xlabel(r'Instant $t$ (secondes)',fontsize=14)
plt.ylabel(r'Température $T$ (K)',fontsize=14)
plt.title(r'Évolution temporelle de la température')
plt.legend()
plt.show()