Wyznaczono pobory mocy w różnych bioreaktorach wielorazowych i jednorazowych o objętościach roboczych od 1 L do 10 L. Szczegóły geometryczne podsumowano w Tabeli 2. W przypadku naczyń jednorazowych, zamontowane u góry porty rurek oraz gniazda wałów mieszadeł musiały zostać usunięte z płyt górnych, aby umożliwić umieszczenie naczyń w uchwycie. Ponadto wbudowane plastikowe wały zostały przytwierdzone do wału ze stali nierdzewnej, który był stosowany wraz z łożyskiem powietrznym, jednak nie wymagano dalszych modyfikacji.
Moment obrotowy mierzono dla prędkości mieszadła od 100 rpm do 300 rpm w naczyniach bez przegród oraz od 100 rpm do 700 rpm w naczyniach z przegrodami, co odpowiada maksymalnym prędkościom liniowym końcówki łopatek wynoszącym odpowiednio 1,13 m·s-1 oraz 1,54 m·s-1 (patrz Równanie 4).
(4)
Zdefiniowane prędkości mieszadła w dolnym zakresie były ograniczone przez dokładność pomiaru czujnika momentu obrotowego oraz względne odchylenie standardowe powtarzalności wynoszące odpowiednio ± 0,2 % i < 0,05 % momentu nominalnego (określone przez producenta36). Ponadto maksymalne prędkości mieszadła były definiowane przez moment nominalny (0,2 Nm), w szczególności dla badanego zbiornika 10 L, oraz przez powstawanie wirów w naczyniach bez przegród. Aby zapobiec uszkodzeniu czujnika, maksymalny moment obrotowy podczas pomiaru określono na poziomie 60 % momentu nominalnego (0,12 Nm), a głębokość wiru ograniczono do około 20 mm na podstawie inspekcji wizualnej.
Przy zastosowaniu stopniowego zwiększania prędkości mieszadła obrotowego, typowy profil momentu obrotowego przedstawiono na Rysunku 2. Sygnał momentu obrotowego wzrastał przy każdym skokowym zwiększeniu prędkości obrotowej, zgodnie z przewidywaniami z równania 1. Wartości szczytowe sygnału momentu obrotowego obserwowano po każdej regulacji prędkości mieszadła, co można wyjaśnić początkowym przyspieszeniem cieczy oraz sterowaniem PID prędkością mieszadła. Pomiary quasi-stabilne uzyskano po około 1 min, w zależności od zastosowanej prędkości obrotowej i rodzaju mieszadła. Pozostałe fluktuacje wokół średniej czasowej wartości momentu obrotowego dla poszczególnych faz wynosiły zazwyczaj około 5 % wartości średniej dla większości badanych mieszadeł i prędkości mieszania.
Do dalszej oceny wykorzystano średnie wartości momentu obrotowego dla danej fazy, natomiast pominięto szczytowy moment obrotowy po każdej korekcie prędkości. Przy częstotliwości pomiaru 2 Hz zmierzone momenty obrotowe (TL) stanowiły średnią z co najmniej 240 punktów danych, co zapewniło wystarczającą pewność statystyczną; względne odchylenia standardowe tych wartości średnich były niższe niż 3 % dla większości punktów pomiarowych, co wskazuje na stabilne sygnały pomiarowe. Co ciekawe, odchylenia standardowe zazwyczaj malały wraz ze wzrostem prędkości mieszania, co sugeruje, że względne znaczenie wspomnianych fluktuacji zmniejsza się przy wyższym intensywnym mieszaniu.
Jak wykazano wcześniej35, moment obrotowy jałowy, tj. moment mierzony bez cieczy w naczyniu, który może wynikać ze strat tarcia w łożysku, uszczelnieniach i napędzie silnika lub z niewielkich odgięć bądź niewyważenia wału mieszadła (zwłaszcza w przypadku jednorazowych wałów plastikowych), może zostać znacznie zredukowany dzięki zastosowaniu łożyska powietrznego. Ogólnie wartości momentu jałowego dla mieszadeł ze stali nierdzewnej były mniejsze niż dla tych wykonanych z tworzywa sztucznego. Można to wyjaśnić większą sztywnością wałów stalowych, co skutkuje mniejszymi oscylacjami podczas obrotu. Dla większości zastosowanych mieszadeł pozostały wartości momentu jałowego z łożyskiem powietrznym były tak niskie jak 0,5 mN·m i w konsekwencji znajdowały się poniżej lub blisko rozdzielczości zastosowanego miernika momentu obrotowego (0,4 mN·m). Najwyższy pozostały moment jałowy zaobserwowano w bioreaktorze #6, w którym na dnie naczynia znajduje się opornik wału mieszadła. Podczas obrotu wał mieszadła kolidował z tym opornikiem, co można zaobserwować również podczas eksperymentów z hodowlą, co prowadzi do dodatkowego tarcia.
Jak widać na Rysunku 3, po obliczeniu mocy doprowadzanej na podstawie efektywnych momentów obrotowych (zgodnie z równaniem 1) i przedstawieniu ich jako funkcji liczb Reynoldsa (równanie 3), otrzymano indywidualne profile dla każdego z badanych mediów modelowych. Na każdej z tych krzywych moc doprowadzana rosła wraz ze wzrostem liczby Reynoldsa, a nachylenia były zbliżone do zależności PL
Re3. Korelację tę można wyprowadzić z równania 2 i równania 3, zakładając stałą liczbę mocy oraz średnicę mieszadła. Wynik ten uzyskano dla wszystkich testowanych mieszadeł przy R2 > 0,99.
Na podstawie uzyskanych eksperymentalnych danych dotyczących momentu obrotowego obliczono charakterystyki mocy wszystkich badanych mieszadeł zgodnie z równaniem 2 (patrz Rycina 4, Rycina 5, Rycina 6). Jako punkt odniesienia wykorzystano standardową turbinę Rushtona, której liczby mocy są dobrze udokumentowane w literaturze1,16,17,18,19,20,21,22,23,24,25. Jak wynika z Ryciny 4a, liczba mocy w mniejszym naczyniu o pojemności 2 L (bioreaktor nr 1) malała przy niskich liczbach Reynoldsa (100 < Re < ≈500) od P0 = 6,3 do P0 ≈ 3,3, po czym ponownie wzrosła powyżej Re ≈ 2000. W warunkach w pełni turbulentnych (Re > 104) uzyskano prawie stałą liczbę mocy P0 = 4,17±0,14. Porównywalną wartość P0 = 4,34±0,22 wyznaczono dla większego naczynia o objętości roboczej 10 L (bioreaktor nr 2), podczas gdy pewne rozbieżności między obiema skalami stwierdzono dla zakresu przejściowego 600 < Re < 104 (patrz Rycina 4a). Niemniej jednak, tendencje jakościowe w obu skalach były w pełni zgodne z danymi literaturowymi1,19, w których wyznaczono odpowiednio moc pobieraną przez pojedynczą turbinę Rushtona w objętościach roboczych 20 L1 oraz 40 L19. Należy zauważyć, że liczby mocy dla zakresu turbulentnego są o maksymalnie 25% niższe niż wartości podane w danych referencyjnych P0 ≈ 4,719 oraz P0 ≈ 5,51. Bezpośrednie porównanie jest jednak często utrudnione ze względu na zastosowanie różnych technik pomiarowych, a także rozbieżności w parametrach geometrycznych, w tym stosunku średnic (d/D), prześwicie nad dnem (zM/D) oraz geometrii dna zbiornika i przegród. Inni badacze wyznaczyli liczby mocy dla turbin Rushtona w naczyniach z przegrodami w zakresie od 3,6 do 5,9, w zależności od zastosowanej geometrii mieszadła i naczynia17,18,21,24,27,29,37,38. Zatem można uznać, że obecne wyniki były zadowalające.
Na Rysunku 4b porównano liczby mocy bioreaktorów nr 3 i nr 4, o objętościach roboczych odpowiednio 1 L i 2 L, dla szerokiego zakresu liczb Reynoldsa. Wartości P0 dwóch geometrycznie podobnych mieszadeł stale spadały w zakresie przejścia i stały się stałe (bioreaktor nr 3: P0 = 3.67±0.06; bioreaktor nr 4: P0 = 4.46±0.05) przy w pełni rozwiniętej turbulencji dla Re > 104, co jest kryterium wcześniej stwierdzonym dla turbiny Rushtona i innych mieszadeł38. Co ciekawe, zaobserwowano niemal stałe przesunięcie między obiema skalami, co można wyjaśnić różnicami w geometrii naczynia i mieszadła. Mimo że konfiguracja mieszadła w obu naczyniach jest podobna, nie można było zachować stałości wszystkich parametrów geometrycznych. Na przykład naczynie o objętości 1 L jest wyposażone tylko w dwa wbudowane przegrody, podczas gdy naczynie o objętości 2 L wyposażono w trzy przegrody. Powszechnie wiadomo, że liczba mocy rośnie wraz ze wzrostem liczby przegród, aż do osiągnięcia krytycznego warunku wzmocnienia38. Ponadto kształt tarczy mieszadła w mniejszym naczyniu musiał zostać zmodyfikowany ze względu na możliwości wykonania, co mogło mieć wpływ na dostarczaną moc. Należy również zauważyć, że zmierzone wartości momentu obrotowego w mniejszym naczyniu wynosiły od 4.2 mN·m do 12.8 mN·m, co stanowi zaledwie do 6 % nominalnego momentu obrotowego użytego miernika momentu. W tym zakresie niewielkie odchylenia sygnału pomiarowego mogą mieć znaczący wpływ na wyniki. Ponieważ nie dysponowano danymi porównawczymi z pomiarów referencyjnych, trudno jest wyciągnąć ostateczne wnioski na temat niezawodności pomiarów w najmniejszej skali zastosowanej w tej badaniu, co wymaga dalszych analiz.
Rysunek 5 przedstawia charakterystyki mocy trzech badanych, dostępnych komercyjnie bioreaktorów jednorazowych. W przeciwieństwie do naczyń z przegrodami, liczby mocy mieszadeł jednorazowych stale spadały w całym badanym zakresie liczb Reynoldsa (100 < Re < 3·104), a wartości stałe nie zostały osiągnięte z powodu postępującego tworzenia się wirów przy wysokich prędkościach mieszania w naczyniach bez przegród. Najwyższe liczby mocy, w zakresie od P0 ≈ 6 do P0 ≈ 1.8, uzyskano dla bioreaktora #5, który jest mieszany za pomocą radjalnego mieszadła łopatkowego oraz aksjalnego mieszadła segmentowego z łopatkami nachylonymi pod kątem 45°.
Zgodnie z oczekiwaniami, dla bioreaktora nr 7, mieszanego za pomocą dwóch mieszadeł segmentowych z łopatkami nachylonymi pod kątem 30°, co skutkuje przepływem głównie osiowym, uzyskano niższe liczby mocy w zakresie od P0 ≈ 5,1 do P0 ≈ 1,1. Powszechnie wiadomo, że mieszadła o przepływie osiowym mają mniejsze liczby mocy niż mieszadła o przepływie promieniowym ze względu na niższe opory przepływu łopatek nachylonych38. Należy zauważyć, że dane eksperymentalne dotyczące dostarczonej mocy w bioreaktorze nr 7, które zgłoszono wcześniej32, są nieco wyższe (np. P0 = 1,9 dla Re = 1,4·104). Jednakże wcześniej opublikowane dane wykazały tę samą zależność P0
Re-0.336, co stwierdzono również w niniejszym badaniu. Za różnice w wartościach bezwzględnych mogą odpowiadać odmienne techniki pomiarowe.
Wśród badanych bioreaktorów jednorazowych bioreaktor nr 6, mieszany przez jedno mieszadło morskie umieszczone blisko dna, wykazał najniższe liczby mocy w zakresie P0 ≈ 0,8 oraz P0 ≈ 0,3 (patrz Rycina 5). Tą niską moc wejściową można wyjaśnić małym kątem nachylenia łopatek mieszadła, mimo że analiza numeryczna dynamiki płynów (CFD) wykazała dość dominującą radialną składową przepływu wokół łopatek mieszadła39. Można stwierdzić dobrą zgodność obecnych wyników z opublikowanymi danymi z modeli CFD39 oraz eksperymentów32.
W końcu układ pomiarowy wykorzystano do zbadania wpływu średnicy mieszadła i kąta łopat w bioreaktorze nr 7. Jak wynika z Rysunku 6, wszystkie krzywe mocy spadają w sposób ciągły w całym zakresie liczb Reynoldsa, zgodnie z oczekiwaniami. Odnotowano istotne różnice pomiędzy dwoma kątami łopat (30° i 45°), gdzie większy kąt łopat wiązał się z wyższym poborem mocy (dla 30°: 1,13 < P0 < 4,25 oraz dla 45°: 1,65 < P0 < 4,46), niezależnie od turbulencji (tj. liczby Reynoldsa). Jest to znane również w przypadku klasycznych mieszadeł z łopatami skośnymi40 i ponownie można to wyjaśnić wyższym oporem przepływu wokół łopat o większym skosie. Co ciekawe, nie wykryto istotnych różnic w liczbach mocy pomiędzy dwiema średnicami mieszadeł. Stwierdzono to również dla mieszadeł z łopatami skośnymi, podczas gdy liczby mocy mieszadeł o przepływie promieniowym zazwyczaj mają tendencję do spadku wraz ze wzrostem stosunku d/D40.

Rycina 1: Schemat układu testowego. Układ składa się z (1) zbiornika mieszalnika, (2) uchwytu naczynia, (3) łożyska z tuleją powietrzną, (4) miernika momentu obrotowego, (5) napędu silnika, (6) przetwornika A/D, (7) jednostki sterującej, (8) komputera PC do akwizycji danych i sterowania. Zgodnie z zaleceniami producenta, do tulei powietrznej doprowadzono sprężone powietrze (5.5 bar). Wskazano również główne wymiary geometryczne zbiornika mieszalnika i mieszadła. Rycina została zmodyfikowana na podstawie opracowania35. Kliknij tutaj, aby wyświetlić większą wersję tej ryciny.

Rysunek 2: Typowy profil pomiarowy ze stopniowym zwiększaniem prędkości obrotowej mieszadła (tj. N1 < N2 < N3) w odstępach 5 min, zaznaczonych pionowymi przerywanymi liniami. Poziome przerywane linie reprezentują 5 % przedział ufności wokół średnich czasowych wartości momentu obrotowego dla odpowiadających im faz (oznaczonych poziomymi liniami ciągłymi). Wartości szczytowe zaobserwowano podczas pierwszej minuty każdego przedziału, co można wyjaśnić początkowym przyspieszeniem cieczy wewnątrz zbiorników oraz sterowaniem prędkością mieszadła w oparciu o regulator PID. Do dalszej ewaluacji wykorzystano jedynie sygnał momentu obrotowego podczas fazy quasi-stabilnej, w której sygnał pomiarowy oscylował wokół średniej wartości w ramach 5 % przedziału ufności. Kliknij tutaj, aby zobaczyć powiększoną wersję tego rysunku.

Rycina 3: Obliczona moc wejściowa w bioreaktorze nr 1 w funkcji liczby Reynoldsa dla różnych mediów modelowych. Osobne profile uzyskano dla każdego z badanych mediów modelowych. Linie ciągłe reprezentują przewidywania modelu przy założeniu P
Re3, co wykazało bardzo dobrą zgodność z danymi eksperymentalnymi (przy R2 > 0,99). Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 4: Wyznaczone liczby mocy w funkcji liczby Reynoldsa w zbiornikach z przegrodami. (a) Porównanie danych dla turbin Rushtona w małym i dużym zbiorniku (z objętością roboczą odpowiednio 2 L i 10 L) pokazuje, że bezwymiarowe liczby mocy dla warunków w pełni turbulentnych są równe dla obu skal. Odnotowano niewielkie odchylenia w zakresie przejściowym przy Re < 104, gdzie liczba mocy rosła wraz ze wzrostem liczby Reynoldsa. (b) Porównanie danych z bioreaktorów nr 3 i nr 4 wykazuje jakościowo podobny spadek liczb mocy wraz ze wzrostem liczby Reynoldsa, aż do uzyskania stabilnych wartości w warunkach w pełni turbulentnych. Liczby mocy dla bioreaktora 1 L wykazują większe fluktuacje w porównaniu z odpowiednikiem 2 L. Dla naczynia 1 L nie uzyskano danych dla liczb Reynoldsa w zakresie 550 < Re < 950 przy zastosowaniu tego samego medium modelowego co w naczyniu 2 L. Ilościowa różnica między skalami może wynikać z różnic w geometrii naczynia i mieszadła lub może być wynikiem czułości czujnika. Wymagane są dalsze badania. Linie ciągłe reprezentują modele regresji wielomianowej. Kliknij tutaj, aby zobaczyć powiększoną wersję tego rysunku.

Rysunek 5: Wyznaczone liczby mocy w funkcji liczby Reynoldsa dla różnych jednorazowych bioreaktorów. Liczby mocy dla każdego z naczyń malały wraz ze wzrostem liczb Reynoldsa. W przeciwieństwie do naczyń z przegrodami, w naczyniach bez przegród nie uzyskano stabilnych liczb mocy ze względu na postępujące tworzenie się wirów przy wysokich prędkościach mieszania. Linie ciągłe reprezentują modele regresji wielomianowej. Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.

Rysunek 6: Wyznaczone liczby mocy w funkcji liczby Reynoldsa dla różnych modyfikacji bioreaktora #7. Uzyskano odrębne profile dla dwóch różnych kątów łopatek wynoszących 30° i 45°, jednak nie stwierdzono istotnych różnic między dwoma stosunkami średnic mieszadła (d/D = 0.43 i d/D = 0.57). Liczby mocy dla wszystkich konfiguracji wykazywały ciągły spadek w całym badanym zakresie liczb Reynoldsa z powodu postępującego tworzenia się wirów przy wysokich prędkościach mieszania w naczyniach bez przegród. Linie ciągłe reprezentują modele regresji wielomianowej. Kliknij tutaj, aby zobaczyć powiększoną wersję tego rysunku.
| Końcowe stężenie sacharozy | Gęstość cieczy ρL | Lepkość cieczy ηL | Liczba Reynoldsa Re |
| (%w/w) | (kg·m-3) | (mPa·s) | (-) |
| 0 | 998.2 | 1 | 11954 |
| 20 | 1081 | 2 | 6486 |
| 30 | 1127 | 3.2 | 4226 |
| 40 | 1176.4 | 6.2 | 2277 |
| 50 | 1231.7 | 15.5 | 954 |
| 55 | 1259.8 | 28.3 | 534 |
| 60 | 1288.7 | 58.9 | 263 |
Tabela 1: Podsumowanie gęstości i lepkości cieczy dla wybranych roztworów sacharozy w 20 °C oraz wynikowa bezwymiarowa liczba Reynoldsa dla mieszadła o średnicy 60 mm i prędkości obrotowej 200 rpm. Liczba Reynoldsa została obliczona zgodnie z równaniem 3.

Tabela 2: Podsumowanie szczegółów geometrycznych badanych bioreaktorów. Kliknij tutaj, aby pobrać ten plik.