Le bilan massique dynamique est linéarisé en remplaçant la relation non linéaire d'équilibre vapeur-liquide par sa forme linéaire locale, puis en regroupant les termes de composition. Pour l'étage (n), le bilan du composant dans la phase liquide est

[ H_n\frac{dx_n}{dt}

Lx_{n-1}-Lx_n+Vy_{n+1}-Vy_n. ]

Si la relation d'équilibre est approximée localement par (y_i=m_i x_i), avec une pente constante (m_i), l'équation de l'étage devient une équation différentielle linéaire couplée aux étages voisins. Cette forme sert de base à la dérivation des fonctions de transfert, à l'étude de la propagation des perturbations et à la conception de systèmes de contrôle de la composition pour les installations pilotes d'absorption et de distillation.

À retenir : La linéarisation transforme le bilan de l'étage en un modèle dynamique du premier ordre en considérant la rétention, les débits et la pente locale d'équilibre comme constants. Le modèle obtenu est généralement un ensemble couplé d'équations du premier ordre, et non une équation indépendante pour chaque étage, car chaque étage échange de la matière avec les étages adjacents.

Commencer par le bilan dynamique de l'étage

Accumulation du composant dans la phase liquide

Pour une rétention liquide constante (H_n), l'accumulation du composant absorbé ou distillé sur l'étage (n) est

[ H_n\frac{dx_n}{dt}. ]

Ici, (x_n) est la fraction molaire du composant dans la phase liquide sur l'étage (n). Les unités doivent être cohérentes : si (H_n) est une rétention molaire, (L) et (V) doivent être des débits molaires.

Contributions des débits

Le flux liquide entrant dans l'étage (n) apporte (Lx_{n-1}), tandis que le flux liquide sortant emporte (Lx_n).

Le flux de vapeur entrant depuis l'étage adjacent apporte (Vy_{n+1}), et la vapeur quittant l'étage (n) évacue (Vy_n). Par conséquent,

[ H_n\frac{dx_n}{dt}

Lx_{n-1}-Lx_n+Vy_{n+1}-Vy_n. ]

Ce bilan suppose que (L) est constant, que (V) est constant, que la rétention est constante, qu'il n'y a pas de réaction et que l'accumulation dans la phase vapeur est négligeable.

Application de l'approximation de l'équilibre local

Relation d'équilibre linéaire

Les compositions de la vapeur et du liquide sont reliées par la fonction d'équilibre

[ y=f(x). ]

Autour d'un point de fonctionnement sélectionné, cette relation peut être approximée comme suit :

[ y_i \approx y_{i,0}+m_i(x_i-x_{i,0}), ]

où

[ m_i=\left.\frac{dy}{dx}\right|{x{i,0}}. ]

Pour un modèle simplifié avec une pente constante de la droite d'équilibre et des compositions mesurées par rapport au point de fonctionnement, cela devient

[ y_i=m x_i. ]

La pente (m) peut représenter une relation locale de volatilité ou une relation de type loi de Henry, selon le système.

Substitution dans le bilan de l'étage

En utilisant (y_n=mx_n) et (y_{n+1}=mx_{n+1}),

[ H_n\frac{dx_n}{dt}

Lx_{n-1}-Lx_n+mVx_{n+1}-mVx_n. ]

En regroupant les termes faisant intervenir (x_n), on obtient

[ \boxed{ H_n\frac{dx_n}{dt}

Lx_{n-1}

(L+mV)x_n + mVx_{n+1} } ]

ou, après division par la rétention,

[ \boxed{ \frac{dx_n}{dt}

\frac{L}{H_n}x_{n-1}

\frac{L+mV}{H_n}x_n + \frac{mV}{H_n}x_{n+1} }. ]

Il s'agit de l'équation dynamique linéarisée standard pour un étage d'équilibre idéal dans l'hypothèse de débits constants.

Signification de la « linéarisation » dans ce modèle

Linéarisation de la relation d'équilibre

La relation d'équilibre initiale peut être non linéaire, notamment dans les systèmes de distillation présentant une variation importante de la volatilité relative. La linéarisation remplace cette relation non linéaire par sa tangente autour du point de fonctionnement stationnaire sélectionné.

L'approximation est fiable pour des variations de composition relativement faibles autour de ce point. Des perturbations importantes peuvent nécessiter la relation complète d'équilibre vapeur-liquide non linéaire.

Linéarisation du bilan matière

Le bilan matière est déjà linéaire en (x) lorsque (H), (L) et (V) sont constants. Des non-linéarités supplémentaires apparaissent lorsque ces grandeurs varient avec la composition, la température, la pression, le niveau ou le débit total.

Pour un modèle plus général, on écrit

[ H_n(x_n)\frac{dx_n}{dt}

L(x_{n-1}-x_n) + V(y_{n+1}-y_n), ]

puis on développe le modèle autour d'un état stationnaire :

[ x_i=x_{i,0}+\delta x_i, \qquad y_i=y_{i,0}+\delta y_i. ]

La perturbation de l'équilibre est approximée par

[ \delta y_i=m_i\delta x_i. ]

Le modèle de perturbation obtenu possède la même structure que l'équation encadrée, mais ses coefficients sont évalués au point de fonctionnement.

Correction importante de l'indexation

La composition de la vapeur voisine doit être conservée

Une erreur fréquente d'indexation consiste à substituer (y_{n+1}=mx_{n-1}). Dans la numérotation habituelle des étages à contre-courant, la vapeur entrant dans l'étage (n) est (y_{n+1}), qui est liée à la composition liquide de l'étage (n+1) :

[ y_{n+1}=mx_{n+1}. ]

Par conséquent, l'équation correctement indexée est

[ H_n\frac{dx_n}{dt}

Lx_{n-1}

(L+mV)x_n + mVx_{n+1}. ]

L'expression

[ H\frac{dx_n}{dt}

Lx_{n-1}-(L+mV)x_n+mVx_{n-1} ]

n'est valide que si la variable notée (x_{n-1}) a été délibérément définie comme la composition équivalente liquide amont de la vapeur entrante. Il ne s'agit pas de la substitution d'équilibre habituelle indexée par étage.

Pourquoi cette distinction est importante

La conservation de (x_{n+1}) produit un modèle de colonne couplé. La composition d'un étage réagit aux deux étages voisins, ce qui est essentiel pour représenter la propagation des perturbations dans une colonne à plateaux ou une colonne à garnissage discrétisée.

Le remplacement du terme correspondant à l'étage voisin par la même variable amont peut produire un modèle plus simple à un seul étage, mais il modifie la connectivité physique et peut conduire à des prédictions dynamiques incorrectes.

Des équations d'étage aux fonctions de transfert

Forme dans le domaine de Laplace

En supposant des perturbations initiales nulles et des coefficients constants, la transformée de Laplace donne

[ H_n sX_n(s)

LX_{n-1}(s)

(L+mV)X_n(s) + mVX_{n+1}(s). ]

Après réarrangement,

[ \boxed{ \left[H_n s+(L+mV)\right]X_n(s)

LX_{n-1}(s)+mVX_{n+1}(s) }. ]

Cette équation constitue une ligne du modèle dynamique de fonction de transfert ou d'espace d'état de la colonne.

Approximation à un seul étage

Si la composition de la vapeur entrante est traitée comme une entrée connue plutôt que comme un état lié à (x_{n+1}), on définit

[ u_n=Lx_{n-1}+Vy_{n+1}. ]

Le bilan devient

[ H_n\frac{dx_n}{dt}

u_n-(L+mV)x_n. ]

La fonction de transfert correspondante est

[ \boxed{ \frac{X_n(s)}{U_n(s)}

\frac{1}{H_n s+L+mV} } ]

avec une constante de temps

[ \boxed{ \tau_n=\frac{H_n}{L+mV} }. ]

Si l'entrée est uniquement la composition liquide amont et que la contribution de la vapeur est négligée ou intégrée séparément, alors

[ \frac{X_n(s)}{X_{n-1}(s)}

\frac{L}{H_n s+L+mV}. ]

Application à l'absorption et à la distillation

Absorption gaz-liquide

Dans l'absorption, le soluté passe généralement de la phase gazeuse à la phase liquide lorsque la composition réelle du gaz dépasse sa valeur d'équilibre :

[ y>y^*(x). ]

La pente d'équilibre (m) décrit la variation de la composition d'équilibre du gaz avec la composition du liquide. Un terme (mV) plus important augmente l'intensité du couplage avec la phase vapeur et modifie la réponse de l'étage.

Pour les colonnes d'absorption à garnissage, le même concept est souvent formulé comme un modèle différentiel distribué plutôt que comme des bilans discrets sur des plateaux. Une discrétisation en volumes finis produit des équations de type étage avec rétention et couplage entre cellules voisines.

Distillation

En distillation, la relation d'équilibre est généralement représentée par

[ y_i=y^*(x_i), ]

où (y^*(x)) peut être non linéaire. Autour d'une composition de fonctionnement,

[ \delta y_i \approx \left.\frac{dy^*}{dx}\right|{x{i,0}}\delta x_i. ]

La pente locale peut varier d'un étage à l'autre. Une colonne de distillation possède donc généralement des coefficients (m_n), (L_n), (V_n) et (H_n), plutôt qu'un unique (m), (L) et (V).

Conditions aux limites

Les équations des étages internes ne suffisent pas à elles seules. Le condenseur, le rebouilleur, l'étage d'alimentation, le reflux, le distillat, le produit de fond et les extrémités de la colonne fournissent les conditions aux limites et les relations avec les variables manipulées.

Par exemple, le reflux modifie le débit liquide entrant dans les étages supérieurs, tandis que le débit de rebouillage modifie le débit de vapeur entrant dans les étages inférieurs. Ces effets aux limites constituent souvent les principales voies de propagation des perturbations de commande.

Extension du modèle à une rétention et des débits variables

Bilan massique total

Lorsque le niveau de liquide ou la rétention varie, un bilan total doit être inclus. Le bilan général de rétention liquide est

[ \frac{dH_n}{dt}

L_{n-1}-L_n+V_{n+1}-V_n, ]

avec les signes et les définitions des flux adaptés à la convention d'étage choisie.

Le bilan du composant doit alors être écrit en termes de rétention du composant :

[ \frac{d(H_nx_n)}{dt}

\text{composant entrant}

\text{composant sortant}. ]

Le développement de la dérivée donne

[ H_n\frac{dx_n}{dt} + x_n\frac{dH_n}{dt}

\text{composant entrant}

\text{composant sortant}. ]

Le terme supplémentaire (x_n,dH_n/dt) est important lorsque les niveaux de liquide sur les plateaux ou les rétentions du rebouilleur et du condenseur varient de manière significative.

Effets de l'énergie et de la pression

Pour des modèles précis d'installations pilotes, la température affecte à la fois l'équilibre et les débits. La pression influence la composition de la vapeur, les conditions de saturation et le débit volumique.

Une linéarisation limitée à la composition est donc principalement adaptée aux perturbations modestes autour d'un point de fonctionnement à température et pression constantes. Des variations plus importantes nécessitent des équations couplées de bilan de matière, d'énergie et de pression.

Comprendre les compromis

Avantages du modèle linéarisé

Le modèle linéarisé est compact, efficace sur le plan du calcul et adapté à l'analyse par fonctions de transfert. Il donne un aperçu direct des constantes de temps, du couplage entre étages, de l'atténuation des perturbations et de la contrôlabilité.

Il est particulièrement utile pour le réglage des contrôleurs et les expériences pédagogiques sur installations pilotes lorsque la colonne fonctionne près d'une condition nominale stable.

Limites d'une pente d'équilibre constante

Une seule valeur de (m) peut être inadéquate lorsque la volatilité relative varie fortement avec la composition ou la température. L'approximation peut également échouer à proximité des azéotropes, dans les mélanges fortement non idéaux ou sur de larges plages de fonctionnement.

Pour ces systèmes, utilisez des pentes locales propres à chaque étage ou conservez la fonction d'équilibre non linéaire.

Limites des débits et de la rétention constants

Supposer que (L), (V) et (H) sont constants exclut les dynamiques de niveau et les effets hydrauliques. Cela peut conduire à des prévisions trop optimistes de la vitesse de réponse et masquer d'importantes interactions entre le débit, la pression et le contrôle de la composition.

Les installations pilotes devraient utiliser des rétentions mesurées ou modélisées indépendamment lorsque le noyage des plateaux, le pleurage, la formation de mousse ou une variation importante du niveau de liquide sont possibles.

Confusion entre les droites de fonctionnement et les bilans dynamiques

La droite de fonctionnement décrit une relation stationnaire de bilan matière entre les compositions des phases. Elle ne décrit pas à elle seule l'évolution des compositions au cours du temps.

Le modèle dynamique nécessite des termes de rétention tels que (H,dx/dt). La suppression de ces termes transforme le modèle en une relation stationnaire plutôt qu'en une relation transitoire.

Faire le bon choix selon votre objectif

Le modèle approprié dépend de la question de savoir si l'objectif est la conception du contrôle, l'analyse conceptuelle ou la prédiction haute fidélité du fonctionnement de l'installation pilote.

Si votre objectif principal est la conception d'un contrôleur : Linéarisez (y^*(x)) autour du point de fonctionnement prévu, conservez le couplage entre les étages voisins et assemblez les équations sous forme d'espace d'état.
Si votre objectif principal est l'étude d'une réponse simple à un seul étage : Traitez la composition de la vapeur entrante comme une entrée connue et utilisez la fonction de transfert du premier ordre (1/[Hs+(L+mV)]).
Si votre objectif principal est l'étude de la propagation des perturbations dans toute la colonne : Conservez (x_{n-1}) et (x_{n+1}) dans chaque équation d'étage et appliquez les conditions aux limites correctes du condenseur, du rebouilleur, de l'alimentation, du reflux et des produits.
Si votre objectif principal concerne des transitoires importants de l'installation pilote : Utilisez un équilibre non linéaire, une rétention variable, des débits variables et des bilans couplés d'énergie et de pression.

Un modèle linéarisé de colonne fiable préserve la connectivité physique des étages tout en simplifiant le comportement local de l'équilibre autour d'un point de fonctionnement clairement défini.

Tableau récapitulatif :

Élément clé	Interprétation du modèle linéarisé
Bilan dynamique	`H_n dx_n/dt = Lx_{n-1} - Lx_n + Vy_{n+1} - Vy_n`
Approximation de l'équilibre	Remplacer `y=f(x)` par `delta y=m delta x` près d'un point de fonctionnement stationnaire
Équation d'étage couplée	`H_n dx_n/dt = Lx_{n-1} - (L+mV)x_n + mVx_{n+1}`
Hypothèses principales	Rétention, débits liquide et vapeur constants, et accumulation dans la phase vapeur négligeable
Forme de la fonction de transfert	`[H_n s + (L+mV)]X_n = LX_{n-1} + mVX_{n+1}`
Constante de temps d'un seul étage	`tau_n = H_n/(L+mV)`
Point critique concernant l'indexation	Le terme de vapeur entrant est lié à `x_{n+1}`, ce qui préserve le couplage entre les étages voisins
Meilleures applications	Réglage des contrôleurs, analyse des perturbations, fonctions de transfert et modélisation en espace d'état
Situations nécessitant des modèles non linéaires	Perturbations importantes, rétention ou débit variables, forte non-idéalité, variations de pression ou effets énergétiques

Vous avez besoin d'une installation pilote pratique pour l'enseignement, la recherche ou le développement des procédés ? Contactez LABPARK pour discuter des installations pilotes d'opérations unitaires éducatives et professionnelles destinées à l'absorption gaz-liquide, à la distillation, aux bioprocédés et aux biotechnologies, ainsi qu'aux applications environnementales et de traitement de l'eau. LABPARK peut aider les universités, les instituts de recherche et les entreprises à sélectionner des équipements adaptés à la formation pratique, aux démonstrations de procédés, au réglage des contrôleurs, aux études de bilans massiques dynamiques et aux travaux expérimentaux évolutifs. Partagez votre procédé, votre capacité, vos conditions de fonctionnement ainsi que vos objectifs pédagogiques ou de recherche, et notre équipe vous recommandera une configuration d'installation pilote adaptée.