IAPWS-IF97 i własności pary wodnej w obiegach jądrowych

Streszczenie

Elektrownia jądrowa nie zamienia ciepła z reaktora w prąd bezpośrednio. W klasycznych blokach energetycznych robi to obieg wodno-parowy: woda odbiera ciepło, para wykonuje pracę w turbinie, skraplacz zamyka obieg, a pompy podnoszą ciśnienie kondensatu. Żeby taki obieg policzyć, trzeba znać własności wody i pary: entalpię, entropię, gęstość, objętość właściwą, ciepła właściwe, temperaturę nasycenia i stopień suchości.^1,2

IAPWS-IF97 jest przemysłowym standardem obliczania tych własności. Zamiast ręcznie odczytywać tabele parowe, program może wyznaczyć punkt termodynamiczny z równań podzielonych na regiony: ciecz sprężona, para przegrzana, gęsty płyn blisko punktu krytycznego, krzywa nasycenia i para wysokotemperaturowa.¹ W tym serwisie ten standard jest podstawą kalkulatora punktu pary wodnej i kalkulatora cyklu Rankine'a.^2,3

Rozszerzenie tematu

Dlaczego elektrownia potrzebuje tabel parowych

W uproszczonym opisie mówi się: reaktor grzeje wodę, para kręci turbiną, generator produkuje prąd. Ten opis jest prawdziwy, ale zbyt ubogi, by policzyć moc.

Do bilansu trzeba wiedzieć:

• ile energii niesie kilogram pary na wlocie turbiny,
• ile energii pozostaje po rozprężeniu,
• ile pracy trzeba włożyć w pompę,
• ile ciepła trzeba doprowadzić w wytwornicy pary albo kotle,
• ile ciepła trzeba odrzucić w skraplaczu,
• jak mokra jest para na końcu turbiny,
• jaki przepływ pary odpowiada danej mocy cieplnej.

Wszystkie te pytania prowadzą do entalpii h, entropii s, objętości właściwej v i punktu nasycenia. Bez spójnej funkcji własności wody nie da się dobrze policzyć nawet prostego Rankine'a.

Entalpia jako waluta obiegu

W obiegach parowych najczęściej liczy się różnice entalpii. Jeśli para w turbinie przechodzi od stanu 1 do stanu 2, idealna praca jednostkowa turbiny jest w pierwszym przybliżeniu różnicą:

w_t = h_1 - h_2

Jeśli woda po pompie i para po podgrzaniu mają stany 4 i 1, ciepło doprowadzone w wytwornicy pary wynosi:

q_in = h_1 - h_4

Sprawność idealizowanego obiegu można więc zapisać jako:

eta = (w_t - w_p) / q_in

To dlatego sama temperatura pary nie wystarcza. Dwa punkty o podobnej temperaturze mogą mieć różne ciśnienia, różną suchość i różną entalpię. Tabela pary jest mapą stanów, nie ozdobą do wykresu.

Entropia i rozprężanie turbiny

Entropia jest potrzebna, bo idealne rozprężanie w turbinie przyjmuje się jako izentropowe, czyli zachodzące przy stałej entropii. W praktyce turbina ma sprawność izentropową: rzeczywisty spadek entalpii jest mniejszy niż spadek idealny.

W kalkulatorze Rankine'a logika jest następująca:

1. z p_1 i T_1 liczymy h_1 oraz s_1,
2. dla ciśnienia skraplacza szukamy punktu izentropowego o s = s_1,
3. dostajemy h_2s,
4. sprawność turbiny przesuwa wynik do realnego h_2,
5. z h_2 i danych nasycenia oceniamy suchość pary na wylocie.³

To jest bardzo dobry przykład, po co potrzebne są funkcje odwrotne albo numeryczne szukanie punktu. Obieg energetyczny nie zawsze pyta "jaka jest entalpia dla p,T". Czasem pyta "jaka jest entalpia dla p,s".

Suchość pary i erozja turbiny

Para mokra jest mieszaniną pary nasyconej i kropelek cieczy. Jej stan opisuje stopień suchości x:

x = 0 oznacza ciecz nasyconą, a x = 1 parę nasyconą suchą.

Dla mieszaniny nasyconej wiele wielkości liczy się liniowo między gałęzią cieczy i pary:

h = h_f + x (h_g - h_f)

gdzie h_f jest entalpią cieczy nasyconej, a h_g entalpią pary nasyconej. Analogicznie można liczyć np. entropię mieszaniny.

W turbinie zbyt mokra para jest problemem mechanicznym. Krople uderzają w łopatki, zwiększają erozję i obniżają sprawność. Lokalny kalkulator ostrzega, gdy końcowa suchość pary spada poniżej progu dydaktycznego 0.85, bo wtedy prosty obieg Rankine'a przestaje być wiarygodnym obrazem realnego bloku i trzeba myśleć o separatorach wilgoci, przegrzewie wtórnym albo innej konfiguracji turbiny.^2,3

Czym jest IAPWS-IF97

IAPWS-IF97 to IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam. Oficjalny release IAPWS zaleca tę formulację do zastosowań przemysłowych i wskazuje, że zastąpiła wcześniejszą przemysłową formulację IFC-67.¹

Słowo "industrial" jest tu kluczowe. IAPWS ma także formulację naukową, zbudowaną pod bardzo szeroki i precyzyjny opis zwykłej wody. IF97 jest natomiast zaprojektowana tak, by była szybka i wygodna w obliczeniach technicznych: turbinach, kotłach, elektrowniach, wymiennikach ciepła i symulatorach procesu.

Oficjalny release dzieli zakres ważności na kilka regionów i dla każdego daje równania podstawowe oraz równania pomocnicze albo wsteczne. Dzięki temu obliczenia mogą być szybkie i stabilne w typowych zadaniach inżynierskich.

Zakres ważności

Release IAPWS-IF97 obejmuje dwa główne zakresy temperatury i ciśnienia:

• 273.15 K <= T <= 1073.15 K przy p <= 100 MPa,
• 1073.15 K < T <= 2273.15 K przy p <= 50 MPa.¹

Rewizja z 2007 roku dotyczyła przede wszystkim rozszerzenia regionu 5: dla temperatur 1073 K do 2273 K górny zakres ciśnienia został rozszerzony z 10 MPa do 50 MPa.¹

Te zakresy są znacznie szersze niż typowy prosty obieg jądrowy, ale ich znajomość jest ważna. Kalkulator powinien ostrzegać, kiedy użytkownik wychodzi poza zakres formulacji albo wchodzi w region, gdzie obliczenie jest numerycznie trudniejsze.

Pięć regionów IF97

IF97 dzieli płaszczyznę p-T na regiony:

• region 1: ciecz sprężona lub podchłodzona,
• region 2: para przegrzana,
• region 3: gęsty płyn w pobliżu punktu krytycznego,
• region 4: krzywa nasycenia ciecz-para,
• region 5: para wysokotemperaturowa.¹

To rozdzielenie nie jest kosmetyką. Woda zachowuje się inaczej jako prawie nieściśliwa ciecz, inaczej jako para przegrzana, inaczej w pobliżu punktu krytycznego i inaczej jako mieszanina nasycona. Jedno proste równanie na całym obszarze byłoby albo powolne, albo niedokładne, albo niewygodne.

Region 1: ciecz sprężona

Region 1 obejmuje ciecz w zakresie od 273.15 K do 623.15 K, przy ciśnieniach od ciśnienia nasycenia do 100 MPa.¹ W obiegu parowym odpowiada to wodzie zasilającej, kondensatowi po pompie i cieczy pod wysokim ciśnieniem.

W regionie 1 IF97 używa równania fundamentalnego dla właściwej energii swobodnej Gibbsa g(p,T) zapisanej bezwymiarowo jako gamma = g/(RT). Z pochodnych tej funkcji oblicza się objętość właściwą, entalpię, entropię, ciepła właściwe i inne wielkości.¹

Dla inżyniera najważniejsze jest to, że ciecz ma małą objętość właściwą. Dlatego praca pompy zwykle jest znacznie mniejsza niż praca turbiny i można ją dobrze przybliżyć przez:

w_p approx v_f (p_boiler - p_cond) / eta_p

To przybliżenie stosuje lokalny kalkulator Rankine'a.³

Region 2: para przegrzana

Region 2 opisuje parę przegrzaną i część pary metastabilnej. W typowym Rankine'ie stan na wlocie turbiny często leży właśnie tutaj albo blisko granicy nasycenia.

Tak jak region 1, region 2 jest opisany równaniem Gibbsa g(p,T), ale w postaci rozdzielonej na część gazu idealnego i część resztową.¹ To dobrze pasuje do pary: część zachowania przypomina gaz, a część wynika z rzeczywistych oddziaływań między cząsteczkami.

W kalkulatorze użytkownik widzi tylko h, s, v i fazę. Pod spodem są pochodne funkcji Gibbsa.

Region 3: okolice punktu krytycznego

Region 3 jest najtrudniejszy. Obejmuje gęsty płyn w pobliżu punktu krytycznego wody:

• T_c = 647.096 K,
• p_c = 22.064 MPa,
• rho_c = 322 kg/m3.¹

Tu własności zmieniają się gwałtownie, a rozdział ciecz/para traci prosty sens. IF97 używa w regionie 3 równania Helmholtza f(rho,T), a nie bezpośrednio g(p,T).¹

To ma konsekwencję programistyczną: jeśli użytkownik podaje p,T, kalkulator musi najpierw znaleźć gęstość rho, dla której równanie daje zadane ciśnienie. Lokalna klasa SteamTables robi właśnie to: numerycznie dobiera gęstość, a potem z pochodnych Helmholtza wyznacza h, u, s, v, c_p, c_v i prędkość dźwięku.²

Region 4: nasycenie

Region 4 nie jest obszarem jednofazowym. To krzywa nasycenia ciecz-para. IF97 daje równania na ciśnienie nasycenia p_s(T) i temperaturę nasycenia T_s(p). Zakres obejmuje linię od okolic punktu potrójnego do punktu krytycznego, a równanie temperatury nasycenia obejmuje zakres ciśnień od około 611 Pa do 22.064 MPa.¹

W praktyce region 4 odpowiada pytaniom:

• przy jakiej temperaturze wrze woda przy danym ciśnieniu,
• jakie są własności cieczy nasyconej,
• jakie są własności pary nasyconej,
• jaka jest entalpia parowania,
• jaka jest suchość mieszaniny.

To kluczowe dla skraplacza, końca turbiny i mokrej pary. Jeżeli po rozprężeniu turbiny punkt ląduje wewnątrz kopuły nasycenia, kalkulator musi przejść z języka p,T na język p,x albo T,x.

Region 5: para wysokotemperaturowa

Region 5 obejmuje wysokotemperaturową parę w zakresie 1073.15 K <= T <= 2273.15 K i p <= 50 MPa.¹ W klasycznych blokach jądrowych PWR/BWR zwykle nie jest to główny region pracy, bo temperatury pary są dużo niższe. Region 5 staje się jednak ważny dla szerszej energetyki cieplnej, wysokotemperaturowych koncepcji reaktorów (HTGR, MSR powyżej 800°C), walidacji kalkulatora i porównania z obiegami niejądrowymi — w tym z obiegiem He turbiny Braytona dla reaktorów gazowych chłodzonych.

Lokalny kalkulator punktu pary obsługuje region 5, ale z jasnym ostrzeżeniem zakresowym. To właściwe podejście: nie każdy region jest równie istotny dla PWR, ale standard powinien mieć spójny obszar działania.

Punkt nasycenia w kalkulatorze

Kalkulator punktu pary przyjmuje p, T i opcjonalną suchość x. Jeśli punkt leży blisko krzywej nasycenia, wykorzystuje region 4 i podaną suchość. W przeciwnym razie wybiera region 1, 2, 3 albo 5 na podstawie temperatury, ciśnienia i granicy nasycenia.²

To rozwiązuje częsty błąd dydaktyczny. Użytkownik wpisuje np. temperaturę bliską wrzenia przy danym ciśnieniu i oczekuje "pary", ale bez informacji o suchości stan nie jest jednoznaczny. Na krzywej nasycenia ta sama para wodna może być mieszaniną o bardzo różnych entalpiach zależnie od x.

Dlatego formularz pyta o suchość, gdy punkt jest nasycony. Jest to fizyczny wymóg jednoznacznego opisu stanu termodynamicznego, nie niedopatrzenie interfejsu użytkownika. Na krzywej nasycenia stan układu opisują dwie niezależne zmienne: ciśnienie (lub temperatura) i suchość x — temperatura nasycenia jest wtedy funkcją ciśnienia i nie dostarcza dodatkowej informacji.

Cykl Rankine'a w serwisie

Kalkulator Rankine'a używa IF97 do prostego obiegu:

1. wlot turbiny: p_1, T_1,
2. wylot turbiny do skraplacza,
3. ciecz nasycona po skraplaczu,
4. ciecz po pompie zasilającej.³

Następnie liczy:

• entalpię wlotową h_1,
• entropię wlotową s_1,
• izentropowy wylot turbiny h_2s,
• rzeczywisty wylot turbiny h_2,
• pracę turbiny,
• pracę pompy,
• ciepło doprowadzone,
• sprawność termiczną,
• moc elektryczną,
• przepływ pary,
• strumień wody chłodzącej dla przyjętego Delta T.

To nie jest projekt turbiny. To jest dydaktyczny bilans energii, który pokazuje, skąd biorą się liczby rzędu 30-37% dla wielu bloków jądrowych.

Dlaczego sprawność jądrowego Rankine'a bywa niższa

Typowy PWR pracuje z temperaturami pary dużo niższymi niż nowoczesny blok gazowo-parowy. Ograniczenia pochodzą z temperatur w obiegu pierwotnym, materiałów, marginesów wrzenia, wytwornic pary, bezpieczeństwa i konstrukcji turbiny na parę nasyconą albo słabo przegrzaną.

Z punktu widzenia termodynamiki oznacza to niższą górną temperaturę obiegu. Granica Carnota:

eta_C = 1 - T_cold / T_hot

nie jest osiągalną sprawnością turbiny, ale pokazuje kierunek. Im niższa temperatura źródła ciepła, tym niższy sufit sprawności. Kalkulator Rankine'a zestawia sprawność termiczną z granicą Carnota po to, by użytkownik widział, że strata nie wynika z "nieudolności elektrowni", tylko z temperatur i nieodwracalności rzeczywistego obiegu.³

PWR, BWR, CANDU i presety

W formularzu Rankine'a są presety odpowiadające orientacyjnym punktom pracy: typowy PWR, BWR, CANDU oraz wariant wysokotemperaturowy jako porównanie.³

Ich rola jest dydaktyczna. Nie są kartą projektową konkretnej elektrowni. Mają pokazać:

• niższe temperatury pary w klasycznej energetyce jądrowej,
• zależność sprawności od skraplacza,
• wpływ sprawności izentropowej turbiny,
• ryzyko zbyt mokrej pary na końcu rozprężania,
• różnicę między mocą cieplną reaktora MWt i mocą elektryczną MWe.

Szczególnie ważne jest rozróżnienie MWt i MWe. Reaktor o mocy 3200 MWt przy sprawności 33% nie daje 3200 MWe, tylko około jednej trzeciej tej wartości.

Skraplacz i ciepło odrzucone

Skraplacz jest termodynamicznie równie ważny jak reaktor. Obniża ciśnienie końca turbiny, zwiększa spadek entalpii i umożliwia zamknięcie obiegu jako ciecz. Ale całe ciepło, którego nie zamieniono na pracę, trzeba odprowadzić do chłodzenia.

W prostym kalkulatorze ciśnienie skraplacza wynika z temperatury nasycenia. Jeśli użytkownik obniży temperaturę skraplacza, poprawia sprawność, ale jednocześnie zakłada lepsze warunki chłodzenia. To dobry sposób na zrozumienie, dlaczego lokalizacja elektrowni, dostęp do wody chłodzącej i warunki pogodowe mają znaczenie energetyczne.

Regeneracja i przegrzew wtórny

Rzeczywisty blok nie jest prostym prostokątem na wykresie T-s. Ma podgrzewacze regeneracyjne, upusty pary, separatory wilgoci, przegrzew wtórny i wielosekcyjną turbinę.

Lokalny kalkulator reheat/regeneration jest jakościowym pomostem. Pokazuje, że:

• regeneracja podnosi temperaturę wody zasilającej i poprawia sprawność,
• przegrzew wtórny ogranicza wilgoć na końcu turbiny,
• upusty pary zabierają część przepływu z głównej ekspansji,
• sprawność netto jest niższa niż sprawność brutto.⁴

Ten kalkulator celowo nie jest pełnym modelem turbiny. Ma pokazać, dlaczego rzeczywisty schemat jest bardziej złożony niż podstawowy Rankine.

Co IF97 daje programiście

Dla programisty kalkulatora IF97 daje trzy korzyści.

Po pierwsze, daje spójne regiony. Kod nie musi mieszać wielu przypadkowych tabel i interpolacji.

Po drugie, daje równania pochodnych. Z funkcji Gibbsa albo Helmholtza można policzyć wiele własności w sposób termodynamicznie spójny.

Po trzecie, daje testy weryfikacyjne. Release zawiera wartości kontrolne dla wybranych punktów, dzięki którym implementację można porównać z dokumentem źródłowym.¹

To ostatnie jest bardzo ważne. Błąd w entalpii rzędu kilku procent może przesunąć sprawność, przepływ pary i wnioski o wilgotności. Kalkulator edukacyjny nie musi być narzędziem projektowym, ale nie powinien liczyć z przypadkowych wzorów znalezionych w internecie — zwłaszcza że IF97 jest standardem wolnodostępnym, dobrze udokumentowanym i weryfikowalnym przez każdego. Kalkulator edukacyjny nie musi być narzędziem projektowym, ale nie powinien liczyć z przypadkowych wzorów znalezionych w internecie.

Co IF97 nie robi

IF97 liczy własności czystej wody i pary. Nie liczy automatycznie:

• strat ciśnienia w rurach,
• sprawności poszczególnych stopni turbiny,
• separacji wilgoci,
• wymienników ciepła,
• sprzężenia z reaktorem,
• chemii obiegu,
• korozji,
• drgań łopatek,
• pełnego bilansu upustów,
• odpowiedzi dynamicznej przy zmianie mocy.

To jest standard własności termodynamicznych, nie symulator elektrowni. Właśnie dlatego artykuł należy do kategorii "Dane i modele": IF97 jest warstwą danych fizycznych, na której dopiero buduje się model procesu.

Niepewności i punkt krytyczny

Oficjalny release zawiera oszacowania niepewności dla wybranych własności, m.in. objętości właściwej, ciepła właściwego, prędkości dźwięku i ciśnienia nasycenia. Dokument zwraca uwagę, że niepewności rosną w pobliżu punktu krytycznego, zwłaszcza dla takich wielkości jak c_p i prędkość dźwięku.¹

To jest ważna lekcja ogólna. Standard przemysłowy nie znaczy "magicznie dokładny wszędzie tak samo". Region blisko krytyczny jest fizycznie trudny, a nie tylko programistycznie niewygodny.

Jak opisywać wyniki kalkulatora

Dobry wynik punktu pary powinien pokazywać:

• region IF97,
• fazę,
• T,
• p,
• h,
• s,
• v albo gęstość,
• u,
• c_p,
• c_v,
• prędkość dźwięku, jeśli jest liczona,
• temperaturę nasycenia dla porównania,
• ostrzeżenia zakresowe.

Dobry wynik Rankine'a powinien dodatkowo pokazywać punkty obiegu i wyraźnie rozdzielać:

• sprawność Carnota,
• sprawność termiczną obiegu,
• moc cieplną,
• moc elektryczną,
• ciepło odrzucone,
• przepływ pary,
• suchość końcową.

Bez tego użytkownik widzi jedną liczbę sprawności, ale nie rozumie, skąd się wzięła.

Ciepło resztkowe po wyłączeniu — dlaczego reaktor nie przestaje grzać

Gdy reaktor zostaje wyłączony (SCRAM), reakcja łańcuchowa natychmiast ustaje, ale ciepło resztkowe (decay heat) nadal się wytwarza przez rozpad produktów rozszczepienia. Wartość ciepła resztkowego zaraz po wyłączeniu wynosi ok. 7% mocy nominalnej i maleje według przybliżonej formuły (norma ANS-5.1):

$$P_{decay}(t) \approx 0,066 \cdot P_0 \cdot \left[t^{-0,2} - (t + t_0)^{-0,2}\right]$$

Ciepło resztkowe musi być odprowadzane do skraplacza lub do systemów awaryjnego chłodzenia rdzenia (ECCS). Tutaj IAPWS-IF97 pojawia się w analizach bezpieczeństwa:

• Jaka jest temperatura pary w primary coolant loop po SCRAM?
• Czy ciśnienie w zbiorniku reaktora (PWR) pozostaje poniżej limitu bezpieczeństwa?
• Czy naturalna cyrkulacja wody reaktora wystarczy do odprowadzenia ciepła bez pomp?

Dla BWR dodatkowy problem: para w reaktorze może kondensować po wyłączeniu, co wpływa na ciśnienie i geometrię cyrkulacji. Obliczenia transientowe dla SCRAM wymagają IF97 jako podrutyny termodynamicznej w kodach bezpieczeństwa (RELAP5, TRACE, MELCOR).

Fukushima a ciepło resztkowe: gdy po tsunami zawiodło zasilanie pomp chłodzenia, ciepło resztkowe (7→3→1% mocy w trakcie kilku godzin) nie miało gdzie uciec. Bez chłodzenia temperatura paliwa UO₂ rosła powyżej granicy integralności pokrycia cyrkonowego — stąd wyzwolenie wodoru i wybuchy. Ciepło resztkowe, obliczane z ORIGEN i zarządzane przez RELAP5/TRACE, jest jednym z kluczowych parametrów bezpieczeństwa każdego reaktora.^1,5

Turbiny na parę mokrą vs przegrzaną — szczegóły projektowe

W klasycznych reaktorach parowych (BWR, PWR) turbina pracuje ze słabo przegrzaną lub mokrą parą, co różni się od turbin elektrowni konwencjonalnych:

Para mokra:

• Krople cieczy (x < 1) uderzają w łopatki — erozja (szczególnie ostatnie stopnie)
• Typowy dopuszczalny limit: x_end ≥ 0,85 na wyjściu z turbiny
• Separatory mokrości między częścią WP (wysokociśnieniową) a NP (niskociśnieniową)

Para lekko przegrzana (BWR/PWR):

• Temperatura odparowania nasycenia + kilka °C przegrzewu
• Redukuje erozję łopatek, ale nie eliminuje całkowicie

Turbiny wielowalcowe:
Duże reaktory mają turbiny z kilkoma walcami (HP turbine + 2-4× LP turbines) i separatorami mokrości między nimi. Para z wylotu HP turbiny przechodzi przez separator wilgoci i ewentualny wtórny podgrzewacz (moisture separator reheater, MSR) przed wlotem do LP turbiny.

Obliczanie tego cyklu w IF97:

1. Punkt wylotowy HP turbiny: z p_HP_out, h_HP_out (znając s i sprawność HP)
2. Po separatorze: tylko sucha para (x = 1 dla par, ciekła odnoga do spływu)
3. Para po MSR: wyższa entalpia dzięki dostarczeniu ciepła
4. LP turbina: dalsze rozprężanie do ciśnienia skraplacza

Dobre modelowanie tego procesu wymaga dwóch dodatkowych par (p,h) — jedno dla wylotu HP i jedno dla wlotu LP po MSR. Lokalny kalkulator reheat/regeneration pokazuje uproszczone podejście.^3,4

Prędkość dźwięku w parze i kryterium choking

Prędkość dźwięku jest własością termodynamiczną obliczalną z IF97. Dla pary w typowych warunkach PWR (280°C, 60 bar): ok. 500–550 m/s.

Dlaczego to ważne w inżynierii:
Zjawisko "choking" (dławienie) pojawia się, gdy prędkość pary w pewnym przekroju osiąga prędkość dźwięku. Dalsze zwiększenie przepływu jest niemożliwe bez zmiany geometrii. Dla turbiny parowej i zaworów bezpieczeństwa:

• Maksymalny przepływ przez dyszę turbiny lub zawór bezpieczeństwa jest ograniczony warunkiem krytycznego przepływu
• Obliczenie krytycznego przepływu wymaga prędkości dźwięku i gęstości pary z IF97

Przypadek Fukushimy:
Zawory bezpieczeństwa i zawory relief reaktora (SRV/SV) działają przez krytyczny przepływ pary. Obliczenie przepustowości tych zaworów podczas transientu wymaga IF97 dla konkretnych ciśnienia i entalpii pary w zbiorniku reaktora. Kody RELAP5 używają IF97 (lub IAPWS-95) do tego obliczenia przy każdym kroku czasu symulacji.^1,5

Związek z energetyką jądrową

IF97 nie jest "jądrowe" w sensie fizyki jądrowej. Jest jednak niezbędne dla energetyki jądrowej, bo większość elektrowni jądrowych jest od strony turbiny klasyczną elektrownią parową.

Reaktor dostarcza ciepło. Cała reszta jest termodynamiką obiegu:

• PWR oddziela obieg pierwotny i wtórny przez wytwornicę pary,
• BWR wytwarza parę bezpośrednio w zbiorniku reaktora,
• CANDU ma własną konfigurację kanałów paliwowych, chłodziwa i wytwornic,
• przyszłe reaktory wysokotemperaturowe mogą przesuwać punkt pracy ku wyższym temperaturom i innym obiegom.

Dlatego student atomistyki powinien umieć przejść od jądra reaktora do turbiny bez utraty bilansu energii. IAPWS-IF97 jest jednym z narzędzi tego przejścia.

Konkretne parametry pracy dla głównych typów reaktorów

Dla studenta obliczającego sprawność elektrowni jądrowej, kluczowe są rzeczywiste parametry obiegów. Dane orientacyjne (nie gwarantowane dla konkretnych instalacji):

PWR (pressurized water reactor — reaktor ciśnieniowy):

• Obieg pierwotny: ~155 bar, 285–320°C
• Para na wlocie turbiny (po wytwornicy): ok. 55–75 bar, ok. 275–285°C (nasycona lub lekko przegrzana)
• Temperatura skraplacza: ~30–35°C (rzeka/morze) lub ~40°C (wieże chłodnicze)
• Sprawność brutto: ok. 33–35%
• Przykłady: EPR (1650 MWe, 450 bar/323°C w obiegu pierwotnym, para ~72 bar), AP1000

BWR (boiling water reactor — reaktor wrzątkowy):

• Para wytwarzana bezpośrednio w reaktorze: ~70–75 bar, ok. 285°C
• Brak wytwornicy pary między reaktorem a turbiną (para z reaktora bezpośrednio do turbiny)
• Turbina musi być zaprojektowana na pracę z parą radioaktywną
• Sprawność brutto: ok. 33–35%
• Przykłady: ABWR, ESBWR

CANDU (Canadian Deuterium Uranium):

• Chłodziwo pierwotne (ciężka woda): ~100 bar, ok. 310°C
• Para w układzie wtórnym: ok. 45–50 bar, ok. 260°C
• Sprawność brutto: ok. 29–32% (niższa ze względu na niższe temperatury)
• Zaleta: możliwość wymiany paliwa bez zatrzymania

HTGR / HTR (wysokotemperaturowe reaktory gazowe):

• Chłodziwo: hel pod ciśnieniem ok. 70 bar, temperatura wlotu/wylotu: ~400/850–950°C
• Możliwy obieg Braytona (turbina gazowa z He) — sprawność >40–50% teoretycznie
• Obieg parowy z wysoką temperaturą pary: możliwość np. 160–200 bar, 510–540°C
• Pebble bed HTR (Chiny, projekt HTR-PM): 250 MWth, sprawność >40%

Zestawienie tych danych pokazuje bezpośrednio, dlaczego klasyczne PWR/BWR mają niższą sprawność niż nowoczesne elektrownie gazowo-parowe: ograniczenie wynika z temperatur chłodziwa i materiałów wytwornic pary, nie z "niedoskonałości" Rankine'a.^1,3

Historia standardów własności wody — od tablic do IF97

Zanim pojawiły się cyfrowe standardy, inżynierowie używali drukowanych tablic pary wodnej. Historia standardyzacji:

• 1934: Pierwsze międzynarodowe tablice pary (Steam Tables 1934) — ręcznie interpolowane, 6 cyfr znaczących
• 1967: IFC-67 (International Formulation Committee) — pierwszy cyfrowy standard dla przemysłu parowego, zastąpił różne krajowe tablice
• 1997: IAPWS-IF97 zastępuje IFC-67 — bardziej spójna matematycznie, szybsza obliczeniowo, lepsza dokładność w pobliżu punktu krytycznego
• 2007: Rewizja IF97 dla regionu 5 — rozszerzony zakres ciśnienia do 50 MPa
• 2012: Aktualna wersja (IAPWS R7-97(2012)) — z drobnymi korektami numerycznymi

Równolegle: IAPWS-95 (standard naukowy):

• Formulacja z 1995 roku dla badań naukowych — jedna funkcja Helmholtza na cały zakres
• Dokładniejsza od IF97, ale wolniejsza obliczeniowo
• Zakres: 251,165 K – 1273 K, do 1 GPa
• IF97 był tworzony jako szybki i dokładny odpowiednik IF95 dla zastosowań przemysłowych

Dla programisty: jeśli kalkulator edukacyjny potrzebuje sprawdzić wynik IF97, IAPWS-95 jest odpowiednikiem "ground truth" (jest wolniejszy, ale bardziej generalny).¹

Superkrytyczna woda — reaktory IV generacji SCWR

Jednym z reaktorów IV generacji zaproponowanych przez Generation IV International Forum (GIF) jest SCWR (Supercritical Water Reactor — reaktor z superkrytyczną wodą). Superkrytyczna woda to woda w stanie powyżej punktu krytycznego:

• T > 647 K (374°C)
• p > 22,064 MPa

W superkrytycznym stanie woda nie wrze — przejście ciecz-gaz znika. Para turbinowa może mieć parametry 25 MPa / 625°C, co daje sprawności >44%.

Termodynamika superkrytyczna a IF97:
Właśnie dla tego zakresu region 3 IAPWS-IF97 jest kluczowy. Superkrytyczna woda leży powyżej punktu krytycznego — brak wyraźnej granicy faz, ale ciągły gradient własności od "gęstej" do "rzadkiej" wody.

Wyzwania materiałowe:
W 625°C z wodą pod 25 MPa żaden znany materiał reaktorowy nie jest w pełni sprawdzony. Stale ferrytyczno-martenzytyczne (HT9, F82H, EUROFER97) są kandydatami, ale wymagają badań odporności na promieniowanie i korozję przez kilkadziesiąt lat.

Polska a SCWR:
Polska nie ma aktywnego programu SCWR, ale NCBJ uczestniczy w europejskim projekcie badawczym CORTEX i innych projektach EU H2020/Horizon dotyczących zaawansowanych reaktorów. Wiedza o termodynamice superkrytycznej jest przydatna dla doktorantów orientujących się w trendach reaktorowych.^1,5

Ciepło skojarzone (CHP) — elektrownia jądrowa jako producent ciepła

Elektrownia jądrowa nie musi produkować tylko prądu. Ciepło odrzucone do skraplacza (stanowiące 60–70% energii reaktora) można zamiast marnować — dostarczać do sieci ciepłowniczych.

Nuclear CHP (Combined Heat and Power — trigeneracja):

• Para z turbiny upuszczona przed skraplaczem i dostarczona do wytwornicy ciepłowniczej
• Woda ciepłownicza (90–120°C) ogrzewa domy, budynki przemysłowe
• Elektrownia traci część mocy elektrycznej, ale zysk w cieple kompensuje straty

Przykłady:

• Loviisa (Finlandia): ciepło dla Helsinek pochodzi ze skraplacza reaktora
• Szwajcaria: dyskusja o ciepłownictwie z Leibstadt NPP
• Chiny: projekt ACP100 (SMR) zaprojektowany jako reaktor ciepłowniczy dla miast

Termodynamika CHP z IF97:
Dla obliczenia CHP kalkulator musi umieć modelować "upust pary" — wylot pośredni turbiny przekierowany do wytwornicy ciepłowniczej. Entalpię i temperaturę takiego upustu oblicza się z IF97 dla odpowiedniego ciśnienia i entropii. Następnie bilans energii rozdziela moc elektryczną i cieplną.^1,5

Woda zasilająca — jakość i chemia obiegu

Własności termodynamiczne wody obliczone przez IF97 zakładają czystą wodę. W rzeczywistości obieg wodno-parowy elektrowni wymaga ścisłej kontroli chemii:

Parametry kontrolowane w wodzie zasilającej:

• pH: 9,0–9,5 (alkaline dla ochrony stali w obiegu wtórnym, 7,5–8,0 dla Inconel wytwornicy)
• Tlen roztworzony: <10 ppb (tlen przyspiesza korozję)
• Żelazo (Fe): <2 ppb
• Miedź (Cu): <2 ppb (szczególnie problem w systemach z mosiądzem)
• Utlenienie: AVT (All-Volatile Treatment) — amoniak/hydrazyna
• Twardość: woda demineralizowana (odsolona)

Dlaczego to nie przeszkadza IF97:
Małe stężenia domieszek (ppb) nie zmieniają własności termodynamicznych wody mierzalnie. IF97 dla czystej wody jest w praktyce dokładne dla wody technicznej elektrowni. Korozja i osady (szczególnie magnetyt Fe₃O₄) wpływają na wymianę ciepła w wytwornicach — to osobny problem (fouling), który IF97 nie modeluje.

Woda zasilająca a reaktory PWR:
W PWR objętość wody w układzie pierwotnym (chłodziwo reaktora) jest odseparowana od obiegu wtórnego. Do pierwotnego dodaje się boranu litu (LiOH + kwas borowy H₃BO₃) jako moderator reaktywności. IF97 nie jest używany dla wody z borem (właściwości termodynamiczne wody borowej różnią się nieznacznie od czystej wody).^1,2

Wykres T-s i h-s dla typowych obiegów jądrowych

Diagram T-s (temperatura-entropia):
Dla prostego obiegu Rankine'a (bez przegrzewu) punkty tworzą charakterystyczną "pętle":

• 1→2: izentropowe rozprężanie w turbinie (pionowa linia w dół)
• 2→3: izobarne skraplanie (pozioma linia na dole, pod kopułą nasycenia)
• 3→4: izentropowe pompowanie (pionowa linia w górę, niewidoczna bo mała)
• 4→1: izobaryczne nagrzewanie i parowanie w wytwornicy (pozioma linia na górze, częściowo pod kopułą)

Stosunek pola powierzchni "wnętrza pętli" do pola pod krzywą górną wizualizuje sprawność termodynamiczną.

Diagram h-s (entalp-entropia, wykres Molliera):
Turbinowe obliczenia są często wykonywane na diagramie Molliera:

• Izentropowe rozprężanie: pionowa linia w dół (h maleje, s stałe)
• Rzeczywiste rozprężanie: linia w prawo i w dół (h i s rosną przez nieodwracalność)
• Kopuła nasycenia: wyraźna granica mieszaniny mokrej od pary suchej

Diagram Molliera jest standardem w inżynierii turbin parowych — narzędzie codziennego projektowania na poziomie każdego etapu turbiny.^1,3

SMR — małe reaktory modułowe a termodynamika pary

Małe reaktory modułowe (SMR, < 300 MWe) są rozwijane przez wiele firm i krajów. Większość z nich to reaktory wodne, ale o innych parametrach niż duże PWR:

NuScale (USA):

• 77 MWe per moduł (do 12 modułów = 924 MWe na jednej lokalizacji)
• Pasywne chłodzenie bez pomp (ciśnienie pierwotne ~12 MPa, temperatury niższe niż typowy PWR)
• Para do turbiny: ~35–40 bar (możliwa sprawność ~28–30%)

BWRX-300 (GE-Hitachi):

• 300 MWe, wrzątkowy, bezpośredni parowanie w reaktorze
• Para ~74 bar, ok. 285°C — parametry typowe BWR

RITM-200 (AAEM Rosja):

• 55 MWe per moduł (lodołamacz Arktika, dwie jednostki)
• Reaktor pierwotny w pełni integruje wytwornicę pary — kompaktowa obudowa

Polskie plany:

• ORLEN Synthos Green Energy (OSGE) planuje zakup reaktorów BWRX-300
• PKN ORLEN i inne firmy badają SMR jako lokalne źródła ciepła i prądu
• NuScale jest analizowany przez NCBJ/PEJ jako opcja dla Kopalni Bełchatów

Wszystkie te reaktory wymagają tych samych obliczeń termodynamicznych z IF97 — parametry różne, ale równania IF97 te same. Dlatego kalkulator Rankine'a z możliwością wprowadzenia dowolnych parametrów jest przydatnym narzędziem do porównywania koncepcji SMR.⁵

Numeryczne aspekty implementacji IF97 — pułapki programistyczne

Implementując IF97, programista może napotkać kilka typowych problemów:

Granice regionów a numeryczna niestabilność:

• Na granicy regionu 1/2 (krzywa nasycenia) oba regiony dają różne wyniki. Lokalna implementacja rozwiązuje to przez wybór regionu po sprawdzeniu granicy nasycenia, nie po samym p,T.
• Region 3 wymaga iteracyjnego szukania ρ przy zadanym p,T — metodą Newtona-Raphsona lub bisekcji. Typowe kryterium zbieżności: |p_obliczone - p_zadane| < 1 Pa.

Funkcje wsteczne:
IF97 zawiera "backward equations" — równania inverse, które z p,h lub p,s wyznaczają T bezpośrednio, bez iteracji. To przyspiesza obliczenia w trybie "podaj h, dostań T" — kluczowym dla Rankine'a i programowania turbin.

Testy weryfikacyjne:
Release IF97 podaje wartości kontrolne dla specyficznych punktów. Prawidłowa implementacja powinna odtworzyć je z dokładnością do ostatniej podanej cyfry (zwykle 6–9 cyfr). Implementacja używana przez lokalny kalkulator punktu pary przeszła weryfikację dla kluczowych regionów.^1,2

Zakres temperatur Celsjusza kontra Kelwiny:
IF97 wymaga temperatury w Kelwinach. Konwersja T[K] = T[°C] + 273,15 musi być stosowana konsekwentnie — błąd o 0,15 K (przez użycie 273 zamiast 273,15) może być znaczący w okolicach punktu potrójnego wody.

Porównanie pracy turbiny parowej i gazowej — wniosek o sprawności

Warto zestawić obieg Rankine'a (para wodna) z obiegiem Braytona (gaz), by zrozumieć, dlaczego sprawności się różnią:

Parametr	Rankine (para)	Brayton (gaz)
Czynnik roboczy	H₂O (fazy)	powietrze/gaz
Temperatura max.	~540°C (para przegrzana)	~1300–1400°C (spaliny)
Temperatura min.	~30–45°C (skraplacz)	~400–500°C (exhaust)
Sprawność Carnota	~62% (przy 540/35°C)	~75% (przy 1350/450°C)
Realna sprawność	~30–38%	~38–42%

Obieg gazowo-parowy (CCGT — combined cycle gas turbine) łączy oba: turbina gazowa odrzuca ciepło do kotła, który generuje parę dla turbiny parowej. Sprawność CCGT sięga 60–62%, co jest niemalże nieosiągalne dla czystego reaktora parowego.

To pokazuje fundamentalną termodynamiczną "cenę" bezpieczeństwa jądrowego: reaktory LWR muszą utrzymać temperatury rdzenia poniżej granic bezpieczeństwa materiałów paliwowych i pokrycia cyrkonowego (~1200°C dla zr-alloy), co ogranicza temperaturę chłodziwa i pary znacznie poniżej temperatury spalin turbiny gazowej.

Reaktory wysokotemperaturowe (HTR, MSR z solą fluorową) mogą zbliżyć się do sprawności CCGT dzięki wyższym temperaturom (600–950°C), ale kosztem wyzwań materiałowych i technologicznych, które wciąż są w fazie R&D.^1,5

Podsumowanie

IAPWS-IF97 jest niewidzialnym fundamentem obliczeń wodno-parowych. Użytkownik widzi w kalkulatorze entalpię, entropię, suchość i sprawność, ale pod spodem działa podział na regiony, równania Gibbsa, równanie Helmholtza i krzywa nasycenia.

Najważniejsza lekcja jest prosta: elektrownia jądrowa jest jednocześnie urządzeniem jądrowym i cieplnym. Bez danych pary wodnej nie policzymy, co realnie dzieje się z ciepłem reaktora po stronie turbiny. IF97 jest mostem między fizyką reaktora a inżynierią turbiny — i zrozumienie, gdzie ten most stoi, jest podstawą dla każdego doktoranta energetyki jądrowej. Elektrownia jądrowa jest jednocześnie reaktorem i elektrownią parową — a dwa standardy danych (ENDF dla fizyki jądrowej i IAPWS-IF97 dla termodynamiki wody) razem opisują pełny łańcuch: od neutronu w paliwie do kilowatoGodziny w sieci elektroenergetycznej. Zrozumienie obu, ich zakresu ważności i ograniczeń, jest wyróżnikiem inżyniera/fizykia zdolnego myśleć o całej instalacji, a nie tylko o wybranej poddziedzinie.

Dodatkowe materiały multimedialne

Ćwiczenia praktyczne

Pierwsze ćwiczenie: dla zadanego punktu p,T określić, czy znajduje się w cieczy sprężonej, parze przegrzanej, regionie nasycenia czy pobliżu punktu krytycznego. Następnie porównać odpowiedź z kalkulatorem punktu pary.

Drugie ćwiczenie: wyjaśnić, dlaczego przy punkcie nasycenia sama para p,T nie wystarcza do określenia entalpii. Trzeba wprowadzić suchość x.

Trzecie ćwiczenie: policzyć ręcznie prosty bilans Rankine'a z podanych wartości h_1, h_2, h_4 i porównać znaczenie pracy turbiny oraz pracy pompy.

Czwarte ćwiczenie: zmienić temperaturę skraplacza w kalkulatorze Rankine'a i opisać, co dzieje się ze sprawnością, ciepłem odrzuconym i przepływem wody chłodzącej.

Piąte ćwiczenie: napisać krótką notę ostrzegawczą dla wyniku, w którym suchość końcowa pary jest zbyt niska. Nota ma wyjaśniać, dlaczego realna turbina potrzebuje separatora wilgoci albo przegrzewu wtórnego.

Szóste ćwiczenie: porównać rolę IF97 i ENDF: oba są standardami danych fizycznych, ale pierwszy dotyczy termodynamiki wody i pary, a drugi danych jądrowych (przekroje czynne, parametry rezonansów, wyniki rozszczepienia). Wskazać podobieństwa w podejściu (regionalne modele, walidacja punktami referencyjnymi, aktualizacje historyczne) i różnice (dziedzina: termodynamika vs. fizyka jądrowa; użytkownicy: inżynierowie obiegu vs. fizycy rdzenia).

Przejdź do ćwiczenia interaktywnego

Powiązane kalkulatory i narzędzia

• Punkt pary wodnej
• Cykl Rankine'a
• Reheat/regeneration cycle
• Bilans mocy bloku
• LCOE cashflow

Powiązane artykuły

• Dane jądrowe ENDF, GNDS i droga od przekrojów czynnych do kalkulatora
• Reaktory energetyczne: typy, moderator, chłodziwo, obieg i logika klasyfikacji
• Infrastruktura przed pierwszą elektrownią jądrową: NEPIO, dozór, kadry, KSE, paliwo, odpady i komunikacja społeczna
• Lokalizacja elektrowni jądrowej: sejsmika, powódź, chłodzenie, zaludnienie, KSE i dlaczego samo dużo miejsca nie wystarcza