Étude d'un Précipitateur Électrostatique (ESP)
📝 Situation du Projet
Vous êtes Ingénieur en Génie des Procédés au sein du bureau d'études "PureAir Solutions". Votre client, la cimenterie "Cimenterie du Nord", fait face à un durcissement des normes environnementales (Arrêté Préfectoral de Mise en Demeure). L'usine rejette actuellement des fumées chargées en particules fines (cendres volantes et poussières de clinker) issues du four rotatif principal. Ces particules, d'un diamètre moyen critique, échappent aux cyclones mécaniques actuels.
La direction technique de l'usine vous mandate pour concevoir et dimensionner un Précipitateur Électrostatique (ESP - Electrostatic Precipitator) à plaques parallèles. Votre objectif est de garantir un taux de captation suffisant pour respecter le seuil de rejet légal, tout en optimisant la consommation énergétique de la haute tension. L'enjeu est double : assurer la conformité réglementaire pour éviter la fermeture administrative et récupérer la matière première (poussières) pour la réinjecter dans le process.
En tant qu'expert en électrostatique appliquée, vous devez dimensionner la surface de collecte totale de l'électrofiltre. Vous devrez calculer le champ électrique, la charge acquise par les particules, leur vitesse de migration, et enfin déterminer la géométrie nécessaire pour atteindre l'efficacité cible de 99,5%.
"Attention, nous travaillons ici avec des tensions de l'ordre de 50 000 Volts (DC Négatif). L'Effet Corona génère de l'ozone. Dans vos calculs, ne négligez pas la permittivité des particules, car c'est elle qui conditionne la charge limite de saturation !"
L'étude repose sur les caractéristiques physico-chimiques des fumées sortant du préchauffeur et sur les contraintes géométriques du site d'implantation.
📚 Référentiel Normatif & Physique
Loi de CoulombLimite de PauthenierÉquation de Deutsch-Anderson| FLUIDE PORTEUR (GAZ DE COMBUSTION) | |
| Débit volumique total | 50 m³/s |
| Température de service | 150 °C |
| Viscosité dynamique du gaz (\(\eta\)) | 2.4 × 10⁻⁵ Pa.s |
| AÉROSOL (PARTICULES) | |
| Diamètre moyen des particules (\(d_{\text{p}}\)) | 2.5 µm |
| Permittivité relative des particules (\(\varepsilon_{\text{r}}\)) | 5.0 (Adimensionnel) |
| Facteur de Cunningham (\(C_{\text{c}}\)) | 1.0 (Négligé ici) |
📐 Géométrie & Électricité
- Tension Appliquée (\(V_0\)): 50 000 V (50 kV)
- Distance inter-plaques (\(2b\)): 30 cm (soit \(b = 15\) cm)
- Permittivité du vide (\(\varepsilon_0\)): 8.854 × 10⁻¹² F/m
⚖️ Objectif de Performance
| Variable | Symbole | Valeur | Unité SI |
|---|---|---|---|
| Tension | \(V_0\) | 50 000 | Volts [V] |
| Demi-distance inter-électrodes | \(b\) | 0.15 | Mètres [m] |
| Diamètre particule | \(d_{\text{p}}\) | 2.5 × 10⁻⁶ | Mètres [m] |
E. Protocole de Résolution
Pour dimensionner correctement cet équipement, nous allons suivre une démarche physique séquentielle : du champ électrique jusqu'à la surface de collecte.
Calcul du Champ Électrique Moyen
Déterminer l'intensité du champ électrostatique régnant entre le fil et la plaque, moteur de la migration.
Détermination de la Charge de Saturation
Calculer la charge électrique maximale (\(q_{\text{s}}\)) qu'une particule peut acquérir par bombardement ionique (Limite de Pauthenier).
Vitesse de Migration (Vitesse de Drift)
Établir le bilan des forces (Coulomb vs Traînée de Stokes) pour trouver la vitesse terminale de la particule vers la plaque.
Dimensionnement de la Surface de Collecte
Utiliser l'équation de Deutsch-Anderson pour calculer la surface de plaques nécessaire pour atteindre 99.5% d'efficacité.
Étude d'un Précipitateur Électrostatique (ESP)
🎯 Objectif Scientifique
L'objectif primordial de cette première étape est de quantifier l'intensité du champ vectoriel électrostatique (\(E\)) qui règne au sein de l'espace inter-électrodes. Ce champ est l'initiateur de tout le processus : c'est lui qui va provoquer l'ionisation du gaz (effet couronne autour du fil) et c'est encore lui qui exercera la force de Coulomb sur les particules une fois chargées. Sans une estimation précise de ce champ, aucun dimensionnement mécanique n'est possible.
📚 Référentiel Théorique
Théorème de Gauss Électrostatique FondamentaleDans la réalité industrielle, la géométrie fil-plaque crée un champ électrique très hétérogène : il est extrêmement intense près du fil (ce qui déclenche l'effet couronne) et diminue en s'approchant de la plaque. Cependant, pour un pré-dimensionnement de niveau Avant-Projet Détaillé (APD), l'intégration rigoureuse de l'équation de Poisson est trop complexe. Nous ferons donc l'hypothèse simplificatrice, mais conservative, d'un champ électrique moyen équivalent. Cette approximation est courante dans l'industrie pour estimer la vitesse de migration moyenne.
Le champ électrique \(E\) (exprimé en Volts par mètre) dérive du potentiel électrique \(V\). Dans une configuration de condensateur plan idéal, le champ est uniforme et vaut \(E = V/d\). Ici, bien que la géométrie soit cylindrique-plane, nous utiliserons la grandeur moyenne \(E_{\text{moy}}\) sur la distance \(b\) séparant le fil de la plaque.
📊 Visualisation : Profil du Champ Électrique \(E(x)\)
Note : L'aire sous la courbe rouge correspond au potentiel \(V_0\). L'approximation \(E_{moy}\) (ligne bleue) conserve cette même aire (rectangle) pour simplifier le calcul.
Démonstration de la formule simplifiée : Le champ électrique dérive du potentiel :
En géométrie 1D simplifiée (condensateur plan) :
En intégrant entre le fil (0) et la plaque (b) :
Si on suppose \(E\) constant (moyenne), alors \(E_{\text{moy}} \times b = V_0\), d'où :
Étape 1 : Données d'Entrée
| Paramètre | Valeur |
|---|---|
| Tension appliquée (\(V_0\)) | 50 000 V |
| Distance fil-plaque (\(b\)) | 0.15 m |
Attention aux unités ! La tension est souvent donnée en kV (kilovolts) et les distances en cm ou mm. Pour obtenir un champ en V/m (Système International), il est impératif de tout convertir : kV en V et cm en m AVANT le calcul.
Calcul du Champ Électrique Moyen
Nous appliquons le rapport simple de la tension sur la distance pour obtenir l'ordre de grandeur du champ moteur.
Application Numérique :Interprétation : Nous obtenons un champ d'environ 3.3 kV/cm. C'est une valeur cohérente pour un ESP industriel. En dessous de 2 kV/cm, l'effet couronne est instable. Au-dessus de 4-5 kV/cm, on risque le claquage diélectrique (arc électrique) qui mettrait l'installation en défaut.
L'ordre de grandeur est validé. Dans l'air sec, le champ disruptif est de 30 kV/cm (3 000 000 V/m). Nous sommes à 3.3 kV/cm, soit un facteur de sécurité d'environ 10 vis-à-vis du claquage théorique, ce qui est correct compte tenu de la poussière et de la température qui réduisent la tenue diélectrique.
Ce calcul simpliste du "champ moyen" ignore deux phénomènes critiques :
1. Hétérogénéité Géométrique : Le champ est en réalité hyperbolique, très intense près du fil (dépassant le seuil d'ionisation de l'air) et plus faible vers la plaque.
2. Charge d'Espace : Le nuage d'ions et de particules chargées modifie le potentiel local (Équation de Poisson), ce qui tend à uniformiser le champ en régime établi.
🎯 Objectif Scientifique
Une particule neutre ne subit aucune force dans un champ électrique statique. Pour être collectée, elle doit acquérir une charge électrique (\(q\)). L'objectif ici est de calculer la charge maximale qu'une particule sphérique de diamètre \(d_{\text{p}}\) peut acquérir lorsqu'elle traverse le champ d'ions créé par l'effet couronne. C'est ce qu'on appelle la "charge par champ".
📚 Référentiel
Loi de Pauthenier Mécanisme de charge par champPour des particules supérieures à 1 µm (ce qui est notre cas avec \(d_{\text{p}} = 2.5\) µm), le mécanisme prédominant est la "charge par champ" (bombardement ionique). Les lignes de champ électrique sont déformées par la particule, attirant les ions vers elle. Au bout d'un certain temps (très court, < 0.1s), la charge accumulée sur la particule repousse les nouveaux ions arrivants : on atteint la saturation. Nous devons calculer cette limite théorique.
Il existe deux modes de charge : la "charge par diffusion" (mouvement thermique des ions, prédominant pour \(d_{\text{p}} < 0.2 \mu m\)) et la "charge par champ" (prédominant pour \(d_{\text{p}} > 1 \mu m\)). La formule de Pauthenier modélise uniquement la charge par champ. Elle montre que la charge limite est proportionnelle au carré du diamètre de la particule et à l'intensité du champ électrique.
⚛️ Visualisation : Mécanisme de Bombardement Ionique
Genèse de l'équation de Pauthenier :
1. Une particule diélectrique dans un champ \(E_0\) déforme les lignes de champ.
2. Le champ électrique radial à la surface de la sphère est renforcé par le facteur diélectrique :
3. La charge \(q(t)\) accumulée crée un champ répulsif \(E_{\text{q}}\).
4. La saturation est atteinte quand le champ répulsif annule le champ externe attractif.
Cette formule lie la charge maximale aux propriétés diélectriques de la particule et à la géométrie.
Où \(\varepsilon_0\) est la permittivité du vide, \(\varepsilon_{\text{r}}\) la permittivité relative de la particule.
Étape 1 : Préparation des Données
| Paramètre | Valeur |
|---|---|
| Champ électrique (\(E_{\text{moy}}\)) | \(3.33 \times 10^5\) V/m |
| Diamètre particule (\(d_{\text{p}}\)) | \(2.5 \times 10^{-6}\) m |
| Permittivité vide (\(\varepsilon_0\)) | \(8.854 \times 10^{-12}\) F/m |
| Permittivité relative (\(\varepsilon_{\text{r}}\)) | 5.0 |
Le terme entre parenthèses \(\frac{3\varepsilon_{\text{r}}}{\varepsilon_{\text{r}} + 2}\) varie entre 1 (pour \(\varepsilon_{\text{r}}=1\), isolant parfait comme l'air) et 3 (pour \(\varepsilon_{\text{r}} \to \infty\), conducteur parfait). Avec une valeur de 5.0, nous sommes proches d'un comportement conducteur, ce qui favorise la prise de charge.
Calculs Détaillés
Nous procédons d'abord au calcul du facteur diélectrique (facteur de Clausius-Mossotti modifié), puis au calcul global de la charge.
1. Calcul du Facteur Diélectrique :Ce terme représente la capacité de la particule à concentrer les lignes de champ. Pour une particule conductrice infinie, ce terme vaudrait 3.
On injecte maintenant toutes les valeurs dans l'équation de Pauthenier.
Interprétation Globale : La charge acquise est extrêmement faible en valeur absolue (\(10^{-16}\) C). C'est normal pour une micro-particule isolée. Cependant, rapportée à sa masse infime, le rapport charge/masse est très élevé, ce qui permettra une accélération significative.
La charge élémentaire d'un électron est \(e \approx 1.6 \times 10^{-19}\) C. Notre particule a donc acquis environ \(1.24 \times 10^{-16} / 1.6 \times 10^{-19} \approx 775\) électrons excédentaires. C'est un ordre de grandeur parfaitement réaliste pour une particule de 2.5 microns dans un champ intense.
Si la résistivité des particules est trop élevée (> \(10^{10} \Omega.cm\)), elles peuvent conserver leur charge une fois arrivées sur la plaque collectrice, créant une couche isolante. Cela peut mener au phénomène de "Back Corona" (contre-émission) qui ruine l'efficacité de l'électrofiltre. Ici, nous supposons que la décharge sur la plaque se fait normalement.
🎯 Objectif Scientifique
Nous cherchons la vitesse transversale (\(w\)) à laquelle la particule se déplace vers la plaque collectrice. Cette vitesse est le résultat d'un équilibre dynamique entre la force électrique qui "tire" la particule et la force de frottement de l'air (traînée) qui la "freine".
📚 Référentiel
Loi de Stokes Force de CoulombLa particule atteint sa vitesse limite quasi-instantanément (temps de relaxation de l'ordre de la milliseconde). Nous pouvons donc écrire l'égalité des forces : \(F_{\text{électrique}} = F_{\text{traînée}}\). Pour des particules fines dans un gaz laminaire, la loi de Stokes s'applique.
Dans un fluide visqueux, une sphère en mouvement subit une force de traînée \(F_{\text{d}} = 3 \pi \eta d_{\text{p}} w\) (Loi de Stokes pour \(Re < 1\)). La force motrice électrique est \(F_{\text{e}} = q_{\text{s}} E\). À l'équilibre, \(F_{\text{e}} = F_{\text{d}}\). Si la particule est très petite (< 1 µm), le milieu n'est plus continu et il faut introduire le facteur de correction de Cunningham \(C_{\text{c}}\). Ici, avec 2.5 µm, nous sommes à la limite mais nous négligeons \(C_{\text{c}}\) (soit \(C_{\text{c}} \approx 1\)).
⚖️ Visualisation : Diagramme de Corps Libre (Équilibre Dynamique)
Manipulation algébrique du bilan des forces : On part du Principe Fondamental de la Dynamique (PFD) :
La phase d'accélération est négligeable (\(t \to 0\)), donc \(\sum \vec{F} = \vec{0}\). Forces en présence : 1. Force Électrique :
2. Force de Traînée (Fluide) :
À l'équilibre scalaire : \(q_{\text{s}} E = 3 \pi \eta d_{\text{p}} w\). On isole \(w\) en divisant par le terme de friction :
L'égalité des forces nous permet d'isoler la vitesse \(w\).
Avec \(\eta\) la viscosité dynamique du gaz.
Étape 1 : Hypothèses & Données
| Paramètre | Valeur |
|---|---|
| Charge (\(q_{\text{s}}\)) | \(1.24 \times 10^{-16}\) C |
| Viscosité (\(\eta\)) | \(2.4 \times 10^{-5}\) Pa.s |
| Diamètre (\(d_{\text{p}}\)) | \(2.5 \times 10^{-6}\) m |
Notez que la vitesse de migration dépend du carré du champ électrique (car \(q_{\text{s}}\) dépend déjà de \(E\)). Doubler la tension quadruple théoriquement la vitesse de capture (tant qu'on ne claque pas). C'est pourquoi le maintien d'une haute tension stable est critique.
Calculs Détaillés
Calculons d'abord le numérateur (Force motrice) puis le dénominateur (Facteur de friction) pour obtenir la vitesse.
1. Force de Coulomb (Motrice) :Interprétation Globale : La particule se déplace transversalement à 7.3 cm/s. Cette vitesse peut sembler faible comparée à la vitesse du gaz (souvent 1 à 2 m/s longitudinalement), mais c'est suffisant si le temps de séjour dans l'appareil est assez long (plusieurs secondes).
Pour des cendres volantes dans des conditions industrielles, les vitesses de migration varient typiquement entre 5 cm/s et 15 cm/s. Notre résultat de 7.3 cm/s est donc parfaitement cohérent et valide le modèle.
Cette vitesse est théorique. Dans la réalité, la turbulence du gaz, le "ré-entraînement" des poussières lors du décolmatage (frappage des plaques) et l'hétérogénéité du champ réduisent la vitesse effective. On applique souvent un facteur empirique de correction (ex: \(w_{\text{k}} = 0.5 \times w_{\text{th}}\)) dans les projets réels pour être conservateur.
🎯 Objectif Scientifique
C'est l'étape finale du design. Connaissant la vitesse de migration et le débit à traiter, nous devons déterminer la Surface de Collecte Totale (\(A\)) nécessaire pour garantir que 99.5% des particules seront capturées avant la sortie de l'appareil. C'est le dimensionnement "process" qui fixera le coût de l'installation.
📚 Référentiel
Loi de Deutsch-AndersonLe modèle de Deutsch-Anderson considère que la turbulence mélange parfaitement les particules dans la section droite du canal. La concentration diminue donc exponentiellement le long du filtre. Plus l'appareil est grand (\(A\)) ou la vitesse de migration élevée (\(w\)), meilleure est l'efficacité. Inversement, un débit élevé (\(Q\)) dilue l'effort de capture.
L'équation de Deutsch-Anderson est la loi fondamentale du dimensionnement des ESP. Elle exprime l'efficacité \(\eta\) comme suit : \(\eta = 1 - \exp(-w \cdot SCA)\), où \(SCA = A/Q\) est la Surface Spécifique de Collecte (Specific Collection Area). C'est une loi statistique basée sur la probabilité de capture.
📈 Visualisation : Loi Exponentielle d'Efficacité
Obtention de la surface A : 1. On part de l'équation d'efficacité :
2. On cherche à isoler \(x\) (qui contient A). On passe le 1 de l'autre côté :
3. On applique le logarithme népérien (fonction réciproque de l'exponentielle) :
4. On remplace \(x\) par \(\frac{A \cdot w}{Q}\).
5. On obtient :
6. On isole A en multipliant par \(-Q/w\) :
L'efficacité \(\eta_{\text{eff}}\) est liée à la surface \(A\) par une exponentielle décroissante.
En inversant la formule, on isole la surface \(A\).
Étape 1 : Données Techniques
| Type | Valeur |
|---|---|
| Débit volumique (\(Q\)) | 50 m³/s |
| Efficacité cible (\(\eta_{\text{eff}}\)) | 0.995 (99.5%) |
| Vitesse de migration (\(w\)) | 0.073 m/s |
L'efficacité est très sensible à l'exposant. Passer de 99% à 99.9% (diviser les rejets par 10) ne demande pas "un peu plus" de surface, mais beaucoup plus (logarithmiquement). Chaque "9" supplémentaire coûte cher en acier !
Calcul de la Surface
On retourne la formule pour exprimer \(A\) en fonction des autres paramètres : \(A = - \frac{Q}{w} \ln(1 - \eta_{\text{eff}})\).
1. Calcul du terme logarithmique (Difficulté de capture) :Interprétation Globale : Il faut déployer environ 3630 m² de tôles collectrices pour traiter ce débit avec cette efficacité. C'est une surface considérable, impliquant un appareil de grande taille, probablement constitué de plusieurs "champs" (sections électriques indépendantes) en série.
Pour un débit de 50 m³/s, une surface de 3630 m² donne une "Surface Spécifique de Collecte" (SCA = A/Q) de \(3630/50 = 72.6\) s/m. C'est une valeur typique pour des électrofiltres haute performance dans l'industrie cimentière (généralement entre 50 et 100 s/m).
L'équation de Deutsch suppose un profil de vitesse plat. Si la distribution du gaz est mauvaise (zones mortes, passages préférentiels), l'efficacité réelle chulera drastiquement. Il est impératif de prévoir des grilles de distribution en entrée d'électrofiltre et de réaliser une simulation CFD (Mécanique des Fluides Numérique) pour valider l'écoulement.
📄 Livrable Final (Rapport APD)
SOLUTIONS
| Ind. | Date | Objet de la modification | Rédacteur |
|---|---|---|---|
| A | 10/10/2024 | Première émission - Hypothèses préliminaires | Ing. J. Dupont |
| B | 24/10/2024 | Validation Client & Dimensionnement Final | Ing. M. Martin |
- Mécanisme de charge : Saturation par champ (Pauthenier) prépondérante (\(d_{\text{p}} > 1\mu m\)).
- Régime hydraulique : Écoulement Piston avec mélange turbulent transverse (Deutsch-Anderson).
- Facteur de correction de vitesse : 1.0 (Théorique pur).
| Débit Gaz (Q) | 50 m³/s |
| Température | 150 °C |
| Tension Service | 50 kV |
Synthèse des grandeurs calculées pour une particule critique de 2.5 µm.
Bureau d'Étude PureAir
Expert Sénior
Ajouter un commentaire