Utilisez les débits de composants comme inconnues. Dans les bilans matières des opérations unitaires et des unités pilotes en génie chimique, formulez chaque inconnue comme un débit massique, un débit molaire ou une quantité de composant individuel. Évitez d'utiliser un débit total de courant inconnu multiplié par une fraction massique ou molaire inconnue, car ce produit crée une équation non linéaire et empêche une résolution directe par les méthodes matricielles linéaires standard.

La formulation linéaire la plus fiable consiste à résoudre d'abord les débits de composants, puis à calculer les débits totaux et les compositions à partir de ces résultats. Cela produit un système d'équations de la forme Ax = b, qui peut être résolu efficacement avec l'algèbre linéaire numérique standard lorsque les équations sont indépendantes et que les degrés de liberté sont nuls.

Pourquoi le choix des inconnues détermine la linéarité

La représentation non linéaire des courants

Un débit de composant est souvent écrit comme suit :

[ \dot{m}_i = \dot{m}w_i ]

où (\dot{m}) est le débit massique total et (w_i) est la fraction massique du composant.

Si (\dot{m}) et (w_i) sont tous deux inconnus, leur produit contient deux variables inconnues. Le bilan résultant est non linéaire, même si la loi de conservation sous-jacente est linéaire.

La représentation linéaire par débit de composant

Définissez le débit du composant lui-même comme l'inconnue :

[ \dot{m}_i = x_i ]

Un bilan de composant peut alors être écrit directement comme suit :

[ \sum_j a_{ij}x_j = b_i ]

où les coefficients (a_{ij}) représentent des relations de courant connues, des facteurs de routage ou des coefficients stoechiométriques. Le vecteur inconnu contient les débits de composants plutôt que des produits de débits totaux et de compositions.

Récupération ultérieure des propriétés des courants

Une fois les débits de composants calculés, déterminez le débit total du courant à partir de :

[ \dot{m}_{\text{total}} = \sum_i \dot{m}_i ]

Calculez ensuite la fraction massique du composant en utilisant :

[ w_i = \frac{\dot{m}i}{\dot{m}{\text{total}}} ]

La même approche s'applique sur une base molaire :

[ \dot{n}_{\text{total}} = \sum_i \dot{n}_i \qquad y_i = \frac{\dot{n}i}{\dot{n}{\text{total}}} ]

Cela sépare le calcul du bilan linéaire du calcul de la composition.

Comment construire le modèle de bilan

Établir la base de calcul

Choisissez une base avant de définir les inconnues. Pour un fonctionnement continu en régime permanent, il s'agit généralement d'un débit par unité de temps ; pour un fonctionnement en discontinu (batch), il peut s'agir d'un lot ou d'une charge spécifiée.

Utilisez le débit massique ou molaire de manière cohérente dans tout le modèle, à moins que les conversions entre les bases ne soient explicitement incluses.

Dessiner le procédé et numéroter les courants

Un schéma de procédé doit identifier chaque courant entrant, sortant, recyclé, dérivé et intermédiaire. Chaque débit de composant associé à un courant inconnu devient une variable candidate.

Par exemple, un courant de sortie à deux composants peut être représenté par :

[ x_1 = \dot{m}{A,\text{out}} \qquad x_2 = \dot{m}{B,\text{out}} ]

plutôt que par la paire composée du débit de sortie total et d'une fraction de composition inconnue.

Écrire des bilans de composants indépendants

Pour une unité en régime permanent sans réaction :

[ \text{entrée de composant} = \text{sortie de composant} ]

Pour une unité réactive, incluez les termes connus de production ou de consommation stoechiométrique. Si les taux de réaction, les expressions de vitesse, les relations d'équilibre ou d'autres quantités sont également inconnus et multiplient d'autres inconnues, le modèle peut ne plus rester linéaire.

Assembler la forme matricielle

Après avoir arrangé chaque équation de manière à ce que les termes inconnus soient à gauche et les termes connus à droite, écrivez :

[ A\mathbf{x}=\mathbf{b} ]

Par exemple :

[ \begin{bmatrix} a_{11} & a_{12}\ a_{21} & a_{22} \end{bmatrix} \begin{bmatrix} x_1\ x_2 \end{bmatrix}

\begin{bmatrix} b_1\ b_2 \end{bmatrix} ]

Ici, (A) contient les coefficients connus, (\mathbf{x}) contient les inconnues de débit de composant, et (\mathbf{b}) contient les débits d'alimentation connus, les spécifications ou les termes sources.

Vérifier si le système est résoluble

Effectuer une analyse des degrés de liberté

Calculez :

[ \text{Degrés de liberté}

\text{nombre d'inconnues}

\text{nombre d'équations indépendantes} ]

Un système avec zéro degré de liberté est exactement spécifié et peut avoir une solution unique si la matrice des coefficients est de plein rang.

Une valeur positive indique que des mesures ou des spécifications supplémentaires sont nécessaires. Une valeur négative indique que le modèle peut contenir des équations redondantes ou des spécifications incohérentes.

Vérifier l'indépendance des équations

Compter les équations est insuffisant car certains bilans peuvent être linéairement dépendants. Vérifiez le rang de (A) avant de résoudre.

Pour un système carré, une matrice de plein rang est requise pour un résultat unique en utilisant un solveur direct tel que :

import numpy as np

x = np.linalg.solve(A, b)

Si la matrice est singulière ou de rang déficient, inspectez le modèle pour détecter les bilans en double, les variables omises ou les courants insuffisamment spécifiés.

Vérifier la conservation après la résolution

Substituez les débits de composants calculés dans chaque bilan. Le vecteur résiduel est :

[ \mathbf{r}=A\mathbf{x}-\mathbf{b} ]

De petits résidus indiquent que la solution numérique satisfait les équations formulées. Ils ne prouvent pas, en soi, que le modèle de procédé ou les mesures sont corrects.

Comprendre les compromis

Les inconnues de composants augmentent le nombre de variables

L'utilisation des débits de composants peut créer plus d'inconnues qu'une formulation basée sur le débit total et la composition. C'est généralement un compromis acceptable car les équations résultantes restent linéaires et sont plus faciles à résoudre et à auditer.

Les variables supplémentaires rendent également la conservation des composants explicite, ce qui est précieux dans l'analyse des unités pilotes et l'enseignement en laboratoire.

Les mesures de composition peuvent rester nécessaires

Les inconnues de débit de composant n'éliminent pas le besoin de mesurer les compositions lorsque le procédé est sous-spécifié. Elles fournissent une meilleure représentation mathématique, mais le modèle nécessite toujours suffisamment de données indépendantes et de spécifications de fonctionnement.

Les débits totaux connus ou les compositions connues peuvent être inclus comme équations linéaires là où c'est approprié.

Certains modèles de procédé sont fondamentalement non linéaires

La conservation de la matière est linéaire par rapport aux débits de composants, mais le modèle complet de l'opération unitaire peut ne pas l'être. L'équilibre liquide-vapeur, les relations pression-débit, la cinétique de réaction, les corrélations de transfert thermique et les calculs de propriétés peuvent introduire des termes non linéaires.

Utilisez un modèle de bilan linéaire là où les informations disponibles le permettent. Ne forcez pas un comportement de procédé intrinsèquement non linéaire dans une formulation linéaire simplement pour utiliser un solveur matriciel.

Les débits calculés négatifs nécessitent une enquête

Un débit de composant négatif signale généralement une erreur de signe, un courant ou un terme de procédé omis, des mesures incohérentes ou une spécification de procédé invalide. Cela peut également indiquer que le régime de fonctionnement supposé ou la structure du modèle est incompatible avec les données fournies.

Vérifiez les directions des courants, les signes des bilans, les unités, le rang des équations et les entrées mesurées avant d'accepter le résultat.

Faire le bon choix pour votre objectif

Utilisez la formulation qui correspond à l'objectif du modèle :

Si votre objectif principal est un solveur linéaire robuste : Définissez les débits massiques ou molaires de chaque composant comme vecteur inconnu et résolvez le système (A\mathbf{x}=\mathbf{b}) résultant.
Si votre objectif principal est de déterminer les compositions des courants : Résolvez d'abord les débits de composants, puis obtenez les débits totaux et les fractions massiques ou molaires par sommation et normalisation.
Si votre objectif principal est le rééquilibrage des données d'unité pilote : Incluez les informations mesurées sur les courants et les compositions comme équations indépendantes, puis vérifiez le rang de la matrice et les résidus des bilans.
Si votre objectif principal est une simulation de procédé détaillée : Conservez les variables de débit de composant pour les bilans de conservation, mais utilisez un solveur non linéaire lorsque l'équilibre, la cinétique ou les relations de propriétés l'exigent.

En choisissant les débits de composants comme inconnues principales, vous rendez les équations de conservation transparentes, préservez la linéarité là où c'est possible et donnez aux solveurs numériques le problème le plus clair à résoudre.

Tableau récapitulatif :

Élément de modélisation	Approche recommandée	Avantage
Inconnues principales	Débits massiques ou molaires de chaque composant	Préserve les équations linéaires
Modèle matriciel	Exprimer les bilans sous forme `A x = b`	Permet l'utilisation de solveurs linéaires standard
Débit total	Sommer les débits de composants après résolution	Sépare le calcul de la composition
Composition	Diviser chaque débit de composant par le débit total	Fournit les fractions massiques ou molaires
Résolubilité	Vérifier les degrés de liberté et le rang de la matrice	Confirme si une solution unique est possible
Validation	Évaluer les résidus de bilan et enquêter sur les débits négatifs	Détecte les problèmes de modélisation ou de données
Comportement non linéaire	Utiliser des solveurs non linéaires pour l'équilibre, la cinétique et les propriétés	Correspond au modèle complet du procédé

Besoin d'un modèle de bilan matière fiable pour votre laboratoire d'enseignement, votre unité pilote ou votre projet de recherche sur les procédés ? Contactez LABPARK pour discuter des unités pilotes d'opérations unitaires éducatives et professionnelles pour le génie chimique, les bioprocédés et la biotechnologie, ainsi que les applications de traitement de l'eau et de l'environnement. LABPARK aide les universités, les instituts de recherche et les entreprises à construire des flux de travail de formation et de test plus clairs grâce à des systèmes d'unités pilotes pratiques, un support de modélisation basé sur les composants et des équipements adaptés à l'enseignement en laboratoire, au développement de procédés et aux études de mise à l'échelle.