ADS5

REPUBLIQUE DEMOCRATIQUE DU CONGO
UNIVERSITE PROTESTANTE DE LUBUMBASHI

FACULTE POLYTECHNIQUE
TRAVAIL PRATIQUE D’ANALYSE DES SYSTEMES CONTINUS
Dirige par : PhD Jacques KAZAKU

Effectue par : SANGWA MVUMVU Eliel
Promotion : BAC 3 GMP (Métallurgie)
ANNEE ACADEMIQUE : 2023-2024

Etude d’un réacteur chimique
1. Modéliser un réacteur batch de volume constant dans lequel se déroule de façon isotherme la
réaction
A −k→ B r = kcA
avec les conditions initiales suivantes:
cA(0) = cA0, cB(0) = cB0
a) Montrer que pour ce système réactionnel comportant deux composants mais une seule
réaction, il est nécessaire d’écrire un bilan molaire pour le composant B, car cB(t) se
calcule directement à` partir de cA0, cB0 et cA(t).
b) Considérer maintenant le système réactionnel isotherme suivant se d´déroulant dans un
réacteur semi-batch avec volume variable :
A + B −k→ 2C r = kcAcB
avec cA(0) = cA0, cB(0) = cB0, cC(0) = cC0, V (0) = v0. B est alimenté avec le d´débit
d’alimentation q de concentration cB,e.
c) Est-il toujours possible d’écrire un seul bilan molaire pour d´écrire entièrement le
système dynamique, c’est-a`-dire pour d´ecrire cA(t), cB(t) et cC(t) ?
d) Généraliser les conclusions précédentes par rapport au nombre de réactions et au type
de réacteur (batch, feb-batch et continu).
SOLUTION
a) Réacteur batch
A−k→ B r=Kca
Conditions initiales : Ca (o) = Ca ,Cb (o) = Cbo
Bilan global
A=E-S+G-D Qm=Qout =0
dm
=Qa-Qb V= Cste
dt
 Bilan molaire pour A
dVCa
= 0-0-KCaV
dt
VCa
= -KCaV
dt
dCa
= -KCa
dt
 Bilan molaire pour B
dVCb
= 0-0-KCbV
dt
VCb
= -KCbV
dt
dCb
= -KCb
dt
Cb Ca
dCa dCb dCb
=- →∫ = - ∫ dCa
dt dt Cbo dt Cao dt
→Cb(t) – Cbo = - Ca(t) – Cao
Cb(t) = (Cao + Cbo) – Ca(t)
b) A + B −k→ 2C rKCaCb : Réacteur semi-batch

Condition initiales : Ca(o) = Cao, Cb (o) = Cbo ; Cc (o)
Modelé dynamique : en considérant le nombre de mole
dV
= Q V(o) = Vo
dt
 Bilan molaire pour A

dna na∗nb
= -K na(o) = nao = VCao
dt V
 Bilan molaire pour B
dnb na∗nb
= Q Cbin – K nb(o) = nbo = VCbo
dt V
 Bilan molaire de C
dc na∗nb
= 2K nc(o) = nco = VCco
dt V
c) Non, il n’est pas toujours possible d’écrire un seul bilan molaire pour décrire entièrement le
système dynamique, car deux équations différentielle pour Ca(t) et V(t) sont nécessaires
pour décrire le système dynamique.
Vo Vo
Cb(t) = Cbin [1 - ]+ (Cbo – Cao ) + Ca(t)
V (t) V (t)
d) Le nombre minimal de variables d’état pour un système réactionnel dans un réacteur r + p
ou r représente le nombre de débits d’alimentation de sorties importantes.
Etude d’un banc de flottation

On considère le banc de flottation donne à` la figure ci-dessous. Comme les deux cellules ont les
mêmes géométries
Figure 1 : Banc de flottation `a deux cellules
On considère que M2 −M1 est proportionnel à la différence de pression au niveau de l’orifice de

communication entres les cellules. On considère les modèles d’écoulement suivant :
Q1 = k1(M1 − M2)
Q2 = k2M2
Donner les fonctions de transfert de M2(s)/Q(s) et M1(s)/Q(s).
SOLUTION
La conservation de la pulpe s’écrit :
dM 1
 =Q−k 1 ( M 1−M 2 )
dt
dM 2
 =k−k 1 ( M 1−M 2 ) −k 2 M 2
dt
En régime permanent on a :
 Q1° = Q2° = Q°
Q°
 M2° =
K2
k 1+ k 2
 M1° = Q°
k1k 2
Si on définit les débits et masses comme des variations autour des régimes permanent, les
valeurs des variables au repos sont nulles, ce qui permet d’écrire les transformées de
LAPLACE des deux équations de conservation :
 {sMsM2=k
1=Q−k 1 M 1+k 1 M 2
1 M 1−( k 1+k 2 ) M 2
dont la résolution algébrique conduit à :
k1
 M 2( s) = 2
Q( s)
s + ( 2 k 1+k 2 ) s +k 1 k 2
s+ k 1+ k 2
 M 1 ( s) = 2 Q( s)
s + ( 2 k 1+k 2 ) s +k 1 k 2
Ces deux fonctions de transfert permettent de calculer M1(t) et M2(t) pour des variations
de débit Q(t) dont on conduit la transformée de LAPLACE (par exemple un échelon de
débit). Elles permettent aussi de calculer le débit Q2(t) qui est donne par :
1
Q 2 ( s )= Q(s)
 1 2 2 k 1+ k 2
s+ s +1
k 1k 2 k1k 1
dont la réponse à un échelon d’amplitude A(Q(t) passe de Qo a Qx + A), est d’après
les tables de la transforme de la LAPLCA
A Q1t Q2t
 Qx ( t )=Qo+ (e −e ) Ou S1 et S2 sont des pôles de la fonction de
S1S 2
transfert
−2 k 1−r 2± √ 4 r 12 +r 2²

2
Etude du bioréacteur´
Considérons une cuve de volume V alimentée par un d´débit volumique d’entrée q, et une
réaction de croissance de la biomasse X par la synthèse d’un substrat S donnée par :
S → X.
La dynamique du système est d´écrite par les équations différentielles non linéaires de bilan
matière suivantes :
Ou µ est le taux spécifique de croissance et Sin est la concentration en substrat dans

l’alimentation. Le taux spécifique de croissance peut être modélise par diverses expressions ;
dans le cadre de cet exercice nous utiliserons le modèle de Haldane donne par l’expression
suivante :
Avec KS la constante d’affinité, KI la constante d’inhibition et µ0 un paramétré qui peut être

associe au taux spécifique de croissance maximum. La figure 2 illustre le modèle cinétique
de Haldane. L’entrée du systèmes est
0.1
0.08
(1/h)
0.06
0.04
0.02
0
0 50 100 150 200 250 300 350 400 450 500 550
S (g/L)
Figure 2: Modelé Haldane – KI = 0.5, µ0 = 4.1, KS = 200, K1 = 2.
u = Sin, la sortie y = S et la perturbation v = V q .

Pour l’application numérique, les valeurs numériques des paramètres sont données a` la
Table 1.
KI 0.5 g/L
µ0 4.1 h−1
KS 200 g/L
k1 2 (sans unités)
Table 1: Valeurs numériques des paramètres pour le bio réacteur
1) Soit u = u⋆ g/L et v = v⋆ h−1, des valeurs fixes connues. Dans le cas v∗ > 0, identifier les
trois points d’´équilibre du système tels que
.
Remarquez que les points d’équilibre doivent être écrits comme des fonctions de u⋆ et
de v⋆. Pour l’application numérique, u∗ = Sin∗ = 500 g · L−1 et v∗ = (V q )∗ = 0.085 h−1.
Identifier les équilibres sur le graphique du modèle cinétique de Haldane (figure 2).
Combien d’équilibres aurez-vous si v∗ = 0?
2) Donner les représentations d’état du système linéarise autour de chaque équilibre
3) Etudier la stabilité des systèmes linéarise. Que pouvez-vous en conclure sur la stabilité
des équilibres ? Etudier et interpréter la réponse indicielle des systèmes linéarises.
4) Calculer les fonctions de transfert correspondantes. Etudier les zéros des systèmes
linéarisent.
5) Refaire les questions 3 et 4 en considérant la biomasse X comme sortie du système.
A titre d’information, la représentation d’état (linéaire) du bio-réacteur est de la forme :

avec

ADS5

Transféré par

Droits d'auteur :

Formats disponibles

ADS5

Transféré par

Informations du document

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

ADS5

Transféré par

Droits d'auteur :

Formats disponibles

REPUBLIQUE DEMOCRATIQUE DU CONGO

UNIVERSITE PROTESTANTE DE LUBUMBASHI

TRAVAIL PRATIQUE D’ANALYSE DES SYSTEMES CONTINUS

Dirige par : PhD Jacques KAZAKU

ANNEE ACADEMIQUE : 2023-2024

cA(0) = cA0, cB(0) = cB0

 Bilan molaire pour A

 Bilan molaire pour B

→Cb(t) – Cbo = - Ca(t) – Cao

Cb(t) = (Cao + Cbo) – Ca(t)

b) A + B −k→ 2C rKCaCb : Réacteur semi-batch

 Bilan molaire pour A

Etude d’un banc de flottation

Figure 1 : Banc de flottation `a deux cellules

On considère que M2 −M1 est proportionnel à la différence de pression au niveau de l’orifice de

La conservation de la pulpe s’écrit :

Ou µ est le taux spécifique de croissance et Sin est la concentration en substrat dans

Avec KS la constante d’affinité, KI la constante d’inhibition et µ0 un paramétré qui peut être

Figure 2: Modelé Haldane – KI = 0.5, µ0 = 4.1, KS = 200, K1 = 2.

u = Sin, la sortie y = S et la perturbation v = V q .

5) Refaire les questions 3 et 4 en considérant la biomasse X comme sortie du système.

A titre d’information, la représentation d’état (linéaire) du bio-réacteur est de la forme :

Vous aimerez peut-être aussi