1 OPERATIONS DE MELANGE

1.1 Mélanges et systèmes d’agitation.

Les systèmes d’agitation permettent le mélange entre phases identiques comme un

mélange de deux liquides ou entre phases différentes comme un mélange entre un
liquide et un solide. Les phases peuvent être miscibles ou immiscibles comme dans
la préparation d’émulsion de deux liquides ou de la dispersion d’un gaz dans un
liquide.
Les systèmes d’agitation sont également utilisés pour optimiser la mise en contact
de réactifs dans le cas de réactions chimiques.
L’efficacité des transferts de chaleurs est également en relation avec le système
d’agitation.

1.2 Systèmes dispersés.

Les systèmes dispersés sont des mélanges de fines particules d’une phase solide,
liquide ou gazeuse non miscibles et dispersées dans une phase dispersante, solide
liquide ou gazeuse. Parmi ces systèmes nous trouvons les émulsions, suspensions
ou mousses.
Les émulsions, suspensions ou mousses sont des mélanges qui macroscopiquement
semblent être des systèmes monophasiques homogènes mais qui
microscopiquement sont des systèmes diphasiques dispersés.

Particules dispersées

Phase dispersante

Vision macroscopique. Vision microscopique.

Figure 1.1 Visions macroscopique et microscopique des émulsions, suspensions ou mousses.

Les particules de la phase dispersée, phase discontinue, sont réparties dans la phase
dispersante, phase continue, comme dans le cas d’une émulsion d’huile dans de
l’eau ou d’eau dans l’huile.

D.Toueix Septembre 2020 1

1.2.1 Système dispersé et état colloïdal.

Lors de la dispersion de la matière au sein d’une phase dispersante on obtient,

suivant les dimensions des particules dispersées, une solution, un état colloïdal ou
une suspension.
Système dispersé

10-9 m 10-6 m

Solution Etat colloïdal Suspension

Figure 1.2 Tailles des particules dans les mélanges.

Les systèmes dispersés colloïdaux sont intermédiaires entre les suspensions,

particules de taille supérieure au micromètre, et les solutions, particules de taille
inférieure au nanomètre.

Un colloïde est donc un mélange d'une ou plusieurs substances, dispersées

régulièrement dans une autre substance, formant un système à deux phases séparées.
Dans un fluide, il forme une dispersion homogène de particules dont les
dimensions peuvent aller du nanomètre au micromètre.

On distingue différentes catégories de colloïdes selon la nature des phases en

présence :

Phase dispersée
Phase dispersante
Gaz Liquide Solide
Gaz totalement Aérosol. Aérosol.
Gaz
miscibles entre Brouillard. Fumée.
eux.
Mousse. Emulsion. Sol.
Liquide
Crème fouettée. Lait. Peinture, encre.
Mousse solide. Emulsion solide. Sol solide.
Solide
Polystyrène Beurre. Verre.
expansé.

Les sols sont des systèmes colloïdaux constitués d’une phase solide dispersée dans
une phase dispersante liquide ou solide.

Les aérosols sont des systèmes colloïdaux constitués d’une phase dispersante
gazeuse et d’une phase dispersée liquide ou gazeuse.

D.Toueix Septembre 2020 2

Les mousses sont des systèmes colloïdaux constitués d’une phase gaz dispersée
dans une phase dispersante liquide ou solide.

Enfin les émulsions sont des systèmes colloïdaux constitués d’une phase liquide
dispersée, phase discontinue, dans une phase dispersante, phase continue, liquide
ou solide.

Les systèmes dispersés sont largement présents dans les produits agroalimentaires
et pharmaceutiques et cosmétologiques.

1.2.2 Etat colloïdale et micelles.

Un état colloïdal peut être caractérisé par la présence de particules dispersées dans
une phase continue dispersante. Ces particules peuvent être des micelles de
dimensions comprises entre 10-9 m et 0,610-7 m. Il peut s’agir de grosses molécules
ou d’agrégats de petites molécules. Ces micelles peuvent êtres des particules solides
cristallisées ou amorphes, gouttes de liquides ou bulles de gaz.

Une micelle est constituée d’un agrégat sphéroïdal de molécules possédant une tête
hydrophile s’orientant vers les phases polaires et une chaîne hydrophobe s’orientant
vers les phases apolaires.

Phase polaire

Chaîne hydrophobe

Phase apolaire

Tête hydrophile
Figure 1.3 Formation de micelles en présence d’une phase polaire ou apolaire.

1.2.3 Stabilité des systèmes dispersés.

Les dispersions sont thermodynamiquement instables. Elles peuvent être stables du

point de vue de la cinétique sur une longue période en particulier lorsque ces
dispersions se trouvent en présence de dispersant comme un émulsifiant dans le cas
des émulsions. L’instabilité des systèmes dispersés est due à des forces
intermoléculaires aux interfaces entre les molécules dispersées et les molécules
dispersantes. Le phénomène est particulièrement important dans le cas de l’eau

D.Toueix Septembre 2020 3

ayant une tension de surface très élevée ainsi que dans beaucoup de milieu bio-
organique dont le lait par exemple.

1.2.3.1 Tension superficielle.

La tension superficielle ou tension interfaciale est une énergie potentielle de

surface ou énergie surfacique. L’énergie interne des molécules d'un fluide
s’exprime vers les molécules voisines constituant une autre phase.

Pour les molécules qui sont dans la

Phase 1 même phase, la résultante totale des
F F forces est globalement nulle, et en
F F
Interface interagissant avec toutes les autres, elles
F F
sont toutes en équilibre.
F F
F F F
A l’interface des deux phases les forces
F
F F F F  
résultantes R1 et R 2 sont la somme des
F F R
forces exercées par l’interaction des
F Phase 2
F différentes molécules.

Figure 1.4 Approche microscopique de la tension superficielle.

Dans le cas spécifique de deux phases, les variables associées sont le volume et la
tension interfaciale. Il est possible d’obtenir une définition de la tension de surface
ou interfaciale (γ) en utilisant le modèle de Gibbs–Duhem généralisée.

dG  VdP  S dT  Ad    ni di (Eq 1.1)

dG, Energie libre est la somme de quatre composantes énergétiques. VdP Energie
due au travail de pression.  S dT Energie due l’entropie du système. Adγ Energie
interfaciale avec A surface de l’interface.  ni di Energie moléculaire avec ni
nombre de moles du composant i et μi potentiel chimique.
A pression, température et composition constante du système.

 G 
dG  Ad    
 A  P ,T ,ni

∆G = A × ∆  (Eq 1.2)

D.Toueix Septembre 2020 4

La tension interfaciale ou de surface, (γ), est une énergie par unité de surface
(J.m−2), qui équivaut aux dimensions à une force par unité de longueur (N.m−1).
C’est cette unité qui est habituellement utilisée pour définir la tension interfaciale
ou de surface.
Pour les interfaces solide/gaz (γ) est de l’ordre de 8 10-2 N·m−1.
Pour les interfaces solide/liquide (γ) est de l’ordre de 2 10-2 N·m−1 comme dans le
cas de la glace et de l’eau.
Pour les interfaces liquide/gaz (γ) est de l’ordre de 10-2 N·m−1 à 10-1 N·m−1.
Par exemple nous avons pour les interfaces :
Eau/air à 20°C (γ) =7,28 10-2 N·m−1. Eau/air à 0°C (γ) =7,56 10-2 N·m−1.
Ethanol/air à 20°C (γ) =2,23 10-2 N·m−1. Eau/huile à 20 °C (γ) = 2 10−2 N·m−1.
(γ) diminue avec la température d’environ 2% tous les 10°C.

1.2.3.2 Différents phénomènes de déstabilisation d’un système dispersé.

Crémage

Floculation
Sédimentation

Rupture de
Phase
Coalescence
Figure 1.5 Mécanisme de déstabilisation des systèmes dispersés.

La déstabilisation d’un système dispersé peut intervenir à la suite d’une

augmentation de la taille des particules ou d’une modification de la densité de la
phase dispersée dans la phase dispersante.
Dans un premier temps il y a floculation les particules de la phase dispersée
formant des agrégats.

D.Toueix Septembre 2020 5

Ces agrégats peuvent modifier la densité de la phase dispersée dans la phase
continue et provoquer, crémage ou sédimentation. Ces phénomènes sont réversibles
l'homogénéité pouvant être facilement rétablie par simple agitation.
Par contre, dans le cas des émulsions, si la floculation est à l’origine de coalescence,
alors il y a rupture complète de phase et irréversibilité.

1.2.4 Stabilisation des systèmes dispersés.

1.2.4.1 Tension interfaciale et stabilité des systèmes dispersés.

La tension interfaciale tend à déstabiliser les suspensions ou à s’opposer à leur

formation.
Considérons la dispersion de 5 cm3 d’une phase discontinue en fines gouttelettes
sphériques de 3,10-4 cm de dans une phase dispersante. Si nous admettons que les 5
cm3 non dispersés représentent une sphère parfaite.

4 3
Surface initiale des 5 cm3 : V    r3  r  3 V  (Eq 1.3)
3 4
3 2
r  3 5 r  1,06 cm . Si = 4  r = 4×3,14×1,062 = 14,14 cm2.
4  3,14

Nombre de gouttelettes dans le système dispersé:

N = volume total / volume d’une gouttelette.

V 5 3
N  ≈4,4.1010 (Eq 1.4)
4
  r 3 4  3,14  3 
4 2

La surface totale de ces sphères est telle que :

2 -4 2 2
Stotale = 4  r N = 4  3,1416 × (3 10 )  4,4 1010  5.10 4 cm (Eq 1.5)

En faisant une dispersion on augmente la surface de séparation entre les phases de

façon considérable, 3500 fois la surface initiale. Cette surface est d’autant plus
grande que les particules dispersées sont fines. On accroît ainsi l'énergie libre du
système, donc son instabilité.
La tension de surface totale, ou mieux l’énergie de surface, a donc été multipliée
par un facteur important conférant ainsi au mélange une énergie interne importante.
C’est un système très instable qui tend à revenir à la situation où coexistent deux

D.Toueix Septembre 2020 6

phases uniques séparées par une zone inter faciale l'énergie étant minimale lorsque
la surface est minimale.
C’est ce qui se passe si on essaie de préparer une vinaigrette constituée
uniquement d’huile et de vinaigre.

1.2.4.2 Les surfactifs et la stabilité des systèmes dispersés.

Les surfactifs, agents de surface ou tensioactifs, ont la capacité généralement de

réduire l’énergie libre interfaciale en diminuant la tension interfaciale totale.
∆G = A × ∆  (Eq 1.2)
Les agents tensio-actifs ont une structure moléculaire caractéristique, connue
comme structure d’amphipathique c’est-à-dire lorsqu'ils possèdent à la fois un
groupe, polaire, hydrophile et un groupe, apolaire, hydrophobe.
<<<<<<<<<
Partie hydrophobe
Partie hydrophile

Figure 1.6 Modèle et exemple d’un agent tensio-actif le dodécyle sulfate de sodium.

1.2.4.3 Divers types de surfactifs.

La nature des groupes de tête et la queue des agents tensio-actifs peut être très
différents. Le groupe hydrophile peut être chargé ou neutre.
-
- Tensio-actif anionique R -SO3

+ Tensio-actif cationique

Tensio-actif non ionique R-O-(CH2CH2O)n-H

la portion tête ne porte aucune charge apparente.
- + Tensio-actif type zwitterionique

Figure 1.6 Modèles et exemples de surfactifs.

Un tensio-actif de type zwitterionique ou amphion est constitué par une molécule

a m p h o t è r e , possédant des charges électriques formelles d’une unité, de signes
opposés et situées en général sur des atomes non adjacents. On trouve ce type de
molécule dans le jaune d’œuf ou le soja comme la « lécitine ».

D.Toueix Septembre 2020 7

 Figure 1.7 Molécule de phosphatidylcholine « lécithine »

Il est intéressant de noter que les agents de surface anioniques et cationiques ne

sont généralement pas compatibles entre eux, tandis que les non ioniques et
zwitterioniques sont compatibles avec tous les autres types.

1.3 Suspensions et sédimentation.

Les suspensions sont des mélanges hétérogènes d’une phase solide dispersée dans
un fluide, phase continue dispersante.
De nombreuses opérations mettent en œuvre des suspensions. C’est le cas dans les
opérations unitaires de cristallisation, de précipitation, de mélanges hétérogènes
solide-liquide et solide-gaz. Dans le domaine de l’environnement c’est le cas des
aérosols, particules fines polluantes dispersées, dans l’atmosphère.

Deux opérations fondamentales sont à étudier dans le cas des suspensions :

- Leur stabilité physique.
- Leur mise en suspension homogène.

1.3.1 Stabilité des suspensions.

Les phénomènes influençant la stabilité physique des suspensions sont la mouillabilité des
particules, la sédimentation ou la flottation et l’homogénéité des particules solides.

1.3.1.1 La mouillabilité.

Lorsque les particules sont recouvertes d’une mince couche d’air, elles ne sont pas
mouillées par la phase dispersante. Ces particules auront tendance à s’agglomérer
et provoquer la sédimentation ou le crémage.
Le mouillage d’un solide par un liquide est un phénomène de surface qui peut être
apprécié par l’angle de contact θ que forment le solide et la tangente à la
gouttelette liquide. Le mouillage est parfait pour θ=0. Il est considéré nul à partir de
θ=90°.
Dans un système tertiaire gaz-solide-liquide le point A est soumis à l’action de trois
tensions interfaciales : liquide-air LA, solide-liquide SL, solide-air SA.

D.Toueix Septembre 2020 8

 A l’équilibre au point (A) le travail
LA d’adhésion entre le solide et le liquide, WSL,
dépend essentiellement de la tension

superficielle du liquide  LA .
SA A SL
WSL   LA  1 cos   (Eq 1.6)

Figure : 1.8 Schéma des tensions interfaciales air-liquide.

Plus la tension superficielle du dispersant sera faible plus la mouillabilité du solide

sera forte. Tout additif, comme un tensio-actif, permettant d’abaisser la valeur de
l’angle  augmente le mouillage du liquide dispersant sur le solide dispersé.

1.3.1.2 Séparation de phase.

Par la différence des masses volumiques du solide et du liquide dispersant les

particules solides de la suspension ont naturellement, sous le champ gravitaire,
tendance à se redéposer ou flotter.
Une taille optimale des particules permet d’obtenir une faible vitesse de sédimentation.
Par exemple dans le cas de compositions pharmaceutiques la taille des particules en suspensions
peuvent aller de 1 à 50 10-6 m pour les suspensions buvables et 3 à 5 10-6 m pour les suspensions
injectables.
La stabilité des suspensions peut être améliorée par augmentation de la viscosité de
la phase dispersante. Toujours dans le cas de produits pharmaceutiques l’utilisation
du glycérol peut, par exemple, permettre d’ajuster la viscosité du milieu.
La modification par ajout d’additif pour modifier la différence de masse volumique
entre la phase solide et liquide est également utilisée pour stabiliser une suspension.
L’idéal est que la masse volumique de la phase dispersante soit égale à la masse
volumique de la phase dispersée.
La concentration de la phase solide influe également sur la vitesse de sédimentation.
Les suspensions concentrées sont en général plus stables que les suspensions
diluées.
La floculation des particules solide permet d’obtenir une structure du sédiment
augmentant sa porosité et diminuant ainsi sa probabilité de redéposition.

1.3.2 Vitesse de chute d’une particule sphérique et sédimentation.

La vitesse de chute d’une particule solide discrète et sphérique est fonction de la résultante des
forces dues à l’attraction terrestre (Fg), à la poussée d’Archimède (Fa) et à la trainée (Ft).

D.Toueix Septembre 2020 9

 4  g  d  (  s  l )
Pour Fa + Ft - Fg = 0 alors u  (Eq 1.7)
3  C x  l

= Volume de la particule.
(Eq 1.8)
= Masse volumique du solide.
Fa
= Accélération due à la gravité.
Ft (Eq 1.9)
= Masse volumique du fluide.
= Coefficient de trainée.
Fg = Surface de la particule.
(Eq 1.10) = Vitesse de sédimentation.
= diamètre de la particule.
Figure 1.9 Vitesse de décantation d’une particule discrète et sphérique dans un fluide statique.

1.3.3 Vitesse de chute et coefficient de trainée.

Le coefficient de trainée (Cx) est fonction du régime d’écoulement créé par la vitesse de chute
d’une particule dans un fluide.

1000
Coefficient de
trainée
100
1 Régime laminaire

10

2 Régime intermédiaire
1 3 Régime
turbulent

0,1
0,1 1 10 100 1000 10000 100000 1000000 10000000
Nombre de Reynolds de grain(Re)

Figure 1.10 Cx de sédimentation d’une sphère en fonction du nombre de Reynolds de grain.

𝑢×𝑑×𝜌𝑙
Le nombre de Reynolds de grain est déterminé par : 𝑅𝑒 = (Eq 1.11)
𝜇𝑙
Avec :
u = Vitesse de sédimentation. d = diamètre de la particule.
 l = Masse volumique du fluide.  l = Viscosité dynamique du fluide.

D.Toueix Septembre 2020 10

L'écoulement est soit laminaire, soit turbulent Re > 1. La chute rapide du
coefficient de traînée correspond à la transition avec une réduction du sillage et
de la traînée.

Laminaire Turbulent Transition

1
2
3

Figure 1.11 Diverses phases de la sédimentation.

1.3.4 Vitesse de chute en régime laminaire.

Dans le cas particulier ou la sédimentation s’effectue en régime laminaire Re ≤ 1

alors le coefficient de trainée est inversement proportionnel au nombre de Re.
24
𝐶𝑥 = 𝑅𝑒 (Eq 1.12)

Dans ce cas la loi de Stokes permet d’apprécier la vitesse de chute d’une particule sphérique dans
un fluide statique.
Avec :
u, vitesse de chute m.s-1.
d, diamètre de la sphère en m.
g, accélération m.s-2
, masse volumique des particules solides en kg.m-3.
(Eq 1.13) , masse volumique des particules solides en kg.m-3.
μ, viscosité du fluide en Pa.s.

1.4 Mélange des systèmes dispersés par mélangeurs discontinus.

1.4.1 Conditions de mise en suspension.

Les mélangeurs discontinus utilisés pour la dispersion de phases sont constitués

d’une cuve, d’un mobile d’agitation avec présence ou non de chicanes.
Les caractéristiques initiales des constituants du mélange, les conditions
opératoires de l’opération et les caractéristiques finales du mélange sont les
données qui orientent le choix d’un mélangeur.

D.Toueix Septembre 2020 11

La géométrie de la cuve associée au mobile d’agitation, la vitesse d’agitation et la
durée du mélange sont des données fondamentales quant aux résultats attendus en
termes d’efficacité, d’homogénéité et de consommation d’énergie.

En général on cherche à combiner action de pompage des matières et mise en

turbulence pour provoquer un brassage des différentes phases dans les mouvements
d’écoulements ainsi créés.

La mise en suspension d’un milieu dispersé dans un milieu dispersant est obtenue par agitation du
système. Les particules solides soumises à des forces gravitaires et ascendantes se trouvent alors
maintenues en suspension dans le fluide par transfert de l’énergie mécanique provoquée par
l’agitation du système.

La vitesse de chute d’une particule sphérique peut permettre d’apprécier la facilité de mettre en
suspension des particules solides dans un fluide. En effet plus la valeur sera faible moins d’énergie
sera consommée pour la mise en suspension et plus la stabilité du mélange sera grande.

1.4.1.1 Vitesse limite minimale de chute d’une particule dans un fluide en milieu turbulent.

Nous pouvons constater que la faciliter à créer une suspension sera d’autant plus facile que la
différence entre les masses volumiques des solides et du fluide sera plus faible, que le rayon des
particules sphériques sera faible et que la viscosité du milieu sera importante.
Pour une particule donnée, soumise à un mouvement turbulent dans un liquide, sa vitesse limite
minimale peut être estimée en considérant qu’il y a dans ces conditions égalité entre l’énergie
cinétique développée dans le liquide et l’énergie potentielle de la particule dans ce liquide.

Energie cinétique dans le fluide Energie potentielle de la particule

l  u 2
l , min = d p  g   s  l
(Eq 1.18) (Eq 1.19)

Avec :
u, vitesse de chute m.s-1.
1/ 2 dp, diamètre de la particule en m.
 d p  g   s  l  g, accélération m.s-2.
u l ,min  



l  s , masse volumique des particules
 
(Eq 1.20) solides en kg.m-3.
l , masse volumique des particules
solides en kg.m-3.

D.Toueix Septembre 2020 12

1.4.1.2 Nombre d’Archimède, faisabilité et stabilité d’un mélange.

Le nombre d’Archimède permet d’évaluer la faisabilité d’obtention d’une suspension solide-fluide.

Le nombre d'Archimède est un nombre sans dimension utilisé pour caractériser le
mouvement d'un solide dans un fluide. C’est le rapport entre les forces
gravitationnelles, les forces d'inertie et les forces visqueuses.
Par expérience plus le nombre d’Archimède est faible plus la faisabilité et la stabilité du mélange
solide-fluide est importante.

d p  g   s  l 
3

Ar  (Eq 1.21)
l   2

Avec :
d, diamètre des particules en m.
g, accélération m.s-2
 s , masse volumique des particules solides en kg.m-3.
l , masse volumique des particules solides en kg.m-3.
μ, viscosité du fluide en Pa.s.

Comme dans le cas de la vitesse de chute d’une particule solide dans un fluide nous pouvons
constater que la faciliter à créer une suspension sera d’autant plus facile que la différence entre les
masses volumiques des solides et du fluide sera plus faible, que le diamètre des particules
sphériques sera faible et que la viscosité du milieu sera importante.

1.4.2 Divers types de mobiles.

Il existe trois familles de mobiles d’agitation se différenciant par le mouvement des

fluides qu’ils engendrent par rapport à l’axe de rotation du mobile :

Type axial Type radial Type tangentiel

D.Toueix Septembre 2020 13

 Figure 1.12 Divers types de mobiles.
1.4.3 Zones d’utilisation des différents types de mobiles.

Ruban
Ancres

Turbines Suivant son type, un mobile est plus

particulièrement adapté pour une
Hélices plage de viscosité déterminée.

10-310-2 10-1 1 10 100 1000 Pa.s

Figure 1.13 Plage d’utilisation de divers mobiles.

1.4.4 Caractéristiques du couple cuve-mobile d’agitation.

Pour assurer un bon écoulement des matières il est important de respecter des
proportions entre le diamètre intérieur de la cuve (Dint), le diamètre de l’agitateur
(dag) et sa hauteur d’implantation (H). Le tableau suivant présente ces principales
caractéristiques.
Type d'agitateur Rapport (dag)/ (Dint)
Agitateur axial 0,2<(dag)/ (Dint)<0,7
Agitateur radial 0,2< (dag)/ (Dint)<0,5
Agitateur raclant 0,95

Pour les mobiles axial ou radial, le rapport entre la hauteur d'implantation du

mobile et le diamètre de la cuve, H/(Dint) est toujours de 1/3.
Si la hauteur de liquide est supérieure à 1,2 H, alors le système est multi-étagé et
comporte deux mobiles ou plus sur l'arbre. La hauteur du fluide dans la cuve doit si
possible être égale au diamètre de la cuve. Parmi les paramètres fonctionnels pour
une cuve et un mobile donnés c’est la viscosité qui conditionne le régime
d’écoulement du fluide dans la cuve et la puissance dissipée.

1.4.5 Paramètres fonctionnels d’un système d’agitation.

L’opération de mélange des systèmes dispersés fait intervenir un certain nombre de

variables comme :
N : la vitesse de rotation de l'agitateur en tours/seconde;
d : le diamètre de l'agitateur en m;
ρ : la masse volumique du fluide en kg.m-3;
P : la puissance dissipée dans le fluide en w;
 : la viscosité dynamique du fluide en Pa.s;

D.Toueix Septembre 2020 14

 
 : la viscosité cinématique;   en en m2.s-1.

Ces différentes variables donnent naissances à des relations permettant d’exprimer
les paramètres fonctionnels souvent sous la forme de nombres adimensionnels.

1.4.5.1 Régime hydrodynamique et nombre de Reynolds.

C’est le nombre permettant de déterminer le régime hydrodynamique d’écoulement

du fluide mélangé. Le nombre de Reynolds exprime le rapport des forces d'inertie
sur les forces de viscosité. Si on a un nombre de Reynolds faible, c'est-à-dire que
les forces de viscosité prédominent, l'écoulement est appelé laminaire. Un nombre
de Reynolds élevé indique que les forces d'inertie prédominent dans l'écoulement,
le régime est dit turbulent.

Lorsque Re < 10 le régime est laminaire.

Lorsque 10 < Re < 104 le régime est

intermédiaire.

Pour les mobiles à débit radial lorsque Re > 104

(Eq 1.22)
le régime est turbulent.

Pour les mobiles à débit axial lorsque Re > 105

le régime est turbulent.
Pour les fortes viscosités les constructeurs de mélangeurs associent
souvent un agitateur radial ou axial à un agitateur raclant.

1.4.5.2 Nombre de puissance ou coefficient de trainée de l'agitateur dans le fluide.

Np est appelé nombre de puissance ou nombre de Newton. Nombre sans dimension,

il caractérise l'énergie transmise au fluide par le système d'agitation.

𝑃
𝑁𝑝 = (Eq 1.23) 𝑃 = 𝑁𝑝 × 𝜌 × 𝑁 3 × 𝑑 5 (Eq 1.24)
𝜌 × 𝑁 3 × 𝑑5

Cette puissance (P) est la puissance transmise au fluide et n’inclut pas la puissance
utilisée au niveau de la partie mécanique.
A titre indicatif, la perte de puissance des roulements et engrenages ainsi qu’au
niveau du moteur compte pour environ 30 à 40% de P.

D.Toueix Septembre 2020 15

A partir de données expérimentales, les courbes de puissance peuvent être établies.
La courbe exprime la variation du nombre de puissance en fonction du nombre de
Reynolds et tous les deux axes sont en coordonnées logarithmiques.
Lorsque le système comporte des chicanes, le nombre de puissance est plus élevé
parce qu'en effet, le rôle principal des chicanes est de casser l'écoulement tangentiel
induit par la rotation. Elles transforment le mouvement tangentiel en un
mouvement axial et c'est pour cela que l'on a une légère augmentation de la
puissance consommée.

Nombre de puissance en fonction du nombre de Re

1000

100
Np 1
Np 2
Np

10
Np 3
1 Np 4
0,1 10 1000 100000 10000000
0,1
Re

Figure 1.14 Diagramme Np en fonction du nombre de Reynolds et du mobile d’agitation.

Mobile Rep. Nb de Nb de H/D d/D

chicanes pales
décollées
Turbine à disque et pales droites. Np1 3 6 1 1/3
Type Rushton
Turbine à pales incurvées. Np2 3 6 1 1/3

Hélice type marine Np3 3 3 1 1/3

Raclant Np4 0 2 1 0,95

Figure 1.15 Nombre de puissance en fonction du nombre de Re pour divers mobiles d’agitation.

Jusqu’à 10 Re, en régime laminaire, le nombre de puissance est proportionnelle au

nombre de Reynolds.
Ces courbes sont aussi utilisables pour les mêmes mobiles mais dans des réservoirs
sans chicanes, si le Nombre de Reynolds est inférieur à 300.

D.Toueix Septembre 2020 16

Pour Re>104, le nombre de puissance est invariant avec Re. On est dans le régime
turbulent, le fluide est peu visqueux. On peut aussi déduire que, dans un système
donné, le nombre de puissance est indépendant de la vitesse de rotation.
Deux facteurs majeurs sont à l'origine d'une consommation énergétique élevée : la
surface entre les pales et le fluide et le profil des pales. Par exemple, la turbine à
disque type Rushton, consomme beaucoup d'énergie parce qu'elle génère un
mouvement du fluide avec un fort phénomène de cisaillement.

Puissance dissipée par unité de volume.

Pour une hauteur de liquide (H) égal au diamètre (D) de la cuve :

𝑃 𝜋×𝐷 2
𝑎𝑣𝑒𝑐 𝑉= × 𝐻 (Eq 1.25)
𝑉 4

La puissance dissipée pour le mélange par unité de volume permet d’évaluer

l’intensité de l’agitation en : Une agitation douce de 0,1 à < 0,4 kW.m-3, une
agitation vigoureuse de 0,4 à < 0,8 kW.m-3 ou une agitation intense de 0,8 à 2
kW.m-3.

Parmi les principales grandeurs caractéristiques des écoulements après le régime

correspondant au nombre de Reynolds on trouve également deux nombres
adimensionnels : le nombre de circulation et le nombre de pompage.

1.4.5.3 Le nombre de pompage.

Le débit de pompage Qp est défini comme le débit de fluide qui traverse le volume
balayé par l'agitateur en rotation. Ce débit de pompage est proportionnel à la
vitesse de rotation, N, du mobile d’agitation et au cube du diamètre du mobile, d,
suivant la relation :

𝑄𝑝
𝑄𝑝 = 𝑁𝑄𝑝 × 𝑁 × 𝑑3 → 𝑁𝑄𝑝 =
𝑁 × 𝑑3
(Eq 1.26) (Eq 1.27)

Le coefficient, NQP, est appelé nombre de pompage. C’est l'expression

adimensionnelle du débit de pompage. Dans le cas d’un régime turbulent ce
nombre de pompage peut être considéré comme constant. Ce nombre est une des
données fournies par le constructeur. A titre d’exemple sont regroupées dans le
tableau 1.2 quelques caractéristiques de mobiles couramment utilisés.

D.Toueix Septembre 2020 17

 Tableau 1.2 Caractéristiques de quelques mobiles usuels.

Mobile Np En régime turbulent NQP

Turbine type Rushton 4,3 à 6 0,75 à 0,83
Turbine 6 pales inclinées à 45° 1,7 0,83
Hélice de type marine 0,3 à 0,35 0,52 à 0,56

1.4.5.4 Le nombre de circulation.

Le débit total de fluide mis en mouvement est supérieur au débit de pompage. Le

débit total est appelé le débit de circulation.
Le débit de circulation Qc est la somme du débit d’entraînement Qe et le débit de
pompage Qp soit Qc=Qe+Qp.
Il est admissible que le rapport Qc/Qp est constant et égal à 1,8. On peut ainsi
définir le nombre de circulation NQc :

𝑄𝑐
𝑄𝑐 = 𝑁𝑄𝑐 × 𝑁 × 𝑑 3 → 𝑁𝑄𝑐 = = 1,8 × 𝑁𝑄𝑝
𝑁 × 𝑑3
(Eq 1.28) (Eq 1.29)

Par corrélation le nombre de circulation et de nombre de pompage permettent

d’estimer le temps de mélange.
De Qp et Qc il est possible, connaissant le volume du mélange V, de calculer les
temps de pompage et de circulation.
𝑉 𝑉
𝑡𝑝𝑜𝑚𝑝𝑎𝑔𝑒 = 𝑄 (Eq 1.30) 𝑡𝑐𝑖𝑟𝑐𝑢𝑙𝑎𝑡𝑜𝑛 = 𝑄 (𝐸𝑞 1.31)
𝑝 𝑐

1.4.5.5 Vitesse d’agitation minimale de mise en suspension et énergie dissipée dans le cas d’un
mélange solide-liquide.

La mise en suspension d’un solide d’un liquide s’effectue dans la majorité des cas par agitation du
liquide et dispersion de solide dans le liquide. L’agitation doit donc permettre aux particules solides
de se disperser dans se redéposer.

h Nous exposerons ici un modèle simplifié d’une cuve à fond plat

de diamètre Dc, d’un mobile d’agitation de diamètre Da et d’une
hauteur du mobile d’agitation par rapport au fond h.
Da
Dc
Figure 1.16 Modèle simplifié d’une cuve à fond plat.

D.Toueix Septembre 2020 18

Le premier modèle corrélé pour déterminer la vitesse minimale de rotation par seconde pour une
mise en suspension d’un solide dans un liquide dans une cuve à fond plat a été proposé par
Zwietering en 1958.
0,2
𝜇 0,1 ×𝑑𝑝 (𝜌𝑠 −𝜌𝑙 ) 0,45
𝑁𝑚𝑖𝑛 = 𝑆𝑍 × × (𝑔 × ) × (∅𝑠 )0,13 (Eq 1.32)
𝐷𝐶0,85 𝜌𝑙

Avec :
N min , en tour(s).s-1.
Dc, diamètre de la cuve en m.
dp, diamètre des particules en m.
g, accélération m.s-2.
 s , masse volumique des particules solides en kg.m-3.
l , masse volumique des particules solides en kg.m-3.
μ, viscosité du fluide en Pa.s.
s , Masse solide/ Masse de suspension.

Cette corrélation se révèle la plus fiable pour des taux de solide < à 5% en masse
des viscosités < à 0,01 Pa.s.
Le paramètre de Zwietering S Z dépend, en outre du rapport du diamètre de la
cuve Dc, au diamètre externe de l’agitateur Da, du rapport du diamètre de la cuve
de Dc, à la distance entre le fond de cuve et l'agitateur h et du type d'agitateur. Ce
paramètre peut prendre des valeurs entre 3 et 15.
Nienow a proposé des tables des valeurs de ce paramètre pour un grand nombre de
géométries de cuve et d'agitateurs.

La puissance par unité de masse PM est déterminée par la relation :

2
4 D  D 
PM   N p ,s  N  Da   a    a  (Eq 1.33)
3 2

  Dc   h 
 s 
avec N p ,s  N p ,s 0  1   (Eq 1.34)
 0,3  N p , 0 
 s 

1.4.5.6 Temps de mélange.

Le temps de mélange tM est la durée nécessaire pour obtenir une composition

homogène à partir du début de l’agitation.
𝐾×𝑅𝑒 𝑥
Pour un système donné, en régime turbulent nous avons : 𝑡𝑀 = (Eq 1.35)
𝑁
La constance K est fonction du mobile utilisé.

D.Toueix Septembre 2020 19

Le temps nécessaire pour réaliser un mélange peut être estimé dans le cas de
matières parfaitement miscibles suivant divers modèles expérimentaux.

La méthode la plus simple est basée, en régime turbulent et sur le diamètre

intérieur de la cuve (d), le diamètre du mobile d’agitation (D) et la vitesse
d’agitation (N) et un type de mobile spécifique. Dans le cas de systèmes à 3 ou 4
chicanes et Re > 104

𝐷 2
6×( 𝑑 )
Pour une hélice marine: 𝑡𝑚é𝑙𝑎𝑛𝑔𝑒 = 𝑁
(𝐸𝑞 1.36)

𝐷 2
4×( 𝑑 )
Pour une turbine de type Rusthon : 𝑡𝑚é𝑙𝑎𝑛𝑔𝑒 = 𝑁
(𝐸𝑞 1.37)

Une autre méthode, toujours en régime turbulent, se base sur le degré de mélange
(A), sur le diamètre intérieur de la cuve (d), le diamètre du mobile d’agitation (D)
et la vitesse d’agitation (N).

𝐴
𝑙𝑛( )
2
𝑡𝑚é𝑙𝑎𝑛𝑔𝑒 = (𝐸𝑞 1.38)
−𝐾
100−𝐷𝑒𝑔𝑟é 𝑑𝑒 𝑚é𝑙𝑎𝑛𝑔𝑒 𝑒𝑛 %
𝐴= (𝐸𝑞 1.39)
100

𝑁 𝑑 2
Pour une hélice marine: 𝐾 = 0,9 × (𝐷) (𝐸𝑞 1.40)

𝑁 𝑑 2/3
Pour une turbine de type Rusthon : 𝐾 = 0,5 × (𝐷) (𝐸𝑞 1.41)

Nous présentons ici une dernière méthode prenant en compte le facteur (ft)
dépendant du régime d’écoulement.

f t  H 1/ 2  D 3 / 2
tM  (Eq 1.42)
N  d 
2 2/3
 g 1/ 6  d 1/ 2

Cette méthode sous évalue le temps de mélange dans le cas d’un régime turbulent.

tmélange proche de t pompage

D.Toueix Septembre 2020 20

Elle présente cependant l’avantage de s’appliquer aux régimes intermédiaire et
laminaire.

ft en fonction du nombre de Re
1000

100
Coefficient ft

10

1
0,1 10 1000 100000
Re

Figure 1.17 Coefficient ft en fonction du nombre de Re.

Pour modéliser les performances de mélange d'une installation donnée on réalise

des expériences en injectant un traceur dans l’installation.

Une corrélation donnant le temps de mélange en relation avec les caractéristiques

du mélangeur peut alors être établie.

Une corrélation proposée par Nienow en 1997 présente, pour des mélangeurs agités
mécaniquement, l’avantage de faire apparaitre dans sa relation le nombre de
pompage comme un des paramètres.

Sonde Traceur % du traceur

𝑁𝑄 1/3 𝐷 2
𝑝
𝑡𝑀 = 5,3 × × (𝑑 ) (Eq 1.43) Temps
𝑁

Figure 1.18 Principe de suivi du temps de mélange.

D.Toueix Septembre 2020 21

1.4.6 Mélanges liquide-gaz.

Les opérations de mélange gaz-liquide dans une cuve agitée provoque la présence
de bulles provoquant une agitation complémentaire. Le plus souvent le gaz ajouté
est un additif indispensable à la mise en œuvre de nombreux procédés comme dans
le domaine du traitement des eaux, des réactions chimiques ou des bio-procédés.

Nous retrouverons les mêmes lois que dans le cas des mélanges liquide-liquide ou
liquide-solide. Cependant la présence des bulles modifie le comportement
hydrodynamique des appareils.

1.4.6.1 Régimes hydrodynamiques dans les mélanges liquide-gaz.

Il existe trois régimes hydrodynamiques dépendant de la vitesse de rotation du mobile

et du débit de gaz injecté dans le réacteur :

- Le régime d’engorgement lorsque les bulles de gaz remontent le long du

mobile d’agitation sans effet de cette agitation.
- Le régime de charge lorsque les bulles de gaz se mélange au-dessus de
l’agitateur.
- Le régime de dispersion lorsque les bulles de gaz sont dispersées dans toute la
cuve.

Agitation croissante

Débit croissant du gaz

Régimes
d’engorgement de charge de dispersion

Figure 1.19 Régimes hydrodynamiques des mélanges liquide-gaz.

Pour un diamètre donné de mobile d’agitation le meilleur régime sera donc un

compromis entre la vitesse de rotation de ce mobile d’agitation et le débit
d’alimentation en gaz.

Ces régimes sont caractérisés par les valeurs de deux nombres sans dimension :

D.Toueix Septembre 2020 22

Le nombre d’aération (𝑁𝑎 ) et le nombre de Froude (𝐹𝑟 ) qui exprime le rapport des
forces d’inerties aux forces gravitaires impliquées dans la formation du vortex.

𝑄𝑔 𝑁 2 ×𝑑
𝑁𝑎 = 𝑁×𝑑3 (Eq 1.44) 𝐹𝑟 = (Eq 1.45)
𝑔

Avec : Qg est le débit du gaz en m3 . s −1 .

N est la vitesse d′ agitation du mobile en m. s −1 .
d est le diamètre du mobile d′ agitation en m.
g est l′ accélération de la pesenteur g = 9,81 m. s −2 .

On observe que la meilleure efficacité est obtenue pour un nombre d’aération >0,1.

Dans le cas où le liquide est de l’eau le régime de dispersion peut s’exprimer par la
relation :
0,5
𝑑 0,5 𝑑 0,5 𝑁 2 ×𝑑
𝑁𝑎 = 0,2 × (𝐷) × 𝐹𝑟 0,5 = 0,2 × (𝐷) × ( ) (Eq 1.46)
𝑔

Par exemple dans le cas d’un fermenteur de diamètre interne égale à D=1 m et
d’une turbine de diamètre d=0,33m on peut estimer la vitesse limite inférieure de
rotation afin d’obtenir une bonne dispersion.
0,5 0,5
𝑑 0,5 𝑁2 × 𝑑 0,33 0,5 𝑁 2 × 0,33
0,2 × ( ) × ( ) = 0,1 ; 0,2 × ( ) ×( ) = 0,1
𝐷 𝑔 1 9,81

Par approximation successive nous obtenons une vitesse limite de 5 tours.s-1 soit
environs 300 rotations par minutes.

1.4.6.2 Puissance en milieu aéré.

La puissance Pg consommée par une cuve aérée est inférieure à la puissance P

consommée par une cuve non aérée de même configuration, agitée à la même
vitesse et renfermant le même liquide.
A partir de nombreuses données expérimentales plusieurs auteurs proposent des
corrélations sous la forme d’une fonction :
𝑃𝑔
𝑓(𝑄𝑔 ) = (1 − ) (Eq 1.47)
𝑃

Pour une cuve standard avec une turbine à disque Cui et Al proposent en 1996 :

D.Toueix Septembre 2020 23

 𝑃𝑔 𝑄𝑔 ×𝑁 0,25 𝑄𝑔 ×𝑁0,25
(1 − ) = 0,52 + 0,62 × ( ) 𝑝𝑜𝑢𝑟 ≥ 0,055 (Eq 1.48)
𝑃 𝑑2 𝑑2

𝑃𝑔 𝑄𝑔 ×𝑁0,25 𝑄𝑔 ×𝑁0,25
(1 − ) = 9,9 × ( ) 𝑝𝑜𝑢𝑟 ≤ 0,055 (Eq 1.49)
𝑃 𝑑2 𝑑2

Avec : Qg est le débit du gaz en m3 . s −1 .

N est la vitesse de rotation du mobile en r. s −1 .
d est le diamètre de l′agitateur en m.
P est la puissance absorbée par le mobile en w sans aération.
Pg est la puissance absorbée par le mobile en w avec aération.

Dans les industries de la fermentation une relation empirique proposée par Michel
et Miller est couramment utilisée pour estimer Pg :

0,45
𝑃 2 ×𝑁×𝑑3
𝑃𝑔 = 0,34 × √𝑛𝑝 × ( ) (Eq 1.50)
𝑄𝑔

np représentant le nombre de pales de la turbine.

𝑃𝑔
Le rapport, 𝑃 , diminue régulièrement avec le nombre d’aération pour devenir
constant. Ce phénomène intervient après que Na soit supérieur à 0,1.
𝑃𝑔
L’hydrodynamique d’une installation peut être étudiée par suivi la variation de 𝑃
en fonction du débit d’air admis à vitesse constante ou de la variation de la vitesse
d’agitation pour un débit constant de gaz.

1.4.6.3 Transfert de matière.

Le coefficient de transfert KLa est lié au système d’agitation, à sa puissance, à la

vitesse d’introduction de l’air et de l’état de dispersion dans le milieu.

Une relation de corrélation permet d’estimer le KLa pour un système donné :

𝑃𝑔 𝛼 𝛽
𝑘𝐿 𝑎 = 𝑘 × (𝜌×𝑉) × (𝑈𝑔 ) (Eq 1.51)

Le coefficient (k) dépend du système utilisé.

𝑃𝑔
est la puissance volumique, en kW.m-3, dissipée par l’agitateur dans le milieu.
𝑉

D.Toueix Septembre 2020 24

(Ug), la vitesse spécifique du gaz, en m.s-1, est le rapport du débit de gaz en m3.s-1
et de la section de la cuve en m2.

(  et  ) dépendent de la dispersion du gaz dans le milieu et de l’agitation dans ce

milieu. Nous trouvons en général (  ) compris entre 0,5 et 0,7 et (  ) entre 0,2 et
0,6.

Dans l’eau le 𝑘𝐿 𝑎 est le coefficient de transfert en d’oxygène. Il est représentatif

de la vitesse de dissolution du dioxygène dans l’eau et est fortement lié à la
température de cette eau.

Middleton, en (1997), a synthétisé de nombreux résultats expérimentaux et

préconise, pour un système air-eau, d’utiliser les valeurs suivantes:

k = 1,2  = 0,7  = 0,6.

Le coefficient volumétrique de transfert de masse côté liquide, kLa peut être

déterminé par une méthode dynamique de désorption-absorption.
Dans un premier temps l’aération de la cuve est remplacée par un bullage modéré
d’azote, sous une agitation vigoureuse. La concentration en oxygène dissous CL
diminue pour atteindre une valeur CL0 proche de zéro. On procède alors de nouveau
à l’aération de la cuve en fixant le débit d’air QG et la vitesse d’agitation N.
La concentration de l’eau en O2 dissous se stabilise à la valeur de saturation CL*,
pour la température de travail.
La phase liquide considérée comme parfaitement mélangée la concentration en
oxygène dissous obéit à l’équation de conservation :

𝑑𝐶𝐿 (𝑡)
= 𝑘𝐿 𝑎[𝐶𝐿 ∗ − 𝐶𝐿 (𝑡)] 𝑞𝑢𝑖 𝑝𝑎𝑟 𝑖𝑛𝑡é𝑔𝑟𝑎𝑡𝑖𝑜𝑛
𝑑𝑡
𝐶𝐿 ∗ − 𝐶𝐿 (𝑡)
→ 𝑙𝑛 [ ∗ ] = −𝑘𝐿 𝑎 × 𝑡 (Eq 1.52)
𝐶𝐿 − 𝐶𝐿 0

1.4.7 Les mousses.

Les mousses trouvent de nombreuses applications dans les domaines de

l’agroalimentaire, de la pharmacie, de la cosmétologie et des matériaux.
Les mousses sont le résultat d’une dispersion d’un gaz dans un liquide à faible
tension superficielle.

D.Toueix Septembre 2020 25

 La mousse est formée de bulles de
gaz entre lesquelles se trouve le
liquide. Les bulles sont plus ou
moins sphériques suivant l’épaisseur
du film. Le canal est formé par
l’intersection de plusieurs films.

Figure 1.20 Base élémentaire d’une mousse liquide.

Le mélange est obtenu par souflage du gaz dans le liquide ou vive agitation du
liquide en présence de gaz. La tension superficielle du liquide peut être modifié par
un agent tensio-actif.
La mousse ainsi formée peut avec le temps évoluer vers un phénomène de
coalescence ou un phénomène de nourissement.
Le phénomène de coalescence est non régulier et génère donc une mousse d’une
structure non homogène.
Dans une mousse le gaz pouvant être plus ou moins soluble dans le liquide qui
l’entoure a la possibilité de diffuser dans les bulles avoisinnantes. Ce phémomène
est le mûrissement et conduit à une diminution globale du nombre de bulles et à
leur grossissement en fonction du temps.

Les mousses ainsi obtenues possèdent donc la particularité d’avoir une faible masse
volumique.

Figure 1.21 Phénomènes de coalescence et de mûrissement.

D.Toueix Septembre 2020 26

Les mousses solides sont obtenues en introduisant des bulles de gaz dans un liquide
en fusion et en solidifiant le tout par refroidissement rapide.

Propriétés des mousses.

Les mousses sont caractérisées par une grande surface gaz/liquide la stabilité du
système étant favorisé par la viscosité de la phase dispersante.

La fraction d’espace poreux dans une mousse est sa porosité. Ø=1- dr.

dr est la densité relative d’une mousse :

m Masse volumique de la mousse.

dr 
 s Masse volumique du solide la composant.

Par exemple, les mousses utilisées dans le cadre de matériaux d’isolation ont une
densité relative, dr, comprise entre 0,05 et 0,2.

1.4.8 Les gels.

Un gel est un réseau tridimensionnel de polymères emprisonnant un liquide. Les

gels, comme les solides, ne s’écoulent pas mais peuvent se déformer, dans certains
cas, sous l'effet de contraintes.

La gélification phénomène de réticulation lie entre eux les polymères indépendants

jusqu'à ce que l'ensemble des polymères forment un réseau tridimensionnel plus ou
moins résistant.

Polymères

Réseau de polymères
emprisonnant le solvant
entre ses mailles.
Solvant

Figure 1.22 Phénomène de gélification par réticulation.

D.Toueix Septembre 2020 27

On distingue généralement deux catégories de gels:

Les gels physiques caractérisés par une réticulation due aux interactions de forces
de Van der Waals réversibles.
Les gels physiques sont en général dépendants de la température et du pH. La
réversibilité des forces de Van der Waals font que les gels physiques sont
thermoréversibles et peuvent se déformer sous l'effet d'une contrainte.

Les gels chimiques caractérisés par une réticulation due aux interactions de
liaisons covalentes irréversibles.
Les points de réticulations des gels chimiques étant des liaisons covalentes ils
empêchent ce dernier de se déformer facilement sous l'effet d'une contrainte. Les
gels chimiques sont aussi thermorésistants.

En pratique n’importe quel liquide ou gaz peut être utilisé comme solvant. Dans le
cas de l’eau on obtient un hydrogel et dans le cas de l’air un aérogel. Un gel
présente, en général, une masse volumique voisine du solvant le constituant.
Dans le cas d’un hydrogel la masse volumique est environ celle de l’eau.

1.4.9 Les aérosols.

Les aérosols sont des mélanges de fines particules liquide ou solide dispersées dans
un gaz représentant la phase dispersante.
La finesse des particules pouvant être inférieures à 2,5 10-6 m et jusqu’à 10 10-6 m
présentent une vitesse de sédimentation négligeable.
Les aérosols sont utilisés dans la préparation de produits cosmétiques et
pharmaceutiques mais aussi, sous forme de sprays, pour des produits d’entretien
ou certaines préparations agroalimentaires.
Les industries polluantes, les fumées issues des combustions de carburants fossiles
et certains phénomènes naturels comme la formation de poussières arrachées au sol
par le vent peuvent être à l’origine d’aérosols.
Il est fondamental de limiter à priori la production d’aérosols et de séparer à
postériori les particules fines des aérosols formés.

D.Toueix Septembre 2020 28

Références et sources.

Site entièrement consacré aux Fondamentaux de la Cristallisation et de la

Précipitation.
IMT Mines Albi ~ Mines Saint-Étienne ~ Mines ParisTech.
https://nte.minesalbi.fr/CristalGemme/co/uc_TempsMelange.html

Process's et un site très complet dédié aux technologies de procédé. Un large

chapitre est consacré à l’agitation des cuves et réacteurs.
http://processs.free.fr/Pages/VersionClassique.php?page=3720

Document p.d.f très complet sur une étude expérimentale du mélange solide-liquide
et caractérisation des suspensions concentrées en milieux visqueux réalisée par
Manon Lessaigne à l’école polytechnique de Montréal.
https://publications.polymtl.ca/1803/1/2015_ManonLassaigne.pdf

Azprocede propose un excellent simulateur dynamique avec cours et exercices

dans le domaine des opérations unitaire du génie des procédés. Un chapitre est
consacré en particulier sur l’agitation.
https://www.azprocede.fr/Cours_GC/agitation_introduction.html

B.E.I Bureau d’études industrielles est spécialisé dans le domaine énergétique et

procédés dont le mélange et les mobiles d'agitation.
http://hmf.enseeiht.fr/travaux/bei/beiep/book/export/html/1518
http://hmf.enseeiht.fr/travaux/bei/beiep/content/g22/caracterisation-mobiles

Ce site regroupe la définition de différentes grandeurs physiques, équations aux

dimensions et unités utilisées entre autre dans le domène du génie des procédés.
Un chapitre comporte les formules physiques en mécanique des fluides.
http://www.formules-physique.com/categorie/1169

Le site regroupe des ressources dans les domaines de la formulation cosmétique et

des émulsions
http://ressources.unisciel.fr/formulation_cosmetique/co/1-3.html

Ce site collaboratif est ouvert à tous. Il concerne le génie alimentaire et les

connaissances qui gravitent autour des sciences de l’alimentation.
https://www.genie-alimentaire.com/spip.php?article286

D.Toueix Septembre 2020 29

Agitation. Mélange - Concepts théoriques de base Techniques de l’ingénieur J3.
Michel ROUSTAN, Jean-Claude PHARAMOND, Alain LINE, 2015.

Agitation. Mélange - Aspects mécaniques Techniques de l’ingénieur J3.

Patrice COGNART, Florent BOUQUET, Michel ROUSTAN, 2015.

Agitation et mélange. Aspects fondamentaux et applications industrielles. Dunod

Catherine Xuereb, Martine Poux, Joël Bertrand, 2006.

Cristallisation et cristallisoirs. Edition ISTE.

Jean-Paul Duroudier, 2017.

Mélangeurs fluide-fluide, fluide-solide, solide-solide. Edition ISTE.

Jean-Paul Duroudier, 2017.

Mixing by agitation of miscible liquids Part I. ELSEVIER.

J.G.van de Vusse, 1955.

Concepts de génie alimentaire (2° Éd.). Procédés associés, applications à la

conservation et transformation des aliments. Lavoisier.
BAZINET Laurent, CASTAIGNE François, 2019.

Génie des procédés alimentaires - 2e édition. Des bases aux applications

technique et ingénierie, Dunod, 2017
Gilles Trystram, Albert Duquenoy, Jean-Jacques Bimbenet.

Émulsions alimentaires et foisonnement

Traité RTA, série Génie des procédés. Lavoisier.
Coordonnateur : LEGRAND Jack Directeur de Collection :
SABONNADIÈRE Jean-Claude. 2013.

D.Toueix Septembre 2020 30

EC 1 Etude de cas.
Opération de décantation à flux horizontal de particules grenues.
Dans une usine d’épuration une opération de décantation a lieu après un traitement
des eaux et une floculation. Les particules en suspension s’accumulent au fond
d’un bassin, d’où on les extrait périodiquement. L’eau ainsi clarifiée est récoltée en
surface par débordement du bassin.
Le bassin de décantation considéré est un bassin de type rectangulaire. Sa longueur
est de 18 m, (L), pour une largeur de 3 m. Sa profondeur est de 1,2 m (H).
Dans ce type de bassins, la condition pour qu'une particule soit retenue et qu'elle ait
le temps d'atteindre le fond doit parcourir dans ce temps, horizontalement la
longueur et verticalement la hauteur du bassin.
A la température de 12°C nous avons pour l’eau à traiter une masse volumique de
1 000,4 kg.m-3 et une viscosité dynamique 1,3 10-3 Pa.s.
A l’une des extrémités on alimente la suspension à un débit de 6 m3.h-1(Qv). Cette
suspension renferme des particules de diamètre comprises de 8 10-6 m et 10-4 m.

La masse volumique d’une particule solide est de 2200 kg.m-3.

L’écoulement du liquide est considéré comme constante sur toute sa section
verticale du bassin de décantation.
Les particules étant grenues elles décantent indépendamment les unes des autres
avec chacune une vitesse de chute constante.

EC 1.1 Diamètre limite d’une particule pour une sédimentation laminaire.

Nous pouvons tirer des relations (Eq 1.10) et (Eq 1.12) une relation permettant
d’estimer le diamètre limite de particule pour une sédimentation laminaire.

𝑢×𝑑×𝜌𝑙 𝑅𝑒×𝜇𝑙
𝑅𝑒 = 𝜇𝑙
𝑢= 𝑑×𝜌𝑙
(Eq 1.14)


d 2  g   s  l  =
Re l
d3 
Re l  18  l
(Eq 1.15)
18  l d  l l  g  (  s  l )

D.Toueix Septembre 2020 31

Nous prenons le nombre de Reynolds égal à 1 donc :

d 3
Re l  18  l
3
 2
1  1,3.103  18
 0,0001372
l  g  (  s  l ) 1000,4  9,81 (2200  1000,4)

d≡1,37 10-4 m.

EC 1.2 Caractéristiques du bassin de décantation.

ul
us
Eau à clarifier
Eau clarifiée

Boues

Figure EC 1.1.1 Schéma d’une décantation à flux horizontal de particules grenues.

La section transversale du bassin est égale à S = 3 × 1,2 = 3,6 m2.

Le volume du bassin de décantation est égal à V= 3,6 × 18 = 64,8 m3.

La vitesse d’une particule entrant dans le bassin à son niveau supérieur a deux
composantes :

ul : vitesse horizontale du fluide égale à Qv/S.

us : vitesse verticale de chute donnée par la loi de Stokes.
uH = Qv/S est appelée vitesse de Hazen ou charge hydraulique superficielle.

Exprimée en m.h–1, uH est indépendante de la profondeur du bassin.

Dans notre cas uH = 6/3,6 = 1,66 m.h-1 soit environ 4,63 10-4 m.s-1.
Toutes les particules ayant des vitesses de sédimentation supérieures à uH seront
théoriquement éliminées.

L'efficacité d'un décanteur horizontal ne dépend, en théorie, que de la vitesse de

Hazen et non de la hauteur ou du temps de rétention (soit de 0,5 à 1,5 m/h).

Le temps de séjour dans le bassin est donné par : t= V/Qv. t= 64,8/6= 10,8 h.

D.Toueix Septembre 2020 32

EC 1.3 Sédimentation sous flux horizontal des particules grenues.

Pour qu’une particule sédimentant se retrouve au fond du bassin au niveau de la

sortie de l’eau clarifiée cette particule doit parcourir dans le même temps la
longueur et verticalement la hauteur du bassin. Sa vitesse de sédimentation doit être
telle que :

H L H  ul
  us  (Eq 1.16)
us ul L

us  18  l
Le diamètre de la particule sera alors égale à : d  (Eq 1.17)
g   s   l 

H  ul 1,2  4,63 104

us    3,08 105 u s = 3,08 105 m.s 1
L 18

us  18  l 3,08 105  18  1,3 103

d   7,82 106 m.
g   s  l  9,81 2200  1000,4

Toutes les particules supérieures à 7,82 106 m. de diamètre seront retenues dans le
bassin. Et pour garantir une sédimentation en régime laminaire le diamètre limite
supérieur des particules sera de 1,37 10-4 m.
Dans ces conditions la vitesse de sédimentation de la particule sera suivant
(𝐸𝑞 1.12) :
𝑑 2 × 𝑔 × (𝜌𝑠 − 𝜌𝑙 ) 1,37 2 × 𝑔 × (𝜌𝑠 − 𝜌𝑙 )
𝑢= 𝑢=
18 × 𝜇 18 × 𝜇

u

d 2  g   s  l


(1,37 104 )2  9,81 (2200  1000,4)
 9,53 103 m.s 1
18   18  0,0013

H L H  ul 1,2  1,66
Soit environ 34,33 m.h-1. Alors  L   0,058 m
us ul us 34,33

Le temps pour déposer au fond du bassin cette particule sera de :

longueur parcourue  temps pour la particule fine 0,058  10,8

t   0,0348 h
Longueur du bas sin 18

D.Toueix Septembre 2020 33

Diamètre de particule déposée au bout une heure de sédimentation.

temps  Longueur du bas sin 1 18

L   1,66 m
temps pour la particule fine 10,8
H L H  ul 1,2  1,66
  us    1,2 m.h 1  3,33.104 m.s 1
us ul L 1,66
us  18  l 3,33 104  18  1,3 103
d   2,56 105 m.
g   s  l  9,81 2200  1000,4

Figure EC 1.1.2 Graphe de décantation à flux horizontal de particules grenues.

D.Toueix Septembre 2020 34

Etude de cas EC 2
Potabilisation de l’eau de mer.

Dans des régions en déficit d’eau douce et proche de la mer une solution pour
obtenir une eau destinée à la consommation humaine est la distillation de l’eau de
mer suivie d’une reminéralisation.

Pour une eau distillée qui ne requiert qu’une reminéralisation, il est possible
d’introduire du CO2 et de l’eau de chaux saturée au moyen de mélangeurs en ligne
en réglant le pH. Le CO2 utilisé est en général celui récupéré lors de la distillation.
Cette reminéralisation s’effectue suivant l’équation de la réaction.

2CO2+Ca (OH)2 = Ca (HCO3)

Eau distillée

M1
Ca(OH)2 CO2

M2

D Eau potable

R1

M3 Hypochlorite de sodium

R2

Eau potable reminéralisée

R3

Figure EC 1.2.1 Schéma de procédé simplifié d’une unité de reminéralisation d’eau distillée.

D.Toueix Septembre 2020 35

Rep Désignation Fonction
M1 Mélangeur statique Mélange de l’eau distillée et du dioxyde de
carbone
M2 Mélangeur Préparation du lait de chaux
D Décanteur, saturateur Séparation des boues et de l’eau de chaux
R1 Réservoir Stockage de l’eau de chaux
M3 Mélangeur statique Mélange de l’eau de chaux et de l’eau
carbonatée
R2 Réservoir Stockage de la solution d’hypochlorite de
sodium
R3 Réservoir Stockage tampon de l’eau potable
reminéralisée

La capacité de traitement de l’unité est de 1700 m3.h-1 d’eau distillée. Le dioxyde

de carbone est ajouté sous une pression de 120 kPa à la température de 20°C dans
le mélangeur M1. L’eau de chaux est, elle, ajoutée dans le mélangeur M2 de façon
à obtenir une au reminéralisée contenant 17,7 mg.L-1 de HCO3- et 5,8 mg.L-1 de
Ca++. L’ajout d’une quantité adaptée d’hypochlorite de sodium permet d’assurer la
potabilité de l’eau ainsi reminéralisée.
L’eau de chaux stockée dans R1 provient du décanteur D permettant de la séparer
des matières en suspension issus du mélange du lait de chaux et de l’eau potable.
Le lait de chaux est lui préparé dans le mélangeur M3 par ajout d’hydroxyde de
calcium dans de l’eau potable.

EC 2.1 Traitement au dioxyde de carbone.

Le dioxyde de carbone dissout dans l’eau doit réagir dans un second temps avec
l’hydroxyde de calcium suivant l’équation de la réaction.

2 CO2 + Ca(OH)2 = 2 HCO3- + Ca2+

Pour obtenir une eau renfermant 17,7 mg.L-1 de HCO3- soit 17,7 g.m-3 le flux de
dioxyde de carbone en moles.h-1mélangé à l’eau sera de :

FEau  CHCO  1700  17,7

FCO2  3
  493,3 moles.h 1
M HCO  61
3

D.Toueix Septembre 2020 36

Ce qui dans les conditions opératoires représentera en utilisant la loi des gaz
parfaits :

n  R  T 493,3  8,32  293

P  V  n  R T  V    10 m3 .h 1
P 120000

Le débit massique de dioxyde de carbone ajouté au flux des 1700 m3.h-1 d’eau
distillée est d’environ 30 kg.h-1. La variation de débit volumique entre l’entre
l’entrée et la sortie du mélangeur M1 sera donc négligeable.

EC 2.2 Traitement à l’eau de chaux

La solubilité à 20°C de l’hydroxyde de calcium dans l’eau est de 1,73 g.L-1.

F2 Eau de chaux CF2,Ca++

F1 Eau + CO2
M3
F3 Eau reminéralisée CF3,Ca++
Figure 1.2.2 Schéma de principe du traitement au lait de chaux.

L’eau de chaux issue du lait de chaux renferme une mole de Ca2+ par mole de
Ca(OH)2.
mCa (OH )2  M Ca   1,73  40
CF 2 ,Ca      0,935 g.L1  935 g.m 3
M Ca (OH ) 2 74
.
L’eau reminéralisée doit contenir 5,8 mg.L-1 de Ca2+ soit 5,8 g.m-3.

F1 + F2 = F3
F2 × CF2,Ca++ = F3 × CF3,Ca++
1700 + F2 = F3
F 2  935
F2 × 935 = F3 × 5,8 F3   161,21 F 2
5,8
1700
1700 + F2 = 161,21 × F2 F2   10,6 m3 .h 1
160,2

D.Toueix Septembre 2020 37

EC 2.3 Préparation du lait de chaux.

Le lait de chaux se prépare dans une cuve de mélange à partir de di hydroxyde de

calcium dans de l’eau.
Chaux
Lait de chaux Ca(OH)2
Eau de chaux 5/100 Ca(OH)2 18,335 kg.h-1
10,6 m3.h-1 M2
Ca2+ 935 g.m-3
Eau traitée + CO2 D 348,365 kg.h-1
1700 m3.h-1
M3
Impuretés Eau reminéralisée
Eau traitée 10,6 m3.h-1
1700 m3.h-1 Ca2+ 5,8 g.m-3
-3
Ca2+ 5,8 g.m-3 HCO3- 17,7 g.m
HCO3- 17,7 g.m-3
Figure 1.2.3 Schéma de principe de la recalcification de l’eau.

L’hydroxyde de calcium de masse volumique de 2240 kg.m-3 et d’une

granulométrie moyenne de 200 microns est introduit avec un taux d’excès de 1,1.
Soit un débit de 20,17 kg.h-1 pour 383,23 kg.h-1 d’eau.

La cuve utilisée pour réaliser le mélange à un

diamètre intérieur (D) de 0,8 m. Pour un
volume utile de 0,4 m3. Le mobile utilisé est
une hélice de type marine de diamètre (d)
égal à 0,264 m et placé à une hauteur (J) de
C H
0,264 m du fond de cuve. Cette cuve est en
outre équipée de 3 chicanes.
d J

D
Figure 1.2.4 Caractéristique de la cuve de mélange eau-chaux.

D.Toueix Septembre 2020 38

Dans ces conditions la vitesse de rotation minimale du mobile d’agitation pour une
mise en suspension est déterminée à l’aide de la relation de Zwietering pour Sz = 7.

0,45
𝜇 0,1 × 𝑑𝑝0,2 (𝜌𝑠 − 𝜌𝑙 )
(Eq 1.32) 𝑁𝑚𝑖𝑛 = 𝑆𝑍 × × (𝑔 × ) × (∅𝑠 )0,13
𝐷𝐶0,85 𝜌𝑙
∅𝑠 représentant le rapport massique du solide dans le mélange nous avons 5/100.

0,45
0,0010,1 × 0,00020,2 (2240 − 1000)
𝑁𝑚𝑖𝑛 = 7× × (9,81 × ) × (0,05)0,13
0,80,85 1000

𝑁𝑚𝑖𝑛 ≈1,6 rotation par seconde. Soit 96,6 rotations par minute.

En utilisant la corrélation proposée par Nienow, (Eq 1.43), nous pouvons estimer le
temps nécessaire pour obtenir la mise en suspension des particules de Ca(OH)2.
Dans ces conditions le nombre de pompage 𝑁𝑄𝑝 = 0,56.

𝑁𝑄 1/3 𝐷 2 0,561/3 0,8 2

𝑝
𝑡 = 5,3 × 𝑁
× (𝑑 ) 𝑡 = 5,3 × 1,6
× (0,264) ≈ 25 𝑠𝑒𝑐𝑜𝑛𝑑𝑒𝑠

EC 2.4 Opération de décantation

L’opération de décantation est réalisée dans un bac de décantation à section

rectangulaire. La sédimentation est assimilée à la décantation de particules grenues
les plus fines étant de 5 10-6 m. Les dimensions du décanteur sont respectivement
une hauteur de 0,5 m, une largeur de 1,5 m et une longueur de 9 m.
Le débit d’alimentation global dans le décanteur est de 10,6 m3.h-1. L’opération
étant réalisées à 25°C la masse volumique du liquide est estimée à 997 kg.m-3 et sa
viscosité dynamique à 0,00094 Pa.s.

EC 2.4.1 Efficacité du décanteur.

La vitesse horizontale de déplacement du liquide, Qv/S, permet de vérifier les

bonnes conditions de la sédimentation.

𝑄𝑣 𝑄𝑣 10,6
𝑢𝑙 = = = ≈ 0,785 𝑚. ℎ−1
𝑆 𝐿 × 𝑙 1,5 × 9

D.Toueix Septembre 2020 39

Cette valeur étant comprise entre 0,5 et 1,5 m.h-1 la sédimentation peut s’effectuer
dans de bonnes conditions.
Le temps de séjour dans ce décanteur sera de :

𝑉 𝐻 × 𝐿 × 𝑙 0,5 × 1,5 × 9
𝜏= = = ≅ 0,64 ℎ𝑒𝑢𝑟𝑒
𝑄𝑣 𝑄𝑣 10,6

Soit approximativement de 38 minutes.

EC 2.4.2 Diamètre des particules les plus fines retenues.

La vitesse de sédimentation dépend de la vitesse horizontale. Les particules les plus

fines retenues seront alors fonction de la vitesse de sédimentation de la viscosité du
milieu ainsi que de la différence entre la masse volumique du solide et celle du
liquide.

𝐻 × 𝑢𝑙 0,5 × 0,785/3600
(Eq 1.16) 𝑢𝑠 = = = 1,2 10−5 𝑚. 𝑠 −1
𝐿 9

𝑢𝑠 × 18 × 𝑢𝑙 1,2 10−5 × 18 × 9,4 10−4

(Eq 1.17) 𝑑=√ = = 4 10−6 𝑚.
𝑔 × (𝜌𝑠 − 𝜌𝑙 ) 9,81 × (2240 − 997)

La contrainte des 5 10−6 𝑚 pour les particules les plus fines est donc respectée.

D.Toueix Septembre 2020 40

Etude de cas EC 3
Etude hydrodynamique et transfert de matière dans un
bioréacteur pilote.
Le bioréacteur pilote est constitué d’une cuve de 0,55 m de diamètre interne pour
une capacité de 120 L utiles. La cuve dispose de 3 chicanes standards et d’un
agitateur de type turbine 6 pales de 0,18 m de diamètre.

EC 3.1 Etude de l’hydrodynamique.

Cette étude est menée le bioréacteur étant rempli de 120 litres d’eau. La vitesse
d’agitation fixée à 150 rotations par minutes on alimente progressivement de l’air
au débit de 0 à 6,5 m3.h-1. L’évolution du rapport Pg/P en fonction du débit
d’alimentation permet d’observer la diminution de ce rapport jusqu’au point de
charge à partir duquel le rapport reste constant.
Par exemple pour 150 rotations par minutes et un débit d’air de 5,5 m3.h-1 nous
avons suivant (Eq 1.48) :

𝑄𝑔 × 𝑁 0,25 5,5/3600 × (150/60)0,25

= ≈ 0,059 ≥ 0,055 𝑑𝑜𝑛𝑐
𝑑2 0,182

𝑃𝑔 𝑄𝑔 × 𝑁 0,25 𝑃𝑔
(1 − ) = 0,52 + 0,62 × ( ) = 0,556 → = 0,44
𝑃 𝑑2 𝑃

𝑄𝑔 5,5/3600
(Eq 1.44) 𝑁𝑎 = 𝑁×𝑑3 = 150/60×0,183 = 0,104

Ce nombre d’aération Na permet une bonne condition de dispersion.

1,2
1
Rapport Pg/P

0,8
0,6
0,4
0,2
0
0 2 4 6 8
Débit d'air en m3/h

Figure 1.27.a Rapport Pg/P en fonction du débit d'air

D.Toueix Septembre 2020 41

 0,15

0,1

Na
0,05

0
0 2 4 6 8
Débit d'air en m3.h-1

Figure 1.27 Etude hydrodynamique d’un mélange eau-air. Na fonction du débit d'air

.
Nous pouvons observer qu’au point de charge, pour un débit d’air supérieur à
5m3.h-1, le nombre d’aération est supérieur à 0,1 donc en accord avec des bonnes
conditions de dispersion de l’air dans l’eau.

EC 3.2 Transfert d’oxygène dans l’eau.

La cuve est emplie d’eau à 25°C avec 𝜌 = 996,7 kg.m-3. On fixe une vitesse
d’agitation constante à 150 rotations par minutes avant d’introduire de l’azote dans
le milieu. Lorsque la concentration en oxygène atteint 0,050 mole.m-3 le débit d’air
est fixé à 5,5 m3.h-1 et l’on suit en fonction du temps l’évolution de la concentration
jusqu’à saturation.

Figure 1.28 Evolution de la concentration en oxygène dans l’eau. Cuve de 120 L,

150 rotations par minutes 5,5 m3.h-1 d’air.

D.Toueix Septembre 2020 42

 0
0 50 100 150 25°C
-0,2 Temps en secondes

-0,4 Concentration en O2
−𝑘𝐿 𝑎 × 𝑡
-0,6 En mole.m-3
-0,8 y = -0,0082x 𝐶𝐿 ∗ = 0,277
-1 R² = 0,993 𝐶𝐿 0 = 0,05
-1,2 𝐶𝐿 ∗ −𝐶𝐿 (𝑡)
𝑙𝑛 [ 𝐶 ∗ ]
𝐿 −𝐶𝐿 0

Figure 1.29 Détermination du KLa

𝑃 𝛼 𝛽
𝑘𝐿 𝑎 = 𝑘 × (𝑉) × (𝑈𝑔 ) (Eq 1.51) et 𝑃 = 𝑁𝑝 × 𝜌 × 𝑁 3 × 𝑑 5 (Eq 1.24)

4 5,5
V = 0,120 m3 et 𝑈𝑔 = 3,14 ×0,552 × 3600 ≈ 0,0064 𝑚. 𝑠 −1

Le régime étant considéré comme turbulent 𝑁𝑝 reste constant et dans les

caractéristiques du système égal à 6,20.

𝑃𝑔
𝑃 = 6,20 × 996,7 × 150/603 × 0,185 ≈ 19 𝑊 ; = 0,44 𝑃𝑔 = 8,36 𝑊
𝑃

𝑃𝑔 𝛼 𝛽
𝑘𝐿 𝑎 = 𝑘 × (𝑉×𝜌) × (𝑈𝑔 ) (Eq 1.51) 𝑎𝑣𝑒𝑐 k = 1,2  = 0,7  = 0,6.
0,7
8,36
𝑘𝐿 𝑎 = 1,2 × ( ) × (0,0064)0,6 ≈ 0,0089
0,120 × 997

D.Toueix Septembre 2020 43