26.07.2013 Views

ELKRAFT Slutrapport for projektet: - Energinet.dk

ELKRAFT Slutrapport for projektet: - Energinet.dk

ELKRAFT Slutrapport for projektet: - Energinet.dk

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!