Obliczenia_wymiennikow.pdf

(323 KB) Pobierz
Obliczenia płytowych wymienników ciepła
Materiał opublikowano z pewnymi
skrótami w czasopiśmie „Instalacje
Sanitarne” Nr 2(7)/2003, str. 27-30
Projektowanie płytowych wymienników ciepła
Lidia Zander, Zygmunt Zander
Uniwersytet warmińsko-Mazurski w Olsztynie
Katedra Inżynierii i Aparatury Procesowej
Wśród użytkowników płytowych wymienników ciepła panuje przekonanie, że aparaty tego
typu mogą być konfigurowane w sposób niemalże nieograniczenie zróżnicowany. Trzeba
przyznać, że w porównaniu z wymiennikami rurowymi istnieje o wiele więcej możliwości
dokonywania przebudowy istniejących urządzeń i dostosowywania ich do aktualnych potrzeb,
jeżeli tylko dysponuje się odpowiednio dużym zestawem płyt. W rzeczywistości jednak nie
można płyt tych zestawiać w sposób dowolny. Nie wystarcza tu przybliżone określenie wielkości
powierzchni wymiany ciepła i odpowiedniej liczby płyt. Taką samą powierzchnię wymiany
ciepła można bowiem uzyskać w wielu kombinacjach przepływów, co w konsekwencji daje
aparaty istotnie różne z punktu widzenia ich wydajności cieplnej i oporów przepływu.
Charakterystyka przepływu strumieni czynników w wymiennikach płytowych
Powierzchnię ogrzewalną płytowych wymienników ciepła zestawia się z tzw. płyt
strumieniowych. Ich powierzchnia jest pofalowana charakterystycznymi wytłoczeniami
zapewniającymi płycie odpowiednią sztywność i, co najważniejsze, pozwalającymi na
rozwijanie odpowiednio wysokiej burzliwości strumienia cieczy przy relatywnie niskich
prędkościach przepływu. Płyta zaopatrzona jest w cztery otwory rozmieszczone w narożnikach,
przy czym dla każdego czynnika płynącego przez aparat zarezerwowana jest jedna, zawsze ta
sama para otworów. Sposób uformowania uszczelki pozwala na to, że po powierzchni płyty od
strony uszczelki może przepływać tylko jeden z czynników – zimny, bądź gorący. Projektant
dysponuje zatem płytami tzw. lewymi i prawymi (rys. 1), które zestawia się na przemian. W ten
sposób otrzymuje naprzemiennie biegnące kanały nośnika ciepła i czynnika podgrzewanego.
Nakładające się otwory w narożnikach płyt tworzą cztery kolektory, za pośrednictwem których
oba czynniki są doprowadzane iodprowadzane z przestrzeni międzypłytowych. Przez
wykorzystanie płyt o zaślepionych otworach możliwa jest zmiana kierunku przepływu każdego
zczynników i wywołanie układu wielobiegowego. Liczba możliwych kombinacji jest
praktycznie nieograniczona, ale wszystkie rozwiązania można tak usystematyzować, że
wyróżnia się trzy główne rodzaje przepływu czynników przez aparat – tzn. układy szeregowe,
równoległe i mieszane (rys. 2).
W układzie szeregowym (rys. 2a) czynnik płynie kolejno przez poszczególne kanały. Droga
przepływu jest długa, ponieważ łączna długość strumienia jest sumą długości wszystkich
kanałów międzypłytowych. Układ ten cechuje się zatem dużym oporem hydraulicznym przy
stosunkowo niewielkim natężeniu przepływu, limitowanym wielkością pola powierzchni
poprzecznego przekroju kanału i dopuszczalną prędkością liniową.
Układ równoległy (rys. 2b) charakteryzuje jednakowy kierunek przepływu wszystkich
strumieni czynnika. Droga przepływu jest krótka, ponieważ odpowiada długości pojedynczego
kanału, zatem i spadek ciśnienia cieczy jest mały. Przez odpowiedni dobór liczby kanałów
równoległych można budować aparaty dostosowane do szerokiego zakresu natężeń przepływu
czynników. Jedynym aspektem limitującym jest średnica otworów wlotowych na płytach -
403572355.003.png 403572355.004.png
najczęściej przyjmuje się, że ze względu na opory hydrauliczne prędkość przepływu cieczy
w kolektorze nie powinna przekraczać 3 m/s. Równoległe połączenie kanałów powinno zawsze
występować w przypadkach stosowania nasyconej pary wodnej w charakterze nośnika ciepła.
Wówczas wlot pary powinien mieć miejsce górnym króćcem, zaś odbiór skroplin – z dolnego
kolektora.
W sytuacjach praktycznych nie zawsze spotyka się czyste układy szeregowe, bądź
równoległe. Zazwyczaj w jednym aparacie konieczne jest uzyskanie żądanego natężenia
przepływu czynników, zakresu zmian temperatur, a także często jako parametr limitujący określa
się dopuszczalny spadek ciśnienia jednego, lub obu czynników. Wówczas stosuje się układy
mieszane (rys. 2c-d), łączące cechy obu układów podstawowych. Odpowiednia liczba
równoległych kanałów składających się na wiązkę strumieni płynących w jednym kierunku
może zapewnić żądane natężenie przepływu, zaś szeregowe połączenie wiązek pozwala na taką
drogę przepływu czynnika, jaka jest niezbędna dla osiągnięcia zamierzonego zakresu zmian
temperatury. We wszystkich przykładach pokazanych na rys 2. powierzchnia wymiany ciepła
jest taka sama, zatem wybór konkretnego rozwiązania wymaga przeprowadzenia obliczeń
projektowych.
Obliczanie konfiguracji zestawu płyt
W obliczeniach projektowych płytowych wymienników ciepła trzeba uwzględnić
rzeczywiste wymiary kanału międzypłytowego traktując go w przybliżeniu jako układ płaski –
rys. 3. Wielkość pola powierzchni poprzecznego przekroju strumienia cieczy wynika
z rozwiniętej szerokości kanału b i średniej odległości między płytami h . W takim układzie ciecz
płynie ze średnią prędkością u , którą oblicza się z równania:
V
=
u
f
(1)
Nie wszyscy producenci wymienników płytowych podają graniczne wartości dopuszczalnych
prędkości u , ale zazwyczaj mieszczą się one w przedziale 0,2< u <0,6 m/s. Dolna granica wynika
z konieczności zachowania ciągłości strumienia przy przepływie z góry na dół, ponieważ przy
zbyt niskich obciążeniach przestrzeni międzypłytowej jest prawdopodobne, że prędkość
swobodnego spadania jest większa od prędkości wynikającej z prawa ciągłości. W takich
przypadkach przy zbyt niskiej prędkości liniowej trzeba liczyć się z ryzykiem rozrywania się
strumieni płynących w dół. Górna granica podanego przedziału ma związek ze znacznym
wzrostem oporów hydraulicznych ze wzrostem prędkości przepływu. Ze względu na sztuczną
turbulizację strumienia i znaczne straty ciśnienia cieczy w płytowych wymiennikach ciepła
prędkości przepływu zbliżające się do 1 m/s stosuje się raczej rzadko. Przy żądanej wydajności
i założonej prędkości przepływu otrzymuje się z równania (1) wielkość pola powierzchni
poprzecznego przekroju strumienia f, niezbędnego dla spełnienia przyjętych założeń. Przy znanej
wartości pola powierzchni przekroju pojedynczego kanału f
0 uzyskuje się liczbę strumieni
połączonych równolegle w jednej wiązce:
n =
f
(2)
f
0
Jeżeli wielkość f 0 nie jest podawana przez producenta, można ją obliczyć jako iloczyn średniej
odległości między płytami h i rozwiniętej szerokości pojedynczej płyty b 0 . Wielkość n
2/10
&
V &
zaokrągla się do najbliższej liczby całkowitej według ogólnie obowiązujących zasad, a następnie
dokonuje się korekty prędkości przepływu:
u
=
V
&
=
V
&
(3)
nf
nb
h
0
0
Skorygowana wg równania (3) prędkość u jest następnie wykorzystywana w dalszych
obliczeniach.
Obliczenie całego zestawu płyt wymiennika wymaga sporządzenia bilansu ciepła
Ciepło potrzebne do podgrzania czynnika:
Q
= &
V
ρ⋅
c
t
(4)
Ciepło przenoszone przez przeponę wymiennika:
Q
=
k
F
t
m
(5)
Równania (4) i (5) jednocześnie określają wydajność cieplną (obciążenie cieplne)
projektowanego aparatu. Proste przekształcenie układu równań (4) i (5) prowadzi do otrzymania
bilansu cieplnego w postaci bezwymiarowej
K
= &
kF
=
t
(6)
t
V
ρ
c
t
m
Wielkość K t zwana kryterium termicznym (ang. NTU = Number of Thermal Units) jest często
wykorzystywana w obliczeniach projektowych płytowych wymienników ciepła, zwłaszcza
w procedurach optymalizacji złożonych układów i stosowania inżynierskich algorytmów
opracowywanych oddzielnie dla każdego typu płyt [W ANG i S UNDEN , 2003]. Jak wynika
z definicji (6) określenie wielkości kryterium K t wymaga znajomości rozkładu temperatur
w projektowanym aparacie.
Niezależnie od sposobu sprecyzowania założeń wyjściowych sumaryczną powierzchnię
wymiany ciepła wyraża równanie
F
=
K
V
&
ρ
c
(7)
t
k
W celu uzyskania informacji o liczbie biegów wiązek strumieni w całym aparacie niezbędne jest
obliczenie sumarycznej długości strumienia L (rys. 3.). Wypada przy tym podkreślić, że ciecz
płynąca w przestrzeni międzypłytowej ograniczona jest ścianami płyt z dwóch stron (rys. 3.),
wobec czego
F
=
2
nb
0
L
(8)
stąd otrzymuje się łączną długość kanału
L
=
K
V
&
ρ
c
(9)
t
2
nb
k
0
3/10
403572355.005.png
Ponieważ jednak długość pojedynczego kanału międzypłytowego L 0 jest określona przez
geometrię płyty, to stosunek
L
i =
L
(10)
wyraża liczbę wiązek strumieni połączonych szeregowo. Ze względu na bezpieczeństwo
obliczeń wielkość i powinna być zawsze zaokrąglana „w górę” do najbliższej liczby całkowitej.
Obliczenia takie wykonuje się oddzielnie dla czynnika gorącego A otrzymując wielkości n A
oraz i A i dla czynnika chłodnego B - wielkości n B oraz i B . W rezultacie dane te składają się na
kompletny układ przepływu obu czynników w aparacie, który można przedstawić za pomocą
ułamka kompozycji
n
A
×
i
A
(11)
n
×
i
B
B
Przykłady stosowania ułamków kompozycji do uproszczonego przedstawienia układu strumieni
w wymienniku płytowym podano na rys. 2.
W zdecydowanej większości przypadków łączna liczba strumieni obu czynników jest taka sama.
W niektórych szczególnych zastosowaniach (np. w przemyśle spożywczym) praktykuje się
zasadę, że łączna liczba strumieni nośnika ciepła jest o 1 większa od sumy wszystkich strumieni
produktu. Liczba płyt składających się na wymiennik ciepła wynosi:
z
=
n
A
i
A
+
n
B
i
B
+
1
(12)
Obliczanie współczynnika przenikania ciepła
Współczynnik przenikania ciepła k jest wielkością kluczową w obliczeniach wymienników
płytowych, a jednocześnie dostępność danych, umożliwiających określenie wartości liczbowej k
dla konkretnego typu płyty jest dość ograniczona. Najczęściej bowiem producenci płyt
dysponują empirycznymi charakterystykami, których raczej się nie publikuje, lecz co najwyżej
są one udostępniane wraz z licencją tym firmom, które dostarczają kontrahentom kompletne
aparaty. Pomimo tego możliwe jest zastosowanie klasycznych metod obliczeniowych.
Obliczenie wartości współczynnika k wymaga znajomości wszystkich oporów cieplnych
występujących w układzie:
k
=
1
1
s
1
(13)
+
+
α
λ
α
A
B
Najmniej problematyczne jest określenie sumarycznego oporu przepony, jednak opory wnikania
ciepła w obu czynnikach wymagają przeprowadzenia obliczeń. Okazuje się bowiem, że w
wymiennikach płytowych intensywność ruchu ciepła jest ściśle związana zoporami
hydraulicznymi (T ARASOV , 1969, 1970). Ma to związek ze specyficzną aranżacją wytłoczeń na
płytach w taki sposób, że burzliwość wywołana w początkowym odcinku kanału nie wygasa,
4/10
403572355.006.png
lecz jest podtrzymywana przez następne. W rezultacie opór hydrauliczny takiego kanału jest
wprawdzie znacznie większy niż w kanale gładkim o tych samych wymiarach, ale jednocześnie
jest on mniejszy od sumy wszystkich oporów miejscowych na płycie.
Na podstawie analizy charakterystyk warunków hydraulicznych i cieplnych w wymiennikach
zbudowanych z płyt o różnej geometrii wytłoczeń zostały opracowane oraz opublikowane
zależności empiryczne, pozwalające na wykonywanie obliczeń projektowych aparatów
płytowych (T ARASOV , 1969, 1970).
Opory hydrauliczne w kanale międzypłytowym charakteryzuje współczynnik oporu
przepływu ζ u (umowny) zdefiniowany równaniem Darcy-Weisbacha
L
u
2
P
=
ζ
ρ
(14)
u
d
2
e
Równanie to może też być wykorzystywane do oszacowania wielkości spadku ciśnienia
w kanałach międzypłytowych.
Wartość współczynnika ζ u jest funkcją liczby Reynoldsa i odległości między płytami
C
0
,
65
+
1
07
lg
h
h
(15)
ζ
=
0
u
Re
n
gdzie wartości stałych C i n wynoszą: C = 340 i n = 0,85, gdy Re 1500 oraz C = 4,23
i n = 0,25, gdy Re 1500. Liczba Reynoldsa Re wprawdzie nie w pełni odzwierciedla charakter
ruchu czynnika w kanale międzypłytowym, ale stanowi ona ważną wielkość w procedurze
obliczeń cieplnych. Jako charakterystyczny wymiar liniowy we wszystkich obliczeniach
przyjmuje się podwójną odległość między płytami
d e =
2
h
(16)
Do obliczenia wartości liczbowych współczynników wnikania ciepła trzeba także wyznaczyć
współczynniki oporu przepływu w kanałach gładkich, korzystając ze znanych zależności:
ζ , gdy Re 2100 lub
96
ζ
=
0
3164
, gdy Re > 2100
(17)
Re
0
Re
0
25
Sztuczna turbulizacja strumienia cieczy wskutek obecności wytłoczeń sprawia, że opór
hydrauliczny kanału międzypłytowego jest większy w porównaniu z oporem kanału gładkiego
przy tej samej wartości liczby Reynoldsa. Ilościowo wyraża to współczynnik wymuszonej
burzliwości strumienia w kanale:
β
=
1
+
0
,
33
0
,
66
h
0
ln
ζ
u
(18)
t
d
ζ
e
0
Oba przeciwstawne zjawiska występujące w kanale międzypłytowym, tzn. sztuczna burzliwość,
wymuszana przez wytłoczenia na powierzchniach płyt i jednoczesne jej tłumienie w wyniku
5/10
,
0 =
,
403572355.001.png 403572355.002.png
Zgłoś jeśli naruszono regulamin