Cours Méthodes Numériques CH3
Cours Méthodes Numériques CH3
Cours Méthodes Numériques CH3
CH3 INSAT
CHERIF Hayet
Maitre assistant Génie Chimique Procédés
2022-2023
expérimentale.
𝑏
L’objectif est d’approcher la valeur de 𝑥𝑑 𝑥 𝑓 𝑎
h= (b- a)/n
𝑏 𝑓 𝑎 +𝑓(𝑏) ∗ 𝑏−𝑎
=𝑥𝑑 𝑥 𝑓 𝑎 2
du temps. 537
568
0.35
0.34
L'Energie spécifique consommée est calculée par l’équation suivante: 612
649
0.33
0.32
𝐼 𝑡 𝑑𝑡 683 0.31
SPC=E 722 0.3
𝑉 752 0.29
787 0.28
821 0.27
Avec: SPC en KWh/m3 854 0.26
898 0.25
928 0.24
964 0.23
V: volume d’eau 1.5 L 1007 0.22
1055 0.21
E: voltage appliqué (10Volt) 1098 0.2
1149 0.19
1193 0.18
1259 0.17
1318 0.16
1366 0.15
Présenter sur un papier millimétré les données expérimentaux et tracer la courbe 1442 0.14
Appliquer les trois méthodes d’intégration pour calculer le SPC en choisissant un pas 1513
1570
0.13
0.12
du temps faible pour approximer l’intégrale 1628
1791
0.11
0.1
1862 0.09
2004 0.08
2174 0.07
2362 0.06
cours méthode numériques
13
Dr cherif hayet
Application 2
• Quand on veut étudier expérimentalement les caractéristiques de mélange dans un réacteur continu, on peut
utiliser un test consistant en une augmentation rapide de la concentration d'entrée pendant une très courte
période de temps (mathématiquement, une impulsion). La concentration de sortie varie alors en fonction du
temps d'une façon qui sera caractéristique du degré de mélange du réacteur. On verra alors si le réacteur se
comporte comme un réacteur « idéal », par exemple du type CSTR ou PFR. On définit la RTD (Residence
Time Distribution) comme:
• En utilisant les fonctionnalités graphiques de MATLAB, localiser la racine positive de l’équation : f(x) = 2sin(x) − x
• Testez ce script sur matlab (bissectin/bissecting)
a=1.5 ;
b=2 ;
c=(a+b)/2;
tol=1e-6;
iter=0;
while abs(2*sin(c)-c) > tol
if (2*sin(a)-a)*(2*sin(c)-c) <0
b=c;
end
if (2*sin(c)-c)*( 2*sin(b)-b)<0
a=c;
end
c=(a+b)/2;
iter=iter+1;
end
c
iter
a=0;
b=0.5;
f=@(x)exp(x)-4.*x
figure(1),fplot(f,[a,b]);
grid on
• 2-2 A l’aide de la méthode de la sécante, trouver une approximation de cette racine en effectuant 8 itérations.
L'équation qui donne la profondeur x en mètres à laquelle la boule est immergée sous l'eau est donnée par :
x0 donné
F(xn)
xn +1 = xn − F '(x0)
On pourra utiliser la valeur fournie par la loi de Boyle-Mariotte ( pv = RT ) comme valeur initiale. On arrêtera le calcul lorsque deux valeurs successives du volume différeront de
moins de 1%.
On aimerait se servir de cette méthode pour calculer le volume molaire, v, du gaz carbonique à 350°C et à 100 atm à partir de la résolution de l’équation de van der Waals suivante :
T(°C) Cp(J/kg°C)
22 4181
42 4179
52 4186
82 4199
100 4217
1)Déterminer la valeur de la chaleur spécifique à C61°C en utilisant l’approximation polynôme du Lagrange du premier ordre.
2) Déterminer la valeur de la chaleur spécifique à C61°C en utilisant l’approximation polynôme du Lagrange du second ordre.
3) Calculer la valeur absolue de l’erreur entre les deux méthodes
4) Répéter les mêmes opérations pour l’ordre trois.
t=[0:0.5:2.5]
% Euler explicite: ue(t)
ue(1)=u0;
for i=2:6
ue(i)=ue(i-1)+h*f(t(i-1),ue(i-1));
end
plot(t,ue)
xlabel('t')
ylabel('biomasse')
title('Problème à valeur initiale')
grid on