EXPERIMENTÁLNÍ OVĚŘENÍ VLASTNOSTÍ PROUDĚNÍ V LOPATKOVÉ KASKÁDĚ STŘEDORYCHLOSTNÍHO TUNELU A POTVRZENÍ VÝSLEDKŮ POMOCÍ CFD SIMULACÍ S OHLEDEM NA VLIV DRSNOSTI POVRCHŮ. SVOČ – FST 2015 Bc. David Fenderl Západočeská univerzita v Plzni, Univerzitní 8, 306 14 Plzeň Česká republika ABSTRAKT Cílem práce bylo získat představu o nárůstu ztrát při expanzi media se zvyšující se drsností lopatky. Vznikající ztráta je úzce spjata s turbulencí. Proudění bylo zkoumáno na experimentální lopatkové mříži. Rovněž byla řešena úloha pomocí CFD. Pro získání okrajových podmínek bylo nutné experimentální ověření vstupní podmínky. Snažení se směřovalo k získání profilových a okrajových ztrát. Po získání dostatečných informací o proudovém poli bylo možné ovlivnit okrajovou ztrátu smykovou vrstvou u povrchu lopatky. ABSTRACT The aim of this article was to get an idea about the losses of increase blade roughness during the expansion. Emerging loss is closely associated with turbulence. The flow was examined on experimental blade cascade. Also, the task was solved using CFD. To obtain boundary conditions was necessary experimental verification of input conditions. Efforts are directed to obtain of profile losses and secondary losses. After obtaining sufficient information on the flow field we may be able to influence secondary losses of shear layer at the surface of the blade. KLÍČOVÁ SLOVA Lopatková mříž, lopatky VS33b a VS33T, drsnost povrchu, intenzita turbulence, rychlostní profil 1.
ÚVOD
Během provozu parní turbíny dochází k erozivnímu a korozivnímu rozrušování lopatek. Ty mají poté vliv na termodynamickou účinnost stupně, objektivně na celý díl parní turbíny a tedy na účinnost celého turbosoustrojí. Zda tento vliv je kladný nebo záporný záleží na typu lopatky a na okrajových podmínkách. Pro tento případ budou řešeny vysokotlaké lopatky založené na profilu VS33. Jedná se o rozváděcí lopatky, kde pracovní médium expanduje v subsonické oblasti. V této publikaci budou řešeny dva typy modelových lopatek s uvažováním nezbytných podobnostních čísel. Hlavní náplní je získat závislost ztrátového součinitele na zvyšující se drsnosti povrchu. Následně získat vstupní okrajovou podmínku na vstupu modelové mříže, a to experimentálně i pomocí CFD. Zde je sledován rychlostní profil, především oblast mezní vrstvy u stěny. Hlavním důvodem ztrát je vznik turbulentní mezní vrstvy na lopatce a sekundárních vírů od stěn lopatkové mříže, ty byly sledovány pomocí aparatury CTA anemometru. Snahou bylo rovněž porovnat experiment se CFD výpočty z pohledu turbulence a vzniku ztrát. K výpočtu byl použit i model rozšířený dvěma transportními rovnicemi pro přechod z laminární do turbulentní mezní vrstvy. V této publikaci se autor zaměří spíše na výsledky, potřebná teoretická část je uvedena ve stejnojmenné diplomové práci.[8]
2.
VSTUPNÍ OKRAJOVÁ PODMÍNKA
Okrajová podmínka na vstupu do lopatkové mříže byla zkoumána ze dvou hlavních důvodů. Prvním z nich bylo nutné ověřit vyrovnanost rychlostního profilu a získat informaci o ztrátě v mezní vrstvě u stěny. Druhým důvodem bylo získat okrajovou podmínku včetně chování turbulizujícího media pro namodelovíní v CFD.
Nejprve byl rychlostní profil proměřen pomocí Pitotovy sondy traverzované do poloviny šířky kanálu. Střední hodnoty rychlostí byly proloženy regresní funkcí. Pro přesnější výsledky u stěny byla aplikována korekční funkce u stěny. Vlivem přiblížení Pitotovy sondy ke stěně dochází k ovlivnění proudového pole a odklonu proudnic od stěny. Následkem tohoto efektu je zkreslení rychlostního profilu. Vše je patrné na Obr. 1, kde je vidět i energetická tloušťka mezní vrstvy.
Obr. 1 Rychlostní profil na vstupu od stěny – Pitototva sonda
Následně byl ten samý rychlostní profil proměřen pomocí měřící aparatury CTA anemometru. Pomocí tohoto zařízení se vzhledem k citlivosti zařízení zjistila určitá nevyrovnanost rychlostního profilu. Hlavním důvodem využití této metody bylo zjištění intenzity turbulence na vstupu. Vše je patrné na Obr. 2. Střední hodnota intenzity turbulence činí 1,856%.
Obr. 2 Rychlostní profil od stěny – CTA aparatura
Pro měření byla použita drátková sonda, viz Obr. 3, pro kterou byla sestavena citlivá traverzovací aparatura. Rovněž bylo nutné vyztužit vstupní kanál a snížit vlastní frekvenci stěn.
Obr. 3 Drátková sonda CTA anemometru
Rovněž byla sledována rychlostní mezní pomocí CFD výpočtu odpovídající experimentu. Při dikretizování výpočtové oblasti byla především sledována jemnost buněk u stěny pro sledování mezní vrstvy, viz Obr. 4.
Obr. 4 Výpočtová oblast pro sledování rychlostního profilu
Po nastavení okrajových podmínek odpovídajících experimentu byl proveden výpočet s užitím dvourovnicového turbulentního modelu k-ω SST. Pro tento model byla rovněž využita intenzita turbulence na vstupu získaná experimentálně. Stupeň turbulence, který je při izotropním proudění roven intenzitě turbulence, je uveden v rovnici (1). Tuto rovnici bylo nutné vložit do programu Fluent jako externí funkci. Fluent využívá pro výpočet intenzity turbulence sofistikovanějšího vztahu, který v mém modelu neodpovídal realitě. (1)
Pole izočar stupně turbulence s definovaným rozsahem je vykreslen na Obr. 5. Výstupní okrajová podmínka je náběžnou hranou lopatek, která byla laděna podle rychlostního profilu ve sledované oblasti proměřovaného rychlostního profilu.
Obr. 5 Stupeň turbulence – definovaný rozsah
Na Obr. 6 je vykreslený rychlostní profil v oblasti Pitotovy a CTA sondy odpovídající experimentu (modře). Červeně je vykreslen vývoj rychlostního profilu u náběžné hrany lopatky, tato ztráta v mezní vzroste na 2,23%. V následující tabulce, Tab. 1, jsou uvedeny ztráty v mezní vrstvě vycházející z energetické tloušťky mezní vrstvy.
Obr. 6 Rychlostní profil na vstupu – CFD
Tloušťky mezních vrstev [mm] Způsob měření Pošinovací tloušťka mezní vrstvy Impulsní tloušťka mezní vrstvy Energetická tloušťka mezní vrstvy Ztráta v mezní vrstvě [%]
Pitotova sonda 0.311504 0.157629 0.280058
CTA aparatura 0.593986 0.272658 0.499441
CFD 0.593944 0.310075 0.522714
0.373411
0.665922
0.696952
Tab. 1 Tloušťky mezních vrstev
3.
INTENZITA TURBULENCE V DEFINOVANÉ ROVINĚ ZA ODTOKOVÝMI HRANAMI LOPATEK
Po vyhodnocení vstupní oblasti lopatkové mříže byly zkoumány dva typy lopatek. Jedná se o lopatky VS33b a VS33T, které jsou specifikovány v [8]. Intenzita turbulence byla zkoumána ve vdálenosti 20% délky tětivy lopatkového profilu za rovinou odtokových hran lopatek. Na následujích dvou grafech, Obr. 7 a Obr. 8, jsou vykreslena intenzita turbulence po rozteči lopatkové mříže s lopatkami VS33T. V prvním grafu je vykreslena intenzita v polovině výšky lopatky a ve druhém 20 mm od stěny, v oblasti maximálních ztrát v sekundárním víru.
Obr. 7 VS33T – Intenzita turbulence ve vzdálenost l/2 Obr. 8 VS33T – Intenzita turbulence ve vzdálenosti 20 mm od stěny
Rovněž byl sledován vliv délky časového záznamu dat na střední hodnotu intenzity turbulence. Jako nejvhodnější se ukázala hodnota 1 sekundy při vzorkovací frekvenci 64 kHz. V Tab. 2 jsou pak uvedeny střední hodnoty turbulence. Vzdálenost l/2 Vzdálenost 20 mm od stěny Střední hodnota intenzity 2.3139 2.8500 turbulence [%] Tab. 2 Střední hodnota intenzity turbulence
4.
ZTRÁTA PŘI EXPANZI VE SLEDOVANÉ ROVINĚNA VÝSTUPU Z LOPATKOVÉ MŘÍŽE – EXPERIMENT
Dalším úkolem bylo sledování ztrát, které by bylo možné následně porovnat s výpočty ze CFD. V následujích dvou tabulkách je uvedeno rozložení ztrátového součinitele při Re2is = 8,8x105 a M2is = 0,24 pro oba typy lopatek. Profilové ztráty jsou patrné z prostředního grafu. Úhrnný výstupní úhel proudu z lopatkové mříže .
Tab. 3 VS33b – Rozložení ztrát v rovině 0.2b, ztráty po rozteči v polovině výšky lopatky, ztráty po rozteči ve vzdálenosti od stěny
Tab. 4 VS33T - Rozložení ztrát v rovině 0.2b, ztráty po rozteči v polovině výšky lopatky, ztráty po rozteči ve vzdálenosti od stěny
5.
ZTRÁTA PŘI EXPANZI VE SLEDOVANÉ ROVINĚ NA VÝSTUPU Z LOPATKOVÉ MŘÍŽE – CFD
Následně byly vytvořeny výpočtové oblasti odpovídající experimentálním lopatkovým mřížím obou lopatek. Kvalita diskretizované oblasti v polovině výšky lopatky je patrná na Chyba! Nenalezen zdroj odkazů.. V polovině výšky lopatky jsou u obou lopatek mříže schodné, proto byly vyhodnoceny současně. Na Obr. 10 je vykreslen stupeň turbulence dle vztahu (1).
Obr. 9 Pohled na síť v okolí profilu v polovině výšky lopatky Obr. 10 Stupeň turbulence v polovině výšky lopatky
Dále byla výpočtová oblast řešena pomocí RANS matenatického modelu s využitím modelu turbulence k-ω SST a rovněž tohoto modelu rozšířeného o gama-Re přechodu do turbulence. Na následujících dvou konturách turbulentní kinetické energie, Obr. 11 a Obr. 12, je patrný rozdíl přechodu z laminární do turbulentní mezní vrstvy na podtlakové straně. Proto při vyhodnocení intenzity turbulence a ztrát na celkovém tlaku je vyhodnocení provedeno pro oba modely.
Obr. 11 Pole izočar turbulentní kinetické energie K, model K-ω SST
Obr. 12 Pole izočar turbulentní kinetické energie K, model K-ω SST včetně γ-Re přechodu
Pole izočar ztrátového součinitele pro jeden mezilopatkový kanál obou lopatkových mříží je vykreslen na Obr. 13 a Obr. 14. Zde je vidět kromě rozložení ztrát i přechod do oblasti malých záporných ztrát, tento fakt způsobený numerikou je třeba brát na vědomí.
Obr. 13 VS33b - izočáry místního ztrátového součinitele celkového tlaku
Obr. 14 VS33T - izočáry místního ztrátového součinitele celkového tlaku
Z výsledků by rád autor ještě uvedl rozložení stupně turbulence po rozteči s užitím obou turbulentních modelů, Obr. 15 a Obr. 16, které jsou porovnatelné z experimentálními daty. Stupeň turbulence získaný numericky je na nižších hodnotách než intenzita získaná experimentálně.
Obr. 15 VS33T - stupeň turbulence Tu na rozteči model K-ω SST
6.
Obr. 16 VS33T - stupeň turbulence Tu na rozteči model K-ω SST včetně γ-Re přechodu
VLIV DRSNOSTI POVRCHU LOPATKY NA PROFILOVÉ A OKRAJOVÉ ZTRÁTY
Zvýšení drsnosti povrchu lopatky bylo uskutečněno pomocí brusiva nanášeného na povrch lopatky pomocí barvy ředitelné nitridovým rozpouštědlem. Pro samotné měření bylo vybráno pět drsností charakteristické velikostí zrna brusiva, viz Tab. 5. Jedna z použitých drsností je patrná na Obr. 17. Střední aritmetické odchylky drsnosti povrchu byla proměřena drsnoměrem. Lopatky bez zdrsnění mají ~Ra 0,4.
Obr. 17 Nanesená drsnost
Norma FEPA F320 F180 F100 F70 F40
Velikost zrna[μm] 45 - 32 90 - 75 150 - 125 250 - 212 500 - 425
~Ra[μm] 13 20 30 36 40
Tab. 5 Úrovně zdrsnění lopatky
Celkový ztrátový součinitel v rovině 0,2b byl proměřen pro všechny výše uvedené úrovně zdrsnění. Na následujících dvou grafech, Obr. 18 a Obr. 19, je vykreslen celkový ztrátový součinitel vážený hmotnostním tokem. Jsou zde vidět profilové ztráty (červeně) a tyto ztráty společné s okrajovými (modře).
Celkový ztrátový součinitel ζM vážený hmotnostním tokem v rovině 0,2b
18,00 17,00 16,00 15,00 14,00 13,00 12,00 11,00 10,00 9,00 8,00 7,00 6,26 6,00 4,41 5,00 4,00 3,00 2,00 1,00 0,00 0
VS33T 13,50 11,10
11,50
10,10
10,00
7,92
8,31
8,78
6,95
6,33
cela_oblast po_rozteci
5
10
15
20
25
30
35
40
45
Ra[μm]
Celkový ztrátový součinitel ζM vážený hmotnostním tokem v rovině 0,2b
Obr. 18 VS33b - celkový ztrátový součinitel ζMC vážený hmotnostním tokem v rovině 0,2b
18,00 17,00 16,00 15,00 14,00 13,00 12,00 11,00 10,00 9,00 8,00 7,21 7,00 6,00 4,47 5,00 4,00 3,00 2,00 1,00 0,00 0
16,20
VS33b 13,10 10,80
13,40
12,70
10,80 9,59
8,73
8,51 cela_oblast
6,30
po_rozteci
5
10
15
20
25
30
35
40
45
Obr. 19 VS33T - celkový ztrátový součinitel ζMC vážený hmotnostním tokem v rovině 0,2b
Ra[μm]
7.
ZÁVĚR
Bylo ověřeno proudové pole experimentální lopatkové mříže VS33b a VS33T včetně vstupní okrajové podmínky a mezních vrstev. Získala se závislost profilových a okrajových ztrát na úrovni zdrsnění povrchu lopatky. Rovněž byly výsledky porovnány s výpočty ze CFD se sledováním turbulence a následných ztrát. V stejnojmenné diplomové práci bude rovněž uvedena závislost profilových a okrajových ztrát na úrovni zdrsnění lopatky užitím turbulentního modelu s gama-Re modelem přechodu. Dále pak ovlivnění a případné snížení okrajových ztrát smykovou vrstvou na povrchu lopatky.
PODĚKOVÁNÍ Zvláštní poděkování patří společnosti DOOSAN ŠKODA POWER s.r.o., ve které byla tato úloha řešena. LITERATURA
[1] ANSYS, „Fluent Theory Guide,“ 2011. [2] M. Hoznedl, J. Mach a J. Matas, „Protokol o zkoušce drsnosti“. [3] K. Sedlák, „Vyhodnocení ztrátových součinitelů přímých lopatkových mříží, měřených na středorychlostním tunelu ŠKODA,“ 2012. [4] G. L. Morrison, M. T. Schobeiri a K. Pappu, „ Five-hole pressure probe analysis technice. Flow Measurement and Instrumentation 9 (1998) 153–158,“ 1998. [5] F. R. Menter, „Two-Equation Eddy-Viscozity Turbulence Models for Engineering Applications,“ AIAA Journal, 1994. [6] J. Linhart, „Mechanika tekutin 2, stručné učební texty, Západočeská univerzita v Plzni,“ 2009. [7] J. Amecke a P. Šafařík, „Data reduction of wake flow measurements with injection of an other gas.,“ 1995. [8] D. Fenderl, „Experimentální ověření vlastností proudění v lopatkové kaskádě středorychlostního tunelu a potvrzení výsledků pomocí CFD simulací s ohledem na vliv drsnosti povrchů.“ 2015