Nouvelle approche de l'étude thermodynamique des systèmes ouverts en PT et écoles d'ingénieurs<br/>par Jean-Noël BEURY Lycée Gustave Eiffel - 33000 BORDEAUX jnbeury@laposte.net _________________ RÉSUMÉ La méthode actuellement enseignée en thermodynamique industrielle pour prendre en compte les irréversibilités dans un système est basée sur des hypothèses implicites qui sont en contradiction avec le premier principe. Nous montrerons les incohérences de cette méthode A dans le cas d'un échangeur thermique et les nombreuses erreurs d'interprétation physique. Nous proposerons une nouvelle approche (méthode B) de l'étude des systèmes ouverts en se recentrant sur les fondamentaux de la thermodynamique. Cette méthode est entièrement conforme aux nouveaux programmes de PT et peut être appliquée dans les écoles d'ingénieurs, BTS, IUT... 1. PRÉSENTATION ET UTILISATION DE LA MÉTHODE A 1.1. Présentation de la méthode A actuellement enseignée En thermodynamique industrielle, on travaille très souvent avec des systèmes ouverts. Les causes d'irréversibilité ont plusieurs origines : détente dans le vide, surcharge brutale sur un piston, diffusion de particules, réactions chimiques, irréversibilité mécanique (frottements du fluide sur la paroi, frottements fluides internes), irréversibilité thermique (transfert thermique entre des systèmes à des températures différentes)... Dans les Classes Préparatoires aux Grandes Ecoles PT et dans les écoles d'ingénieurs, l'irréversibilité en thermodynamique industrielle est très souvent traitée en utilisant qf . On rencontre dans la littérature les présentations suivantes pour introduire qf : 1.1.1. Première présentation (référence [1]) On envisage une transformation irréversible de l'état 1 à l'état 2 d'un système ouvert en régime permanent d'écoulement, à une entrée et une sortie. On définit (voir annexe 1) : ·  w i = travail indiqué massique (travail reçu par unité de masse de fluide de la · · part des parties mobiles de la machine).  qe = transfert thermique massique reçu (échangé avec l'extérieur) à travers toute la frontière du système ouvert.  qf = quantité de chaleur massique dissipée par frottement ( qf  0 )<br/><br/>Page 1<br/><br/>Jean-Noël BEURY<br/><br/><br/>·<br/><br/> w f = travail massique des forces de frottement ( w f  0 )<br/><br/>Le travail massique des forces intérieures de viscosité (forces de frottement) est transformé intégralement en chaleur massique de frottement :  w f +  qf = 0 . On remplace la transformation irréversible par une transformation réversible respectant les conditions suivantes : · même loi de variation que la transformation réelle. · chaleur massique de frottement  qf supposée fournie à température variable (sans irréversibilité thermique) par une source fictive pour respecter la condition précédente. Dans ces conditions, la chaleur reçue par le fluide est :  q =  qe +  qf . La variation d'entropie massique pour une transformation réversible est : ds =<br/><br/>q<br/>T<br/><br/>=<br/><br/> qe +  q f<br/>T<br/><br/>.<br/><br/>dh + dec + dep =  w i +  qe<br/><br/>On a 3 équations : dh = Tds + vdp  q +  qf ds = e et  qf  0 T La première équation est l'expression du premier principe de la thermodynamique pour un système ouvert en régime permanent, à une entrée et une sortie (démonstration Annexe 1). La deuxième équation est l'identité thermodynamique qui est universelle : elle ne traduit que les propriétés thermodynamiques du fluide dans ses différents états. 1.1.2. Deuxième présentation (référence [2])<br/><br/>Le long d'une transformation réversible, on a par définition : ds =<br/><br/>. T Le long d'une transformation irréversible, on a recours à un « tour de passe-passe » en assimilant la chaleur générée par les frottements à un apport de l'extérieur  qf .<br/>On pose ds =<br/><br/>q<br/><br/> qe +  qf<br/>T<br/><br/>=  se +  sf avec  qf &gt; 0   sf &gt; 0 .<br/><br/>1.1.3. Questions découlant immédiatement de ces deux présentations<br/><br/>· ·<br/><br/> qf est fourni par une source extérieure. Pourquoi ne pas rajouter  qf dans le<br/>premier principe alors qu'on le rajoute dans le second principe ? Pourquoi  qe garderait-il la même valeur pour la transformation irréversible et la transformation réversible alors que  qe dépend du chemin suivi ?<br/><br/>Dans l'épreuve de PT 2002 (voir référence [8]), question 4 : « Calculer la création d'entropie massique due à l'irréversibilité de l'évolution dans le compresseur ». L'étudiant q doit donc calculer  sc = f . T<br/><br/>Nouvelle approche de l'étude thermodynamique des systèmes ouverts en PT.<br/><br/>Page 2<br/><br/><br/>On ne peut pas identifier l'entropie échangée massique à<br/><br/>puisque pour une T transformation monotherme (système en contact avec une source de chaleur à la  qe température TS) elle vaut . De gros risques de confusions apparaissent !!! De TS  qf même, l'entropie créée massique ne peut pas être identifiée à . T 1.1.4. Ces deux présentations sont incorrectes !!! Elles ne respectent pas les principes de base de la thermodynamique : le transfert thermique dépend du chemin suivi. On devrait tenir compte de tous les transferts thermiques échangés avec l'extérieur pour écrire le premier principe de la thermodynamique. Pourquoi ne pas rajouter  qf<br/><br/> qe<br/><br/>· ·<br/><br/>·<br/><br/>dans le premier principe ? Il est donc faux d'affirmer avec les présentations précédentes que « tout se passe comme si le fluide recevait un transfert thermique d'une source fictive extérieure correspondant à l'échauffement du fluide » puisqu'il faudrait rajouter  qf dans le premier principe de la thermodynamique.<br/>1.2. Utilisation de la méthode A<br/><br/>1.2.1. Echangeur thermique calorifugé On considère un échangeur thermique parfaitement calorifugé à contre courant en régime permanent. On note T1 et T1' les températures, supposées connues, des fluides à l'entrée de l'échangeur, T2 et T2' sont respectivement les températures à la sortie. Le système ouvert étudié est un élément de longueur dx de la conduite 12. On note T(x) la température à l'abscisse x. La température à l'abscisse x de la conduite 1'2' est notée T'(x). On fait l'étude en régime permanent. Il n'y a pas de partie mobile de la machine :  w i = 0 . On néglige les variations d'énergie cinétique et potentielle.<br/><br/>T2' T1 Dm, cp O<br/><br/>Dm', cp'<br/><br/>T1' T2 L<br/><br/>x<br/><br/>x+dx<br/><br/>x<br/><br/>Les 3 équations de la méthode A s'écrivent :<br/><br/>Page 3<br/><br/>Jean-Noël BEURY<br/><br/><br/>  qe  qf dh  qf  dh + dec + dep =  w i +  qe ds = T + T = T + T   On a donc :  . dh = Tds + vdp dh v ds =  - dp  q +  qf  ds = e  TT T  vdp On en déduit que  qf = -vdp et  sc = - T Dans un échangeur thermique (voir annexe 2), on rencontre deux types d'irréversibilité : · irréversibilité mécanique (due à la viscosité du fluide). Elle se traduit par une baisse de pression (pertes de charge). On peut montrer en mécanique des fluides que la viscosité des fluides peut s'interpréter par un phénomène de diffusion de quantité de mouvement engendré par un gradient de vitesse. · irréversibilité thermique (diffusion de chaleur) : écart entre les températures des deux conduites. Le transfert thermique se fait spontanément du corps le plus chaud vers le corps le plus froid. Ces deux types d'irréversibilité complètement indépendants sont dus aux phénomènes diffusifs (diffusion de quantité de mouvement et diffusion de chaleur). En pratique, la baisse de pression n'est pas négligeable par rapport aux autres pertes de charge dans un circuit hydraulique mais le transfert thermique mis en jeu est négligeable par rapport au flux de chaleur échangé entre le fluide froid et le fluide froid. On a également une irréversibilité thermique (voir annexe 2 où il faudrait avoir un échangeur de très grande longueur). La méthode A donne  qf = 0 et  sc = 0 alors qu'on a une réversibilité mécanique<br/>et une irréversibilité thermique !!! 1.2.2. Utilisation avec la méthode A de deux chemins pour aller de l'état 1 à l'état 2 Actuellement, on utilise avec la méthode A deux chemins différents pour aller du même état initial au même état final : Dans la référence [7] (Banque PT 1998) Question 6 : « Pour caractériser une compression adiabatique, on définit son rendement indiqué polytropique, noté ikc<br/><br/>comme étant le rapport entre le travail indiqué massique, noté wikc, de la compression assimilée à une évolution polytropique réversible (non adiabatique) et le travail indiqué massique de la compression déterminé lors de la question 4a. »<br/>La transformation fictive réversible est donc déjà utilisée avec la méthode A pour calculer le travail de transvasement (voir paragraphe 3.5).<br/><br/>Nouvelle approche de l'étude thermodynamique des systèmes ouverts en PT.<br/><br/>Page 4<br/><br/><br/>1.3. Méthode A confrontée aux principes de la thermodynamique 1.3.1. Définition de qf On considère un système ouvert en régime permanent d'écoulement, à une entrée et une sortie en contact avec une source de chaleur à la température TS : la transformation est monotherme puisqu'on a une seule source de chaleur. a) Le premier principe de la thermodynamique pour un système ouvert s'écrit (voir annexe 1 pour la démonstration) : (u2 - u1 ) + ( ep 2 - ep1 + ec 2 - ec1 ) = w i + p1v1 - p2v 2 + qe<br/><br/>En notation infinitésimale, on a : du + dec + dep =  w i - d( pv ) +  qe<br/>En posant h = u + pv , on obtient : dh + dec + dep =  w i +  qe avec<br/><br/> w i = travail indiqué massique (travail massique reçu par le fluide de la part des parties<br/><br/>mobiles de la machine).  qe = transfert thermique massique algébriquement reçu (échangé avec l'extérieur) à travers toute la frontière du système ouvert. b) Le théorème de l'énergie cinétique s'écrit : dec =  w ext +  w int TRAVAIL MASSIQUE DES FORCES EXTÉRIEURES : · Le travail massique du fluide amont et aval (voir démonstration annexe 1 avec p1v1 ­ p2v2) est égal à ­ d(pv). · Le travail indiqué massique (travail massique reçu par le fluide de la part des parties mobiles de la machine) est égal à  w i . Le travail des forces extérieures dérivant d'une énergie potentielle est à égal à ­dep. TRAVAIL MASSIQUE DES FORCES INTÉRIEURES :  w int = travail massique des forces intérieures = travail massique des forces intérieures de pression (pdv : voir les références [3] et [4] pour la démonstration de l'expression des forces intérieures de pression) + travail massique des forces intérieures de viscosité ( w f  0 ). On a donc : dec = {-d( pv ) +  w i - dep } +  w int = - pdv - vdp +  w i - dep + pdv +  w f ·<br/><br/>{<br/><br/>}<br/><br/>c) L'identité thermodynamique s'écrit : dh = Tds + vdp dh + dec + dep =  w i +  qe (1)  On obtient finalement : dh = Tds + vdp (2) de + de = -vdp +  w +  w (3) p i f c<br/><br/>Des équations (1) et (3), on a : dh - vdp +  w i +  w f =  w i +  qe En utilisant (2), on a : Tds + vdp +  w f -vdp =  qe Soit Tds +  w f =  qe , d'où ds = (4) T T Le bilan d'entropie pour une transformation monotherme s'écrit :<br/>Page 5 Jean-Noël BEURY<br/><br/> qe<br/><br/>-<br/><br/> wf<br/><br/><br/>ds =  se +  sc =<br/><br/> qe<br/>TS<br/><br/>+  sc (5)<br/><br/>En identifiant les équations (4) et (5), on a : (T - T ) -  w f  q  q w  sc = e - e - f =  qe S T TS T T TS T<br/><br/>  q   (TS - T ) -  w f   Le bilan d'entropie s'écrit donc : ds = { se } + { sc } =  e  +  qe  (6) T TS T   TS   <br/>Interprétation physique : · le premier terme traduit l'entropie échangée : transfert thermique avec une source de chaleur à la température TS. · Le deuxième terme traduit l'entropie créée qui est bien la somme de deux (T - T )  0 irréversibilité  wf contributions : -  0 irréversibilité mécanique et  qe S T T TS<br/><br/>thermique à cause d'un transfert thermique entre deux corps à des températures différentes. On définit  qf le transfert thermique massique dû aux frottements comme l'opposé<br/>du travail des forces intérieures de viscosité :  qf = - w f . L'interprétation physique de  w f est évidente (travail des forces intérieures de viscosité) alors que l'interprétation physique de  qf et de<br/><br/> qf<br/>T<br/><br/>est beaucoup moins évidente comme on<br/><br/>va le voir... L'équation (6) devient alors :   (T - T )  +   qf  =   qe  +   qf   q ds = { se } + { sc } =  e +  qe S     T TS   T   T   T     TS<br/><br/>de pression) = - ( -d ( pv ) + pdv ) = pdv + vdp - pdv = vdp .<br/><br/>On définit le travail de transvasement massique par :  w t = ­ (travail massique du fluide amont et aval + travail massique des forces intérieures<br/><br/>Interprétation physique dans le cas d'un compresseur avec dec = dep = 0. Le théorème de l'énergie cinétique s'écrit : dec + dep = -vdp +  w i +  w f , d'où  w i =  w t -  w f , soit<br/><br/> w i   w t . Le travail de transvasement correspond donc au travail minimum à fournir au<br/>compresseur pour vaincre le travail des forces de pression amont et aval, ainsi que le travail des forces intérieures de pression. On verra au paragraphe 3.5 une interprétation physique plus simple du travail de transvasement<br/>1.3.2. Erreurs d'interprétation fréquemment rencontrées q q q q  se = e et  sc = f . Le premier terme ds = e + f =  se +  sc avec T T T T représenterait l'entropie échangée et le deuxième terme l'entropie créée. C'est faux : le<br/><br/>·<br/><br/>Nouvelle approche de l'étude thermodynamique des systèmes ouverts en PT.<br/><br/>Page 6<br/><br/><br/>·<br/><br/>·<br/><br/>représente l'irréversibilité uniquement d'origine mécanique alors que le T  qe terme n'a pas d'interprétation physique simple (il contient en fait l'entropie T échangée et l'irréversibilité d'origine thermique). Les énoncés de concours et les présentations de la méthode A (voir paragraphe 1.1) suggèrent cette relation fausse : - Voir référence [7] PT 1998 question 5.c : « Calculer qf la quantité de chaleur massique due à l'irréversibilité ». Il faudrait plutôt écrire : « calculer qf la quantité de chaleur massique due à l'irréversibilité mécanique ». - Voir référence [8] PT 2002 question 4 : « Calculer la création d'entropie massique due à l'irréversibilité de l'évolution dans le compresseur ». L'énoncé demande donc q à l'étudiant de calculer  sc = f . T Il faudrait plutôt écrire : « Calculer la création d'entropie massique due à l'irréversibilité mécanique de l'évolution dans le compresseur ». On rencontre souvent  qf = 0  transformation réversible . Cette équivalence est fausse puisque  qf ne mesure que l'irréversibilité d'origine mécanique. On peut donc avoir une transformation irréversible avec  qf = 0 ! Il n'est donc plus étonnant d'avoir qf = 0 dans l'échangeur thermique du paragraphe 1.2.1 avec dp = 0 mais par contre, il est faux d'écrire que l'entropie créée est nulle. En réalité, dans cet exemple, on a :  q   (T - T ) -  w f  =  qe = dh . L'entropie créée vaut   ds = { se } + { sc } =  e  +  qe S  T TS T T T  TS     donc :  sc =  qe<br/><br/>terme<br/><br/> qf<br/><br/>(TS - T )<br/>T TS<br/><br/>!!! Ce résultat est complètement masqué avec la méthode A.<br/><br/>·<br/><br/>·<br/><br/>·<br/><br/>1.3.3. Problèmes pédagogiques liés à l'utilisation de qf La méthode A est très délicate à introduire correctement car la démonstration faite dans le paragraphe 1.3.1 est hors de propos pour les élèves : le travail des forces de pression intérieure (avec risques de confusion entre p dv et ­p dv ) et le théorème de l'énergie cinétique pour un système ouvert sont hors programme en PT. On se contente souvent des présentations faites dans le paragraphe 1.1 qui contiennent de nombreuses contradictions. Les étudiants ont l'impression que la thermodynamique est un tour de passe-passe alors qu'au contraire elle est rigoureuse. Sans précaution, on peut très facilement faire des erreurs d'interprétation physique. La rédaction même des problèmes de concours (cf références [7] et [8]) laisse supposer q q que  se = e et  sc = f . T T<br/><br/>Page 7<br/><br/>Jean-Noël BEURY<br/><br/><br/>· ·<br/><br/>Utiliser simultanément l'entropie créée et  qf pose problème puisque  qf masque complètement l'entropie créée d'origine thermique. On ne peut pas interpréter qf par la phrase : « Tout se passe comme si le fluide recevait un transfert thermique massique qf correspondant à l'échauffement du fluide<br/><br/>·<br/><br/>à cause des frottements fluides internes (forces de viscosité) » car elle fait référence au paragraphe 1.1 qui est faux. Cette interprétation est pourtant facilement compréhensible par les étudiants !!! Un autre problème pédagogique apparaît avec l'utilisation de  qf : Les étudiants se trouvent confrontés à deux relations pour écrire le bilan d'entropie : ds =<br/>ds =  se +  sc =<br/><br/> qe<br/>T<br/><br/>+<br/><br/> qf<br/>T<br/><br/>et<br/><br/> qe<br/>TS<br/><br/>+  sc (pour une transformation monotherme). Laquelle choisir<br/><br/>dans un exercice ?<br/>2. QUELQUES REMARQUES CONCERNANT LA PROGRAMMES DE THERMODYNAMIQUE EN CPGE RÉFORME DE 1995 DES<br/><br/>2.1. Présentation de la réforme de 1995 L'enseignement de la thermodynamique en MP, PC et PSI a beaucoup évolué en 1995 suite à la réforme des programmes inspirés par Hubert Gié, Inspecteur Général de q l'Education Nationale : la relation ds = rev a été remplacée par l'identité T thermodynamique pour les calculs d'entropie. L'entropie échangée et l'entropie créée se trouvent au coeur du bilan d'entropie. La compréhension des phénomènes physiques s'en trouve facilitée pour les étudiants. Les équations de base sont devenues le premier principe de la thermodynamique, le second principe de la thermodynamique et l'identité thermodynamique (cf référence [6]). dh + dec + dep =  w i +  qe  Pour un système ouvert, on a donc : ds =  se +  sc dh = Tds + vdp <br/><br/>Prenons le cas d'un gaz parfait pour illustrer le bilan d'entropie qui se fait très simplement : 1) Calcul de la variation d'entropie : dh = Tds + vdp = c p dT , d'où ds = c p<br/>T  p  dT v dT dp et s = s2 - s1 = c p ln  2  - r ln  2  - dp = c p -r TT T p  T1   p1 <br/><br/>2) Calcul de l'entropie échangée. Elle est nulle si la transformation est adiabatique. Pour q une transformation monotherme, on a  se = e . TS 3) On en déduit immédiatement l'entropie créée massique :  sc = ds -  se avec son interprétation physique immédiate : l'entropie créée est strictement positive pour une transformation irréversible et nulle pour une transformation réversible.<br/><br/>Nouvelle approche de l'étude thermodynamique des systèmes ouverts en PT.<br/><br/>Page 8<br/><br/><br/>2.2. Programmes de PT 2004 On trouve dans les nouveaux programmes de PT 2004 (cf référence [9]) les phrases suivantes : · « Formulation du second principe pour un système ouvert simple en régime permanent, à une entrée et une sortie : ds =  se +  sc ».<br/><br/>« La représentation de l'entropie générée par les frottements n'est pas exigible ». q q La relation ds = e + f ne figure donc pas explicitement au programme. T T On retrouve donc la même approche dans la filière PT que dans les filières MP, PC et PSI. La connaissance de qf et de l'entropie générée par les frottements est importante pour la description d'une machine adiabatique irréversible en thermodynamique industrielle.<br/>Le programme de PT met au premier plan l'entropie échangée et l'entropie créée pour les systèmes ouverts. Nous allons voir que l'on peut se contenter de la relation q ds =  se +  sc = e +  sc et ne pas aborder la discussion délicate avec TS q q ds = e + f pour traiter les systèmes ouverts. On pourra même reprendre T T l'interprétation physique de la méthode A !!! 3. MISE EN PLACE DE LA MÉTHODE B<br/><br/>·<br/><br/>Soit un système ouvert en régime permanent d'écoulement, à une entrée et une sortie.<br/>3.1. Rappel des équations de la méthode A La méthode A, actuellement utilisée est basée sur les équations suivantes : dh + dec + dep =  w i +  qe<br/><br/>dh = Tds + vdp  q +  qf ds = e T  w t = vdp<br/>Le bilan d'entropie pour tous les systèmes ouverts se fait donc avec la relation : q  q +  qf ds = rev = e qui n'est pas conforme au programme PT 2004 (cf référence T T [9]) puisqu'on utilise systématiquement une transformation réversible ! On perd également l'entropie échangée et l'entropie créée. La technique utilisant une transformation réversible était utilisée systématiquement pour les systèmes fermés avant la réforme de 1995 !<br/><br/>Page 9<br/><br/>Jean-Noël BEURY<br/><br/><br/>3.2. Mise en place de la méthode B 3.2.1. Cas général Pour la transformation réelle (réversible ou irréversible), on peut écrire le premier principe (cf annexe 1 pour la démonstration), le second principe de la thermodynamique, l'identité thermodynamique et le travail de transvasement massique. dh + dec + dep =  w i +  qe<br/>ds =  se +  sc dh = Tds + vdp<br/><br/> w t = vd p<br/><br/> w i est le travail indiqué massique pour la transformation réelle (travail massique<br/>reçu par le fluide de la part des parties mobiles de la machine).  qe est le transfert thermique massique reçu (échangé avec l'extérieur) à travers toute la frontière du système ouvert. Dans l'annexe 1, on a vu que l'on peut se ramener à un système fermé. Dans ce cas, l'entropie échangée est le transfert thermique divisé par T ' la température de la surface<br/><br/>d'échange (c'est-à-dire la surface (  ' ) du volume de contrôle) à travers laquelle se fait le transfert thermique :<br/><br/> Se = - <br/>'<br/><br/>S j .dSext dt T '<br/><br/>L'avantage de cette définition est d'avoir une seule expression de l'entropie échangée, même si les deux systèmes qui interagissent ne sont pas à l'équilibre thermique. On suppose simplement l'équilibre thermodynamique local réalisé au voisinage de la surface. On rencontre très souvent les cas particuliers suivants : · Pour une transformation adiabatique :  Se = 0 et donc  se = 0 .<br/>·<br/><br/>Le système échange de la chaleur avec une source de chaleur (voir référence [6] : système fermé n'échangeant aucun travail et capable de fournir un transfert thermique sans que sa température TS varie. La transformation est monotherme.  Qe L'entropie échangée vaut :  Se = . En grandeurs massiques, on a : TS<br/><br/> se =<br/>·<br/><br/> qe<br/>TS<br/><br/>Si le système échange de la chaleur avec un autre système dont la température varie, on peut utiliser la notion de pseudo-source de chaleur (voir référence [6]).<br/><br/>L'entropie créée est notée  Sc . En grandeurs massiques, on définit  sc l'entropie créée massique. L'entropie créée massique est calculée entre deux états d'équilibre du système par sa définition :  sc = ds -  se . Interprétation physique : ·  sc &gt; 0 pour une transformation irréversible<br/>·<br/><br/> sc = 0 pour une transformation réversible.<br/>Page 10<br/><br/>Nouvelle approche de l'étude thermodynamique des systèmes ouverts en PT.<br/><br/><br/>En thermodynamique industrielle, on a très souvent dec et dep nuls ou négligeables sauf pour les gaz en sortie de tuyère. q  q +  qf On peut remarquer que la relation ds = rev = e a été remplacée par l'identité T T thermodynamique pour les calculs d'entropie (cf paragraphe 2.1). L'enthalpie est un paramètre énergétique essentiel en thermodynamique industrielle. L'entropie échangée et l'entropie créée se trouvent au coeur de l'étude du système (conforme au programme PT 2004 cf référence [9]).<br/>Ces 4 relations suffisent pour presque toutes les applications réversibles ou irréversibles rencontrées en thermodynamique industrielle (avec les hypothèses : système ouvert en régime permanent d'écoulement, à une entrée et une sortie). Dans certaines transformations adiabatiques irréversibles (compresseurs, turbines...), on tient compte des irréversibilités mécaniques dues aux forces de frottements internes (forces de viscosité). Elles se traduisent par un travail des forces intérieures de viscosité ( w f  0 ) . Ce travail des forces intérieures se manifeste par une dissipation d'énergie et donc un échauffement du fluide. On définit  qf le transfert thermique massique dû aux frottements comme l'opposé du travail massique des forces intérieures de viscosité  qf = - w f . On cherchera donc dans ces applications à calculer qf  0 et l'entropie générée par les frottements. On va voir une technique permettant de répondre à ces deux questions. 3.2.2. Cas particulier d'une transformation adiabatique irréversible où le calcul de qf est nécessaire On ne connaît pas en général le chemin effectivement suivi par la transformation réelle. On cherche parfois à la modéliser par une autre transformation (on verra que c'est souvent une polytropique) avec deux chemins différents pour aller du même état initial au même état final. L'objectif est de calculer qf et l'entropie générée par les frottements. Dans ce cas, on mettra en évidence deux chemins (transformation réelle irréversible et transformation FICTIVE réversible) pour aller du même état initial au même état final. On a les mêmes dh, dec, dep, ds car h, ec, ep et s sont des fonctions d'état et leurs variations ne dépendent pas du chemin suivi.<br/><br/>·<br/><br/>Pour la transformation REELLE irréversible, on écrit (comme dans le paragraphe 3.2.1) le premier, le second principe de la thermodynamique ainsi que l'identité thermodynamique dh + dec + dep =  w i +  qe ds =  se +  sc dh = Tds + vdp<br/><br/>( w t<br/><br/>= vdp )<br/><br/>Page 11<br/><br/>Jean-Noël BEURY<br/><br/><br/>Il est très souvent inutile d'écrire le travail de transvasement massique car on ne connaît pas en général le chemin effectivement suivi par la transformation réelle.<br/>·<br/><br/>Pour la transformation FICTIVE réversible, on écrit l'identité thermodynamique, le premier principe, le second principe de la thermodynamique et éventuellement (si nécessaire) le travail de transvasement fictif. Pour cette transformation FICTIVE réversible, on définit un travail indiqué massique réversible  w i rev noté  w if et un transfert thermique massique réversible  qe rev noté  qef . On note :<br/><br/> w if = travail indiqué massique fictif (calculé sur la transformation FICTIVE réversible)  w tf = travail de transvasement massique fictif (calculé sur la transformation FICTIVE<br/>réversible)  qef = transfert thermique massique fictif (calculé sur la transformation FICTIVE<br/>réversible).<br/><br/>dh + dec + dep =  w if +  qef dh = Tds + vdp q ds = ef T  w tf = vdp Souvent la transformation fictive est une polytropique, on remplacera l'indice f par l'indice k comme c'est l'usage en thermodynamique industrielle. Cette notation a déjà été utilisée avec la méthode A : voir paragraphe 1.2.2.<br/>3.3. Avantages d'utiliser la méthode B On a une présentation unitaire de la thermodynamique pour les systèmes fermés et ouverts en gardant les fondamentaux de la thermodynamique : premier principe, second principe, identité thermodynamique, entropie échangée et entropie créée. Les étudiants peuvent comprendre très facilement les notions introduites. La démonstration du premier principe de la thermodynamique pour un système ouvert en régime permanent à une entrée et une sortie est au programme (cf annexe 1). La transformation fictive n'est utilisée que lorsque c'est absolument nécessaire (modélisation par une polytropique par exemple). Les notations employées pour la transformation fictive sont transposées de la transformation réelle. Il n'y a pas de confusion possible. On dissocie bien la transformation réelle de la transformation FICTIVE. 3.4. Interprétation physique 3.4.1. Interprétation physique de qf Rappelons que l'on applique cette transformation fictive réversible uniquement pour des transformations réelles adiabatiques :<br/><br/>·<br/><br/>·<br/><br/>·<br/><br/>Nouvelle approche de l'étude thermodynamique des systèmes ouverts en PT.<br/><br/>Page 12<br/><br/><br/>T T Les deux grandeurs  qf et  qef sont donc égales puisque la variation d'entropie ne<br/><br/>La méthode A donne : ds =<br/><br/> qe +  qf<br/><br/>. La méthode B donne : ds =<br/><br/> qef<br/><br/>dépend pas du chemin suivi.<br/>Pour répondre aux questions des problèmes de concours concernant qf, il suffit donc : · de calculer qef avec la méthode B · et d'écrire à la fin des calculs que qf = qef . INTERPRETATION PHYSIQUE DE qf = qef : tout se passe comme si le fluide recevait un transfert thermique massique qf correspondant à l'échauffement du fluide à cause des frottements fluides internes (forces de viscosité).<br/><br/>On a vu que interprétation était incorrecte avec la méthode A alors qu'elle devient correcte avec la méthode B !<br/>3.4.2. Quelle est l'entropie générée par les frottements ?<br/><br/>Pour la transformation réelle, on a : ds =  se +  sc . Pour la transformation fictive réversible, on a : ds =<br/><br/> qef<br/>T<br/><br/>.<br/><br/>On a donc qf = qef &gt; 0 et une entropie créée strictement positive. C'est tout à fait cohérent. La seule cause d'irréversibilité ici est d'origine mécanique puisque la transformation est adiabatique.<br/>3.5. Application de la méthode B à un compresseur réel adiabatique irréversible. 3.5.1. Compresseur adiabatique réversible Un compresseur est un organe mécanique destiné à augmenter par un p2 procédé uniquement mécanique la pression d'un fluide. On veut faire passer la pression du fluide de p1 à p2 avec p2 &gt; p1. Les compresseurs sont utilisés dans de nombreuses applications industrielles : réfrigération, compresseur à air. Les compresseurs peuvent être p1 Schéma d'un regroupés en deux grandes catégories : compresseur · les compresseurs volumétriques : le fluide est emprisonné dans un volume fermé que l'on réduit pour réaliser la compression. · les turbocompresseurs ou compresseurs à aubages. On a un mouvement continu du fluide. La compression est obtenue en convertissant l'énergie cinétique communiquée au fluide en pression par des aubages mobiles. On a donc un système ouvert (machine avec transvasement). Le système considéré sera un gaz.<br/><br/>Page 13<br/><br/>Jean-Noël BEURY<br/><br/><br/>On étudiera dans cet article uniquement les turbocompresseurs. Ce sont des machines compactes avec des surfaces d'échange faibles et des coefficients d'échange thermique faibles (pour un gaz). On suppose donc très souvent la compression adiabatique. Dans ce paragraphe, on suppose la transformation réversible. On néglige les variations d'énergie cinétique et d'énergie potentielle. On appelle par définition le travail de transvasement massique : w t =<br/>def<br/><br/><br/><br/>2<br/><br/>1<br/><br/>v dp .<br/><br/> dh + dec + dep =  w i +  qe   Les relations fondamentales s'écrivent : dh = Tds + vdp  ds =  qe   T dh =  w i On a donc :   d h = vd p<br/><br/>On en déduit immédiatement que  w i = vdp =  w t et h = w i = w t =  vdp<br/>1<br/><br/>2<br/><br/>Pour un compresseur : wi &gt; 0. Le travail de transvasement correspond au travail que l'on doit fournir au compresseur réversible pour faire passer 1 kg de fluide de la pression p1 à la pression p2.<br/><br/>Interprétation graphique dans le diagramme (T, s) : La transformation est adiabatique réversible, donc isentropique. On définit le point 1' tel que h1' = h1, c'est-à-dire T1' = T1 pour un gaz parfait. Le travail indiqué massique wi est noté wis : travail indiqué massique isentropique. Sur l'isobare p2 : dh = Tds + v dp , d'où h2 is - h1' =  Tds = aire A<br/>1' 2 is<br/><br/>T p p2 2is 1' p1<br/><br/>1<br/><br/>s aire A = h = wis<br/><br/>On a donc : w is = h = h2 is - h1 = w ts = h2is - h1' = aire A<br/>3.5.2. Présentation du compresseur réel adiabatique Dans un compresseur réel, il faut tenir compte des forces de viscosité (frottements fluides internes). On a donc une irréversibilité mécanique.<br/>p T<br/>isobare p2 2 2is isobare p1<br/><br/>1 a) On considère une compression adiabatique réversible (donc isentropique) de la même pression initiale p1 à la s même pression finale p2. L'état final est noté 2is. On a une augmentation de la température lors de la compression, donc T2is &gt; T1. On se déplace sur une verticale dans le diagramme (T, s).<br/><br/>Nouvelle approche de l'étude thermodynamique des systèmes ouverts en PT.<br/><br/>Page 14<br/><br/><br/>b) La transformation réelle est irréversible : on la représente en traits pointillés. Le fluide passe de la pression p1 à la pression p2. A cause des forces de viscosité (frottements fluides internes), il y a dissipation d'énergie et donc échauffement du fluide. La température finale T2 (pour la transformation réelle) est plus grande que la température T2is (pour la transformation isentropique) : T2 &gt; T2is. On envisage 2 chemins pour aller de l'état 1 à l'état 2 : · transformation réelle adiabatique irréversible. · transformation FICTIVE réversible 1  2 qui a le même état initial 1 (T1, p1) et le même état final 2 (T2, p2). Le gaz reçoit de l'extérieur qef &gt; 0 : transfert thermique massique reçu pour la l'échauffement du fluide. transformation FICTIVE réversible correspondant à<br/><br/>INTERPRETATION PHYSIQUE : tout se passe comme si le fluide recevait un transfert thermique massique qf = qef correspondant à l'échauffement du fluide à cause des frottements fluides internes (forces de viscosité).<br/><br/>h = w i + qe  Transformation réelle :  (qe = 0 et se = 0 car adiabatique) s = se + sc  dh + dec + dep =  w if +  qef  dh = Tds + vdp  Transformation FICTIVE, réversible :   qef ds = T  2 w tf = vdp 1  h = w if + qef  donc   qef &gt;0 ds =  T Bilan : On a le même h et s pour la transformation réelle et la transformation FICTIVE réversible puisqu'on a des variations de fonctions d'état avec le même état initial (T1, p1) et le même état final (T2, p2). h = w i = w if + qef = w tf + qef  h = w i &gt; 0 , d'où  2 qef = 1 T ds  h = h2 - h1 peut être calculé avec le modèle du gaz parfait ou avec des tables,<br/>diagrammes thermodynamiques.<br/>On va voir 3 méthodes pour calculer qf .<br/><br/>Page 15<br/><br/>Jean-Noël BEURY<br/><br/><br/>3.5.3. Calcul de qf avec une transformation FICTIVE réversible, polytropique, gaz parfait. Dans de nombreux problèmes de concours PT (1998, 2002...), on modélise la transformation réelle irréversible par une transformation FICTIVE réversible, polytropique avec un gaz parfait. On connaît l'état initial (T1, p1) et l'état final (T2, p2). On peut mesurer facilement les températures et les pressions avec des capteurs. Pour calculer k, il suffit<br/><br/>d'écrire la relation Tp<br/><br/>1- k k<br/><br/>= cte pour l'état 1 et pour l'état 2.<br/><br/>p  ln  2   p1  On peut démontrer (voir Annexe 3) que k = T  p  ln  1  + ln  2   T2   p1 <br/><br/>On a : k &gt;  pour une compression et k &lt;  pour une<br/><br/>détente. Démonstration dans le cas particulier d'une compression adiabatique réversible avec un gaz parfait pour la transformation 1  2is<br/><br/>p p2 2is 2<br/><br/><br/><br/>Graphiquement, on représente dans le diagramme (p, v) la transformation réversible 1  2is et la<br/><br/>k<br/><br/>transformation polytropique (d'indice k) réversible p1 1 1  2 (qui modélise la transformation réelle). On a vu v que T2 &gt; T2is. Comme p2 = p2is, on a donc v2 &gt; v2is (car pv = rT pour un gaz parfait). La pente de la transformation 1  2 est donc plus grande que la pente de la transformation 1  2is et donc k &gt;  .<br/>L'indice f peut être remplacé par k pour rappeler que l'on travaille sur une transformation polytropique. a) Calcul du travail de transvasement polytropique : La transformation polytropique est caractérisée par :<br/><br/>pv k = cte = p1v1k = p2v 2k , donc v = Ap k avec A une constante.<br/>1  - k +1  1 1 2 2 p  = Ak  p - k +1 - p - k +1  . On peut rentrer A dans la w tk =  vdp =  Ap dp = - A 1 2  1 1 1  - +1  k - 1  k  1   k parenthèse, on obtient alors : w tk = ( p2v 2 - p1v1 ) . k -1 kr Finalement, on a : w tk = (T2 - T1 ) . k -1 b) Calcul de wi : w i = h = c p (T2 - T1 ) 1 - k 2<br/><br/>-<br/><br/>1<br/><br/>c) On en déduit : qef = w i - w tk<br/><br/>Nouvelle approche de l'étude thermodynamique des systèmes ouverts en PT.<br/><br/>Page 16<br/><br/><br/>D'où qf = qef = w i - w tk =<br/><br/>r kr (T2 - T1 ) - (T2 - T1 ) (1ère méthode de calcul de qf) k -1  -1<br/><br/>On a vu dans le paragraphe 3.5.2 que qf &gt; 0 ainsi que son interprétation physique. La relation entre le travail indiqué et le travail de transvasement figure explicitement au programme de PT avec son interprétation physique : w i &gt; w tk (cf référence [9]).<br/><br/>Remarque : le calcul de qf = qef peut se faire également en utilisant la relation Tds = cdT (voir problème de la banque PT 1998 Première partie cf référence [7]). On en déduit immédiatement : qf = qef =  Tds = c (T2 - T1 ) (2ème méthode de calcul de qf)<br/>1 2<br/><br/>Cette relation ne figure évidemment pas au programme PT mais peut être abordée dans un problème de concours.<br/>L'intérêt de la polytropique vient justement de la facilité à calculer qf.<br/><br/>Remarques concernant la polytropique pour les gaz parfaits : Selon les auteurs, la polytropique est une transformation adiabatique irréversible, tandis que pour d'autres elle est par définition réversible. On peut montrer (cf référence [5]) que ces deux approches sont rigoureusement équivalentes pour les gaz parfaits.<br/>Prenons l'exemple de la question 5 du problème de la banque PT 1998 (cf référence [7]) : on travaille avec une transformation polytropique adiabatique irréversible d'exposant k. L'énoncé demande de calculer la quantité de chaleur massique due à l'irréversibilité : qf. · La méthode A donne : dh = Tds + vdp =  qf +  w t . On en déduit : qf = h - w tk Il suffit là encore d'appliquer la méthode B : on travaille sur une transformation fictive polytropique réversible. On obtient immédiatement : h = qek + w ik = qek + w tk car dh = Tds + vdp =  qef +  w t On obtient le même résultat : qf = h - w tk .<br/>3.5.4. Évaluation graphique des transferts thermiques et de la création d'entropie T Dans certaines applications, on ne modélise pas p explicitement par une polytropique mais on s'intéresse à 2p une représentation graphique dans le diagramme (T, s). 2 p1 On définit le point 1' tel que h1' = h1, c'est-à-dire T1' = T1 pour un gaz parfait. 1' On a vu dans le paragraphe 3.5.2 que : 1<br/><br/>·<br/><br/>On a : qef =  Tds = aire A2 et w if = w tf .<br/>1<br/><br/>2<br/><br/>h = w i = w tf + qef<br/>aire A1 aire A 2<br/><br/>s<br/><br/>Page 17<br/><br/>Jean-Noël BEURY<br/><br/><br/>Or w i = h = h2 - h1 = h2 - h1' = aire sous la courbe 1'  2 = aire A1 + A2 donc A1 + A2 = w tf + A2 , d'où w tf = aire A1<br/>Conformément au programme PT (cf référence [9]), on peut linéariser la partie T +T 12 : qf = qef = aire A2 = aire du trapèze = ( s2 - s1 ) 1 2 (3ème méthode de calcul de qf) 2 2qf L'entropie générée par les frottements vaut donc : s = sc = s2 - s1 = T1 + T2 3.5.5. Comment justifier graphiquement que wi réel &gt; wis ? On représente très souvent sur les diagrammes la compression isentropique 12is. On a vu que : · w is = h2 is - h1 = h2 is - h1' = aire hachurée A3<br/><br/>·<br/><br/>w i réel = h2 - h1 = h2 - h1' = aire sous la courbe 1'  2 = aire A3 + A4<br/>T p 2 2is p2 p1<br/><br/>1'<br/><br/>1<br/><br/>s aire A3 aire A 4<br/><br/>w is = aire hachurée A3 On peut donc conclure que : w i réel &gt; w is :  w i réel = aire A3 + A4<br/><br/>Interprétation physique : Il faut fournir plus de travail au compresseur fonctionnant entre deux pressions fixées avec une transformation irréversible qu'avec une transformation isentropique : le travail de transvasement isentropique apparaît donc comme le travail minimal à fournir à un compresseur fonctionnant entre deux pressions fixées.<br/><br/>Remarque : On pourrait envisager une compression isotherme (en refroidissant au cours de la transformation) mais cette compression n'existe pas en pratique. On peut être amené à refroidir quand on passe d'un étage de compression à un autre étage de compression. Lorsque cette compression isotherme intervient dans les exercices, on la considère comme réversible (liquéfaction de l'azote par le procédé Linde-Hampson par exemple).<br/><br/>Nouvelle approche de l'étude thermodynamique des systèmes ouverts en PT.<br/><br/>Page 18<br/><br/><br/>3.6. Application de la méthode B à l'échangeur thermique possédant des irréversibilités thermiques. On reprend l'échangeur à contre courant étudié dans le paragraphe 1.2.1. La méthode B est appliquée simplement pour montrer qu'elle donne une interprétation physique aisée alors que l'interprétation physique avec la méthode A était délicate. Le fluide dans la conduite 12 passe de l'état 1 (p1, T1) à l'état 2 (p2, T2). On néglige la baisse de pression (p1 = p2) et on étudie le cas particulier où il se réchauffe. L'échange thermique étant irréversible, la température T'(x) (température à l'abscisse x de la conduite 1'2') est nécessairement supérieure à T(x) (température à l'abscisse x de la conduite 12). On considère un élément de longueur dx de la conduite 12 à l'abscisse x :<br/>T2' T1 Dm, cp O Dm', cp' T1' T2 L<br/><br/>x<br/><br/>x+dx<br/><br/>x<br/><br/>Pour la transformation REELLE irréversible, on a : dh + dec + dep =  w i +  qe    qe +  sc (  w i = 0 car pas de partie mobile de la machine) ds =  se +  sc = T'  dh = Tds + v dp  2 dh 2 v dp T2 dh - = On en déduit la variation d'entropie massique : s =  1T 1 T1 T T T2 dh L'entropie échangée massique vaut se =  T1 T ' T2 dh T2 dh T2 1 1 L'entropie créée massique vaut : sc = s - se =  - =  dh  -  T1 T T1 T ' T1 T T '  Comme T' &gt; T et dh &gt; 0 , on a donc sc &gt; 0. Remarque : On obtient le même résultat avec un refroidissement du fluide dans la conduite 12 avec T' &lt; T et dh &lt; 0 . L'application de la méthode B donne donc des résultats tout à fait cohérents alors que la méthode A aboutissait à qf = 0. Remarque : Si on a une irréversibilité thermique et une irréversibilité mécanique, on a  T2  1 1    2 -vdp  alors sc = s - se =   dh  -   +   . T T '   1 T   T1 Interprétation physique : L'entropie créée est égale à la somme de deux termes : T2 1 1 · Le premier terme est positif  dh  -  . Il traduit l'irréversibilité thermique T1 T T '  (diffusion de chaleur).<br/><br/>Page 19<br/><br/>Jean-Noël BEURY<br/><br/><br/>Le deuxième terme est positif : l'irréversibilité mécanique (due à la viscosité du fluide) se traduit par une baisse de pression. Conclusion : en pratique, la mise en équation des échangeurs se fait en étudiant simultanément les évolutions de l'état des deux fluides (voir annexe 2). L'objet de ce paragraphe est simplement de montrer que les problèmes d'interprétation rencontrés avec la méthode A n'existent plus avec la méthode B. On retrouve bien les deux sources d'irréversibilité dues aux phénomènes diffusifs (diffusion de chaleur et diffusion de quantité de mouvement).<br/>3.7. Bilan de la méthode B · La méthode B peut être comprise par les étudiants en PT puisque la démonstration (cf annexe 1) du premier principe dans le cas d'un système ouvert en régime permanent d'écoulement, à une entrée et une sortie est au programme PT (cf référence [9]). · Le seul résultat à admettre pour les étudiants en PT (voir démonstration paragraphe 3.4) est que pour une transformation adiabatique : qf = qef .<br/><br/>·<br/><br/>·<br/><br/>·<br/><br/>puisqu'il suffit de T calculer l'entropie créée, l'étude concernant exclusivement des transformations adiabatiques. D'ailleurs, on trouve la phrase suivante dans les programmes de PT 2004 (cf référence [9]): « La représentation de l'entropie générée par les frottements n'est pas exigible ». Tous les problèmes de concours PT peuvent se traiter sans aucune difficulté en utilisant la méthode B.<br/><br/>Avec la méthode B, on n'a plus besoin de définir le terme<br/><br/> qf<br/><br/>CONCLUSION Nous avons vu dans cet article que l'utilisation en PT de la relation q q ds = e + f peut entraîner de nombreuses erreurs d'interprétation et pose de T T nombreux problèmes pédagogiques :  qf . · On ne peut pas identifier l'entropie créée à T q q · L'entropie créée d'origine thermique est masquée dans la relation ds = e + f . T T · Les présentations proposées aux étudiants pour introduire cette relation suscitent beaucoup de questions : pourquoi  qf n'est pas rajouté dans le premier principe puisqu'il vient d'une source fictive extérieure et pourquoi  qe est le même dans la<br/><br/>·<br/><br/>transformation fictive de la méthode A ? q q Pourquoi utiliser ds = e + f pour étudier les systèmes ouverts en PT et T T ds =  se +  sc pour étudier les systèmes fermés (machines dithermes par exemple) étudiés en PTSI ?<br/><br/>Nouvelle approche de l'étude thermodynamique des systèmes ouverts en PT.<br/><br/>Page 20<br/><br/><br/>·<br/><br/>La méthode A n'est pas compatible avec le programme PT 2004 qui insiste sur  qf . l'entropie créée que l'on ne peut pas identifier à T<br/><br/>La méthode A correspond à un tour de passe-passe et risque d'entraîner un problème de compréhension globale de la thermodynamique. Nous avons donc proposé une méthode B qui se recentre sur les fondamentaux de la thermodynamique : premier principe, second principe, identité thermodynamique, entropie échangée et entropie créée. Tous les systèmes rencontrés en thermodynamique industrielle en PT peuvent s'étudier avec ces équations fondamentales. En pratique, on utilise qf = qef uniquement pour étudier l'irréversibilité (due aux forces de viscosité) d'un compresseur et d'une turbine adiabatiques. Cette nouvelle approche de l'étude des systèmes ouverts en PT ne remet pas en cause la réforme de la thermodynamique effectuée en 1995 qui insiste sur l'identité thermodynamique, l'entropie échangée et l'entropie créée. Au contraire, elle va dans le même sens puisqu'elle remet l'entropie échangée et l'entropie créée au coeur de l'étude des systèmes ouverts comme pour les systèmes fermés depuis 1995. La méthode B permet de répondre aux problèmes de concours PT avec une approche rigoureuse. Je termine la rédaction de la partie thermodynamique industrielle dans le nouveau livre « Physique ­ Tout-en-un ­ 2ème année PT Cours et exercices corrigés ­ Collection J'intègre ­ DUNOD (parution prévue juin 2008) ». Cet ouvrage couvre, en un seul volume, la totalité des programmes de physique de 2ème année de la filière PT. Les différents systèmes rencontrés en thermodynamique industrielle sont traités en détail avec la méthode B : différents types de compresseurs, turbine, tuyère, détente de Joule-Thomson ou laminage, systèmes ouverts à plusieurs entrées et plusieurs sorties, échangeur thermique, soutirage... avec de nombreux exercices, problèmes de concours corrigés et applications industrielles (cogénération, centrale nucléaire...), contraintes environnementales liées au réchauffement climatique, exercices utilisant le logiciel Thermoptim. Dans cet ouvrage, qf = qef est appliqué uniquement pour étudier l'irréversibilité du compresseur et de la turbine. L'étude des bilans énergétiques est très souvent complétée par des bilans exergétiques (voir référence [5]) qui sont hors programme en PT. Une analyse plus détaillée de l'entropie se fait dans certaines spécialités en écoles d'ingénieurs. La création interne d'entropie est la conséquence de toutes les dissipations liées aux gradients de pression, de concentration, de potentiel chimique, de température, de vitesses... Il faut alors comparer les échelles de temps de ces phénomènes à celles de l'écoulement. Il y a transformation locale d'énergie « noble » en énergie interne. Les spécialistes essayent de calculer cette source d'entropie en traitant finement les phénomènes dissipatifs avec des logiciels de mécanique des fluides utilisant des maillages très fins et des équations microscopiques.<br/><br/>Page 21<br/><br/>Jean-Noël BEURY<br/><br/><br/>REMERCIEMENTS Tous mes remerciements chaleureux à Stéphane OLIVIER (professeur de physique en PC* au lycée Louis le Grand à Paris), Renaud GICQUEL (professeur de thermodynamique à l'école des Mines de Paris et auteur de l'excellent logiciel Thermoptim qui est utilisé dans les classes de PT, PT* en cours, en séance d'informatique et en TIPE), Gilles BERNARD (professeur de thermodynamique à l'ENSAM de Cluny), Michael HEYER (professeur de thermodynamique à l'ENSAM de Metz) ainsi qu'à mes collègues du lycée Gustave Eiffel de Bordeaux : Gérard DAREAU, Jean-Pierre JORRE (calculs de l'annexe 2), Georges ROUSSELET, Pierre SALLES, pour les nombreux échanges constructifs et la relecture attentive de cet article.<br/><br/>BIBLIOGRAPHIE [1] OUZIAUX, Cours de thermodynamique. [2] DESCAMPS Didier, ENSAM Lille http://didierdescamps.free.fr/thermo.pdf [3] OUZIAUX Roger, « Mécanique des fluides appliquées », Dunod, 2004. [4] KIRILLIN, SYTCHEV, SHEINDLIN « Thermodynamique technique, technique soviétique », éditions Mir Moscou, 1976. [5] GICQUEL Renaud, « Systèmes énergétiques, tome 1 et 2 », Ecole des Mines de Paris Les Presses, 2001. [6] OLIVIER Stéphane et GIE Hubert, « Thermodynamique 1ère et 2ème année », éditions Tec et Doc collection de sciences physiques, 1996. [7] Problème de concours - Banque PT 1998 ­ Epreuve IIA. [8] Problème de concours - Banque PT 2002 ­ Epreuve IIA. [9] Programme des Classes Préparatoires aux Grandes Ecoles, PT ­ programme 2004 : B.O. hors série n°6 du 16 septembre 2004 ­ Annexe 6 page 1734.<br/><br/>Jean-Noël BEURY Professeur en Sup PCSI-SI Lycée Gustave Eiffel 33000 BORDEAUX<br/><br/>Nouvelle approche de l'étude thermodynamique des systèmes ouverts en PT.<br/><br/>Page 22<br/><br/><br/>ANNEXE 1 : EXPRESSION DU PREMIER PRINCIPE ET DU SECOND PRINCIPE POUR UN SYSTÈME OUVERT EN RÉGIME PERMANENT D'ÉCOULEMENT, À UNE ENTRÉE ET UNE SORTIE<br/>On considère un fluide qui s'écoule de l'état 1 (p1, v1, T1) à l'état 2 (p2, v2, T2). La démonstration suivante est très générale et suppose que l'on est en régime permanent d'écoulement. On note en minuscules les grandeurs massiques.<br/>a) Comment se ramener à un système fermé ? · Si on considère le système A'B'CD (appelé partie commune PC ou volume de contrôle) avec une frontière fixe, on a un système ouvert (échange de matière). On ne peut pas appliquer le premier principe de la thermodynamique. · On se ramène à un système fermé (  ) de la façon suivante :<br/><br/>- système fermé (  ) à t : partie commune (PC) + masse dm1 (de volume dV1) qui rentre pendant dt. - système fermé (  ) à t+dt : partie commune (PC) + masse dm2 (de volume dV2) qui sort pendant dt.<br/>masse dm2 qui sort pendant dt Système fermé à t C partie commune (PC) D A P1, v1, T1 B B' A' D' Système fermé à t + dt C' P2, v2, T2<br/><br/>masse dm1 qui rentre pendant dt<br/><br/>b) Régime permanent d'écoulement On dit aussi que le régime est stationnaire ou indépendant du temps. Toutes les fonctions d'état du système ouvert sont les mêmes à t et à t+dt : UPC ( t ) = UPC ( t + dt ) ; Em PC ( t ) = Em PC ( t + dt ) ; SPC ( t ) = SPC ( t + dt ) et<br/>MPC ( t ) = MPC ( t + dt )<br/><br/>c) Conservation de la masse Le masse du système fermé est la même à t et à t+dt (définition d'un système fermé).  M ( t ) = MPC ( t ) + dm1  M ( t + dt ) = MPC ( t + dt ) + dm2 <br/><br/>Comme M  ( t ) = M ( t + dt ) , on en déduit immédiatement que dm1 = dm2 = dm et<br/>dm1 dm2 = . On a conservation du débit massique : Dm1 = Dm 2 dt dt<br/><br/>Page 23<br/><br/>Jean-Noël BEURY<br/><br/><br/>d) Expression du 1er principe Le premier principe pour le système fermé (  ) s'écrit :<br/><br/>dU + dEm =  Wextérieur +  Qextérieur · · · ·<br/><br/>dU = U ( t + dt ) - U ( t ) = {UPC ( t + dt ) + dm u2 } - {UPC ( t ) + dm u1} = dm ( u2 - u1 )<br/><br/>= dm ( ep 2 - ep1 + ec 2 - ec1 )<br/><br/>dEm = Em ( t + dt ) - Em ( t ) = {Em PC ( t + dt ) + dm ec 2 + dm ep 2 } - {Em PC ( t ) + dm ec1 + dm ep1}<br/><br/> Qextérieur = qe dm = transfert thermique reçu de l'extérieur (échangé avec l'extérieur)<br/>à travers toute la frontière du système ouvert.  Wext =  Wpression amont +  Wpression aval +  Wi = travail des forces extérieures ne dérivant<br/><br/>pas d'une énergie potentielle (compris dans dEm) W  pression amont = F1.dl 1 = + p1S1dl1 = + p1dV1 = p1v1dm &gt; 0 car force motrice (le volume de la masse qui rentre est exprimé avec le volume massique v1) On peut imaginer qu'un piston pousse du côté 1 pour faire rentrer la masse dm1. Le travail de cette force est donc positif. W  pression aval = F2 .dl 2 = - p2S2dl 2 = - p2dV2 = - p2v 2dm &lt; 0 car force résistante (le volume de la masse qui sort est exprimé avec le volume massique v2).  Wi = w i dm = travail indiqué = travail reçu par le fluide de la part des parties<br/>mobiles de la machine.<br/><br/>Finalement, on obtient : ( u2 - u1 ) + ( ep 2 - ep1 + ec 2 - ec1 ) = qe + p1v1 - p2v 2 + w i 12 12 On pose h2 = u2 + p2v 2 et h1 = u1 + p1v1 et on a : h2 - h1 + ep 2 - ep1 + c2 - c1 = w i + qe 2 2 En grandeurs massiques, le premier principe pour un système ouvert en régime permanent s'écrit : h + ep + ec = w i + qe Si on l'écrit en notation infinitésimale (états 1 et 2 très proches) : dh + dep + dec =  w i +  qe<br/><br/>Tout se passe comme si une masse de 1 kg était transvasée de l'entrée (ABA'B') à la sortie (CDC'D') : on s'est donc ramené à un système fermé. On dit qu'on a une machine avec TRANSVASEMENT puisqu'un fluide s'écoule en passant de l'état 1 à l'état 2. Tout se passe comme si une masse dm était « transvasée » de l'état 1 à l'état 2.<br/><br/>wi est appelé le travail indiqué massique. Pour un compresseur : wi &gt; 0 et pour une turbine : wi &lt; 0. qe est le transfert thermique massique reçu de l'extérieur à travers toute la frontière<br/>du système ouvert.<br/><br/>Nouvelle approche de l'étude thermodynamique des systèmes ouverts en PT.<br/><br/>Page 24<br/><br/><br/>Remarque : La démonstration précédente a été faite en supposant le régime permanent. On peut remplacer cette hypothèse par l'hypothèse d'une machine cyclique : c'est une machine qui revient dans le même état au bout du temps dt. Le dt choisi est le commun multiple de toutes les périodes des différents composants du système.<br/>e) Expression du second principe Le deuxième principe pour le système fermé (  ) s'écrit : dS =  Se +  Sc<br/><br/>avec dS = S ( t + dt ) - S ( t ) = {SPC ( t + dt ) + dm s2 } - {SPC ( t ) + dm s1} = dm ( s2 - s1 ) = dm s<br/><br/>On pose  Se = se dm et  Sc = sc dm On a donc : s = s2 - s1 = se + sc L'entropie échangée massique se est nulle si la transformation est adiabatique. L'entropie créée massique sc est strictement positive si la transformation est irréversible. L'entropie créée massique est nulle si la transformation est réversible. Si on l'écrit en notation infinitésimale (1 et 2 très proches) : ds =  se +  sc<br/><br/>ANNEXE 2 : THERMIQUE<br/><br/>CRÉATION<br/><br/>D'ENTROPIE<br/><br/>DANS<br/><br/>UN<br/><br/>ÉCHANGEUR<br/><br/>On considère un échangeur thermique parfaitement calorifugé, de longueur L, parcouru par deux fluides à contre-courant :<br/><br/>T2' T1 Dm, cp O<br/><br/>Dm', cp'<br/><br/>T1' T2 L<br/><br/>x<br/><br/>x+dx<br/><br/>x<br/><br/>On note T1 et T1' les températures, supposées connues, des fluides à l'entrée de l'échangeur, T2 et T2' sont respectivement les températures à la sortie. On admet que les températures ne dépendent que de x et on néglige la conduction thermique le long de l'axe Ox. On note g la conductance thermique par unité de longueur de la paroi échangeuse. Le premier principe de la thermodynamique pour chacun des fluides en régime permanent s'écrit : Dm [ h( x + dx ) - h( x )] = gdx [T '( x ) - T ( x )] ; Dm ' [ h '( x ) - h '( x + dx )] = gdx [T ( x ) - T '( x )] Soit : dh dh ' = g (T '- T ) ; Dm ' = g (T '- T ) dx dx On assimile les deux fluides à des gaz parfaits : dT dT ' Dm c p = g (T '- T ) ; Dm ' c p ' = g (T '- T ) dx dx<br/>Dm<br/><br/>Page 25<br/><br/>Jean-Noël BEURY<br/><br/><br/> g g dT dT ' = T '- T (eq.1) ;  ' = T '- T (eq.2) On obtient :  dx dx · En faisant la différence des deux équations précédentes, on obtient : dT dT '  -' = 0 , on a T -  'T ' = T1 -  'T2 ' (eq.3) et dx dx T -  'T ' = T2 -  'T1 ' (eq.4)<br/>On divise (eq.1) par  et (eq.2) par  ' . En faisant la différence : d (T '- T ) dT dT ' 1 1 1 1 = -k (T '- T ) - = T '  -  - T  -  , soit dx dx dx  '  ' On en déduit : T '- T = (T2 '- T1 ) e - kx (eq.5) et T '- T = (T1 '- T2 ) e k ( L - x ) (eq.6) En combinant (eq.3) et (eq.5), on a : T '- T1 - T = T1 + 2 (1 - e-kx ) (eq.7) ; T ' = T2 '+ T2k' T1 (1 - e -kx ) (eq.8) k '<br/>·<br/>En combinant (eq.4) en x = 0 et (eq.7) en x = L, on obtient :  '-  (T '- T ) (eq.9) T2 '- T1 =  ' -  e - kL 1 1 En combinant (eq.6) en x = 0 et (eq.9), on a :  '-  (T '- T ) (eq.10) T1 '- T2 =  ' e kL -  1 1 On en déduit immédiatement (avec eq.7 et eq.9) :  ' 1 - e - kx T = T1 + (T1 '- T1 ) (eq.11) ;  ' -  e - kL  e - kL - e - kx puis avec (eq.5), (eq.11) et (eq.9) on a : T ' = T1 '+ (T1 '- T1 )  ' -  e - kL<br/><br/>On pose :  =<br/><br/>Dm c p<br/><br/>; ' =<br/><br/>Dm ' c p '<br/><br/>et k =<br/><br/>1<br/><br/>-<br/><br/>1 . '<br/><br/>(<br/><br/>)<br/><br/>(<br/><br/>)<br/><br/>Remarque : Dans le cas très particulier où  =  ' ( k = 0 ) : T '- T = T2 '- T1 = T1 '- T2 En combinant (eq.1) et (eq.3), on a T = T1 +<br/>T2 '- T1<br/><br/><br/><br/>x (eq.11) ; T ' = T2 '+<br/><br/>T2 '- T1<br/><br/><br/><br/>x<br/><br/>(eq.12) Dans ce cas, et ce cas seulement, l'écart des températures est constant et vaut : T '- T1 (en combinant eq.11 et eq.12 en x = L) T '- T = 1 1+ L /  Cet écart devient nul si L &gt;&gt;  . Un bilan entropique sur l'échangeur complet donne :  Sc dS = Dm (s2 - s1 ) + Dm '(s2 '- s1 ') = dt dt dh dT , donc : Or, pour une isobare dans le cas particulier d'un gaz parfait, ds = = cp T T<br/><br/>Nouvelle approche de l'étude thermodynamique des systèmes ouverts en PT.<br/><br/>Page 26<br/><br/><br/>T  T '  = Dm c p ln  2  + Dm ' c p ' ln  2  dt  T1   T1 '   T ' - T  ' 1 - e - kL   T '- T  e - kL - 1  SC 1 1 1  + g  ' ln  1 + 1 = g  ln  1 + - kL    dt T1  '-  e T1 '  '-  e - kL    Un développement limité à l'ordre un donne : - kL - kL  SC T1 '- T1  ' 1 - e T1 '- T1  e - 1 = g + g ' dt T1  ' -  e - kL T1 '  '- e - kL<br/><br/> SC<br/><br/>(<br/><br/>)<br/><br/>(<br/><br/>) <br/> <br/><br/>(<br/><br/>)<br/><br/>(<br/><br/>)<br/><br/> SC<br/><br/>g  ' (T1 '- T1 ) 1 - e - kL  1 1  = -  dt  ' -  e - kL  T1 T1 ' <br/>2<br/><br/>(<br/><br/>)<br/><br/> SC<br/>dt<br/><br/>=<br/><br/>g  ' (T1 '- T1 ) 1 - e - kL T1T1 '  '- e - kL<br/><br/>(<br/><br/>(<br/><br/>) = g  ' (T '- T ) ( e - e ) - e ) T T '( ' e )<br/>2 -L/  ' -L/  1 1 -L/  ' -L/  11<br/><br/>On obtient une expression symétrique en    ' . Si on suppose que  ' &gt;  on vérifie que tous les termes de cette expression sont positifs et donc que l'entropie créée est bien strictement positive. Elle est nulle si et seulement si T1 = T1 ' . On remarque que l'entropie créée n'est que du second ordre en T1 '- T1 , ce qui signifie que pour les petites différences de températures l'échange devient réversible. On peut aussi donner les formes asymptotiques : 2  SC gL (T1 '- T1 ) · Si kL &lt;&lt; 1 alors  dt T1T1 '<br/><br/> dt T1T1 ' Si l'échange n'est pas isobare il faut rajouter les termes de pression : T  T '  p   p '  SC = Dm c p ln  2  + Dm ' c p ' ln  2  + r ln  1  + r ln  1  dt  T1   T1 '   p2   p2 '  · Si kL &gt;&gt; 1 alors Les deux termes rajoutés sont positifs car p1 - p2 &gt; 0 et p1 '- p2 ' &gt; 0 (ce sont les pertes de charge dues à la viscosité des fluides).<br/>L'entropie créée apparaît donc comme la somme de l'entropie créée par l'irréversibilité thermique et de l'entropie créée par l'irréversibilité mécanique, c'està-dire par l'ensemble des phénomènes diffusifs (diffusion de chaleur et diffusion de quantité de mouvement).<br/><br/> SC<br/><br/>g  (T1 '- T1 )<br/><br/>2<br/><br/>ANNEXE 3 : COMPLEMENTS SUR LA POLYTROPIQUE<br/>Pour le compresseur et la turbine, on modélise souvent la transformation réelle irréversible par une transformation fictive réversible polytropique avec un gaz parfait.<br/>a) Introduction de la polytropique pour un gaz parfait Cela revient à considérer que les irréversibilités sont réparties uniformément tout au long de la compression, ce qui sur le plan technologique, correspond assez bien au cas des machines multi étagées (voir référence [5]). Ce modèle se traduit par la relation suivante :<br/><br/>Page 27<br/><br/>Jean-Noël BEURY<br/><br/><br/>T ds = c dT (s entropie massique, T température absolue en Kelvin) où c est la capacité thermique massique vraie de cette évolution (c est une constante, indépendante de la température). Il ne faut pas confondre pour le gaz parfait c, cv et cp. On peut écrire l'identité thermodynamique et la relation entre h et T pour un gaz parfait : dT dT dp . dh = c p dT = Tds + vdp . On en déduit que c p =c +v T T T dT dT dp (eq.1). Pour trouver une relation entre p et v, il Comme pv = rT , on a c p =c +r T T p dp dv dT faut éliminer T et dT : pv = rT  + = (eq.2). En combinant les équations 1 et p v T  dp dv  dp dp dv + = 0 , d'où ( c - c p + r ) + ( c - cp ) = 0 . Or 2, on obtient : ( c - c p )  +r v p p v p dp dv r = c p - cv , on a donc c - c p + r = c - cv . On en déduit : ( c - cv ) + ( c - cp ) =0 p v En intégrant : ( c - cv ) ln p + ( c - c p ) ln v = cte1 , d'où ln p + ln v c -cv = cte2 ; Cette évolution est donc définie par pv k = cte avec k = arrive à Tp<br/>1- k k<br/><br/>c -c p<br/><br/>c - cp c - cv<br/><br/>. En utilisant pv = rT , on<br/><br/>= cte .<br/><br/>k est appelé le coefficient polytropique de la transformation.<br/>b) Calcul de k et de c Pour calculer k, on connaît très souvent dans les problèmes l'état initial (T1, p1) et l'état final (T2, p2). La température et la pression s'obtiennent très facilement avec des capteurs de température et de pression. On utilise très souvent en thermodynamique industrielle le jeu de variables (T, p).<br/><br/>Il suffit d'écrire la relation Tp<br/>1- k k 11 1- k k<br/><br/>1- k k<br/><br/>= cte pour l'état 1 et pour l'état 2.<br/>1- k<br/><br/> T  1 - k  p2  T p  k d'où ln  1  = Tp = T2 p2 ; d'où 1 =  2  ln   k T2  p1   T2   p1   T   p  p  T  p    d'où k ln  1  = (1 - k ) ln  2  et k ln  1  + ln  2   = ln  2    T2   p1    p1   T2   p1    p  ln  2   p1  . Finalement, on obtient : k =  T1   p2  ln   + ln    T2   p1 <br/><br/>Connaissant k, cp et cv, on peut en déduire c =<br/><br/>kcv - c p k -1<br/><br/>.<br/><br/>Nouvelle approche de l'étude thermodynamique des systèmes ouverts en PT.<br/><br/>Page 28