These
These
N° d’ordre :
Série :
THESE
Présentée pour l’Obtention du Diplôme de Doctorat en Science
cience en Génie des
Procédés
Par
Nadjiba BENMAKHLOUF
Année Universitaire
2020-2021
Université de Constantine 3
N° d’ordre :
Série :
THESE
Présentée pour l’Obtention du Diplôme de Doctorat en Science
cience en Génie des
Procédés
Par
Nadjiba BENMAKHLOUF
Année Universitaire
2020-2021
REMERCIEMENTS
Avant tout je remercie Dieu le tout puissant de m’avoir donné la santé et le courage,
pour terminer ce modeste travail.
Je dédie ce travail à :
Ma sœur ;
Mes frères ;
iv
Table des matières
v
Table des matières
vi
Liste des figures
vii
Liste des figures
viii
Liste des figures
ix
Liste des Tableaux
x
Liste des abréviations
xi
Nomenclature
NOMENCLATURE
xii
Nomenclature
xiii
RESUME
اﻷكسدة الحرارية المائية في الشروط فوق الحرجة هي تقنية لمعالجة النفايات تعمل على أكسدة المواد
العضوية باستخدام الماء في الشروط فوق الحرجة ،تم تطويرها كتقنية بديلة من أجل الحد من مخاطر التلوث الثانوي.
الهدف من هذه الدراسة هو تطوير أداة محاكاة في الحالة المستقرة لمفاعل أنبوبي لﻸكسدة في المياه فوق الحرجة
) ،(SCWOتابع ومتواجد بجامعة قادس في إسبانيا ،وقد تمت دراسة هذا المفاعل باستخدام نهج النمذجة ثنائية اﻷبعاد
بتطبيق النموذج المضطرب ،k-εباﻹضافة إلى صيغ ) IAPWSالرابطة الدولية لخصائص المياه والبخار( لحساب
الخصائص الحرارية والفيزيائية للماء في ظل الظروف فوق الحرجة ،تم افتراض أن ﺧصائص وسط التفاعل عبارة
عن ﺧليط من اﻷنواع الموجودة ،كما تم تحسين معامل نقل الحرارة.
النموذج ثنائي اﻷبعاد طور باستخدام برنامج " ،"COMSOL Multiphysics 5.2حيث سمح هذا اﻷﺧير بتحسين
نتائج النموذج أحادي اﻷبعاد بشكل كبير ،ولذلك تم استخدامه لنمذجة سلوك المفاعل لتحلل الملوث النموذجي :الفينول
والملوثات الحقيقية :زيت القطع Biocut35ونفايات زيت الزيتون.
النتائج التي تم الحصول عليها متوافقة بشكل جيد مع القيم التجريبية ،لذلك تم إستخدامها لدراسة حدودية لتأثير ظروف
التشغيل على أداء المفاعل باستخدام منهجية تصميم التجارب "." DOE
م ّكن البرنامج الذي تم تطويره ،ﻷول مرة ،من رسم منطقة أمان للمفاعل فوق الحرج لكل ملوث في ظروف تشغيل
مختلفة ،مما يضمن التشغيل ﻷقصى قدر من الكفاءة مع احترام حدود اﻷمان للمفاعل.
الكﻠمات المفتاحية :المياه فوق الحرجة ،اﻷكسدة ،النمذجة ،مفاعل أنبوبي ،منهجية تصميم التجارب ،منطقة اﻷمان.
Introduction générale
INTRODUCTION GENERALE
1
Introduction générale
2
Introduction générale
3
Chapitre I Synthèse Bibliographique
CHAPITRE I
SYNTHESE BIBLIOGRAPHIQUE
1.1. Introduction
4
Chapitre I Synthèse Bibliographique
Le principe de l'oxydation hydrothermale est basé sur l'oxydation des eaux usées
dans des conditions sous et supercritique de l'eau. L’oxydation dans les conditions sous
critiques est une oxydation en voie humide (Wet Air Oxidation : WAO), qui se produit à
une pression de 2 à 15 MPa, sous une température de 150 à 320 °C et l’oxydation dans les
conditions supercritiques est une oxydation supercritique (Super Critical Water Oxidation :
SCWO) se produit à des températures et des pressions supérieures au point critique de
l’eau (P_c = 22.1 MPa, T_c = 374 °C) [2],[3].
L’eau dans ces conditions supercritique spéciales devient un excellent milieu réactionnel,
et en raison de ses nouvelles propriétés, l’oxydation présente les avantages suivants: une
seule phase réactionnelle (pas de limitation de transfert de masse), des rendements de
réaction très élevés, une oxydation complète [4] et en plus, l'effluent liquide final n'est pas
toxique et peut être évacué directement dans l’environnement sans traitement [5], [6],
l'effluent gazeux est également propre, le monoxyde de carbone ne dépasse pas quelques
parties par million [7], aucun des composés NOx ne se forme en raison des températures
relativement basses des procédés SCWO (400-650 °C) comparé au procédé d'incinération
(900-1300 ° C) [8], [9] [10]. La réaction d'oxydation dans l'eau supercritique est hautement
exothermique, et la quantité d'énergie libérée par la réaction peut être récupérée et offre des
opportunités intéressantes pour la réalisation de procédé en autosuffisance énergétique
[11].
Le procédé d’oxydation hydrothermale est constitué de quatre étapes principales: la
préparation de la matière première pour l’alimentation et la pressurisation, la réaction, la
séparation du sel, la récupération de la chaleur et enfin la dépressurisation [1]. Il est
évident que certains éléments peuvent varier d’un système à l’autre, en fonction des
caractéristiques spécifiques des déchets à traiter et des objectifs spécifiques de chaque
traitement, et aussi du procédé lui même.
Dans l’eau supercritique, les propriétés varient fortement et sont très différentes de
celles de l’eau liquide normale ou de sa vapeur. La masse volumique change rapidement et
son évolution présente un grand gradient au point critique; elle prend une valeur
intermédiaire entre celle de l'eau à l'état liquide (1 g / cm3) et celle de l'eau à l'état gazeux
5
Chapitre I Synthèse Bibliographique
(0,001 g / cm3) [2]. Ce point est particulièrement important car il affecte la solubilité des
composés organiques et des gaz dans le milieu réactionnel [12].
L'eau supercritique a une viscosité d'un ordre de grandeur inférieure à celle de l'eau liquide
et, par conséquent, les coefficients de diffusion sont d'un ordre de grandeur supérieur [13].
La constante diélectrique statique chute de 78,5 F/m à 25 °C à 5 F/m prés de la région
critique jusqu’à l’intervalle entre 1 et 2 F/m à 450 °C [3]. Il en résulte que l'eau
supercritique (SCW) agit comme un gaz dense non polaire et que ses propriétés de
solvatation ressemblent à celles d'un solvant organique de faible polarité [6]. Elle montre
une miscibilité totale avec les gaz tels que l'azote, l'oxygène et le dioxyde de carbone [14],
en revanche, les sels inorganiques sont presque insolubles dans le SCW [13].
En ce qui concerne les propriétés thermiques, elles varient fortement au passage critique, la
valeur de la conductivité thermique diminue du fait de la rupture des liaisons hydrogène, ce
qui, du point de vue du processus, indique que les échanges thermiques seront moins
efficaces que pour l'eau dans des conditions ambiantes [2]. La capacité calorifique de l'eau,
est infinie au point critique, ceci est dû à la contribution énergétique nécessaire pour
réaliser la forte dilatation thermique subie par l'eau au voisinage du point critique [12].
Dans cet état, les propriétés physiques, thermiques et chimiques de l’eau sont évidemment
modifiées [15].
La combinaison de la solvatation et des propriétés physiques fait de l’eau supercritique un
milieu idéal pour l’oxydation des matières organiques. Lorsque les composés organiques et
l'oxygène sont dissous dans l'eau au-dessus du point critique, ils sont immédiatement mis
en contact moléculaire intime dans une seule phase homogène à haute température, sans
limitation de transport à l'interface et, pour des températures suffisamment élevées, la
cinétique est rapide et la réaction d'oxydation s'achève rapidement [16].
6
Chapitre I Synthèse Bibliographique
State) capable de décrire les systèmes SCWO [22]. Très peu des EoS cubiques classiques
donnent des résultats précis dans cette région. La corrélation de Peng et Robinson,
modifiée par Farrokh-Niae et al, a été utilisée par A. Leybros et al [23] avec une déviation
entre les valeurs calculées de toutes les propriétés de l’eau et les valeurs expérimentales de
2%.
Pour le calcul de la viscosité et de la conductivité thermique, Sengers et al [24] ont
proposé les équations de la conductivité thermique et de la viscosité en fonction de la
température et de la densité. Les corrélations proposées par Chung et al [25] peuvent être
utilisées pour calculer à la fois la viscosité et la conductivité thermique [7].
La formulation la plus utilisée pour calculer les propriétés de l'eau à l'état supercritique
avec précision est celle de l'association internationale des propriétés de l'eau et de la vapeur
(IAPWS-IF 97 : International Association of Properties of Water and Steam), pour le
calcul de la densité, la viscosité, la conductivité thermique, l’enthalpie, la capacité
calorifique et la constante diélectrique. Cette formulation calcule ces propriétés en fonction
de la température et de la pression, elle peut être appliquée pour des pressions allant jusqu'à
100 MPa et des températures allant jusqu'à 2273,15 K [26]. Elle a été utilisée dans les
conditions supercritiques par A. Fourcault et al [18], J. Mercadier et al [20], Jan A. M.
Withag et al [27], J. P. Serin et al [28], N. Benmakhlouf et al [29] et elle a donné de bons
résultats.
Afin d'expliquer la variation des propriétés de l'eau au point critique, nous avions
développé un programme sur le logiciel Mathcad [1] ; à partir de la formulation de
l'AIPWS. Les variations des propriétés : masse volumique, viscosité dynamique, capacité
calorifique, constante diélectrique, conductivité thermique et enthalpie, au voisinage du
point critique en fonction de la température à des pressions de 22.1 MPa, 25 MPa et 30
MPa sont présentées sur la Figure (1. 1) suivante :
7
Chapitre I Synthèse Bibliographique
1000
0.0016
0.0004
200
0 0.0000
0 200 400 600 800 0 200 400 600 800
100
80
Constante diélectrique(F.m-1)
60
P=22.1 MPa
P=25 MPa
P=30 MPa
40
20
0
0 200 400 600 800
Température (°C)
0.8
4000
Conductivité thermique (W.m .K )
-1
-1
0.6
3000
Enthalpie (KJ/Kg)
0.2 1000
0.0 0
0 200 400 600 800 0 200 400 600 800
8
Chapitre I Synthèse Bibliographique
9
Chapitre I Synthèse Bibliographique
10
Chapitre I Synthèse Bibliographique
destruction totale de l'azote. Les températures optimales étaient de 400 °C pour le carbone
organique total et de 500 °C pour l'azote total, et le temps optimal et l'excès d'oxygène
étaient respectivement de 30 min et 100%.
En général, certains des composés traités sont efficacement oxydés dans le SCW, tandis
que d'autres sont plus réfractaires dans ce processus, sauf si des conditions sévères sont
utilisées pour atteindre un niveau élevé de destruction. D’autre molécules et déchets traités
par SCWO avec leurs taux de destruction sont présentés sur L’annexe (A).
11
Chapitre I Synthèse Bibliographique
avec restes d’eau), une phase aqueuse et une phase solide [51]. En ajoutant un courant
d’eau pour laver la phase huileuse, après en passant dans une nouvelle centrifugeuse,
générant une eau de lavage qui ressemble à une dilution : l'Alpechin classique.
Les déchets d’huile d’olive Alpechin, contiennent également de fortes concentrations de
composés phytotoxiques et inhibiteurs microbiens, tels que les composés phénoliques et les
acides gras à longue chaîne [52] , donc un traitement biologique sera difficile et n'est pas
l'option la plus appropriée, aussi la concentration de ce type de déchet est trop faible pour
une incinération économique, d’où la recommandation d’un traitement d'oxydation surtout
à l’eau supercritique [47]. Le déchet d’huile d’olive utilisé dans ce travail provient de la
coopérative viticole La Purísima, située dans la province de Córdoba en Espagne.
D’un autre coté, l'oxygène peut réagir directement avec les substances organiques pour
former un radical organique (R•) et des radicaux (HOO•), cette réaction se produit à des
températures élevées où l'oxygène réagit facilement avec les substances organiques,
ensuite le radical organique réactif (R•) réagit avec l'oxygène pour générer un radical
12
Chapitre I Synthèse Bibliographique
13
Chapitre I Synthèse Bibliographique
14
Chapitre I Synthèse Bibliographique
J. Shànchez -Oneto
Oneto et al [59] ont proposé un modèle cinétique pour les deux fluides
de coupe Servol et Biocut , dans une plage de température comprise entre 400 °C et 5500 °C
et une pression de 25 MPa. Dans le cas de Servol, un pseudo
pseudo-premier
premier modèle cinétique a
été utilisé. Par ailleurs,, un modèle mathématique à deux paramètres comportant deux
étapes, une réaction rapide suivie d'une réaction lente, était nécessaire pour décrire la
cinétique de Biocut et pour calculer les constantes cinétiques. Les deux modèles sont
présentés dans le Tableau (11. 1).
Ce mécanisme a été utilisé, avec succès, pour modéliser l'oxydation SCWO de l'huile de
transformateur minérale « mineral transformer oil » par V. Marulanda [60], dont la
composition est très similaire à celle de
des déchets d'huile de coupe.. Il a aussi été utilisé par
V. Vadillo et al [43] pour la simulation de Biocut dans le cas des concentrations élevées.
Le même groupe, J. Shànchez -Oneto et al [42], a proposé aussi un modèle cinétique global
pour les fluides de coupe, dans la m
même plage de température et la même pression,
pression ce
modèle dépendait de la concentration en oxygène
oxygène, dont la réaction était de l’ordre de 0.58
[42]. Ce mécanisme a aussi été utilisé avec succès pour modéliser l'oxydation SCWO de
l'huile de coupe « Biocut » par J.M. Benjumea et al [61].
La cinétique des réactions de la matière organique dans l’eau supercritique est
fréquemment considérée comme une cinétique du premier ordre ou pseudo-premier
pseudo ordre,
onn prend en considération la concentr
concentration des composés organiques [7].
[7] Le phénol était le
polluant le plus utilisé par les chercheurs pour la détermination de ses paramètres
cinétiques (Gopalan et Savage [62], Krajnc et Levec [63], Koo et al [64],
[64] Portela et al [65],
YoshitoOshima et al [66].. D. Mateo et al [67], TD Thornton et al [33]).
Le Tableau (1. 1) suivant présente quelques exemples d'études cinétiques réalisées par
plusieurs groupes de recher
recherche.
15
Chapitre I Synthèse Bibliographique
16
Chapitre I Synthèse Bibliographique
35 250
30
200
P=22.1 MPa
C a p ac ité c alo rifiq u e (J/(g .K ))
P= 22.1 MPa
25 P=24 MPa
P= 24 MPa
P=26 MPa
N o m bre de P randtl
P= 26 MPa
20 P=30 MPa 150
P= 30 MPa
15
100
10
50
5
0 0
320 360 400 440 320 360 400 440
Température (°C) Température (°C)
(a) (b)
Figure 1. 4: Capacité calorifique et le nombre du Prandtl en fonction de la température et de
la pression pour SCW (notre programme [1]).
Il y a eu plusieurs travaux sur les méthodes expérimentales ainsi que les méthodes
numériques réalisés dans le but d’estimer les coefficients de transfert de chaleur dans des
17
Chapitre I Synthèse Bibliographique
conditions supercritiques [80], l’estimation de ces paramètres a été utilisé beaucoup plus
dans le domaine « supercritical water-cooled nuclear reactor concepts » [81], [82].
Plusieurs corrélations sont proposées pour le calcul du coefficient de transfert de chaleur
interne en SCW, la corrélation de Dittus-Boelter (1930) est la plus largement utilisée à des
pressions sous-critiques pour la convection forcée [81]:
Nu = C. Re Pr . F (1. 6)
18
Chapitre I Synthèse Bibliographique
Figure 1. 5: Comparaison des valeurs de coefficient de transfert interne calculées avec des
corrélations proposées selon S. Mokrya et al [81].
Sur la Figure (1. 5), les valeurs du coefficient de transfert interne sont calculées à partir
des corrélations de : Dittus et Boelter (1930), Bioshop et al (1964), Swenson et al. (1965)
and Jackson (2002). Le Tableau (1. 2) montre ces corrélations, vérifiées pour le cas de
transfert de chaleur dans l’eau supercritique dans des conditions opératoires différentes.
Dans ce qui suit, l'indice « b » signifie que la propriété est calculée à la température
moyenne de l'écoulement, et l'indice « w » signifie que la propriété est calculée à la
température de l'écoulement à surface interne de la paroi.
Pour choisir la meilleure corrélation, il faut comparer la prédiction numérique avec les
données expérimentales obtenues dans des conditions similaires aux conditions de
fonctionnement [83].
19
Chapitre I Synthèse Bibliographique
flux de flux de
Temperatu Pression
masse chaleur Corrélation
re (°C) (MPa)
(kg/m2.s) (MW/m2)
. .
<=350 23.3-26.3 / 0.58-4.65 Nu = 0.1 Re Pr
Domin
(1963) μ
. .
250-350 23.3-26.3 / 0.58-4.65 Nu = 0.036Re Pr ∙
μ
.
Bioshop ρ D
Nu = 0.0069Re . Pr .
. 1 + 2.4.
et al 282–527 22.6 -27.5 651–3662 0.31–3.46 ρ L
(1964)
.
Swenson . .
ρ
Nu = 0.00459Re Pr .
et al. 75–576 22.8-41.4 542–2150 / ρ
(1965)
.
Jackson ρ Cp
/ / / / Nu = 0.0183Re .
Pr .
. .
(2002) ρ Cp
Nu
Feng et al 600-1200 ρ
. .
/ 23-30 0.1-0.6 µ
(2005) = 0.02269Re .
Pr .
. .
ρ µ
Sara ρ
.
. .
Mokry et 320-406 24 200-1500 0.07-1.25 Nu = 0.0061Re Pr .
ρ
al (2011)
20
Chapitre I Synthèse Bibliographique
Ce réacteur a prouvé son efficacité pour la destruction des composés organiques simple et
complexe comme : l’acide acétique, le méthanol, le phénol [67], les matériaux
énergétique : RDX (Héxogène), HMX (Octogène) et CL20 (2, 4, 6, 8, 10, 12-hexanitro- 2,
4, 6, 8, 10, 12 –hexaazaisowurtzitane) [85].
Un réacteur de ce type a été développé à l’échelle semi-industrielle par le HOO (Option
d’Oxydation Hydrothermale) en France.
21
Chapitre I Synthèse Bibliographique
interne est d'environ 12,32 mm. Ce pilote a été utilisé de manière satisfaisante lors de
plusieurs études antérieures sur le SCWO pour la destruction de nombreux polluants; il a
démontré son efficacité pour le phénol [18], [34], huile de coupe [43], [44], éthanol,
déchets de l’huile d’olive [18], eaux usées industrielles inflammables [86].
Initialement ce réacteur a été utilisé avec une seule injection d’oxygène, et afin d’améliorer
l'efficacité de l'oxydation de l'eau supercritique, en particulier l'efficacité de la destruction
des polluants et du contrôle de la température dans le réacteur, des injecteurs d’eau de
refroidissement et d’oxydants séparés ont été installés [44], les deux systèmes sont
capables de refroidir le réacteur. La Figure (1. 7) suivante présente le pilote complet :
Il s’agit d’un réacteur de forme Y, comportant une section inclinée à gauche et une
section verticale à droite, la section gauche est la zone de la réaction et la section droite est
la zone de séparation.
La section gauche est un réacteur à trois parois constitué d'un mur extérieur
supportant la pression, d'un mur de distribution intermédiaire et d'un mur poreux intérieur
en acier inoxydable 316L.
L’air est injecté dans le réacteur par la seule entrée située du côté de la paroi sous
pression, il est distribué par des trous situés sur la paroi de distribution et finit par traverser
la paroi poreuse dans la zone de réaction. L'air agit en tant que fluide transpirant qui forme
22
Chapitre I Synthèse Bibliographique
un film de gaz protecteur sur la surface intérieure de la paroi poreuse pour empêcher la
corrosion et le dépôt de sel, il agit également en tant qu'oxydant pour le SCWO, de plus,
une conduite d’alimentation indépendante en oxygène pur (O2) est ajoutée en tant que
source principale d’oxydant.
Les produits de la partie gauche du réacteur sont séparés dans la partie droite, et
cela dans des conditions idéales, en un écoulement ascendant propre (eau et gaz sans
particules solides) et une saumure à écoulement descendant (sels et particules solides avec
eau), qui sont libérés du réacteur par des sorties en haut et en bas, respectivement.
Les deux sections gauche et droite du réacteur sont équipées avec deux fours de
chauffage électriques indépendants, les fours sont conçus pour fournir de la chaleur au
démarrage et maintenir la température dans le réacteur pendant les expériences [87], ce
réacteur est présenté dans la Figure (1. 8) suivante :
23
Chapitre I Synthèse Bibliographique
Dans ce même laboratoire, un autre pilote SCWO a été développé, il s’agit d’un réacteur
tubulaire en acier inoxydable de 2.5 m, et de diamètre interne de 5.2 mm. L'eau et les
composés organiques sont amenés par l'intermédiaire de deux pompes, avant d'être
préchauffés à 673 K à partir des colliers chauffants en céramiques entourent le tube. Le
mélange des deux lignes eau/déchets et oxydant s'effectue dans le réacteur maintenu à 773
K par un bain de sable fluidisé. Après, un refroidisseur abaisse la température du mélange
réactionnel à une température comprise entre 293 et 303 K, ensuite, l'effluent est
dépressurisé à travers un déverseur. Un schéma de principe de cette installation est
présenté sur la Figure (1. 10).
24
Chapitre I Synthèse Bibliographique
25
Chapitre I Synthèse Bibliographique
Toutefois, dans la technologie de conception des réacteurs SCWO, des progrès et des
développements supplémentaires sont nécessaires pour le contrôle de la corrosion des
réacteurs, l'amélioration de la durabilité, l'amélioration de l'efficacité de la réaction et la
facilité de maintenance.
26
Chapitre I Synthèse Bibliographique
récupération d'énergie dégagée par les réactions. Une grande variété de logiciels de CFD
est disponible maintenant pour effectuer la modélisation numérique d'un processus
d’oxydation dans l'eau supercritique. Parmi les plus utilisés par les chercheurs et
l'industrie, on retrouve le COMSOL [60],[27] et le FLUENT [19].
Les études réalisées dans le coté simulation, étaient uni, bi et tridimensionnelles, basées sur
des hypothèses simplificatrices relatives en premier lieu, les propriétés du milieu
réactionnel (où la majorité suppose que les propriétés thermo-physiques du milieu
réactionnel sont celles de l'eau dans l'état supercritique), l’estimation du coefficient de
transfert de chaleur, le polluant organique à éliminer, et aussi le régime d'écoulement dans
les réacteurs. L’état stationnaire était le plus utilisé par les modèles étudiés et très peu de
travaux se sont intéressés à l’état transitoire.
27
Chapitre I Synthèse Bibliographique
28
Chapitre I Synthèse Bibliographique
D’un autre coté, J.R. Portela et al [91], avaient étudié l'oxydation du méthanol dans un
réacteur multi-injection d'oxygène développé au niveau de l'ICMCB (Figure (1. 6)).
Différentes expériences avaient été menées dans des conditions quasi adiabatiques, avec
une température d'alimentation de 430 ºC à l'entrée du réacteur avec une concentration du
méthanol de 6.67 g/L. La simulation était réalisée avec le logiciel PROSIM plus, pour
obtenir les bilans de matière et d’énergie, les propriétés thermo-physiques du milieu
réactionnel étaient calculées à partir du mélange, les paramètres cinétiques du méthanol
utilisés étaient ceux trouvés expérimentalement par D. Mateos et al [67] dans le même
réacteur. Afin de tester la sensibilité du modèle à l'enthalpie de la réaction d'oxydation
hydrothermale du méthanol, trois simulations, avec ∆H égal à : -635 kJ /mol, -685 kJ /mol
et -735 kJ /mol ont été effectuées et comparées avec le profil de la température comme
montré sur la Figure (I. 14).
En ce qui concerne le réacteur à paroi froide, M.J. Cocero et al [56] avaient utilisé un
modèle pour le réacteur développé au niveau de l’université de Valladolid (Figure (1. 11)),
pour l’alcool isopropylique et le dimethylsulfoside, ce modèle basé sur un bilan de masse
et d’énergie, et le modèle cinétique généralisé proposé par Li et al [53] présenté sur la
Figure (1. 2). Le réacteur est divisé en trois parties : la partie annulaire externe (Coque
externe du réacteur), bobine intermédiaire (Intermediate coil), et chambre de réaction. Les
29
Chapitre I Synthèse Bibliographique
températures obtenues par la simulation numérique sont en accord avec les mesures
expérimentales comme montré sur la Figure (1. 15).
Figure 1. 15: Profil de la température avec 8wt% alcool isopropylique et 5000 ppm
dimethylsulfoside (Tmax=780 °C, tr=29 s), M.J. Cocero et al [56].
Dans le cas d’un régime transitoire, J.M.Benjumea et al [61] avaient développé un modèle
mathématique unidimensionnel résolu en utilisant Matlab, pour le réacteur tubulaire
développé à l’université de Cádiz pour les fluides de coupe. Le modèle avait basé sur des
bilans de conservation de masse, de l’énergie et de quantité de mouvement et la variation
des différentes espèces. Les paramètres cinétiques utilisés étaient ceux de J. Shánchez –
Oneto et al [42], et les propriétés utilisées étaient celles du mélange. Les résultats du
programme Matlab étaient en bon accord avec les résultats expérimentaux comme montré
sur la Figure (1. 16). Ce modèle montre aussi la variation temporelle en fonction des
changements des paramètres de fonctionnement tels que la température, la concentration
ou les pertes de chaleur le long du réacteur.
30
Chapitre I Synthèse Bibliographique
Dans ces cas, les travaux ont fait recours aux logiciels de simulation pour la
résolution des modèles à utiliser. Les modèles de turbulence ont aussi été appliqués dans le
cadre des simulations 2D et 3D [12], le tenseur des contraintes devrait inclure le tenseur
des contraintes de Reynolds et un modèle approprié pour le calculer [18]. Le modèle le
plus couramment utilisé est le modèle k-ɛ, car il a été utilisé avec succès pour simuler
l’écoulement dans l'eau supercritique [80] .
S. Moussiere et al [12] avaient établi et validé un modèle basé sur la géométrie
bidimensionnelle et tridimensionnelle, en utilisant le simulateur CFD Fluent, pour le flux
réactif monophasique dans le réacteur qui a été développé au CEA (Figure (1. 10)), et cela
pour le Dodécane et le méthanol. Deux modèles de turbulence ont été choisis pour décrire
les phénomènes de turbulence: le modèle k-ω et le modèle k-ε. Les résultats montraient
que le modèle k-ε offre une meilleure description de l'écoulement comparé au modèle k-ω
dans le cas sans réaction chimique, tandis que dans le cas avec réaction chimique, le
modèle k-ω offrait une meilleure description de l’écoulement, et le profil de la température
expérimental dans ce cas était très proche du profil prédit, comme montré sur la Figure (1.
17).
31
Chapitre I Synthèse Bibliographique
(a) (b)
Figure 1. 17: Simulation du profil de la température avec les modèles k-ω et k-ε, dans le
cas sans réaction chimique (a) et avec réaction chimique (b) [19].
La comparaison des modèles 2D et 3D est appliquée sur le tube interne de réacteur double
enveloppe présenté sur la Figure (1.9), cette comparaison montre que la simulation 3D
n'apporte pas de précision supplémentaire sur la prédiction des transferts thermiques ; bien
que les temps de calcul soient largement supérieurs. Ainsi, les profils de température
calculés à partir des modèles 3D sont identiques aux profils calculés par le modèle 2D et
dans tous les cas les températures sont sur-estimées par rapport aux températures
expérimentales. La Figure (1. 18) présente les profils de température avec les deux
modèles 2D et 3D.
Figure 1. 18: Profils de température prédits le long du tube interne du réacteur par le
modèle 2D et le modèle 3D [12].
32
Chapitre I Synthèse Bibliographique
Dans un autre travail, V. Marulanda [60] avait utilisé le modèle de turbulence k-ε sous le
simulateur CFD Comsol Multiphysics 3.3, avec une représentation 2D, où les propriétés du
milieu réactionnel étaient supposées celles de l’eau pure. L’huile minérale de
transformateur, fortement contaminés par les PCB (polychlorinated biphenyls) était
éliminée avec le peroxyde d'hydrogène comme agent oxydant, les paramètres cinétiques
utilisés étaient aussi ceux trouvés par J. Shánchez–Oneto et al [59]. Le modèle proposé
avait pour objectif de trouver les conditions opératoires optimales pour atteindre un
fonctionnement idéal du pilote dont le réacteur est tubulaire. Les simulations avaient
montré que le fonctionnement du pilote proposé n’était pas possible ; en raison de la
température élevée à l’intérieur du réacteur qui résultait des réactions d'oxydation
fortement exothermique et qui pouvait atteindre 1000 K pour une température à l’entrée de
800 K, comme montré sur la Figure (1. 19). Ils recommandaient, donc d'utiliser une
température d'entrée inférieure afin d'éviter des conditions de fonctionnement de
température dangereuses, tout en maintenant les objectifs de destruction de la charge
organique et de PCB identiques.
Figure 1. 19: Le profil de température avec le taux de conversion de PCB pour différente
valeurs de température à l’entrée, V. Marulanda [60].
Dans le cas d’une simulation 3D, N. Zhou et al [92] avaient réalisé une représentation dans
un réacteur réservoir de l’oxydation supercritique du méthanol ; avec l’utilisation du code
commercial CFD-ACE+. Le modèle standard k-ɛ a été utilisé pour simuler la turbulence de
l'écoulement, le réacteur est considéré totalement adiabatique, et parfaitement agité ; et les
propriétés physiques du fluide étaient ceux du mélange et obtenues à l'aide du logiciel
ASPEN PLUS.
33
Chapitre I Synthèse Bibliographique
Figure 1. 20: Résultats CFD de la conversion du méthanol avec Simulation 3D, N. Zhou et
al [92].
34
Chapitre I Synthèse Bibliographique
L'idée de l’eau supercritique pour le traitement des eaux usées et de déchets est
apparue à la fin des années 1970, et l'étude s'est accélérée durant les années 1980 [98], où
elle était devenu un support intéressant pour un grand nombre d'applications,
notamment dans le domaine du traitement des eaux usées et les déchets aqueux en général.;
35
Chapitre I Synthèse Bibliographique
36
Chapitre I Synthèse Bibliographique
Mitsubishi Heavy Industries (MHI), qui avait mis en œuvre en 2005 son premier procédé à
grande échelle au Japon pour le traitement des déchets contenant des PCB, d’une capacité
de 6 tonnes / jour.
La société française Hydrothermale Oxidation Option (HOO), en collaboration avec le
CNRS, avaient développé deux installations de démonstration industrielle (toutes les deux
100 kg / h) destinées à être commercialisées dans le traitement des déchets industriels.
En France, Innoveox avait exploité en 2011 une usine SCWO pour les déchets industriels,
notamment de pétrochimie et de la pharmacie.
Le tableau (I. 3) suivant donne ces installations commerciales SCWO à grande échelle
avec d’autre, [100], [99], [51], [54].
37
Chapitre I Synthèse Bibliographique
SuperWater
Solutions Orlando, Floride Boues d'épuration 2013-présent
38
Chapitre I Synthèse Bibliographique
1.13.1. La corrosion
H O ⇔ H + OH (1. 7)
Cette réaction a un caractère endothermique et, par conséquent, à des températures élevées,
l'équilibre se déplace vers le côté droit, c'est-à-dire qu'à une pression de 25 MPa et à des
températures élevées, les concentrations de H+ et de OH- sont d'environ trois ordres de
grandeur au-dessus des valeurs dans l’eau ambiante, de sorte que l’eau puisse être
considérée à la fois acide et alcaline [100].
En fait, la corrosion des réacteurs est la principale raison qui retarde la commercialisation
du procédé SCWO [51]. Plusieurs auteurs avaient réalisé des études pour tenter de
résoudre ce problème de corrosion. M.Hayward et al [104] avait étudié la corrosion d'un
réacteur en acier inoxydable pour l'oxydation de l'eau supercritique; l'effet de l'oxygénation
sur la corrosion a été examiné en utilisant de l'eau dégazée, de l'eau oxygénée et une
solution de H2O2 à 3%. Pour l'utilisation d'eau dégazée et d'eau oxygénée, aucune variation
39
Chapitre I Synthèse Bibliographique
significative du degré de corrosion n'a été observée, mais lorsque H2O2 était ajouté dans le
flux, la corrosion augmentait considérablement, à des températures proches de la
température critique, et le taux de corrosion moyen le plus élevé était de 0,125 mm / an à
375 ° C. Ils ont prouvé aussi que l'exposition du réacteur à l'air libre entre les expériences
est un facteur contribuant à la corrosion du matériau du réacteur. D'autre part, BR Foy
[105] avaient étudié la corrosion du réacteur en titane et le peroxyde d'hydrogène est
l'agent oxydant. La vitesse de corrosion la plus élevée observée était de 14 mm / an.
Yuzhen Wang et al [106], avaient réalisé la comparaison entre la corrosion de l'acier
inoxydable 316L et de l'alliage 625 « alloy 625 », la réaction a été testée dans un procédé
d'oxydation à l'eau supercritique pour le traitement de l'eau de cokéfaction «coking waste
water » , ils avaient constaté une corrosion par piqûre à la surface d’acier inoxydable 316L
avec une vitesse de corrosion de 0,334 mm / an. Bien qu’aucune corrosion remarquable
n’ait été constatée à la surface de l’alliage 625, sa vitesse de corrosion était de 0,052
mm/an, d’après cette étude, la résistance à la corrosion de l'alliage 625 est beaucoup plus
élevée que celle de l'acier inoxydable 316L dans les conditions supercritiques.
En ce qui concerne l'économie et la résistance à la corrosion, les aciers inoxydables et les
alliages à base de nickel ont d’avantage attiré l'attention des chercheurs que les autres
matériaux [107]. Jusqu'à présent, aucun alliage spécifique n'a été développé pour une
utilisation dans le SCW. Les principaux matériaux métalliques couramment utilisés dans
les systèmes des fluides supercritiques comprennent les aciers doux, les aciers inoxydables
et les alliages à base de nickel [108].
En raison des propriétés uniques de l'eau dans des conditions supercritiques, les
substances organiques et l'oxygène ont une solubilité complète, tandis que les espèces
polaires telles que les sels inorganiques sont insolubles et précipitent [5]. Ce phénomène
est dû essentiellement à la faible constante diélectrique de l'eau supercritique et peut aller
jusqu'à l’obturation du réacteur [54]. S'ils ne sont pas contrôlés, les sels déposés peuvent
bloquer le réacteur, l'échangeur de chaleur, les lignes de transport ou d'autres zones, et
même provoquer une explosion, entraînant la défaillance où l'arrêt coûteux des systèmes.
Le dépôt de sel à long terme provoque également un encrassement et une corrosion plus
rapide du réacteur.
40
Chapitre I Synthèse Bibliographique
Par conséquent, mis à part les méthodes de nettoyage conventionnelles qui sont coûteuses,
des conceptions spéciales de réacteurs et/ou des approches opérationnelles pratiques sont
nécessaires pour surmonter le défi du dépôt de sel [99]. Chematur dans ses applications
commerciales avait procédé à augmenter la vitesse du fluide pour maintenir les particules
en suspension et éviter le dépôt [82]. L’augmentation de la pression de fonctionnement
permet d’augmenter la densité de la solution, et par la suite la solubilisation des sels. Le
problème de cette solution est que la solubilité de la couche d'oxyde protectrice augmente
également, ce qui favorise la corrosion [100]. P. Kritzer et E. Dinjus [109] avaient souligné
que la densité de la solution dans le réacteur devait être inférieure à 200 kg / m3 afin de
minimiser la corrosion, et ils ont indiqué que la meilleure solution pour éviter la
précipitation de sel est de réduire la quantité de sel présente dans l’alimentation, grâce à
des dispositifs de séparation solide-fluide. Ces dispositifs peuvent être utilisés avant ou
après l'étape de réaction.
1.14. Conclusion
41
Chapitre I Synthèse Bibliographique
conceptions de réacteur capables d'assurer une bonne protection des matériaux, évitant
ainsi la corrosion et le dépôt de sel tout en maximisant la récupération d'énergie.
42
Chapitre II Etude expérimentale et modélisation de l’oxydation hydrothermale
CHAPITRE II
HYDROTHERMALE
2. 1. Introduction
Les simulations numériques ont pris une importance grandissante ces dernières
années, tant pour l’étude des phénomènes physiques que pour le calcul d'écoulements
industriels. Moins coûteuses que les expériences sur les bancs d'essai, elles permettent de
modifier facilement les paramètres de l'écoulement et la géométrie. Pour cela, il est
nécessaire de disposer d’un outil numérique performant, qui s'adapte facilement aux
géométries complexes rencontrées dans l'industrie, et qui prenne en compte les
comportements complexes d'un fluide en situation réelle (turbulence, compressibilité,
ect...). Il faut toutefois noter que l'utilisation d'un modèle numérique" n'exclut pas
l'expérience, mais constitue une approche complémentaire permettant de remplacer à
moindre coût de trop nombreux essais.
La simulation numérique à l'aide de Computational Fluid Dynamics (CFD) est considérée
comme une approche utile pour le développement fiable des corrélations de transfert de
chaleur ; et de fournir une meilleure compréhension des mécanismes réactionnels et de
transfert de chaleur dans des géométries et des conditions opératoires pour les réacteurs
d’oxydation dans l’eau supercritique.
CFD est un ensemble de méthodes numériques appliquées pour obtenir des solutions
approchées de problèmes de dynamique des fluides, et de transfert de chaleur dans l'espace
multidimensionnel. Selon cette définition de CFD, elle est un moyen d'appliquer les
méthodes d'une discipline (analyse numérique) à un autre (de la dynamique des fluides et
de transfert de chaleur).
Ce chapitre illustrera la partie expérimentale et de simulation numérique de cette
recherche. Concernant la partie expérimentale, cela inclut une description complète des
équipements utilisés et du fonctionnement du système dans le cas d’oxydation dans des
43
Chapitre II Etude expérimentale et modélisation de l’oxydation hydrothermale
Cette section décrit les équipements dans lesquels le travail expérimental a été
développé, à savoir, l'installation de l'oxydation dans les conditions supercritiques à
l'échelle pilote appartenant au groupe de recherche du département de génie chimique et
technologie alimentaire de la faculté des sciences, à l’université de Cádiz en Espagne.
L’installation pilote est conçue pour fonctionner en continu avec un débit nominal de 25
L/h, elle est constituée essentiellement d’un réacteur tubulaire et des éléments
supplémentaires, comprenant des réservoirs de préparation de l’alimentation, des pompes à
impulsion, un compresseur à air, des échangeurs de chaleur pour préchauffer les flux
d'alimentation, le refroidisseur de l'effluent final et des vannes de régulation de pression.
La Figure (2. 1) représente une photo du pilote, et sur la Figure (2. 2) un schéma
représentatif de la configuration du pilote pour faciliter l'explication ultérieure du
fonctionnement de l'équipement.
44
Chapitre II Etude expérimentale et modélisation de l’oxydation hydrothermale
45
Chapitre II Etude expérimentale et modélisation de l’oxydation hydrothermale
maintenue à 250 bars à l’aide d’une soupape de régulation. Ensuite, l'alimentation est
introduite dans un échangeur de chaleur à tubes concentriques à contre-courant, pour la
préchauffer. Cependant, pour démarrer l’installation, il existe un système de résistances
enroulées sur le tube d’alimentation. D’un autre coté, l'air utilisé comme oxydant, est aussi
conduit par un compresseur (voir Figure (2. 4)) vers deux bouteilles où il est stocké, et
après il conduit vers un autre échangeur de chaleur également en tubes concentriques à
contre-courant, pour le préchauffage. Les deux fluides sont mélangés à l'entrée du réacteur
tubulaire et la température du mélange recommandée (environ 400 °C) pour commencer
l'oxydation du résidu dans des conditions supercritiques. Comme cette réaction est très
exothermique, la chaleur produite peut être utilisée pour préchauffer partiellement les
alimentations en liquide et en air. Ce courant est alors utilisé pour alimenter les échangeurs
de chaleur destinés au préchauffage des courants d'entrée. A la sortie des échangeurs,
l’effluent reste à une température élevée (≈200 °C), un refroidissement est donc nécessaire
pour atteindre des conditions acceptables de température (T <40°C) avant son entrée dans
le séparateur final gaz-liquide. Pour cela, on utilise un refroidisseur à plusieurs étages, qui
utilise l'eau du réseau comme réfrigérant.
Une fois refroidi, le courant est dépressurisé à l'aide de l'électrovanne de régulation de
pression, passe dans une conduite dirigée vers le séparateur gaz-liquide, où le courant est
séparé en une phase gazeuse et une phase liquide.
Un autre réservoir contient juste l’eau, est utilisé pour le nettoyage du système après les
expériences.
46
Chapitre II Etude expérimentale et modélisation de l’oxydation hydrothermale
a) Le Réacteur
Il s’agit d’un réacteur tubulaire constitué de trois sections de tubes en acier inoxydable
(AISI 316L), de 9.52 m de longueur, de 12.32 mm de diamètre intérieur et 19.05 mm de
diamètre extérieur. Les sections ont une longueur de 2920, 2960, 3000 mm,
respectivement, elles sont placées horizontalement et sont reliées par des pinces « Hoke »
pour les applications à haute pression, avec des sections verticales de longueurs 69,3 et
69,9 mm. Au total, le réacteur a un volume de 1230 cm3, et le long du réacteur, il y a 24
thermocouples soudés à des colliers métalliques et placés sur le tube, comme indiqué dans
le Tableau (2. 1) suivant :
Tableau 2. 1 : Les thermocouples placés sur le réacteur.
Thermocouple Longueur du réacteur (m) Thermocouple Longueur du réacteur (m)
1 0.05 13 3.63
2 0.22 14 4.65
3 0.32 15 4.82
4 0.48 16 5.04
5 0.64 17 6.16
6 1.13 18 6.26
7 1.63 19 6.36
8 2.11 20 7.39
9 2.62 21 7.79
10 3.12 22 9.01
11 3.3 23 9.27
12 3.43 24 9.52
47
Chapitre II Etude expérimentale et modélisation de l’oxydation hydrothermale
L’isolation thermique
Comme indiqué précédemment, l'installation pilote est conçue pour que le processus soit
autonome du point de vue énergétique, de sorte que la chaleur dégagée lors du processus
d'oxydation soit suffisante pour préchauffer l'alimentation, en liquide et en gaz, à la
température requise à l’entrée du réacteur (380-400 °C). Pour minimiser les pertes
thermiques, il est essentiel d'isoler thermiquement à la fois le réacteur et les échangeurs de
chaleur des courants d'entrée, afin de maximiser la chaleur dégagée. Le matériau isolant
utilisé est une couverture thermique de la maison de commerce Unifrax, plus précisément,
le Fiberfrax Durablanket S avec une composition de 43-47% d’AL2O3 et de 53-57% de
SiO2. Le Tableau (2. 2) suivant indique les épaisseurs d'isolation utilisées dans chacun des
équipements :
Tableau 2. 2 : Les épaisseurs d'isolation du réacteur, échangeurs et préchauffeurs
électrique.
Unité Epaisseurs (mm)
Réacteur 45.73
Échangeur de liquide 32.99
Échangeur de gaz 32.78
Préchauffeur électrique 45.73
48
Chapitre II Etude expérimentale et modélisation de l’oxydation hydrothermale
e) Refroidisseur multi-étages
Le courant de l’effluent de sortie du réacteur est utilisé pour alimenter les échangeurs de
chaleur destinés au préchauffage des courants d'entrée, à la sortie des échangeurs, l’effluent
reste à une température élevée (≈200 °C), valeur bien supérieure à celle nécessaire à son
rejet, pour refroidir ledit flux du produit à une température de 40 °C, l'installation dispose
d'un dernier échangeur. Ce refroidisseur est un échangeur à plusieurs étages, et présentant
des dimensions extérieures très compactes, il est alimenté en eau de robinet.
49
Chapitre II Etude expérimentale et modélisation de l’oxydation hydrothermale
Pour mesurer la composition de la phase gazeuse produite lors des différents tests réalisés
dans l’installation pilote, il existe un analyseur portable de CO, CO2, O2, H2 et CH4,
modèle 79000 P-5, fabriqué par NOVA (USA). Pour effectuer les différentes mesures, il
est muni d’un détecteur infrarouge de CO, CO2 et CH4, capteur électrochimique O2, cellule
de thermo-conductivité pour H2.
50
Chapitre II Etude expérimentale et modélisation de l’oxydation hydrothermale
Pour démarrer l’installation pilote de manière sûre et stable, il faut suivre les étapes
suivantes:
1. Allumer l’interrupteur principal et l’ordinateur.
2. Charger le programme Vijeo Look sous le logiciel du système de contrôle, et passer
à l'état « Run ». Les variables seront affichées et le contrôle de projet du pilote
ouvert, en sélectionnant l'état Démarrage manuel : allumer la pompe, ouvrir la
vanne manuelle des réservoirs du liquide… etc.
3. La pression commence à monter dans le système jusqu’à atteindre
automatiquement une valeur de référence, introduite dans le PID correspondant, par
exemple 50 bars. Une fois stabilisée, la pression augmente à nouveau en
51
Chapitre II Etude expérimentale et modélisation de l’oxydation hydrothermale
introduisant une nouvelle valeur de référence dans le PID, par exemple 75 bars, et
ainsi de suite jusqu'à atteindre environ 250 bars.
4. Au même temps, le compresseur d'air peut être activé en appuyant sur le bouton
d'alimentation. La vanne de régulation du débit d'air est alors ouverte en cliquant
sur l'icône correspondante sur l'écran de l'ordinateur et sur l'écran PID
correspondant.
5. À ce stade, le contacteur des résistances est activé sur l'écran PID correspondant, on
place un signal de référence pour la température et il est modifié en mode
automatique, de sorte qu'à la jonction des différents courants à l'entrée du réacteur
atteint la température de l’entrée demandée.
6. Au moment où toutes les conditions optimales pour l'essai en question sont
atteintes, la conduite du déchet à traiter est ouverte, puis la conduite d'eau est
fermée et le test démarre.
7. Attendre la stabilité des valeurs des variables mesurées, pendant tout ce temps, les
données de pression, les débits et les températures sont collectés dans le serveur de
données activé au début du test.
8. Une fois l'état d'équilibre atteint, prélever plusieurs échantillons de liquide à la
sortie du système.
9. Une fois qu’un nombre suffisant d’échantillons a été prélevé, il est possible de
modifier certaines variables d’opération afin de réaliser un nouveau test, pour
lequel il sera nécessaire d’attendre à nouveau que l’état d’équilibre soit atteint.
10. Une fois la série de tests terminée, le flux de déchets est remplacé par de l'eau pour
commencer le nettoyage du système. La vanne correspondante du réservoir d'eau
est ouverte et le réservoir du déchet est fermé. Le préchauffage est éteint et la
température dans le système diminue progressivement. Le système reste sous
pression jusqu'à ce que la température baisse suffisamment et que l'eau nettoie le
système.
11. Diminuez ensuite progressivement la pression à des intervalles de 50 bars, par la
suite, la conduite d'air est désactivée.
12. Attendre jusqu’à ce que de l'eau propre sorte du système, sans aucun résidu,
ensuite, éteindre la pompe et arrêter l'eau de refroidissement.
13. Fermer le programme Vijeo Look et déconnecter l’ordinateur.
52
Chapitre II Etude expérimentale et modélisation de l’oxydation hydrothermale
53
Chapitre II Etude expérimentale et modélisation de l’oxydation hydrothermale
CFD, telles que les équations de Navier-Stokes et la théorie des modèles de turbulence. Les
relations utilisées pour le calcul des propriétés thermo physiques du milieu réactionnel et le
coefficient de transfert de chaleur sont présentés, ainsi qu’un aperçu sur les méthodes
numériques utilisées pour la résolution des modèles obtenus et la méthodologie pour la
récupération d’énergie.
∂(ρ. u⃗)
+ ∇. (ρ. u⃗. u⃗) = −∇P + ∇. (τ) + ρ. g⃗ (2. 2)
∂t
Avec : P : la pression (Pa), 𝑔⃗ : le vecteur gravité (m/s2).
𝜏̿ exprime le tenseur des contraintes visqueuses et il est donné par la relation suivante :
2
𝜏̿ = μ. (∇u⃗ + ∇u⃗ ) − ∇. u⃗ 𝐼 (2. 3)
3
54
Chapitre II Etude expérimentale et modélisation de l’oxydation hydrothermale
∂ ̇ . u⃗)
(ρ. E) + ∇(ρ. E. u⃗) = −∇. q − ∇(P. u⃗) + Q̇ − ∇(τ (2. 4)
∂t
Avec : E : désigne l’énergie totale par unité de masse, q : vecteur du flux de chaleur, Q̇ :
source de chaleur volumique.
Conservation des espèces
L’équation de conservation des espèces est écrite comme :
∂
ρy + ∇. ρ. u⃗. y = −∇. J⃗ + R + S (2. 5)
∂t
Avec : J⃗ : le flux de diffusion des espèces, yj: la fraction massique des espèces j, Rj : est le
taux de production nette des espèces j par réaction chimique, Sj : est le terme source de la
phase dispersée et n’importe quelle source définie par l'utilisateur [110].
Modèles de turbulence applicables au flux SCW
La turbulence est une propriété du champ d'écoulement et se caractérise par une large
gamme d'échelles d'écoulement: les échelles les plus grandes qui dépendent de la
géométrie, les échelles les plus petites, et toutes les échelles intermédiaires [111].
Dans le cas d'un écoulement turbulent, il est nécessaire de résoudre les équations (2. 1), (2.
2), (2. 4) et (2. 5) sur des cellules de petites tailles afin de représenter correctement toutes
les échelles de la turbulence [12].
L'analyse théorique et la prévision de la turbulence a été, et jusqu'à ce jour, le problème
fondamental de la dynamique des fluides ; notamment ; de la Dynamique des Fluides
Computationnel (CFD). Les principales difficultés sont dues à la nature aléatoire ou
chaotique des phénomènes de turbulence. Aussi le problème réside dans le fait qu’il n’y a
pas de séparation entre les grandes et les petites échelles : on ne peut pas traiter le
comportement à grande échelle sans tenir compte explicitement des mécanismes mis en jeu
à petites échelles, et inversement. Les deux attitudes, sont donc fortement liées l’une à
l’autre.
Il existe principalement trois méthodes pour modéliser les écoulements turbulents, la
simulation numérique directe (Direct Numerical Simulations : DNS), la simulation des
55
Chapitre II Etude expérimentale et modélisation de l’oxydation hydrothermale
grandes échelles (Large Eddy Simulations : LES), et les méthodes statistiques (Reynolds
Averaged Navier-Stokes : RANS) [110].
La méthode RANS avec des modèles de turbulence appropriés est largement utilisée pour
les applications industrielles [112], et c’est celle qui est la plus utilisée par les codes CFD
commerciaux, pour cette raison, dans la présente étude nous nous attarderons tout
particulièrement sur l’approche RANS.
Équations gouvernantes: Reynolds Moyenne Navier-Stokes (RANS)
L’approche RANS implique l'application des propriétés d'écoulement moyen pour la
résolution de l'équation de Navier Stokes, et pour la fermeture de toute équation de
transport utilisée pour modéliser la turbulence [110].
Une fois l'écoulement est devenu turbulent, toutes les quantités fluctuent dans le temps et
dans l'espace. Alors que la moyenne de Reynolds est appliquée pour les écoulements
incompressibles, une méthode plus pratique pour l'écoulement compressible est la
moyenne de Favre [27]. Dans ce cas la moyenne pondérée en fonction de la densité est
souvent utilisée pour tenir compte de fortes fluctuations de densité, tout en conservant la
forme relativement simple des équations obtenues par l'approche classique de moyenne de
Reynolds [113]. Bien que la décomposition de Favre et de Reynolds soient utilisées,
l'équation obtenue sera appelée équation de « Reynolds Averaged Navier-Stokes équation:
RANS » [114].
Sur la représentation moyenne de Reynolds, chaque variable est décomposés en deux
termes, un terme moyen et un terme fluctuant lié au petites échelles à une valeur de
Reynolds moyen nulle, ce qui donne alors pour n’importe quelle quantité scalaire du flux ɸ
[111]:
ф=ф+ф (2. 6)
Où ф représente la partie moyenne et ф représente la partie fluctuante.
Par l’application de la décomposition de Reynolds, des termes supplémentaires
apparaissent dites tenseurs de Reynolds, et le système devient ouvert. Pour la fermeture de
ce système, un modèle de turbulence sera utilisé, ce modèle va donner une description plus
réaliste des phénomènes d’interaction turbulente en suivant l’évolution de chaque
paramètre turbulent par des équations de transport.
Différents modèles de turbulence invoquent différentes hypothèses sur les échelles non
résolues résultant en différents degrés de précision, pour différents cas d'écoulement [111],
56
Chapitre II Etude expérimentale et modélisation de l’oxydation hydrothermale
tel que modèles algébriques Zero-Equation, les modèles à une-Equation, les modèles à
deux équations, les modèles de fermeture de second ordre.
Les deux modèles turbulents connus pour donner des résultats raisonnablement précis sont
le k-ε et le k-ω de modèles à deux équations.
Dans ce travail, le modèle utilisé est le k-ε de Launder et Spalding, qui résout deux
équations de transport supplémentaires pour l'énergie cinétique turbulente (k) et le taux de
dissipation de la turbulence (ε) [113]. L'énergie cinétique turbulente est la mesure de
l'intensité de la turbulence dans un écoulement de fluide, donc c'est la mesure de l'énergie
cinétique contenue dans la région turbulente d'un écoulement de fluide [80]. D'autre part, le
taux de dissipation de turbulence est défini comme le taux auquel l'énergie cinétique
turbulente est dissipée sous forme d'énergie interne dans le système, telle que l'énergie
thermique. Ces deux paramètres sont à la base de la modélisation de la turbulence k-ε, ils
seront présentés selon le logiciel COMSOL Multiphysics version 5.2 dans la partie (2.3.3).
57
Chapitre II Etude expérimentale et modélisation de l’oxydation hydrothermale
fortement présente dans les réacteurs supercritiques tubulaires, dont le Reynolds peut
atteindre les 10 000. Compte tenu de la complexité des modèles de turbulence et du cas
bidimensionnel, un logiciel recommandé par la littérature pour ce type d’étude a été
utilisé : le COMSOL Multiphysics [60], [27].
COMSOL est un logiciel qui a été largement utilisé pour la modélisation et la simulation
physique, c’est un environnement interactif puissant, permettant de modéliser et de
résoudre tous types de problèmes scientifiques, et d’ingénierie basés sur des équations aux
dérivées partielles (EDP). L’avantage de ce logiciel, c’est qu’on peut coupler et combiner
plusieurs phénomènes physiques lors d’une même simulation numérique, ce qui permet
d’étendre des modèles conventionnels pour un type de physique à des modèles
Multiphysiques résolvant des phénomènes couplés simultanément : écoulement de fluide,
transfert de chaleur, réaction chimique, électrochimie, mécanique des structure,…etc.
Le module COMSOL CFD permet de modéliser les écoulements turbulents à l'aide de
modèles physiques déjà mis en œuvre. Les interfaces de la physique des écoulements des
fluides de ce module sont basées, comme dans tout code numérique pour la dynamique des
flux et le transfert de chaleur, sur les lois de conservation de quantité de mouvement, de la
masse, de l'énergie et de la concentration des différentes espèces dans le milieu réactionnel.
Ces lois de la physique sont traduites en équations aux dérivées partielles et sont résolues
avec les conditions initiales et aux limites spécifiées [111].
Le processus de modélisation et de simulation numérique sous le logiciel COMSOL
Multiphysics 5.2 comporte les étapes suivantes:
Sélection du mode d'application approprié dans le navigateur de modèle, et la prise
en compte des différents phénomènes physiques pouvant exister dans le problème
considéré.
La définition globale des paramètres et des variables liés à notre modèle.
La définition de la géométrie.
Le choix des matériaux et des propriétés du fluide utilisé.
Le maillage.
La résolution du problème à l'aide de l'un des différents solveurs.
Visualisation des résultats en mode post-traitement.
a) Les physiques choisies
La première étape pour réaliser un modèle sur le logiciel COMSOL est de choisir les
physiques appropriées avec les phénomènes physiques et chimiques relatives à notre
58
Chapitre II Etude expérimentale et modélisation de l’oxydation hydrothermale
b) Définition de la géométrie
Il est possible sur COMSOL Multiphysics de créer des modèles 1D, 2D, 2D
axisymétriques et 3D, dans lesquels la géométrie, la physique et le maillage peuvent être
définis, le tout sans aucun logiciel supplémentaire.
Du fait de l’augmentation du coût de calcul entre 2D et 3D, à la fois en mémoire et en
temps de traitement, nous avons considéré une modélisation bidimensionnelle (2D) pour la
résolution de notre problème de la simulation du réacteur tubulaire du pilote de l’UCA
présenté dans la partie (2.2.1), avec une géométrie axisymétrique, les directions examinées
sont la longueur (L) égale à 9.52 m, et le rayon (r) dont le diamètre intérieur D_in, égal à
12,32 mm, comme montré dans la Figure (2. 9) suivante.
59
Chapitre II Etude expérimentale et modélisation de l’oxydation hydrothermale
60
Chapitre II Etude expérimentale et modélisation de l’oxydation hydrothermale
Chaque région de la figure ci-dessus est représentée par une équation. À partir de cette
formulation, nous avons développé un programme sous le logiciel Mathcad; ce programme
détermine la variation de toutes les propriétés thermodynamiques de l'eau en fonction de la
température et de la pression dans l'intervalle qui couvre les cinq régions (densité,
viscosité, conductivité thermique et capacité thermique). Nous avons comparé nos résultats
avec ceux des valeurs tests de W. Wagner et H.J. Kretzschmar [26], et nous avons trouvé
un très bon accord. La formulation IAPWS utilisée ainsi que le programme développé et sa
validation sont bien détaillés dans notre thèse de Magistère [1].
A partir des résultats du programme développé, les données sont transférées dans le
programme principal COMSOL et une interpolation linéaire est réalisée pour chaque
propriété en fonction de la température avec une pression fixe (25 MPa). La température
est variée entre 298,15 et 1000 K et la fonction d'interpolation est définie par un fichier
contenant les valeurs de la fonction en points discrets, comme montré sur la Figure (2.11)
suivante.
61
Chapitre II Etude expérimentale et modélisation de l’oxydation hydrothermale
Dans ce travail, nous avons considéré les deux cas de figures : le premier, lorsque l’eau
pure est considérée comme milieu réactionnel pour lequel les propriétés IAPWS
interpolées sont utilisées, et le deuxième en considérant le milieu réactionnel comme
mélange de toutes les espèces présentes dans le milieu réactionnel dont plus de 70% sont
de l’eau. Pour les propriétés de l’azote, l’oxygène, le CO2 elles proviennent à une pression
constante de 25 MPa de la base de données NIST [115] (le polluant étant considéré
négligeable).
Pour le calcul de ces propriétés, le mélange des différentes espèces est ensuite supposé
idéal, c’est-à-dire que chaque propriété du mélange est calculée par la somme pondérée à
partir d’une moyenne massique de chaque composé selon l’expression suivante [86] ,
[34], [12] :
62
Chapitre II Etude expérimentale et modélisation de l’oxydation hydrothermale
∇. (ρ. 𝐮) = 0 (2. 8)
2 2
ρ. (𝐮. ∇)𝐮 = ∇ −PI + (μ + μ )(∇𝐮 + (∇𝐮) ) − (μ + μ )(∇𝐮)I − ρkI (2. 9)
3 3
μ
ρ. (𝐮. ∇)k = ∇ μ + ∇k + P − ρε (2. 10)
σ
μ ε ε
ρ. (𝐮. ∇)ε = ∇ μ+ ∇ε + C P −C ρ (2. 11)
σ k k
k
μ = ρC (2. 12)
ε
2 2
P =μ ∇𝐮: (∇𝐮 + (∇𝐮) ) − (∇𝐮) − ρk. ∇𝐮 (2. 13)
3 3
63
Chapitre II Etude expérimentale et modélisation de l’oxydation hydrothermale
Avec: r la vitesse de réaction (mol/(m3.s)) et ∆H (T) est la chaleur de la réaction (kJ / mol).
64
Chapitre II Etude expérimentale et modélisation de l’oxydation hydrothermale
Avec:
Di représente le coefficient de la diffusion (m2/s) de l’espèce i, et ri représente la vitesse de
réaction (mol/ (m3.s)), ci la concentration en espèces i (i = polluant, O2, H2O et CO2
(mol/m3)).
Les trois caractéristiques physiques de ce modèle: écoulement turbulent, transfert de
chaleur dans les fluides, chimie et transport d'espèces diluées sont couplées pour obtenir le
modèle global. En effet, les équations de transport donnent le champ de vitesse; le transfert
de chaleur représenté par l'équation (2. 14) utilise ce champ de vitesse et donne la variation
de température, et la réaction représentée par l'équation (2. 17) donne la distribution des
espèces en utilisant le champ de vitesse et la température.
65
Chapitre II Etude expérimentale et modélisation de l’oxydation hydrothermale
66
Chapitre II Etude expérimentale et modélisation de l’oxydation hydrothermale
D D D
1 1 D . ln D D . ln D D . ln D D (2. 20)
= + + + +
H hi λ1 λ2 λ3 D . he
Avec :
H : représente le coefficient de transfert de chaleur global (W/m2.K).
hi, he : représentent respectivement les coefficients de transfert de chaleur interne et
externe (W/m2.K).
Din, D2, D3 et Dext : représentent respectivement le diamètre intérieur du réacteur, Fiberfrax
Durablanket S, la coque en aluminium et le diamètre extérieur (m).
67
Chapitre II Etude expérimentale et modélisation de l’oxydation hydrothermale
Le
Débit de
P Débit de coefficient
Comsopsés T (°C) Concentration l’air Référence
(MPa) l’eau (kg/h) d’excès en
(kg/h)
O2
A.Fourcaul
25 430 5 g/L 12.5 3.96 1
t et al [18]
B. G
Phénol
25 374.3 12.36 g/L 10.18 3.49 3 .Jarana
[34]
25 428 19 g O2/L 10.6 4.07 4.2 V.Vadillo
Biocut 35
25 388 94.4 g O2/L 10.1 4.22 1 et al [43]
Déchet J.M.
d’huile 25 405 44.2 g O2/L 14.4 4.54 1.5 Benjumera
d’olive [51]
68
Chapitre II Etude expérimentale et modélisation de l’oxydation hydrothermale
La formule proposée par J. S-Oneto et al [41] pour l’huile de coupe « Biocut 35 » , à partir
de cette analyse est C6H17O.
- Pour la formule chimique de déchet d’huile d’olive, nous avons utilisé l’analyse
élémentaire réalisée par Hafidi et al [120], qui donne :
69
Chapitre II Etude expérimentale et modélisation de l’oxydation hydrothermale
Concernant l'oxydation de Biocut "35", les paramètres cinétiques utilisés dans ce travail
sont ceux développés par J. Shànchez -Oneto et al [59], où ils avaient constaté que
l'évolution de la DCO suivait un schéma en deux étapes: la première étape est une
élimination rapide où la plupart des substances organiques sont oxydées (95%) et la
deuxième étape est une oxydation lente du reste des composés.
Au niveau de notre programme bidimensionnel principal sous COMSOL, une expression
conditionnelle a été utilisée, en fonction du taux d’élimination du Biocut 35.
f) Le Maillage
La création du maillage est une étape cruciale dans le contexte de la méthode des éléments
finis, que COMSOL utilise pour résoudre l'ensemble des équations différentielles.
Sur COMSOL, on peut soit utiliser la valeur par défaut, via l'utilisation d'une option
appelée "maillage contrôlé par la physique", soit créer notre propre maillage pour l'adapter
70
Chapitre II Etude expérimentale et modélisation de l’oxydation hydrothermale
Figure 2. 12: Le maillage de calcul avec des éléments triangulaires sur le domaine
intérieur et des éléments quadrilatéraux sur les parois.
71
Chapitre II Etude expérimentale et modélisation de l’oxydation hydrothermale
2.3.4.
.3.4. Optimisation de certains paramètres du modèle utilisé
72
Chapitre II Etude expérimentale et modélisation de l’oxydation hydrothermale
et de préciser des contraintes. Les données d'entrée du modèle, qu'il s'agisse des
dimensions géométriques, des formes d'une pièce, des propriétés d’un matériau ou des
distributions spatiales des matériaux, peuvent être traitées comme variables de contrôle.
Tout type de données de sortie du modèle peut être utilisé pour définir la fonction
objective. Dans notre programme, nous l’avions utilisé pour optimiser la valeur du
coefficient de transfert de chaleur.
Dans ce cas, nous avons choisi de créer un objectif des moindres carrés (Figure (2. 14)),
avec la méthode de Levenberg-Marquardt avec une tolérance de 10-6, chaque
caractéristique de l’objectif du moindre carré correspond à une expérience dans laquelle les
mesures ont été obtenues. Par rapport au profil de la température; nous cherchons par
exemple la valeur du coefficient de transfert de chaleur qui correspond le mieux aux
données expérimentales du profil de la température. Cette démarche comprend les étapes
suivantes :
D'abord, nous définissons la fonction objective : cette fonction permet de décrire le
comportement du système simulé en rapport à un objectif précis, la fonction
objective est la somme des carrés de déviations entre le profil de température
obtenu à partir de notre simulation et les valeurs expérimentales de température.
Ensuite, nous définissons un ensemble de variables de contrôle, c'est-à-dire les
données d'entrée du modèle qu’on veut modifier, en lien avec l'objectif.
Troisièmement, nous précisons un ensemble de contraintes : les limites des
intervalles de valeurs des variables de contrôle ou des conditions de
fonctionnement qui doivent être remplies.
Enfin, nous utilisons le module Optimisation pour optimiser notre design en
ajustant ces variables.
73
Chapitre II Etude expérimentale et modélisation de l’oxydation hydrothermale
74
Chapitre II Etude expérimentale et modélisation de l’oxydation hydrothermale
L’influence de ces trois paramètres a été étudiée sur les deux variables dépendantes
suivantes:
L’efficacité de traitement ou le taux d’élimination,
La chaleur à récupérer dans le réacteur.
Les réponses du DOE (Design Of Experiments) considéré ont été obtenues en remplaçant
les expériences par des exécutions numériques du programme 1D développé.
Deuxièmement, et une fois notre modèle bidimensionnel sous COMSOL Multiphysics a
été validé, une étude paramétrique similaire a été réalisée en utilisant un plan d’expérience
factoriel complet, en étudiant l’influence dans ce cas de trois variables indépendantes:
La température à l’entrée du réacteur,
La concentration de la charge exprimée en DCO,
Le débit de l’effluent à l’entrée du réacteur.
L’influence de ces trois paramètres a été étudiée sur les trois variables dépendantes:
L’efficacité d'élimination de la DCO,
La température maximale dans le réacteur,
La chaleur récupérée dans le réacteur.
Les réponses du DOE considéré ont été obtenues en remplaçant toujours les expériences
par des exécutions numériques du programme bidimensionnel développé sous COMSOL
pour les déchets réels : huile de coupe Biocut 35 et le déchet d’huile d’olive.
Le traitement numérique du DOE a été effectué dans le cas du phénol, huile de coupe
Biocut 35 et le déchet d’huile d’olive à l'aide du logiciel Minitab16.
2.5. Conclusion
75
Chapitre II Etude expérimentale et modélisation de l’oxydation hydrothermale
Pour le calcul des propriétés, on avait étudié les deux cas de figures : considérer les
propriétés du milieu comme étant celles de l’eau pure et ensuite le considérer comme un
mélange. L’optimisation de quelques paramètres a aussi été introduite dans le but
d’améliorer le modèle 2D étudié.
76
Chapitre III Résultats et discussion
CHAPITRE III
RESULTATS ET DISCUSSION
3.1. Introduction
Dans ce chapitre seront présentés les résultats obtenus ainsi que leur discussion.
Dans l’objectif d’une maitrise du procédé de l’oxydation hydrothermale pour une
éventuelle application industrielle en Algérie, nous avons essayé d’étudier tous les aspects
de la SCWO, dans la limite de la disponibilité des données et des équipements, allant du
transfert de chaleur, modélisation et simulation du procédé, récupération énergétique et
stabilité du réacteur. Nous commencerons par la présentation des propriétés de l’eau dans
les conditions supercritiques sous le logiciel COMSOL, et la présentation des propriétés en
cas du mélange eau_air, ensuite seront présentés les résultats du programme basé sur le
modèle bidimensionnel, ainsi que la comparaison avec le modèle unidimensionnel
développé sous le logiciel Mathcad durant le travail de magister. Les résultats de l’étude du
coefficient de transfert de chaleur seront ensuite présentés, avec une nouvelle corrélation
ajustée basée sur la corrélation de Dittus-Boelter.
Une fois le programme validé par des valeurs expérimentales de traitement de molécule
modèle et des déchets réels, il sera exploité pour étudier l’effet des conditions de
fonctionnement telles que la température, la concentration d'alimentation en polluants et le
débit d'entrée du réacteur sur les profils de la température et de la concentration. A la base
de ces profils, une analyse utilisant un plan d’expérience factoriel complet (DOE) sera
réalisée afin d’étudier l’effet des variables indépendantes, telles que la température
d’entrée, la concentration initiale en polluant et le débit sur des variable dépendantes telles
que le rendement, la température maximale et la chaleur récupérée. Ces résultats seront
exploités par la suite dans l'étude de sécurité du fonctionnement du réacteur.
77
Chapitre III Résultats et discussion
3.2. Modélisation à 2D
Pour la simulation, il est nécessaire que les propriétés thermo physiques soient
décrites le plus précisément possible, et comme l’état supercritique n’est pas inclut
directement dans les variations des propriétés de l’eau dans la base des données du
COMSOL, nous avons commencé par inclure les données à pression de fonctionnement et
procédé à une interpolation linéaire à partir de nos résultats du programme Mathcad, ces
propriétés sont donc inclues dans le programme en utilisant l’option « interpolation » dans
la rubrique « définition ». Les propriétés de l’eau sur le logiciel COMSOL sont présentées
dans la Figure (3.1), pour des températures entre 300 K et 1000 K et une pression de 22.1
MPa, 25 MPa et 30 MPa.
Notons que ces résultats peuvent être utilisés dans n’importe quel travail de
modélisation dans les conditions supercritiques et non supercritiques de l’eau dans ce
logiciel, dans l’intervalle de température entre 273.15 K et 2273.15 K et une pression
inférieure à 100 MPa, ils peuvent donc être utilisés comme un code de calcul pour générer
les propriétés de l’eau dans cet intervalle de température et de pression sur le logiciel
COMSOL.
78
Chapitre III Résultats et discussion
(a) (b)
(c) (d)
79
Chapitre III Résultats et discussion
températures expérimentales, elles ont aussi été inclues par interpolation à la base de
données du COMSOL.
Les résultats montrent la variation de toutes les propriétés thermo-physiques utilisées dans
le programme bidimensionnel COMSOL, à savoir: la masse volumique, la viscosité
dynamique, la conductivité thermique et la capacité calorifique en fonction de la
température et pour différents pourcentages de l’eau.
(a) (b)
(c) (d)
80
Chapitre III Résultats et discussion
pour la densité car il affecte la solubilité des composés organiques et des gaz dans le milieu
réactionnel [12]. Aussi, la différence entre les propriétés de l’eau pure et celles du mélange
est fortement influencée par le pourcentage surtout à l’état sous-critique et ce sont la masse
volumique et la conductivité thermique qui en sont le plus influencées. Après le point
critique la différence est négligeable. Le pic de la capacité calorifique diminue en fonction
de la fraction de l’air, et la masse volumique augmente avec la fraction en eau et diminue
avec la température pour des températures inférieures à 718 K, pour des températures
supérieures à 718 K, les valeurs de la masse volumique diffèrent peu quelle que soit la
composition eau/air, ceci peut s’expliquer par l’évolution de la masse volumique de l’eau
pure relativement à celle de l’air : bien que l’eau soit nettement plus dense que l’air à
température ambiante, cette tendance s’inverse au-delà d’une température d’environ 718 K.
81
Chapitre III Résultats et discussion
D’après la Figure (3. 3) et dans le cas sans réaction chimique, le modèle bidimensionnel
développé sous COMSOL donne un meilleur résultat et représente mieux les résultats
expérimentaux avec un coefficient de corrélation R2 égal à 0.96, contre un R2 de 0.7 pour
le modèle 1D. On remarque que la déviation entre le modèle 1D et les valeurs
expérimentales à la sortie du réacteur (trouvée aussi par le modèle 1D de Foucault et al
[18]) n’apparait plus avec le modèle 2D.
En présence de la réaction chimique, les Figures (3. 4) et (3. 5) suivantes présentent une
comparaison des profils de la température et de la concentration du phénol expérimentaux,
modèle 1D et modèle 2D.
82
Chapitre III Résultats et discussion
D’après la Figure (3.4), on peut voir que le profil de température simulé en 1D est en
accord avec le profil de température expérimental. Cependant, il y a une déviation pour des
points du réacteur surtout à l’entrée et à la sortie, par contre le profil du modèle 2D est très
proche du profil expérimental notamment à l’entrée et à la sortie du réacteur avec un
coefficient de corrélation de R2=0.952 contre un R2= 0.79 pour le modèle 1D.
Concernant le profil de la concentration du phénol présenté sur la Figure (3.5), pour lequel
nous n’avons que la dernière valeur expérimentale à la sortie du réacteur et non le profil
complet (la conception du pilote ne permet pas une mesure de la concentration à différents
points), le modèle 1D nous donne un résultat meilleur que le modèle 2D, toutefois la
différence n’est pas aussi significative que pour le profil de température (déviation de 0.5
% pour le modèle 1D et 4.63% pour le modèle 2D) et n’est basée que sur un seul point.
Le modèle 2D apporte une précision supplémentaire sur le profil de la température dans le
réacteur par rapport au modèle 1D. On note que c’est le premier travail de comparaison
entre modélisation 1D-2D dans les réacteurs supercritiques. Ce résultat est très important et
permet d’atteindre le premier objectif de cette thèse et qui est d’améliorer le résultat obtenu
avec le modèle 1D, pour lequel, les déviations à l’entrée et à la sortie du réacteur étaient
attribuées en partie au modèle. Ceci valide aussi le programme 2D et permet de l’exploiter
pour l’étude de la SCWO.
Après avoir testé le modèle bidimensionnel avec l’hypothèse de l’eau pure comme milieu
réactionnel, on présente les profils dans le cas ou le milieu réactionnel est supposé un
83
Chapitre III Résultats et discussion
mélange : eau, N2, O2 et CO2. Le programme calcule toutes les propriétés du mélange à
partir d’une moyenne massique de chaque composé selon l’expression (2. 7).
La Figure (3. 6) présente la comparaison entre le profil de la température dans les deux cas
avec les résultats expérimentaux. Les résultats montrent une légère amélioration surtout à
l’entrée du réacteur entre le profil de la température dans le cas ou les propriétés sont
égales à celles de l’eau pure ou le coefficient de corrélation R2 est égal à 0.952, et celles du
mélange réactionnel avec un R2 égal à 0.966.
Par rapport à la concentration, la Figure (3. 7) présente une comparaison pour le profil de
la concentration et le résultat montre aussi une légère amélioration : le taux d’élimination
du phénol dans le cas du mélange est plus proche du résultat expérimental avec une
déviation de 2.7%, au lieu de 4.63% dans le cas des propriétés de l’eau pure.
Figure 3. 6 : Comparaison entre les profils de la température dans le cas des propriétés
égales à celles de l’eau pure et de mélange.
84
Chapitre III Résultats et discussion
Figure 3. 7 : Comparaison entre les profils de la concentration dans le cas des propriétés
égales à celles de l’eau pure et du mélange.
Le fait que l’amélioration ne soit pas très significative, s’explique par le fait que le milieu
réactionnel dans ce cas soit constitué de plus de 70% de l’eau.
Donc d’après ces résultats, dans la section suivante, on va considérer que les propriétés du
milieu réactionnel sont celles du mélange, et on va prendre en considération les propriétés
des espèces présentes dans le milieu réactionnel.
Par rapport au temps de calcul, l’exécution du programme 1D sous le logiciel Mathcad est
très lente, le premier point est obtenu après environ 3 mn et la courbe de variation de
température est obtenue après 130 mn, par contre le programme 2D développé sous le
logiciel COMSOL, l’exécution est très rapide, les résultats sont obtenus pour tous les
points après 168 s, sachant que le processeur du micro ordinateur utilisé est de 2.4 GHz.
Le modèle 2D sera donc utilisé dans cette partie pour l’étude de l’oxydation dans
l’eau supercritique du composé modèle, qui est le phénol et les deux déchets réels qui sont
l’huile de coupe Biocut 35 et le déchet de l’huile d’olive.
a) Le phénol
Deux simulations ont été réalisées en utilisant le modèle bidimensionnel à partir des
conditions de fonctionnement des deux expériences déjà présentés au Tableau (2.4), pour
l'étude expérimentale de A. Fourcault et al [18] et B.G. Jarana et al [34], obtenus à partir
85
Chapitre III Résultats et discussion
du réacteur de pilote de l’université de Cádiz. Dans ce cas nous avons utilisé directement la
valeur expérimentale du coefficient de transfert de chaleur qui est de 30 W/m2·K [18], [34]
et nous avons considéré que les propriétés du milieu réactionnel sont celles du mélange
réactionnel dont 77% est de l’eau dans la première expérience et 74 % pour la deuxième.
Sur les Figures (3. 8) et (3. 9), nous présentons une comparaison entre les profils de
température et de concentration expérimentaux et les résultats de simulation de notre
modèle 2D avec les deux conditions opératoires différentes.
R2=0.966 R2=0.935
86
Chapitre III Résultats et discussion
Comme on peut voir sur les Figures (3. 8) et (3. 9), le programme 2D donne un bon accord
pour les profils de température pour les deux conditions d’entrée (T=430 °C, C=5 g/L et
T=374.3 °C, C=12.36 g/L) avec un coefficient de corrélation R2 égal à 0.966 et 0.935
respectivement.
Dans la première expérience, en raison de la faible concentration du phénol (5 g/L) et les
pertes thermiques importantes (H=30 W/m2 K), la température diminue au dessous de
370°C à la sortie du réacteur, avec une efficacité d'élimination de 79.2% (une déviation de
2.7 % avec les résultats expérimentaux). Dans la deuxième expérience, malgré le fait que
les pertes thermiques sont importantes, et la température d'entrée est relativement faible
(374.3°C), la température diminue légèrement en raison de la concentration d'entrée élevée
du phénol (12.36 g/L), avec une efficacité d'élimination de phénol de 68.93% proche de la
valeur expérimentale qui est 62%.
Les champs de la température et de la concentration à l’intérieur du réacteur simulés par le
programme 2D, pour les deux conditions d’entrée sont présentés sur la Figure (3. 10) et (3.
11) respectivement, et les Figures (3. 12) et (3. 13) montrent le changement de profil de la
température et de la concentration dans la direction radiale à différentes positions axiales
de longueur 1 m jusqu’à 9 m, et cela pour la première expérience.
87
Chapitre III Résultats et discussion
88
Chapitre III Résultats et discussion
(a) (b)
(a) (b)
Figure 3. 13 : Profil de la concentration en fonction du rayon du réacteur (a), avec un plan
de coupe à la sortie (b), « expérience une, phénol ».
Comme on peut voir sur les Figures (3. 10) et (3. 12), la température augmente en fonction
de la longueur lorsque la réaction a lieu, en raison de la chaleur dégagée par l’oxydation
exothermique, et avec les pertes thermiques dans la paroi du réacteur, la température dans
les couches de fluide adjacentes à cette dernière diminue. Dans la dernière section du
réacteur, la température à la paroi baisse en dessous de point critique, par conséquent la
89
Chapitre III Résultats et discussion
densité du milieu réactionnel, qui est composé de plus de 70 % d’eau, augmente près de la
paroi, c’est pourquoi la solubilité du polluant est faible dans cette région.
On note qu’aucun travail relatif à la modélisation des réacteurs d’oxydation supercritique
n’avait présenté les profils de la température et de la concentration en fonction du diamètre
ou rayon du réacteur.
90
Chapitre III Résultats et discussion
R2=0.8 R2=0.95
Figure 3. 14 : Comparaison entre les profils de température expérimentaux et simulés de
SCWO de Biocut 35.
Les résultats montrent un bon accord pour le profil de la température principalement pour
les conditions : température d'entrée = 428 °C et DCO d'entrée = 19 g O2/L. Pour une
température plus basse et une DCO plus élevée, le modèle est moins précis, notamment à
l'entrée du réacteur.
91
Chapitre III Résultats et discussion
92
Chapitre III Résultats et discussion
(a) (b)
Figure 3. 18 : Profils de la température (a) et de la concentration (b) en fonction du rayon
du réacteur suivant différentes longueurs « Biocut 35 ».
Les résultats montrent une dégradation rapide du Biocut 35 au début du réacteur et à des
températures élevées à l'axe du réacteur qui diminuent en fonction de la longueur du
réacteur.
La température est minimale à l’entrée, et elle augmente rapidement après environ le
deuxième mètre du réacteur, ou elle peut atteindre une valeur maximale de 511 °C, due à la
93
Chapitre III Résultats et discussion
chaleur dégagée par les réactions d'oxydation fortement exothermique de l’huile de coupe
Biocut 35, ensuite elle diminue en fonction de la longueur.
L’utilisation du modèle 2D nous permet de connaitre exactement la position dans le
réacteur ou la réaction a lieu, et d’éviter les points chauds dans ce dernier, qui diffèrent en
fonction du type de déchets et des conditions opératoires, ainsi la géométrie
bidimensionnelle permet d'affiner la configuration du réacteur avec, par exemple, plusieurs
injections d’oxydants, ou injections de l’eau froide dans le but de contrôler la température
dans le réacteur [122].
10Mpa
1800
15Mpa
1600
22.1MPa
1400
25MPa
1200
30MPa
1000
800
600
400
200
0
550 600 620 630 640 647 650 660 670 680 700 750 800
Température (K)
94
Chapitre III Résultats et discussion
D’après ce résultat, le coefficient de transfert de chaleur calculé n’est pas constant et varie
légèrement le long du réacteur entre 8.57 et 8.89 W/m2·K, ce résultat concorde avec
certains travaux de la littérature, où il a été prouvé que le coefficient de transfert était en
fonction de la température et du débit [34], toutefois la valeur moyenne a une erreur
relative de 0.108 par rapport à la valeur expérimentale.
Dans la Figure (3. 21), nous avons présenté le profil de la température en utilisant les
données d'entrée présentées dans le Tableau (2.3) [43], [93], et ce, en utilisant le modèle
2D avec un coefficient de transfert ayant : la valeur expérimentale de V. Vadillo [93] et
ceux obtenus par notre calcul.
95
Chapitre III Résultats et discussion
Figure 3. 21: Comparaison du profil de la température dans le cas sans réaction chimique
avec CTC calculé, sa valeur expérimentale et le résultat expérimental.
96
Chapitre III Résultats et discussion
(a) (b)
Figure 3. 22: Comparaison entre le profil de température avec CTC calculé, expérimental
et optimisé en fonction de L (a) et en fonction du D (b) Tracé de surface.
Comme on peut voir sur la Figure (3. 22), le profil de la température avec CTC optimisé
est superposé au profil de la température dans le cas où le CTC est calculé, cependant la
déviation par rapport au dernier point persiste, cela peut être dû au fait de travailler avec un
coefficient de transfert constant le long du réacteur, or la variation de la température peut
avoir une influence sur son estimation, d’autant plus que les variations des propriétés
physiques sont importantes.
Dans l’objectif d’améliorer encore plus le résultat, nous avons procédé à l’optimisation des
constantes de la corrélation de Dittus-Boelter présentée dans l’équation (2. 18). En effet,
ces constantes, sont modifiées ou bien ajustées selon les résultats expérimentaux par
plusieurs auteurs [83], certaines de ces modifications ont été présentées au Tableau (1. 2).
Les résultats de l’optimisation en utilisant notre programme Comsol sont:
a=0.3158, b=0.17007, c= 0.35023.
La nouvelle corrélation devient donc:
Nu = 0.3158. Reb0.17007. Prb0.35023 (3. 2)
La Figure suivante (3. 23) présente le profil de la température obtenu en utilisant la
corrélation ajustée, avec une comparaison avec le profil de la température dans le cas de
CTC expérimental, calculé à partir de corrélation de Dittus-Boelter, calculé à partir de
corrélation de Dittus-Boelter ajustée et CTC optimisé.
97
Chapitre III Résultats et discussion
Le module d'optimisation permet d'avoir une corrélation ajustée pour le calcul de CTC qui
a donné un profil de température plus précis avec un coefficient de corrélation de 0.975 et
une diminution de l’écart à la sortie; et sera donc à utiliser dans la modélisation de
l’oxydation supercritique dans ce réacteur. La précision des résultats du modèle dépend
fortement de la valeur du coefficient de transfert de chaleur.
Après avoir obtenu une corrélation ajustée pour le calcul du coefficient du transfert
de chaleur, nous avons appliqué notre programme 2D pour validation avec les résultats
expérimentaux du déchet d’huile d’olive, obtenus par le groupe de recherche de
l’université de Cádiz, sur le réacteur tubulaire du pilote, et pour lesquels l’étude de
simulation et le calcul du CTC n’ont pas été réalisés.
Concernant les propriétés du milieu réactionnel, nous les avons considérées comme celles
du mélange (eau, oxygène, dioxyde de carbone et azote); les contributions de déchet de
l'huile d’olive peuvent être considérées comme négligeables, car sa concentration dans
l'eau est inférieure à 0,03% en poids et le pourcentage de l’eau est de 78%. La chaleur
globale de la réaction a été obtenue à partir de la formule moléculaire du déchet de l’huile
d’olive : C5H8O2 obtenue par analyse élémentaire, et en utilisant l'expression empirique (3.
1) trouvée par M. J. Cocero [121].
98
Chapitre III Résultats et discussion
99
Chapitre III Résultats et discussion
100
Chapitre III Résultats et discussion
(a) (b)
Figure 3. 27 : Profils de la température (a) et de la concentration (b) en fonction du rayon
du réacteur pour différentes longueurs « déchet d’huile d’olive ».
Après validation du programme 2D, une simulation numérique a été réalisée afin de
mieux comprendre le comportement du réacteur. Pour cela, le programme 2D a été utilisé
pour étudier l'effet de certains paramètres tels que la température, la concentration en
polluant et le débit à l'entrée du réacteur sur le profil de la température et le taux
d’élimination du phénol, l’huile de coupe Biocut 35 et le déchet d’huile d’olive.
a) Phénol
Les Figures (3. 28) et (3. 29) montrent l'effet de la température initiale, de la concentration
initiale et du débit d’effluent à l’entrée du réacteur sur le profil de la température (Figures
3. 28 : a, b, c) et le taux d'élimination (Figures 3. 29 : a, b, c). La température à l’entrée
varie entre 375 et 460°C, la concentration initiale de 10 à 30 g/L et le débit de l’effluent de
9 à 25 kg/h.
101
Chapitre III Résultats et discussion
102
Chapitre III Résultats et discussion
Figure 3. 30: Effet de la température initiale (a), de la concentration initiale (b) et le débit
(c) sur le profil de la température de Biocut 35.
103
Chapitre III Résultats et discussion
Figure 3. 31 : Effet
et de la température initiale (a
(a), de la concentration initiale (b)
(b et le débit
(c)) sur le profil du taux d’élimination du Biocut 35.
104
Chapitre III Résultats et discussion
Ces résultats relatifs aux profils de température sont en accord avec les résultats de
simulation unidimensionnel réalisés par V. Vadillo et al [43] pour l’huile de coupe Biocut
35.
c) Le déchet d’huile d’olive
Nous avons aussi étudié l’effet de la température, la DCO et le débit à l’entrée sur les
profils de la température et de la concentration de déchet d’huile d’olive le long du
réacteur, les résultats sont présentés sur les Figures (3. 32) et (3. 33) suivantes.
La température à l’entrée varie toujours entre 375 et 420°C, la concentration initiale de 35
à 71 gO2/L et le débit de l’effluent de 9 à 25 kg/h.
105
Chapitre III Résultats et discussion
106
Chapitre III Résultats et discussion
Le programme élaboré 1D, a été validé à l'aide des données expérimentales réalisé
au niveau du pilot de SCWO de l’université de Cádiz, comme on a déjà motionné, par A.
Faucault et al [18], avec R2 de 0.79 pour la température, et valeur proche pour la
concentration à la sortie avec une déviation de 0.5%. Nous l'avons exploité pour effectuer
d'autres conditions de fonctionnement afin d'étudier l'effet de la température, de la
concentration et du débit d'effluent à l'entrée du réacteur sur la réaction d’oxydation
supercritique du phénol en utilisant un plan d’expérience factoriel complet 23; remplaçant
les expériences sur pilote, (difficiles à réaliser toutes expérimentalement avec des
variations importantes sur chaque facteur), par des exécutions numériques du programme
1D selon la matrice DOE utilisée sous le logiciel Minitab.
Les trois facteurs indépendants considérés sont la température d'entrée (Tin), la
concentration du phénol à l’entrée (Cin) et le débit d'alimentation (Qin), avec deux niveaux
pour chaque facteur correspondant au minimum (-1 en valeur codée) et au maximum (+1
107
Chapitre III Résultats et discussion
Les flux froid et chaud du procédé ppeuvent être schématisés sur la Figure
igure (3. 34) suivante :
La chaleur requise pour refroidir le flux de sortie et qui peut d'être réutilisée pour
préchauffer le flux d'entrée, est calculée selon l’expression suivante :
108
Chapitre III Résultats et discussion
D’après les résultats du Tableau (3. 2), deux expériences permettent d’atteindre un taux
d’élimination presque de 100%, à savoir l’expérience 2 et l’expérience 4, avec une
température à l’entrée maximale et un débit minimal, concernant la chaleur à récupérer,
elle est maximale pour la température et la concentration initiale les plus élevées.
109
Chapitre III Résultats et discussion
Les corrélations mathématiques de premier ordre reliant les réponses étudiées aux facteurs
ont été obtenues comme suit :
110
Chapitre III Résultats et discussion
(a) Effets des facteurs principaux (b) Graphique des effets standardisés de Pareto (c) Effets d'interaction des facteurs
111
Chapitre III Résultats et discussion
le Tableau (3.4), l’effet des facteurs significatifs est beaucoup plus marqué pour la chaleur
à récupérer que pour le taux d’élimination, les facteurs température, concentration initiale
et leur interaction ont les valeurs du T-Student de 320.57, 226.19 et 222.83 respectivement.
L’effet de l’interaction dans ce cas est très important.
(a) Effets des facteurs principaux (b) Graphique des effets standardisés de Pareto (c) Effets d'interaction des facteurs
L’analyse de variance (ANOVA) de ce modèle est présentée dans le tableau (3. 4) suivant :
112
Chapitre III Résultats et discussion
L'importance de l'application d’un plan des expériences dans cette partie du travail est
d'étudier l'effet des facteurs cités ci-dessus, et aussi de trouver un optimum. Comme
présenté ci-dessus, les trois facteurs étudiés n'ont pas eu la même influence sur les deux
réponses: taux d'élimination, et la quantité de la chaleur à récupérer. Pour cela, nous
utilisons l’Optimizer de Minitab pour trouver les conditions optimales afin de maximiser
les deux réponses au même temps. Les résultats de cette étape sont présentés sur les
Figures (3.37.a) et (3.37.b):
(a) (b)
Figure 3. 37 : Optima codés pour maximiser l'efficacité d'élimination (a) et la quantité de
la chaleur à récupérer (b).
Dans cette partie, nous avons exploité le programme 2D élaboré sous le logiciel
COMSOL Multiphysics pour étudier l’effet des conditions opératoires, température,
113
Chapitre III Résultats et discussion
La chaleur requise qui peut être réutilisée pour préchauffer le flux d'entrée à partir de flux
de sortie du réacteur est calculée selon l’expression suivante :
Q_rec_
Q_rec_polluant = m·Cp (T_out-T_ex) (3. 7)
114
Chapitre III Résultats et discussion
L'objectif de cette partie était de quantifier les effets, les interactions et d'obtenir des
corrélations entre les variables dépendantes et indépendantes de l'étude, et cela dans le cas
réel et avec les deux déchets réels, l’huile de coupe Biocut 35 et le déchet d’huile d’olive.
L'analyse des données du plan factoriel a été réalisée à l'aide du logiciel Minitab 16.
A. Huile de coupe Biocut 35
Pour l’huile de coupe Biocut 35, les trois facteurs indépendants considérés, comme nous
avons mentionné, sont la température d'entrée, la concentration d'entrée exprimé en DCO
et le débit d'effluent, avec deux niveaux pour chaque facteur correspondant au minimum (-
1 en valeur codée) et au maximum (+1 en valeur codée).
Les réponses du DOE étaient le taux d’élimination, la température maximale à l’intérieur
du réacteur et la chaleur récupérée.
La limite de température de sécurité du réacteur étant fixée entre 550 ° C et 600 °C [94],
[123]. A partir des résultats présentés sur la Figure (3. 30), la température peut atteindre
des valeurs non autorisées à l’intérieur du réacteur, pour cette raison, nous avons décidé
d’exploiter les résultats du plan d’expérience pour délimiter une «zone d'opérations de
sécurité».
Le Tableau (3. 5) présente les valeurs codées et réels pour chaque facteur. Les valeurs de la
concentration sont estimées à partir des valeurs minimale et maximale déjà réalisées sur ce
pilote.
Tableau 3. 5: Variables codées et réelles utilisées dans le plan factoriel
complet « Biocut 35».
Variable -1 +1
Température initiale (°C) 375 420
Concentration initiale DCO (g O2/L) 19 50
Le débit de l’effluent (kg/h) 9 25
La matrice du plan factoriel complet utilisée et les réponses pour les huit expériences
correspondantes au DOE 23 est présentée sur le Tableau (3. 6). Par conséquent, 8
exécutions numériques du programme développé 2D en fonction des conditions de
fonctionnement de la matrice DOE nous donnent la température à la sortie (T_out), le taux
d’élimination (Y_DCO), la température maximale (T_max) et la chaleur récupérée
(Q_rec_Bio) dans le procédé comme suit :
115
Chapitre III Résultats et discussion
Q_rec_Bio =16067.5+1990 T+997.5 DCO+577.5 Q+530 T DCO+65 T Q+117.5 DCO Q (3. 10)
116
Chapitre III Résultats et discussion
(a) Effets des facteurs principaux (b) Graphique des effets standardisés de Pareto (c) Effets d'interaction des facteurs
117
Chapitre III Résultats et discussion
b) Température maximale
La Figure (3. 40) présente les principaux effets, les graphes d'interactions et le diagramme
de Pareto pour la réponse température maximale dans le réacteur. La température s'est
avérée être l'effet positif le plus important, suivi par la DCO, par contre, le débit d’effluent
n’a presque aucun effet. L'interaction température-DCO est la seule interaction
significative. Comme confirmé par graphique de Pareto Figure (3. 40. b) et l'analyse de la
variance présentée sur le Tableau (3. 8), la température, la concentration et leur interaction
sont significatives avec une valeur de p < 0,05. Cependant, le débit d’effluent et ses
interactions n’ont pas d’effet significatif avec une valeur de p> 0,05.
(a) Effets des facteurs principaux (b) Graphique des effets standardisés de Pareto (c) Effets d'interaction des facteurs
118
Chapitre III Résultats et discussion
c) La Chaleur récupérée
La Figure (3. 41) présente les principaux graphes des effets et des interactions ainsi que le
graphe de Pareto pour la réponse chaleur récupérée. Pour cette réponse, et comme montré
sur le diagramme de Pareto de la Figure (3. 41.b), la température est également l'effet le
plus important avec un effet positif, suivi de la DCO initiale et le débit, et comme confirmé
par l'analyse de la variance présentée sur le Tableau (3. 9), la température est le seul facteur
significatif avec p <0,05, les interactions à deux niveaux ne sont pas significatifs avec une
valeur de p> 0,05.
(a) Effets des facteurs principaux (b) Graphique des effets standardisés de Pareto (c) Effets d'interaction des facteurs
119
Chapitre III Résultats et discussion
SE
Term Effect Coef T p
Coef
Constant 16067.5 135.0 119.02 0.005
T 3980 1990 135.0 14.74 0.043
DCO 1995 997.5 135.0 7.39 0.086
Q 1155 577.5 135.0 4.28 0.146
T*DCO 1060 530 135.0 3.93 0.159
T*Q 130 65 135.0 0.48 0.714
DCO*Q 235 117.5 0.7563 0.87 0.544
Dans ce qui suit, nous utilisons le logiciel Minitab pour trouver les conditions optimales
afin de maximiser les deux réponses étudiées : le taux d’élimination et la chaleur récupérée
dans le réacteur, au même temps. Les résultats de cette étape sont présentés sur la Figure
(3.42) ci-dessous:
120
Chapitre III Résultats et discussion
Tableau 3. 10: Variables codées et réelles utilisées dans le plan factoriel complet « déchet
d’huile d’olive».
Variable -1 +1
Température initiale (°C) 375 420
Concentration initiale DCO (g O2/L) 35 71
Le débit de l’effluent (kg/h) 9 25
La matrice du plan factoriel complet utilisée et les réponses pour les huit expériences
correspondantes au DOE 23 sont présentées sur le Tableau (3. 11). Par conséquent, 8
exécutions numériques du programme développé 2D en fonction des conditions de
fonctionnement de la matrice DOE nous donnent le taux d’élimination, la température
maximale, et la chaleur récupérée dans le procédé.
121
Chapitre III Résultats et discussion
Tableau 3. 11: Matrice du plan factoriel complet et réponses « Déchet d’huile d’olive ».
D’après les résultats du plan factoriel présentés sur le Tableau (3. 11), dans ce cas aussi, les
deux expériences 4 et 8 ont conduit au dépassement important de la température limite à
l’intérieur du réacteur, pour lesquelles la température et la DCO à l’entrée sont maximales,
ce sont aussi les deux réponses qui permettent d’avoir le maximum de la chaleur à
récupérer et le taux d’élimination le plus élevé.
Les corrélations mathématiques du premier ordre reliant les réponses étudiées aux facteurs
ont été obtenues comme montrées sur les équations suivantes :
Y_DCO=95.611+2.844 T+1.231 DCO-2.111 Q-0.686 T DCO +1.906 T Q+ 1.019 DCO Q (3. 11)
T_max=516.954+139.481 T+99.806 DCO-0.819 Q+97.554 T DCO-0.801 T Q-0.766 DCO Q (3. 12)
Q_rec_DHO=18913+4555 T+3535 DCO+3360 Q+3058 T DCO+2787 T Q+ 2892 DCO Q (3. 13)
Les coefficients de détermination correspondants sont de 0,966, 0.999 et 0,891 pour le taux
d’élimination, la température maximale et la chaleur récupérée respectivement.
a) Le taux d’élimination
La Figure (3. 43) présente les principaux graphes des effets et des interactions ainsi
que le graphe de Pareto pour la réponse taux d’élimination dans le cas de déchet d’huile
d’olive. La température initiale a l'effet positif le plus important; suivi du débit d'effluent
qui a un effet négatif et enfin la DCO initiale. Selon le diagramme de Pareto de la Figure
122
Chapitre III Résultats et discussion
(3. 43. b) et l’analyse de variance présentée sur le Tableau (3.12), tous les facteurs et leurs
interactions à deux niveaux ne sont pas significatifs avec une valeur de p> 0,05.
(a) Effets des facteurs principaux (b) Graphique des effets standardisés de Pareto (c) Effets d'interaction des facteurs
Tableau 3. 12: Analyse de la variance pour le taux d’élimination « déchet d’huile d’olive».
SE
Term Effect Coef T p-value
Coef
Constant 95.611 0.8138 117.49 0.005
Température 5.687 2.844 0.8138 3.49 0.177
DCO 2.463 1.231 0.8138 1.51 0.372
Q -4.222 -2.111 0.8138 -2.59 0.234
Température*DCO -1.372 -0.686 0.8138 -0.84 0.554
Température*Q 3.813 1.906 0.8138 2.34 0.257
DCO*Q 2.037 1.019 0.8138 1.25 0.429
123
Chapitre III Résultats et discussion
b) La température maximale
Sur la Figure (3. 44), les principaux effets, les graphes d'interactions et le diagramme de
Pareto pour la réponse température maximale sont présentés dans le cas de déchet d’huile
d’olive. La température présente l'effet positif le plus important, suivi par la concentration,
et le débit n’a aucun effet sur la température maximale. Dans ce cas et suivant le graphique
Pareto et l'analyse de la variance présentée sur le Tableau (3. 13), la température, la
concentration et leur interaction sont significatives avec une valeur de p < 0,05, cependant,
le débit d’effluent et leurs interactions ne sont pas significatifs avec une valeur de p> 0,05
et l'interaction température-DCO est la seule interaction significative. Donc pour une
question de sécurité du réacteur il faut éviter de travailler avec un maximum de
température et de DCO à l’entrée.
(a) Effets des facteurs principaux (b) Graphique des effets standardisés de Pareto (c) Effets d'interaction des facteurs
124
Chapitre III Résultats et discussion
c) La chaleur récupérée
La Figure (3. 45) présente les principaux graphes des effets et des interactions ainsi que le
graphe de Pareto pour la réponse chaleur récupérée dans le cas de déchet d’huile d’olive.
Pour cette réponse, et comme montré sur le diagramme Pareto présenté sur la Figure (3.
45.b), la température a également l'effet positif le plus important, suivi de la concentration
initiale et du débit d’effluent. Comme confirmé par l'analyse de la variance présentée sur le
Tableau (3. 14), tous les facteurs et leurs interactions à deux niveaux ne sont pas
significatifs avec une valeur de p> 0,05.
(a) Effets des facteurs principaux (b) Graphique des effets standardisés de Pareto (c) Effets d'interaction des facteurs
125
Chapitre III Résultats et discussion
Tableau 3. 14: Analyse de la variance pour la chaleur récupérée « déchet d’huile d’olive».
SE
Term Effect Coef T p-value
Coef
Constant 18913 2920 6.48 0.098
Température 9110 4555 2920 1.56 0.363
DCO 7070 3535 2920 1.21 0.440
Q 6720 3360 2920 1.15 0.455
Température*DCO 6115 3058 2920 1.05 0.485
Température*Q 5575 2787 2920 0.95 0.515
DCO*Q 5785 2892 2920 0.99 0.503
Dans ce qui suit nous utilisons l’Optimizer de Minitab pour trouver les conditions
optimales afin de maximiser les deux réponses : taux d’élimination et chaleur récupérée, au
même temps. Les résultats de cette étape sont présentés sur la Figure (3.46) ci-dessous:
126
Chapitre III Résultats et discussion
D’un autre coté, on note que d’un point de vue économique, l'augmentation de la
température initiale implique des coûts d'investissement et d'exploitation plus élevés;
tandis que l’augmentant de la DCO initiale n'a pas, en principe, de coûts supplémentaires
associés, car de nombreuses eaux usées industrielles sont déjà très concentrées.
Par conséquent, par rapport à la température maximale, la température et la concentration
d'alimentation sont des paramètres critiques de fonctionnement, en raison de la réactivité
élevée de l'huile de coupe Biocut 35 et de déchet d’huile d’olive, la température augmente
rapidement au niveau de la première section du réacteur. Ainsi, pour maintenir des
conditions de sécurité pour le pilote, la concentration et la température à l’entrée du
réacteur doivent être contrôlées pour éviter de dépasser la limite de température du
matériau.
127
Chapitre III Résultats et discussion
(a) (b)
Figure 3. 47 : Zone de sécurité pour éviter une valeur maximale de 550 ° C (a) et 500 ° C
(b) dans le réacteur et assurer un taux d'élimination de DCO de 96% et 375 ° C comme
température de sortie « cas de l’huile de coupe Biocut 35 ».
Les résultats présentés dans cette étude, relatifs au déchet réel l’huile de coupe Biocut 35,
ont fait l’objet d’une publication scientifique dans le journal international: « International
Journal of Chemical Reactor Engineering » [124].
En ce qui concerne le déchet d’huile d’olive, les résultats suivants (Figure (3. 48))
présentent la zone de sécurité, et cela toujours pour un débit de 17 kg/h, température
maximale du réacteur limitées à 500 ° C et 550 °C, une efficacité de traitement minimale
de 96%, et une température minimale à la sortie de 375 ° C.
128
Chapitre III Résultats et discussion
(a) (b)
Figure 3. 48 : Zone de sécurité pour éviter une valeur maximale de 550 ° C (a) et 500 ° C
(b) dans le réacteur et assurer un taux d'élimination de DCO de 96% et 375 ° C comme
température de sortie « cas de déchet d’huile d’olive ».
La zone de sécurité se réduit pour des contraintes plus sévères (température maximale de
500°C). Pour une température maximale ne dépassant pas 500 °C, la concentration de
déchet d’huile d’olive ne doit pas dépasser 45 g O2/L, cette concentration augmente à
environ 57 gO2/L si on augmente la température maximale jusqu’à 550 °C. La surface de la
zone de sécurité dans les mêmes conditions change d’un polluant à un autre.
Dans la partie qui suit, on va fixer la température maximale à l’intérieur du réacteur à une
limite de 550 °C, et étudier l’effet du débit sur la zone de sécurité, en maintenant
l’efficacité de traitement minimale à 96 %, et la température minimale à la sortie de 375
°C.
Les résultats sont présentés sur les Figures (3. 49) et (3. 50) pour l’huile de coupe Biocut
35 et le déchet d’huile d’olive successivement, pour différentes valeurs de débit (9, 15, 20
et 25 kg/h).
129
Chapitre III Résultats et discussion
(a) (b)
(c) (d)
Figure 3. 49 : Effet du débit sur la zone de sécurité « cas de l’huile de coupe Biocut 35 »,
(a) : 9 kg/h, (b) : 15 kg/h, (c)=20 kg/h et (d) : 25kg/h.
130
Chapitre III Résultats et discussion
(a) (b)
(c) (d)
Figure 3. 50 : Effet du débit sur la zone de sécurité « cas de déchet d’huile d’olive », (a) :
9, (b) : 15 kg/h, (c)=20 kg/h et (d) : 25kg/h.
Les résultats montrent que la zone de sécurité dans le cas de l’huile de coupe Biocut 35 est
légèrement influencée par le débit, par contre dans le cas du déchet d’huile d’olive cette
zone dépend significativement du débit et on a plus de sécurité pour un débit de 15 kg/h.
131
Chapitre III Résultats et discussion
3.4. Conclusion
Les résultats présentés dans ce chapitre peuvent être résumés comme suit :
-Le modèle 2D a permis d’améliorer les résultats obtenus avec le modèle 1D, et d’autres
améliorations ont été apportées en ajustant le coefficient de transfert de chaleur. En effet ce
paramètre est très important dans le modèle et son estimation avec précision affecte
visiblement le résultat du modèle.
- La comparaison réalisée avec le modèle 2D entre le cas où le milieu réactionnel est
supposé de l’eau pure et le cas mélange des espèces présentes dans ce milieu, le deuxième
cas s’est avéré plus approprié.
-Une fois le programme validé, le modèle dans sa forme adoptée, a été utilisé pour étudier
le comportement du réacteur en réalisant une étude paramétrique classique qui a montré
l’effet des conditions opératoires : température, concentration du polluant et débit de
l’effluent à l’entrée, sur les profils de la température et de la concentration. Cette étude a
montré clairement des zones de température maximale dans le réacteur, qui risque parfois
et de dépasser les limites de sécurité.
- Une étude paramétrique utilisant la méthodologie des DOE a été réalisée en exploitant le
modèle 1D pour le phénol, afin de quantifier l’effet des conditions opératoires étudiées sur
les deux réponses : taux d’élimination et la chaleur à récupérer, les résultats ont montré que
la température avait un effet positif significatif sur le taux d’élimination, et avait avec la
concentration un effet positif significatif sur la chaleur à récupérer.
-Une autre étude paramétrique utilisant la méthodologie des DOE a ensuite été réalisée en
exploitant cette fois le modèle 2D afin de quantifier l’effet des conditions opératoires
étudiées sur les réponses précédentes en incluant en plus, le taux d’élimination, la
température maximale dans le réacteur et la quantité de la chaleur récupérée. Les résultats
ont montré, que la température avait un effet positif sur le taux d’élimination, il était
significatif dans le cas de l’huile de coupe Biocut 35. Pour le cas de déchet d’huile d’olive,
la température et la concentration ainsi que leur interaction avaient un effet significatif sur
la température maximale.
Concernant la quantité de chaleur récupérée ; la température avait un effet positif
significatif dans le cas de l’huile de coupe Biocut 35 mais pas pour le déchet d’huile
d’olive.
132
Chapitre III Résultats et discussion
Enfin, une étude de superposition des corrélations obtenues avec limites imposées a permis
de délimiter des zones de sécurité, pour la première fois, pour un réacteur SCWO et aussi
de montrer la variation de cette zone en fonction du débit et de la variation des contraintes.
133
Conclusion générale et perspectives
134
Conclusion générale et perspectives
135
Références
REFERENCES
136
Références
137
Références
138
Références
139
Références
140
Références
141
Références
142
Références
143
Annexes
Taux
Température Pression Type de
Composants d’élimination Référence
(°C) (MPa) réacteur
(%)
Huile de coupe 422 25 Tubulaire 96.17 [61]
Lexivat
400 25 Batch 100 [48]
d’enfouissement
PCB 539 24.1 Tubulaire 99.6 [60]
Phénol 430 25 Tubulaire 76.5 [18]
Eaux usées
Industrielles 500 25 Tubulaire 90 [125]
Inflammables
Ethanol 405 25 Tubulaire 96.2 [51]
cokéfaction des
eaux usées
650 25 Tubulaire 99 [126]
« coking
wastewater »
processus de
fabrication
d'affichage à
cristaux
615 25 / 99.99 [48]
« liquides
(LCD) liquid
crystal
display »
RDX 500 25 Batch < 70 [85]
HMX 500 25.7 Batch 95.8 [85]
NH4ClO4
85%+PBHT 500 25.5 Batch 94.6 [85]
15%
144
Annexes
145
Annexes
Résultat de
Régio P Valeur test de FluidMat Lib Waterstea
T(K) notre
n (MPa) référence 1 97 2 m pro 3
Programme
1 298.15 0.1 8.90022551.10-4 8.90022551.10-4 8.90022551.10-4 8.87.10-4
2 873.15 20 3.39743835.10-5 3.39743835.10-5 3.39743835.10-5 3.398.10-5
3 673.15 60 7.26093560.10-5 7.2609356.10-5 / 7.257.10-5
146
Annexes
Résultat de
P Valeur test de FluidMat Water
Région T (K) notre
(MPa) référence 1 Lib 97 2 steam pro 3
Programme
1 298.15 5 78.04662743 78.5907250 / /
2 873.15 10 1.12620467 1.12620970 / /
3 673.15 40 10.30623674 10.3126058 / /
1 2 3
valeurs tests du livre [26], FluidMat Lib [127] , Waterstream pro [128]
147
Annexes
148
i n t e r n a t i o n a l j o u r n a l o f h y d r o g e n e n e r g y 4 2 ( 2 0 1 7 ) 1 2 9 2 6 e1 2 9 3 2
ScienceDirect
Article history: The hydrothermal oxidation is an efficient and clean way for wastewater treatment; it is
Received 2 August 2016 also an energy self sufficient process under certain conditions, by the mean of the highly
Received in revised form exothermic reaction of oxidation.
3 November 2016 In this present work, we studied the modeling of hydrothermal oxidation of phenol in
Accepted 7 November 2016 polluted water using supercritical conditions. First we apply a simple mathematical model
Available online 1 January 2017 for tubular reactor, resulting from mass, momentum and energy balances, which is often
difficult to solve in its general form. According to adequate simplifications, the profiles of
Keywords: temperature and concentration of chemical species are computed along the reactor. Finite
Supercritical water oxidation differences numerical method was used to solve the obtained system of equations and a
Wastewater treatment MathCAD program was developed for this purpose. The results are close to experimental
Energy recovery works from the literature.
Process modeling To study the influence of some operating conditions on the efficiency of the treatment
Design of experiments and the amount of heat generated by the reaction and which will be recovered to heat the
inlet flows, we exploited the modeling results to perform a full factorial design of experi-
ment (DOE). The results allow fixing the optimal conditions for the studied factors on the
efficiency and energy recovery of the process treatment.
© 2016 Hydrogen Energy Publications LLC. Published by Elsevier Ltd. All rights reserved.
* Corresponding author.
E-mail address: [email protected] (O. Nawel).
http://dx.doi.org/10.1016/j.ijhydene.2016.11.046
0360-3199/© 2016 Hydrogen Energy Publications LLC. Published by Elsevier Ltd. All rights reserved.
i n t e r n a t i o n a l j o u r n a l o f h y d r o g e n e n e r g y 4 2 ( 2 0 1 7 ) 1 2 9 2 6 e1 2 9 3 2 12927
Species conservation
Continuity
vu$r$yj
vru Mj $V j $r ¼ 0 (11)
¼0 (3) vx
vx
V j stands for the stoichiometric coefficient of species j
r stands for the density of mixture in (kg/m3) and u for its
(phenol, O2 and CO2), Mj for the molar weight of species j.
velocity in (m/s).
Operating conditions and parameter values used
Momentum conservation
The simulation was performed with an initial temperature of
vuz r vP vtxy waste to about 450 C and 25 MPa, a flowrate of waste is 12.5 l/h,
¼ þ (4)
vx vx vx
the concentration of phenol is 5 g/l. Air also enters the reactor
where P is the pressure in (Mpa), txy is the stress tensor in in excess at a pressure of 25 MPa, a temperature of 200 C [7].
direction x, linked to a gradient of velocity in the direction y. The kinetic parameters used in this model are shown in the
following Table 1:
vux vuy 2 vuk
txy ¼ meff $ þ $ $dxy (5)
vy vx 3 vxk
Numerical approach
For turbulent flow, the effective viscosity meff is not equal to
dynamic viscosity m as for laminar flow meff ¼ m þ mT where mT is The governing equations of mass, momentum, energy and
the turbulent viscosity [10]. species conservation (3, 4, 10 and 11), which form a system of
In turbulent mode of operation, the stress tensor should differential equations, are solved numerically. Finite differ-
include the Reynolds's stress tensor and an appropriate model ence method with upwind scheme is used to get a purely
to compute it (such: keε and keu models) [11], this complicates algebraic system [13]. For example, the continuity differential
the mathematical description of the process and often need equation becomes in the algebraic form:
CFD software programs. Thus, for mathematical simplifica-
vu$r ui ri ui1 $ri1
tion, and in order to get a one-dimensional problem, the ¼ (12)
vx Dx
divergence of the stress tensor is assumed to be as the local
pressure drop over the control volume under consideration [7]. Dx is the length step, ui the velocity at point i and ui1, the
In this study, this computation has been done according to velocity at the previous mesh point.
Blasius and McAdams correlations. The resulting balance The final system contained the non-linear equations ob-
equation is then: tained after discretization of the differential system, plus the
equations of thermodynamic properties “IAPWS” of pure water:
vu2 r vP density r (T, P), enthalpy h (T, P) and viscosity m (T, P). This
¼ Ck $r$u2 (6)
vx vx algebraic model is solved using an appropriate numerical
The coefficient Ck in the momentum conservation Eq. (6) is method for system of linear and non-linear equations in a
given by: Mathcad code which was implemented with IAPWS properties
f calculations.
Ck ¼ (7)
2$Dh
Heat recovery
Where Dh is the hydraulic diameter of the channel in (m),
and f is the friction factor, which can be determined using the
As mentioned above, the oxidation reaction is highly
Blasius and McAdams correlations for a smooth tube [12],
exothermic and the heat produced can be recovered to rise the
given by:
effluent temperature at the inlet until a desired operating
value. The cold and hot streams in the process can be sche-
0:316:Re0:25 Re < 30000
f¼ (8)
0:184:Re0:20 30000 < Re < 106 matized as follows (see Fig. 2):
r$u$Dh
Re ¼ (9)
m
Table 1 e Kinetics parameters in Eq. (1) for phenol
m is the dynamic viscosity in (Pa s). oxidation [7].
Kinetics parameters
Energy conservation
K ((mol l)1mn s1) Ea (kJ/mol) m n
Using the same assumption as above the energy conservation
1.34 ± 0.77
equation is as follows: 10 39.2 ± 10.07 1 0
i n t e r n a t i o n a l j o u r n a l o f h y d r o g e n e n e r g y 4 2 ( 2 0 1 7 ) 1 2 9 2 6 e1 2 9 3 2 12929
QR is the required heat load for cooling the outlet reactor Some results of the developed program for some needed
stream in (J), m is the effluent mass in (kg), Cp is the specific properties for the modeling step are presented in Table 2,
heat in (J/kg K),Tr and Tc the temperatures at the outlet of the which shows a good agreement with the recommended test
reactor and the outlet of the cooler respectively in K. values for regions 1, 2, 3 and 5. Region 4, represents a satu-
The required heat is estimated in order to be reused as a ration pressure, and the program results are compared in
heating agent of the inlet stream; it depends on inlet Table 3 to the recommended test values [8].
temperature, initial concentration of pollutant and of the Next, the results of modeling the supercritical oxidation
effluent flowrate. Therefore, the influence of these factors will tubular reactor were obtained using the operating data pre-
be studied on the heat production and also on treatment ef- sented above, to validate the model and the developed pro-
ficiency (expressed by the elimination rate Y) by using a gram. These results are compared with the experimental
factorial design of experiments in Minitab software. results of Fourcault et al. [7]. Temperature profile and species
concentration variation along the reactor are presented in
Y ¼ Cini Cf Cini (14) Figs. 3 and 4 respectively.
Cini and Cf are the initial and final phenol concentration in The results agree well with experimental ones, except for
(g/l). the temperature profile at the outlet of the reactor, which may
The outlet temperature of the reactor and final phenol be due to the use of constant heat transfer coefficient to es-
concentration are estimated, for the given operating condi- timate the heat lost with the reactor surroundings. Indeed,
tions, using the above model and the Mathcad program some works related the overall heat transfer coefficient to
developed for simulating the process. temperature and flowrate [14]. For phenol concentration pro-
file (Fig. 4), we have only the experimental value at the outlet
of the reactor which agrees with the modeling results.
To study the influence of inlet temperature, inlet phenol
Results and discussion
concentration and effluent flowrate on treatment efficiency
and also on the heat generated by the reaction, we perform a
First, we present the results obtained from the Mathcad pro-
full factorial design of experiments using these three factors.
gram for physical properties of water and steam using IAPWS
The responses of the considered DOE are obtained by replac-
formulation. To validate the developed program to be used by
ing experiments by numerical executions of the developed
the modeling program, the results were compared with the
Mathcad program. Using operations conditions fixed by the
test values recommended by Wagner [8].
Table 3 e Results of saturation pressure of region 4 using Table 4 e Factors intervals and design of experiments
IAPWS developed program. matrix.
T ( K) Calculated Ps (MPa) Ps (MPa) Test value [8] Factors Low level High level
3 3 1 þ1
300 3.53658941 10 3.53658941 10
500 2.63889776 2.63889776 A: Initial temperature ( K) 673 873
600 12.34431458 12.3443146 B:Initial Concentration (g/l) 5 15
C: Effluent Flowrate (l/h) 10 30
references
[1] Jarana BG, Oneto JS, Portela Migue lez JR, Nebot Sanz E,
Martenez de la Ossa EJ. Simulation of supercritical water
Fig. 11 e Coded optimums for elimination efficiency oxidation with air at pilot plant scale. Int J Chem React Eng
maximization 2010;8. Article A58.
[2] Bambang V, Jae-Duck K. Supercritical water oxidation for the
destruction of toxic organic wastewaters. J Environ Sci
2007;19:513e22.
[3] Aymonier C. Traitement hydrothermal des de chets
industriels spe ciaux : donnees pour le dimensionnement
d’installations industrielles et concepts innovants de
acteurs sonochimique et e
re lectrochimique. Ecole Doctorale
des Sciences Chimiques, Universite de bordeaux 1; The se
Doctorat N 2212. 2000.
[4] Jin Hui, Guo Simao, Guo Liejin, Cao Changqing. A
mathematical model and numerical investigation for
glycerol gasification in supercritical water with a tubular
reactor. J. Supercrit Fluids JULY 2015. http://dx.doi.org/
10.1016/j.supflu.2015.06.028.
[5] Vielcazals S, Mateos D, Mercadier J, Marias F, Bottreau M,
Cansell F. Modelling and simulation of
hydrothermal oxidation of organic compounds. AIChE J
2006;52:1e8.
[6] Mancini F, Cansell F, Marias F, Mercadier J. Hydrothermal
Fig. 12 e Coded optimums for heat generation oxidation of waste from biomass to generate energy. Oral
maximization Commun WasteEng 2005. Albi 17e19 Mai 2005.
[7] Fourcault A, Garcia-Jarana B, Sanchez-Oneto J, Marias F,
Portela JR. Supercritical water oxidation of phenol with air.
From this optimization study, we have to operate at a Exp results Model Chem Eng J 2009;152:227e33.
maximum of initial temperature, maximum of initial con- [8] Wagner W, Kretzschmar H-J. International Steam Tables,
centration and a minimum flowrate, to achieve a 100% pre- properties of water and steam based on the industrial
dicted rate of phenol elimination and a 2.571 105 J of predicted Formulation IAPWS. 2nd ed. Germany: Springer; 1998.
heat to be recovered by the process. [9] Minkowycz WJ, Sparrow EM, Murthy JY. Handbook of
numerical heat transfer. 2nd ed. John Wiley, Sons, Inc;
2006.
[10] Moussie re S, Joussot-Dubien C, Turc H-A, Boutin O,
Guichardon P, Charbit G, et al. Stirred double shell reactor for
Conclusion waste treatment by supercritical water oxidation. In:
Moussie re S, Joussot-Dubien C, Guichardon P, Boutin O,
A mathematical model for the oxidation of phenol in super- Turc H-A, Charbit G, Fournel B, editors. Proceeding of 20th
critical water based on the conservation of mass, energy, AIRAPT- 43th EHPRG, June 27- July 1, Karlsruhe/Germany; 2005.
momentum and chemical species, varying physical properties [11] Abhijeet J, Sharawage J. CFD analysis of supercritical water
reactor for flow back and produced water treatment. Master
according to IAWPS formulation was performed. The obtained
of Science (MS) Mechanical Engineering -Engineering and
temperature and species profiles were compared with exper- Technology of Ohio University. 2014.
imental and numerical results from the literature, the results [12] Prashant Jain K, Rizwan-uddin. Numerical analysis of
validated the developed programs. supercritical flow instabilities in a natural circulation loop.
A factorial design of experiments with Minitab software USA. Nucl Eng Des 2008;238:1947e57.
was also used to study the effect of initial temperature, initial [13] Holmes MH. Introduction to numerical methods in
differential equations. In: Marsden JE, Sirovich L, Antman SS,
pollutant concentration of the effluent and its flowrate on the
editors. Texts in applied mathematics, 52. Springer; 2007.
degree of pollutant elimination and on the amount of the heat
[14] Dhanuskodi R, Arunagiri A, Anantharaman N. Analysis of
to be recovered. The DOE applications allowed us to estimate variation in properties and its impact on heat transfer in sub
the operation conditions to maximize the treatment efficiency and supercritical conditions of water/steam. Int J Chem Eng
and the quantity of heat to be recovered. Appl October 2011;2(No. 5).
Int. J. Chem. React. Eng. 2020; ▪▪▪(▪▪▪): 20200047
attention of researchers more than other materials; some oils in order to study the effect of a highly chlorinated waste
others have suggested inhibiting corrosion by adding basic on reactor corrosion and metal concentration in liquid
compounds to neutralize the corrosive acid species (Tang output effluent. Their experiments were performed at 490 °
et al. 2016). Recently, Voisin, Erriguible, and Aymonier C with a residence time of 15 s and a 88.8% COD reduction
(2020) proposed a new solvent system composed of a and weak metal concentrations in liquid output effluent
molten salt in pressurized water, called hydrothermal had been achieved.
molten salt (HyMoS) that changed the paradigm of the Portela et al. (2001) studied hydrothermal oxidation of
solubility of inorganic compounds in supercritical water. cutting oil waste in normal and promoted hydrothermal
Indeed, their study showed the possibility to precipitate a oxidation in a continuous flow system, when hydrogen
low melting temperature salt, at a temperature above its peroxide had been used as oxygen and free radicals source
melting one, leading to the instantaneous formation of the when introduced into the reactor with or without previous
HyMoS. The molten salt is then capable of dissolving a thermal decomposition, respectively. They observed that
large amount of inorganic salt. organic material was easily oxidized, obtaining more than
Moreover, as the oxidation reactions are exothermic, 90% Total chemical oxygen (TOC) reduction in less than
autothermal operation and even heat production can be 10 s at 500 °C. At lower temperatures, the use of promoters
achieved in the SCWO process (Abelleira et al. 2013) to clearly enhanced the oxidation process.
reduce its operational costs that are strongly influenced by The determination of kinetic parameters for SCWO of
the energy recovery in the process, the throughput, and the cutting oil under oxygen excess conditions had been pro-
heating value of the waste (Vadillo et al. 2011). posed by Sánchez-Oneto et al. (2007), at a temperature
Several research studies have been carried out to test range of 400 to 500 °C and a pressure of 25 MPa, specifically
the effectiveness of this technology for the degradation of for two commercial cutting oils (Servol and Biocut). In the
many organic compounds; these tests have been carried case of Servol, a pseudo-first kinetic model was used.
out on laboratory and/or pilot scale including different However, a two-parameter mathematical model involving
geometries and materials for the reactors, operating pro- two steps (a fast reaction followed by a slow one) was
cedures, different sets of conditions and mathematical needed to describe the Biocut SCWO kinetics and to
treatment of the data, etc. calculate the kinetic constants. For the same range of
Due to the high cost of building and operating of these temperature and pressure, another kinetic model was also
reactors, the development of simulations is of great interest proposed by Sánchez-Oneto et al. (2008) for Biocut when
(Bermejo et al. 2007). The development of mathematical no oxygen excess conditions are used obtaining a 0.58
models is necessary to better understand the phenomena order for oxygen concentration.
that take place within the, to optimize the operating con- Using a one dimensional model and basing on kinetic
ditions and the design and to assess the process scaling up parameters proposed by Sánchez-Oneto et al. (2007),
(Loppinet-Serani, Aymonier, and Cansell 2010). Vadillo et al. (2001) presented a simulation based on the
The present work considers the modeling of SCWO of energy balance of highly concentrated cutting oil in a
cutting oil “Biocut 35” which is mixed with water in an continuous flow tubular reactor. They also carried out 15
appropriate proportion obtaining the corresponding cut- experiments in order to validate the model at different
ting fluid used to provide cooling and lubrication to metal values of initial temperature, feed concentration of cutting
removal processes, to prevent corrosion and facilitate the oil and excess of oxygen (air), and the comparisons showed
ejection of cut metal from machinery processes. The a very good agreement. In the other hand, Benjumea et al.
chemistry, delivery, and disposal of exhausted cutting (2017) developed a mathematical model based on mo-
fluids may have a negative impact on the environment mentum, mass, species, and energy balances, of cutting oil
(Skerlos 2013). The treatment of this pollutant waste must in tubular reactor, in steady and stationary states. When
meet the strict regulatory framework to avoid any risk of the system was considered using a two-dimensional
pollution of the environment and ensure the safety of op- model, very often, the resolution required a Computa-
erators. It is also noted that other technologies were used tional Fluid Dynamics (CFD) software to obtain tempera-
for this pollutant, like electro coagulation and Fenton’s ture and concentrations profiles along the studied reactor.
process for removal of chemical oxygen demand (COD) The response of Supercritical water reactor (SCWR) based
(Patel and Pandya 2016). in two-dimensional had been scarcely studied in literature.
The study of treatment of cutting oil by SCWO had been Marulanda (2011) used the k-ε turbulence model, with a
carried out by several authors like Cansell, Beslin, and two-dimensional representation for a tubular reactor with a
Berdeu (1998), who performed an experiment with cutting considerable length (66 m) to the oxidation of mineral
N. Benmakhlouf et al.: Modeling and simulation of hydrothermal oxidation 3
transformer oil. The work studied only the effect of the 2 Theory
temperature at the inlet on the increase of the temperature
in the reactor and the abatement of total organic carbon. 2.1 Mathematical modeling
The used model was based on the assumption that the
properties of the reaction medium are those of pure water. For the modeling step, the data from a tubular reactor
However the results showed that it was better to use a lower developed in the University of Cádiz in Spain, were used.
inlet temperature in order to avoid hazardous operating The pilot plant had a 25 kg/h capacity and used air as
temperature conditions (inlet temperature lower than oxidant, it had been satisfactorily used in several previous
650 K). These simulations validated the kinetic models studies on SCWO for the degradation of many pollutants
proposed by Sánchez-Oneto et al. (2007), (2008) for cutting such as phenol (Fourcault et al. 2009), and cutting oil
fluids. (Benjumea et al. 2016; Jimenez-Espadafor et al. 2011;
For parametric studies, the design of experiments Vadillo et al. 2011). The tubular reactor of the pilot plant
(DOE) had shown to be a powerful tool to quantify the was made of stainless steel 316 L tubing; with a special wall
effects of operating factors and their interactions, allowing thickness of 3.365 mm to increase its pressure resistance,
a better understanding of the process. This technique had the internal diameter was about 12.32 mm and the total
been widely used in the chemical industry (Sahu and volume of the reactor was 1229.553 cm3. The reactor was not
Andhare 2018), and also for supercritical water oxidation totally adiabatic, so a priori it was necessary to estimate the
studies. value of the heat transfer coefficient between the reactor
Miller et al. (2015) showed that the feed concentration and the external environment.
of fecal sludge had a significant impact on the increase of A two dimensional axisymmetric computational fluid
the temperature, whereas the reactor pressure and oxygen dynamics model of SCWO reactor was used, with initial
excess had the most significant impact on heat transfer operating conditions of inlet temperature about 400 °C,
efficiency. In the other hand, Cardona, Buitrago, and pressure of 25 MPa and effluent mass flow rate of above
Acosta (2017) found that the temperature-time interaction 10 kg/h. In these operating conditions, the turbulence in-
was significant for total organic carbon and total nitrogen side the reactor is expected to be fully developed
destruction efficiencies of landfill leachate, when the (Re > 10,000), so the simulation including the modeling of
optimal temperatures were 400 °C for total organic carbon turbulence in steady-state was used.
and 500 °C for total nitrogen, and optimal time was 30 min. The used COMSOL Multiphysics 5.2 CFD (computa-
In this study, a two dimensional model of hydro- tional fluid dynamics) module provided the ability to
thermal oxidation in supercritical water was developed model turbulent flows with already implemented physical
for a tubular reactor that used air as oxidant. First, the models. The Fluid Flow physics interface in this module are
heat transfer coefficient was calculated, from the tem- based on the conservation of momentum, mass, energy
perature profile in the absence of reaction. Then, the and concentration of different species in the reaction me-
computer code was validated using the experimental re- dium, expressed in terms of partial differential equation
sults of Vadillo et al. (2011) for cutting oil with and without (Module Users Guide 2015). In this work, the choice of the
reaction. Finally, the reactor behavior with respect to the turbulence model is based on the RANS (Reynolds Aver-
operating conditions (feed concentrations, reactor’s inlet aged Navier–Stokes) model because it is one of the most
temperature, and waste water flow rates) was studied widely used in the field of CFD and it is less Central pro-
applying numerical simulations and a full factorial cessing unit (CPU) intensive than LES (Large Eddy simu-
experimental design to quantify main factors and in- lation) or DES (Detached Eddy Simulation) (Breuer, Jovičić,
teractions effects on the outlet temperature, COD removal and Mazaev 2003). Also the standard k-ε model because it
and maximal temperature. The analysis of variance is the most commonly used to simulate supercritical water
(ANOVA) was used to test the obtained first order models (Marulanda 2011; Shrawage 2014). However some works
and the significance effects. Finally, the developed such that of Moussiere et al. (2007) compared the use of two
simulation program and DOE were coupled to generate models for modeling supercritical water reactor. The k-ε
the safety zone, which was obtained by imposing con- had shown better results to simulate the temperature pro-
strains on the studied responses: maximum temperature file of non-reactive flow while the k-ω showed an
in the reactor to ensure reactor safety, COD removal yield improvement of the results for the reactive one.
to ensure treatment efficiency and outlet temperature for For its closure this model introduces two additional
energy recovery possibilities. variables in the Navier–Stokes equation compared to the
4 N. Benmakhlouf et al.: Modeling and simulation of hydrothermal oxidation
laminar flow model, namely the turbulent kinetic energy k sources other than viscous dissipation (W/m3) and Qvd the
and the dissipation of turbulent energy ε. One major reason viscous dissipation (W/m3) was not set by user but calcu-
that makes the k-ε model computationally inexpensive is lated directly by Comsol.
that it uses wall functions to describe the flow close to walls The convective heat flux on the boundaries in contact
instead of resolving the very steep gradients in these regions with the fluid is then modeled as being proportional to the
(Module Users Guide 2015). temperature difference across a fictitious thermal bound-
For two-dimensional channel flow, the steady-state ary layer. Mathematically, the heat flux is expressed by the
mass, momentum, energy conservation and species in the following equation:
reaction medium, can be written as follows:
−n・(−K ∇ T) H(T ext − T) (8)
2.1.1 Mass equation where H is a heat transfer coefficient and Text the temper-
ature of the external fluid.
In addition, the heat source term Q is given by the heat
∇.ρ.u 0 (1)
of the reaction ΔH (T ) and reaction rate r as:
Q −ΔH(T)・r (9)
2.1.2 Momentum equation
where r is the reaction rate (mol/(m3・s)) and ΔH (T) the heat
of reaction (kJ/mol).
ρ.(u.∇)u ∇ − PI + μ + μT ∇ u + (∇u)T The heat of reaction ΔH (kJ/mol) was calculated from
the relationship proposed by Cocero et al. (2002) at tem-
2 2 perature equal to 400 °C and pressure of 25 MPa for model
− μ + μT (∇u)I − ρkI (2)
3 3 compounds:
μT −ΔH 415・ n + 107・m − 193・f (10)
ρ.(u.∇)k ∇μ + ∇ k + P k − ρε (3)
σk
where n represents the number of carbon atoms, m: num-
μT ε ε 2 ber of hydrogen atoms and f: number of oxygen atoms.
ρ.(u.∇)ε ∇μ + ∇ ε + Cε1 P k − C ε2 ρ (4) Since the heat capacity and enthalpy of formation of
σε k k
cutting oil are not available, this is the only analytic
k2 method to calculate the heat of reaction.
μT ρC μ (5)
ε Therefore, this relationship has been used by many
authors and has proven useful as for phenol by Garcia
2 2
P k μT ∇ u : ∇ u + (∇u)T − (∇u)2 − ρk. ∇ u (6) Jarana et al. (2010), for cutting oil by Vadillo et al. (2011).
3 3
2.2 Physical properties the only considered reaction was the overall oxidation
one for cutting oil and is expressed as:
Because of the high water content in the reaction medium 4 C6 H17 O + 39 O2 → 24 CO2 + 34 H2 O (13)
(more than 70%), some studies assumed that the thermo-
dynamic and transport properties of the reacting medium Cutting oil is a real waste water that is a complex
are equal to that of pure water, and they calculated them mixture of several compounds, and in this case, the reac-
from the International Association for the Properties of tion mechanism and the identification of all the involved
Water and Steam (IAPWS) formulation (Fourcault et al. intermediate compounds is a difficult task. Usually, for
2009; Vielcazals et al. 2006). In the present study, it has complex mixtures, the kinetic study, that is necessary for
been considered that the properties in the reaction medium SCWO reactor simulation, is performed based on COD
are equal to the properties of the mixture (water, oxygen, measurements.
carbon dioxide, and nitrogen) and calculated as an average Therefore, developing generalized kinetic models is
with respect to the weight percentage of each pure com- very useful for SCWO studies. Li, Chen, and Gloyna (1991)
pound (Benjumea et al. 2017; García Jarana et al. 2010; Oh, suggested a simplified mechanism in which acetic acid was
Kochan, and Beller 1997; Vadillo et al. 2011): the most typical intermediate compound. Thereafter, Por-
tela, Nebot, and de la Ossa (2001) proposed a new gener-
∑mi ⋅ Propi (T, P) alized kinetic model considering the carbon monoxide as a
Prop(T, P) (12)
∑mi refractory intermediate compound for cases in which the
with Prop (T, P) represents the property calculated from formation of acetic acid is not pertinent. This model was
each substance “i”, evaluated at temperature T and pres- satisfactory used in case of cutting oil (Vadillo et al. 2014).
sure P (25 MPa), mi represents the mass of substance i. First, a fast elimination of most of the organics (95%)
The contributions of cutting oil were considered as takes place, followed by a second slow step to oxidize the
negligible with a concentration in water of less than 0.032 wt remaining compounds (more refractory compounds as
% and moreover its properties in supercritical conditions short chain organic acids, mainly acetic acid and gases as
were not readily available. The supercritical water proper- CO).
ties were obtained from the IAPWS formulation (Wagner – The kinetics for Biocut 35 used in this work were those
and Kretzschmar 2008) by means of a developed computing developed by Sánchez-Oneto et al. (2007) (See Table 1),
code running under Mathcad software (Nadjiba, Nawel, and the evolution of COD with time was found to follow a
Abdesslam Hassen 2017). The obtained water physical two steps kinetic model. First, a fast elimination of most
properties were expressed as functions of temperature using of the organics (95%) followed by a second slow step to
Comsol linear interpolation. These functions were input as oxidize the remaining compounds. The kinetics for the
data into a Comsol code devised for a temperature range of COD concentration for Biocut 35 is expressed as follows:
298.15 to 1000 K and a pressure of 25 MPa for each property. −dCOD −E a
r k ° exp · COD (14)
The other fluid properties were obtained from the NIST dτ RT
Database (NIST|National Institute of Standards and Tech-
with k° is the pre-exponential factor of Arrhenius law (1/s),
nology, n.d.) at 25 MPa and experimental temperatures.
Ea: the activation energy of reaction (J/mol) and R: the
universal gas constant J/(mol・K).
Table : Operating parameters in no reaction case. Table : Range of operating parameters in real and coded values.
and axial geometry is considered. The generation of a mesh boundary flow was handled by using wall functions, and
is a crucial step in the context of finite element method that no-slip condition was set at the walls.
Comsol used to resolve the set of the differential equations
model. Since the considered geometry was simple, Comsol
used a free triangular mesh (2214 elements) within the flow 2.6 Statistical design
stream away from the walls; while boundary layers were
added close to the walls consisting of quadrilateral ele- In a supercritical water oxidation process, it was very
ments (5290 elements) tightly packed in the direction important to estimate the temperature profile in the
normal to the walls. reactor, its value at the outlet of the reactor gave an idea on
As the computational model is time independent, a the amount of heat that can be recovered to preheat the
stationary solver is considered to solve, the non-linear effluent at the inlet, and the maximal temperature in the
equations, using the Damped Newton’s method. For each reactor should not be higher than the temperature limit of
iteration, the non-linear system of equations is linearized the reactor material, to work in safe conditions.
and solved using a linear solver, by default, and for most A parametric study of the process was performed using
simulations, Comsol uses the direct linear solver, the tur- full factorial design of experiments, studying the influence
bulence quantities were solved separately by using a of three independent variables: temperature, feed con-
segregated solver. The convergence criterion of 10−3 for all centration and flow rate at the inlet of reactor, on three
variables is taken for all the iterations. dependent process variables: the outlet temperature, COD
removal efficiency and the maximal temperature in the
reactor. Once the Comsol program was validated, the re-
2.5 Boundary conditions sponses of the considered DOE were obtained by replacing
experiments by numerical executions of the developed
At first, it was necessary to confirm the value of heat program. The levels of real and coded values for the
transfer coefficient between the experimental apparatus operating parameters are given in Table 4. The DOE nu-
and the surroundings, so the model was first used in the merical treatment was performed using Minitab16
case of no reaction, only water and air were introduced in software.
the reactor without the pollutant (a blank experiment).
Table 2 presents the operating parameters used for this
simulation case. 3 Results and discussion
Second, the model was used with chemical reaction,
and Table 3 presents the used operating parameters for 3.1 Properties of the reaction medium
simulating the process in this case.
Concerning the turbulent kinetic energy k and the As mentioned earlier, IAPWS Mathcad code was used to
dissipation of turbulent energy ε, at the inlet they were find the properties of water, with interpolation in terms of
calculated directly in Comsol. For the turbulence model, temperature on Comsol Multiphysics for a fixed pressure of
Compound Inlet pressure Inlet tempera- Initial concentration Waste water Flow rate of Oxidant excess Ref.
(MPa) ture (°C) (gO/L) Flow(kg/h) air (kg/h) coefficient
25 MPa. The results show all the thermal physical proper- losses between the experimental apparatus and the sur-
ties used in code namely: density, viscosity, thermal con- roundings and its heat transfer coefficient play an impor-
ductivity, and heat capacity versus temperature for a tant role to solve the energy balance presented in equations
mixture. The variation of physical properties assuming a (7). In the present work, since the reactor shows important
mixture water-air composed of 70 wt% water is shown heat losses when work at high temperature, the heat
Figure 1. The results show the great variation of the prop- transfer coefficient was often estimated experimentally.
erties versus temperature in the water-air mixture. This For that estimation, temperature profile along the reactor
point is particularly important for density because it affects in a blank experiment was used, when pure water was
the solubility of organic compounds and gases in the re- introduced into the reactor with an oxidant stream in the
action medium (Moussiere 2006). The accuracy of the absence of chemical reaction.
simulation results depends strongly on the accuracy of the Thus, from the work of Vadillo (2012), the heat transfer
estimation of all properties. coefficient for the pilot plant facility at University of Cádiz
was experimentally estimated to be equal to 9.83 W/(m2・K).
Using the input data of blank experiment of the present
3.2 Optimization of heat transfer coefficient work shown in Table 2 and the experimental heat transfer
estimation coefficient cited above, the execution of the Comsol code
(without reaction module) gave the results shown in
The estimation of thermal losses is necessary to develop a Figure 2.
model that can predict accurately the experimental tem- The results show a good agreement with experimental
perature profile along the reactor and the COD removal values particularly in the first section of the reactor, but
efficiency when the reactor behavior is not adiabatic. Heat some deviation can be seen in the last meter of the reactor,
since the obtained profile was much more accurate and agreement mostly for the second conditions (inlet tem-
very close to the experimental case. Therefore, its use was perature = 428 °C and initial COD = 19 g O2/L), for lower
adopted in the modeling of the SCWO with chemical re- temperature and higher COD, the model was less accurate,
action for which the results were very dependent on the especially at the inlet of the reactor.
reliability of the heat transfer coefficient. As shown in Figure 4, the temperature profiles were
modeled with a reasonable accuracy with R2 equal to 0.8 and
0.95 for Experiments 1 and 2, respectively. Due to the heat
3.3 Temperature and concentrations released by the exothermic oxidation reactions of Biocut 35,
profiles of SCWO the reactor temperature increased rapidly in the reactor inlet
especially in case one, when it reached 430 °C because of the
The supercritical oxidation of cutting fluid “Biocut 35” was high concentration of pollutant. This value was acceptable
simulated using the developed Comsol code with the for a safe pilot plant operation which should not exceed the
optimized heat transfer coefficient and operating condi- maximum value of 550 °C at the outlet of reactor and 600 °C
tions of the work of Vadillo et al. (2011). A comparison inside the reactor (Benjumea et al. 2016).
between simulated and experimental temperature profiles In the other hand, the supercritical water oxidation ef-
are presented in Figure 4. The results show a good ficiency for the degradation of Biocut 35 is presented in
N. Benmakhlouf et al.: Modeling and simulation of hydrothermal oxidation 9
Figure 5 where the plot shows the species concentration 3.4 SCWO numerical results
profiles along the reactor but the comparison with experi-
mental can be performed only for the Biocut 35 concentra- After the validation of the computing code, a numerical
tion at the reactor outlet. The predicted COD removal was of simulation was carried out in order to better understand
the order of 96.66 and 95.47% for operating conditions of the behavior of the reactor. For that, the code was used to
experiments one and two, respectively, compared to the study the effect of certain parameters such as temperature,
respective experimental values of 85.2 and 80.8%. COD concentration and flow rate at the inlet of the reactor
Figure 6 shows a 2D representation of the temperature on the profile of both temperature and COD removal along
and Biocut 35 concentration profiles. The results show a the reactor. Figures 7 and 8 show the effect of initial tem-
fast degradation at the entrance of the reactor and at high perature in the range 375–420 °C, initial waste water con-
temperatures at the reactor axis but which decreased ac- centration in terms of COD concentration from 19 to 50 g O2/
cording to the reactor length. As can be seen most of the L and feed flow rate of waste water from 9 to 25 kg/h, on
reaction takes place in the first two meters of the reactor both temperature profiles (a, b, c) and COD removal (d, e, f)
and therefore, in the last seven meters, the temperatures for cutting oil Biocut 35.
decrease because there is no more energy generation by As can be seen in Figure 7(a) and (b), the temperature
reaction and therefore there are important heat losses. and the concentration of pollutant at the inlet had an
Figure 7: Effect of initial temperature (a), initial waste water concentration (b) and feed flow rate of waste water (c) on the temperature profile
for SCWO of Biocut 35.
important effect on the temperature profile of the reactor needed to complete the overall oxidation reactions was
and, the reaction evolution. The increase of initial tem- achieved in a lower reactor volume, as can be observed in
perature and initial concentration increased the maximum Figure 8(f); moreover the outlet temperature decreased,
reactor achieved temperature and therefore the mean since globally, less organic matter was oxidized in the unit
temperature of the reactor, because the reaction rate was of time and therefore a lower quantity of heat was released
very sensitive to temperature changes and the reaction was in the oxidation reactions.
highly exothermic. As a direct consequence, a higher COD It is important to point out that from a practical point
elimination was achieved in a lower residence time that of view, increasing the initial temperature involves
was in shorter reactor length as can be seen in Figure 8(d) higher investment and operation costs; while increasing
and (e). the initial COD does not have, in principle, higher asso-
The increase of the inlet temperature led to a faster ciated costs, since many industrial wastewaters are
initiation of the reaction, according to Figure 7(a), the already highly concentrated. Therefore, is very inter-
oxidation with an inlet temperature of 420 °C was slightly esting to increase the initial COD as far as possible, since
faster than for 410 and 400 °C, so the maximum tempera- it would provide a higher temperature profile inside the
ture achieved was higher as the inlet temperature reactor, and the efficiency of the process (rate of organic
increased, obtaining a higher COD removal efficiency material treated/time) will improve. However, to deter-
closer to that at the reactor inlet, with a lower residence mine the best initial temperature an economic optimiza-
time, leading to a lower reactor volume, hence a lower tion would be necessary. Regarding the feed flow rate to
investment cost in the corresponding industrial applica- be treated and as expected, a calculated reactor volume
tion with the same COD removal. (hence a reactor length) to obtain an adequate COD re-
Regarding to feed flow rate of waste water, as can be movals is necessarily depending on the feed flow rate. It
seen in Figure 7(c), its influence on the temperature profile is convenient to select the lowest ratio Volume/flow rate
in the studied range was lower than that of initial tem- needed for a desired COD removal, but overdesigning the
perature and initial concentration; as soon as the flow rate reactor is not convenient because it would decrease the
decreased, the maximum temperature could be reached to temperature at the reactor outlet even if a small increase
be closer to the reactor inlet because the residence time in COD removal is achieved.
12 N. Benmakhlouf et al.: Modeling and simulation of hydrothermal oxidation
Figure 8: Effect of initial temperature (d), initial waste water concentration (e) and feed flow rate of waste water (f) on the rate of elimination for
SCWO of Biocut 35.
3.5 The full factorial design According to that, the levels of each factor were fixed to
avoid that limit of temperature.
Since the elaborated code was validated using the avail- Table 5 presents the used full factorial design matrix
able experimental data, it was exploited to test different and the responses for the eight experiments corresponding
values of operating conditions in order to study the effect of to a 23 DOE. Therefore, eight numerical executions of the
temperature, concentration and effluent flow rate at the developed code under the operating conditions of the DOE
inlet of reactor on the SCWO reaction using factorial matrix leading to the outlet temperature, percentage of
experimental. In all the cases, the pressure was 25 MPa and COD removal and maximal temperature achieved along the
the air flow rate was 4 kg/h. The objective of this part was to reactor.
quantify the effects, the interactions and to get a mathe-
matical relationship between dependent and independent 3.5.1 Outlet temperature, percentage of COD removal
variables of the study. and maximal temperature analysis
The three considered independent factors were inlet
temperature, inlet COD concentration and feed flow rate The data analysis of the factorial design was performed
with two levels for each factor corresponding to the mini- using Minitab16 software. First, mathematical correlations
mum (−1 in coded value) and a maximum (+1 in coded of first order relating the studied responses to factors were
value) as presented in Table 4. The full factorial DOE re- obtained as follows:
sponses were outlet temperature, COD removal and T_out 377.399 + 18.234 T + 9.171 COD + 5.674 Q
maximal temperature in the reactor.
Based on the results the results shown in Figure 7, as + 4.696 T COD + 0.789 T Q − 1.396 COD Q (15)
the temperature limit of a stainless steel reactor was fixed
Y_COD 95.8438 + 0.7588 T + 0.6162 COD − 0.0612 Q
at maximum of 550 °C, it is possible to identify a safety
operation range for the variables shown in Table 4. + 0.5812 T COD − 0.0163 T Q + 0.0112 COD Q (16)
N. Benmakhlouf et al.: Modeling and simulation of hydrothermal oxidation 13
Run Inlet temperature Inlet COD Feed flow Outlet temperature Percentage of COD Maximal temperature
(°C) (gO/L) rate (kg/h) (T_out) (°C) removal (Y_COD) (%) (T_max) (°C)
T_max 432.119 + 55.631 T + 34.119 COD − 0.994 Q by the initial COD and finally the feed flow rate. The in-
teractions were not significant. The Pareto chart of the
+ 32.631 T COD − 0.756 T Q − 0.994 COD Q (17)
standardized effects shown in Figure 9(b) classified all the
The corresponding coefficients of determination were main factors and interactions effects. As confirmed by the
0.994, 0.999 and 0.993 for the outlet temperature (°C), analysis of variance of Table 6 all the factors and their two
percentage of COD removal (%) and maximal temperature level interactions are no significant with a p-value >0.05.
(°C) models respectively. For a heat recovery purpose, the outlet temperature
must be high, thus, the inlet temperature and initial COD
3.5.1.1 Outlet temperature concentration had to be higher, and moreover feed flow
Figure 9 presents the main effects and interactions plots and rate had to be higher, since they have a positive effect on
also the Pareto chart for the outlet temperature response. In the response. Although, their influence do not seems to be
this case and depending in the studied ranges, the initial significant because the increase in temperature and COD
temperature was the most important positive effect followed would not induce very high outlet temperature. This could
Figure 9: Outlet temperature response, (a) Main factors effects, (b) Pareto Chart of the standardized effects, (c) Interactions factors effects.
14 N. Benmakhlouf et al.: Modeling and simulation of hydrothermal oxidation
Table : Analysis of variance for Outlet temperature. Table : Analysis of variance for COD removal.
Term Effect Coef SE Coef T p-value Term Effect Coef SE Coef T p-value
Constant . . . . Constant . . . .
T . . . . . T . . . . .
COD . . . . . COD . . . . .
Q . . . . . Q −. −. . −. .
T*COD . . . . . T*COD . . . . .
T*Q . . . . . T*Q −. −. . −. .
COD*Q . . . . . COD*Q . . . . .
be explained by the fact that the reaction was rapid and shown in Figure 10(b) classified all the main factors. As
therefore mainly took place in these cases in the first two confirmed by the analysis of variance of Table 7, the tem-
meters of the reactor, as can be appreciated in Figure 5. In perature, COD and the interaction of the two factors were
this case the tubular reactor was long enough to give time the most significant factors p-value <0.05. Also, from the
to important heat losses, and finally the outlet temperature t-student values the temperature and the concentration
was quite similar in all the cases. had the most important effect.
Therefore, for an efficient pollutant degradation,
3.5.1.2 COD removal which corresponded to a maximum rate of elimination of
Figure 10 presents the main effects and interactions plots COD, the inlet temperature might be higher.
and also the Pareto chart for the percentage of COD removal
response. For this response, the temperature was also the 3.5.1.3 Maximal temperature
most important effect and it had a positive effect, followed Figure 11 presents the main effects, interactions plots and
by initial COD. On the other hand, the feed flow rate had no Pareto chart for the maximal temperature in the reactor
effect on the pollutant degradation. In this case, only the response. The temperature was found to be the most
interaction between Temperature and COD showed a sig- important positive effect followed by the COD, and finally
nificant effect. The Pareto chart of the standardized effects the lowest effect corresponds to the flow rate. The
Figure 10: COD removal response, (a) Main factors effects, (b) Pareto Chart of the standardized effects, (c) Interactions factors effects.
N. Benmakhlouf et al.: Modeling and simulation of hydrothermal oxidation 15
Figure 11: Maximal temperature response, (a) Main factors effects, (b) Pareto Chart of the standardized effects, (c) Interactions factors effects.
interaction temperature COD was the only significant one. Table : Analysis of variance for Maximal temperature.
The Pareto chart of the standardized effects shown in
Figure 11(b) classified all the main factors and interactions Term Effect Coef SE Coef T p-value
effects, as confirmed by the analysis of variance of Table 8, Constant . . . .
the temperature, the concentration and their interaction T . . . . .
were significant with p-value <0.05, however flow rate and COD . . . . .
Q −. −. . −. .
their interactions were not significant with a p-value >0.05.
T*COD . . . . .
Therefore, the temperature and feed concentration T*Q −. −. . −. .
were critical operating parameters, because of high reac- COD*Q −. −. . −. .
tivity of cutting oil, the temperature increased rapidly at
the first section of the reactor. Thus, feed concentration and
temperature at the reactor inlet might be controlled to means of the overlaid graphic of Minitab and by imposing
avoid exceeding the material temperature limit, and to constraints on these variables.
maintain safe conditions for working plant. Indeed, for operating safety reasons, two maximum
reactor temperature values limited to 500 and 550 °C were
3.5.1.4 Safety zone considered, a minimum treatment efficiency of 96%, and a
Supercritical water oxidation is a severe process in terms of minimum outlet temperature of 375 °C for energy recovery
temperature and pressure, so it is crucial to delimit a safety reasons.
zone where the reactor material would be protected from The obtained safety zone is shown in Figure 12 for a
an eventual dangerous operating environment, particu- flow rate of 17 kg/h, as a white band of the graph delimited
larly towards its structural and mechanical properties by the imposed constraints. As shown in the two graphs of
would be concerned. The obtained models using the full Figure 12, the more severe are the constraints, the more the
factorial design were also exploited to define a safety zone area is reduced. Indeed, from a temperature limit of 550 to
by superposition of the dependent variables limits by 500 °C the safety zone has been clearly reduced.
16 N. Benmakhlouf et al.: Modeling and simulation of hydrothermal oxidation
Li, L., P. Chen, and E. F. Gloyna. 1991. “Generalized Kinetic Model for Wet of an Industrial Wastewater.” Chemical Engineering Journal 144
Oxidation of Organic Compounds.” AIChE Journal 37 (11): 1687–97. (3): 361–7.
Loppinet-Serani, A., C. Aymonier, and F. Cansell. 2010. “Supercritical Sahu, N. K., and A. Andhare. 2018. “Design of Experiments Applied
Water for Environmental Technologies.” Journal of Chemical to Industrial Process.” In Statistical Approaches with
Technology & Biotechnology 85 (5): 583–9. Emphasis on Design of Experiments Applied to Chemical
Marulanda, V. 2011. “Reacting Flow Simulations of Supercritical Water Processes. London, UK: IntechOpen.
Oxidation of PCB-Contaminated Transformer Oil in a Pilot Plant Sánchez-Oneto, J., J. R. Portela, E. Nebot, and E. M. de la Ossa. 2007.
Reactor.” Brazilian Journal of Chemical Engineering 28 (2): 285–94. “Hydrothermal Oxidation: Application to the Treatment of
Miller, A., R. Espanani, A. Junker, D. Hendry, N. Wilkinson, D. Bollinger, Different Cutting Fluid Wastes.” Hazardous Materials 144 (3):
J. M. Abelleira-Pereira, M. A. Deshusses, E. Inniss, and W. Jacoby. 639–44.
2015. “Supercritical Water Oxidation of a Model Fecal Sludge Shrawage, Abhijeet J. 2014. CFD Analysis of Supercritical Water
without the Use of a Co-fuel.” Chemosphere 141: 189–96. Reactor for Flow Back and Produced Water Treatment.
Module Users Guide, COMSOL. CFD. 2015. Comsol Multiphysics 5(2). Ohio: Russ College of Engineering and Technology of Ohio
(n.d.). USA: COMSOL. University.
Moussiere, S., C. Joussot-Dubien, P. Guichardon, O. Boutin, H. A. Turc, Skerlos, S. J. 2013. “Cutting Fluids and Their Environmental Impact.” In
A. Roubaud, and B. Fournel. 2007. “Modelling of Heat Transfer Encyclopedia of Tribology, edited by Q. J. Wang and Y.-W. Chung,
and Hydrodynamic with Two Kinetics Approaches during 655–60. USA: Springer.
Supercritical Water Oxidation Process.” The Journal of Tang, X., S. Wang, L. Qian, M. Ren, P. Sun, Y. Li, and J. Qiao Yang. 2016.
Supercritical Fluids 43 (2): 324–32. “Corrosion Properties of Candidate Materials in Supercritical
Moussiere S. 2006. Etude par simulation numérique des écoulements Water Oxidation Process.” Journal of Advanced Oxidation
turbulents réactifs dans les réacteurs d’oxydation Technologies 19 (1): 141–57.
hydrothermale : Application à un réacteur agité double Vadillo, V., M. B. García-Jarana, J. Sánchez-Oneto, J. R. Portela, and
enveloppe. Paul Cezanne Aix-Marseille III. E. J. Martínez de la Ossa. 2011. “Simulation of Real Wastewater
Nadjiba, B., O. Nawel, and M. Abdesslam Hassen. 2017. “Modeling Supercritical Water Oxidation at High Concentration on a Pilot
and Optimization of Phenol Hydrothermal Oxidation in Plant Scale.” Industrial & Engineering Chemistry Research 50
Supercritical Water.” International Journal of Hydrogen Energy (22): 12512–20.
42 (17): 12926–32. Vadillo, V., J. Sánchez-Oneto, J. R. Portela, and E. J. M. de la Ossa.
NIST. National Institute of Standards and Technology (NIST). (n.d.). 2014. “Supercritical Water Oxidation for Wastewater Destruction
[Text]. Also available at https://www.nist.gov/national-institute- with Energy Recovery.” In Supercritical Fluid Technology for
standards-and-technology (Retrieved August 23, 2019). Energy and Environmental Applications. 181–90. Boston:
Oh, C. H., R. J. Kochan, and J. M. Beller. 1997. “Numerical Analysis and Elsevier.
Data Comparison of a Supercritical Water Oxidation Reactor.” Vadillo M. V. 2012. Estudio del proceso de oxidation en agua
AIChE Journal 43 (6): 1627–36. supercritica para su escalamiento industrial : implantacion de
Patel, J. K., and A. P. M. J. Pandya. 2016. “To Compare Technologies nuevas soluciones tecnologicas, simulation y optimization.
for Treatment of Different Types of Metal Cutting Fluid.” Cadiz: Faculdad de Ciencias, Departamento de Ingenieria
International Journal of Advance Research and Innovative Ideas Quimica y technologia de alimentos, Universidad de Cadiz.
in Education 2 (3): 528–35. Vielcazals, S., J. Mercadier, F. Marias, D. Mateos, M. Bottreau,
Portela, J. R., E. Nebot, and E. M. de la Ossa. 2001. “Generalized F. Canselle, and C. Marraud. 2006. “Modeling and Simulation of
Kinetic Models for Supercritical Water Oxidation of Cutting Oil Hydrothermal Oxidation of Organic Compounds.” AIChE Journal
Wastes.” The Journal of Supercritical Fluids 21 (2): 135–45. 52 (2): 818–25.
Portela, J. R., J. López, E. Nebot, and E. M. de la Ossa. 2001. Voisin, T., A. Erriguible, and C. Aymonier. 2020. “A New Solvent
“Elimination of Cutting Oil Wastes by Promoted Hydrothermal System: Hydrothermal Molten Salt.” Science Advances 6 (17):
Oxidation.” Journal of Hazardous Materials 88 (1): 95–106. eaaz7770.
Sánchez-Oneto, J., F. Mancini, J. R. Portela, E. Nebot, F. Cansell, and Wagner, W., and H.-J. Kretzschmar. 2008. International Steam Tables -
E. J. Martínez de la Ossa. 2008. “Kinetic Model for Oxygen Properties of Water and Steam based on the Industrial
Concentration Dependence in the Supercritical Water Oxidation Formulation IAPWS-IF97, 2nd ed. Verlag: Springer.
Nom et Prénom : BENMAKHLOUF Nadjiba
Titre : Etude du Procédé d’Oxydation Hydrothermale dans l’eau supercritique
pour le Traitement des Eaux Polluées
Thèse en vue de l’obtention du Diplôme de Doctorat en sciences en Génie des
Procédés
Résumé
L’oxydation hydrothermale dans l'eau supercritique est une technique de traitement des déchet, qui consiste à
oxyder des matières organiques en utilisant de l'eau dans des conditions supercritiques. Elle a été développée comme
technique alternative afin de limiter les risques de pollution secondaire.
Le but de cette étude est de développer un outil de simulation à l'état stationnaire, pour le réacteur tubulaire d'oxydation en
Le travail a été réalisé en utilisant une approche de modélisation bidimensionnelle appliquant le modèle de turbulence k-ε
et la formulation de l'Association Internationale pour les Propriétés de l'eau et de la vapeur (IAPWS-IF97) pour calculer
les propriétés de l’eau dans les conditions supercritique. Les propriétés du milieu réactionnel ont été supposées celles du
Le modèle 2D développé en utilisant le COMSOL Multiphysics, a été nettement amélioré par rapport au modèle 1D et a
donc été utilisé pour modéliser le comportement du réacteur supercritique pour la dégradation du polluant modèle : le
phénol et les polluants réels : l’huile de coupe Biocut35 et le déchet de l’huile d’olive.
Les résultats obtenus sont en bon accord avec les valeurs expérimentales, et ont aussi été exploités pour une étude
paramétrique de l’effet des conditions opératoires sur les performances du réacteur, en utilisant la méthodologie des plans
d’expériences.
Le programme développé a permis, pour la première fois, de tracer une zone de sécurité du réacteur supercritique pour
chaque polluant, et pour des conditions opératoires différentes, ce qui assure un fonctionnement pour un maximum
Mots clés : Eau supercritique, Oxydation, Modélisation, réacteur tubulaire, Plan d'expériences, Zone de sécurité.