Analyse Vrais Stirling1 PDF
Analyse Vrais Stirling1 PDF
Analyse Vrais Stirling1 PDF
Le but de cette tude est de fournir un instrument de travail mathmatique aux dveloppeurs de moteurs Stirling. Par dfinition, il rsulte de ces dveloppement des ralisations imparfaites. Le moteur contient des volumes morts incontournables et il est quip de rgnrateurs dont l'efficacit varie de 0 !. "abituellement, la littrature propose des thories supposant une absence de volumes morts et la perfection impossible du rgnrateur. Parfois ces anomalies constructives sont traites sparment. Schmidt suppose des volumes morts sans tenir compte du rgnrateur et #artini explore le rgnrateur en tenant compte des volumes morts, seulement partir du moment ou le rgnrateur atteint une efficacit de $0! la mo%enne logarithmique des tempratures ces extrmits est infinie. La mo%enne arithmtique de ces deux tempratures donne une bonne approximation, c'est donc cette approche que &'utilise dans ce document. 'ous allons maintenant anal%ser ce moteur imparfait en nous appu%ant sur les lois de la thermod%namique classique et extraire des r(gles importantes qui nous guiderons dans l'tude de notre futur moteur en contribuant ainsi fournir l'nergie de demain. Serge )lutchen*o Symboles utiliss : Prfixes et suffixes + . / P reprsente l'nergie thermique, exprime en ,oules reprsente un volume exprim en m(tres cubes. reprsente une temprature exprime en )elvin. reprsente un travail exprim en ,oule. reprsente une valeur de pression exprime en 'e0ton par m1
Les suffixes numriques sont utiliss dans les quations d'tat. 2ls font rfrence au diagramme p-v de la figure 3. Par exemple T 4 correspond la temprature la plus chaude du moteur alors que T 5 correspond la plus froide. 6ans le m7me esprit V 4 et V 5 rf(rent aux volumes correspondants. Les suffixes numriques exprims sous la forme 583, 384, 489 etc.. sont utiliss pour exprimer des transitions d'tats. Les suffixes c, r et f concernent respectivement l'espace chaud, le rgnrateur, et l'espace froid. ,'ai parfois, pour des raisons de clart, a&out ces suffixes dans les formules. Par exemple pour des expressions telles que V C3 et V C4 il faut lire -olume chaud sur le point 4 du diagramme et -olume chaud sur le point 9. L'ensemble des autres s%mboles utiliss est donn dans le tableau suivant. : ;ttention ) et *<
CV
e
Chaleur spcifique volume constant (J/Kg K) Efficacit du rgnrateur. Efficacit thermique du moteur Stirling. acteur dfini par l!quation (") #apport de chaleur spcifique.
Et
K k
k Mc k Mr k Mf k MT k MDP
m p
V Mc V MT
V Mr V MT V Mf V MT
V MT VT
V MT V D V P
(asse totale du fluide contenu dans le moteur ( Kg). )ression a&solue ( */m+). )ression mo'enne effective du c'cle complet (*/m+) Chaleur totale a,oute depuis une source e-trieure (Joules). Chaleur re,ete vers l!e-trieur par le c'cle. Constante du ga. de travail (J//g K) 0emprature du ga. de travail dans l!espace chaud. ( K) 0emprature du ga. de travail la sortie du rgnrateur. ( K) 0emprature du ga. de travail dans l!espace froid. ( K ) 0emprature du ga. de travail l!entre du rgnrateur.( K ) 0emprature du ga. de travail contenu dans le rgnrateur.( K ) 0emprature du refroidisseur. ( K ) 0emprature du rchauffeur. ( K ) 1olume total du ga. de travail contenu dans le moteur. ( 1olume mort chaud ( rchauffeur2 liaisons2 c'lindre) ( 1olume mort froid ( refroidisseur2 liaisons2 c'lindre) ( 1olume mort du rgnrateur.( 0otal des volumes morts. (
pm Qe Qs
R
m4
) )
m4 m4
m4
)
m4
m4
m4
#chauffeur
3e (oteur
3s #efroidisseur
)uissance mcanique
igure 45
P
3 3'
4 2 1 1'
V
igure 65 7iagrame pression volume
La figure 5 reprsente un schma simplifi d'un moteur stirling quip d'un rgnrateur imparfait. Son diagramme p-v est reprsent sur la figure 3. Pour une rgnration idale, la chaleur totale re&ete pendant la transformation 985 serait absorbe par un rgnrateur parfait et rendue au fluide de travail pendant la transformation 384. ;utrement dit, ce rgnrateur parfait aurait une surface infinie ou bien il serait ncessaire de disposer d'un temps de rgnration infini. Pour un rgnrateur imparfait, la temprature du fluide de travail l'entre et la sortie du rgnrateur correspond .5' et .4' respectivement. 6es changeurs externes doivent alors augmenter la temprature de .4' .4 et la diminuer de .5' vers .5. Le dveloppeur de moteur Stirling doit pouvoir valuer l'efficacit de son rgnrateur et d'anal%ser cette imperfection pour calculer le rendement final de sa ralisation. Le volume mort correspond la somme des volumes du fluide de travail contenus dans les changeurs, le rgnrateur ainsi que dans les raccordements entre ces diffrents lments. 2l est vident qu'un moteur Stirling contient des volumes morts irrductibles. =n pratique, on trouve des moteurs avec un volume mort reprsentant $>! du volume total.
Analyse du cycle de Stirling. Les volumes morts. 6isons que les volumes morts contenus dans le rchauffeur,le rgnrateur et le refroidisseur sont exprims en m(tres cubes et sont respectivement exprims par V Mc , V Mr et V Mf alors le volume mort total s'exprime par?
V MT =V McV Mr V Mf = k Mc k Mr k Mf V MT
:5<
ou
k Mc= V Mc V MT
V Mr reprsente le rapport du volume mort du rgnrateur sur le V MT V Mf volume mort total et k Mf = qui est le rapport du volume mort froid sur le V MT volume mort total.
k Mr=
@eprsentons maintenant le rapport entre le total des volumes morts et le volume V MT total de gaA contenu dans le moteur par k MT = alors le total des volumes VT morts peut 7tre exprim sous la forme?
V MT =k MT V T = k MT V MT V D V P
:3<
ou V D et V P reprsentent les volumes en m(tres cubes dplacs par le piston dplaceur et le piston de travail. Le volume mort d'un moteur stirling est souvent exprim en rfrence au volumes total de gaA dplac, cela nous donne alors?
V MT =k MDP V DV P
:4<
=n consquence le rapport du volume mort sur le volume total est celui du volume mort sur le total des volumes bala%s. 2l est exprim par?
k MT = k MDP k MT ou k MDP= 5 k MDP 5 k MT
:9<
:$<
La valeur de e est de 5 pour 500! d'efficacit, cela correspond une rgnration idale. Bne valeur de Aro correspond une absence totale de rgnrateur ou pire, une efficacit de 0!. La temprature du gaA de travail la sortie du rgnrateur peut 7tre calcule partir de l'efficacit de celui8ci. Cela nous donne alors?
T 4 ' = T 5e T 4T 5
:$a<
Pour un rgnrateur a%ant une m7me efficacit de chauffage et de refroidissement? Q 34 ' =Q 95 ' la temprature du gaA de travail l'entre sera?
T 5 ' =T 4e T 5T3 = T 4 e T 4 T 5
:D<
Pour un moteur Stirling contenant un volume mort notable, il est important d'avoir une temprature correcte du fluide de travail au niveau du rgnrateur. La temprature effective du gaA contenu dans le volume mort du rgnrateur peut 7tre dtermin par une simple mo%enne arithmtique?
T R= T 4 ' T 5 ' 3
:E<
:Ea<
'ous observons qu'en utilisant la mo%enne arithmtique, la temprature effective du rgnrateur devient indpendante de son efficacit. !uations d"tat: 'ous commenFons par l'quation d'tat correspondant la transformation isotherme 583. 'ous avons les espaces chaud et froid, respectivement V c et V f ou les tempratures du fluide de travail correspondent respectivement T 4 et T 5 . Sans oublier la temprature du rgnrateur ? T R . 6ans cette quation nous tenons compte des volumes morts V Mc , V Mf et V Mr Cela nous donne? mR mR p= = V c V Mc V Mr V Mf V f Vc Vf :>< K T4 T4 Tr T5 T5 T4 T5 ou
K= V Mc V Mr V Mf T4 Tr T5
: <
: a<
2l devient clair maintenant que pour des tempratures donnes du fluide de travail le facteur ) est fonction des volumes morts. La compression isot#erme. $ %-&' 6urant la transformation correspondant la compression isotherme, le fluide de travail contenue dans l'espace froid est comprim depuis V f1 =V D V P vers V f2 =V D . V f1 et V f2 reprsentent les volumes correspondant aux positions 5 et 3 du diagramme p-v de la figure 3. 6onc la chaleur re&ete pendant la transformation isotherme 583 est?
Q 53=W 53= pdV
V f1 V
f2
f2
G mRT 5
dV f / V f KT 5
V
f1
V f2 KT 5 V f1 KT 5 V D KT 5 G mRT 5 ln V D V P KT 5
G mRT 5 ln
:50<
2l est bon de noter que le travail de compression dpend seulement du facteur ( , fonction des volumes morts. Le c#auffage isoc#ore. $ &-)' =n principe, la chaleur a&oute durant le chauffage isochore :384< est?
Q 34= mC V T 4T 3 G mC V T 4T 5
:55<
ou C V est la chaleur spcifique volume constant : &H*g )< suppose 7tre une constante. Sans rgnrateur, cet apport de chaleur est obtenu depuis une source externe. ;vec un rgnrateur idal, cet apport proviendrait uniquement de celui8 ci. La chaleur restitue par un rgnrateur imparfait durant cette transformation est?
Q 34 ' = mC V T 4T 4' G emC V T 4 T 5
:53<
:54<
In voit bien que pendant cette phase du c%cle, la chaleur apporte dpend de l'efficacit du rgnrateur seulement. L"expansion isot#erme $ )- * ' 6urant la transformation isotherme : 489< le volume chaud contenant le gaA de travail varie de V C3 =V D V C4 =V D V P Le volume contenant le gaA froid V f est suppos 7tre Aro. Pendant cette transformation, la chaleur apporte est?
Q 49 = W 4 9 = m pDV c
V c3 V c4
dV c V V c KT 4 V C4 KT 4 G mRT 4 ln V C3 KT 4 V DV P KT 4 G mRT 4 ln V D KT 4
G mRT 4
c3
V c4
:59<
2l est maintenant vident que le travail fourni est dpendant des volumes morts. Le refroidissement isoc#ore: $*-%' La chaleur re&ete pendant la transformation isochore :985<, phase de refroidissement du fluide est?
Q 9 5 =mC V T 5T 9 G mC V T 4T 5
:5$<
Sans rgnrateur, cette chaleur est dissipe dans le refroidisseur. ;vec un rgnrateur parfait elle serait compl(tement absorbe par celui8ci. ;vec une rgnration imparfaite la chaleur absorbe par le rgnrateur est alors exprime par la formule suivante?
Q9 5 ' = mC V T 5 ' T 9 G emC V T 4T 5
:5D<
:5E<
'ous vo%ons ici que le transfert de chaleur pendant la phase de refroidissement dpend de l'efficacit du rgnrateur.
+#aleur totale a,oute. Pour une rgnration imparfaite, la chaleur totale a&oute depuis une source extrieure dans le c%cle est? Q e =Q 44 ' Q 49 :5>< -o%ons cela plus en dtails?
Qe = mC V T 4T 4 ' k 5 T 4 ln V c4 KT 4 V c3 KT 4
V D V P KT 4 :5> a< V D KT 4 2l est bon de rappeler que - reprsente le rapport de chaleur spcifique. 2l paraJt maintenant vident que la chaleur totale introduite dans le moteur dpend de la qualit du rgnrateur et de la quantit de volumes morts prsente. Sans rgnration ou avec un rgnrateur d'une efficacit nulle, la chaleur apporte au moteur depuis une source extrieure serait
G mC V 5 e T 4T 5 k 5 T 4 ln
:5 < :30<
+#aleur totale re,ete vers l"extrieur. ;vec un rgnrateur imparfait la chaleur extraite du c%cle est
Q s = Q 5 5 ' Q 5 3 cela donne?
:35<
V f1 KT 5 V f2 KT 5
Q s=mC V T 5T 5' k 5 T 5 ln
V D V P KT 5 :35a< V D KT 5 'ous avons ici une forte ressemblance entre chaleur a&oute et chaleur extraite. La dpendance de ces phnom(nes la qualit du rgnrateur et de la prsence de volumes morts devient indiscutable.
G mC V 5 e T 4T 5 k 5 T 4 ln
:33<
:34<
2l paraJt vident qu'une s%mtrie de fonctionnement du rgnrateur cot chaud et cot froid apparaJt ici aussi comme indispensable.
.ravail effectif. Le surplus d'nergie des deux transformations isothermes 583 et 489 est transform en travail mcanique utile. Ce travail, avec un moteur muni d'un rgnrateur imparfait et comprenant des volumes morts peut 7tre calcul comme ceci?
W net = Q G Q e Q s G Q 4 4 ' Q 4 9 Q 5 5 ' Q 5 3 G Q4 9 Q 5 3 V c4 KT 4 V f1 KT 5 ]T 5 ln [ ]} G mR {T 4 ln [ V c3 KT 4 V f2 KT 5
:39<
2l faut noter que V c4 : -olume chaud en position 9 sur le diagramme p8v< est V f1 :isochore< donc gal V D V P=V T . =t que gal V c3 =V c2 =V D=V 3 on en dduit?
W net =mR {T 4 ln [ V DV P KT 4 V D V P KT 5 ]T 5 ln [ ]} V D KT 4 V D KT 5
:39 a<
2l est vident que le travail fourni par le c%cle dpend des volumes morts. Si ils taient rduit nant cela nous donnerait ceci?
W net =mR ln V5 T 4T 5 V3
:3$<
Cette quation est celle que l'on trouve dans les ouvrages de thermod%namique. Constatons aussi qu'en l'absence de volumes morts le travail fourni par le c%cle est totalement indpendant de l'efficacit du rgnrateur K
/ression effective moyenne. Le travail mcanique peut 7tre dtermin par la pression mo%enne effective du c%cle : P m < et des changements de volumes qui se rsument finalement au volume bala% par le piston de travail?
V c4 V C3 = V f1 V
f2
= V 5V V2 = V p
cela donne ?
W net = P m V P
:3D<
:3E<
'otons encore que la pression effective dpend des volumes morts. Si il n'% en avait aucun, cela nous donnerait ceci,. en utilisant la loi des gaA parfaits ainsi qu'en se souvenant que T 5=T 3 et ue V 5=V 4 ?
P m=[ P 4 P 3 V 4 V5 ] ln V 5 V 3 V3
P 4 P 3 V5 ln V5 V3 G V 35 5
:3><
Bne autre solution consiste calculer la pression aux quatres L coins M du c%cle et d'en faire la mo%enne. Cette solution passe par le calcul de la densit du gaA par rapport au volume spcifique. C'est le principe adopt dans la mise en pratique expose plus loin.
0fficacit t#ermi!ue. L' efficacit thermique d'un moteur Stirling W net Et= Qe cela donne pour notre moteur imparfait?
E t peut 7tre vu sous la forme ?
:3 <
V D V P KT 4 V D V P KT 5 {T 4 ln { }T 5 ln { }} V D KT 4 V D KT 5 Et= :3 a< V DV P KT 4 5 e {T 4 ln [ ] T 4T 5 } V D KT 4 k 5 In remarque que l'efficacit thermique du moteur Stirling dpend la fois de la qualit du rgnrateur et de la proportion de volumes morts prsents.
:40<
#7me lorsque une notation diffrente est emplo%e comme dans l'quation :40<, les rsultats sont identiques en utilisant la thermod%namique classique ou la thermod%namique en temps fini. =n l'absence de rgnrateur :e 1 2<, ce qui reprsente le pire des cas pour l'efficacit du c%cle du moteur Stirling. L'efficacit thermique serait alors?
E t =T 4T 5 /[ T 4 T 4T 5 / k 5 ln V5 ] V3
:45<
In remarque que les quations :40< et :45< sont identiques en l'absence de volumes morts. Si maintenant nous disposions du rgnrateur parfait e G 5 et d'une absence totale de volumes morts, ce qui serait idal pour un moteur Stirling, l'efficacit thermique serait?
Et= T 4T 5 T4
:43<
Ce qui nous fait dire que le c%cle de Stirling poss(de l'efficacit du c%cle de Carnot oprant entre deux sources thermiques T 4 et T 5 respectives. 2l est maintenant vident qu'en thorie, le moteur Stirling peut 7tre une machine tr(s rentable pour convertir de la chaleur en travail mcanique. L'efficacit d'une conception et d'une fabrication qui exclurait les volumes morts tout en incluant un rgnrateur quasi8parfait quivaudrait raliser une machine de Carnot. ;utrement dit ce serait la ralisation du moteur thermique le plus efficace qu'il soit en termes de rendement nergtique.
3ise en prati!ue: #uni de nos belles formules, nous allons examiner les effets des volumes morts et l'effet de l'efficacit du rgnrateur sur un moteur Stirling L rel M non pressuris. 'ous allons prendre un taux de compression optimis que nous appellerons C t VD T1 pour cet exemple, les param(tres initiaux de notre moteur gal T4 seront les suivants? .emprature chaude? .emprature froide? -olume de gaA dplac par le piston de travail? -olume de gaA dplac par le dplaceur avec -olume mort morts cOt chaud -olume mort du rgnrateur -olume mort cOt froid
T 4 G 34 ) : D$0 NC<. T 5 G 44> ) : D$ NC<. V P G 3E 9D9 cc : optimis < V D G E$000cc C t = 0.4DD V Mc G 300 cc V Mr G 3$0 cc V Mf G 5$0 cc
La procdure de calcul d'un moteur Stirling imparfait peut maintenant 7tre ralise comme suit? 6e l'quation :Ea< pour une temprature du gaA de travail maintenu ) et T 5 G 44> ), la temprature effective du rgnrateur est?
T R= T 4 T 5 3
T 4 G 34
G D40,$ )
Iccupons nous maintenant des volumes morts. Le volume mort total contenu dans le moteur V MT =V Mc V Mr V Mf = 300 P 3$0 P 5$0 G D00 cc G D.0=89 m4 6e l'quation :9< le rapport des volumes morts sur les volumes bala%s est?
k MDP = V MT = 0.00$>$$ V D V P
:44<
:49<
-o%ons maintenant le rapport de chacun des volumes morts sur le volume mort total suivant l'quation :3<
K Mc = V Mc =0,444 V MT
:4$<
K Mr =
V Mr =0,95D V MT VM = 0,3$0 V MT
:4D<
K Mf =
:4E<
:4><
:4 < :90<
@emplissons notre moteur avec de l'air pression atmosphrique. Comme il n'est pas tablis de position de remplissage normalis nous choisissons de placer notre moteur en position 4. Le volume correspond alors au volume total maximum. V T = V D V PV MT . La temprature de remplissage tant fixe 30 N C : 3 4.0 )< , elle est valable pour les trois volumes. 6onc ?
m = VT p 505900 0.5040 = = 0.5393E R 3 4,0 3 4 3 4 k!
:95<
'ous allons maintenant calculer les diffrentes transitions d'tat pour un c%cle complet en tenant compte de l'efficacit du rgnrateur, fix arbitrairement E$! : e G 0,E$<, et des volumes morts prsents. 6e l'quation :50< , le travail de compression isotherme583 est?
W 53= mRT 5 ln [ V D KT 5 ] G 84E9D.9DD3 ,oules V D V P KT 5
:93<
:94<
'ous savons que le travail fournit par le moteur est la somme de ces deux derniers rsultats. Les deux autres transformation nous servant calculer le rendement nergtique de notre moteur.
W net = W 5 3W 4 9 = D953.D330 ,oules
:99<
Cependant si nous utilisons la formule :39a< pour calculer ce travail le rsultat est identique la somme voque plus haut.
W net =mR {T 4 ln [ V DV P KT 4 V D V P KT 5 ]T 5 ln [ ]} V D KT 4 V D KT 5
:9$<
Calculons le rendement nergtique des deux isochores en tenant compte de l'efficacit du rgnrateur. =st il besoin de rappeler que C " reprsente la Chaleur spcifique volume constant :,H)g )<. Sa valeur est de E5> pour l'air. 2l est bon de connaJtre maintenant les tempratures de sortie du rgnrateur 5' et 4', dpendantes de l'efficacit du rgnrateur dfinit par e. 'ous utilisons les formules :$a< et :D<.
T 5 ' = T 4e T 5T 4
8
9>9.3$ ) )
:9D< :9E<
T 4 ' = T 5 e T 4 T 5
8 EED.E$
Le chauffage isochore 384 est la sommes des chaleurs apportes par le rgnrateur et par celle apporte depuis une source extrieure. 6e l'quation :55< la chaleur totale apporte pendant cette transformation est?
Q34= mC V T 4 T 5 G $3305.03 5 ,oules
:9><
6e l'quation :54< la chaleur apporte depuis une source extrieure pendant la transformation 4'84 est?
Q44 ' = 5e mC V T 4T 5 G 540$0.3$E ,oules
:9 <
Le refroidissement isochore 985 est la somme de la chaleur re&ete dans le rgnrateur et celle extraite vers le milieu extrieur. =n se servant de l'quation :5$< cela nous donne?
Q9 5 =mC V T 5 T 4 G 8$3305.03 5 ,oules
:$0<
6e l'quation :5E< la chaleur extraite vers le milieu extrieur pendant la transformation 585' est?
Q55 ' = mC V T 5T 5 ' G 8540$0.3$E3 &oules
:$5<
E3.D D5 ,oules
Le travail net par c%cle , que nous avons d& calcul, peut aussi 7tre vue comme la somme de Q e et de Q s . :$9<
:$$<
L'efficacit de Carnot suivant la formule :44< est de 0.D44> ce qui montre que nous pouvons encore amliorer notre moteur.
Pour obtenir la pression aux quatre L coins M du moteur le plus simple est de procder ainsi. Calculons la pression sur le point 5 du diagramme P -. Le volume est gal est alors ?
V DV PV Mt et la temprature gale T5
la pression
:$D<
Le point 3 maintenant ?
5 V D V Mt H=8$ G 5.$ 9D *g H cm1 P 3 = RT 5 m
:$E<
Le point 4?
5 V D V Mt H=8$ G 9.4$9E *g H cm1 P 4 = RT 4 m
:$><
et enfin le point 9?
5 V DVP V Mt H =8$ G 4.5 93E *g H cm1 P 9 = RT 4 m
:$ <
La pression effective mo%enne s'obtient ici en faisant la mo%enne de ces 9 pressions. Cela nous donne 3.$E>4$ *g H cm1. : remplissage 5.059 *g H cm1<. :D0<
0n conclusion: 6e cette tude nous pouvons conclure ceci? 5. Pour un moteur Stirling avec un volume mort donn, un rgnrateur inefficace n'affectera pas le travail effectif si la temprature effective du rgnrateur est une mo%enne arithmtique des tempratures extr7mes du gaA de travail. Cependant un moteur Stirling muni d'un mauvais rgnrateur demandera plus de chaleur et ncessitera un meilleur refroidissement que si il tait quip d'un rgnrateur efficace. 3. Les volumes morts vont dgrader la fois l'efficacit nergtique et augmenteront les apports de chaleurs externes. 'otons que la construction d'un moteur Stirling entraJne fatalement la prsence de volumes morts incontournables. 4. Bn faible travail mcanique peut quand m7me 7tre obtenu si le moteur contient un taux de volumes morts importants. 9. Pour atteindre de hautes performances , un moteur Stirling doit 7tre quip du rgnrateur le plus efficace qu'il soit. Si ce n'tait pas le cas, il vaut mieux ne pas avoir de rgnrateur du tout. Ce qui nous fait l'conomie de son volume mort.
Serge (lutc#en-o.