Formulation du problème
Points de réflexion :
Données
\(i \) : indice des courants chauds , \(i\in C={1,NC}\), \(NC\) est le nombre de courants chauds
\( j\) : indice des courants froid , \(j\in F={1,NF}\), \(NF\) est le nombre de courants froid
\(k\) : indice des étapes, \(k\in E={1,NE}\), \(NE\) est le nombre d'étapes
\(FCp_i\) : Produit du débit et de la capacité calorifique du courant chaud \(i,\ [kW\ K^{-1}]\)
\(FCp_j\) : Produit du débit et de la capacité calorifique du courant froid \(j,\ [kW\ K^{-1}]\)
\(T_i^e\) : Température d'entrée du courant chaud \(i\)
\(T_i^s\) : Température de sortie du courant chaud \(i\)
\(T_j^e\) : Température d'entrée du courant froid \(j\)
\(T_j^s\) : Température de sortie du courant froid \(j\)
\(T^e_{eau}\) : Température d'entrée de l'utilité froide (eau)
\(T^s_{eau}\) : Température de sortie de l'utilité froide (eau)
\(T_{vap}\) : Température d'entrée et de sortie de l'utilité chaude (vapeur saturée)
\(qc_i^{tot}=\left(T_i^e-T_i^s \right)FCp_i\) : Puissance totale échangée par le courant chaud \(i\)
\(qf_j^{tot}=\left(T_j^s-T_j^e \right)FCp_j\) : Puissance totale échangée par le courant froid \(j\)
\(h_i\) : Coefficient de transfert de chaleur du courant chaud \(i,\ kW\ m^{-2}\ K^{-1}\)
\(h_j\) : Coefficient de transfert de chaleur du courant froid \(j,\ kW\ m^{-2}\ K^{-1}\)
\(h_{eau}\) : Coefficient de transfert de chaleur de l'utilité froide \(kW\ m^{-2}\ K^{-1}\)
\(h_{vap}\) : Coefficient de transfert de chaleur de l'utilité chaude \(kW\ m^{-2}\ K^{-1}\)
\(U_{i,j}=\frac{1}{\frac{1}{h_i}+\frac{1}{h_j}}\) : Coefficient global d'échange entre \(i\) et \(j\)
\(UF_{i}=\frac{1}{\frac{1}{h_i}+\frac{1}{h_{eau}}}\): Coefficient global d'échange entre \(i\) et l'utilité froide
\(UC_{j}=\frac{1}{\frac{1}{h_j}+\frac{1}{h_{vap}}}\): Coefficient global d'échange entre \(j\) et l'utilité chaude
\(C_{aire}\) : Coût de l'air des échangeurs
\(C_{fix}\) : Coût fixe des échangeurs
\(C_{eau}\) : Coût de l'utilité froide
\(C_{vap}\) : Coût de l'utilité chaude
\(dT^{min},dT^{max}\) : Valeur minimale et maximale des différences de température
Variables
Variables libres :
\(Z\) : Fonction objectif
Variables positives [Borne inf ; Borne sup] :
\( q_{i,j,k}\) : Puissance échangée, à l’étape \(k\), entre le courant chaud \(i\) et le courant froid \(j\)
\(quf_i\) : Puissance échangée entre le courant chaud \(i\) et l’utilité froide
\(quc_j\) : Puissance échangée entre le courant froid \(j\) et l’utilité chaude
\(T^e_{i,k}\) : Température d’entrée du courant chaud \(i\) dans l’étape \(k\)
Intervalle de définition \(\left[T_i^s ;T_i^e\right]\)
\(T^s_{i,k}\) : Température de sortie du courant chaud \(i\) dans l’étape \(k\)
Intervalle de définition \(\left[T_i^s ;T_i^e\right]\)
\(T^e_{j,k}\) : Température d’entrée du courant froid \(j\) dans l’étape \(k\)
Intervalle de définition \(\left[T_j^s ;T_j^e\right]\)
\(T^s_{j,k}\) : Température de sortie du courant froid \(j\) dans l’étape \(k\)
Intervalle de définition \(\left[T_j^s ;T_j^e\right]\)
\(dT^g_{i,j,k}\) : Différence de température, à gauche de l’étape \(k\), entre le courant chaud \(i\) et le courant froid \(j\)
Intervalle de définition \(\left[dT^{min} ;dT^{max}\right]\)
\(dT^d_{i,j,k}\) : Différence de température, à droite l’étape \(k\), entre le courant chaud \(i\) et le courant froid \(j\)
Intervalle de définition \(\left[dT^{min} ;dT^{max}\right]\)
\(dTuf_i\) : Différence de température, à gauche de l’échangeur, entre le courant chaud \(i\) et l’utilité froide
Intervalle de définition \(\left[dT^{min} ;dT^{max}\right]\)
\(dTuc_j\) : Différence de température, à droite de l’échangeur, entre le courant froid \(j\) et l’utilité chaude
Intervalle de définition \(\left[dT^{min} ;dT^{max}\right]\)
Variables binaires :
\(y_{i,j,k}\) : Couplage, à l’étape \(k\), entre le courant chaud \(i\) et le courant froid \(j\)
\(yuf_i\) : Couplage entre le courant chaud \(i\) et l’utilité froide
\(yuc_j\) : Couplage entre le courant froid \(j\) et l’utilité chaude
Contraintes
Bilan global sur chaque courant :
Courants chauds
\(qc_i^{tot} = \underset{k\in E}{\sum} \underset{j\in F}{\sum} q_{i,j,k} + quf_i\) avec \(i\in C\)
Courants froids
\(qc_j^{tot} = \underset{k\in E}{\sum} \underset{i\in C}{\sum} q_{i,j,k} + quc_j\) avec \(j\in F\)
Bilan par étape et par courant :
Courants chauds
\(\left( T^e_{i,k}-T^s_{i,k}\right)FCp_i = \underset{j\in F}{\sum} q_{i,j,k}\) avec \(i\in C, k\in E\)
Courant froids
\(\left( T^s_{j,k}-T^e_{j,k}\right)FCp_j = \underset{i\in C}{\sum} q_{i,j,k}\) avec \(j\in F, k\in E\)
Bilan sur les utilités :
Utilités froides
\(\left(T^s_{i,NE}-T^s_i \right)FCp_i=quf_i\) avec \(i\in C\)
Utilités chaudes
\(\left(T^s_j-T^s_{j,1}\right)FCp_j=quc_j\) avec \(j\in F\)
Connexion des températures d'entrée :
Courants chauds
\(T^s_{i,k} = T^e_{i,k+1}\) avec \(i\in C, k=[1 ; NE-1]\)
Courants froids
\(T^e_{j,k} = T^s_{j,k+1}\) avec \(j\in F, k=[1 ; NE-1]\)
Affectation des températures d'entrée :
Courants chauds
\(T^e_{i,1}=T^e_i\) avec \(i\in C\)
Courants froids
\(T^e_{j,NE}=T^e_j\) avec \(j\in F\)
Couplage :
Chaud/Froid
\(q_{i,j,k} \le y_{i,j,k} Max\left( qc_i^{tot}, qf_j^{tot}\right)\) avec \(i\in C, j\in F, k \in E\)
Chaud/utilité
\(quf_i\le yuf_i.qc_i^{tot}\) avec \(i\in C\)
Froid/utilité
\(quc_j\le yuc_j.qf_j^{tot}\) avec \(j\in F\)
Topologie - un seul couplage par étape :
courants chauds
\(\underset{j\in F}{\sum}y_{i,j,k}\le 1\) avec \(i\in C, k\in E\)
courants froids
\(\underset{i\in C}{\sum}y_{i,j,k}\le 1\) avec \(j\in F, k\in E\)
Différence de température à gauche (cas général et utilité froide) :
\(dT^g_{i,j,k} \le T^e_{i,k}-T^s_{j,k}+dT^{max}(1-y_{i,j,k})\) avec \(i\in C, j\in F, k\in E\)
\(dTuf_i\le T^s_{i,NE}-T^s_{eau}+dT^{max}(1-yuf_i)\) avec \(i\in C\)
Différence de température à droite (cas général et utilité chaude) :
\(dT^d_{i,j,k} \le T^s_{i,k}-T^e_{j,k}+dT^{max}(1-y_{i,j,k})\) avec \(i\in C, j\in F, k\in E\)
\(dTuc_j\le T_{vap}-T^s_{j,1}+dT^{max}(1-yuc_j)\) avec \(j\in F\)
Connexion des différences de températures E/S entre les étapes ()
\(dT^d_{i,j,k}=dT^g_{i,j,k+1}\) avec \(i\in C, j\in F, k\in E\)
Fonction objectif :
\(z=C_{fix}\left(\underset{i\in C}{\sum} \underset{j\in F}{\sum}\underset{k\in E}{\sum}y_{i,j,k}+\underset{i\in C}{\sum}yuf_i+\underset{j\in F}{\sum}yuc_j \right)\\+C_{eau}\underset{i\in C}{\sum}quf_i+C_{vap}\underset{j\in F}{\sum}quc_j\\+C_{aire}\frac{\underset{i\in C}{\sum} \underset{j\in F}{\sum}\underset{k\in E}{\sum}q_{i,j,k}}{U_{i,j}LMTD\left(dT^g_{i,j,k},dT^d_{i,j,k}\right)}\\+C_{aire}\frac{\underset{i\in C}{\sum}quf_i}{U_{i,eau}LMTD\left(dTuf_i,T_i^s-T^e_{eau}\right)}\\+C_{aire}\frac{\underset{j\in F}{\sum}quc_j}{U_{j,vap}LMTD\left(T^e_{vap}-T^s_j,dTuc_j\right)}\)
avec :
\(LMTD\left(dT_1,dT_2\right)\approx \left[dT_1dT_2\frac{(dT_1+dT_2)}{2} \right]^{1/3}\)
Initialisation des variables
\(q_{i,j,k}=Min\left(qc_i^{tot}, qf_j^{tot}\right)\) avec \(i\in C, j\in F, k\in E\)
\(T^e_{i,k}=T^s_{i,k}=T^e_i\) avec \(i\in C, k\in E\)
\(T^e_{j,k}=T^s_{j,k}=T^e_j\) avec \(j\in C, k\in E\)
\(dT^g_{i,j,k}=dT^d_{i,j,k}=\left(T^e_i-T^e_j\right)\) avec \(i\in C, j\in F, k\in E\)
\(dTuf_i = \left(T^e_i-T^e_{eau} \right)\) avec \(i\in C\)
\(dTuc_j = \left(T^e_{vap}-T^e_j \right)\) avec \(j\in F\)