La variation locale de la pression transmembranaire modifie à la fois le flux de perméat et le champ d’écoulement axial. Dans une unité pilote de filtration membranaire en ligne, la vitesse locale de perméation est approximativement proportionnelle à la pression transmembranaire locale (TMP) par l’intermédiaire de la perméabilité hydraulique de la membrane. Lorsque le perméat quitte le canal d’alimentation, le débit et la vitesse d’écoulement axiaux varient le long de la membrane, ce qui modifie à son tour la perte de charge par frottement ainsi que la TMP disponible plus en aval. La perte de charge et la collecte du perméat doivent donc être modélisées comme un problème couplé et spatialement variable, plutôt que comme des calculs indépendants entre l’entrée et la sortie.

À retenir : Une TMP locale élevée augmente le flux de perméat, mais le retrait de fluide qui en résulte modifie le profil de vitesse axiale et la perte de charge. Une modélisation précise d’une unité pilote nécessite de résoudre conjointement le transport membranaire et la perte de charge dans le canal le long de la direction d’écoulement.

Pourquoi la TMP locale est importante

La TMP détermine le flux local de perméat

Le flux transmembranaire local peut être représenté par une relation fondée sur les résistances :

[ J(z) = \frac{\Delta P_{\mathrm{TMP}}(z)} {\mu \left[R_m + R_g(z) + R_f(z)\right]} ]

où (J) est le flux local de perméat, (\Delta P_{\mathrm{TMP}}) est la TMP locale, (\mu) est la viscosité, et (R_m), (R_g) et (R_f) représentent les résistances de la membrane, de la couche de gel et de l’encrassement.

Si les termes de résistance sont définis de manière à inclure la viscosité, le facteur de viscosité est omis. La relation essentielle est que le flux dépend de la différence de pression locale et de la résistance locale, et non simplement de la pression moyenne à travers l’unité.

La TMP varie le long de la membrane

Dans un canal membranaire, la TMP locale est déterminée par la différence de pression entre le fluide côté alimentation et le fluide côté perméat :

[ \Delta P_{\mathrm{TMP}}(z)

P_f(z) - P_p(z) ]

Les deux pressions peuvent varier selon la position axiale. La pression côté alimentation diminue généralement en raison de la perte de charge par frottement, tandis que la pression côté perméat peut également varier en raison de la résistance du canal de perméat, de la géométrie de collecte ou de la contre-pression.

Par conséquent, une même membrane peut subir des forces motrices sensiblement différentes à son entrée et à sa sortie.

La variation de pression crée une sollicitation non uniforme de la membrane

Lorsque la TMP est plus élevée, la membrane tend à produire davantage de perméat, à condition que l’encrassement et la polarisation de concentration ne compensent pas l’augmentation de pression. Lorsque la TMP est plus faible, le flux local diminue et moins de perméat est collecté.

Cela produit une distribution non uniforme du flux. Un modèle fondé uniquement sur la TMP moyenne peut donc représenter de manière erronée l’encrassement local, la polarisation de concentration, la récupération de perméat et la surface membranaire requise.

Comment la perméation modifie la perte de charge

L’écoulement axial diminue lorsque le perméat est évacué

L’extraction du perméat retire du fluide du canal d’alimentation. Le débit volumique axial local diminue donc avec la distance le long de la membrane :

[ \frac{dQ}{dz} = -wJ(z) ]

où (Q(z)) est le débit axial d’alimentation, (w) est la largeur de membrane exposée à la filtration, et (J(z)) est le flux local de perméat.

Pour un canal de section transversale approximativement constante, la vitesse axiale moyenne diminue également :

[ u(z) = \frac{Q(z)}{A} ]

Il s’agit de l’analogue en phase liquide d’un écoulement à vitesse variable dans d’autres opérations unitaires pilotes : le retrait local modifie les conditions d’écoulement utilisées pour calculer la perte de charge en aval.

Une vitesse axiale plus faible réduit généralement les pertes par frottement

Dans l’hypothèse d’un écoulement de lubrification, le gradient de pression est lié au débit axial local. Pour un canal simple de type fente, la relation générale s’écrit

[ -\frac{dp}{dz} \propto \frac{\mu Q(z)}{h^3 w} ]

où (h) est l’espace du canal et (w) est la largeur du canal.

Lorsque (Q(z)) diminue en raison de la perméation, le gradient de pression axial local diminue également. Par conséquent, un canal perméable peut présenter une perte de charge axiale totale inférieure à celle d’un canal imperméable par ailleurs identique, transportant le débit d’entrée initial sur toute sa longueur.

La perte de charge est couplée au profil de flux

Le couplage fonctionne dans les deux sens :

La TMP locale détermine (J(z)).
(J(z)) détermine la vitesse à laquelle le débit axial diminue.
Le débit axial détermine le gradient local de pression par frottement.
Le profil de pression détermine la TMP plus en aval.

Cette rétroaction signifie que le profil de pression ne peut généralement pas être défini indépendamment du profil de flux membranaire.

Comment modéliser une unité pilote en ligne

Définir les variables d’état locales

Un modèle unidimensionnel pratique suit généralement :

La pression côté alimentation, (P_f(z))
La pression côté perméat, (P_p(z))
Le débit axial d’alimentation, (Q(z))
La TMP locale, (\Delta P_{\mathrm{TMP}}(z))
Le flux local de perméat, (J(z))
Les résistances de la membrane, du gel et de l’encrassement

Le modèle doit également inclure la viscosité, les dimensions du canal, la perméabilité de la membrane, les conditions d’entrée, les conditions de sortie et toute résistance hydraulique côté perméat.

Avancer le long de la membrane

Le canal peut être divisé en petits segments axiaux (\Delta z). À chaque segment, le modèle doit :

Calculer la TMP locale à partir des pressions d’alimentation et de perméat.
Calculer la résistance locale et le flux de perméat.
Retirer le débit de perméat correspondant du flux d’alimentation axial.
Recalculer la vitesse axiale.
Mettre à jour la pression côté alimentation à l’aide du gradient hydraulique local.
Mettre à jour la pression côté perméat si les pertes dans le canal de perméat sont importantes.

Une procédure numérique d’avancement, telle que la méthode d’Euler explicite, est souvent suffisante pour un modèle initial d’unité pilote :

[ Q_{i+1}

Q_i - wJ_i\Delta z ]

[ P_{f,i+1}

P_{f,i} + \left(\frac{dp_f}{dz}\right)_i\Delta z ]

La taille des segments doit être suffisamment petite pour résoudre les variations significatives de pression, de flux et de résistance.

Utiliser les conditions aux limites de manière cohérente

Le modèle peut être défini à l’aide de la pression d’entrée de l’alimentation et de la pression de sortie du perméat, ou à l’aide de conditions de débit et de pression appliquées à différentes limites. Ces choix déterminent si le problème constitue un simple calcul de problème à valeur initiale ou s’il nécessite une résolution itérative.

Par exemple, si les pressions d’alimentation et de rétentat sont toutes deux fixées tandis que la pression côté perméat est également imposée, le modèle peut devoir ajuster un débit d’entrée inconnu ou une distribution de flux afin de satisfaire simultanément toutes les conditions.

Vérifier la conservation

Le débit de perméat intégré doit être égal à la différence entre le débit d’alimentation à l’entrée et le débit de rétentat :

[ Q_{\mathrm{perm,total}}

Q_{\mathrm{feed,in}}

Q_{\mathrm{retentate,out}} ]

Cette vérification du bilan massique est essentielle. Un profil de pression peut sembler physiquement cohérent tout en produisant un débit total de perméat incorrect si l’intégration du flux local n’est pas cohérente.

Ce que le modèle révèle sur les performances de l’unité pilote

La collecte du perméat peut être non uniforme

Un maximum local de TMP ne garantit pas automatiquement la production utile de perméat la plus élevée. Un flux plus important peut accroître la polarisation de concentration ou l’encrassement, en augmentant (R_g) ou (R_f) et en réduisant la réponse effective du flux.

Le débit de perméat mesuré doit donc être interprété conjointement avec la pression axiale, le débit d’alimentation et les données de résistance.

Le débit de rétentat contrôle les conditions en aval

Au fur et à mesure de la filtration, le débit de rétentat diminue. Cela modifie le nombre de Reynolds du canal, les conditions de cisaillement et la perte de charge par frottement le long du module.

Dans un système à écoulement tangentiel, la diminution du débit axial peut également affecter l’action de balayage qui limite l’accumulation de la couche de gel et de l’encrassement. Le bénéfice hydraulique d’une perte de charge plus faible peut donc s’accompagner d’une modification des conditions de transfert de matière.

La géométrie des fibres et des canaux est importante

Dans les systèmes à fibres creuses, la perméation à travers la paroi réduit la vitesse axiale moyenne. Selon les hypothèses simplifiées de propriétés constantes et de perméation uniforme à travers la paroi, la vitesse peut être représentée par

[ u(z) = u_0 - \frac{2v_w z}{R} ]

où (u_0) est la vitesse d’entrée, (v_w) est la vitesse normale à la paroi et (R) est le rayon de la fibre.

Ce résultat illustre le principe général suivant : le flux à travers la paroi modifie la vitesse axiale, et la vitesse axiale modifie la perte de charge. L’équation exacte dépend de la géométrie et des hypothèses ; un modèle de canal plan ne doit donc pas être directement transposé à un module à fibres creuses sans modification.

Comprendre les compromis

Une TMP plus élevée n’améliore pas toujours la productivité

L’augmentation de la TMP accroît la force motrice immédiate de la perméation, mais elle peut également accélérer le blocage des pores, la formation du gâteau, la croissance de la couche de gel ou l’encrassement irréversible.

La plage de fonctionnement utile est donc déterminée par le flux durable et la qualité du produit, et non par le flux initial le plus élevé.

La TMP moyenne peut masquer des problèmes locaux

L’utilisation d’une seule TMP moyenne simplifie les calculs, mais peut dissimuler une région de sortie à faible flux ou une région d’entrée surchargée. Cela peut conduire à des estimations incorrectes de la surface membranaire, de la vitesse d’encrassement et de la récupération de perméat.

Les modèles fondés sur la pression moyenne sont les plus défendables lorsque la variation axiale de pression est manifestement faible par rapport à la TMP de fonctionnement.

Une longueur excessive du module peut réduire la validité du modèle

Si la perméation retire une quantité trop importante de débit axial, la vitesse peut approcher zéro près de l’extrémité aval. Pour la mise à l’échelle des fibres creuses décrite dans les références, le maintien d’un régime d’écoulement axial significatif nécessite une longueur nettement inférieure à

[ L \ll \frac{u_0R}{2v_w} ]

Lorsque cette condition n’est pas satisfaite, des hypothèses telles qu’une échelle de vitesse presque uniforme, un écoulement simple de lubrification ou une faible variation axiale peuvent ne plus être valables.

La TFF et la NFF se comportent différemment

Dans la filtration à écoulement normal, la matière est dirigée directement vers la membrane, ce qui peut favoriser une formation rapide du gâteau et une diminution du flux. La filtration à écoulement tangentiel utilise un écoulement transversal pour balayer la surface de la membrane et limite généralement l’accumulation de la couche de gel et de l’encrassement.

Un modèle qui inclut la TMP mais ignore la résistance dépendante de l’écoulement transversal peut donc surestimer le flux à long terme, particulièrement lors d’un fonctionnement en TFF.

Faire le bon choix selon votre objectif

Utilisez un modèle couplé à résolution spatiale lorsque la variation de pression, la résistance côté perméat, l’encrassement ou l’échelle du module est importante.

Si votre objectif principal est de prévoir la perte de charge : Couplez le bilan du débit axial à une relation d’écoulement de lubrification ou de frottement dans le canal, en mettant à jour la vitesse locale après chaque étape de perméation.
Si votre objectif principal est la récupération de perméat : Intégrez le flux local à partir des pressions locales d’alimentation et de perméat plutôt que de multiplier la surface membranaire par un flux moyen.
Si votre objectif principal est le comportement d’encrassement : Incluez la résistance de la couche de gel et de l’encrassement en fonction du flux local et des conditions d’écoulement transversal.
Si votre objectif principal est la mise à l’échelle de l’unité pilote : Mesurez la pression axiale, le débit d’alimentation, le débit de perméat et le flux de filtrat afin de comparer le modèle spatial aux données d’exploitation réparties.

Un modèle fiable de membrane en ligne considère la pression, le débit axial, la perméation et la résistance comme des variables interdépendantes plutôt que comme des calculs séparés.

Tableau récapitulatif :

Facteur clé	Effet sur le système	Implication pour la modélisation
TMP locale	Détermine le flux de perméat à chaque position de la membrane	Calculer le flux à partir des pressions locales d’alimentation et de perméat
Retrait du perméat	Réduit le débit et la vitesse axiaux d’alimentation	Mettre à jour le débit et la vitesse le long de la membrane
Perte de charge axiale	Réduit la pression côté alimentation en aval	Coupler le frottement dans le canal à la perméation
Résistance côté perméat	Modifie la TMP locale par la variation de pression	Inclure les pertes dans le canal de perméat lorsqu’elles sont importantes
Résistance due à l’encrassement et au gel	Réduit le flux effectif et augmente la résistance	Modéliser la résistance en fonction des conditions locales
Géométrie du module	Contrôle le comportement hydraulique et les profils de vitesse	Utiliser des équations spécifiques à la géométrie pour les canaux plans ou les fibres creuses
Conditions aux limites	Déterminent si le débit ou la pression d’entrée doit être résolu de manière itérative	Appliquer de manière cohérente les conditions d’alimentation, de rétentat et de perméat
Vérification de la conservation	Confirme que le perméat total est égal à la perte de débit d’alimentation	Vérifier le débit de perméat intégré par rapport au bilan massique

Vous planifiez une unité pilote de filtration membranaire en ligne ou souhaitez améliorer un procédé existant ? Contactez LABPARK pour découvrir nos unités pilotes d’opérations unitaires éducatives et professionnelles, conçues pour les universités, les instituts de recherche et les entreprises. Nos systèmes de génie chimique, de bioprocédés et biotechnologies, ainsi que de traitement environnemental et de l’eau, vous aident à évaluer la TMP locale, la perte de charge, la récupération de perméat, le comportement d’encrassement et les performances de mise à l’échelle grâce à des données pratiques d’unité pilote. LABPARK peut vous aider à sélectionner la configuration de filtration membranaire adaptée, à intégrer des points de mesure pour les pressions d’alimentation et de perméat, et à établir une base expérimentale fiable pour le développement des procédés, l’enseignement, la recherche et la validation industrielle.