TZ 5 TEPELN ZAZEN Tepeln a hydraulick vpoet

TZ 5 TEPELNÁ ZAŘÍZENÍ Tepelný a hydraulický výpočet výměníků tepla a dimenzování Dimensování S&T výměníků. Tepelný výpočet výměníků tepla metodou LMTD. Rudolf Žitný, Ústav procesní a zpracovatelské techniky ČVUT FS 2010

TZ 5 Shell & Tube – přenos tepla Dali

TZ 5 Shell & Tube – přenos tepla W 2 T 2 ’ Q=k. S T W 1 T 1 ’ W 2 T 2’’ Qloss Tepelná kapacita proudu W [W/K] pro jednofázový tok W 1 T 1’’ Entalpická bilance celého výměníku Entalpická bilance jednoho proudu Tři neznámé při kontrolním výpočtu: výstupní teploty T 1’’, T 2’’ a Q (tepelný výkon). Tepelné ztráty Qloss zpravidla zanedbáme nebo je lze snadno odhadnout. Při návrhovém výpočtu jsou naopak teploty i výkon dané, počítá se teplosměnná plocha. .

TZ 5 Shell & Tube – přenos tepla Entalpická bilance jsou 2 rovnice, ta třetí je kde k je celkový součinitel prostupu tepla počítaný ze seriového řazení termických odporů Střední teplotní diference vyžaduje výpočet teplotních profilů.

TZ 5 Shell & Tube – přenos tepla Tekutina voda pára Stav kapalina var <5 Bar kondenzace 1 Bar W/(m 2. K) 5000 - 7500 3000 - 10000 - 15000 uhlovodík kapalina 0. 5 -2. 5 cp var 0. 5 -2. 5 cp kondenzace 0. 5 -2. 5 cp 750 -1500 1000 - 3500 1500 - 4000 plyn 1 Bar 10 Bar 80 - 125 250 - 400

Shell & Tube – fouling R TZ 5 Tekutina voda pára Stav kapalina var <5 Bar kondenzace 1 Bar R (m 2. K/W) 0. 0001 - 0. 00025 0. 0001 - 0. 0002 0 - 0. 0001 uhlovodík kapalina 0. 5 -2. 5 cp var 0. 5 -2. 5 cp kondenzace 0. 5 -2. 5 cp 0. 0002 - 0. 001 0. 0001 - 0. 0003 plyn 1 Bar 10 Bar 0 - 0. 0001

TZ 5 Shell & Tube – přenos tepla Pouze v případě dvou paralelních toků je možné počítat střední teplotní diferenci metodou LMTD (logarithmic mean temperature difference) T 1’ Souproudý výměník T 1’’ ’=T 1’-T 2’ =T 1 -T 2 T 2’’ T 2’ s ’’=T 1’’-T 2’’ S T 1’’ ’=T 1’-T 2’ =T 1 -T 2 T 2’’ T 2’ Qs ’’=T 1’’-T 2’’ Q ’ ’’ jsou teplotní diference na obou koncích výměníku

TZ 5 Shell & Tube – přenos tepla Důkaz plyne z lineární závislosti mezi teplotou a entalpickým tokem Qs T 1’’ ’=T 1’-T 2’ =T 1 -T 2 Což je obyčejná diferenciální rovnice kterou lze integrovat s výsledkem: T 2’’ T 2’ Qs ’’=T 1’’-T 2’’ Q

TZ 5 Shell & Tube – přenos tepla Přesně totéž platí i protiproudý výměník tepla = důkaz je identický T 1’ T 1 ’=T 1’-T 2’’ atd. , se stejným výsledkem: T 1’’ =T 1 -T 2 T 2’’ T 2 Qs ’’=T 1’’-T 2’ Q Možná to trochu mate, souproudý výměník při stejných vstupních teplotách, stejných průtocích a stejné teplosměnné ploše (k. S) má výrazně menší výkon Q než protiproudý. Jak to, když je vzoreček úplně stejný? Odpověď je v tom, že na obou koncích výměníků nejsou stejné teplotní diference ’ a ’’. Ale hlavně: Jak postupovat, když není jasné co je to „konec výměníku“ a jak potom definovat teplotní diference ’ a ’’? Odpověď je na další stránce…

TZ 5 Shell & Tube – přenos tepla Především je třeba předefinovat pojem středního logaritmického spádu a vyjádřit ho ne rozdílem teplot na koncích výměníku, ale jen rozdíly teplot vstupních a výstupních proudů (tato definice pak bude zcela univerzální) T 1’ ’=T 1’-T 2’’ T 1 T 1’’ =T 1 -T 2 ’’=T 1’’-T 2’ Takto definované LMTD odpovídá protiproudému výměníku tepla. Pro jakýkoliv jiný (neideální) výměník tepla je třeba tuto hodnotu zmenšit, a tomu se říká F-korekce je návrhová metoda, kdy jsou známé všechny teploty i výkon, tedy i P a R. Počítá se teplosměnná plocha k. S.

TZ 5 Shell & Tube – přenos tepla Korekce TLM je funkcí dvou parametrů. Pokud je proud číslo 2 slabší (W 1>W 2, R<1) nazývá se parametr P efektivita (a označuje se symbolem ). Je to poměr tepelného výkonu analyzovaného výměníku k tepelnému výkonu ideálního (protiproudého) výměníku s nekonečně velkou teplosměnnou plochou – v tom případě se totiž výstupní teplota proudu číslo 2 přiblíží vstupní teplotě proudu číslo 1 (T 1’=T 2’’) T 1 ’ T 2’’ T 1 ’ W 1<W 2 T 1’’ T 2’’ T 1’’ T 2 ’ W 1>W 2 T 2 ’

TZ 5 Shell & Tube – přenos tepla F-korekce LMTD R Pro shell & tube výměníky se dvěma a více tahy v trubkách T 1 -teploty v plášti, T 2 -teploty v trubkách

Shell & Tube – přenos tepla TZ 5 F-korekce LMTD Překřížení teplot B A T 1 ’ T 2’’ T 1’’ PA T 2’’ T 1’’ PB FTA T 2 ’ Varianty A, B téhož výměníku se dvěma tahy v trubkách a jedním tahem v plášti lišící se jen obrácením směru proudu v plášti. Varianta B je CHYBA, projeví se překřížením teplotního profilu (ve druhém tahu teplota ohřívaného proudu dokonce klesá). Na efektivitu výměníku to však kupodivu nemá žádný vliv.

TZ 5 Shell & Tube – přenos tepla Excelovský program výpočtu teplotních profilů ve výměníku se dvěma tahy v trubkách Definuj teploty v trubkách (modrý proud) a v plášti (červený proud) NTU 1=k. S/W 1 - tubes, NTU 2=k. S/W 2 - shell

TZ 5 Shell & Tube – tlakové ztráty Tlaková ztráta v trubkách z Darcy Weissbach rovnice Třecí součinitel f závisí na Re a na relativní drsnosti

TZ 5 Deskové VT Femina Tepelné a hydraulické výpočty Použita metodika vyvinutá Martinem Holgerem FEMINA umožňuje modelovat proudění, přenos tepla i hmoty v potrubních sítích. Systém je tvořen propojením 1 D konečných elementů typu trubka (PIPE), výměník (HEX), čerpadlo (PUMP), nádoba (TANK). Power Q [k. W] 80 70 60 50 40 30 20 10 20 30 40 50 60 70 80 chevron angle [deg] 90

TZ 5 Shell & Tube Femina Tepelné a hydraulické výpočty Použita metodika Bell vyvinutá na univerzitě Delaware (Jiří Taborek!)

TZ 5 Bell Delaware (viz TEMA) A-průsak štěrbinou v přepážce B-příčné obtékání C-bypass vně trubkového svazku E-průsak štěrbinou mezi přepážkou a pláštěm J-faktor (Colburn)

Bell Delaware (viz TEMA) TZ 5 Max p plášť Průměr pláště (počet trubek) Požadovaný výkon Příliš nízké rychlosti proudění v trubkách i plášti: riziko foulingu, nízké přestupy tepla Max p trubky Minimální rychlost v trubkách (fouling) Max. rychlost v trubkách (eroze) Vibrace Příliš vysoké rychlosti proudění v trubkách i plášti: vysoké tlakové ztráty, vibrace, eroze Délka trubek

Bell Delaware (viz TEMA) Průměr pláště TZ 5 Nejmenší průměr a nejkratší trubky -nejlevnější řešení Min. Max. Délka trubek

TZ 5 Shell & Tube - dimenzování

Shell & Tube - dimenzování TZ 5 HEDH (Heat Exchanger Design Handbook kapitola 4) a TEMA popisují poměrně detailně strukturní analýzu konstrukčních uzlů výměníku Shell & Tube. Metodika je založena na kategorizaci napětí dle ASME (Boiler and Pressure Vessel Code): idea je v tom, že výpočty stačí provést jen v oblasti lineární pružnosti (to umí každý konečněprvkový program a pro většinu uzlových prvků existují i analytická řešení), a plasticita či únava materiálu se zahrne do vhodných koeficientů bezpečnosti. I. Membránová napětí konstantní po průřezu (např. napětí v plášti od vnitřního přetlaku) II. Ohybová napětí (např. ohybová napětí v trubkovnici zatížené přetlakem – i když dojde k překročení meze kluzu neznamená to nekontrolovatelný nárůst deformace) III. Napětí vyvolaná omezením deformací nebo tepelným zatížením (přechod válcové a kulové skořepiny, plastizace opět neznamená kolaps) IV. Špičková napětí (singularity v ostrých rozích, závitech šroubocých spojů ap) V.

TZ 5 Shell & Tube - dimenzování Z tabulek materiálových parametrů se pro zvolený materiál stanoví Návrhové napětí S=min( u/3, mez pevnosti y/1. 5, mez průtažnosti 2% R/1. 5, T/1. 1) časová mez pevnosti 100000 h mez průtažnosti 1% Vypočtená napětí (lineární analýza) pak musí být dle jejich kategorizace menší než určitý násobek (1, 1. 5, 3) návrhového napětí S, schematicky I<S II<1. 5 S III<3 S IV<SD Únavové napětí SD I+ II<1. 5 S I+ III <3 S

TZ 5 Shell & Tube - dimenzování Membránová napětí patří do kategorie I – nesmí překročit návrhové napětí S Membránové napětí v trubce Membránové napětí v kulové skořepině p s R

TZ 5 Shell & Tube - dimenzování Tato ohybová napětí patří do kategorie II – nesmí překročit napětí 1. 5 S Trubkovnice je kruhová deska s pravidelně rozmístěnými otvory (jejich geometrii určuje rozteč f a průměr d). Počítá se jako plná deska stejné tloušťky, ale se sníženým modulem pružnosti E a modifikovanou Poissonovou konstantou, viz obr. 1 * E*/E 0. 2 d 0 f R s 0. 1 0. 2 0. 3 0. 5 (f-d)/f 1

TZ 5 Shell & Tube - dimenzování Tato napětí patří do kategorie III – nesmí překročit napětí 3 S Přechodové napětí (trubka v trubkovnici) Dosah přechodových napětí p s R

TZ 5 Shell & Tube - dimenzování Téměř všechny předchozí vztahy lze odvodit z diferenciální rovnice deformace nosníku na pružném podkladě (spojité zatížení)

TZ 5 Shell & Tube - dimenzování Ztráta stability vnějším přetlakem trubek 1. Plastický kolaps (napětí je třeba držet pod mezí kluzu) 2. Elastický kolaps U plovoucí hlavy může dojít i ke ztrátě osové stability (vzpěr) F L

TZ 5 Shell & Tube - dimenzování Teplotní roztažnost materiál (1/0 C) ocel 12 E-6 až 19 E-6 Al 23 E-6 Cu 14 E-6 mosaz 18 E-6 Si 2 E-6 cín 27 E-6

TZ 5 Shell & Tube - dimenzování Kompenzátory v plášti d R Maximální napětí odpovídající prodloužení jedné vlny Tuhost kompenzátoru

TZ 5 Shell & Tube - vibrace Vibrace se týkají především trubek výměníků a jsou způsobeny ØPeriodickým odtrháváním vírů v úplavu při příčném obtékání trubek (Karmánova vírová stezka, vortex shedding) ØTurbulentními fluktuacemi (energetické turbulentní víry) ØAkustickými kmity stlačitelného media v plášti (stojaté akustické vlnění) Problémy vznikají, když se frekvence těchto budicích sil přiblíží vlastní frekvenci kmitání trubek vetknutých do trubkovnic a s omezeným pohybem v místě průchodu přepážkami (rezonance). Důsledkem vibrací je hluk (až 150 db), který ovšem ještě sám o sobě nezpůsobuje poruchu (to je případ akustických vibrací). Nebezpečnější je vzájemný kontakt kmitajících trubek nebo trubek a pláště (otěr), stejně jako otěr materiálu v místě přepážek. Kmity trubek způsobují i uvolnění trubek z trubkovnice (a netěsnosti). Vibrace trubek jsou i příčinou únavových lomů.

TZ 5 Shell & Tube - vibrace Při výpočtu výměníku z hlediska vibrací se postupuje takto 1. Stanoví se nejmenší vlastní frekvence příčného kmitání trubek f 0 (s ohledem na geometrii, f 0 klesá s kvadrátem délky trubky, uplatní se i závislost na axiální síle, f 0 trubek natahovaných teplotní dilatací roste) 2. Určí se základní frekvence budících sil (odtrhávání vírů v úplavu, frekvence turbulentních vírů, frekvence akustických kmitů). S výjimkou akustických kmitů je tato frekvence přímo úměrná rychlosti proudění tekutiny v mezitrubkovém prostoru. 3. Posoudí se riziko vzniku rezonancí (shody vlastních frekvencí a frekvence budících sil) 4. Odhadne se amplituda kmitů trubek (tento problém zatím není uspokojivě řešen a problematika vibrací stále zůstává jedním z klíčových a otevřených problémů navrhování výměníků – spolu s problematikou foulingu).

TZ 5 Shell & Tube - vibrace První vlastní frekvence f příčného kmitání trubky délky L s momentem setrvačnosti průřezu J je dán řešením diferenciální rovnice kmitání nosníku: U trubek výměníku je třeba tento základní vztah modifikovat, protože trubka je na několika podpěrách (přepážky), např. L L Cn n=3 10 n=počet stejně dlouhých sekcí uprostřed n>5 1 Hmotnost trubky+tekutiny+ekvivalentní hmotnost tekutiny vně trubky vztažená na 1 m délky 0. 1 1 10

TZ 5 Shell & Tube - vibrace Předchozí vztahy se týkaly první (a nejnebezpečnější) vlastní frekvence trubky, která není axiálně zatížená. V reálu ale zatížená je, především teplotními dilatacemi. Pokud je teplota trubek vyšší než teplota pláště budou trubky trubkovnicí stlačovány (a obráceně, chladnější trubky jsou natahovány). Korekce základní vlastní frekvence na osovou sílu F (kladná při natahování): Při natahování trubek se tedy frekvence vlastních kmitů zvýší (každý, kdo hraje na kytaru, to ví). Tato změna vlastní frekvence je typicky 40% a rozhodně ji nelze zanedbat.

TZ 5 Shell & Tube - vibrace Téměř všechny předchozí vztahy lze odvodit z diferenciální rovnice kmitajícího nosníku

TZ 5 Shell & Tube - vibrace Při příčném obtékaní trubky (svazku trubek) dochází k periodickému odtrhávání vírů, která začíná už při hodnotách Reynoldsova čísla Re>100. Tyto víry vyvolávají síly, které jsou kolmé na směr proudění. Jejich frekvence je úměrná rychlosti nabíhajícího proudu u St je Strouhalovo číslo, které je režimu plně vyvinuté turbulence a pro osamocenou trubku přibližně konstantní St=0. 21. U svazku trubek je hodnota Strouhalova čísla (a tedy i frekvence) vyšší, HEDH ji vyjadřuje jako funkci rozteče trubek grafem St 0. 8 0. 6 0. 4 0. 2 1 1. 5 2 s/D

TZ 5 Shell & Tube - vibrace Turbulentní víry (turbulent buffeting) mají široké spektrum pulzací. Tyto víry interagují s oscilacemi trubek, což vše ještě víc komplikuje. Na základě experimentů s turbulentním obtékání svazku trubek plyny, doporučuje HEDH korelaci pro dominantní frekvenci energetických turbulentních vírů D u st sl

TZ 5 Shell & Tube - vibrace Frekvence pulzací stlačitelného media v rezonanční dutině Průměr pláště Rychlost zvuku (ve vzduchu cca 300 m/s)

TZ 5 CFD Ardon

TZ 5 CFD Petr Fiala: Tepelně hydraulický a pevnostní výpočet výměníku tepla Diplomová práce FS VUT Brno, 2010 Ukázka diplomové práce jejímž jediným cílem bylo provést výpočet shell&tube výměníku tepla použitím 1. Standardního programu HTRI Xchanger Suite 5 (metodika TEMA) 2. Konečněprvkovou analýzu systémem ANSYS (proudění v hlavách programem ANSYS CFX, a pevnostní analýza sestavy hlav, trubkovnic a hrdel programem ANSYS Workbench 11)

TZ 5 CFD Diplomová práce FS VUT Brno, 2010 Zadané parametry výměníku (viz předchozí stránka) – dva tahy v trubkách, jeden tah v plášti. V plášti je horká voda 1250 C (musí být održována pod tlakem aby se nevypařila), která ohřívá chladnou vodu v trubkách (vstupní teplota 330 C). 125 70. 8 33 62. 4

TZ 5 CFD Diplomová práce FS VUT Brno, 2010 Ve svazku je cca 3000 trubek. Práce se ani nesnažila modelovat jednotlivé trubky, jen hlavy výměníku (modelář Solid. Works). Modelovala se jen polovina, vzhledem k symetrii geometrie. CFD (Computer Fluid Dynamics) programem ANSYS CFX se pak modelovalo turbulentní proudění. Cílem bylo zmapovat teplotní pole v kapalině a stanovit součinitele přenosu tepla (vyšly až podezřele vysoké 7000 -9000 W/m 2 K)

TZ 5 CFD Diplomová práce FS VUT Brno, 2010 Vypočtené teploty pak byly spolu se zvoleným tlakem vody použity ve statickém výpočtu tlakové nádoby (konečněprvkový systém ANSYS) nebyl to celý výměník, jen hlavy a kus pláště. Bylo hodnoceno napětí v místě hrdel, pláště v blízkosti plovoucí hlavy a v plovoucí hlavě. Není mi jasné jaké elementy byly použity (pravděpodobně skořepinové) ani kolik jich bylo (z obrázků je patrné, že síť měla cca 70000 uzlů). Trubkovnice byla nahrazeno plnou deskou bez otvorů jen s modifikovaným modulem pružnosti. Opravdu nevím, zda se dá výsledkům věřit.