PRŮRAZ VZDUCHOVÉ MEZERY MEZI ELEKTRODAMI GENERÁTORU NÍZKOTEPLOTNÍ PLAZMY (reg.číslo GAČR 101/05/0643)
Ing. Libor Fiala
Generátor nízkoteplotní plazmy je nejmodernější zařízení pro roztápění a stabilizaci uhelných energetických bloků v uhelných elektrárnách s granulačním či výtavným ohništěm. Pro prvotní inicializaci plazmové pochodně jiskrovým výbojem se používá vysokonapěťový a vysokofrekvenční ionizační zdroj. Tento zdroj má velmi nepříznivé vlastnosti ovlivňují funkci celé technologie. Z tohoto důvodu byla vytvořena idea nahrazení ionizačního zdroje při prvotní ionizaci generátoru nízkoteplotní plazmy metodou přibližování elektrod k sobě až do okamžiku vzniku elektrického oblouku vlivem průrazu vzduchu v mezielektrodovém prostoru intenzitou elektrického pole. Pro ověření této metody ionizace byla provedena matematická simulace elektrického pole v mezielektrodovém prostoru.
Generator of low-temperature plasma is state-of-the-art equipment for smelting and stabilization of coal generating blocks in coal power stations with granulating or smelting fireplace. The high-voltage and high-frequency ionisation source has been used for primary initialization of plasma cresset with spark discharge. This source has very adversely features influencing functions of the whole technology. On this account it was created the idea of supplying ionisation source for primary ionizing plasma cresset, method approximation electrode on till moment rise electric arc owing to electric intensity becomes disruptive discharge air in space between electrode. The mathematic simulation of electric field in space between electrode was performed for checking of those method of ionisation. 13. ANSYS Users’ Meeting, 21. – 23.září 2005 Přerov -1-
Plazmová technologie představuje nový trend ve způsobu iniciace hoření uhelné směsi při roztápění kotlů ze studeného či teplého stavu, dále je využitelná pro stabilizaci kotle. Tento způsob roztápění kotlů je založen na iniciaci hoření uhelné směsi proudem nízkoteplotního plazmatu využívaného jako počátečního aktivačního či stabilizačního zdroje. J e d i n ý m p a l i v e m využívaným pro roztápění kotle touto nejmodernější technologií je uhelný prášek aktivovaný proudem nízkoteplotního plazmatu vytékající z anodové části plazmatronu, který šetří provozní náklady na roztápění kotle oproti mazutu či zemnímu plynu. Generátor nízkoteplotní plazmy je elektrické zařízení principielně se skládající z dvou elektrodových částí (katodové a anodové) a oddělovací i z o l a č n í č á s t i . K a ž d á e l e k t r o d o v á č á s t s e s k l á d á z v l a s t n í mě d ě n é elektrody a z vlastního ocelového obalu. Potřebnými medii pro provoz plazmového generátoru jsou tlakový vzduch a voda pro chlazení měděných elektrod. Každá elektroda generátoru plazmy má vlastní chladící okruh integrovaný do katodového a anodového obalu a chladící voda je dodávána z externího zdroje vody. T l a k o v é v z d u c h y m a j í z a ú k o l f o r mo v a t o b l o u k v o s e p l a z m a t r o n u j a k o nosné médium a dále unášet elektrický oblouk přes výstupní anodu difuzorového provedení do aktivační a iniciační komory. Plnící tlakový vzduch je do generátoru nízkoteplotní plazmy přiváděn ve dvou rozdělených větvích, a to na vzduch proudící do katody přes tří drážkový zavířovač s osovou tryskou a na vzduch proudící do anody přes osmi cestný tryskový zavířovač integrovaný do oddělovací izolační vložky. Na obrázku č.1 je zobrazen osový řez generátorem nízkoteplotní plazmy.
Obr.1 Osový řez generátorem nízkoteplotní plazmy 13. ANSYS Users’ Meeting, 21. – 23.září 2005 Přerov -2-
Nízkoteplotní plazma vzniká působením stejnosměrného elektrického oblouku na zavířenou soustavu plnících tlakových vzduchů v prostoru vymezeném měděnými elektrodami. Plazmatron je napájený stejnosměrným stabilizovaným proudem a dosahuje elektrického výkonu v regulačním rozsahu 160 – 320 kW. Stěžejními podmínkami pro správný chod generátoru nízkoteplotní plazmy a stabilní proces hoření elektrického oblouku v mezielektrodovém prostoru jsou stejnosměrný stabilizovaný proud vztažený k průtočnému množství plnících tlakových vzduchů. Na obrázku č.2 je zobrazen boční pohled na plazmovou pochodeň vystupující z plazmového generátoru.
Obr.2 Boční pohled na plazmovou pochodeň vystupující z generátoru nízkoteplotní plazmy Pro prvotní ionizaci mezielektrodového prostoru je používán vysokonapěťový a vysokofrekvenční ionizační zdroj (oscilátor) s parametry výstupního napětí 15 kV a frekvenci 2 MHz. Tento způsob inicializace elektrického oblouku má několik zásadních nevýhod. Mezi tyto nevýhody patří vysoká hladina elektromagnetického rušení, která má vliv na veškerá elektrická zařízení v okolním prostoru. Vliv elektromagnetického rušení lze eliminovat pouze do stavu, že ionizační zdroj ruší v omezené míře pouze vlastní elektrické zařízení plazmové technologie. Z důvodu, aby toto rušení neovlivňovalo systémy měření a regulace celé technologie byly vyvinuty nezávislé ochrany, které oddělují citlivé elektrické okruhy od zbylých zarušených okruhů. Mezi další nevýhody patří citlivost na kvalitu povrchu elektrod v místech vzniku ionizační elektrické jiskry vyvolané oscilátorem. S touto nevýhodou úzce souvisí nepříliš velká vlastní spolehlivost ionizačního zdroje. Z následujícího plyne, že s uvedeným způsobem inicializace plazmové pochodně elektrickou jiskrou oscilátoru, nelze zaručit 100 % spolehlivost při startu generátoru nízkoteplotní plazmy. P r á v ě z t ě c h t o d ů v o d ů b y l a o b j e v e n a i d e a z mě n y z p ů s o b u i o n i z a c e elektrodového prostoru. Nový způsob zapalování plazmového generátoru 13. ANSYS Users’ Meeting, 21. – 23.září 2005 Přerov -3-
spočívá v přibližování měděných elektrod k sobě, čímž se zmenšuje vzduchová mezera a roste intenzita elektrického pole. V určité vzdálenosti elektrod od sebe dojde k elektrickému průrazu mezielektrodového prostoru a tím ke vzniku elektrického výboje mezi katodou a anodou. Protože plazmová technologie je nasazená na elektrárně Prunéřov I. a tudíž jsme limitováni elektrickými parametry zde nasazeného výkonového tyristorového usměrňovače. Tento usměrňovač dává při zatížení naprázdno napětí max. 560 V. Tato hodnota napětí je přivedena na katodovou elektrodu, druhá anodová elektroda je přizemněna. Při analytickém ověření možnosti použití metody přibližování elektrod pro inicializaci elektrického oblouku jsme vyšli ze známé Paschenovy křivky. Z Paschenovy křivky uvedené na obrázku č.3 a z přibližných výpočtu dle dostupných analytických rovnic pro homogenní pole jsme došli k výsledku, že k průrazu elektrické pole a vzniku elektrického oblouku dojde při vzdálenosti vzduchové mezery přibližně 0,20 mm. Z analytického přiblížení řešeného problému pro homogenní pole jsme došli k přibližné vzduchové mezeře potřebné k přeskoku elektrického oblouku definovaným stejnosměrným napětím. Pro hodnotu tlakových poměrů vzduchu mezi elektrodami a pro hodnotu napětí 560V mezi měděnými elektrodami je přeskoková vzdálenost přibližně 0,2 mm pro homogenní elektrické pole.
-
100
experimentální závislost vypočítaná závislost
10
Up
Up [kV] 1
0,1 0,01
0,1
1
10
100
1000
pd [kPa cm] Obr.3 Závislost přeskokového napětí v homogenném poli pro vzduch
13. ANSYS Users’ Meeting, 21. – 23.září 2005 Přerov -4-
Protože elektrické pole v oblasti elektrod generátoru nízkoteplotní plazmy nelze považovat za homogenní, bylo nutné pro zpřesnění analýzy přistoupit k řešení simulace elektrického pole numerickými metodami. Pro takové řešení analýzy byl použit software ANSYS Emag, který využívá numerické metody konečných prvků. Pro reálnost analýzy byl vytvořen přesný geometrický numerický model oblasti mezi elektrodami plazmatronu. Pro urychlení výpočtu analýzy bylo s výhodou použito osové geometrické souměrnosti a byl tedy vytvořen dvourozměrný model s definovanou osou symetrie (na obrázku č.4 je zobrazen geometrický model oblasti mezi elektrodami plazmatronu, obrázek č.5 přestavuje detailní zobrazení oblasti vzduchového prostředí mezi měděnými elektrodami). Výpočet numerickými metodami provedeme z hlediska spolehlivosti průrazu mezielektrodového prostoru plazmového generátoru pro snížené napětí na katodě o hodnotě 510V. 1
AREAS
JUL 13 2005 09:03:17
TYPE NUM
Y Z
X
Obr.4 Zobrazení vytvořené geometrie modelu prostoru mezi elektrodami generátoru nízkoteplotní plazmy
13. ANSYS Users’ Meeting, 21. – 23.září 2005 Přerov -5-
1
AREAS
JUL 13 2005 09:04:02
TYPE NUM
Obr.5 Detailní zobrazení geometrie modelu prostoru mezi elektrodami generátoru nízkoteplotní plazmy Pro přesnost výpočtu je velmi důležitá kvalitně vytvořená výpočetní síť. Při pokrývání modelu sítí byl dán důraz na vytvoření kvalitní mapované sítě a dále na vhodné nastavení hustoty elementů v různých místech modelu dle předpokládaného rozložení elektrického pole a to podle amplitudy pole a především změn směru pole. V místech s předpokládaným vysokým gradientem pole byla nastavena hustší síť, v místech s malým gradientem jemnější síť. Důraz byl rovněž kladen na síť představující okolní vzduch, vnější plocha vytvořeného modelu byla pokryta pro zpřesnění výpočtu sítí INFIN zohledňující prostup elektrického pole dále do okolního prostoru. Provedení modelu sítě konečných prvků je ukázáno na obrázku č.6.
13. ANSYS Users’ Meeting, 21. – 23.září 2005 Přerov -6-
1
ELEMENTS
JUL 13 2005 09:01:30
Obr.6 Zobrazení sítě konečných prvků řešeného modelu Řešení této numerické analýzy bylo provedeno pro 12 předem definovaných vzdáleností měděných elektrod generátoru nízkoteplotní plazmy. Tyto vzdálenosti byli vybrány dle analytického výpočtu a představují postupné přibližování elektrod dle kroku, který se zmenšuje s klesající vzduchovou mezerou mezi elektrodami. V okamžiku, kdy p ř e k r o č í i n t e n z i t a e l e k t r i c k é h o p o l e h o d n o t u 3 M V / m mů ž e m e d l e teoretických předpokladů předpokládat, že nastane průraz mezielektrodového prostředí intenzitou elektrického pole v důsledku překročení průrazné pevnosti vzduchu. Zobrazení průběhu elektrické intenzity mezi elektrodami ve vybraných velikostech vzduchové mezery jsou uvedeny na obrázcích č.7, 9, 11, 13 a 15. Dále jsou vyobrazeny grafické výstupy závislosti intenzity elektrického pole na vzdálenosti vzduchové mezery v nejužším místě mezielektrodového prostoru mezi měděnými elektrodami generátoru nízkoteplotní plazmy. Tyto závislosti jsou uvedeny na obrázcích č.8, 10, 12, 14 a 16.
13. ANSYS Users’ Meeting, 21. – 23.září 2005 Přerov -7-
1
1
NODAL SOLUTION
JUL 13 2005 12:56:48
POST1 STEP=1 SUB =1 TIME=1 PATH PLOT NOD1=682 NOD2=702 EFSUM
intenzita el. pole [V/m]
STEP=1 SUB =1 TIME=1 EFSUM (AVG) RSYS=0 SMN =2.492 SMX =330626
MX
JUL 13 2005 13:24:39
(x10**2) 3306.256 3282.598 3258.943 3235.288 3211.633 3187.978 3164.323 3140.668 3117.013 3093.358 (x10**-3)
3069.703 0
.324 .162
.648 .486
.972 .81
1.296 1.134
1.623 1.458
vzduchová mezera [m]
2.492
36164
72326
108488
144650
180812
216974
253136
289298
330626
Obr.7 Kontury intenzity el. pole Obr.8 Závislost intenzity el. pole na v mezielektrodovém prostoru pro vzdálenosti vzduchové mezery o vzduchovou mezeru 1,62 mm velikosti 1,62 mm 1
1
NODAL SOLUTION
JUL 13 2005 13:06:26
POST1 STEP=1 SUB =1 TIME=1 PATH PLOT NOD1=690 NOD2=731 EFSUM
intenzita el. pole [V/m]
STEP=1 SUB =1 TIME=1 EFSUM (AVG) RSYS=0 SMN =2.586 SMX =508780
MX
JUL 13 2005 14:13:04
(x10**2) 5069.934 5044.650 5019.362 4994.074 4968.786 4943.498 4918.210 4892.922 4867.634 4842.346 (x10**-3)
4817.058 0
.208 .104
.416 .312
.624 .52
.832 .728
1.043 .936
vzduchová mezera [m]
2.586
55650
111298
166945
222593
278240
333888
389535
445183
508780
Obr.9 Kontury intenzity el. pole Obr.10 Závislost intenzity el. pole v mezielektrodovém prostoru pro na vzdálenosti vzduchové mezery o vzduchovou mezeru 1,04 mm velikosti 1,04 mm 1
1
NODAL SOLUTION
JUL 13 2005 13:32:59
POST1 STEP=1 SUB =1 TIME=1 PATH PLOT NOD1=710 NOD2=792 EFSUM
intenzita el. pole [V/m]
STEP=1 SUB =1 TIME=1 EFSUM (AVG) RSYS=0 SMN =2.657 SMX =.113E+07
MX
JUL 13 2005 14:20:20
(x10**3) 1119.911 1117.246 1114.577 1111.908 1109.239 1106.570 1103.901 1101.232 1098.563 1095.894 (x10**-4)
1093.225 0
.926 .463
1.852 1.389
2.778 2.315
3.704 3.241
4.632 4.167
vzduchová mezera [m]
2.657
123974
247945
371916
495887
619859
743830
867801
991772
.113E+07
Obr.11 Kontury intenzity el. pole Obr.12 Závislost intenzity el. pole v mezielektrodovém prostoru pro na vzdálenosti vzduchové mezery o vzduchovou mezeru 0,46 mm velikosti 0,46 mm 13. ANSYS Users’ Meeting, 21. – 23.září 2005 Přerov -8-
1
1
NODAL SOLUTION
JUL 13 2005 13:48:58
POST1 STEP=1 SUB =1 TIME=1 PATH PLOT NOD1=710 NOD2=792 EFSUM
intenzita el. pole [V/m]
STEP=1 SUB =1 TIME=1 EFSUM (AVG) RSYS=0 SMN =2.676 SMX =.237E+07
MX
JUL 13 2005 14:28:47
(x10**3) 2324.346 2321.005 2317.666 2314.327 2310.988 2307.649 2304.310 2300.971 2297.632 2294.293 (x10**-4)
2290.954 0
.454 .227
.908 .681
1.362 1.135
1.816 1.589
2.265 2.043
vzduchová mezera [m]
2.676
258764
517526
776287
.104E+07
.129E+07
.155E+07
.181E+07
.207E+07
.237E+07
Obr.13 Kontury intenzity el. pole Obr.14 Závislost intenzity el. pole v mezielektrodovém prostoru pro na vzdálenosti vzduchové mezery o vzduchovou mezeru 0,23 mm velikosti 0,23 mm 1
1
NODAL SOLUTION
JUL 13 2005 14:00:52
POST1 STEP=1 SUB =1 TIME=1 PATH PLOT NOD1=690 NOD2=730 EFSUM
intenzita el. pole [V/m]
STEP=1 SUB =1 TIME=1 EFSUM (AVG) RSYS=0 SMN =2.68 SMX =.335E+07
MX
JUL 13 2005 14:02:52
(x10**3) 3346.364 3334.938 3323.508 3312.078 3300.648 3289.218 3277.788 3266.358 3254.928 3243.498 (x10**-4)
3232.068 0
.31 .155
.62 .465
.93 .775
1.24 1.085
1.550 1.395
vzduchová mezera [m]
2.68
366438
732874
.110E+07
.147E+07
.183E+07
.220E+07
.257E+07
.293E+07
.335E+07
Obr.15 Kontury intenzity el. pole Obr.16 Závislost intenzity el. pole v mezielektrodovém prostoru pro na vzdálenosti vzduchové mezery o vzduchovou mezeru 0,15 mm velikosti 0,15 mm Na uvedených obrázcích č.7-16 je názorně vidět, jak se zvyšuje intenzita elektrického pole při snižující se vzdálenosti vzduchové mezeře mezi měděnými elektrodami. Pro větší přehlednost závislosti intenzity elektrického pole na velikosti vzduchové mezery bylo provedeno grafické zhodnocení této závislosti pro všechny vypočtené vzdálenosti vzduchové mezery. V grafickém zobrazení závislosti je zobrazena oblast teoretického průrazu mezielektrodového prostoru překročením průrazné pevnosti vzduchového prostředí mezi elektrodami intenzitou elektrického pole. Tato závislost intenzity elektrického pole na velikosti vzduchové mezery je zobrazena na obrázku č.17.
13. ANSYS Users’ Meeting, 21. – 23.září 2005 Přerov -9-
Intenzita elektrického pole mezielektrodového prostoru v závislosti na velikosti vzduchové mezery
3500 3346
Oblast teoretického průrazu mezielektrodového prostoru intenzitou elektrického pole 3000
Intenzita el.pole - max. hodnoty [kV/m]
2974
2500 2324
2000 1911
1500 1411
1000
1119
795 500
619 507
430
374
330
0 0
0,2
0,4
0,6
0,8
1
1,2
1,4
1,6
1,8
velikost vzduchové mezery [mm]
Obr.17 Závislost intenzity elektrického pole mezielektrodového prostoru na velikosti vzduchové mezery Ze závislosti intenzity elektrického pole na velikosti vzduchové mezery uvedené na obrázku č.17 vyplývá, že pro tvar elektrod plazmového generátoru má tato závislost tvar hyperbolické křivky a hranici teoretického průrazu mezielektrodového prostoru dosahuje při hodnotě vzdálenosti vzduchové mezery 0,17 mm. Tuto definovanou vzdálenost elektrod můžeme považovat za bezpečnou hodnotu vzduchové mezery, při které spolehlivě dojde k přeskoku. Protože jsme pro výpočet numerické analýzy použili snížené napětí, bude tedy minimální vzdálenost vzduchové mezery větší a dle experimentálním pokusům provedených na zařízení plazmové technologie je přeskoková vzdálenost v rozmezí vzdálenosti 0,20 – 0,25 mm. Přeskoková vzdálenost je závislá nejen na použitém napětí mezi elektrodami, ale také na tlaku vzduchového média v mezielektrodovém prostoru ( viz. Paschenova křivka – obrázek č.3). Vzhledem k tomu, že v mezielektrodovém prostoru proudí tlakový vzduch z rozváděcích kanálů katody a zejména anody je zřejmé, že ke zvýšené přeskokové vzdálenosti mezi elektrodami plazmového generátoru oproti numerické analýze přispívá nejenom zvýšené napětí, ale také právě proudění tlakové vzdušiny v mezielektrodovém prostoru.
Příspěvek zpracován v rámci řešení projektu GAČR 101/05/0643 „Využití plazmové technologie v uhelné energetice“.
13. ANSYS Users’ Meeting, 21. – 23.září 2005 Přerov - 10 -