1re année Informatique

Calcul de la puissance d’un moteur

thermique
Remarque : Bien qu’il reste utile de comprendre la démarche de modélisation et de simulation, il
n’est pas nécessaire de comprendre les équations de thermodynamique pour répondre aux questions :
les équations à résoudre numériquement sont clairement données dans les questions et seule l’écriture
d’un programme est demandée. Par ailleurs, toutes les questions sont indépendantes.

I Introduction à la simulation du problème physique

I.1 Présentation du problème physique
Un moteur thermique est un organe de conversion de l’énergie chimique en énergie mécanique.
Un mélange d’air et d’essence est introduit dans la chambre de combusion lorsque le piston est au
plus bas (au point mort bas PMB), il est compressé par le piston jusqu’au point le plus haut (le point
mort haut PMH), puis une étincelle déclenche la combusion qui a pour effet d’augmenter fortement la
température du gaz presque instantanément (et par conséquence la pression). Le piston amorce alors
la descente jusqu’au point mort bas, et capte l’énergie mécanique développée par la pression.

Figure 1 : Compression et détente du gaz dans un moteur thermique. De gauche à droite, position
PMB à θ = 0°, phase de compression, position PMH à θ = 180° et phase de détente.

Du fait de l’inertie des pièces mécaniques, le moteur tourne à vitesse constante à la vitesse de
θ̇ = 3000 tr/min. La loi cinématique liant la position angulaire θ du moteur et la hauteur λ de la
chambre de combustion délimitée par le piston est décrite par :
p
λ(θ) = L + H + R cos θ − L2 − R2 sin2 θ (1)
où L = 0,2 m, R = 0,04 m et H = 0,05 m. Le diamètre du piston vaut d = 0,05 m.
La loi de comportement du gaz est souvent approchée par la loi des gaz parfait pV = nRg T
mais pour une pression de plusieurs dizaines de bars, une loi plus proche de la réalité est choisie,
correspondant au modèle de Van der Waals :

n2 a
(p + )(V − nb) = nRg T (2)
V2
où p est la pression, V est le volume de la chambre, n est le nombre de moles du gaz dans la chambre,
Rg est la constante des gaz parfaits, T est la température et a et b sont des constantes propres au gaz.
Par ailleurs, le gaz échangeant très peu de chaleur avec les parois durant le cycle, les lois thermo-
dynamiques imposent la relation :
T (V − nb)γ−1 = cste (3)

Calcul de la puissance d’un moteur thermique 1/5

 1re année Informatique

où la constante dépend des conditions initiales et γ est appelé coefficient de Laplace des gaz parfaits
γ = 1,4.
À l’entrée dans la chambre, le gaz est à la pression atmosphérique de 1 bar = 105 Pa, une tempé-
2
rature de T0 = 27 °C = 300 K et occupe le volume V0 = Sλ(0) = (R + H) πd4 = 1,8 × 10−4 m3 .
Lorsque le piston remonte, la pression monte et absorbe de l’énergie mécanique. La puissance
mécanique s’exprime comme le produit de la force par la vitesse, c’est-à-dire sous la forme P = pS λ̇
(où p > 0 et λ̇ < 0).
Lorsque le piston est au point mort haut (PMH), la bougie déclenche par une étincelle électrique
la combustion du gaz, réalisant la réaction chimique exothermique :

2C8 H18 + 25O2 −→ 16CO2 + 18H2 O

Lorsque le piston redescend, la pression fournit de l’énergie au piston. La puissance garde la même
expression P = pS λ̇ (où p > 0 et λ̇ > 0).
L’énergie fournie par le gaz au piston est l’intégrale de la puissance sur l’intervalle de temps
correspondant à 1 tour.

I.2 Simulation du problème

On souhaite élaborer un programme permettant de calculer la puissance du moteur à partir de
l’étude thermodynamique du cycle réel.
La démarche de simulation est la suivante :
1. déterminer la courbe de vitesse du piston λ̇ pour 1 tour de moteur,
2. déterminer les conditions thermodynamiques initiales du gaz,
3. déterminer la courbe de pression au cours d’un cycle de montée puis descente du piston,
4. déterminer la courbe de puissance transmise au piston,
5. calculer l’énergie reçue par le piston durant 1 tour de moteur.
Le programme principal est organisé de façon synthétique et fait appel à des fonctions, que nous
allons créer par suite, pour chaque étape du calcul :
1 from numpy import *
2 from matplotlib . pyplot import *
3 # ######################
4 # ## Dimensions du Moteur
5 H = 0.05 # hauteur initiale de la chambre
6 R = 0.04 # rayon de l ’ excentration du vilebrequin
7 L = 0.2 # longueur de la bielle
8 d = 0.05 # diametre du piston
9 Nmot = 3000 # vitesse du moteur en tr / min
10 S = pi * d **2/4 # calcul de la surface du piston
11 theta = linspace (0 ,2* pi ,100) # generation d ’ une serie de position
angulaire
12 w = Nmot * pi / 30 # passage des tr / min en rad / s
13 t = theta / w # generation du vecteur de temps
14 # ## Parametres des gaz
15 a0 = 135.8 e -3 # constantes du melange gazeux avant explosion
16 b0 = 0.0364 e -3
17 a1 = 400 e -3 # constantes du melange gazeux apres explosion
18 b1 = 0.04 e -3
19 Rg = 8.314462 # constante des gaz parfait
20 gamma = 1.4 # constante de Laplace
21 T0 = 300 # temperature initiale ( d ’ admission )
22 P0 = 1 e5 # pression initiale ( d ’ admission )
23 T1 = 2300 # temperature a l ’ explosion
24 # #########################################
25 # ### DEBUT DU PROGRAMME DE SIMULATION

Calcul de la puissance d’un moteur thermique 2/5

 1re année Informatique

26 close ( ’ all ’)
27 # calcul de la geometrie
28 lambda1 = geometrie ( theta , R ) # Q1
29 figure (1)
30 plot ( theta , lambda1 )
31 # calcul de la vitesse
32 vitesse = derivee (t , lambda1 ) # Q2
33 figure (2)
34 plot ( theta , vitesse )
35 # calcul du nombre de moles initiales
36 volume = lambda1 * S
37 n = moles ( P0 , volume [0] , T0 ) # Q4
38 print ( ’ nb moles : ’ ,n )
39 # calcul de la pression
40 pression = vanderwaals ( theta , volume ) # Q5
41 figure (3)
42 plot ( theta , pression )
43 # calcul de la puissance et de l ’ energie
44 puissance = pression * S * vitesse
45 figure (4)
46 plot ( theta , puissance )
47 energie = integration (t , puissance ) # Q6
48 print ( ’L energie mecanique recue est : ’ , energie , ’ Joules ’)

II Écriture du programme de simulation

Les différentes fonctions que vous écrirez pourront faire appel aux variables globales
définies au début du programme. Pour vous aidez, la ligne où est utilisée la fonction à
réaliser sera donnée.

II.1 Loi géométrique

Q1. Écrire une fonction geometrie(theta,R) (ligne 28) prenant en argument le tableau des angles
theta et le rayon R (qui variera lors des dernières questions) et qui retourne un tableau contenant les
hauteurs λ, grâce à l’équation (1).
II.2 Calcul de la vitesse
La puissance étant le produit de la vitesse du piston par la force du gaz sur le piston (P = F λ̇ =
pS λ̇), il est nécessaire de déterminer la vitesse du piston, qui est la dérivée de la position λ par rapport
au temps : λ̇ = dλ dt .
On s’assurera d’avoir un tableau de vitesses de même taille que le tableau de position lambda1 en
prenant la vitesse initiale nulle à t = 0.
Q2. Écrire une fonction derivee(temps,lambda1) (ligne 32) prenant en argument le tableau des
valeurs du temps (temps) et le tableau des positions λ (lambda1) et retourne la dérivation de la
position par rapport au temps avec la méthode de dérivation numérique de votre choix que vous
expliciterez. Vous veillerez à ce que le tableau retourné soit de la même taille que les tableaux temps
et lambda1.

II.3 Calcul du nombre de moles dans la chambre

La relation de comportement du gaz (la loi de Van der Waals) dépend du nombre de moles n dans
la chambre de combustion. n est constant et fixé par les conditions d’admission du gaz (p0 , V0 , T0 )
2
avec p0 = 1 bar, V0 = (R + H) πd4 = 1,8 × 10−4 m3 et T0 = 27 °C = 300 K. Il est donc la solution de
l’équation :
n2 a0
(p0 + )(V0 − nb0 ) = nRg T0 (4)
V02

Calcul de la puissance d’un moteur thermique 3/5

1re année Informatique

Cette équation va être résolue numériquement, par dichotomie. Pour obtenir une estimation relati-
vement proche de la solution, on utilise la loi des gaz parfait pV = nRg T , soit n0 = Rp0gVT00 . L’intervalle
de recherche de la solution sera pris égal à [0,8n0 , 1,2n0 ].
Q3. Définir une fonction f en Python dont vous cherchez le zéro permettant de résoudre l’équation (4).
Vous préciserez bien la nature du ou des argument(s).

Q4. Définir une fonction moles(p0,V0,T0) (ligne 37) prenant en argument la pression p0, le volume
T0 et la température T0 et qui retourne le nombre de moles n incluses dans le piston par dichotomie.
Cette fonction moles devra faire appel à la fonction f et on supposera qu’il y a une unique solution
sur l’intervalle de recherche.

II.4 Calcul de la pression dans la chambre

Une fois la vitesse calculée, il faut déterminer la pression pour obtenir la force du gaz sur le piston.
La pression dépend de la loi de comportement thermodynamique du gaz, la loi de Van der Waals, et
des conditions de transformation, ici traduites par la loi T (V − nb)γ−1 = cste.
En calculant la constante à partir des conditions initiales et en éliminant la température, la pression
est donc, pour le volume V imposé, la solution de l’équation :

n2 a0
(p + )(V − nb0 )γ = nRg T0 (V0 − nb0 )γ−1 (5)
V2
2
où T0 = 300 K, V0 = (R + H) πd4 = 1,8 × 10−4 m3 et pour l’air, a0 = 135,8 × 10−3 Pa · m6 /mol2 et
b0 = 0,0364 × 10−3 m3 /mol.
La formule précédente n’est valable que pour la phase de compression (θ ∈ [0 , 180°]). Après 180°,
elle doit être modifiée pour tenir compte de la combustion du gaz lorsque le piston arrive au point
mort haut : le gaz se transforme essentiellement en H2 O et en CO2 (ce qui modifie les constantes a et
b) et augmente fortement en température pour atteindre plus de 2000 °C.
Tous calculs faits, la pression est, pour le volume V imposé, la solution de l’équation :

n2 a1
(p + )(V − nb1 )γ = nRg T1 (V1 − nb1 )γ−1 (6)
V2
2
où T1 = 2300 K, V1 = (H −R) πd4 = 0,2 × 10−4 m3 et pour les gaz brulés, a1 = 400 × 10−3 Pa · m6 /mol2
et b1 = 0,04 × 10−3 m3 /mol.
Q5. Écrire une fonction vanderwaals(theta,volume) (ligne 40) prenant en argument le tableau des
angles theta et le tableau des volumes volume (défini à la ligne 36) et qui retourne le tableau de
pressions dans la chambre à partir de l’équation de Van der Waals donné par les équations (5) et (6).

II.5 Calcul de l’énergie mécanique reçue par le piston sur un tour

L’énergie est l’intégrale de la puissance fournie au piston sur l’intervalle de temps correspondant
à un tour.
Q6. Écrire une fonction integration(temps,puissance) (ligne 47) prenant en argument le tableau
des valeurs du temps (temps) et le tableau contenant la puissance calculé pour chacun de ces instants.
Cette fonction un nombre réel qui est l’intégrale de la puissance par rapport au temps. Vous expliciterez
la méthode d’intégration choisie.

II.6 Dimensionnement du moteur

On souhaite avoir un moteur dont la puissance moyenne est de 10 kW (la puissance moyenne est
la valeur de l’énergie mécanique obtenue précédemment multipliée par deux et divisée par le temps
mis pour faire un tour, car il y a quatre cylindres, un cycle tous les deux tours).
Le seul paramètre qui peut être modifié dans la conception est la valeur du rayon R.

Calcul de la puissance d’un moteur thermique 4/5

1re année Informatique

Q7. Définir une fonction fobj(R), faisant appel aux fonctions définies précédemment (geometrie,
derivee,moles, vanderwalls et integration) qui permet dans un premier temps de calculer la
puissance moyenne à partir d’un rayon R passé en argument, puis dans un second temps renvoie une
quantité qui devra être nulle si la puissance moyenne est de 10 kW.

Q8. Donner une estimation en octets de la quantité de mémoire nécessaire pour réaliser une évaluation
de la fonction fobj en supposant que tout les tableaux contiennent des nombres réels codés en double
précision. Vous rappelerez tout d’abord le nombre d’octets nécessaire pour représenter un réel en
double précision. On rappelle que theta est un vecteur de 100 éléments.

Q9. Définir une fonction newton2 dont vous préciserez les arguments permettant de trouver le zéro
d’une fonction dont la dérivée sera approchée numériquement.

La fonction fobj précédente est utilisée pour déterminer le rayon optimum mais en fonction de la
valeur initiale la réponse obtenue est soit divergente, soit non physique, soit correcte.
On donne ci-dessous la figure représentant la fonction objectif (celle dont on cherche le zéro pour
répondre au problème posé) en fonction du rayon R ainsi que la courbe de la puissance en fonction
du rayon R.

2500
2000
1500
Fonction objectif en W

1000
500
0
500
1000
1500
2000
0.030 0.032 0.034 0.036 0.038 0.040 0.042 0.044 0.046 0.048
Rayon R en mm

15

10
Puissance en kW

10
0.00 0.01 0.02 0.03 0.04 0.05
Rayon R en mm
Q10. Expliquer pourquoi il y a des problèmes de convergence à partir de ces courbes. Comment
paramétrer la résolution pour s’assurer d’atteindre la solution recherchée ?

Calcul de la puissance d’un moteur thermique 5/5