ELKRAFT Slutrapport for projektet: - Energinet.dk
ELKRAFT Slutrapport for projektet: - Energinet.dk
ELKRAFT Slutrapport for projektet: - Energinet.dk
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
PSO-FU2204 1 9 kW Stirlingmotor<br />
<strong>ELKRAFT</strong><br />
<strong>Slutrapport</strong> <strong>for</strong> <strong>projektet</strong>:<br />
Komponentoptimering - Stirlingmotorer til biobrændsler<br />
PSO-projekt FU2204<br />
Revideret slutrapport, oktober 2006<br />
Henrik Carlsen<br />
Inst. For Mekanik Energi og Konstruktion<br />
Danmarks Tekniske Universitet
PSO-FU2204 2 9 kW Stirlingmotor<br />
Indhold<br />
1 Indledning .......................................................................................................................5<br />
1.1 Formål med <strong>projektet</strong> ...............................................................................................7<br />
2 Udvikling af stirlingmotor til <strong>for</strong>søg...................................................................................8<br />
2.1 Grundlæggende konstruktion af stirlingmotorerne til biobrændsler...........................8<br />
2.2 Valg af <strong>for</strong>søgsmotor..............................................................................................10<br />
2.3 Design af 9 kW <strong>for</strong>søgsmotor.................................................................................11<br />
2.4 Cylinder med topstykke, heder, regenerator og køler .............................................16<br />
2.5 Krumtapmekanisme ...............................................................................................18<br />
2.6 Krumtaphus............................................................................................................21<br />
3 Laboratorie<strong>for</strong>søg..........................................................................................................23<br />
3.1 Resultater fra laboratorieafprøvning af SM5A.........................................................26<br />
3.2 Afprøvning af andre regeneratortyper.....................................................................28<br />
3.3 Forbedring af målinger og databehandling .............................................................30<br />
3.4 Afprøvning af 9 kW motor på Fachhochschule Reutlingen .....................................32<br />
3.5 Maskine til afprøvning af regeneratorer..................................................................33<br />
4 Model til simulering af stirlingmotorer ............................................................................40<br />
4.1 Modifikation af eksisterende simuleringsmodel ......................................................41<br />
4.2 Sammenligning af resultater fra <strong>for</strong>søg med resultater fra simulering.....................47<br />
4.3 Særlige krav til stirlingmotorer til biomasse ............................................................50<br />
5 Undersøgelse af særlige komponenter og driftsbetingelser...........................................51<br />
5.1 Dellast....................................................................................................................51<br />
5.2 Ventiler...................................................................................................................52<br />
5.3 Lejetætninger og smørefedt ...................................................................................55<br />
5.4 Brint som arbejdsgas .............................................................................................57<br />
5.5 Simulering af arbejdsåg med indbyggede rulle og nålelejer....................................58<br />
6 Evaluering af teknologi og marked ................................................................................63<br />
6.1 Små fjernvarmecentraler til boligopvarmning..........................................................63<br />
6.2 Virksomheder <strong>for</strong>trinsvis i træindustrien .................................................................64<br />
6.3 Landbrug................................................................................................................64<br />
6.4 Støtteordninger i Europa ........................................................................................64<br />
6.5 Det oversøiske marked ..........................................................................................65<br />
6.6 Konkurrerende produkter .......................................................................................65<br />
7 Konklusion og <strong>for</strong>slag til videre arbejde .........................................................................66<br />
Referencer ...........................................................................................................................68<br />
Bilag 1. Stirlingmotorens virkemåde<br />
Bilag 2: Praktiske aspekter ved anvendelse af brint som arbejdsgas.<br />
Bilag 3. Tegningningsliste og enkelte detailtegninger<br />
Bilag 4: 9 kW stirlingmotor SM5A<br />
Bilag 5: Resultater fra <strong>for</strong>søg med SM5A-002 stirlingmotoren<br />
Bilag 6: Stirling engines and biomass –what is the problem?<br />
Bilag 7: ANSYS Finite element raport vedrørende arbejdsåg.
PSO-FU2204 3 9 kW Stirlingmotor<br />
Abstract (english)<br />
Based on the experience from a <strong>for</strong>mer development project, a new single cylinder,<br />
hermetically sealed Stirling engine with an electric power output of 9 kW has been designed<br />
and manufactured. The new engine is a complete redesign of the <strong>for</strong>mer design and many<br />
improvements have been implemented.<br />
One of the engines from the <strong>for</strong>mer project has been utilized <strong>for</strong> a comprehensive laboratory<br />
test of engine per<strong>for</strong>mance and new components and materials. Results from tests have<br />
been compared with results from an improved Numerical Simulation Program (NSP) <strong>for</strong><br />
prediction of Stirling engine per<strong>for</strong>mance which showed that the NSP was able to reproduce<br />
the results from laboratory with a satisfactory accuracy.<br />
A novel machine <strong>for</strong> testing of heat transfer in regenerators has been designed,<br />
manufactured and utilized <strong>for</strong> testing of different regenerator designs.<br />
Although the actual project differs from the original plan, most of the tasks <strong>for</strong>eseen have<br />
been made, and the results have already contributed to the development of the Stirling<br />
engines.
PSO-FU2204 4 9 kW Stirlingmotor<br />
Resumé<br />
Projektets <strong>for</strong>mål var at fremstille og instrumentere en 9 kW stirlingmotor, der alene skulle<br />
anvendes som <strong>for</strong>søgsmotor til detaljerede målinger i laboratoriet samt afprøve nye<br />
komponenter og materialer med henblik på at <strong>for</strong>bedre stirlingmotorernes effekt og<br />
virkningsgrad samt levetid. Formålet var yderligere at <strong>for</strong>bedre en simuleringsmodel, der var<br />
udviklet i <strong>for</strong>bindelse med tidligere gennemførte projekter.<br />
Forsøgsmotoren skulle bygge på erfaringerne fra et tidligere projekt, hvor en én-cylindret<br />
hermetisk stirlingmotor med en el-effekt på 9 kW var udviklet. Desværre viste det sig, at der<br />
var mange problemer med denne motors konstruktion, og der blev brugt betydelige<br />
ressourcer og en del tid på at få de eksisterende motorer til at fungere tilfredsstillende. Det<br />
medførte, at revisionen af den én-cylindrede motors design blev mere omfattende end<br />
beregnet. Samtidig medførte andre problemer, at design og fremstillingen af komponenter til<br />
den nye motor blev væsentlig <strong>for</strong>sinket, og den nye motor var der<strong>for</strong> først samlet og klar til<br />
afprøvning ved <strong>projektet</strong>s afslutning.<br />
En væsentlig del af <strong>projektet</strong>s indhold var at afprøve <strong>for</strong>søgsmotoren i laboratoriet og<br />
afprøve <strong>for</strong>skellige nye <strong>for</strong>slag til komponenter og delløsninger. En af de to første 9 kW<br />
motorer er der<strong>for</strong> benyttet til at gennemføre mere end 75 <strong>for</strong>søg i <strong>for</strong>skellige driftstilstande,<br />
samt til at få afprøvet <strong>for</strong>skellige delkomponenter. Motoren er blandt andet benyttet til<br />
afprøvning af ventiler samt nye typer smørefedt. En afprøvning med brint som arbejdsgas er<br />
af sikkerhedsmæssige grunde erstattet af en afprøvning med nitrogen.<br />
En ny type fiberregenerator er afprøvet i <strong>for</strong>søgsmotoren og resultaterne sammenlignet med<br />
den oprindelige regenerator fremstillet på DTU. Resultaterne viste, at regeneratorerne ikke<br />
var bedre end den eksisterende type, der er fremstillet af tynd tråd, men at de nye<br />
fiberregeneratorer meget vel kan vise sig at være billigere at fremstille.<br />
Det er vanskeligt at måle varmeovergangstallet <strong>for</strong> regeneratorerne i en stirlingmotor, <strong>for</strong>di<br />
det ikke kan gøres under stationære <strong>for</strong>hold. Der<strong>for</strong> er der udviklet en særlig maskine til at<br />
afprøve 9 kW motorens regeneratorer under periodisk strømning. Dette er ikke set før, og<br />
trods problemer med den elektriske varmeveksler er resultaterne meget interessante. Der er<br />
desuden gennemført stationære målinger af trykfald, og resultaterne er benyttet til at<br />
verificere de empiriske algoritmer i et numerisk simuleringsprogram til beregning af<br />
stirlingmotorers egenskaber. Resultater fra dette program er efterfølgende sammenlignet<br />
med resultater fra <strong>for</strong>søg, som viste at der er en god overensstemmelse mellem beregnede<br />
og målte resultater.<br />
Selv om <strong>for</strong>løbet af <strong>projektet</strong> afviger noget fra den oprindelige plan, er de fleste af <strong>projektet</strong>s<br />
delmål opnået, og resultaterne har allerede bidraget væsentlig til udviklingen af<br />
stirlingmotorerne.
PSO-FU2204 5 9 kW Stirlingmotor<br />
1 Indledning<br />
Distribueret el-produktion ved hjælp af små decentrale anlæg bliver med stor sandsynlighed<br />
en vigtig del af fremtidens el-<strong>for</strong>syning. Man <strong>for</strong>estiller sig, at mange små anlæg i fremtiden<br />
samles i netværker og danner et virtuelt kraftværk, der styres centralt. Markedet <strong>for</strong> små elproducerende<br />
anlæg ventes der<strong>for</strong> at vokse kraftigt i de kommende år.<br />
Kombineres el- og varmeproduktion, således at spildevarmen fra el-produktionen (kraftproduktionen)<br />
udnyttes til opvarmnings<strong>for</strong>mål, er der store miljømæssige og økonomiske<br />
<strong>for</strong>dele at hente. Disse kraft-varme-anlæg må nødvendigvis placeres i nærheden af de<br />
steder, hvor der er behov <strong>for</strong> den producerede varme. Det medfører, at der først og fremmest<br />
er behov <strong>for</strong> små kraftvarmeanlæg til biomasse, <strong>for</strong>di biomassen typisk anvendes i små<br />
varmeanlæg med en varmeeffekt på 200 – 5000 kW.<br />
Som en følge af behovet af små kraftvarmeanlæg til biobrændsler har man i Tyskland og<br />
Østrig vedtaget særdeles gunstige regler <strong>for</strong> sådanne anlæg. Således afregnes el produceret<br />
ved hjælp af biomasse på små anlæg til 0,215 EURO i Tyskland og 0,16 EURO i Østrig.<br />
El-producerende anlæg, der anvender biomasse som brændsel, er typisk baseret på<br />
dampturbineteknologi med en el-effekt på mellem 2 MW og 50 MW. Som følge af den høje<br />
politiske prioritering af små kraftvarmeanlæg til biomasse i Tyskland og Østrig er de første<br />
små anlæg med en el-effekt under 1 MW ved at komme på markedet. Disse anlæg er<br />
baseret på <strong>for</strong>gasning af biomasse og anvendelse af <strong>for</strong>gasningsgassen i en konventionel<br />
<strong>for</strong>brændingsmotor. Anlæggene, der <strong>for</strong>ventes at blive kommercielt konkurrencedygtige, har<br />
en el-effekt på mindst 100 kW, og mange mener, at el-effekten skal være højere, før det er<br />
muligt at opnå en tilfredsstillende anskaffelsespris. Kraftvarmeanlæg til biobrændsler med en<br />
el-effekt under 150 kW kan der<strong>for</strong> med <strong>for</strong>del baseres på en anden teknologi.<br />
I en stirlingmotor <strong>for</strong>egår <strong>for</strong>brændingen ikke inde i cylindrene som i en almindelig motor,<br />
men uden<strong>for</strong> cylindrene. Varmen fra <strong>for</strong>brændingen overføres til motoren i en varmeveksler.<br />
Det gør stirlingmotoren velegnet til faste brændsler, <strong>for</strong>di faste brændsler som flis eller andre<br />
<strong>for</strong>mer <strong>for</strong> biobrændsler <strong>for</strong>brændes på samme måde som i en almindelig varmekedel. Et lille<br />
decentralt kraftvarmeværk til fast brændsel kan dermed opbygges ved at kombinere en<br />
stirlingmotor med en fastbrændselskedel.<br />
I 1989 indledtes et samarbejde mellem det daværende Vølund R & D og Henrik Carlsen,<br />
DTU, om udvikling af store stirlingmotorer. I Energistyrelsen var man på det tidspunkt ikke<br />
interesseret i <strong>projektet</strong>, og <strong>projektet</strong> kom alene i gang, <strong>for</strong>di <strong>ELKRAFT</strong> stillede midler til<br />
rådighed <strong>for</strong> udviklingen. Arbejdet var i perioden indtil 1993 koncentreret om udvikling af en<br />
stirlingmotor med en el-effekt på 149 kW. Motoren skulle anvendes til decentral kraftvarme<br />
og anvende flis som brændsel. I 1993 blev aktiviteterne udvidet med endnu et projekt, der<br />
sigtede på udvikling, fremstilling og afprøvning af en stirlingmotor med en el-effekt på 35 kW,<br />
der ligeledes anvendte flis som brændsel. Dette projekt fik tilskud fra Energistyrelsen til at<br />
fremstille en prototype, der blev afprøvet i laboratoriet i 1996.<br />
I samarbejde med Naturgas Midt-Nord blev der yderligere igangsat et projekt med det <strong>for</strong>mål<br />
at udviklet en mindre stirlingmotor med en el-effekt på 9 kW. Hensigten var, at denne motor<br />
skulle anvendes i små kraftvarmeanlæg, der anvendte biogas eller naturgas som brændsel.<br />
De danske projekter vedrørende udvikling af stirlingmotorer til decentral kraftvarme med<br />
biobrændsler omfattede dermed ved <strong>projektet</strong> start udvikling og afprøvning af følgende to<br />
motortyper:
PSO-FU2204 6 9 kW Stirlingmotor<br />
1. Fire-cylindret stirlingmotor/generator med en el-effekt på 35 kW til direkte fyring med flis<br />
og andre biobrændsler<br />
2. Én-cylindret stirlingmotor/generator med en el-effekt på 9 kW til fyring med gas<strong>for</strong>mige<br />
biobrændsler (biogas, <strong>for</strong>gasningsgas) og naturgas<br />
Ved <strong>projektet</strong>s start i sommeren 2002 var der afprøvet 3 <strong>for</strong>skellige modeller af den 4cylindrede<br />
motor, SM3A, SM3B og SM3C. Desuden var der fremstillet to identiske motorer af<br />
den én-cylindrede motor type, SM5A. Status <strong>for</strong> udviklingen og afprøvningen af den éncylindrede<br />
motor er behandlet i de følgende kapitler. Status <strong>for</strong> udviklingen, produktion og<br />
afprøvningen af den fire-cylindrede motor medio 2005 er, at der er fremstillet 9 motorer af<br />
den 4-cylindrede type, idet to af disse motorer er bygget sammen til en 8-cylindret enhed<br />
med fælles generator. Følgende motortyper og anlæg er produceret:<br />
Antal Samlet Med flis Tilskud, DK aktiviteter<br />
Model stk. driftstimer driftstimer Driftsafprøvning<br />
SD3A 1 2400 1400 Skarp Salling, DK Danske Energistyrelse<br />
SD3B 1 1600 800 Ansager, DK Danske Energistyrelse<br />
SD3C 1 9500 9000 Oberlech, A Østrigske Energistyrelse<br />
SD3D 1+2x1 4500 4500 Lund, S og Lingenau, A Svensk og EU<br />
SD3E 3 100 100 Feldbach & Bielefeld, A, Hjortshøj, DK ABC og Energistyrelsen<br />
I alt 9 18100 15800<br />
Status <strong>for</strong> afprøvningen af den fire-cylindrede motor er dermed følgende<br />
• Afprøvningen af den første motor, type SM3A, blev afbrudt i 1999 på grund af<br />
mekaniske problemer<br />
• De øvrige motorer er i drift eller ved at blive sat i drift efter installation eller revision<br />
• En motor (SM3C), som p.t. er i drift på et hotel i Oberlech i Østrig, har kørt ca. 9.000<br />
timer med flis som brændsel.<br />
• Motorerne har tilsammen kørt i alt ca. 18.000 timer, hvoraf ca. 16.000 timer er med flis<br />
som brændsel<br />
• På de seneste versioner af den 4-cylindrede motor er effekten i praksis over den<br />
ønskede el-effekt på 35 kW og virkningsgraden er væsentligt <strong>for</strong>bedret, så varmen til<br />
kølevandet i dag udgør ca. 90 kW ved fuld last mod tidligere helt op til 110 kW.<br />
• Konstruktionen er løbende <strong>for</strong>bedret, så det er væsentlig nemmere at samle og<br />
servicere den nye motor<br />
• Driftsafprøvninger har vist, at anlæggene kan køre ubemandet i fuldautomatisk drift.<br />
• Der er både udviklet fyringssystemer til direkte fyring med flis på en rist og til <strong>for</strong>gasning<br />
af flis og efterfølgende afbrænding af den urensede gas i en brænder på stirlingmotoren.<br />
Resultaterne fra nærværende projekt har haft stor indflydelse <strong>for</strong> de <strong>for</strong>bedringer, som er<br />
gennemført ved videreudviklingen af type SM3C til de følgende modeller SM3D og SM3E.<br />
I 2004 stiftedes selskabet Stirling.<strong>dk</strong> ApS (SD), idet DTU-Innovation investerede midler i den<br />
videre udvikling og markedsføring. SD har indgået en licens- og samarbejdsaftale med det<br />
østrigske firma MAWERA, som har installeret 4 anlæg i Østrig og yderligere solgt to anlæg af<br />
to motorer hver til henholdsvis en tysk og en østrigsk kunde. Et andet resultat af dette<br />
samarbejde er et projekt vedrørende afprøvning af nogle konkrete <strong>for</strong>slag til <strong>for</strong>bedringer af<br />
den fire-cylindrede motor, der finansieres af det østrigske Austrian Bioenergy Centre.<br />
SD er ved at installere en af de tre sidst producerede motorer (type SM3E) i Hjortshøj, hvor<br />
anlægget skal producere el og varme til ca. 100 boliger. Samtidig er SD ved at revidere
PSO-FU2204 7 9 kW Stirlingmotor<br />
anlægget i Ansager, der består af en modstrøms<strong>for</strong>gasser og en af de fire-cylindrede<br />
stirlingmotorer, således at det kan sættes i drift i starten af 2006.<br />
Udviklingsprojekterne er gennemført i samarbejde mellem skiftende parter med DTU som det<br />
centrale omdrejningspunkt, og aktiviteterne har fået tilskud fra Energistyrelsen, <strong>ELKRAFT</strong>,<br />
Naturgas Midt-Nord, EU, Austrian Bioenergy Centre m.fl.<br />
Sammenlignes resultaterne fra ovennævnte projekter med aktiviteterne andre steder i<br />
verden, er der næppe tvivl om, at den danske udvikling af stirlingmotorer til små<br />
kraftvarmeværker, der anvender biomasse som brændsel, er førende trods de vanskelige<br />
<strong>for</strong>hold <strong>for</strong> finansiering og hjemmemarked i de seneste år.<br />
1.1 Formål med <strong>projektet</strong><br />
Projektets <strong>for</strong>mål var at fremstille og instrumentere en 9 kW stirlingmotor, der alene skulle<br />
anvendes som <strong>for</strong>søgsmotor til detaljerede målinger i laboratoriet samt afprøve nye<br />
komponenter og materialer med henblik på at <strong>for</strong>bedre stirlingmotorernes effekt og<br />
virkningsgrad samt levetid. Formålet var yderligere at <strong>for</strong>bedre en simuleringsmodel, der var<br />
udviklet i <strong>for</strong>bindelse med tidligere gennemførte projekter.<br />
Ønsket om at fremstille en <strong>for</strong>søgsmotor skal ses i lyset af, at de tidligere projekter alle<br />
sigtede på at udvikle kraftvarmeanlæg til biomasse, der efterfølgende blev installeret hos en<br />
vært, hvor anlægget skulle sættes i normal drift. Dermed var anlægget ikke til rådighed <strong>for</strong><br />
<strong>for</strong>søg, og der var oftest lang vej fra det team på DTU, der stod <strong>for</strong> udviklingen af<br />
stirlingmotorerne, til de værter, hvor anlæggene var installerede. Hensigten med at fremstille<br />
en stirlingmotor til laboratorie<strong>for</strong>søg var der<strong>for</strong> at fokusere på udvikling og afprøvning af<br />
<strong>for</strong>skellige <strong>for</strong>hold af betydning <strong>for</strong> stirlingmotorernes ydelse, virkningsgrad og drift.<br />
Resultater fra <strong>for</strong>søgsmotoren skulle yderligere anvendes til at evaluere<br />
simuleringsprogrammet ved bl.a. at afprøve specielle komponenter, der alene er fremstillet<br />
med henblik på at eftervise modellens godhed. Desuden ønskes en evaluering af<br />
stirlingmotoren som "prime mover" i små kraftvarmeinstallationer.<br />
Følgende emner ønskedes undersøgt:<br />
• Dellast. Undersøgelse af de tekniske muligheder <strong>for</strong> effektregulering af motorerne<br />
undersøges inklusive måling af <strong>for</strong>søgsmotorens egenskaber under dellast.<br />
• Test af ventiler. Undersøgelse af mulighederne <strong>for</strong> at <strong>for</strong>bedre ventilernes pålidelighed<br />
og levetid.<br />
• Test af slidmaterialer. Undersøgelse af de smøringsfri stempelringe og –tætninger, som<br />
anvendes i stirlingmotoren.<br />
• Lejetætninger. Undersøgelse af <strong>for</strong>skellige typer smørefedt til de nåle- og rullelejer, som<br />
benyttes i stirlingmotorerne, samt undersøgelse af lejetætningerne, der holder<br />
smørefedtet inde i lejerne.<br />
• Reduktion af støj. Undersøgelse af mulighederne <strong>for</strong> at reducere klaringerne i lejerne<br />
med henblik på reduktion af mekanisk støj fra motoren.<br />
• Brint i stedet <strong>for</strong> Helium. Undersøge af de materiale- og håndteringsmæssige<br />
problemer ved anvendelse af brint i <strong>for</strong>søgsmotoren.
PSO-FU2204 8 9 kW Stirlingmotor<br />
Forbedring af simuleringsmodeller. Forbedring af den eksisterende simuleringsmodel,<br />
der er udviklet på Inst. <strong>for</strong> Energiteknik, DTU. Der sigtes mod <strong>for</strong>bedring af beregningen<br />
af tryktab i motoren samt udvikling af bedre delmodeller af de enkelte komponenter.<br />
• Optimering af design-parametre. Optimering af effekt og virkningsgrad <strong>for</strong><br />
stirlingmotorer.<br />
• Vurdering af stirlingmotorteknologien. Ud fra resultaterne fra laboratorie<strong>for</strong>søg og<br />
beregninger vurderes stirlingmotorernes potentiale som basismodul i små<br />
kraftvarmeanlæg til biobrændsler.<br />
Budgettet taget i betragtning er listen over opgaver og emner, der ønskedes undersøgt,<br />
omfattende. Som følge af store mekaniske problemer med de to 9 kW motorer, der er<br />
udviklet i <strong>for</strong>bindelse med et tidligere projekt, har det heller ikke været muligt at gennemføre<br />
<strong>projektet</strong> som planlagt. Først langt inde i <strong>projektet</strong> var der skabt et overblik over årsagen til<br />
problemerne, og det var der<strong>for</strong> ikke muligt at nå at udvikle en ny motor, fremstille den og<br />
anvende den til de planlagte <strong>for</strong>søg. Projekt<strong>for</strong>løbet blev der<strong>for</strong> ændret, således at en af<br />
motorerne fra det tidligere projekt blev bragt i køreklar stand, og <strong>for</strong>søgene blev derefter<br />
gennemført med denne motor, samtidig med at den nye motor blev detailkonstrueret og<br />
komponenterne fremstillet.<br />
Projektet har yderligere været påvirket af, at der ikke i samme omfang som planlagt har<br />
været projektmidler til rådighed, så det har været muligt at fastholde en fast projektstab med<br />
9 kW motoren som fokus. Der har der<strong>for</strong> været en stor udskiftning af medarbejdere under<br />
projekt<strong>for</strong>løbet, hvilket har medført betydelige vanskeligheder <strong>for</strong> gennemførelsen af<br />
<strong>projektet</strong>.<br />
2 Udvikling af stirlingmotor til <strong>for</strong>søg<br />
Stirlingmotoren er baseret på en lukket termodynamisk proces, hvor det er den samme<br />
arbejdsgas, der gennemfører en cyklisk proces. Stirlingmotorens arbejdsprincip er nærmere<br />
beskrevet i bilag 1. Normalt benyttes Helium som arbejdsfluid, da denne gas har lav friktion<br />
og gode varmetransport egenskaber. Varmen overføres til og fra processen ved hjælp af<br />
varmevekslere. Forbrændingen <strong>for</strong>egår dermed ikke inde i motorens cylindre, men udvendigt<br />
som i en almindelig kedel, således at <strong>for</strong>brænding og mekanik er fuldstændig adskilt. Dette<br />
er en stor <strong>for</strong>del ved anvendelse af faste brændsler og andre vanskelige brændsler.<br />
2.1 Grundlæggende konstruktion af stirlingmotorerne til biobrændsler<br />
Adskillelsen af <strong>for</strong>brænding og mekanik gør det muligt at fremstille en hermetisk<br />
stirlingmotor, hvor generatoren er indbygget i et tryksat krumtaphus på samme måde som<br />
elmotoren er indbygget i en hermetisk kølekompressor.<br />
Sammenlignet med andre teknologier har den hermetiske stirlingmotor til små, decentrale<br />
kraftvarmeværker følgende <strong>for</strong>dele:<br />
• Simpel opbygning<br />
• Høj virkningsgrad i <strong>for</strong>hold til størrelsen.<br />
• Lave emissioner uden anvendelse af katalysatorer<br />
• Lang levetid med at lavt behov <strong>for</strong> service<br />
• Ingen reduktion af effekt og virkningsgrad med tiden<br />
• Mulighed <strong>for</strong> at anvendelse af mindre rene brændstoffer som f.eks. fyringsolie og<br />
biomasse
PSO-FU2204 9 9 kW Stirlingmotor<br />
Stirlingmotorens eksterne varmetilførsel gør den særligt velegnet til vanskelige brændsler.<br />
Varmen overføres fra <strong>for</strong>brændingen til motoren gennem en varmeveksler ved ca. 700 C, og<br />
partikler i røggasserne er der<strong>for</strong> ikke et stort problem. Det er der<strong>for</strong> en oplagt mulighed at<br />
anvende stirlingteknologien i små biomassefyrede kraftvarmeværker.<br />
Det særlige ved de dansk-udviklede stirlingmotorer er, at de i modsætning til andre<br />
stirlingmotorer er konstrueret direkte til at anvende brændsler, som er vanskelige at anvende<br />
i almindelige motorer. Det kan f.eks. være faste biobrændsler som flis eller træpiller eller<br />
gas<strong>for</strong>mige brændsler som biogas og <strong>for</strong>gasningsgas fra <strong>for</strong>gasning af træ. Det betyder, at<br />
varmebelastningen er lav i motorens varmeoverførende flader, og at der ikke <strong>for</strong>ekommer<br />
snævre passager, hvor partiklerne kan sætte sig og tilstoppe varmeveksleren.<br />
Motorernes øvrige komponenter er tilpasset de kraftige rør i kedeldelen og det store indre<br />
volumen. Det medfører, at motorernes cylindre har et <strong>for</strong>holdsvis stort cylindervolumen og et<br />
lavt middeltryk sammenlignet med stirlingmotorer, der er beregnet til naturgas eller dieselolie.<br />
Motorerne er yderligere konstrueret som hermetisk lukkede enheder, hvor alle de mekaniske<br />
dele sammen med generatoren er indbygget i et lukket hus, og den eneste <strong>for</strong>bindelse fra<br />
motorens indre til omgivelserne er kabel<strong>for</strong>bindelsen til el-nettet. På den måde undgås<br />
bevægelige tætninger, som har givet mange problemer i stirlingmotorer med krumtaphuse, der<br />
har <strong>for</strong>bindelse til omgivelserne.<br />
Stempelringe og -tætninger er fremstillet af smøringsfri PTFE-materialer. En helt ny type<br />
krumtapmekanisme sørger <strong>for</strong>, at sliddet på stempelringene er lille. Olie i motoren er helt<br />
undgået ved at anvende fedtsmurte rulle- og nålelejer. Der skal dermed hverken skiftes olie<br />
eller tændrør i motoren, og behovet <strong>for</strong> service er der<strong>for</strong> reduceret til et minimum.<br />
Som tidligere beskrevet er der udviklet to motortyper. Den største motor har 4<br />
dobbeltvirkende cylindre og en nominel el-effekt på 35 kW. Den integrerede 6-polede<br />
asynkron generator har en nominel el-effekt på 37 kW, hvilket medfører et omdrejningstal lidt<br />
over 1000 omdrejninger per min. Arbejdsprocessen, der anvender Helium som arbejdsgas,<br />
har et middeltryk på 4,5 MPa. De fire cylindre er placeret i hvert hjørne af et kvadrat, hvor de<br />
fire hederpaneler, der udgør den varme varmeveksler, <strong>for</strong> hver af de fire cylindre danner et<br />
kvadratisk kammer, der er placeret <strong>for</strong> enden af fyringssystemets bræn<strong>dk</strong>ammer. Der er til<br />
dato produceret 9 motorer i 5 generationer. Den seneste model, SM3E, er fremstillet i 3<br />
eksemplarer, og en af motorerne har nu kørt over 4000 timer med særdeles godt resultat.<br />
Den maksimale el-effekt ved nominel driftstilstand er 40 kW, og den interne virkningsgrad er<br />
34% (den interne virkningsgrad defineres som motorens akseleffekt i <strong>for</strong>hold til den<br />
varmeeffekt, der tilføres i motorens varme varmeveksler, hederen). Det er en <strong>for</strong>bedring på<br />
15% i <strong>for</strong>hold til de tidligere modeller, hvilket først og fremmest er muliggjort som følge af<br />
resultaterne fra dette projekt.<br />
Den mindre motor, SM5A, har en enkelt cylinder og en nominel el-effekt på 9 kW. Motoren er<br />
konstrueret efter de samme grundlæggende principper som den oven<strong>for</strong> beskrevne firecylindrede<br />
motor, og mange af erfaringerne fra konstruktion og drift kan udveksles mellem de<br />
to motortyper. Der blev fremstillet to prototyper af den første udgave af 9 kW-motoren, SM5A<br />
(se ref. 1), og det var meningen at disse to motorer skulle langtidsafprøves med henholdsvis<br />
biogas og naturgas som brændsel. Denne målsætning viste sig imidlertid at være <strong>for</strong><br />
optimistisk <strong>for</strong> de første prototyper af en hel ny motortype, <strong>for</strong>di problemer med motorerne<br />
<strong>for</strong>hindrede gennemførelsen af det planlagte testprogram.<br />
Ifølge projektbeskrivelsen bestod en del af nærværende <strong>projektet</strong> i at gennemføre en<br />
revision af 9 kW motorens design, så den bliver så enkel som mulig at samle og adskille, og<br />
efterfølgende at fremstille og afprøve en ny motor, som alene skulle anvendes til <strong>for</strong>søg.
PSO-FU2204 10 9 kW Stirlingmotor<br />
Neden<strong>for</strong> er de én-cylindrede motorers opbygning beskrevet i detaljer samtidig med, at<br />
revisionen af motorens konstruktion er behandlet.<br />
2.2 Valg af <strong>for</strong>søgsmotor<br />
Valget af 9 kW-el motoren som udgangspunkt <strong>for</strong> <strong>for</strong>søgsmotoren var betinget af, at det er<br />
den billigste løsning både med hensyn til afprøvning af nye komponenter og med hensyn til<br />
energi<strong>for</strong>brug ved driftsafprøvninger. De større motorer med flere cylindre er mere<br />
tidskrævende at samle og adskille, og det vil ofte være nødvendigt at fremstille nye<br />
<strong>for</strong>søgskomponenter i fire eksemplarer i stedet <strong>for</strong> én. Yderligere kan resultater fra <strong>for</strong>søg<br />
med den én-cylindrede stirlingmotor med en el-effekt på 9 kW direkte overføres til<br />
udviklingen af den fire-cylindrede stirlingmotor med en el-effekt på 35 kW.<br />
• Det er vanskeligt at få den samme driftstilstand på alle cylindre på samme tid i en motor<br />
med mere end én cylinder, og det er der<strong>for</strong> svært at opnå en veldefineret driftstilstand<br />
med en fire-cylindret motor. En én-cylindrede motor er der<strong>for</strong> mere velegnet til<br />
laboratorie<strong>for</strong>søg end en firecylindret motor<br />
• Stempelringe og -tætninger samt stempelstangstætninger har samme ud<strong>for</strong>mning og er<br />
fremstillet af de samme PTFE-materialer i begge motorer. Yderligere er<br />
middelstempelhastighed og den trykvariation, som tætningerne er udsat <strong>for</strong>, identiske <strong>for</strong><br />
de 2 motortyper. Driftsresultater vedrørende slid af PTFE-ringene kan der<strong>for</strong> direkte<br />
overføres fra den én-cylindrede motor til den fire-cylindrede motor.<br />
• Der benyttes de samme én-vejsventiler og bypass-ventiler i de to motortyper. I den éncylindrede<br />
motor udsættes én-vejsventilerne <strong>for</strong> det dobbelte maksimaltryk i <strong>for</strong>hold til<br />
den firecylindrede motor, og ventilerne udsættes dermed <strong>for</strong> en <strong>for</strong>m <strong>for</strong> accelereret test.<br />
Pålidelighed og levetid <strong>for</strong> disse komponenter er af stor betydning <strong>for</strong> 35 kW-motorens<br />
drift.<br />
• Topstykker, hederrør og stempler er fremstillet af de samme varmefaste materialer.<br />
Erfaringer vedrørende krybning som følge af termiske spændinger samt "low cycle<br />
fatigue" kan der<strong>for</strong> overføres fra den ene motor til den anden.<br />
• Begge krumtapekanismerne er opbygget efter de samme grundlæggende principper, og<br />
der anvendes den samme type og størrelse nålelejer i begge motortyper på de<br />
leje<strong>for</strong>bindelser, hvor der <strong>for</strong>ekommer vippe bevægelser. Hovedlejerne, som roterer med<br />
konstant vinkelhastighed, er <strong>for</strong>skellige, men i begge motorer er der tale om rullelejer<br />
med samme type læbetætninger til at holde fedtet inde i lejerne. Desuden er<br />
problemstillingen med fedtsmøring af nåle- og rullelejer i et hermetisk krumtaphus med<br />
Helium under tryk det samme i begge motortyper.<br />
Forsøg med nye PTFE-ringe til stempel- og stempelstangstætninger gennemføres også<br />
bedst på den én-cylindrede motor, <strong>for</strong>di det er nemmere at adskille og samle motoren med<br />
henblik på inspektion af de testede komponenter.<br />
Det var hensigten at anvende erfaringerne fra de to første prototyper af den én-cylindrede<br />
motor til at udvikle en ny motor med samme grundlæggende ud<strong>for</strong>mning, hvor der var taget<br />
hensyn til den oprindelige konstruktions mindre hensigtsmæssige delløsninger. For at undgå<br />
problemer med den nye motor var det ikke hensigtsmæssigt at fastlægge konstruktionen af<br />
denne motor før de to eksisterende motorer kørte nogenlunde problemfrit. Der tog imidlertid<br />
betydeligt længere tid end ventet at opnå dette, og det var først i efteråret 2003, at årsagerne<br />
til problemer med de eksisterende motorer var kortlagt og nye løsninger identificeret.<br />
Årsagen til en væsentlig del af problemerne var en uhensigtsmæssig ud<strong>for</strong>mning af de<br />
eksisterende motorers krumtapmekanisme, der førte til høje belastninger på de<br />
komponenter, der fikserer hovedleje og svingarmenes akseltappe i motorblokken. De høje
PSO-FU2204 11 9 kW Stirlingmotor<br />
belastninger samt fejl under fremstillingen af motorernes komponenter var den direkte årsag<br />
til adskillige havarier under driftsafprøvningen.<br />
Motorens kompakte design medførte andre problemer, idet de <strong>for</strong>holdsvis store<br />
trykvariationer i krumtaphuset medførte, at smørefedt blev pumpet ud af lejerne. Det var dels<br />
et problem <strong>for</strong> lejerne, men det var også et problem <strong>for</strong> stangtætningen, der blev ødelagt af<br />
fedt, der ramte stempelstangen. Der har yderligere været problemer med regulering af<br />
trykket i krumtaphuset samt med slaglodningen af komponenter i krumtapmekanismen. Også<br />
svejseprocessen, der er anvendt til fastgøring af hederrørene til topstykket, har givet<br />
problemer.<br />
Mange af problemerne har ført til havarier, og det har været tidskrævende at finde nye<br />
løsninger og fremstille nye komponenter, <strong>for</strong> derefter at afprøve dem i motoren. Erfaringerne<br />
fra afprøvningen af de eksisterende motorer har ført til følgende konklusioner:<br />
Krumtapmekanisme:<br />
• Belastningerne på bagpladen i krumtaphuset er uhensigtsmæssigt store<br />
• Momentet på de akseltappe, der fikserer svingarmene til bagpladen, skal reduceres<br />
eller elimineres.<br />
• Som følge af motorens kompakte design er trykvariationen i krumtaphuset under drift<br />
<strong>for</strong>holdsvis stort, hvilket medfører, at trykvariationerne pumper fedtet ud af lejerne<br />
• Når lejerne, der <strong>for</strong>binder åget med arbejdsstemplets plejlstang, er udluftet ind mod<br />
stempelstangen, medfører trykvariationerne i krumtaphuset, at smørefedtet bliver<br />
”skudt ud” på stempelstangen, således at stangtætningerne bliver ødelagt.<br />
• Den eksisterende styring af krumtaphusets tryk i <strong>for</strong>hold til middeltrykket i cylinderen<br />
fungerer ikke efter hensigten<br />
• Slaglodning af dele i mekanismen bør undgås af hensyn til styrken<br />
• Motorens svinghjul er <strong>for</strong> lille, hvilket medfører at motoren får en ujævn gang, hvilket<br />
igen fører til en lav cos og en høj maksimalværdi <strong>for</strong> strømmen.<br />
• For at opnå optimal fasevinkel mellem stempel positioner er der monteret en exentrik<br />
på krumtapakslen. Denne konstruktion øger komplexiteten af mekanismen væsentligt.<br />
Topstykke og brænder<br />
• Den proces, der er benyttet til svejsning af hederrørene til topstykket, er ikke velegnet.<br />
• Luft<strong>for</strong>varmerens ud<strong>for</strong>mning er uhensigtsmæssig, da kapacitetsstrømmen af<br />
luftstrømmen på den udvendige side af varmeveklserens rør ikke er jævnt <strong>for</strong>delt over<br />
alle rørene, mens <strong>for</strong>delingen af røggas på rørenes indvendige side er mere jævnt<br />
<strong>for</strong>delt. Dette fører til tab ef effektivitet <strong>for</strong> lufoen<br />
Generator<br />
• Strømgennemføring til generatoren skal flyttes, så der bliver mere plads til montage og<br />
ledningsføring.<br />
2.3 Design af 9 kW <strong>for</strong>søgsmotor<br />
Den nye motor, der er udviklet i <strong>for</strong>bindelse med nærværende projekt, er som oven<strong>for</strong> nævnt<br />
baseret på erfaringerne fra de første 9 kW motorer (typebetegnelse SM5A). Den nye motor<br />
betegnes i det følgende SM5B. Det er en én-cylindret stirlingmotor af ¡ -typen, hvor<br />
arbejdsstempel og <strong>for</strong>trængerstempel er placeret i samme cylinder (se bilag 1). Motoren er<br />
som de tidligere motorer konstrueret som en hermetisk enhed, hvor generatoren og alle<br />
stirlingmotorens mekaniske dele er bygget ind i en lukket beholder, som udgør motorens<br />
krumtaphus. Den eneste <strong>for</strong>bindelse fra motorens indre til omgivelserne er kabel<strong>for</strong>bindelsen<br />
til el-nettet. På den måde undgås bevægelige tætninger, som har givet mange problemer i
PSO-FU2204 12 9 kW Stirlingmotor<br />
stirlingmotorer med konventionelle åbne<br />
krumtaphuse. Den integrerede asynkron<br />
generator har 6 poler svarende til et<br />
omdrejningstal på ca. 1000 min-1, og<br />
generatorens nominelle el-effekt er 11 kW.<br />
Figur 2.1 viser princippet i den hermetiske<br />
stirlingmotor af ¡ -typen.<br />
For at undgå problemer med olie, der<br />
trænger op i arbejdsvolumenerne fra<br />
krumtaphuset, anvendes fedtsmurte rulle-<br />
og nålelejer. Stempelringe og<br />
stangtætninger er fremstillet af PTFE med<br />
<strong>for</strong>skellige <strong>for</strong>mer <strong>for</strong> fyldstoffer. Motoren<br />
anvender Helium som arbejdsgas ved et<br />
middeltryk på 8 MPa. Motoren er udviklet til<br />
først og fremmest at anvende gas<strong>for</strong>mige<br />
brændsler, og dette er baggrunden <strong>for</strong>, at<br />
middeltrykket er noget højere end i den 4cylindrede<br />
35 kW motor.<br />
Den nye 9 kW motors grundlæggende ud<strong>for</strong>mning er ikke ændret i <strong>for</strong>hold til de første<br />
motorer. Det betyder at motorens boring og slaglængde <strong>for</strong> både arbejdsstempel og<br />
<strong>for</strong>trænger er uændret, ligesom der heller ikke er ændret på antal og dimensioner af rørene i<br />
køler og heder samt regeneratorens dimensioner. Som følge af de i <strong>for</strong>egående afsnit<br />
nævnte konklusioner er der imidlertid lavet en omfattende revision af motoren.<br />
Konstruktionen er ændret på følgende punkter:<br />
Krumtapmekanisme:<br />
Figur 2.1 Principskitse af hermetisk en-cylindret<br />
motor.<br />
• Krumtapmekanismen er redesignet med udgangspunkt i en ny idé til en ændret<br />
geometri. Der er opnået følgende <strong>for</strong>bedringer:<br />
1. Den reviderede krumtapmekanisme medfører, at den ønskede fasevinkel<br />
mellem <strong>for</strong>trænger- og arbejdsstempel er en direkte følge af mekanismens<br />
geometri, hvilket medfører en væsentlig simplifikation af mekanismen både<br />
med hensyn til antallet af komponenter, men også montage og demontage.<br />
2. Den reviderede krumtapmekanisme medfører en væsentlig reduktion af<br />
kræfterne på mekanismens komponenter og dermed også på bagpladen.<br />
3. Momentet på akseltappene til svingarmene er elimineret, og akseltappene er<br />
ikke længere monteret på bagpladen, men direkte i krumtaphuset. Dette<br />
reducerer yderligere kræfterne på bagpladen.<br />
• Alle lejer er udluftet i omdrejningscentret, og lejerne, der <strong>for</strong>binder åget med<br />
arbejdsstemplets plejlstang, udluftes ikke ind mod stempelstangen men direkte til<br />
krumtaphuset. Dermed er problemerne med at smørefedt pumpes ud af lejerne som<br />
følge af trykvariationer i krumtaphuset løst. Desuden kan fedtet ikke blive ”skudt ud” på<br />
<strong>for</strong>trængerens stempelstang.<br />
• En ny styring af trykket i krumtaphuset i <strong>for</strong>hold til arbejdsvolumenernes middeltryk er<br />
udviklet og afprøvet på de eksisterende motorer. Styringen består af en kontrollerbar<br />
læk, der <strong>for</strong>binder de to volumener.<br />
• Åg og svingarme fremstilles af højstyrkestål (Weldox 700), hvor den udvendige kontur<br />
vandstråleskæres, så lodning helt undgås.<br />
• Den nye motors svinghjul er mere end dobbelt så stort som svinghjulet på de<br />
oprindelige motorer.
PSO-FU2204 13 9 kW Stirlingmotor<br />
Topstykke og brænder:<br />
• En ny metode til montering af hederrør på det varme topstykke er identificeret og<br />
implementeret på den nye motor. Dermed er svejseprocessen simplificeret, så<br />
vanskeligheder med svejserevner ved svejsning af rørene til topstykket kan undgås.<br />
• En helt ny luft<strong>for</strong>varmer er udviklet, hvor strømningen af både luft og røggas er jævnere<br />
<strong>for</strong>delt.<br />
Generator:<br />
• Strømgennemføringen til generatoren er redesignet på basis af de gode erfaringer fra<br />
strømgennemføringerne på de større 35 kW motorer.<br />
Mens den første motor, SM3A, var konstrueret og detailtegnet ved hjælp af 2D-CADprogrammet<br />
AutoCad, er hele konstruktionsarbejdet vedrørende den nye motor gennemført i<br />
3D-CAD-programmet AutoDesk Inventor. Det har givet nogle udvidede muligheder <strong>for</strong> at lave<br />
en hensigtsmæssig konstruktion, men det har samtidig været et større arbejde end <strong>for</strong>udset<br />
at lave hele tegnearbejdet fra grunden. Det er nemlig ikke muligt at overføre tegninger fra<br />
AutoCad til Inventor, <strong>for</strong>di der ligger en helt anden tankegang til grund <strong>for</strong> at arbejde i 3D-<br />
CAD-programmet. I alt er der lavet ca. 200 tegninger til fremstillingen af motoren, og meget<br />
få af disse tegninger er uændrede i <strong>for</strong>hold til den tidligere type SM5A-motor (se tegningsliste<br />
og detailtegningen af enkelte hove<strong>dk</strong>omponenter i bilag 3)<br />
Figur 2.2 og 2.3 viser den nye motors opbygning. Neden<strong>for</strong> er den nye konstruktion<br />
kommenteret i detaljer: Motorens grundlæggende data fremgår af tabel 2.1. Sammenlignet<br />
med andre stirlingmotorer i effektområdet fra 2 kW til 20 kW fremstår den nye konstruktion<br />
som både mere enkel, mere kompakt og mere robust. Neden<strong>for</strong> er motorens konstruktion<br />
beskrevet i detaljer.
Luft ind<br />
PSO-FU2204<br />
Gas<br />
14 9 kW Stirlingmotor<br />
Forvarmet<br />
Luft<br />
Røggas<br />
Forbrænding<br />
Røggas<br />
ud<br />
9<br />
10<br />
11<br />
12<br />
13<br />
14<br />
15<br />
1<br />
0 100 200 300 400 500 600 700<br />
Figur 2.2: Tegning af SM3B <strong>for</strong>fra.<br />
2<br />
3<br />
4<br />
5<br />
6<br />
7<br />
8<br />
16<br />
18<br />
Pos. Description<br />
1 Gas brænder<br />
2 Hederrør<br />
3 Expansionsvolumen<br />
4 Fortrængerstempel<br />
5 Topstykke<br />
6 Regenerator<br />
7 Køler<br />
8 Bypassventil<br />
9 Luft<strong>for</strong>varmer<br />
10 Kølevand ud<br />
11 Kompressionsvolumen<br />
12 Arbejdsstempel<br />
13 Bagplade til mekanisme<br />
14 Arbejdsåg<br />
15 Fortrængeråg<br />
16 Arbejdssvingarm<br />
17 Hermetisk signalgennemføring<br />
18 Fortrængersvingarm<br />
19 Dæksel<br />
20 Svinghjul<br />
21 6 polet asynkron generator<br />
22 Kølevand ind<br />
23 Krumtap<br />
24 Balancevægt<br />
25 Krumtaphus<br />
26 Hermetisk effektudtag<br />
17<br />
mm
19<br />
PSO-FU2204 15 9 kW Stirlingmotor<br />
Tabel 2.1: Specifikationer <strong>for</strong> 9 kWel<br />
stirlingmotor, SM3B<br />
Boring, mm 114<br />
Slag, mm 54<br />
Fasevinkel, o 79<br />
Antal cylindre 1<br />
Omdrejningshastighed, rpm 1015<br />
Middeltryk, MPa 8.0<br />
Arbejdsgas Helium<br />
Generatorvirkningsgrad, % 89<br />
Driftsdata ved nominel driftstilstand:<br />
Hederrørstemperatur, o C 740<br />
Kølevandstemperatur, frem, o C 55<br />
Akseleffekt, kW 10<br />
El-effekt, kW 9<br />
Vægt inkl. brænder, kg 600<br />
Længde, m 1,0<br />
Bredde, m 0,6<br />
Højde, m 1,75<br />
20 21 22 23 24 25<br />
0 100 200 300 400 500 600 700 800 900 1000<br />
Figur 2.3: Snittegning af SM3B fra siden. Bemærk at snittet er <strong>for</strong>sat, således at man ser centeret i cylinderen<br />
samtidig med centeret i krumtappen.<br />
26<br />
mm
PSO-FU2204 16 9 kW Stirlingmotor<br />
2.4 Cylinder med topstykke, heder, regenerator og køler<br />
Motorens interne varmevekslere består af en heder, hvor varmen tilføres, en regenerator,<br />
hvor varmen lagres, og en køler, hvor den overskydende varme fra processen fjernes. De tre<br />
varmevekslere er anbragt som en rotationssymmetrisk ring om <strong>for</strong>trængerstemplets cylinder,<br />
således at kølevandet i køleren samtidig sørger <strong>for</strong> at holde temperaturen på <strong>for</strong>trængerens<br />
cylinderoverflade på et lavt niveau.<br />
Figur 2.4: Snit I topstyke som det er på de to første motorer. Bemærk at kuplen er lavet som en tolags<br />
konstruktion, og at hederrørene er svejset direkte i kupelmaterialet.<br />
Figur 2.5: Motorens topstykke består af en kuplet endebund I 2 lag, hvor der er påsvejset en flange<br />
og boret huller til montering af hederrørene. I det nye topstykke samles 4 hederrør i et samlestykke til<br />
et fælles tilslutningsrør, således at de 32 hederør kun medfører 8 svejste rørtilslutninger i topstykket.
PSO-FU2204 17 9 kW Stirlingmotor<br />
De U-<strong>for</strong>mede hederrør <strong>for</strong>binder det varme ekspansionsvolumen over <strong>for</strong>trængeren med<br />
regeneratoren, der sidder i en ring om <strong>for</strong>trængercylinderen. Hederen består af 32 rør, ø8 x<br />
2,5 mm, i materialet Alloy 800H, der er bukket og monteret på topstykket ved at svejse<br />
rørenderne indvendigt. Før svejsningen er der monteret finner. På grund af problemer i den<br />
oprindelige konstruktion med svejserevner i svejsningen af hederrørene til topstykket er<br />
hederen <strong>for</strong>søgt ud<strong>for</strong>met på en måde, så dette problem reduceres. Figur 2.5 viser den nye<br />
konstruktion og til sammenligning er den gamle konstruktion vist på fig. 2.4. Som det fremgår<br />
af figuren er der på den nye konstruktion indført nogle samlestykker, hvor fire hederrør<br />
samles til et større rør. Mellemstykke, finner og hederrør loddes sammen i vakuum med<br />
nikkel som lod. Dermed er der<br />
kun 2 x 8 rørender, der skal<br />
svejses til topstykket, og<br />
rørendernes godstykkelse er<br />
større. Finnerne er fælles <strong>for</strong><br />
alle fire rør i et sæt, således at<br />
hederen er 8-kantet, når den<br />
er samlet.<br />
Figur 2.6: Køler til 9kW motor. De to huller på <strong>for</strong>siden er ind-<br />
og udløb til vand. De mange små huller der ses på oversiden er<br />
enderne på Heliumrørene, rørne kan også ses igennem<br />
hullerne på <strong>for</strong>siden. Fortrængerstemplet tætner og styrer<br />
indvendig i den store center boring.<br />
Topstykkets kuplede<br />
endebund er koldpresset af 2<br />
lag plade. Efter påsvejsning af<br />
flange, hederrør og<br />
indercylinder, varmebehandles<br />
topstykket, hvorefter det<br />
bearbejdes til de færdige mål.<br />
Der er yderligere svejset en<br />
plade fast på flangen til<br />
montering af gasbrænderen.<br />
Regeneratoren er fremstillet af<br />
tynd rustfri tråd, der er<br />
behandlet i en særlig proces<br />
og presset til den endelige<br />
<strong>for</strong>m i et presseværktøj.<br />
Køleren er af ”tube and shell”<br />
typen med Helium<br />
arbejdsgassen indvendigt i<br />
rørene, og kølevandet<br />
strømmende på tværs uden<br />
på rørene (se figur 2.6). Der er<br />
brugt ressourcer på at prøve<br />
at finde enklere løsninger end<br />
tynde rustfri rør, der er vacuumloddet i et rustfrit hus, men dette har ikke medført store<br />
<strong>for</strong>bedringer.<br />
Køleren er monteret i en mellemblok, således at det er muligt at demontere topstykke med<br />
brænder og køler ved blot at løsne 4 bolte. Når hele cylinderen med topstykke, køler og<br />
brænder er demonteret, kan man umiddelbart skifte stempeltætninger og stangtætninger.<br />
Det betyder, at et service er nemt og hurtigt at lave, hvilket er vigtigt, <strong>for</strong>di det <strong>for</strong>ventes at<br />
det er nødvendigt at skifte disse tætninger efter ca. 8.000 timers drift. Hvis en almindelig<br />
husholdningsbil kører 8.000 timer i en tilsvarende driftstilstand (fuld last og moderat<br />
omdrejningstal) svarer det til ca. 1 mill. km mellem service.
PSO-FU2204 18 9 kW Stirlingmotor<br />
2.5 Krumtapmekanisme<br />
Krumtapmekanismen i den én-cylindrede motor er baseret på to fire-ledsmekanismer, der<br />
vender henholdsvis opad og nedad som vist på figur 2.7. Hver af de to fireledsmekanismer<br />
består af et trekantet åg og to drejelige <strong>for</strong>bindelser til krumtaphuset. De tre hjørner på et åg<br />
er <strong>for</strong>bundet til hhv. krumtaphuset via en arm med en fast længde, krumtapakslen, og en<br />
<strong>for</strong>bindelse til stemplet. De to mekanismer er ens, idet de er spejlet omkring et vandret plan.<br />
Fortrænger<br />
Arbejds<br />
stempel<br />
Figur 2.7: Den samlede krumtapmekanismes<br />
struktur .<br />
(Xc,Yc)<br />
Figur 2.8 Grundmekanismens struktur.<br />
Den generelle ud<strong>for</strong>mning af<br />
fireledsmekanismen er vist på figur 2.8.<br />
Trekanten med siderne b og R og vinklen<br />
t er <strong>for</strong>bundet til to faste<br />
omdrejningspunkter O2 og O4 ved hjælp<br />
af armene a og c. Når vinklen alfa<br />
ændres, beskriver koble-punktet (xc, yc)<br />
et kurve<strong>for</strong>løb. Forbindes koble-punktet<br />
(xc, yc) til stemplet, mens armen a er<br />
krumtapbugten, der roterer om O2, er det<br />
muligt at bestemme længden af a, b, c, d<br />
og R således, at punktet (xc, yc) bevæger<br />
sig tilnærmelsesvis på en ret linje. Hvor<br />
lille fejlen er, afhænger af størrelsen af<br />
åget i <strong>for</strong>hold til slaglængden. I<br />
designfasen må man altså beslutte hvor<br />
stort krumtaphuset må blive, <strong>for</strong> at<br />
<strong>for</strong>bedre liniariteten af stempelføringen.<br />
Mekanismens dimensioner er optimeret<br />
således, at punktet (xc, yc), der er<br />
<strong>for</strong>bundet til stemplet, bevæger sig<br />
inden<strong>for</strong> et rektangel med en bredde der<br />
er 200 gange længden. Slaglængden er<br />
ca. 2 2 gange radius på den roterende<br />
arm, dvs. krumtappens <strong>for</strong>sætning = a.<br />
Fasevinklen mellem de to stemplers position<br />
er givet ved summen af vinklerne t.
PSO-FU2204 19 9 kW Stirlingmotor<br />
Figur 2.9: Gammel (øverst) og ny mekanisme.<br />
Figur 2.10:<br />
Mekanismens fysiske<br />
ud<strong>for</strong>mning<br />
Gamle mål:<br />
a 19<br />
b 95,96<br />
c 95,9<br />
d 134,8<br />
R 136,03<br />
t 45°<br />
G 95<br />
Slag 54,144<br />
Side 0,2798<br />
Nye mål:<br />
a 20,7<br />
b 97<br />
c 116,72<br />
d 151<br />
R 125,96<br />
t 39,6367°<br />
G 79<br />
Slag 54,079<br />
Side 0,2713<br />
Figur 2.9 (øverst) viser<br />
geometrien i de to første<br />
motorers krumtapmekanisme<br />
sammenlignet med den nye<br />
mekanisme. På den øverste,<br />
gamle, geometri ses det at<br />
vinklen t er 45°. Ser man på<br />
begge ågs bevægelse, er<br />
fasevinklen mellem<br />
topdødpunktet <strong>for</strong> stemplerne<br />
90 o . Dette er lidt i overkanten<br />
af, hvad der er<br />
hensigtsmæssigt, og den ene<br />
mekanisme er der<strong>for</strong> <strong>for</strong>skudt<br />
14 grader på krumtapbugten,<br />
med en eksentrik, således at<br />
den faktiske fasevinkel<br />
mellem <strong>for</strong>trænger og<br />
arbejdsstempel er 76 o .<br />
Den nye mekanisme har en<br />
værdi af ”t” på 39,63°,<br />
svarende til en fasevinkel på<br />
79,27° uden eksentrik. Den<br />
termodynamiske konsekvens<br />
af ændringen er lille, men de<br />
konstruktionsmæssige <strong>for</strong>dele<br />
er meget store. Bemærk at<br />
værdien af ”G” er reduceret<br />
fra 95 til 79mm. Det betyder at<br />
cylinderaksen kommer tættere på krumtaphusets<br />
symmetriakse. Det giver en mere harmonisk motor med<br />
mindre vibrationer. Slaglængde og sidebevægelse er næsten<br />
uændrede. Karrakteristisk <strong>for</strong> den nye mekanisme er at<br />
svingarmene skal lejres i større afstand fra krumtapakslens<br />
center, end i den gamle mekanisme. Det har medført en mere<br />
kompleks bearbejdning af krumtaphuset. Tidligere er nævnt<br />
problemet med stivheden af svingarmens befæstning på<br />
bagpladen. Dette problem elimineres ved at fæstne<br />
svingarmsakslen direkte i krumtaphuset. Således er der god<br />
sammenhæng i det nye design.<br />
Koblepunktets bevægelse tilnærmelsesvis på en ret linje<br />
medfører, at det er muligt at undgå krydshoveder i motoren til<br />
at optage de kræfter vinkelret på bevægelsesretningen, der er<br />
en følge af, at kræfterne ikke overføres i samme retning som<br />
cylinderens centerlinje. For en oliefri krumtapmekanisme som<br />
den her anvendte er det en særlig <strong>for</strong>del at undgå<br />
krydshoveder, <strong>for</strong>di disse kræver en eller anden <strong>for</strong>m <strong>for</strong><br />
smøring.
PSO-FU2204 20 9 kW Stirlingmotor<br />
Figur 2.11: Foto af ny mekanisme.<br />
Det kan ses af figur 2.11 at lejetapperne til krumtapmekanismens to svingarme ikke er<br />
monteret på bagpladen som i den tidligere motor, men direkte i krumtaphuset, idet der er<br />
bearbejdet to ”øjer” direkte i krumtaphuset.<br />
Figur 2.9 viser resultatet af optimeringen af fire-ledsmekanismen med en fasevinkel på 2 x<br />
39 o = 78 o og den ønskede slaglængde på 54 mm. Koblepunktet W’s bevægelse vinkelret på<br />
stempelbevægelsens længderetning er mindre end ± 0,15 mm, hvilket svarer til en<br />
tværbevægelse, der udgør mindre end 0,5 o /oo af slaglængden. Den mindre vinkel i ågets<br />
spidse trekant medfører, at <strong>for</strong>holdet mellem slaglængden af koblepunktet og<br />
krumtapbugtens diameter bliver mindre, og maksimalkraften på hovedlejerne tilsvarende<br />
mindre. Tabellen neden<strong>for</strong> viser en sammenligning af den gamle og den nye mekanisme:<br />
Gamle mek. Nye mek.<br />
Slaglængde, mm 54 54<br />
Tværbevægelse, mm 0,28 0,27<br />
Koblepunktets fase<strong>for</strong>skydning, o 45 39,6<br />
Maksimalkraft på stempel, kN 15 15<br />
Maksimalkraft på plejlstangsleje, N 21,2 19,3<br />
Akseltappe<br />
Som det fremgår af tabellen er maksimalkraften på plejlstangslejerne reduceret med 10%, og<br />
bortset fra lejerne på arbejdsstemplets <strong>for</strong>bindelse til åget er kræfterne på alle lejer reduceret<br />
med 10%. For de anvendte nåle- og rullelejer betyder en 10% reduktion af belastningen en<br />
<strong>for</strong>øgelse af levetiden på næsten 40%.<br />
For at opnå den ønskede slaglængde på 54 mm <strong>for</strong> den nye krumtapmekanisme, er<br />
krumtapbugtens <strong>for</strong>sætning øget fra 19 mm til 20,7 mm, hvilket bl.a. har medført, at det var<br />
nødvendigt at anvende et større hovedleje end tidligere.
PSO-FU2204 21 9 kW Stirlingmotor<br />
Figurene 2.3, 2.10 og 2.11 viser den praktiske ud<strong>for</strong>mning af krumtapmekanismen i den nye<br />
motor. Sammenlignes med de første motorer som vist på figur 2.9 ses, at cylinderens<br />
centerakse er flyttet nærmere til krumtappens omdrejningspunkt, mens de to svingarme, der<br />
<strong>for</strong>binder åget med krumtaphuset, er blevet længere. Netop de længere svingarme er den<br />
vigtigste grund til at ændre ud<strong>for</strong>mningen af krumtapmekanismen. Som det ses er<br />
svingarmenes lejring nu en fast del af krumtaphuset i modsætning til den gamle ud<strong>for</strong>mning,<br />
hvor svingarmenes lejring var boltet fast på bagpladen. Dermed er hovedårsagen til de<br />
mange havarier elimineret.<br />
I den oprindelige krumtapmekanismes konstruktion var der inkluderet en balancemasse, der<br />
skulle reducere vibrationerne som følge af stemplernes massekræfter i lodret retning.<br />
Balancemassen var placeret på en <strong>for</strong>længelse af det åg-par, som <strong>for</strong>bandt arbejdsstemplet<br />
til krumtappen. En nærmere analyse af de dynamiske kræfter i konstruktionen viste, at<br />
balancemassen skulle være 5 gange tungere <strong>for</strong> at bidrage væsentligt til reduktion af<br />
vibrationerne. Det er der slet ikke plads til i krumtaphuset, og det blev der<strong>for</strong> besluttet helt at<br />
fjerne den. Dermed var endnu en mulig årsag til havari fjernet, da balancemassen var faldet<br />
af under langtidsafprøvningen af den motor, der kørte med biogas som brændsel.<br />
2.6 Krumtaphus<br />
Krumtaphuset er ændret betydeligt i <strong>for</strong>hold til den første model. Dels krævede den nye<br />
mekanisme, at afstanden mellem cylinderens centerlinje og krumtappens centerlinje blev<br />
reduceret, og dels er svingarmene, der <strong>for</strong>binder ågene med et fast punkt, nu lejret direkte i<br />
krumtaphuset.<br />
Svinghjulet er monteret på krumtappen mellem hovedlejet ved kumtapbugten og rotoren til<br />
generatoren. For at <strong>for</strong>bedre kvaliteten på den producerede strøm, er svinghjulets tykkelse<br />
mere end <strong>for</strong>doblet, og krumtaphuset er <strong>for</strong>længet tilsvarende. For at reducere vægten af<br />
krumtaphuset og <strong>for</strong> at undgå problemer med de<strong>for</strong>mationer af endebunden, er den flade<br />
endebund erstattet med en kuplet endebund og en separat fikstur til kuglelejet <strong>for</strong> enden af<br />
krumtappen. Samtidig er nettilslutningens gennemføring placeret i midten af endebunden,<br />
således et der er god plads til ledningsføringen på indersiden af krumtaphusets trykbærende<br />
væg.<br />
Figur 2.13 viser.et foto af den nye 9 kW motor, SM5B, monteret i en ramme med gasrampe,<br />
blæser og regulering, og figur 2.12 viser et print fra 3D-CAD programmet fra samme vinkel.<br />
Resultatet af udviklings<strong>projektet</strong> er en gennemgribende rekonstruktion af motorens design,<br />
som <strong>for</strong>ventes at være væsentlig <strong>for</strong>bedret på en række punkter.
PSO-FU2204 22 9 kW Stirlingmotor<br />
Figur 2.12. Den nye 9 kW motor<br />
modeleret ved hjælp af 3D-CAD<br />
programmet Inventor<br />
Figur 2.13. Den nye 9 kW motor<br />
fotograferet efter samling.
PSO-FU2204 23 9 kW Stirlingmotor<br />
3 Laboratorie<strong>for</strong>søg<br />
Som tidligere nævnt var det<br />
hensigten at udvikle en ny<br />
motor til <strong>for</strong>søg baseret på<br />
erfaringerne fra de to første<br />
prototyper af den éncylindrede<br />
9 kW motor,<br />
SM5A. Den oprindelige<br />
motors konstruktion viste sig<br />
imidlertid ikke at være særlig<br />
driftssikker, og det var<br />
nødvendigt at reparere<br />
motorerne efter adskillige<br />
havarier. For at undgå at<br />
gentage grundlæggende fejl i<br />
ud<strong>for</strong>mningen af den nye<br />
motor blev<br />
konstruktionsarbejdet<br />
udskudt i længere tid, indtil<br />
alle problemerne med de<br />
eksisterende motorer var<br />
undersøgt og løsninger på<br />
problemerne identificaret. Det<br />
medførte, at konstruktionen<br />
af den nye motor dels blev<br />
væsentligt <strong>for</strong>sinket og dels<br />
mere omfattende end<br />
oprindelig planlagt.<br />
På grund af de mange<br />
problemer med de første<br />
motorer var det kun den ene<br />
motor, SM5A-001, der blev<br />
installeret hos en anlægsvært<br />
og afprøvet over en længere<br />
periode. Den anden motor,<br />
Figur 3.1: Foto af den ene af de første motorer. Denne motor<br />
er benyttet til laboratorie<strong>for</strong>søgene.<br />
SM5A-002 blev der<strong>for</strong> installeret i laboratoriet og anvendt til at gennemføre de planlagte<br />
<strong>for</strong>søg i dette projekt i stedet <strong>for</strong> den nye motor, som ikke kunne nå at blive færdig. Figur 3.1<br />
viser den ”gamle” 9 kW motor (SM5A) monteret på en ramme med luftblæser og<br />
gasregulering. I bilag 4 er denne motor beskrevet i detaljer.<br />
Efter at motorens krumtapmekanisme og varme topstykke var repareret, gennemførtes en<br />
indledende afprøvning. Ved denne afprøvning blev det konstateret, at trykket i krumtahuset<br />
varierede langsomt med en periodetid på 10 – 15 minutter og en amplitude på 3 – 5 bar.<br />
Dette er ikke et problem under drift i en almindelig kraftvarmeinstallation, men det er<br />
uhensigtsmæssigt, når motoren skal bruges til <strong>for</strong>søg, <strong>for</strong>di det ikke er muligt at opnå en<br />
veldefineret stationær tilstand.
PSO-FU2204 24 9 kW Stirlingmotor<br />
Tabel 3.1<br />
Arbejdsgas Cylindertryktryk, max Bar 95<br />
Cylindertryktryk, min Bar 58<br />
Cylindertryk, krumtaphus Bar 72<br />
Røggas Temperatur efter lufo<br />
o<br />
C 400<br />
Temperatur før lufo 638<br />
Luft Temperatur efter lufo 511<br />
Effektivitet, lufo 0,7928<br />
Kølevand Temperatur før generator<br />
o<br />
C 30,3<br />
Temperatur før køler/motor<br />
o<br />
C 31,9<br />
Temperatur efter køler/motor<br />
o<br />
C 43,1<br />
Flow m 3 /h 1,467<br />
Cooling incl. Generator kW 21,839<br />
Cooling, engine ex. Generator kW 19,109<br />
Hat Temperatur<br />
o<br />
C 681<br />
Rør Temperatur<br />
o<br />
C 755<br />
Gas Volumenflow, start m 3<br />
1,46<br />
Volumenflow, slut m 3<br />
1,765<br />
Tid s 300<br />
Flow m 3 /h 3,66<br />
Forbrug kW 39,65<br />
Luft Volumenflow, start m 3<br />
70<br />
Volumenflow, slut m 3<br />
73,35<br />
Tid s 310<br />
Flow m 3 /h 38,903<br />
Luft/brændstof 10,629<br />
Generator Eleffekt, ekstern måler kW 9,49<br />
Cos fi, ekstern måler 0,657<br />
0<br />
Motor Virkningsgrad 0,2393<br />
Forskellige løsninger til dette<br />
problem blev diskuteret og det<br />
blev besluttet at erstatte det<br />
eksisterende aktive system til<br />
kontrol af <strong>for</strong>holdet mellem<br />
cylinderens middeltryk og<br />
krumtaphusets tryk med et<br />
passivt system bestående af en<br />
variabel læk, der kunne<br />
reguleres, således at<br />
middeltrykket ikke ændrede sig<br />
under <strong>for</strong>søgene. En<br />
afprøvning af systemet viste, at<br />
det var muligt at indstille<br />
lækken således, at trykket ikke<br />
ændrede sig måleligt, og uden<br />
at motorens effekt blev<br />
påvirket. Figur 3.2 viser<br />
blokken med Helium ventiler og<br />
manometre<br />
Der er anvendt naturgas som<br />
brændsel under <strong>for</strong>søgene.<br />
Tabel 3.1 viser de<br />
gennemsnitlige måleresultater<br />
efter en én-times testperiode i<br />
referencetilstanden.<br />
Energi<strong>for</strong>brug og virkningsgrad<br />
er baseret på nedre<br />
brændværdi af naturgassen,<br />
som er opgivet til 39 MJ/kg.<br />
Ved vurderingen af<br />
virkningsgraden skal man tage<br />
hensyn til, at luft<strong>for</strong>varmerens effektivitet kun er 75%. Luft<strong>for</strong>varmeren til den nye motor er<br />
væsentlig <strong>for</strong>bedret, og det <strong>for</strong>ventes, at virkningsgraden er ca. 1,5% point højere <strong>for</strong> den nye<br />
motor. Dataopsamlingssystemet er baseret på en PLC, der registrerer alle temperaturer og<br />
øvrige ”langsomme” signaler, der anvendes til styring og overvågning af motoren. En PC<br />
læser efterfølgende de digitaliserede måleværdier i PC-ens lager. Hurtige signaler<br />
registreres direkte i PC-en ved hjælp af det installerede A/D-kort med en maximal<br />
samplingsfrekvens på 120 kHz. Data fra målinger af kølevandsflow og temperaturer til<br />
beregning af den afsatte effekt i kølevandet overføres ved hjælp af en RS232 seriel<br />
<strong>for</strong>bindelse. Figur 3.3 viser det nye dataopsamlingssystem.
PSO-FU2204 25 9 kW Stirlingmotor<br />
Figur 3.2: Forsøgsmotoren med focus på ventilblokken med manometer og haner.<br />
Figur 3.3: Styring og dataopsamling.
PSO-FU2204 26 9 kW Stirlingmotor<br />
3.1 Resultater fra laboratorieafprøvning af SM5A<br />
Som en del af <strong>projektet</strong> er der gennemført et større måleprogram, hvor egenskaberne <strong>for</strong> 9<br />
kW-motoren, SM5A-002, er blevet kortlagt. Måleresultaterne skal blandt andet anvendes til<br />
at sammenligne med resultater fra en model til simulering af stirlingmotorens egenskaber. I<br />
løbet af måleserien er der i alt lavet 76 målinger af stationære driftstilstande. Forsøgene i<br />
DTU’s laboratorium omfatter måling af den stationære tilstand <strong>for</strong> <strong>for</strong>skellige kombinationer af<br />
følgende parametre:<br />
• middeltryk<br />
• hederrørenes temperatur<br />
• kølevandets temperatur<br />
En tabeloversigt over målingerne kan ses i bilag 5, hvor de vigtigste resultater er vist, og de<br />
følgende figurer illustrerer de målte data. Figur 3.4 viser motorens elektriske effekt som<br />
funktion af hederrørenes temperatur. Alle de på figurerne viste <strong>for</strong>søg er gennemført med<br />
Helium som arbejdsgas. Teksten på figurerne henviser til:<br />
R0: Regenerator type 0, som er standard regeneratoren fremstillet på DTU<br />
P50, P65, P80: Middeltryk, 50 bar, 65 bar, 80 bar<br />
K40, K55, K79: Kølevandstemperatur, 40°C, 55°C, 70°C<br />
Figur 3.5 viser den elektriske virkningsgrad som funktion af hederrørenes temperatur <strong>for</strong> de<br />
samme måledata, som oven<strong>for</strong>. Virkningsgraden er beregnet som <strong>for</strong>holdet mellem den<br />
producerede elektriske effekt og <strong>for</strong>bruget af naturgas beregnet på basis af den nedre<br />
brændværdi. Resultaterne viser, at den maksimale el-effekt er målt til næsten 11 kW<br />
eksklusive eget<strong>for</strong>brug til blæser og styring og den maksimale virkningsgrad til 26%.<br />
Eget<strong>for</strong>bruget udgør ca. 300 W ved fuld last.<br />
Som <strong>for</strong>ventet stiger motorens effekt og virkningsgrad som funktion af hederrørenes<br />
temperatur, mens stigende kølevandstemperatur medfører et fald i effekt og virkningsgrad.<br />
En nærmere analyse viser også, at motorens effekt er nærmest proportional med<br />
middeltrykket, mens trykket har mindre betydning <strong>for</strong> virkningsgraden.<br />
Som anført I projektansøgningen var det hensigten at gennemføre <strong>for</strong>søg med brint som<br />
arbejdsmedie. Som beskrevet i afsnittet om brint, afsnit 5.4, er <strong>for</strong>delene ved anvendelse af<br />
brint i den én-cylindrede motor begrænsede, og sikkerhedsproblemerne taget i betragtning<br />
blev det der<strong>for</strong> besluttet ikke at gennemføre brint-<strong>for</strong>søgene som planlagt. Af hensyn til<br />
sammenligning af måleresultater med resultater fra modellen til simulering af stirlingmotorers<br />
egenskaber, blev det besluttet at gennemføre en række <strong>for</strong>søg med nitrogen som<br />
arbejdsmedie. Viskositeten <strong>for</strong> nitrogen er højere end <strong>for</strong> Helium og varmeledningsevnen er<br />
lavere, hvilket medfører, at motorens egenskaber ændres væsentligt. Det er ikke nogen<br />
<strong>for</strong>del under almindelig drift, men når man alene skal anvende resultaterne til sammenligning<br />
med resultater fra simulering, er der en stor <strong>for</strong>del. Det skyldes, at de høje strømningstab og<br />
den lave varmeledningsevne stiller større krav til nøjagtigheden af de modeller, der<br />
anvendes i simuleringsprogrammet. Resultaterne med nitrogen som arbejdsmedie er<br />
nærmere beskrevet i afsnit 4.2 vedrørende simuleringsmodellen.
PSO-FU2204 27 9 kW Stirlingmotor<br />
El-power, kW<br />
12,00<br />
11,00<br />
10,00<br />
9,00<br />
8,00<br />
7,00<br />
6,00<br />
5,00<br />
Regenerator test<br />
4,00<br />
500 550 600 650 700 750<br />
Finned tube temp., °C<br />
R0-P50-K40<br />
R0-P50-K55<br />
R0-P50-K70<br />
R0-P65-K40<br />
R0-P65-K55<br />
R0-P65-K70<br />
R0-P80-K40<br />
R0-P80-K55<br />
R0-P80-K70<br />
Figur 3.4 El-effekt som funktion af hederrørendes temperatur ved trykkene (P) 50 bar, 65 bar og 80<br />
bar samt middelkølevandstemperatur på (K) 40 o , 55 o , 70 o <strong>for</strong> regeneratoren R0<br />
Efficiency (el)<br />
0,270<br />
0,250<br />
0,230<br />
0,210<br />
0,190<br />
0,170<br />
Regenerator test<br />
0,150<br />
500 550 600 650 700 750<br />
Finned tube temp., °C<br />
R0-P50-K40<br />
R0-P50-K55<br />
R0-P50-K70<br />
R0-P65-K40<br />
R0-P65-K55<br />
R0-P65-K70<br />
R0-P80-K40<br />
R0-P80-K55<br />
R0-P80-K70<br />
Figur 3.5 El-virkningsgrad som funktion af hederrørendes temperatur ved trykkene (P) 50 bar, 65 bar<br />
og 80 bar samt middelkølevandstemperatur på (K) 40 o , 55 o , 70 o <strong>for</strong> regeneratoren R0.<br />
Virkningsgraden er beregnet som <strong>for</strong>holdet mellem produceret el-effekt og gas<strong>for</strong>brug.
PSO-FU2204 28 9 kW Stirlingmotor<br />
3.2 Afprøvning af andre regeneratortyper<br />
En af målsætningerne <strong>for</strong> nærværende projekt var at afprøve nye løsninger og<br />
delkomponenter og sammenligne resultaterne med resultater fra de ”kendte” løsninger.<br />
Regeneratoren er en særdeles vigtig komponent i en stirlingmotor, <strong>for</strong>di blot små tab i<br />
regeneratoren har stor indflydelse på motorens effekt og virkningsgrad. Effektiviteten <strong>for</strong> en<br />
veldesignet regenerator i en stirlingmotor er over 99%, så kravene er store. Regeneratoren<br />
er imidlertid også en dyr komponent, og det er der<strong>for</strong> interessant at finde alternativer, der er<br />
bedre og billigere.<br />
Regeneratoren i de motorer, som er udviklet på DTU, består af tynd rustfri tråd, der er<br />
behandlet ved hjælp af en særlig proces og derefter presset til den endelige facon og<br />
porøsitet. Funktionen af regeneratorerne fremstillet på denne måde adskiller sig kun<br />
marginalt fra funktionen af regeneratorer fremstillet på konventionel måde ved at stable lag af<br />
trådnet til den ønskede højde og porøsitet. Data <strong>for</strong> regeneratoren til 9 kW motoren, der i det<br />
følgende betegnes R0, er:<br />
Indvendig diameter: 118 mm<br />
Udvendig diameter: 168 mm<br />
Længde: 61 mm<br />
Tråddiameter: 0,06 mm<br />
Porøsitet: 78 %.<br />
En ny type fiberregeneratorer udviklet af firmaet Beakert i Belgien kunne udgøre en mulighed<br />
<strong>for</strong> at reducere prisen på regeneratorerne og <strong>for</strong>bedre egenskaberne. Der<strong>for</strong> er to <strong>for</strong>skellige<br />
fiberregeneratorer fra Beakert afprøvet. De er fremstiillet ved hjælp af en særlig proces<br />
udviklet af Beakert, som er beregnet til produktion af filtre. Fibre i rustfrit stål sprøjtes og<br />
bratkøles, hvorefter de sintres i den ønskede <strong>for</strong>m. Data <strong>for</strong> de to regeneratorer, der i det<br />
følgende betegnes R1 og R2, er:<br />
R1 SF 444/65/80; d = 0,065 mm; Porosity: 80%; Mat.: SS 444<br />
R2 SF 316/22/85; d = 0,022 mm; Porosity: 85%; Mat.: SS 316<br />
Efter sintringen er regeneratorerne <strong>for</strong>sigtigt drejet i en drejebænk til de ønskede mål, så de<br />
kunne monteres i det svøb, som anvendes <strong>for</strong> at undgå blow by på regeneratorens yderside.<br />
Regeneratorerne kunne derefter monteres i 9 kW motoren på samme måde som den<br />
originale regenerator.<br />
Forsøgene med fiberregeneratoren R2 med de tynde fibre blev påvirket af, at løse fibre flød<br />
rundt i motoren med arbejdsgassen og satte sig i de PTFE-baserede stempelringe.<br />
Tætningerne virker dårligere med fibre indlejret i overfladen, hvilket indvirkede på<br />
målingerne. Desuden satte der sig fibre i den ventil, der regulerede motorens middeltryk i<br />
<strong>for</strong>hold til trykket i krumtaphuset. Problemet med løse fibre må nødvendigvis løses, inden det<br />
er muligt at anvende fiberregeneratorerne i praksis.<br />
Efter at motoren var blevet adskilt og renset, blev der gennemført en <strong>for</strong>søgsrække med<br />
begge fiberregeneratorer. Figur 3.6 og 3.7 viser en sammenligning af resultater fra <strong>for</strong>søg<br />
med den originale regenerator R0 og de to nye fiberregeneratorer R1 og R2<br />
Figur 3.6 viser den målte elektriske effekt som funktion af hederrørenes temperatur ved et<br />
middeltryk på 80 bar en kølevandstemperatur på 40 o C. Tilsvarende viser figur 3.7 den<br />
elektriske virkningsgard <strong>for</strong> det samme sæt parametre. Resultaterne viser, at<br />
virkningsgraden ikke afviger særligt meget <strong>for</strong> de 3 <strong>for</strong>skellige regeneratorer, og at <strong>for</strong>skellen<br />
på den originale regenerator R0 og fiberregeneratoren R1 med de ”grove” fibre er marginal.
PSO-FU2204 29 9 kW Stirlingmotor<br />
Derimod er el-effekten <strong>for</strong> regeneratoren R2 med de tynde fibre reduceret i <strong>for</strong>hold til de to<br />
øvrige regeneratorer.<br />
I kapitel 4 er resultaterne fra <strong>for</strong>søgene yderligere behandlet I <strong>for</strong>bindelse med<br />
sammenligninger med resultater fra simulering.<br />
Resultaterne viser, at regeneratoren R1 med de ”grove” fibre kunne være et alternativ til den<br />
originale regenerator, da det <strong>for</strong>ventes at fremstillingsprisen er væsentlig mindre. Der er dog<br />
ikke observeret nogen <strong>for</strong>bedring af motorens effekt og virkningsgrad som følge af<br />
anvendelsen af fiberregeneratorer.<br />
El-power, kW<br />
12,00<br />
11,00<br />
10,00<br />
9,00<br />
8,00<br />
7,00<br />
6,00<br />
5,00<br />
Comparison of Regenerators, Tc = 40 C<br />
4,00<br />
500 550 600 650 700 750<br />
Finned tube temp., °C<br />
R0-P80-K40<br />
R1-P80-K40<br />
R2-P80-K40<br />
Figur 3.6 El-effekt som funktion af hederrørendes temperatur <strong>for</strong> regeneratorerne R0, R1 og R2<br />
Efficiency (el)<br />
0,270<br />
0,250<br />
0,230<br />
0,210<br />
0,190<br />
0,170<br />
Comparison of regenerators, Tc = 40 C<br />
0,150<br />
500 550 600 650 700 750<br />
Finned tube temp., °C<br />
Reg. no R0<br />
Reg. no R1<br />
Reg. no R2<br />
Figur 3.7 El-virkningsgrad som funktion af hederrørendes temperatur <strong>for</strong> regeneratorerne R0, R1 og<br />
R2. Virkningsgraden er beregnet som <strong>for</strong>holdet mellem produceret el-effekt og gas<strong>for</strong>brug.
PSO-FU2204 30 9 kW Stirlingmotor<br />
3.3 Forbedring af målinger og databehandling<br />
En analyse af energibalancen <strong>for</strong> de oven<strong>for</strong> beskrevne målinger var ikke tilfredsstillende.<br />
Ved høje temperaturer og tryk og dermed ved høj motoreffekt var der ca. 3 kW, der ikke<br />
kunne gøres rede <strong>for</strong>, hvilket svarede til ca. 7% af den indfyrede effekt. Det blev der<strong>for</strong><br />
besluttet at <strong>for</strong>søge at reducere de tab, der ikke kunne gøres rede <strong>for</strong>.<br />
Varmetabet fra det uisolerede<br />
krumtaphus var et af de steder, der blev<br />
vurderet til at udgøre en væsentlig<br />
fejlkilde. Der<strong>for</strong> blev hele krumtaphuset<br />
isoleret med ”Airflex” isoleringsplader,<br />
der består af en tæt cellestruktur af<br />
kunststof. For at undgå overhedning af<br />
krumtaphuset, der normalt afgiver<br />
varmen fra tabene i mekanismen direkte<br />
gennem krumtaphusets væg til<br />
omgivelserne, var det nødvendigt at<br />
fremstille en van<strong>dk</strong>appe til montering på<br />
krumtaphusets dæksel. Van<strong>dk</strong>appen<br />
blev tilsluttet kølevandssystemet på en<br />
måde, så det var muligt at måle de<br />
Figur 3.8. Isolering af krumtaphus med van<strong>dk</strong>ølet<br />
dæksel<br />
enkelte varmetab individuelt. Figur 3.8<br />
viser den isolerede motor med de to<br />
isolerede rørtilslutninger til<br />
dækselkølingen.<br />
Det varme røggasrør er isoleret med en mineralisoleringsmåtte. Der er indsat en skive med<br />
vinger i røggaskanalen <strong>for</strong> at mikse strømningen, således at temperaturen er den samme<br />
overalt i tværsnittet ved målingen af røggastemperaturen. Dermed er det muligt at bestemme<br />
en mere korrekt energistrøm i den røggas, der <strong>for</strong>lader motoren.<br />
Som vist på figur 3.9 er der installeret et nøjagtigt manometer, der ved hjælp af nogle ventiler<br />
kan måle både det maksimale og det minimale tryk i cylindervolumnet. Trykket i<br />
krumtaphuset kan aflæses på et separat manometer.<br />
Figur 3.9: Præcisionsmanometer med to tilslutninger til<br />
registrering af det maksimale og det minimale<br />
cylindertryk.
PSO-FU2204 31 9 kW Stirlingmotor<br />
Efter ombygninger af testfaciliteterne er der gennemført et måleprogram, hvor<br />
fiberregeneratoren R2 var monteret i motoren. Som ved de tidligere målinger er motorens eleffekt<br />
samt gas<strong>for</strong>brug og varmeeffekten afsat i kølesystemet målt som funktion af<br />
driftstemperaturer og middeltryk. Det <strong>for</strong>bavsende var, at isoleringen af motoren ikke havde<br />
<strong>for</strong>bedret varmebalancen nævneværdigt. Der manglede stadig ca. 3 kW i den totale<br />
energibalance ved fuld last, hvilket indikerer at den anvendte metode ved de tidligere <strong>for</strong>søg<br />
til at estimere varmetabet fra krumtaphuset må have været god. Fejlen må der<strong>for</strong> komme fra<br />
en eller flere af de andre målinger, men det har ikke været muligt at undersøge dette<br />
inden<strong>for</strong> <strong>projektet</strong>s rammer.<br />
I praksis er det ikke muligt at variere en enkelt parameter under <strong>for</strong>søgene mens de andre<br />
parameter er fuldstændigt uændrede. Ændres en temperatur eller et tryk ændres de øvrige<br />
driftsparametre også lidt, inden der igen er opnået en stationær tilstand. Det gør det<br />
vanskeligt at sammenligne resultater fra <strong>for</strong>skellige målinger og fra beregninger.<br />
For bedre at kunne præsentere måleresultaterne og sammenligne dem med resultater fra<br />
f.eks. simuleringer af processen er de målte data tilnærmet med følgende polynomium:<br />
F = A1 TH 2 + A2 TC 2 + A3 P 2 + A4 TH TC + A5 TH P + A6 TC P + A7 TH + A8 TC + A9 P + A10<br />
hvor F Funktionsværdi (effekt, varme til kølevand etc.)<br />
TH Hederrørenes temperatur<br />
TC Kølevandets temperatur<br />
P Arbejdsgassens middeltryk<br />
A1 – A10 er konstanter, som findes ved at minimere residual-funktionerne, der indeholder alle<br />
målepunkterne. Figur 3.10 viser el-virkningsgraden beregnet på basis af målte data <strong>for</strong><br />
motoren monteret med fiberregeneratoren R2.<br />
Electric efficiency [-]<br />
0,25<br />
0,24<br />
0,23<br />
0,22<br />
0,21<br />
0,2<br />
0,19<br />
0,18<br />
0,17<br />
0,16<br />
0,15<br />
450 500 550 600 650 700 750<br />
Heater temperature [°C]<br />
Tcold=38°C<br />
Tcold=45°C<br />
Tcold=52°C<br />
Tcold=60°C<br />
Figur 3.10. Brutto elvirkningsgrad (ex. El-<strong>for</strong>brug <strong>for</strong> blæser og styring) som funktion af<br />
hederrørstemperaturen ved <strong>for</strong>skellige kølevandstemperaturer præsenteret ved hjælp af tilpassede<br />
polynomier.
PSO-FU2204 32 9 kW Stirlingmotor<br />
3.4 Afprøvning af 9 kW motor på Fachhochschule Reutlingen<br />
Der en betydelig international opmærksomhed om de danske aktiviteter inden<strong>for</strong><br />
stirlingmotorer til biomasse. En mangeårig kontakt til Prof. Bernd Thomas, Fachhochschule<br />
Reutlingen, førte til, at en af 9 kW motorerne blev udlånt til laboratorie<strong>for</strong>søg. Det var særligt<br />
interessant, <strong>for</strong>di der er opbygget en omfattende prøvestand i Reutlingen, hvor<br />
omdrejningstallet også kan varieres.<br />
Figur 3.11 viser motorens el-effekt, afsat varme i kølesystemet og gas<strong>for</strong>brug som funktion af<br />
motorens omdrejningstal ved nominele tryk og temperatur. Tilsvarende viser figur 3.12<br />
motorens el-virkningsgrad samt total energiudnyttelse (summen af el- og varmeproduktion i<br />
<strong>for</strong>hold til gas<strong>for</strong>brug) også som funktion af omdrejningstallet ved nominelle tryk og<br />
temperaturer. Som det ses af figurerne 3.11 og 3.12 vokser effekt og varmeproduktion som<br />
<strong>for</strong>ventet nærmest proportionalt med omdrejningstallet, men motorens el-virkningsgrad har et<br />
maksimum ved ca. 800 rpm. En total energiudnyttelse på 85% er tilfredsstillende, selvom<br />
dette kunne være bedre. Man skal her være opmærksom på, at motorens brændersystem er<br />
designet til biogas, hvor røggstemperaturen ikke må være mindre end 130 o C <strong>for</strong> at undgå<br />
problemer med kondensation af svovl<strong>for</strong>bindelser. Der er dermed en øvre grænse <strong>for</strong>, hvor<br />
høj en total energiudnyttelse kan blive.<br />
.<br />
P / Q<br />
[KW]<br />
40,0<br />
35,0<br />
30,0<br />
25,0<br />
20,0<br />
15,0<br />
10,0<br />
5,0<br />
0,0<br />
.<br />
.<br />
QQ_Brenngas Q_BHKW P_el<br />
Brenngas<br />
QBHKW Pel 350 400 450 500 550 600 650 700 750 800 850 900 950 1000 1050<br />
Drehzahl [Upm]<br />
Figur 3.11: Motorens el-effekt, afsat varme i kølesystemet og gas<strong>for</strong>brug som funktion af motorens<br />
omdrejningstal ved nominele tryk og temperatur.
PSO-FU2204 33 9 kW Stirlingmotor<br />
ηh<br />
[%]<br />
100,0<br />
90,0<br />
80,0<br />
70,0<br />
60,0<br />
50,0<br />
40,0<br />
30,0<br />
20,0<br />
10,0<br />
0,0<br />
eta_el η el<br />
η eta_ges ges<br />
350 400 450 500 550 600 650 700 750 800 850 900 950 1000 1050<br />
Drehzahl [Upm]<br />
Figur 3.12: Motorens el-virkningsgrad samt total energiudnyttelse som funktion af omdrejningstallet<br />
ved nominelle tryk og temperaturer.<br />
3.5 Maskine til afprøvning af regeneratorer<br />
Som tidligere nævnt er regeneratoren en særdeles vigtig komponent, når man designer en<br />
stirlingmotor. Det er der<strong>for</strong> også særdeles vigtigt at modeller til simulering af stirlingmotorers<br />
egenskaber indeholder en god analytisk model af trykfald og varmetransmission i<br />
regeneratoren. Særligt er varmetransmission i regeneratoren vanskeligt at måle, <strong>for</strong>di<br />
varmetransporten ikke kan måles under stationære <strong>for</strong>hold. Regeneratoren er et varmelager,<br />
og det er kun muligt at opnå en periodisk stationær tilstand.<br />
Som oven<strong>for</strong> beskrevet kan man afprøve <strong>for</strong>skellige regeneratorer i motoren og måle en<br />
eventuel <strong>for</strong>skel fra en regenerator til en anden, men det er ikke muligt at isolere<br />
tabsmekanismer, så man kan verificere algoritmer til beregning af regeneratorens tryktab og<br />
varmetransmission. Det er systemet <strong>for</strong> sammenhængende til, og desuden har det<br />
varierende tryk i motoren også stor indflydelse på de målinger, man kan lave.<br />
Med henblik på at undersøge de i litteraturen kendte algoritmer til beregning af<br />
varmetransmision og trykfald i regeneratorerne blev det der<strong>for</strong> besluttet at opbygge en<br />
prøvestand baseret på en <strong>for</strong>søgs-stirlingmotor, der blev fremstillet i perioden 1992 – 1993.<br />
En stirlingmotor består af to variable volumener, kompressionsvolumenet og<br />
ekspansionsvolumenet, som er <strong>for</strong>bundet gennem tre varmevekslere, køleren (den kolde<br />
varmeveksler), regeneratoren og hederen (den varme varmeveksler). Systemets egenskaber<br />
kan måles, men det er meget svært at isolere måledata <strong>for</strong> regeneratoren, bl.a. <strong>for</strong>di<br />
kompressionsvarmen og ekspansionsvarmen fra cylindervolumenerne som følge af<br />
trykvariationerne i motoren er meget større end varmetabene fra regeneratoren.<br />
For at komme uden om dette problem er regeneratormaskinen ud<strong>for</strong>met på en måde, så et<br />
stempel i en cylinder flytter arbejdsgassen frem og tilbage gennem regeneratoren med<br />
tilnærmelsesvis konstant tryk. Dette er opnået ved at anbringe to regeneratorer og to kølere<br />
symmetrisk på hver sin side af en elektrisk opvarmet varmeveksler som vist på figur 3.13 og<br />
3.14. Resultatet er, at to variable volumener med tilsammen konstant masse og temperatur<br />
er <strong>for</strong>bundet gennem fem varmevekslere, således at temperaturen af gassen er nærmest<br />
uændret under drift overalt i maskinen. Gassen flyttes mellem det kolde volumen nr. 1 over
PSO-FU2204 34 9 kW Stirlingmotor<br />
stemplet og det kolde volumen nr. 2 under stemplet, men temperaturen af disse variable<br />
volumener ændres ikke uanset temperatur<strong>for</strong>skellen mellem kølevand og heder. Yderligere<br />
er volumenet under stemplet <strong>for</strong>bundet til krumtaphusets volumen, som er meget stort i<br />
<strong>for</strong>hold til cylindervolumenet, hvilket medfører at trykket i dette volumen er konstant under<br />
afprøvningen. Tabet i regeneratoren kan måles ved at måle den varme, der afsættes i<br />
køleren umiddelbart efter regeneratoren.<br />
I en stirlingmotor er det flere varmetab, der sammen med regeneratortabet medvirker til at<br />
reducere motorens virkningsgrad. Et af dem er ”shuttle” tabet, der er en følge af, at<br />
<strong>for</strong>trængerstemplets væg bevæger sig frem og tilbage langs en cylinder med en stor<br />
temperaturgradient. Et givet punkt på <strong>for</strong>trængerens cylindriske del er skiftevis placeret ud<br />
<strong>for</strong> et varmere og et koldere sted på cylinderen, hvorved der transporteres varme via<br />
<strong>for</strong>trængeren fra den varme ende til den kolde. I regeneratortestmaskinen er dødvolumen<br />
snarere en <strong>for</strong>del, så trykket kan holdes konstant, og der<strong>for</strong> består stemplet, der flytter<br />
gassen mellem de to cylindervolumener, i praksis af to stempler <strong>for</strong>bundet med en stang.<br />
Rundt om stangen er cylinderen isoleret, så varmetab ved ledning reduceres og ”shuttle” tab<br />
elimineres.<br />
Tabet fra de to regeneratorer afsættes som nævnt i de to kølere, som sidder umiddelbart i<br />
<strong>for</strong>længelse af regeneratorerne. Det er imidlertid ikke kun regeneratortabet, der afsættes i<br />
disse kølere. Også varmeledning gennem yder- og indercylinder fra den varme til den kolde<br />
del ender i køleren, og det er der<strong>for</strong> vigtigt at minimere dette varmetab. I<br />
regeneratortestmaskinen er ydercylinderen der<strong>for</strong> lavet som en tynd væg, der ikke er<br />
trykbærende, og som er isoleret fra en ydre kold væg, der er trykbærende. Denne metode til<br />
at reducere varmetabet kan kun lade sig gøre ved elektrisk varmetilførsel, og den kan ikke<br />
anvendes i en motor med almindelig brænder.<br />
Figur 3.13 Skematisk fremstilling af regeneratortestmaskinen.
PSO-FU2204 35 9 kW Stirlingmotor<br />
Figur 3.14 Tegning af komplet test maskine <strong>for</strong> regeneratorer
PSO-FU2204 36 9 kW Stirlingmotor<br />
Figur 3.14 viser en snittegning af regeneratortestmaskinen og figur 3.15 viser et foto.<br />
Testmaskinen er designet til et maksimalt driftstryk på 80 bar. Regeneratorerne, der<br />
anvendes i testmaskinen, er de samme som de regeneratorer, der anvendes i 9 kW motoren,<br />
og det er dermed muligt at undersøge 9 kW motorens regeneratorer under de samme <strong>for</strong>hold<br />
med hensyn til arbejdsgas, tryk og temperatur, som <strong>for</strong>ekommer i motoren.<br />
Figur 3.15 Samling af regeneratortestmaskinen. Den røde del er fra en eksisterende hermetisk<br />
stirlingmotor med tryksat generator.<br />
Figur 3.16 viser et foto af den elektriske heder til regeneratortestmaskinen. Den består af 9<br />
steatitter med seks huller med spiralviklet modstnadstråd i hver. Den elektriske heder var<br />
designet til en effekt på 4 kW ved en spænding på 400 Volt og en Heliumtemperatur på 600 o .<br />
Desværre viste det sig, at kølingen af modstandstråden på den elektriske heder ikke var<br />
tilstrækkelig på de steder, hvor strømningshastigheden var lav. Når belastningen kom over<br />
ca. 3 kW brændte tråden over i løbet af nogle minutters drift, og <strong>for</strong>søgene måtte afbrydes.
PSO-FU2204 37 9 kW Stirlingmotor<br />
Figur 3.16 Elektrisk varmeelement til regenerator test maskine.<br />
Først blev den øverste regenerator nr. 1 længst fra krumtaphuset valgt som den regenerator,<br />
som blev benyttet til detaljerede målinger. Dermed var det nemmere at skifte regenerator,<br />
når en ny regenerator skulle monteres. Der er imidlertid et betydeligt tryktab over de to<br />
regeneratorer, som fører til et pdV-arbejde i det øverste cylindervolumen nr.1. Dette pdVarbejde<br />
ender som varme i køleren, der også benyttes til at måle tabet fra regeneratoren.<br />
Ved 1000 rpm var tabet ca. 2,5 kW, hvilket er ca. 3 gange større end regeneratortabet. Det<br />
var der<strong>for</strong> nødvendigt at ændre opstillingen.<br />
Hele cylinderen med varmevekslere, regeneratorer og heder blev der<strong>for</strong> vendt på hovedet,<br />
så måle-regeneratoren nu sad nærmest cylindervolumen nr. 2 ved krumtaphuset. Et andet<br />
problem var, at den elektriske heder er konstrueret på basis af hovedmålene fra hederen på<br />
9 kW stirlingmotoren, da det ikke er muligt at behandle en maskine med 5 varmevekslere<br />
ved hjælp af den eksisterende simuleringsmodel. I modsætning til en almindelig stirlingmotor<br />
er massestrømmen imidlertid den samme i alle 5 varmevekslere, hvilket medfører at<br />
volumenstrømmen i testmaskinen elektriske heder er ca. 3 gange større end i en<br />
stilringmotor med det samme cylindervolumen. Trykfaldet over hederen var der<strong>for</strong> også<br />
større end <strong>for</strong>udsat.<br />
Der blev lavet I alt 11 stationære målinger med regeneratortestmaskinen, inden den<br />
elektriske heder gav op. Figur 3.17 viser varmetabet som funktion af den gennemsnitlige<br />
temperatur i gassen i hederen <strong>for</strong> tre <strong>for</strong>skellige tryk, 30 bar, 50 bar og 70 bar. Målingerne er<br />
vanskelige at gennemføre, <strong>for</strong>di varmetabene i regeneratoren er små i <strong>for</strong>hold til de termiske<br />
masser i kølere, cylindre m.m. Det tager der<strong>for</strong> meget lang tid (ca. 3 timer) før en stationær<br />
tilstand er opnået. Desuden er det nødvendigt at korrigere <strong>for</strong> varmetabene i cylindervægge.<br />
Dette gøres ved at stoppe maskinen umiddelbart efter, at målingen er registreret. Når<br />
maskinen står stille, indstiller der sig i løbet af de næste 1 – 2 minutter en ny tilstand, hvor<br />
det kun er de stationære tab tab tilbage. Disse registreres også og trækkes fra værdierne,<br />
som er registreret mens maskinen kørte.
PSO-FU2204 38 9 kW Stirlingmotor<br />
Q [kJ/s]<br />
1,20<br />
1,00<br />
0,80<br />
0,60<br />
0,40<br />
0,20<br />
0,00<br />
Korrigeret regeneratortab<br />
100 150 200 250 300 350 400<br />
Middel hedertemperatur, C<br />
50 bar<br />
30 bar<br />
75 bar<br />
Figur 3.17 Resultater fra testmaskine med regenerator R2. Testen viser korrigeret varmetab som<br />
function af heder temperatur.<br />
Q_tab, W<br />
2000<br />
1800<br />
1600<br />
1400<br />
1200<br />
1000<br />
800<br />
600<br />
400<br />
200<br />
0<br />
Samlet regeneratortab<br />
0 100 200 300 400<br />
Heder middeltemperatur, C<br />
Kühl<br />
ILK Dresden<br />
G & W felts original<br />
G & W felts (v. T & P)<br />
G & W screens (original)<br />
G & W screens (v. T & P)<br />
Tong/London<br />
Tanaka<br />
Miyabe<br />
Målinger<br />
Figur 3.18. Beregnede og målte vværdier <strong>for</strong> regeneratortabet som funktion af<br />
middeltemperaturen i hederen ved et middeltryk på 50 bar Helium.
PSO-FU2204 39 9 kW Stirlingmotor<br />
Formålet var at vurdere de <strong>for</strong>skellige algoritmer til beregning af Nusselts tal <strong>for</strong><br />
regeneratorer. På figur 3.18 er måleresultaterne fra et tryk på 50 bar og Helium som<br />
arbejdsgas sammenlignet med 9 <strong>for</strong>skellige algoritmer <strong>for</strong> Nusselt tallet kendt fra litteraturen.<br />
Målingerne er vist med mørkeblåt markeret med firkanter på spidsen. Som det ses af figuren,<br />
er der en betydelig spredning mellem resultater fra beregning af varmetabet med de<br />
<strong>for</strong>skellige korellationer <strong>for</strong> Nusselt tallet. Desuden er det målte regeneratortab lidt mindre<br />
end de laveste beregnede værdier. Resultaterne bør undersøges nærmere, og der er der<strong>for</strong><br />
ikke indført ændringer i simuleringsmodellen på basis af de målte værdier <strong>for</strong> Nusselt tallet i<br />
regeneratoren.<br />
Der er også gennemført en måleserie af tryktabet over 9 kW motorens regeneratorer under<br />
stationære <strong>for</strong>hold. Ved hjælp af en ringkammerblæser er luft blæst gennem regeneratorerne<br />
ved relevante Reynolds tal, og flow og tryktab er registreret. Figur 3.19 viser et eksempel,<br />
hvor måledata <strong>for</strong> den originale regenerator, R0 (tråddiam.=0,06 mm, porøsitet=78%) er<br />
sammenlignet med <strong>for</strong>skellige empiriske udtryk kendt fra samme kilder i litteraturen som<br />
oven<strong>for</strong> ved undersøgelsen af Nusselt tallet <strong>for</strong> varmeovergang. Her viser resultaterne, at<br />
måledata ligger i midten tæt på de empiriske udtryk <strong>for</strong>eslået af Blas og ILK Dresden.<br />
Algoritmerne i simuleringsprogrammet til beregning af tryktab i regeneratorerne er der<strong>for</strong><br />
ændret i henhold til måleresultaterne.<br />
Trykfald, Pa<br />
30000<br />
25000<br />
20000<br />
15000<br />
10000<br />
5000<br />
0<br />
Trykfald ved stationær tilstand, tråd 304/60/78<br />
G & W Felt<br />
G & W Screens<br />
Kühl<br />
ILK Dresden<br />
Måledata<br />
Blas<br />
Tanaka<br />
Tong & London<br />
Miyabe<br />
0 10 20 30 40 50 60<br />
Reynolds tal<br />
Figur 3.19: trykfald over regeneratorer.
PSO-FU2204 40 9 kW Stirlingmotor<br />
4 Model til simulering af stirlingmotorer<br />
Det Numeriske SimuleringsProgram (NSP), der hidtil er anvendt ved design og optimering af<br />
de oven<strong>for</strong> beskrevne motorer, er udviklet på Inst. <strong>for</strong> Energiteknik, DTU. I modellen er<br />
stirlingmotorens hove<strong>dk</strong>omponenter opdelt i fem delvolumener, ekspansionsvolumen,<br />
hedervolumen, regeneratorvolumen, kølervolumen og kompressionsvolumen. I heder- og<br />
kølervolumenerne er der yderligere inkluderet eventuelle <strong>for</strong>bindelsesvolumener, som<br />
omfatter de rør eller kanaler, der <strong>for</strong>binder cylindervolumenet med henholdsvis<br />
køleren/hederen og regeneratoren. Figur 4.1 viser diskretiseringen af den simplificerede<br />
motorproces.<br />
Sym. akse<br />
Kompresions<br />
volumen<br />
Tilslutnings vol.<br />
Køler<br />
Regenerator<br />
Spalte afstand<br />
Stempel<br />
Heder- og kølervolumen er yderligere delt i 10 – 20 delvolumener (kan specificeres efter<br />
behov), mens regeneratoren er beskrevet som et enkelt volumen med en<br />
temperaturgradient, der varierer med tiden. De termodynamiske ligninger (energiligningen,<br />
kontinuitet m.m.) <strong>for</strong> hvert volumen danner til sammen et sæt af differentialalgebraiske<br />
ligninger, således at antallet af uafhængige variable i alt udgør 5 + 2 *N, hvor N er antallet af<br />
delvolumener i de to varmevekslere (køler og heder). Dette ligningssystem omfatter også<br />
empiriske ligninger <strong>for</strong> beregning af Nusseltstallet som funktion af Reynoldstallet <strong>for</strong><br />
varmetransmissionen i alle 5 delvolumener. I den eksisterende model er disse ligninger løst<br />
ved at antage, at trykket i hele maskinen til et givet tidspunkt er det samme, og ved samtidig<br />
at linearisere ligningerne i en omegn om det aktuelle punkt. Dermed omdannes<br />
ligningssystemet til et sæt sammenhørende lineære differentialligninger, som løses direkte<br />
som et begyndelsesværdiproblem ved hjælp af en standardmetode (f.eks. 4. ordens Runge<br />
Kutta).<br />
En stationær driftstilstand <strong>for</strong> en maskine med periodisk varierende <strong>for</strong>hold i<br />
arbejdsvolumenerne er karakteriseret ved, at startværdierne <strong>for</strong> tryk og temperaturer i alle<br />
delvolumener er identiske med slutværdierne efter en hel omdrejning.<br />
Ligningssystemet løses ved at gætte et sæt startværdier <strong>for</strong> tryk og temperaturer og<br />
sammenligne de gættede værdier med slutværdierne efter en omdrejning. Herefter<br />
korrigeres startværdierne, og processen gentages indtil <strong>for</strong>skellen mellem start- og<br />
Heder<br />
Figur 4.1 Diskretiseringen af den simplificerede motorproces<br />
Eksp. volumen<br />
Spalte dybde
PSO-FU2204 41 9 kW Stirlingmotor<br />
slutværdier er tilstrækkelig lille. I den aktuelle simuleringsmodel kræves typisk mellem 12 og<br />
24 iterationer <strong>for</strong> at finde en stationær løsning, hvilket på en moderne PC med en clockfrekvens<br />
på 2,5 MHz tager mellem 3 og 6 sekunder..<br />
Tabel 4.2 viser en liste over de 41 input parametre, der beskriver motorens konstruktion. Se<br />
også figur 4.4 der viser parametrene på en tegning af motoren. Tabel 4.3 viser et eksempel<br />
på resultatet fra en simulering, hvor data fra SM5B-motorens er anvendt. Der, hvor der er<br />
angivet to værdier, beskriver den første værdi relevante data <strong>for</strong> ekspansionsvolumenet og<br />
den anden værdi relevante data <strong>for</strong> kompressionsvolumenet.<br />
4.1 Modifikation af eksisterende simuleringsmodel<br />
Fordelen ved den skitserede model er, at den er meget hurtig, <strong>for</strong>di ligningerne kan opstilles<br />
på eksplicit <strong>for</strong>m. Det har stor betydning, når modellen benyttes sammen med programmer til<br />
optimering. Antagelsen om at trykket i hele motoren til et givet tidspunkt er det samme, er<br />
imidlertid ikke velegnet til at behandle motorer, hvor trykfaldet i f.eks. hederen er stort. Det er<br />
netop tilfældet i motorer til biomasse, hvor hederen nødvendigvis skal have en stor udvendig<br />
overflade <strong>for</strong> at gøre det muligt at overføre varmen fra røggasserne til arbejdsgassen.<br />
Der blev der<strong>for</strong> lavet en vurdering af omfanget af arbejdet <strong>for</strong>bundet med at inkludere trykket<br />
i hvert delvolumen i det samlede sæt af differentialligninger. Det blev imidlertid hurtigt klart, at<br />
det resulterende system af differentialalgebraiske ligninger dermed ikke kunne løses<br />
eksplicit, hvilket kolliderede med hele filosofien bag det eksisterende<br />
simuleringsprogram. Dermed ville det være nødvendigt at starte helt <strong>for</strong>fra og udvikle et nyt<br />
simuleringsprogram baseret på iterativ løsning af det differentialalgebraiske ligningssystem.<br />
Dette arbejde vurderedes til alene at kræve ca. to mandeår, og det var der<strong>for</strong> ikke muligt<br />
inden<strong>for</strong> rammerne af nærværende projekt.<br />
Arbejdet blev i stedet koncentreret om, hvordan beregningen af det lokale tryk i motorens<br />
delvolumener kunne <strong>for</strong>bedres i den eksisterende simuleringsmodel uden at ændre den<br />
grundlæggende løsningsmetode. Dette kompromis var særligt attraktivt, <strong>for</strong>di det blev muligt<br />
at sætte et nyt Ph.D. projekt i gang samtidig med dette projekt, der fokuserede på iterativ<br />
løsning af ligningssystemet til beskrivelse af <strong>for</strong>holdene i stirlingmotorens delvolumener<br />
inklusive trykfald og bevægelsesmængde i de enkelte delvolumener. Resultatet fra denne<br />
simuleringsmodel kan efterfølgende sammenlignes med resultater fra den modificerede<br />
udgave af det eksisterende simuleringsprogram.<br />
Det eksisterende simuleringsprogram blev der<strong>for</strong> ændret således, at trykket i midten af<br />
regeneratoren blev valgt som referencetryk. Trykket i motorens cylindervolumener kan<br />
dermed udtrykkes ved:<br />
pv = pR + ¡ pv, i , i = 1 – n<br />
Hvor pv = trykket i ekspansions- eller kompressionsvolumenet<br />
pR = trykket i midten af regeneratoren<br />
¡ pv, i = trykfaldet over delvolumen i<br />
n = antal delvolumener mellem regeneratorens midte og cylinderen
PSO-FU2204 42 9 kW Stirlingmotor<br />
Tabel 4.2: Indput til simuleringsprogram.<br />
Stirling engine simulation program<br />
Version 2.1<br />
Filename: SM5B.DAT<br />
Date: 2005-10-26<br />
BO Bore, mm 114 114<br />
S Stroke, mm 54 54<br />
PS Piston rod, mm 12<br />
PH Phase angle, deg 79<br />
CON Configuration 2 (beta-type)<br />
PMID Mean pressure, MPa 8<br />
N0 Number of rev., rpm 1020<br />
JG Working gas 2 (Helium)<br />
MS Piston mass, kg 5 5<br />
GR Displacer gap, mm 0.8<br />
TS Displ. wall thickness, mm 2.5<br />
MATF Displacer material 2 (Stainless steel)<br />
DS Piston clearance, mm 2 2<br />
F0 Piston friction, N 80 150<br />
NYM Mechanical efficiency 0.9<br />
LPC Displ. cyl. lenght, mm 35<br />
TPC Displ. cyl. thickness, mm 1.5 0<br />
MAT Displ. cyl. Material 2 (Stainless steel)<br />
LRC Reg. cyl. lenght, mm 35<br />
TRC Reg. cyl. thickness, mm 15 15<br />
MATS Reg. cylinder mat. 4 (Heat resistant stainless steel)<br />
TRY Regenerator with, mm 25.5<br />
DIR Reg. inner diameter, mm -1<br />
LR Regenerator length,mm 61<br />
FF Filler factor 0.22<br />
TT Tread diameter, mm 0.06<br />
MATR Regenerator material 2 (Stainless steel)<br />
T0 Temperature, C 720 60<br />
VDA Dead volume, cm**3 10 58<br />
NV Number of tubes 32 261<br />
LV Exchanger tube length, mm 395 132<br />
DLV Inactive tube lenght, mm 180 8<br />
D Exch. tube diameter, mm 5 2.5<br />
TLV Tube wall thickness, mm 2.5 0.4<br />
FVT Entrance losses 3 3<br />
MATV Exchanger tube material 4 2<br />
CH Heattransfer const. 0 1.0E-4<br />
DVM Manifold tube diameter, m 50 45<br />
LVR Reg. manifold lenght, mm 60 14<br />
LVC Cyl. manifold lenght, mm 0 25<br />
DAH Heatconduction const. 0.5 0.5
PSO-FU2204 43 9 kW Stirlingmotor<br />
Tabel 4.3: Eksempel på output fra simuleringsprogram.<br />
Expansion Compression<br />
Ind. work, W 29234 -14159<br />
Regenerator loss, W 1507.9<br />
Shuttleconduction, W 314.9<br />
Pumping loss, W 661.3<br />
Heatconduction, W 1392.1<br />
Heat I/O, W 32139 -20532<br />
Ind. poweroutput, W 15075<br />
Regen. powerloss, W -466.6 -276.7<br />
Heatex. powerloss, W -491.1 -224.3<br />
Manifold powerloss, W -20.6 -86.0<br />
Pistonring loss, W -146.9 -275.4<br />
Mechanical loss, W -1507.5<br />
Poweroutput, W 11580<br />
Efficiency/COP 0.360<br />
Disp. cyl./wall, mm 1.5 1.5<br />
Rege. cyl. thn., mm 15.0 15.0<br />
Pressure, max/min 8.00 1.72<br />
Max pres. & ampl., MP 10.12 4.23<br />
Bearing <strong>for</strong>ce, N 1672.6 15353.4<br />
Det er alene korrektionen af trykket i de to variable cylindervolumener, der er interessante,<br />
<strong>for</strong>di det er her, at effekten afsættes. Arbejdet (med <strong>for</strong>tegn) produceret i de to variable<br />
cylindervolumener kan dermed udtrykkes ved:<br />
Wv = (pR + ¡ pv, i, i = 1 – n) dVv = pR dVv + (¡ pv, i, i = 1 – n) dVv<br />
Hermed adskilles pdV arbejdet fra processen uden tryktab, og tryktabets indflydelse <strong>for</strong> de<br />
enkelte komponenter kan adderes efterfølgende. Det er samtidig muligt at identificere de<br />
steder i maskinen, der bidrager mest til tab af effekt på grund af trykfald.<br />
Ud over de oven<strong>for</strong> beskrevne ændringer er de enkelte delmodeller til beregning af<br />
varmetransmission og trykfald i regeneratoren revurderet. Dette er nærmere behandlet i<br />
afsnittet om regeneratortestmaskinen (afsnit 3.5). Desuden er de empiriske udtryk kendt fra<br />
litteraturen til beregning af spaltetabet undersøgt nærmere ved hjælp af den avancerede<br />
stirlingmotor simuleringsmodel. Spaltetabet opstår i spalten mellem <strong>for</strong>trængeren og<br />
cylindervæggen. Som vist på figur 4.5 er spalten <strong>for</strong> oven åben til det varme<br />
ekspansionsvolumen, mens spalten <strong>for</strong> neden ender ved <strong>for</strong>trængerens stempeltætning.<br />
Temperaturen stiger fra stempeltætningen, hvor temperaturen i tæt på kølevandets<br />
temperatur på 50 o C – 80 o C, til ekspansionvolumenet, hvor temperaturen under drift er i<br />
nærheden af 700 o C. Både cylindervæggen og den cylindriske del af <strong>for</strong>trængerens væg har<br />
dermed en temperaturgradient, hvilket medfører, at når <strong>for</strong>trængeren er i den øverste stilling<br />
ser et givet punkt på <strong>for</strong>trængerstemplet en varmere del af cylinderen, end når <strong>for</strong>trængeren<br />
er i bunddødpunkt. Fortrængerens vægtemperatur vil der<strong>for</strong> indstille sig i en<br />
ligevægtstilstand, således at der overføres varme fra væggen til <strong>for</strong>trængeren, når cylinderen<br />
befinder sig over sin middelposition, mens der overføres varme den anden vej fra<br />
<strong>for</strong>trængervæg til cylindervæg, når cylinderen befinder sig under sin middelposition. På den<br />
måde transporteres der varme fra den varme ende til den kolde ende som følge af<br />
stempelbevægelsen og temperaturgradienten. Dette tab kaldes også ”shuttle conduction”.
PSO-FU2204 44 9 kW Stirlingmotor<br />
Figur 4.4: Illustrering af de fysiske størrelser der benyttes i simuleringsprogrammet.
PSO-FU2204 45 9 kW Stirlingmotor<br />
Der er endnu et tab, der er <strong>for</strong>bundet til spalten mellem <strong>for</strong>trænger og cylinder. Når trykket<br />
varierer i motoren, vil der strømme arbejdsgas (Helium) ind i spalten, når trykket stiger, og ud<br />
af spalten, når trykket falder. Den faldende temperatur fra indløbet til spalten medfører, at<br />
den indstrømmende gas afkøles, og den opvarmes igen, når den strømmer tilbage.<br />
Resultatet er en varmetransport, der ligesom shuttletabet afhænger af spaltens længde og<br />
bredde og temperaturgradienten. Dette tab omtales ofte som spaltepumpetabet.<br />
Mens spaltepumpetabet vokser med spalteafstanden, falder shuttleconduction-tabet med<br />
spalteafstanden. Det medfører, at der er et optimum <strong>for</strong> spalteafstanden, hvor også<br />
dødvolumenet i spalten har betydning <strong>for</strong> motorens effekt og virkningsgrad.<br />
Det oven<strong>for</strong> omtale program til simulering af stirlingmotorer er som nævnt velegnet til<br />
simulering og optimering af motorernes konstruktion, med det egner sig ikke til detaljerede<br />
studier af særlige fænomener som f.eks. spaltetabene. Et andet simuleringsprogram, som er<br />
udviklet i <strong>for</strong>bindelse med et Ph.D. projekt, er blevet udviklet med henblik på detaljerede<br />
undersøgelser af specifikke fænomener i stirlingmotorens <strong>for</strong>skellige volumener. Dette<br />
program kræver ca. 100 gange så lang regnetid, og det er dermed ikke operationelt på<br />
samme måde som det eksisterende program, hvis man vil lave en interaktiv deloptimering af<br />
en given konstruktion. De to simulerings-modeller supplerer der<strong>for</strong> hinanden på bedste<br />
måde.<br />
Indløb til<br />
spalte<br />
Figur 4.5: Spaltpumpetab og<br />
shuttelconduction.<br />
Ekspansionsvolumen<br />
Spalte mellem<br />
Stempel og top<br />
Tætningsring<br />
Traditionelt er de to<br />
tabsmekanismer, som<br />
spaltetabet er sammensat af,<br />
modelleret på basis af<br />
simplificerede analytiske<br />
betragtninger. Det er meget<br />
lidt <strong>for</strong>skning inden<strong>for</strong><br />
spaltetab i stirlingmotorer, og<br />
det var der<strong>for</strong> oplagt at<br />
anvende det nye avancerede<br />
simuleringsprogram til at<br />
beregne en mere nøjagtig<br />
værdi <strong>for</strong> spaltetabet i en<br />
given motor.<br />
Simuleringsmodellen blev<br />
der<strong>for</strong> udviddet med en<br />
diskretisering, der <strong>for</strong>tsatte fra<br />
ekspansionsvolumenet ned i<br />
spalten til <strong>for</strong>trængerens<br />
stempeltætning. Figur 4.6<br />
viser resultatet af beregninger<br />
af spaltetabet som funktion af<br />
spalteafstanden <strong>for</strong> 9 kW<br />
motoren. Figuren viser dels<br />
den oftest anvendte algoritme<br />
<strong>for</strong> henholdsvis shuttle<br />
conduction tabet (kurven<br />
markeret med trekanter) og<br />
spaltepumpetabet (kurven markeret med firkanter) præsenteret af Urieli og Berchwitch (U&B,<br />
ref. 2). Desuden vises summen af de to beregnede værdier (kurven markeret med cirkler),<br />
der umiddelbart kan sammenlignes med resultater fra den nye simuleringsmodel. Som det<br />
ses af figuren underestimerer U&B det samlede spaltetab indtil spalten når en værdi på ca.
PSO-FU2204 46 9 kW Stirlingmotor<br />
1,2 mm, hvorefter det analytiske udtryk <strong>for</strong> spaltepumpetabet stiger voldsomt, mens den<br />
detaljerede analyse viser en faldende tendens.<br />
Resultatet er overraskende, <strong>for</strong>di det viser sig, at den valgte spalteafstand på 0,8 – 1 mm er<br />
netop der, hvor det samlede spaltetab har sin maksimale værdi, hvilket er uheldigt. I den<br />
eksisterende NSP til hurtige beregninger af stirlingmotorens ydelse og virkningsgrad er der<br />
der<strong>for</strong> indført en korrektionsfaktor, således at det af Urieli og Berchowitch angivne udtryk<br />
multipliceres med 1,5 samtidig med, at det beregnede spaltetab <strong>for</strong> voksende spalteafstand<br />
begrænses til en værdi, der svarer til det dobbelte af det minimale spaltetab.<br />
Rate of energy transport [kW]<br />
5<br />
4<br />
3<br />
2<br />
1<br />
0<br />
0.10<br />
0.09<br />
0.08<br />
0.07<br />
0.06<br />
0.05<br />
H&B Shuttle+Pump<br />
Quadratic Shuttle+Pump<br />
U&B Shuttle+Pump<br />
U&B Shuttle<br />
U&B Pump<br />
0.20<br />
0.30<br />
0.40<br />
1.00<br />
0.90<br />
0.80<br />
0.70<br />
0.60<br />
0.50<br />
Clearance gap size [mm]<br />
Figur 4.6. Resultater fra beregninger af spaltetabene ved hjælp af det detaljerede simuleringsprogram<br />
sammenlignet med resultater fra de analytiske udtryk, som oftest anvendes ved simulering af<br />
stirlingmaskiner. U&B er en <strong>for</strong>kortelse af Urieli og Berchowitz (ref. 2).<br />
2.00<br />
3.00<br />
4.00<br />
5.00
PSO-FU2204 47 9 kW Stirlingmotor<br />
4.2 Sammenligning af resultater fra <strong>for</strong>søg med resultater fra simulering<br />
Resultater fra beregninger med det oven<strong>for</strong> omtalte Numeriske SimuleringsProgram (NSP),<br />
der er blevet modificeret med hensyn til beregning af regeneratortab og beregning af tryktab,<br />
er sammenlignet med resultater fra laboratorie<strong>for</strong>søgene med stirlingmotoren SM5A-002 (se<br />
afsnit 3.1).<br />
Virkningsgraden <strong>for</strong> en motor er en velegnet størrelse til at vurdere, hvor effektiv motoren er.<br />
Virkningsgraden er imidlertid ikke hensigtsmæssig som parameter til vurdering af<br />
nøjagtigheden af en simulering med et NSP, <strong>for</strong>di det er en beregnet størrelse, der afhænger<br />
af flere <strong>for</strong>skellige målelige størrelser. Til verifikation af NSP’s evne til at beregne<br />
stirlingmotorens egenskaber er det mere hensigtsmæssigt at anvende motorens akseleffekt<br />
og den varmeeffekt, der afsættes i motorens køler. Begge disse størrelser er nemme at<br />
måle, og tilsammen er de også bestemmende <strong>for</strong> motorens virkningsgrad.<br />
Figur 4.7 viser en sammenligning af målte og beregnede resultater <strong>for</strong> akseleffekt og<br />
varmeeffekt afsat i køleren som funktion af hederrørstemperaturen ved henholdsvis 50 bar<br />
og 84 bar middeltryk. Som det fremgår af figuren er der god overensstemmelse mellem<br />
målte og beregnede værdier, og man kan konkludere, at det anvendte stirling-NSP med de<br />
benyttede data <strong>for</strong> 9 kW motoren kan beregne aksel- og køleeffekt inden<strong>for</strong> ca. 10%. I<br />
betragtning af, at simuleringsmodellen alene er baseret på en beskrivelse af termodynamik,<br />
fluiddynamik og varmetransmission uden nogen <strong>for</strong>m <strong>for</strong> tilpasning til den aktuelle opgave<br />
eller beregning, må det anses <strong>for</strong> at være et særdeles tilfredsstillende resultat.<br />
Shaft power, kW<br />
14<br />
13<br />
12<br />
11<br />
10<br />
9<br />
8<br />
7<br />
6<br />
5<br />
SM5A-002, Tcold = 55 C<br />
4<br />
500 550 600 650 700 750<br />
Finned tube temp., °C<br />
8<br />
R 0-P50 bar-calc<br />
500 550 600 650 700 750<br />
Figur 4.7: Sammenligning af målte (meas) og beregnede (calc) resultater <strong>for</strong> akseleffekt<br />
(Shaft power) og varmeeffekt i køleren (Heat from cooler) som funktion af<br />
hederrørstemperaturen ved henholdsvis 50 bar og 84 bar middeltryk.<br />
For yderlige at afprøve modellens <strong>for</strong>udsætninger er der lavet en sammenligning mellem<br />
målte og beregnede data <strong>for</strong> <strong>for</strong>søgsmotoren afprøvet med nitrogen som arbejdsgas. Figur<br />
4.8 viser en sammenligning af målte og beregnede resultater <strong>for</strong> akseleffekt og varmeeffekt<br />
afsat i køleren som funktion af hederrørstemperaturen ved henholdsvis 65 bar og 82 bar<br />
middeltryk med Nitrogen som arbejdsgas. Den eneste <strong>for</strong>skel mellem disse beregninger og<br />
de <strong>for</strong>egående beregninger er, at Helium arbejdsgassen er erstattet med Nitrogen og<br />
Heat from cooler, kW<br />
22<br />
20<br />
18<br />
16<br />
14<br />
12<br />
10<br />
R 0-P84 bar-m eas<br />
R 0-P84 bar-calc<br />
R 0-P50 bar-m eas<br />
Finned tube tem p., °C
PSO-FU2204 48 9 kW Stirlingmotor<br />
driftstilstandene tilpasset de aktuelle <strong>for</strong> <strong>for</strong>søgene. I betragtning af, hvor store krav der stilles<br />
til modellen ved disse beregninger, må resultatet anses <strong>for</strong> at være meget tilfredsstillende.<br />
Akseleffekten beregnes med meget stor nøjagtighed, og varme afsat i køleren beregnes med<br />
tilfredsstillende nøjagtighed.<br />
Af resultaterne ses yderligere, at effekt og intern virkningsgrad falder betragteligt, når Helium<br />
erstattes med Nitrogen. Dette skyldes den højere viskositet og den lavere<br />
varmeledningsevne <strong>for</strong> Nitrogen i <strong>for</strong>hold til Helium. Den interne virkningsgrad er <strong>for</strong>holdet<br />
mellem akseleffekten og den til hederen tilførte effekt, og den kan beregnes ud fra følgende<br />
betragtninger om energibalance:<br />
Qheater=Pshaft+ Qcooler<br />
internal= Pshaft / Qheater = Pshaft / ( Pshaft+Qcooler)<br />
Den interne virkningsgrad ved anvendelse af henholdsvis Helium og Nitrogen kan dermed<br />
beregnes til:<br />
internal, Helium<br />
internal, Nitrogen<br />
Shaft power, kW<br />
8<br />
7<br />
6<br />
5<br />
4<br />
3<br />
= 12kW / ( 12kW + 20kW ) * 100% = 37,5%<br />
= 7,5kW / ( 7,5kW + 22,5kW ) * 100% = 25%<br />
SM5A-002-R0 Nitrogen<br />
2<br />
500 550 600 650 700 750<br />
Finned tube temp., °C<br />
Figure 4.8. Nitrogen som arbejdsgas.<br />
Heat from cooler, kW<br />
26<br />
24<br />
22<br />
20<br />
18<br />
16<br />
14<br />
12<br />
P = 82 bar-meas<br />
P = 82 bar-calc<br />
P = 65 bar-meas<br />
P = 65 bar-calc<br />
10<br />
500 550 600 650 700 750<br />
Finned tube temp., °C<br />
For yderligere at afprøve det modificerede stirling NSP er der lavet en sammenligning af<br />
måle- og beregningsresultater <strong>for</strong> den ene fiberregenerator R2 (SF 22/316/85) med de fine<br />
fibre (fiberdiameter 0,022 mm, porøsitet 84,6%), der er produceret af firmaet Beakert (se<br />
afsnit 3.2). Den anden fiberregenerator har stort set samme data som den oprindeligt<br />
monterede regenerator R0, så den er ikke taget med i sammenligningen.
PSO-FU2204 49 9 kW Stirlingmotor<br />
Figur 4.9 og 4.10 viser en sammenligning af målte og beregnede resultater <strong>for</strong> henholdsvis<br />
akseleffekt og varmeeffekt afsat i køleren som funktion af hederrørstemperaturen ved et<br />
middeltryk på 84 bar, en kølevandstemperatur på 40 o C og med Helium som arbejdsgas. Som<br />
det fremgår af figurerne, er der god overensstemmelse mellem resultater fra<br />
laboratorie<strong>for</strong>søg og simulering. Som beskrevet I afsnit 3.2 er virkningsgraden næsten den<br />
same <strong>for</strong> de to regeneratorer, mens effekten er lidt lavere med fiberregeneratoren, hvilket<br />
simuleringsresultaterne efterviser. Man kan dermed konkludere, at det ikke er en <strong>for</strong>del at<br />
anvende meget tynde tråddiametre i regeneratoren i denne motor.<br />
Shaft power, kW<br />
14<br />
13<br />
12<br />
11<br />
10<br />
9<br />
8<br />
7<br />
SM5A-002, Tcold = 40 C, Pmean = 84 bar<br />
Reg-2-calc<br />
6<br />
500 550 600 650 700 750<br />
Mean tube temp., °C<br />
Reg-0-meas<br />
Reg-0-calc<br />
Reg-2-meas<br />
Figur 4.9. Sammenligning af målt og beregnet akseleffekt <strong>for</strong> original regenerator til 9<br />
kW motor R0 og kort fiber regenerator R2 (d = 0,022 mm, porositet = 84,6%)<br />
Heat from cooler, kW<br />
21<br />
20<br />
19<br />
18<br />
17<br />
16<br />
SM5A-002, Tcold = 40 C, Pmean = 84 bar<br />
Reg-2-calc<br />
15<br />
500 550 600 650 700 750<br />
Mean tube temp., °C<br />
Reg-0-meas<br />
Reg-0-calc<br />
Reg-2-meas<br />
Figur 4.10. Sammenligning af målt og beregnet varme fra køler <strong>for</strong> original<br />
regenerator til 9 kW motor R0 og kort fiber regenerator R2.
PSO-FU2204 50 9 kW Stirlingmotor<br />
4.3 Særlige krav til stirlingmotorer til biomasse<br />
En stirlingmotor består i praksis af 2 cylindervolumener, tre varmevekslere hvoraf den ene er<br />
en regenerator samt nogle dødvolumener, der <strong>for</strong>binder de fem aktive volumener. Det er<br />
vanskeligt at designe disse komponenter, <strong>for</strong>di det er et sammenhængende system, hvor<br />
man ikke kan ændre en parameter, uden at det har indflydelse på de andre. Varmevekslerne<br />
er yderligere komplicerede, <strong>for</strong>di det ikke kun er et spørgsmål om tilstrækkelig overflade og<br />
tværsnitsareal <strong>for</strong> at overføre den nødvendige varme, <strong>for</strong>di det indre volumen i<br />
varmevekslerne også har betydning <strong>for</strong> motorens tryk<strong>for</strong>hold og dermed indflydelse på effekt<br />
og virkningsgrad.<br />
Hvis der der<strong>for</strong> er særlige krav til varmetransmissionen i f.eks. hederen, vil dette ikke alene<br />
have indflydelse på hederens design, men på hele motorens design og grundlæggende data.<br />
Det særlige ved de danske stirlingmotorer er, at de er udviklet specifik til anvendelse af faste<br />
brændsler som flis og andre <strong>for</strong>mer <strong>for</strong> biomasse. Konsekvensen af dette er, at disse motorer<br />
er større og har har et lavere middeltryk end stirlingmotorer beregnet til naturgas eller andre<br />
højværdige brændsler. Det kan sammenlignes med, at der er stor <strong>for</strong>skel på de<br />
varmeoverførende flader i en kedel til træpiller og en væghængt gaskedel, og det er ikke<br />
umiddelbart muligt at ombygge den ene, så den kan erstatte den anden på lige vilkår.<br />
Ved anvendelse af biobrændsler i stirlingmotorer er det ikke alene stirlingmotorens<br />
ud<strong>for</strong>mning der kræver særlig tilpasning, men også <strong>for</strong>brændingssystemet, <strong>for</strong>di det er<br />
nødvendigt at anvende <strong>for</strong>varmet luft til <strong>for</strong>brændingen, hvis man vil opnå en tilfredsstillende<br />
virkningsgrad. Disse <strong>for</strong>skelligartede krav til stirlingmotorens ud<strong>for</strong>mning inklusive den<br />
<strong>for</strong>varmede lufts betydning <strong>for</strong> virkningsgraden har været genstand <strong>for</strong> en særlig<br />
undersøgelse, som er præsenteret ved konferrence i Durham. Efter aftale med værten <strong>for</strong><br />
konferencen er artiklen vedlagt rapporten som bilag (se bilag 6).
PSO-FU2204 51 9 kW Stirlingmotor<br />
5 Undersøgelse af særlige komponenter og driftsbetingelser<br />
I projektbeskrivelsen var det <strong>for</strong>udsat, at der blev lavet en undersøgelse af specifikke emner<br />
af betydning <strong>for</strong> stirlingmotorens konstruktion og drift. Undersøgelsen af mulighederne <strong>for</strong> en<br />
reduktion af frigangen i lejerne ved anvendelse af Finite Element Analysis (FEA) er<br />
behandlet i afsnittet om konstruktionen af den nye motor. Neden<strong>for</strong> er en række andre emner<br />
behandlet<br />
5.1 Dellast<br />
De danske stirlingmotorer er beregnet til små kraftvarmeværker, der anvender biobrændsler.<br />
Regulering af motorernes effekt har været indskrænket til on/off regulering på samme måde,<br />
som gasmotorerne i små kraftvarmeværker reguleres. I det fremtidige energimarked er<br />
kraftvarmeanlægenes evne til at modulere vigtigt, <strong>for</strong>di det er nødvendigt med en hurtig<br />
tilpasning til variationerne i el-produktionen fra vindmøllerne og andre <strong>for</strong>mer <strong>for</strong> ureguleret<br />
produktion. Hvis motorens skal tilsluttes et lokalt net, hvor motoren også skal styre frekvens<br />
og spænding på nettet, er effektregulering en nødvendighed.<br />
Det er ikke så enkelt at styre en stirlingmotors effekt som det er at styre effekten <strong>for</strong> en<br />
konventionel motor. Der er flere tekniske muligheder <strong>for</strong> at reducere effekten til en<br />
dellasttilstand, og de <strong>for</strong>skellige muligheder er der<strong>for</strong> undersøgt. Der er følgende muligheder<br />
<strong>for</strong> en regulering af en Stirlingmotors effekt:<br />
• Reduktion af middeltryk<br />
• Større dødvolumen<br />
• Ventil til kontrol af tryk<strong>for</strong>løbet i cylinderen<br />
• Reduktion af rørtemperatur i hederen<br />
• Ændring af fasevinkel mellem ekspansions- og kompressionsvolumen<br />
• Reduktion af volumenvariationens størrelse <strong>for</strong> arbejdsstempel eller <strong>for</strong>trænger<br />
Neden<strong>for</strong> er de <strong>for</strong>skellige metoder behandlet i detaljer. Undersøgelsen er baseret på<br />
resultater fra den simuleringsmodel, som er verificeret andre steder i nærværende rapport.<br />
Simuleringsmodellen behandler alene <strong>for</strong>hold i motoren, og der er i undersøgelsen ikke taget<br />
højde <strong>for</strong> de tab, som stammer fra fyringssystemet og overføringen af varmen til<br />
stirlingmotoren. Disse tab er individuelle <strong>for</strong> de enkelte fyringssystemer og type af brændsler,<br />
og resultaterne bør der<strong>for</strong> korrigeres, hvis de skal anvendes til vurdering af et konkret<br />
stirlingmotoranlæg.<br />
I kapitel 3 er resultaterne fra målingerne af 9 kW motorens effekt og virkningsgrad beskrevet<br />
inklusive målinger, hvor hederørenes temperatur og middeltrykket i cylinderen er varieret.<br />
Disse data har naturligvis relevans <strong>for</strong> undersøgelsen af motorens egenskaber under dellast,<br />
men det er valgt alene at benytte disse data til at verificere og <strong>for</strong>bedre simuleringsmodellen,<br />
og derefter anvende simuleringsmodellen til undersøgelsen af motorens egenskaber under<br />
dellast. Se i øvrigt kapitel 3.2 og bilag 5 <strong>for</strong> målinger af <strong>for</strong>skellige driftstilstande.
PSO-FU2204 52 9 kW Stirlingmotor<br />
Virkningsgrad<br />
0,4<br />
0,35<br />
0,3<br />
0,25<br />
0,2<br />
0,15<br />
0,1<br />
0,05<br />
0<br />
Dellast<br />
0 0,2 0,4 0,6 0,8 1<br />
Specifik effekt<br />
Figur 5.1: Virkningsgradskurver <strong>for</strong> <strong>for</strong>skellige dellastløsninger.<br />
Temperatur<br />
Slaglængde<br />
Fase vinkel<br />
Middeltryk<br />
Figur 5.1 viser en sammenligning af de <strong>for</strong>skellige metoder til dellastregulering. Data til<br />
simuleringen er hentet fra 9 kW motoren, selvom ikke alle de undersøgte muligheder<br />
umiddelbart ville kunne anvendes på denne motor.<br />
5.2 Ventiler<br />
Der er to typer ventiler i motoren. I den én-cylindrede motor er der to én-vejsventiler, der<br />
sørger <strong>for</strong> at <strong>for</strong>syne to små beholdere med henholdsvis maksimumstryk og minimumstryk til<br />
pneumatisk styring af den anden type ventil, der består af en tvangsstyret ventil, der benyttes<br />
ved opstart og ved strømsvigt, så motoren ikke løber løbsk. Disse ventiler er identiske med<br />
de ventiler, der anvendes i den 4-cylindrede motor, idet der i denne motor også sidder en énvejsventil<br />
i hver af de fire cylindre, således at alle fire cylindervolumener arbejder med<br />
samme minimumstryk.
PSO-FU2204 53 9 kW Stirlingmotor<br />
5<br />
1<br />
2<br />
4<br />
3<br />
Figur 5.1: Eksploderet tegning af envejsventil.<br />
Figur<br />
Figur 5.1:<br />
5.2: 5.1:<br />
Eksploderet<br />
Eksploderet tegning af envejsventil.<br />
tegning<br />
af<br />
af<br />
envejsventil.<br />
envejsventil.<br />
Særligt én-vejsventilerne har givet mange problemer. Figur 5.2 viser en tegning af énvejsventilerne,<br />
der består af 1) ventilhus, 2) ventillegeme, 3) ventilsæde, 4) PTFE-tætning, og<br />
5) returfjeder. Når trykket er størst på den side, hvor ventillegemet sidder, trykkes den<br />
kegle<strong>for</strong>mede del af ventillegemet mod tætningen indtil det rammer ventilsædet.<br />
Problemerne med tætningerne var dels slid på PTFE-delene og utilstrækkelig præcision ved<br />
fremstillingen af tætningerne, og dels slid på ventillegeme og -sæde.<br />
Der gik <strong>for</strong>holdsvis lang tid fra ansøgningen af nærværende projekt til at <strong>projektet</strong> kunne<br />
startes, og der var allerede inden <strong>projektet</strong>s start fundet en ny løsning til PTFE-tætningen i<br />
én-vejsventilerne. En ny leverandør, der leverede PTFE-tætninger med høj præcision og<br />
finish, havde en standard tætning på programmet med tilnærmelsesvis de samme mål som<br />
den, der hidtil blev anvendt. Ventilkeglen i ventilens bevægelige del blev ændret, så den<br />
passede til den nye geometri, og samtidig blev der valgt et vacuum-hærdningdstål til<br />
ventilkegle og –sæde i stedet <strong>for</strong> det hidtil anvendte stål. Siden at de gamle type énvejsventiler<br />
er blevet erstattet af de nye, har der ikke været problemer med én-vejsventilerne.
PSO-FU2204 54 9 kW Stirlingmotor<br />
Figur 5.3 viser bypassventilen, der i den én-cylindrede motor <strong>for</strong>binder det aktive<br />
cylindervolumen med krumtaphusets volumen, og i de fire-cylindrede motorer <strong>for</strong>binder de<br />
aktive cylindervolumener. Når ventilen er åben, medfører <strong>for</strong>bindelsen mellem<br />
cylindervolumenet og et volumen med konstant tryk, at motoren ikke producerer nogen<br />
effekt. Når ventilen er lukket, er der ingen <strong>for</strong>bindelse mellem cylindervolumenet og<br />
volumenet med konstant tryk, vil motoren producere effekt under <strong>for</strong>udsætning af, at<br />
temperaturen på hederen er tilstrækkelig høj. Det betyder, at det er muligt at stoppe motoren,<br />
selvom hedertemperaturen er høj. Dermed undgår man, at motoren løber løbsk, hvis nettet<br />
falder ud f.eks. som følge af et lynnedslag.<br />
1<br />
2<br />
3<br />
4<br />
Figur 5.3: Eksploderet tegning af bypass ventilen.<br />
5<br />
a<br />
Som vist på figur 5.3 består<br />
servoventilen af 1) ventilhus, 2)<br />
stempel med stempelring, 3)<br />
returfjeder, 4) sæde, og 5)<br />
PTFE-tætning. Når der tilføres<br />
Helium styretryk fra et<br />
højtrykssystem ved a) presser<br />
stemplet ventillegemet mod<br />
sædet, og ventilen er lukket.<br />
Fjernes trykket, returnerer<br />
fjederen stemplet, og ventilen<br />
er åben.<br />
Ventilsædets vinkel har<br />
betydning <strong>for</strong> slid på<br />
ventiltætningen og ventilens<br />
tæthed. Der er der<strong>for</strong> afprøvet<br />
<strong>for</strong>skellige vinkler mellem 0 o og<br />
45 o . Det blev konkluderet, at det<br />
bedste resultat blev opnået<br />
med et ventillegeme med en<br />
vinkel på 45 o som vist på figur<br />
5.3.
PSO-FU2204 55 9 kW Stirlingmotor<br />
5.3 Lejetætninger og smørefedt<br />
I de hermetiske motorer smøres lejerne med fedt, der holdes inde i lejerne ved hjælp af<br />
læbetætninger. Også her er der et betydeligt friktionstab og motorens virkningsgrad og<br />
akseleffekt kan <strong>for</strong>bedres ved at nedbringe friktionen i disse tætninger. Sammenhængen<br />
mellem tætningsdesign og eftersmøringsinterval vil også blive undersøgt. Smøringen i de<br />
oscillerende lejer er kritisk og atypisk. Der er behov <strong>for</strong> et smøremiddel med gode<br />
grænsesmøringsegenskaber i en Heliumatmosfære ved 80°C. Forsøg med specielle<br />
smøremidler kan vise sig at være væsentlige.<br />
I de hermetiske motorer benyttes rulle- og nålelejer, der smøres med fedt. Fedtet holdes inde<br />
i lejerne ved hjælp af læbetætninger. Denne konstruktion benyttes <strong>for</strong> at undgå problemer<br />
med, at smøreolie fra hydrodynamisk smurte lejer trænger ind i arbejdsvolumenerne med<br />
risiko <strong>for</strong>, at olie omdannes til tjæreagtige <strong>for</strong>bindelser i den varme del af cylinderen og<br />
tilstopper varmeveklser og regenerator.<br />
Figur 5.4: Eksploderet tegning af DOB tætning<br />
med indbygningsdele.<br />
I den nuværende konstruktion<br />
anvendes to typer tætninger ved<br />
lejerne. Den ene tætning er standard<br />
tætninger beregnet til montage<br />
sammen med nålelejer. Disse<br />
tætninger er nemme at montere og<br />
de fungerer særdeles godt, men de<br />
kan kun fås i nogle ganske få<br />
diametre over ø50 mm. Den anden<br />
type er nogle specielle tætninger af<br />
typen DOB (se figur 5.4), som<br />
anvendes til de større diametre,<br />
leveres i Danmark af Ingeniørfirmaet<br />
GunnarHaagensen A/S. Disse<br />
tætninger består af en membran af<br />
gummi, der er spændt ud over en<br />
stålskive med et hul, hvor selve<br />
tætningslæben sidder. Fordelene<br />
ved denne type tætning<br />
er, at indbygningsbredden er lille, at<br />
det er muligt at få dem i mange<br />
<strong>for</strong>skellige størrelser, og at friktionen<br />
er lille. Tætningerne er <strong>for</strong> så vidt<br />
også udmærkede, men der er nogle problemer i <strong>for</strong>bindelse med indbygningen. Tætningen<br />
kræver en støtteplade på siderne, og konstruktionen kræver der<strong>for</strong> fremstilling af ekstra<br />
komponenter. Desuden er der en risiko <strong>for</strong>, at læbetætningen kommer til at vende ”den<br />
<strong>for</strong>kerte vej” under montage.<br />
For at imødegå de oven<strong>for</strong> nævnte mindre hensigtsmæssige egenskaber <strong>for</strong> de eksisterende<br />
specialtætninger fra Haagensen, er mulighederne <strong>for</strong> at finde andre løsninger til tætning af<br />
lejerne undersøgt. Trods en nøje gennemgang af markedet var det ikke muligt at finde andre<br />
løsninger end de traditionelle simmerringe med fjederaktivering samt tætninger fremstillet af<br />
PTFE. Simmerringene kan ikke bruges, <strong>for</strong>di de er 8 - 12 mm brede, hvilket der ikke er plads<br />
til, og <strong>for</strong>di fjederaktiveringen giver unødigt høje friktionstab.<br />
PTFE er ikke et elastisk materiale som gummi. Tætningerne fremstillet af fyldt PTFE er<br />
der<strong>for</strong> også fjederaktiverede, og funktionen er baseret på plastiske de<strong>for</strong>mationer som<br />
kompensation <strong>for</strong> slid. Tætningerne fremstilles på en særlig drejebænk af PTFE-emnerør, og
PSO-FU2204 56 9 kW Stirlingmotor<br />
det er der<strong>for</strong> muligt at få præcis den geometri og diameter, der ønskes. Tætningerne er<br />
imidlertid dyre (3 – 4 gange prisen <strong>for</strong> gummitætningerne), og der bliver meget nemt ødelagt<br />
under montage. Et <strong>for</strong>søg med anvendelse af PTFE-tætninger i 9 kW motoren nåede der<strong>for</strong><br />
kun til montagen, hvor den blev så skadet, at <strong>for</strong>søget blev opgivet.<br />
Det er muligt at bestille tætninger fremstillet efter særlige mål, men det kræver, at man<br />
bestiller mange tusinde af gangen. Konklusionen er der<strong>for</strong>, at det <strong>for</strong>eløbigt er nødvendigt at<br />
klare sig med de tætninger fra Haagensen, som anvendes i dag.<br />
Smørefedt.<br />
Valg af smørefedt til lejerne er også problematisk. Dels skal olien i fedtet have et lavt<br />
damptryk <strong>for</strong> at undgå, at der kommer oliedampe op i cylindervolumenerne, og dels skal<br />
viskositeten ved driftstemperaturen være tilstrækkelig til at bevare den ønskede smørefilm<br />
under drift. Den tryksatte Heliumatmosfære i krumtaphuset får i øvrigt flere typer fedt til at<br />
dekomponere, så olien bliver presset ud af sæbematrixen.<br />
Hidtil er der anvendt to <strong>for</strong>skellige typer fedt i lejerne, idet der anvendes én type fedt i de<br />
lejer, der roterer en hel omgang under drift, og en anden type fedt til de lejer i mekanismen,<br />
der kun rokker frem og tilbage ved lav vinkelhastighed. Årsagen er, at hastighederne i lejet er<br />
afgørende <strong>for</strong>, om fedten bliver ”fanget” i kontaktzonen eller den slipper væk. Den<br />
nødvendige baseolieviskositet kan bestemmes ud fra et krav om, at filmtykkelsen i<br />
kontaktzonen er tilstrækkelig til at sikre ruhedstopadskillelse. Vælges en fedttype med<br />
unødvendig høj baseolieviskositet opnås ingen levetids<strong>for</strong>længelse, kun øget friktion og<br />
egenopvarmning. Følgende fedt er anvendt <strong>for</strong> at optimere begge lejetyper:<br />
LGEP2 til rullelejerne på krumtapakslen, der roterer en fuld omdrejning pr. omdrejning.<br />
Baseolieviskositet 200 mm 2 /s.<br />
LGEM2 til ågenes og svingarmenes nålelejer, der kun roterer +/- 12 o – 17 o pr. omdrejning<br />
Baseolieviskositet 500 mm 2 /s<br />
Generelt er disse fedttyper dekomponeret i relativ høj grad, når en motor har været i drift. Det<br />
er normalt, at der er lille sø af olie i bunden af krumtaphuset ved demontering, selvom alt var<br />
tørt ved samling. Problemet er lidt mindre i de fire-cylindrede motorer, <strong>for</strong>modentlig <strong>for</strong>di de<br />
kører med lavere tryk i krumtaphuset, 35 bar imod 80 bar på de en-cylindrede motorer.<br />
I <strong>for</strong>bindelse med afprøvningen af de 4-cylindrede stirlingmotorer i Østrig viste det sig, at de<br />
høje temperaturer i returledningen fra fjernvarmesystemerne førte til høje temperaturer i<br />
krumtaphuset og dermed også i lejerne. Det kunne det anvendte fedt ikke tåle, hvilket førte<br />
til, at lejer og krumtap blev ødelagt. Der blev der<strong>for</strong> sat en undersøgelse i gang i samarbejde<br />
med den østrigske partner og firmaet Klüber, der er specialister i smørefedt. De har<br />
efterfølgende <strong>for</strong>eslået følgende fedttyper:<br />
Klüberalfa BHR 53-402, Højtemperatur langtidssmørefedt. Pris ca. 3400 kr/kg<br />
Klübersynth BH 72-422, Langtids og højtemperatur fedt. Pris ca. 1200 kr/kg<br />
Det særlige stof der findes i <strong>for</strong>skellige koncentrationer i de to fedttyper er per-floureredepolyether-olie,<br />
også kendt som PFPE olie. PFPE olie bliver brugt mange steder i industrien<br />
som problemknuser ved tribologiske problemer.<br />
Klüber har <strong>for</strong>etaget beregninger, der viser, at vi kan <strong>for</strong>vente fedtlevetid omkring 20000<br />
timer ved en krumtaphustemperatur på 85°C. Til sammenligning havde vi tidligere omkring<br />
4000 timer. Det valgte fedt er nu under afprøvning, og trods prisen ser det ud til at være en<br />
attraktiv løsning.
PSO-FU2204 57 9 kW Stirlingmotor<br />
5.4 Brint som arbejdsgas<br />
Stirling motorer arbejder ofte med Helium som arbejdsmedie. En del stirlingmotorer med høj<br />
specifik ydelse anvender dog brint som arbejdsgas, mens motorer med lav specifik ydelse<br />
ofte anvender luft som arbejdsmedie. For at opstille et grundlag til at vurdere brint i <strong>for</strong>hold til<br />
Helium, må virkningen af de termodynamiske <strong>for</strong>dele ved brint demonstreres.<br />
Fordelen ved anvendelse af brint (H2)i stedet <strong>for</strong> Helium (He) skyldes først og fremmest to<br />
<strong>for</strong>hold. Det første er en følge af, at brints molekylvægt kun er 2, mens Heliums molekylvægt<br />
er 4. Samtidig er brint en to-atomig gas, mens Helium er en mono-atomig gas. Ud fra den<br />
statistiske termodynamik kan den specifikke varmekapacitet cv <strong>for</strong> en idealgas i det relevante<br />
temperaturområde beregnes ud fra gaskonstanten og molekylernes frihedsgrader:<br />
^ ^<br />
Specifik varmekapacitet, Helium (mono-atomig gas) cp, He = 5/2 R/M<br />
Specifik varmekapacitet, Brint (to-atomig gas) cp,H2 = 7/2 R/M<br />
^ ^<br />
Hvor R : universalgaskonstanten<br />
M : Molekylvægten <strong>for</strong> den pågældende gas<br />
Molekylvægten <strong>for</strong> Helium er dobbelt så stor som molekylvægten <strong>for</strong> brint, hvilket medfører<br />
at <strong>for</strong>holdet mellem de to gassers specifikke varmekapacitet er cv, He / cv, H2 = 5/14 = 0,357. En<br />
beregning ud fra gasdata viser, at det faktiske <strong>for</strong>hold ved 20 o C og 1 bar er cv, He / cv, H2 =<br />
0,364, hvilket er en god overensstemmelse, når man tager simplifikationerne i betragtning.<br />
Ud fra definitionen af Prandl’s tal fås, at Pr = cp/k, hvor er viskositeten og k er<br />
varmeledningsevnen. Prandl’s tal kan <strong>for</strong> en idealgas tilnærmes med Pr = 0,7 uafhængig af,<br />
hvilken gas, der er tale om, hvilket medfører, at <strong>for</strong>holdet mellem viskositeten og<br />
varmeledningsevnen /k stiger med stigende molekylevægt og falder med stigende antal<br />
frihedsgrader. Molekylvægten er den dominerende faktor, så stigende molekylvægt medfører<br />
stigende <strong>for</strong>hold /k. Når viskositeten stiger, stiger tryktabet i en given kanalstrømning eller<br />
varmeveksler, og når varmeledningsevnen falder, stiger temperatur<strong>for</strong>skellen mellem væg og<br />
gas. Man kan dermed konkludere, at en let gas fører til højere effekt og virkningsgrad i en<br />
stirlingmotor.<br />
Et andet <strong>for</strong>hold af betydning er hysteresetabet, som opstår som følge af, at arbejdsgassens<br />
temperatur stiger under kompression og falder under ekspansion. Når temperaturen i gassen<br />
er højere en vægtemperaturen, overføres varme fra gas til væg, mens det omvendte er<br />
tilfældet, når gastemperaturen er lavere. Den cykliske varmeovergang mellem gas og væg er<br />
årsag til et tab, som benævnes hysteresetabet. Dette tab er relativt stort i den fire-cylindrede<br />
motor på grund af de store dødvolumener i den komplekse heder.<br />
Som oven<strong>for</strong> nævnt er Helium er en mono-atomig ædelgas i modsætning til brint, der er en<br />
di-atomig gas. Som det fremgår af ovenstående udredning, er isentropkoefficienten k = cp/cv<br />
= cp/(cp – R) højere <strong>for</strong> monoatomige gasser (ca. 1,7), end <strong>for</strong> di-atomige gasser (ca. 1,4).<br />
Ved at benytte brint i stedet <strong>for</strong> Helium kan motorens hysteresetab reduceres, <strong>for</strong>di<br />
temperatur og trykændringerne <strong>for</strong> en given volumenændring er mindre. Det betyder, at<br />
motorens virkningsgrad bliver højere og at belastningen på motorens mekaniske dele bliver<br />
mindre.<br />
Ud fra en teoretisk betragtning kan man dermed konkludere, at brint er en bedre arbejdsgas<br />
til stirlingmotorer end Helium. Fordelen ved at anvende brint afhænger imidlertid af den<br />
aktuelle motors konstruktion. For at belyse dette er der lavet en beregning af effekt og<br />
virkningsgrad <strong>for</strong> den nye 9 kW motor SM5B med både Helium og brint som arbejdsgas ved
PSO-FU2204 58 9 kW Stirlingmotor<br />
samme driftstilstand. Resultaterne er sammenlignet med tilsvarende beregninger <strong>for</strong> 35 kW<br />
motoren SM5E (se tabel 5.4)<br />
Tabel 5.4<br />
Effekt pr. cylinder, kW Virkningsgrad, %<br />
Helium Brint Helium Brint<br />
SM5B (9 kW) 11,74 12,64 36,2 38,9<br />
SM3E (35 kW) 11,17 14,60 32,3 40,4<br />
Tabellen viser den beregnede akseleffekt pr. cylinder og motorens virkningsgrad.<br />
Virkningsgraden er defineret som akseleffekt i <strong>for</strong>hold til den varme, der tilføres hederen og<br />
kan dermed betragtes som en virkningsgrad <strong>for</strong> selve motoren uden fyringssystem. Af<br />
tabellen ses, at <strong>for</strong>delen ved at anvende brint i stedet <strong>for</strong> Helium i den én-cylindrede 9 kW<br />
motor (SM5B) er væsentlig mindre end <strong>for</strong> 35 kW motoren (SM3E). Dette skyldes dels den<br />
lille volumenfasevinkel i 35 kW motoren, som er en følge af den dobbeltvirkende<br />
firecylindrede konstruktion, og dels en følge af, at strømningstabene i hederen til biomasse<br />
er væsentlig større.<br />
Der er imidlertid andre <strong>for</strong>hold af betydning <strong>for</strong> valg af brint i stedet <strong>for</strong> Helium i en<br />
stirlingmotor. Prisen taler også til <strong>for</strong>del <strong>for</strong> brint, men der er også alvorlige praktiske<br />
problemer <strong>for</strong>bundet med brint. Dels kan brint ødelægge de materialer, der indgår i<br />
stirlingmotorer især i de varme dele, og dels er de sikkerhedsmæssige aspekter ved<br />
håndteringen af brint mere komplicerede. Disse <strong>for</strong>hold er undersøgt tidligere, og resultatet<br />
er vedlagt i bilag 2.<br />
Af sikerhedsmæssige årsager og <strong>for</strong>di 9 kW motoren ikke vil opføre sig væsentlig anderledes<br />
med brint i stedet <strong>for</strong> Helium som arbejdsgas, blev det besluttet ikke at gennemføre<br />
laboratorie<strong>for</strong>søg med brint. I stedet er der gennemført <strong>for</strong>søg med Nitrogen som arbejdsgas,<br />
hvilket ud fra et udviklingsmæssigt synspunkt var væsentlig mere interessant. Disse <strong>for</strong>søg<br />
er beskrevet nærmere under afsnittet vedrørende simuleringsmodellen.<br />
5.5 Simulering af arbejdsåg med indbyggede rulle og nålelejer.<br />
Kræfterne virker på akslerne i midten af lejerne. For at simulere kræfternes påvirkning af<br />
åget er det nødvendigt at simulere lejernes stivhed. Der<strong>for</strong> er finite element modellen<br />
<strong>for</strong>synet med segmenter af elementer, der <strong>for</strong>binder lejehullernes periferi med akslen.<br />
Kræfterne kan derefter lægges på akslerne.<br />
Lejesegmenternes mål er givet ved lejets radier, rullernes længde, og en deling der giver et<br />
korrekt antal ruller. Lejesegmenterne i de <strong>for</strong>skellige lejer defineres med <strong>for</strong>skellige<br />
materialeparametre, således at materialenes elastisitetsmodul tilpasses således, at et<br />
segment stivhedsmæssigt svarer til en rulle i lejet.<br />
Kontaktarealet mellem en rulle og en bane har <strong>for</strong>m som et meget smalt rektangel. Bredden<br />
afhænger af belastningen og rullens diameter. Rullens flytning skyldes primært<br />
kompressionen af materialet lokalt i kontaktzonen og ikke en generel de<strong>for</strong>mation af rullen.<br />
Inderbane og yderbane belastes samtidig, og den samlede flytning kan dermed udtrykkes<br />
som 2 x u, hvor u er flytningen som følge af en kontaktzonede<strong>for</strong>mation. Flytningen u er<br />
udtrykt ved en empirisk <strong>for</strong>mel, der giver flytningen af centrene som funktion af rullernes<br />
længde og størrelsen af belastningen på en rulle:<br />
(Reference: 3: s.260 ligning 13/8)
PSO-FU2204 59 9 kW Stirlingmotor<br />
u ≅ 3,<br />
98*<br />
10<br />
u = 3,<br />
98*<br />
10<br />
u = 1,<br />
94μm<br />
−5<br />
−5<br />
2*<br />
u = 3,<br />
9μm<br />
*<br />
*<br />
0,<br />
925 −0,<br />
85<br />
P * l [ mm]<br />
0,<br />
925<br />
−0,<br />
85<br />
( 1000N<br />
) * ( 19mm)<br />
[ mm]<br />
Da rullerne i rullelejerne på krumtapakselen har samme længde som nålene i nålelejerne,<br />
kan en bjælke bestemmes, der passer til begge lejetyper. Rullernes længde er 19 mm.<br />
Når et sæt liniekontakter med en længde på 19 mm belastes med en kraft på 1 kN, flyttes<br />
emnerne 3,9 μm mod hinanden. Dette respons skal modelleres med en bjælke med<br />
rektangulært tværsnit og en længde på 22,5 mm <strong>for</strong> krumtaplejet, og 6,5 mm <strong>for</strong><br />
vippelejerne. Bjælkens tøjning, dvs. den relative <strong>for</strong>længelse, betegnes e og beregnes som:<br />
e,krumtapleje = 3,9 μm / 22,5 mm = 0,000173<br />
e,vippelejer=3,9µm/6,5mm=0,0006<br />
Idet der gælder : F = e*E*A => E = F / (A*e),<br />
hvor F: kraften (sat til 1 kN)<br />
A: tværsnitsareal<br />
A,vippelejer = 19mm x 6,3mm = 119,7mm 2<br />
A,krumtap = 19mm x 12,2mm = 231,8mm 2<br />
E: elastisitetsmodulet (= 205 000 N/mm 2 <strong>for</strong> stål)<br />
E,vippelejer = 1000N / (119,7mm 2 x 0,0006) = 13923,7 MPa = 14 GPa<br />
E,krumtap = 1000N / (231,8mm 2 x 0,000173) = 24936,8 MPa = 25 GPa<br />
Det antages i modellen, at lejerne er slørfri, dvs. at alle ruller er i kontakt med både<br />
inderbane og yderbane ved montage, og alle ruller er ubelastede, når lejet i øvrigt er<br />
ubelastet. Når lejet belastes, vil rullerne i den ene side af lejet bære lasten, og rullerne i den<br />
anden side vil slippe banen. I modellen fjernes bjælkerne i den ubelastede side, da der ellers<br />
vil være trækspændinger i dem i modsætning til rullerne i det virkelige leje.<br />
At lokalisere og fjerne de trækbelastede bjælker er en itterativ proces, <strong>for</strong>di de<strong>for</strong>mationen af<br />
lejehuset har en indflydelse på, hvilke bjælker der belastes, og de<strong>for</strong>mationen afhænger af<br />
antallet af bjælker i modellen.<br />
Figur 5.5 viser finite element modellen af arbejdsåget med indspændinger. Geometrien er<br />
importeret fra 3d tegneprogrammet og tilpasset ANSYS analysen. Indspændingerne ses som<br />
en cylindrisk fixering af krumtapakslen nederst, dvs. at positionen er fast med selve akslen<br />
kan rotere. Det øverste leje midt<strong>for</strong> er i motoren <strong>for</strong>bundet til en svingarm. I modellen er<br />
indersiden af akslen fastholdt i x-aksens retning.
PSO-FU2204 60 9 kW Stirlingmotor<br />
Figur 5.5: Finite element model af arbejdsåg med indspændinger.<br />
På figur 5.6 er belastningen lagt på. Stempelkraften topper ved 20 kN og bliver <strong>for</strong>delt på to<br />
symetriske åg. Maksimallasten på et åg er 10kN. For at simulere et centralt angrebspunkt på<br />
den asymetriske aksel er lasten delt op i 3500 N + 6500 N og lagt på hver ende af akselen.<br />
Figur 5.6: Belasning på åg.
PSO-FU2204 61 9 kW Stirlingmotor<br />
Figur 5.7 viser det automatisk genererede net. Nettet er valgt som hexagonalt domineret.<br />
Modellen er fra 3d tegneprogrammet en samlingstegning af flere komponenter. Disse<br />
komponenter er <strong>for</strong>tsat afgrænset som delelementer af den samlede model. Denne opdeling<br />
gør det muligt at definere et materiale med lavere E-modul <strong>for</strong> lejesegmenterne der simulerer<br />
rullerne i de fysiske lejer. Opdelingen giver også mulighed <strong>for</strong> at tilføje og fjerne<br />
lejesegmenter uden at definere modellen <strong>for</strong>fra, idet alle lejesegmenter er med i modellen,<br />
nok til at fylde lejerne helt. De segmenter der er udsat <strong>for</strong> trækspændinger er undertrykt så<br />
de ikke indgår i beregninger eller grafik. Omkostningen ved denne sammensatte<br />
modelopbygning er komplikationer i beregningerne idet der må defineres kontaktzoner og<br />
bestemmes belastninger hvor delelementerne mødes. Programmet håndterer denne<br />
problematik automatisk, men beregningstiden og stabiliteten påvirkes. Kontaktzonernes<br />
fysiske egenskaber er defineret med en række parametre, der skal kontrolleres og vælges<br />
manuelt.<br />
Figur 5.7: Model af åg med automatisk genereret net.<br />
Figur 5.8 viser resultatet af spændingsberegningerne. Spændinger er vist med farveskala og<br />
de<strong>for</strong>mationer er vist idet modellens de<strong>for</strong>mation er <strong>for</strong>stærket 200 gange således at den<br />
de<strong>for</strong>merede facon bliver tydelig. Simuleringen er gennemført med et konvergenskrav på<br />
10%, dvs. at efter første simulering <strong>for</strong>finer programmet nettet og gentager simuleringen indtil<br />
spændingsniveauet afviger mindre end 10% fra <strong>for</strong>rige simulering.
PSO-FU2204 62 9 kW Stirlingmotor<br />
Maksimalspændingen er 56MPa – hvilket svarer til en sikkerhedsfaktor omkring 4,3 idet<br />
udmattelsesgrænsen <strong>for</strong> ck45 er ca. 240MPa.<br />
De<strong>for</strong>mationen af lejehusene er væsentlig. Dels er faconen vigtig <strong>for</strong>di lejet skal understøttes<br />
af en jævn cirkel<strong>for</strong>m <strong>for</strong> at undgå lokale belastningskoncentrationer. Dette fænomen er<br />
årsag til at ågets profil er valgt som vist med en lodret side til højre, i stedet <strong>for</strong> direkte<br />
tangentiel <strong>for</strong>bindelse. Dels giver de<strong>for</strong>mationen anledning til generering af lejestøj når<br />
lejehuset de<strong>for</strong>meres, så der opstår klaring i lejerne. Klaring i det store leje på krumtappen er<br />
ikke problematisk, idet belastningens retning følger en roterende vektor. Derimod er klaring i<br />
vippelejerne meget mere støjgenererende, <strong>for</strong>di belastningen kun har to modsat rettede<br />
retninger. Det vurderes, at ågets stivhed og styrke er tilfredsstillende. Se i øvrigt bilag 7 <strong>for</strong><br />
uddybende materiale.<br />
Figur 5.8: Resultat af finite element analyse. Viser ækvivalente Von Mise spændinger og de<strong>for</strong>mationer<br />
med en <strong>for</strong>størrelse på 200X.
PSO-FU2204 63 9 kW Stirlingmotor<br />
6 Evaluering af teknologi og marked<br />
Et kraftvarmeanlæg dimensioneres til at dække varmebehovet hos kunden, som f.eks. kan<br />
være en fjernvarme<strong>for</strong>ening, et boligselskab, en højskole eller et landbrug. Samtidig med, at<br />
anlægget producerer varme til opvarmning af boliger, erhvervslokaler eller andet, producerer<br />
anlægget også el, der enten anvendes af kunden selv eller sælges til nettet.<br />
Fordelen ved at anskaffe et kraftvarmeanlæg sammenlignet med blot at anskaffe en<br />
varmekedel er, at strøm er dyrere end varme. Man kan dermed spare udgifter til strøm<br />
alternativt sælge strøm til nettet <strong>for</strong> en højere pris end merudgiften til brændsel.<br />
Konkurrence<strong>for</strong>delene <strong>for</strong> de udviklede stirlingmotorer til kraftvarme er:<br />
• Motorerne gør det muligt at etablere decentrale kraftvarmeværker med faste brændsler<br />
på de små varmemarkeder. Det er <strong>for</strong>eløbigt det eneste produkt, som kan<br />
markedsførers inden<strong>for</strong> en overskuelig fremtid<br />
• Høj energiudnyttelse<br />
• Lave emissioner og ingen miljøskadelige restprodukter<br />
• Lavt behov <strong>for</strong> service og dermed også lave udgifter til drift og vedligehold<br />
• Fuldautomatisk. Anlægget vil normalt køre i ubemandet drift. Det vil konkurrerende<br />
teknologier næppe opnå <strong>for</strong>eløbig<br />
• Mulighed <strong>for</strong> at afpasse driften til døgnvariationer i el-priser. Anlæggene kan køre om<br />
dagen, når el-prisen er høj, og stoppes om natten, når el-prisen er lav. Automatisk start<br />
og stop en gang i døgnet bliver meget vanskeligt at gennemføre <strong>for</strong> konkurrerende<br />
teknologier.<br />
Den firecylindrede stirlingmotor er designet som et modul, så to motorer kan kobles sammen<br />
til en enkelt motor med den dobbelte effekt. Man kan altså umiddelbart sælge både en 4cylindret<br />
35 kW-el og en 8-cylindret 70 kW-el motor, når udviklingsarbejdet er tilendebragt.<br />
Kunderne i Europa og Canada kan <strong>for</strong>deles på følgende grupper:<br />
1. store luksusejendomme, hoteller, sammenslutninger af parcelhuse i små klynger, …<br />
2. små fjernvarmecentraler til boligopvarmning<br />
3. virksomheder <strong>for</strong>trinsvis i træindustrien<br />
4. landbrug og gartnerier<br />
I Sydøst Asien og Sydamerika er kunderne f.eks.<br />
1. el-<strong>for</strong>syning i landområderne, hvor stirlingmotoranlæggene f.eks. anvender risskaller<br />
som brændsel til drift af rismøller, småindustri m.m.<br />
2. træindustrien<br />
Neden<strong>for</strong> er de enkelte markeder nærmere behandlet.<br />
6.1 Små fjernvarmecentraler til boligopvarmning<br />
Et lille kraftvarmeanlæg baseret på en 35 kW stirlingmotor vil typisk erstatte en oliekedel eller<br />
en kedel til flis eller træpiller med en varmeeffekt på 400 – 1000 kW. Stirling
PSO-FU2204 64 9 kW Stirlingmotor<br />
kraftvarmeanlægget vil ikke alene sælge varme lavet på basis af flis eller træpiller men også<br />
el, der sælges til nettet eller direkte til kunderne, der er tilknyttet fjernvarmesystemet.<br />
I Danmark er der ca. 2,4 mill. boliger, hvoraf ca. halvdelen er parcelhuse, rækkehuse eller<br />
stuehuse på landet. Ca. 10% af disse boliger kunne med <strong>for</strong>del tilsluttes et stirlinganlæg, der<br />
producerede varme til opvarmning af boligerne og el til nettet. Det svarer til ca. 120.000<br />
boliger, hvilket svarer til ca. 2.400 anlæg med en el-effekt på 35 kW eller 10.000 anlæg med<br />
en el-ydelse på 9 kW.<br />
I Tyskland, Østrig, Schweiz, Storbritanien og Skandinavien er der ca. 80 mill. boliger, hvilket<br />
med de samme antagelser svarer til ca. 80.000 anlæg med en el-effekt på 35 kW eller<br />
240.000 anlæg med en el-ydelse på 9 kW.<br />
Det internationale marked vurderes til at være meget stort. I Tyskland, Østrig og Schweiz er<br />
interessen allerede nu meget stor, og potentialet i lande som f.eks. Brasilien, Argentina,<br />
Filippinerne og Indonesien skønnes at være endog meget stort.<br />
EU-kommissionen har anslået, at man alene i EU i de kommende år vil investere 1000<br />
milliarder DKr. i anlæg til produktion af el og varme på basis af sol, vind, biomasse og<br />
geotermi. Her anses biomasse at blive det største område.<br />
6.2 Virksomheder <strong>for</strong>trinsvis i træindustrien<br />
Næste alle virksomheder i træindustien har et betydeligt behov <strong>for</strong> el til drift af maskinerne og<br />
varme til opvarmning af produktionslokaler og administration samt til tørring af<br />
råmaterialerne. Varmen produceres ofte ved hjælp af en kedel, der <strong>for</strong>brænder høvlspåner,<br />
træsmuld eller andre affaldsprodukter fra produktionen. El købes fra el-nettet.<br />
I en undersøgelse af markedet finansieret af EU angives det, at der i Europa er ca. 224 000<br />
virksomheder inden<strong>for</strong> træindustrien. 94% af disse virksomheder beskæftiger mindre end 20<br />
mennesker. Kun 5% har 100 ansatte eller flere. Den samlede omsætning er ca. 1000<br />
milliarder DKr. Mange af disse virksomheder vil med <strong>for</strong>del kunne installere en el-generator<br />
eller et kraftvarmeanlæg, der anvender spildprodukter fra produktionen som brændsel.<br />
Størrelsen <strong>for</strong> langt hovedparten af virksomhederne vil være 10 – 100 kW-el. Alene i Østrig<br />
skønnes det, at der er et behov <strong>for</strong> 11 000 anlæg med en effekt på 100 kW.<br />
6.3 Landbrug<br />
I Danmark er der 110.000 personer beskæftiget i landbruget, skovbruget og på gartnerierne.<br />
Mange af disse landbrug vil med <strong>for</strong>del have en lille 9 kW stirlingmotor til at producere el og<br />
varme til produktionen og boligen.<br />
6.4 Støtteordninger i Europa<br />
I de første 5 – 10 år vil det være nødvendigt at el-prisen er støttet politisk, hvis økonomien i<br />
anlæggene skal være tilstrækkeligt attraktiv. Det <strong>for</strong>ventes imidlertid, at serieproduktion i<br />
større serier på 1000 – 5000 stk. pr. år vil kunne bringe prisen så langt ned, at anlæggene<br />
ikke vil kræve politiske bestemte støtteordninger <strong>for</strong> at være rentable.<br />
I Tyskland har man indført regler, der medfører at man får mindst 0,19 EURO pr. kWh <strong>for</strong> el<br />
produceret med biomasse som brændsel på anlæg med en el-effekt under 50 kW. Der er<br />
ikke andre teknologier, der bare er i nærheden af kommercialisering, så det er særdeles<br />
vigtigt at sætte fart på kommercialiseringen, mens der er en meget stor interesse blandt<br />
potentielle kunder.
PSO-FU2204 65 9 kW Stirlingmotor<br />
I Østrig er der indført regler, der er næsten lige så lukrative som i Tyskland, om end reglerne<br />
er noget mere komplicerede.<br />
I Danmark er reglerne <strong>for</strong> tilskud til kraftvarmeanlæg, der anvender bio-brændsler, desværre<br />
ikke så attraktive. Ved salg af el til nettet får man 60 øre pr. kWh, men hvis anlægget har en<br />
el-effekt under 150 kW kan man producere el til eget <strong>for</strong>brug og dermed erstatte<br />
in<strong>dk</strong>øbsprisen <strong>for</strong> el, som <strong>for</strong> både private og institutionelle <strong>for</strong>brugere er meget høj (op til 1,5<br />
kr pr. kWh).<br />
6.5 Det oversøiske marked<br />
På det oversøiske marked skal prisen helt ned på ca. 12 000 DKr pr. installeret kW-el, hvis<br />
anlægget skal være salgbart i store styktal. Firmaet Shell har udvidet sine<br />
<strong>for</strong>retningsområder med et femte inden<strong>for</strong> biomasse. Shell skønner, at markedet i EU er lille<br />
sammenlignet med Sydamerika og Sydøst Asien. Der<strong>for</strong> har Shell købt skov gennem de<br />
seneste år, og Shell er i dag mellem de 10 største skovejere i verden. Også Statoil og Hydro-<br />
Texaco har iværksat aktiviteter inden<strong>for</strong> bio-energi, og Statoil er interesseret i at bidrage til<br />
finansiering af opbygningen af et stirlinganlæg, der anvender træpiller som brændsel.<br />
Et detaljeret kendskab til markedet <strong>for</strong> stirlinganlæg er vanskeligt at opbygge, <strong>for</strong>di markedet<br />
ikke eksisterer i dag, og <strong>for</strong>di energimarkedet er domineret af politisk styring. Der er der<strong>for</strong> i<br />
høj grad tale om pionerteknologi.<br />
6.6 Konkurrerende produkter<br />
Den nærmeste konkurrent til stirlingmotoren i små kraftvarmeanlæg til bio-brændsler er<br />
<strong>for</strong>gasning og efterfølgende rensning af gassen, så den kan anvendes i en almindelig<br />
<strong>for</strong>brændingsmotor. Kommercielle anlæg er ved at komme på markedet, men anlæggene er<br />
<strong>for</strong>holdsvis komplicerede. Behovet <strong>for</strong> personale til løbende vedligehold er større<br />
sammenlignet med Stirling-anlæggene, og anlæggene skal serviceres betydeligt oftere.<br />
Yderligere kan der <strong>for</strong> visse af teknologierne være problemer med bortskaffelse af de<br />
restprodukter, der ofte er en følge af gasrensningsprocessen.<br />
En anden konkurrerende teknologi er dampturbiner eller dampmotorer. Ved el-effekter under<br />
2000 kW er disse anlægstyper dyre og virkningsgraden er lav. Desuden kræves uddannet<br />
personale til at passe dampkedelen, så udgifterne til bemanding og service er høje<br />
sammenlignet med stirlingmotorerne.<br />
Det konkluderes, at der er et betydeligt marked <strong>for</strong> små stirlingkraftvarmeværker til<br />
biobrændsler med en el-effekt fra 2 kW til 150 kW, og at der ikke ser ud til at komme<br />
konkurrence fra andre producenter af stirlingmotorer eller andre teknologier på kort sigt.
PSO-FU2204 66 9 kW Stirlingmotor<br />
7 Konklusion og <strong>for</strong>slag til videre arbejde<br />
En omfattende revision af en eksisterende 9 kW stirlingmotor er gennemført, og en ny motor<br />
er fremstillet. Resultatet af udviklings<strong>projektet</strong> er en gennemgribende rekonstruktion af<br />
motorens design, som <strong>for</strong>ventes at være væsentlig <strong>for</strong>bedret på en række punkter. Bl.a. er<br />
maksimalkraften på lejerne reduceret med 10%, bortset fra lejerne på arbejdsstemplets<br />
<strong>for</strong>bindelse til åget. For de anvendte nåle- og rullelejer betyder en 10% reduktion af<br />
belastningen en <strong>for</strong>øgelse af levetiden på næsten 40%. Den nye motor er først afprøvet efter<br />
den <strong>for</strong>melle afslutning af <strong>projektet</strong>, og resultaterne fra afprøvningen vedlægges som bilag 8.<br />
Der er opbygget et dataopsamlingssystem, som er anvendt til afprøvning af en af de<br />
eksisterende 9 kW motorer, der blev renoveret til <strong>for</strong>målet. De nye laboratoriefaciliteter og<br />
den nye motor vil fremover blive anvendt til laboratorie<strong>for</strong>søg med henblik på at undersøge<br />
<strong>for</strong>skellige fænomener samt komponenter og materialer.<br />
Et numerisk simuleringsprogram til analyse af stirlingmotorers egenskaber er modificeret og<br />
<strong>for</strong>bedret med hensyn til beregning af tab i regeneratoren samt tryktab i motorens <strong>for</strong>skellige<br />
varmevekslere. Parallelt med arbejdet i nærværende projekt er der udviklet en meget<br />
detaljeret simuleringsmodel til analyse af detaljerne i stirlingmotorens funktion. Denne<br />
udvikling er gennemført som et Ph.D. projekt, og de samlede resultater fra de to<br />
simuleringsmodeller har ført til en væsentlig bedre <strong>for</strong>ståelse af tabsmekanismer i<br />
stirlingmotorer.<br />
Ved sammenligning med måledata er simuleringsprogrammerne verificeret, og resultaterne<br />
viser en meget tilfredsstillende overensstemmelse mellem beregninger og laboratorie<strong>for</strong>søg.<br />
I modsætning til de fleste andre modeller til simulering af stirlingprocessen er der i disse<br />
programmer ingen ”tilpasninger” til måledata, men alene løsning af de termodynamiske<br />
ligninger, der beskriver processerne, samt empiriske udtryk til beregning af<br />
varmetransmission og trykfald hentet fra litteraturen eller fra laboratorie<strong>for</strong>søg. De to<br />
modeller udgør i dag henholdsvis det hurtigste og det mest detaljerede<br />
simuleringsprogrammer til detaljerede analyser og optimering af stirlingmotorer. Fremtidigt<br />
arbejde inden<strong>for</strong> simulering af stirlingmotorer kan med <strong>for</strong>del udvides med CFD-beregninger<br />
(Computational Fluid Dynamics) af strømningen i lokale områder af stirlingmotoren <strong>for</strong><br />
optimering af delkomponenter og overgange mellem komponenter<br />
De nye og reviderede simuleringsværktøjer er bl.a. anvendt til undersøgelse af tabene i<br />
spalten mellem cylindervæg og stempel i stirlingmotorens varme del. Undersøgelsen viste, at<br />
spalten mellem cylinder og <strong>for</strong>trængerstempel hidtil er beregnet på et <strong>for</strong>kert grundlag, og at<br />
spalten enten skal være større eller mindre <strong>for</strong> at opnå et mindre spaltetab.<br />
Måledata fra laboratorie<strong>for</strong>søg med bl.a. Nitrogen som arbejdsgas har givet interessante<br />
in<strong>for</strong>mationer om motorernes fuldlast og dellast egenskaber, samtidig med at <strong>for</strong>skellige<br />
typer smørefedt til smøring af rulle- og nålelejerne samt slid på de smøringsfri materialer er<br />
undersøgt.<br />
Indledende afprøvning af et smørefedt baseret på en ny type syntetisk olie (PFPE) har<br />
<strong>for</strong>eløbigt vist lovende resultater, og det <strong>for</strong>ventes at det er muligt at <strong>for</strong>længe intervallet<br />
mellem eftersmøring med en faktor 5. Også et af de nye smøringsfri stempelringsmaterialer,<br />
der er afprøvet, har væsentlig <strong>for</strong>bedrede egenskaber i <strong>for</strong>hold til det tidligere anvendte.<br />
Forsøgene med en regenerator fremstillet ved sintring af tynde fibre viste, at<br />
fiberregeneratoren havde tilnærmelsesvis samme egenskaber som de regeneratorer, som
PSO-FU2204 67 9 kW Stirlingmotor<br />
allerede anvendes i motorerne. Undersøgelsen blev dog påvirket af, at løse fibre flød rundt i<br />
motoren med arbejdsgassen og satte sig i de PTFE-baserede stempelringe. Tætningerne<br />
virker dårligere med fibre indlejret i overfladen, hvilket indvirkede på målingerne. Desuden<br />
satte der sig fibre i den ventil, der regulerede motorens middeltryk i <strong>for</strong>hold til trykket i<br />
krumtaphuset. Problemet med løse fibre må nødvendigvis løses, inden det er muligt at<br />
anvende fiberregeneratorerne i praksis.<br />
Alternativer til de eksisterende PTFE-baserede tætningers materiale og ud<strong>for</strong>mning er<br />
undersøgt, og <strong>for</strong>slag til <strong>for</strong>bedringer identificeret og afprøvet. Tætninger til at holde<br />
smørefedt inde i lejerne er også undersøgt, men der er ikke fundet bedre alternativer selvom<br />
behovet er betragteligt.<br />
Et syntetisk smørefedt fra Klüber er undersøgt, og potentialet <strong>for</strong> <strong>for</strong>bedring af driftstid<br />
mellem service samt lejelevetid er betydeligt. Den eneste ulempe ser ud til at være den høje<br />
pris, men da mængden af smørefedt i motoren er lille, er det ikke en hindring <strong>for</strong> anvendelse<br />
af de syntetiske fedt-typer.<br />
I en artikel præsenteret ved en konference i Durham, UK, er de særlige problemer ved<br />
anvendelse af biomasse som brændstof til stirlingmotorer behandlet. Det kan konkluderes, at<br />
der er en sammenhæng mellem den hydrauliske diameter <strong>for</strong> den varme varmevekslers<br />
ud<strong>for</strong>mning på røggassiden og stirlingmotorens middeltryk. Ved direkte <strong>for</strong>brænding af<br />
biomasse er det maksimale middeltryk således ca. 50 bar.<br />
En undersøgelse af <strong>for</strong>delen ved at anvende brint som arbejdsgas i stedet <strong>for</strong> Helium er<br />
gennemført. Ud fra en teoretisk betragtning kan man konkludere, at brint er en bedre<br />
arbejdsgas til stirlingmotorer end Helium. Fordelen ved at anvende brint afhænger imidlertid<br />
af den aktuelle motors konstruktion. Sikkerhedsproblemer og usikkerhed vedrørende brints<br />
eventuelle skadevirkninger på lejerne medfører imidlertid, at det ikke kan anbefales på kort<br />
sigt at anvende brint i stedet <strong>for</strong> Helium.<br />
En vurdering af stirlingteknologien og markedet viser, at der er et betydeligt potentiale <strong>for</strong><br />
afsætning af små anlæg til produktion af el og varme på basis af flis og andre <strong>for</strong>mer <strong>for</strong><br />
biomasse.<br />
Der er <strong>for</strong>tsat et stort behov <strong>for</strong> videreudvikling af teknologien, hvis den danske førerposition<br />
inden<strong>for</strong> stirlingmotorer til små kraftvarmeanlæg med biomasse skal bevares. Der er behov<br />
<strong>for</strong> at udbygge anvendelsen af de udviklede simuleringsprogrammer til optimeringer af<br />
stirlingmotorerne funktion, og et mere systematisk undersøgelse af smøringsfri materialer til<br />
stempelringe er højt prioriteret. Anvendelse af FEM (Finite Element Methode) til analyse af<br />
spændinger og termiske <strong>for</strong>hold i <strong>for</strong>bindelse med krumtapmekanismens lejer samt CFD<br />
(Computational Fluid Mechanics) til analyse af strømnings<strong>for</strong>holdene i stirlingmotorens<br />
komponenter <strong>for</strong>ventes også at kunne bidrage til bedre per<strong>for</strong>mance, længere levetid og<br />
reduceret støj.
PSO-FU2204 68 9 kW Stirlingmotor<br />
Referencer<br />
1. Henrik Carlsen, Jonas Bovin, ”9 kW Stirling Engine <strong>for</strong> Biogas and Natural Gas”,<br />
Energistyrelsen Jr.nr. 5116/96-0011<br />
2. Urieli, I., Berchowitz,”Stirling Cycle Engine Analysis”, Adam Hilger LtD, 1984,<br />
ISBN 0-85274-435-8<br />
3. G. Niemann: Maschinen-elemente. Springer-Verlag 1981. B1