L3 – ME-GT – 2022/2023
Résolution des problèmes énergétiques
Travaux Pratiques
Chargé de TP : Thibaut Ménard
Nom et Prénoms : Zigui Evan Jean Christ
PROJET 3 : RESOLUTION D’EQUATIONS DIFFERENTIELLES
ORDINAIRES – APPLICATION A LA SIMULATION
DYNAMIQUE THERMIQUE DU BATIMENT
Introduction
L'analyse numérique est l'étude d'algorithmes permettant résoudre des problèmes continus par
la méthode de discrétisation numérique. Cela signifie qu'il s'agit principalement de répondre à
des questions numériques avec des variables réelles ou complexes comme l'algèbre linéaire
numérique, de trouver des solutions numériques à des équations différentielles (EDO).
L’objectif de ce projet est d’étudier la température de l’air à l’intérieur ainsi qu’à l’extérieur
d’une pièce selon 3 cas en se fixant sur la méthode de résolution d’équation différentielle.
I. Aide à la programmation
Ø Nous allons en premier temps définis les constantes liées aux problèmes pour mieux
comprendre la suite des enjeux
On importe les modules numpy, matpotlib pour effectuer les calculs numériques.
Ø Programmation de la fonction f(t, y) où y = (Ta, Tm, E)
Cette notion est l'une des plus importante en programmation. L'idée est de dire que si telle
variable a telle valeur alors faire cela sinon cela.
La fonction des EDO dy/dt = f(t,y) avec y = (Ta, Tm, E) = (y[0], y[1], y[2] avec Ta :
Température de l'air
Tm : température du mur
E : Energie consommées.
On définit les constantes des équations différentielles.
def f(t, y):
E = A0*(y[1]-y[0]); puis = P(t,E,y[0]) # MODIF CHAUFFAGE
f0 = E + puis # MODIF CHAUFFAGE
f1 = A1*(y[0]-y[1])+B1*(Te(t)-y[1])
f2 = puis*rcva # MODIF CHAUFFAGE
return np.array([f0, f1, f2]) # MODIF CHAUFFAGE
Ø Déterminons la puissance énergétique à travers une trosième
fonction P
2
les conditions if, else et elif sont des notions les plus importante en programmation. L'idée est
de dire que si telle variable a telle valeur alors faire cela sinon cela.
Il est possible de donner des instructions quel que soit les choix possibles avec le mot clé else.
Il est possible d'ajouter autant de conditions précises que l'on souhaite en ajoutant le mot clé
elif , contraction de "else" et "if", qu'on pourrait traduire par "sinon".
Les mots clé if, elif et else cherchent à savoir si ce qu'on leur soumet est True . En anglais
True signifique "Vrai". Donc si c'est la valeur est True, les instructions concernant la
condition seront exécutées.
Ø Par la suite, on prendra comme température extérieure, celle de Rouen en 2013, à
l’aide du fichier distribué. Celui- ci contient deux colonnes. La première contient le
temps en heure, la seconde la température correspondante. Une méthode
d’interpolation devra être employée.
Ø
On defintit la température exterieur :
#Utiliser premiers cas sans prendre exemple de rouen
#def Te(t):
#return TM-dT*np.cos(w*t)
3
Ø Tracer de courbe
Nous pouvons observer sur le tableau de console deux graphiques, ci-dessous :
4
Commentaire : le nombre de points auquel le résultat est évalué n'est pas, au moins
directement lié à la précision des calculs internes.
La consommation énergétique durant ces 10 jours est de 211.0 KWh
Conclusion
Au terme de notre projet on a puis paramétré les différentes constantes d’un bâtiment
thermique via la méthode de résolution d’EDO (équation différentielle) puis à exploiter
numériquement pour en déduire à travers 3 cas a température ainsi que la puissance. Pour
notre part nous avons appris à comment utiliser python pour exploitée ce type de probléme.
5
6
