1 Modèle de simulation
Dans cette partie, nous allons expliquer notre modèle de simulation. Notre
modèle est basé sur les paramètres suivants :
• V f0 : Volume initial du réservoir d’alimentation (en m3 ).
• Cw : Concentration de l’eau (en kg/m3 ).
• Sa : Surface de la membrane (en m2 ).
• dP : Différence de pression (en kg/m/h2 ).
• Aw : Coefficient de perméabilité de l’eau (en h/m).
• Bs : Perméabilité au soluté (en m/h).
• ϕ : Constante de passage du soluté (en m2 /h2 ).
• Cf0 : Concentration initiale dans le réservoir d’alimentation.
Les équations de notre modèle sont les suivantes :
Sa · Aw · dP
a1 = (1)
Cw
Sa · Aw · ϕ
a2 = (2)
Cw
Aw · dP
a3 =1+ (3)
Bs · Cw
Aw · ϕ
a4 = (4)
Bs · Cw
a5 = V f0 (5)
a6 = Cf0 (6)
Sa · Aw
a7 = (7)
Cw
dCf dCpav
Les équations pour calculer dt
et dt
sont données par :
2
2 ·Cf
dCf (a1 · Cf − a2 · Cf 2 + a3a−a 4 ·Cf
) · (1 − a3 −a14 ·Cf )
= 2 ·Cf ·t
dt a5 − a1 · t + a2 · Cf · t − a3a−a 4 ·Cf
a7 ·Cf Cf
dCpav (a 1 − a 2 · Cf + a3 −a4 ·Cf
) · a3 −a4 ·Cf
− Cpav
=
dt a · Cf −a 6
5 Cf −Cpav
1
L’équation pour calculer le flux d’eau (Jw ) est donnée par :
Jw = Aw · (∆P − ϕ · (Cf − Cp))
L’équation pour calculer la concentration du perméat à partir de la con-
centration de l’alimentation est donnée par :
Cf
Cp =
a3 − a4 · Cf
2 Suivi des paramètres au cours de la simu-
lation
Les paramètres à surveiller au cours de la simulation sont les suivants :
• L’évolution de la concentration d’alimentation dans le réservoir d’alimentation.
• L’évolution de la concentration du perméat à la sortie du module
d’osmose inverse.
• L’évolution de la concentration moyenne du produit dans le réservoir
de stockage.
• L’évolution du flux massique.
2
3 Résultats de la simulation
Les résultats de la simulation sont présentés sous forme de graphes ci-dessous.
Figure 1: Évolution de la concentration d’alimentation (réservoir
d’alimentation) Cf
Le graphe de l’évolution de la concentration d’alimentation (Cf) peut être
divisé en trois phases distinctes, chacune caractérisée par une dynamique
particulière.
Dans la première phase, qui s’étend jusqu’à environ t = 400 heures, la
concentration d’alimentation augmente légèrement et de manière presque
linéaire. Cela peut s’expliquer par le fait que le système est encore en train de
s’équilibrer, et la faible augmentation de la concentration reflète l’ajustement
progressif du réservoir d’alimentation.
La deuxième phase commence après les 400 heures et est marquée par
une accélération de l’augmentation de la concentration d’alimentation. Cette
hausse rapide peut être due à une augmentation de la perméabilité du module
d’osmose inverse ou à une diminution de l’efficacité de l’épuration.
Enfin, dans la troisième phase, la concentration d’alimentation atteint
un plateau ou une saturation. Cela pourrait résulter d’une stabilisation du
système ou de la saturation des capacités d’épuration du module d’osmose in-
verse. Cette phase de saturation indique que le système a atteint un équilibre
stable et que la concentration d’alimentation ne varie plus significativement.
3
Figure 2: Évolution de la concentration du perméat Cp et de la concentration
moyenne Cpav
L’évolution de la concentration du perméat (Cp) suit une tendance simi-
laire à celle de la concentration d’alimentation (Cf). On observe une augmen-
tation modérée jusqu’à t = 400 heures, puis une accélération de l’augmentation
au-delà de ce temps. En revanche, la concentration moyenne dans le réservoir
de stockage (Cpav) est nettement plus faible que la concentration du perméat.
4
Figure 3: Évolution du flux en fonction du temps
La diminution du flux en fonction du temps est expliquée par la diminu-
tion de la force motrice suite à la saturation progressive de la membrane.
5
Figure 4: Évolution de la concentration d’alimentation (réservoir
d’alimentation) Cf pour differentes valeurs de dp
Dans ce graphe, nous observons l’évolution de la concentration de l’alimentation
(Cf ) pour différentes valeurs de dp. Il est notable que l’augmentation de dp
entraı̂ne une augmentation de Cf dans la plage de valeurs entre 2.02 × 1013 et
6.02 × 1013 . Cependant, pour la valeur de 6.02 × 1013 , Cf augmente jusqu’à
atteindre un pic à 920 heures, après quoi elle devenir nulle. Ceci suggère que
dp ne devrait pas dépasser cette valeur.
6
Figure 5: Évolution de la concentration Cpav pour differentes valeurs de dp
7
Figure 6: Évolution de la concentration Jw pour differentes valeurs de dp