Name: Szybkie obrazowanie cząstek w pobliżu powierzchni (Wideo) | JoVE Journal
Uploaded: 2013-06-24T14:00:00.000+00:00
Duration: 11 min 59 s
Description: Poznaj metodę PIV do analizy przepływów płynów w pobliżu powierzchni przy użyciu obrazowania cząstek. Technika ta umożliwia nieinwazyjne pomiary wielowymiarowych pól prędkości z wysoką rozdzielczością czasową i przestrzenną.

ai_agent_IP

ai_agent

maintenance

instantTrialDurationMinutes

isBiopharma

enableNewMyJoveSectionDesign

educationFeatures

seoMetaApi

ai_copilot_institution_blacklist

ai_chat_institution_blacklist

freeAccessWindow

ssoRoleCapture

scheduledMaintenanceBanner

joveLabsFeatures

disable_autoplay_for_bots

extended_trial_duration

5 minutes

playerTrialDesktop

Ten materiał jest oparty na pracy wspieranej przez Amerykańską Narodową Fundację Nauki w ramach grantu nr. CBET-1032930 i praca wykonana w Laboratorium Ilościowej Diagnostyki Laserowej Uniwersytetu Michigan.

1. Bezpieczeństwo w laboratorium
<ol>
	<li>Przed uruchomieniem lasera zapoznaj się z materiałami dotyczącymi bezpieczeństwa pracy z laserami i upewnij się, że wymogi dotyczące szkolenia zostały spełnione.</li>
	<li>Zdobądź odpowiedni sprzęt ochronny do pracy z laserami. Każda osoba powinna nosić okulary ochronne, które blokują długość (lub długości) fali emisji lasera. Umieść znak ostrzegawczy przed laboratorium, aby informować innych o pracy lasera. Rozwieś zasłony ochronne wokół stołu optycznego, aby odizolować go od innych współpracowników w wspólnej przestrzeni laboratoryjnej.</li>
	<li>Zdejmij wszystkie zegarki i biżuterię podczas pracy z laserami.</li>
	<li>Uwzględnij drogę wiązki podczas ustawiania sprzętu: rozmieść urządzenia w taki sposób, aby dokonywanie regulacji nie wymagało sięgania nad lub pod wiązkę.</li>
	<li>Przeczytaj instrukcję obsługi lasera, aby dowiedzieć się, jak bezpiecznie go eksploatować.</li>
	<li>Utrzymuj poziom oczu poza płaszczyzną wiązki lasera!</li>
</ol>
2. Przygotowanie stanowiska pracy
<ol>
	<li>Wyznaczyć powiększenie wymagane do danej aplikacji i dobrać odpowiednią soczewkę. Powiększenie (M) można wyznaczyć, dzieląc długość matrycy kamery przez odpowiadającą jej długość pola widzenia (FOV). W tym przykładzie długość matrycy kamery wynosi 17.6 mm, a odpowiadająca jej długość FOV to 2.4 mm. Zatem M = 17.6 mm / 2.4 mm = 7.33. Aby uzyskać tak małe FOV, zastosowano mikroskop dalekosiężny.
	<ol>
		<li>Przeprowadzić wstępne obliczenia przewidywanych prędkości w obszarze przyściennym. Na podstawie tych szacunków wyznaczyć parametry rejestracji, takie jak liczba klatek na sekundę oraz opóźnienie czasowe, zgodnie z praktycznymi wytycznymi dla PIV 1,2. Wyznaczyć czas, w jakim cząsteczka pokona odległość 8 pikseli. To wyznaczy opóźnienie czasowe między kolejnymi impulsami lasera (dt). W PIV serii czasowej wartość 1/dt wyznaczy wymaganą liczbę klatek na sekundę kamery i musi być ona mniejsza niż maksymalna liczba klatek dopuszczalna przez kamerę. W późniejszym etapie mogą być konieczne niewielkie korekty tych parametrów w celu optymalizacji rejestracji przepływu, aby uzyskać wysokiej jakości dane o prędkości. Jeśli wymagana liczba klatek na sekundę przekracza maksymalną częstotliwość repetycji lasera, można użyć dwóch laserów do przeprowadzenia PIV w trybie rozdzielania klatek (frame-straddling). W tym przykładzie liczba klatek na sekundę (5 kHz) nie przekracza maksymalnej częstotliwości repetycji lasera, zatem do przeprowadzenia PIV w trybie serii czasowej wymagany jest tylko jeden laser.</li>
	</ol>
	</li>
	<li>Wypoziomować laser względem stołu
	<ol>
		<li>Ustawić głowicę lasera na jednym z końców wypoziomowanego stołu optycznego. Umieścić pochłaniacz wiązki bezpośrednio na drodze wiązki na drugim końcu stołu.</li>
		<li>Umieścić szynę optyczną pomiędzy głowicą lasera a pochłaniaczem wiązki. Przymocować tarczę do blokera wiązki, przymocować bloker wiązki do uchwytu, a uchwyt umieścić na szynie.</li>
		<li>Ustawić prąd lasera na niską wartość – wystarczającą do emisji, ale niewystarczającą do przepalenia kartki papieru. Włączyć laser i przesuwać uchwyt w przód i w tył. Dokonywać niewielkich korekt pozycji lasera, aż centrum wiązki laserowej pozostanie w jednym punkcie podczas przesuwania uchwytu. Przymocować laser do stołu optycznego.</li>
		<li>Zmierzyć wysokość centrum wiązki laserowej za pomocą kątownika kombinowanego. Wyłączyć laser.</li>
	</ol>
	</li>
	<li>Zainstalować optykę do kształtowania arkusza laserowego
	<ol>
		<li>Zdjąć szynę, ale pozostawić bloker wiązki z tarczą przed pochłaniaczem wiązki. Włączyć laser i starannie zaznaczyć miejsce, w którym centrum wiązki uderza w tarczę. Umieścić w drodze lasera optykę do kształtowania arkusza, którą w tej demonstracji jest homogenizator wiązki (BH) zawierający również teleskop do kształtowania arkusza, aby utworzyć arkusz laserowy. Wysokość arkusza laserowego musi być większa niż FOV. Dostosować pozycję BH tak, aby wycentrować wysokość i szerokość arkusza laserowego względem znaku na tarczy oraz zapobiec powrotowi odbić do wnęki lasera. Jeśli jest to konieczne, umieścić przysłonę między głowicą lasera a BH, aby uniknąć odbić wstecznych. Wyłączyć laser. Arkusz światła w tej demonstracji miał wysokość 8 mm i grubość 0.5 mm oraz energię impulsu 0.4 mJ/impuls.</li>
		<li>Jeśli miejsce na stole optycznym jest ograniczone, umieścić lustro o wysokim współczynniku odbicia pod kątem 45°, aby obrócić arkusz światła laserowego o 90°. Przymocować kolejną tarczę do blokera wiązki, przymocować bloker wiązki do uchwytu i umieścić uchwyt na szynie. Umieścić zespół szynowy za lustrem. Włączyć laser. Dokonywać niewielkich korekt ustawienia lustra, aż centrum arkusza światła pozostanie w jednym miejscu na tarczy podczas przesuwania jej wzdłuż szyny.</li>
		<li>Ustawić częstotliwość repetycji lasera zgodnie z liczbą klatek na sekundę dla pomiarów (5 kHz dla omawianego przykładu) i ustawić prąd lasera na wartość maksymalną. Umieścić szynę między BH a tarczą. Przymocować drugi bloker wiązki do uchwytu i umieścić zespół na szynie. Włączyć laser. Przesuwać uchwyt w przód i w tył, aby wyznaczyć położenie punktu ogniskowego względem BH. Zaznaczyć położenie punktu ogniskowego względem BH. W przypadku użycia lustra pomiar wykonać względem lustra. Zmierzyć przybliżoną wysokość arkusza laserowego w punkcie ogniskowym. Wyłączyć laser.</li>
	</ol>
	</li>
	<li>Zamontować i wyregulować mikroskop dalekosiężny oraz kamerę
	<ol>
		<li>Wyznaczyć poziome i pionowe linie środkowe apertur mikroskopu dalekosiężnego (LDM) i kamery za pomocą kwadratu centrującego i kątownika kombinowanego. Zmierzyć odległość między stołem a poziomymi liniami środkowymi LDM i kamery.</li>
		<li>Przymocować LDM i kamerę do uchwytów, używając podkładek lub nakrętek, tak aby poziome linie środkowe LDM i kamery znajdowały się na tej samej wysokości. Zamocować LDM i kamerę na szynie. Połączyć LDM i kamerę za pomocą odpowiednich adapterów. Dostosować wysokość zespołu tak, aby poziome linie środkowe znajdowały się w tej samej odległości nad stołem co centrum arkusza światła.</li>
		<li>Zamocować stół translacyjny przed zaznaczonym punktem ogniskowym wiązki. Ruch stołu translacyjnego będzie równoległy do propagacji wiązki. Przymocować szynę z zespołem kamery do stołu translacyjnego tak, aby cały zespół był prostopadły do arkusza światła. Wycentrować zespół kamery, wyrównując pionowe linie środkowe LDM i kamery z punktem ogniskowym.</li>
		<li>Podłączyć kamerę do komputera i kontrolera wysokiej prędkości (HSC). Podłączyć laser do HSC. Pozostawić osłonę zespołu kamery i przeprowadzić kalibrację intensywności w programie PIV (LaVision DaVis 7.2).</li>
		<li>W programie ustawić kamerę w trybie ciągłego przechwytywania (continuously grab) i zdjąć osłonę zespołu kamery. Umieścić kątownik kombinowany w punkcie ogniskowym. Przesuwać kamerę i LDM wzdłuż szyny, aż uzyska się wyraźny obraz linijki. Kontynuować przesuwanie kamery i LDM wzdłuż szyny i ustawić ostrość obrazu za pomocą pręta ostrzącego LDM, aż matryca kamery obejmie pożądane pole widzenia (2.4 x 1.8 mm2 odpowiadające matrycy 800 x 600 pikseli).</li>
	</ol>
	</li>
	<li>Przymocować płytkę do uchwytu tak, aby była równoległa do stołu i umieścić ją w punkcie ogniskowym. Podnieść płytkę tak, aby była widoczna na obrazach na komputerze. Wyłączyć ciągłe przechwytywanie i założyć osłonę zespołu kamery. Włączyć laser i upewnić się, że arkusz światła laserowego styka się z powierzchnią płytki na całej jej długości.</li>
</ol>
3. Konfiguracja przepływu
<ol>
	<li>W tej demonstracji PIV jest wykonywane poprzez rejestrację obrazów światła rozproszonego na kropelkach oleju silikonowego. Kropelki oleju są wytwarzane za pomocą atomizera oleju. Podłącz następujące elementy do zasilania powietrzem: filtr cząsteczkowy, filtr oleju, regulator ciśnienia, przepływomierz masowy oraz atomizer oleju. Wyjście atomizera podłącz do stalowej rurki. Użyj uchwytu i zacisku, aby przymocować stalową rurkę do stołu optycznego, unieś rurkę nad stół i skieruj ją w stronę płyty.</li>
	<li>Otwórz zawór zasilania powietrzem. Ustaw ciśnienie wsteczne na regulatorze ciśnienia na &gt; 140 kPa, aby zapewnić odpowiedni przepływ przez układ.</li>
	<li>Włącz przepływ i dostosuj gęstość zasiań za pomocą dysz atomizera oraz zaworów obejściowych w atomizerze.</li>
</ol>
4. Optymalizacja konfiguracji
<ol>
	<li>Wprowadź częstotliwość klatek w programie komputerowym. Sprawdź, czy HSC wysyła do lasera sygnał wyzwalający zgodny z częstotliwością klatek. W zasilaczu lasera ustaw częstotliwość repetycji oraz prąd (w tym przykładzie odpowiednio 5 kHz i 15,5 A). Ustaw laser w trybie zewnętrznym. Przed przełączeniem na tryb zewnętrzny laser musi stale otrzymywać z HSC sygnał wyzwalający zgodny z ustawioną w laserze częstotliwością repetycji, w przeciwnym razie laser ulegnie przegrzaniu.</li>
	<li>Ustaw kamerę w tryb ciągłego przechwytywania obrazu (continuously grab), włącz laser oraz atomizer. Użyj pręta do ustawiania ostrości w LDM, aby upewnić się, że obrazy cząstek są ostre. Upewnij się również, że intensywność obrazów cząstek nie powoduje nasycenia kamery. Jeśli tak się dzieje, zmniejsz prąd lasera – wpłynie to na położenie punktu ogniskowania! W przypadku zmiany prądu lasera powtórz kroki 2.3.3 i 2.4.3. Wyłącz tryb przechwytywania obrazu po uzyskaniu ostrych obrazów cząstek.</li>
	<li>Rejestracja, analiza i dostosowanie parametrów w celu uzyskania poprawnych danych o prędkości
	<ol>
		<li>Zarejestruj kilkaset obrazów przepływu. Po zakończeniu rejestracji sprawdź zapisane obrazy, aby upewnić się, że przesunięcie cząstek nie przekracza 8 pikseli, gęstość zasiań jest rzędu 8-10 cząstek na okno analizy o wymiarach 32 x 32 piksele oraz aby zweryfikować ostrość obrazów. Powtarzaj kroki 4.3.1-4.3.4 aż do spełnienia powyższych kryteriów.</li>
		<li>Jeśli przesunięcie cząstek przekracza 8 pikseli, zmniejsz wartość dt pomiędzy dwoma impulsami lasera PIV, aby uzyskać maksymalne przesunięcie wynoszące 8 pikseli. Jeśli przesunięcie cząstek jest znacznie mniejsze niż 8 pikseli, odpowiednio zwiększ dt. W systemach PIV z jednym laserem wartość dt reguluje się poprzez zmianę częstotliwości klatek, a w konsekwencji częstotliwości repetycji lasera. W systemach PIV z dwoma laserami dt jest opóźnieniem czasowym między impulsem z pierwszego lasera a impulsem z drugiego lasera. Jeśli regulacja dt nie rozwiązuje problemu, można najpierw dostosować częstotliwość klatek i repetycji lasera, a następnie ponownie precyzyjnie ustawić dt.</li>
		<li>Jeśli śledzenie grup cząstek w serii obrazów jest utrudnione, może to wynikać ze zbyt dużego ruchu poza płaszczyzną obrazowania. Istnieje kilka sposobów na rozwiązanie tego problemu: a) przesuń zespół kamery względem punktu ogniskowania tak, aby kamera obrazowała grubszą płaszczyznę światła; b) zwiększ odległość roboczą między kamerą a zespołem płaszczyzny światła (i ustaw ostrość za pomocą pręta do ogniskowania), aby uzyskać większą głębię ostrości, jednak zmniejszy to rozdzielczość przestrzenną.</li>
		<li>Jeśli gęstość zasiań jest zbyt mała lub zbyt duża, zwiększ lub zmniejsz liczbę dysz atomizera.</li>
	</ol>
	</li>
</ol>
5. Przeprowadzenie eksperymentu
<ol>
	<li>Przeprowadź kalibrację intensywności kamery z założoną osłoną zespołu kamery, aby ustawić punkt odniesienia dla intensywności. Po zakończeniu kalibracji zdejmij osłonę.</li>
	<li>Ustaw laser na zoptymalizowaną częstotliwość repetycji i prąd. Przed przełączeniem lasera w tryb zewnętrzny upewnij się, że laser otrzymuje ciągły sygnał wyzwalający zgodny z ustawioną częstotliwością. Włącz laser.</li>
	<li>Nagraj sekwencję obrazów tła, na których widoczna jest jedynie płaszczyzna światła muskająca powierzchnię płytki. Zapisz te obrazy.</li>
	<li>Włącz przepływ i poczekaj, aż się ustabilizuje.</li>
	<li>Ustaw kamerę w trybie ciągłego przechwytywania obrazu (continuous grab) i zweryfikuj, czy kamera rejestruje ostre obrazy cząstek. Wyłącz tryb ciągłego przechwytywania.</li>
	<li>Wprowadź żądaną liczbę obrazów, a następnie naciśnij przycisk nagrywania.</li>
	<li>Po zakończeniu nagrywania wyłącz przepływ i laser. Przejrzyj sekwencję obrazów i sprawdź przesunięcie cząstek, gęstość zasiewu oraz ostrość obrazu cząstek. Jeśli wynik jest satysfakcjonujący, zapisz nagranie; w przeciwnym razie powtórz kroki 5.4-5.7.</li>
	<li>Powtórz kroki 5.4-5.7, aby zebrać więcej serii pomiarowych.</li>
	<li>Zwiększ czas ekspozycji kamery (czas zbierania obrazów dla każdej klatki).</li>
	<li>Umieść wzorzec kalibracyjny w płaszczyźnie światła i upewnij się, że styka się on z płytką. Oświetl wzorzec od tyłu za pomocą źródła światła (np. latarki). Przy włączonym trybie ciągłego przechwytywania obrazu ustaw wzorzec tak, aby zarejestrowany obraz był ostry i niedodstrojony. Upewnij się, że punkt styku między płytką a wzorcem jest widoczny na obrazie – jest to kluczowe dla określenia położenia płytki na obrazach.</li>
	<li>Nagraj 10 obrazów wzorca kalibracyjnego. Powtarzaj kroki 5.9-5.11 za każdym razem, gdy zostanie zmieniony zespół kamery lub ustawienie ostrości.</li>
</ol>
6. Przetwarzanie danych
<ol>
	<li>Program PIV wykorzystany w tej demonstracji to LaVision DaVis 8.1. Należy uśrednić każdy zestaw obrazów tarczy kalibracyjnej. Uzyskany obraz należy wykorzystać w procedurze kalibracji w celu określenia rzeczywistych wymiarów nabytych obrazów.</li>
	<li>Należy zastosować każdą kalibrację do odpowiadającego jej zestawu obrazów.</li>
	<li>Należy określić położenie płyty na skalibrowanych obrazach. Informacja ta jest niezbędna do stworzenia maski geometrycznej (opisanej w sekcji 6.6).</li>
	<li>Należy uśrednić obrazy tła. Należy określić, czy odbicia laserowe od powierzchni znacząco wpływają na szum tła, porównując wartości intensywności uśrednionego obrazu tła z wartościami intensywności cząsteczek znacznika. Jasne odbicia laserowe w pobliżu ścianki będą miały intensywność wyższą niż intensywność cząsteczek. Będzie to niekorzystnie wpływać na korelacje PIV w pobliżu ścianki i ograniczy położenie pierwszego wiarygodnego wektora najbliższego ścianki. W tym przykładzie odbicia laserowe nie miały znaczącego wpływu na tło.</li>
	<li>Obrazy przepływu po kalibracji należy poddać wstępnej obróbce za pomocą filtra górnoprzepustowego (filtr odejmowania przesuwającego się tła), aby usunąć duże fluktuacje intensywności tła, takie jak odbicia laserowe. Sygnały cząsteczek charakteryzują się niewielkimi fluktuacjami intensywności i przejdą przez filtr.</li>
	<li>Należy zdefiniować maskę geometryczną – należy użyć maski prostokątnej, aby wyłączyć obliczanie wektorów w miejscu, w którym płyta znajduje się na obrazach. Uwaga: program DaVis oferuje dwie opcje masek geometrycznych: jedną, która włącza korelacje PIV w określonym obszarze, oraz drugą, która wyłącza korelacje PIV w określonym obszarze. W tej demonstracji zastosowano maskę włączającą algorytm PIV w określonym obszarze.</li>
	<li>W menu „Advanced mask settings” należy upewnić się, że maska jest zastosowana odpowiednio (tzn. używane są tylko piksele wewnątrz maski).</li>
	<li>Należy określić procedurę obliczania wektorów: w tym przykładzie zastosowano procedurę wieloprzebiegową ze zmniejszającym się rozmiarem okna – 2 pierwsze przejścia z oknami przesłuchującymi o rozmiarze 64 x 64 pikseli i 50% zakładką, a następnie 3 przejścia z oknami przesłuchującymi o rozmiarze 32 x 32 pikseli i 50% zakładką.</li>
	<li>Pola wektorów prędkości w tej demonstracji zostały poddane postprocesowaniu za pomocą pięciu procedur w celu poprawy jakości wyników korelacji wzajemnej: a) trwałe zastosowanie maski; b) usunięcie wektorów z współczynnikiem szczytu (Q) &lt; 1,1; c) zastosowanie filtra medianowego; d) usunięcie grup zawierających &lt; 5 wektorów; e) zastosowanie uzupełniania wektorów (vector fill-up). Współczynnik szczytu (Q) definiuje się jako <img alt="Wzór równowagi Q=((P1−min)/(P2−min)); równanie matematyczne; użytek edukacyjny." fo:src="/files/ftp_upload/50559/50559eq1highres.jpg" src="/files/ftp_upload/50559/50559eq1.jpg" data-altupdatedat="1787630408000" data-alt="Równanie">, gdzie P1 i P2 to odpowiednio pierwszy i drugi najwyższy szczyt korelacji, a min to wartość minimalna w płaszczyźnie korelacji. Q jest metryką służącą do oceny jakości wektora. Q porównuje najwyższy szczyt korelacji, który daje najlepszy wektor, z powszechnym tłem korelacji reprezentowanym przez drugi najwyższy szczyt korelacji. Wektory z wartością Q bliską 1 wskazują, że najwyższy szczyt korelacji jest szczytem fałszywym. Następnie filtr medianowy wyznacza wektor medianowy (umedian, vmedian) grupy wektorów oraz odchylenie wektorów sąsiednich (urms, vrms). Filtr medianowy odrzuca wektor środkowy (u,v), jeśli nie spełnia on następujących kryteriów: umedian - urms ≤ u ≤ umedian + urms oraz vmedian - vrms ≤ v ≤ vmedian + vrms. Ponadto, jeśli przy obliczaniu wektorów prędkości określono dużą zakładkę, mogą powstać grupy wektorów pozornych. Dlatego możliwe jest usunięcie grup zawierających mniej niż określoną liczbę wektorów. Po usunięciu wektorów pozornych można zastosować uzupełnianie wektorów, aby wypełnić puste przestrzenie wektorami interpolowanymi, wyznaczonymi z niezerowych wektorów sąsiednich. Na koniec trwałe zastosowanie maski usunie wszystkie wektory znajdujące się poza nią.</li>
	<li>Należy ocenić jakość wyników: a) Czy wyniki mają sens fizyczny? (tzn. mniejsze prędkości w pobliżu granicy, wzrost prędkości wraz ze wzrostem odległości od ścianki; kierunek wektorów zgodny z ogólnym kierunkiem przepływu, itp.); b) Czy wynikowe pole wektorowe składa się w przeważającej części z wektorów pierwszego wyboru (wskazanych przez oprogramowanie PIV). Zazwyczaj zaleca się, aby udział wektorów pierwszego wyboru był wyższy niż 95%. Szerszy zakres kroków postprocesowania opisano w literaturze, np. 1,2.</li>
</ol>

<ol><li>Raffel, M. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aRaffel%2C+M+%22Particle+Image+Velocimetry%3A+A+Practical+Guide%22">Particle Image Velocimetry: A Practical Guide</a>. , Springer. Berlin. (2007).</li>
<li>Adrian, R. J., Westerweel, J. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aAdrian%2C+RJ+%22Particle+image+velocimetry%22">Particle image velocimetry</a>. , Cambridge University Press. New York. (2011).</li>
<li>Sick, V. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((High+speed+imaging+in+fundamental+and+applied+combustion+research%5BTitle%5D)+AND+%22Proceedings+of+the+Combustion+Institute%22%5BJournal%5D)">High speed imaging in fundamental and applied combustion research.</a> Proceedings of the Combustion Institute. 34 (2), 3509-3530 (2013).</li>
<li>Prasad, A. K. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Particle+image+velocimetry%5BTitle%5D)+AND+%22Current+Science+%28Bangalore%29%22%5BJournal%5D)">Particle image velocimetry.</a> Current Science (Bangalore). 79 (1), 51(2000).</li>
<li>Stitou, A., Riethmuller, M. L. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Extension+of+PIV+to+super+resolution+using+PTV%5BTitle%5D)+AND+%22Measurement+Science+%26+Technology%22%5BJournal%5D)">Extension of PIV to super resolution using PTV.</a> Measurement Science & Technology. 12 (9), 1398-1403 (2001).</li>
<li>Alharbi, A. Y., Sick, V. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Investigation+of+boundary+layers+in+internal+combustion+engines+using+a+hybrid+algorithm+of+high+speed+micro-PIV+and+PTV%5BTitle%5D)+AND+%22Experiments+in+Fluids%22%5BJournal%5D)">Investigation of boundary layers in internal combustion engines using a hybrid algorithm of high speed micro-PIV and PTV.</a> Experiments in Fluids. 49 (4), 949-959 (2010).</li>
<li>Meinhart, C. D., Wereley, S. T., Santiago, J. G. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+PIV+Algorithm+for+Estimating+Time-Averaged+Velocity+Fields%5BTitle%5D)+AND+%22Journal+of+Fluids+Engineering%22%5BJournal%5D)">A PIV Algorithm for Estimating Time-Averaged Velocity Fields.</a> Journal of Fluids Engineering. 122 (2), 285(2000).</li>
<li>Jainski, C., Lu, L., Dreizler, A., Sick, V. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((High-Speed+Micro+Particle+Image+Velocimetry+Studies+of+Boundary-Layer+Flows+in+a+Direct-Injection+Engine%5BTitle%5D)+AND+%22International+Journal+of+Engine+Research%22%5BJournal%5D)">High-Speed Micro Particle Image Velocimetry Studies of Boundary-Layer Flows in a Direct-Injection Engine.</a> International Journal of Engine Research. , (2012).</li>
<li>White, F. M. Fluid mechanics. , McGraw-Hill. New York. 864(2008).</li>
<li>Pope, S. B. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aPope%2C+SB+%22Turbulent+Flows%22">Turbulent Flows</a>. , Cambridge University Press. Cambridge. 771(2000).</li>
<li>Most, D., Leipertz, A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Simultaneous+Two-Dimensional+Flow+Velocity+and+Gas+Temperature+Measurements+by+use+of+a+Combined+Particle+Image+Velocimetry+and+Filtered+Rayleigh+Scattering+Technique%5BTitle%5D)+AND+%22Applied+Optics%22%5BJournal%5D)">Simultaneous Two-Dimensional Flow Velocity and Gas Temperature Measurements by use of a Combined Particle Image Velocimetry and Filtered Rayleigh Scattering Technique.</a> Applied Optics. 40 (30), 5379(2001).</li>
<li>Omrane, A., Petersson, P., Aldén, M., Linne, M. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Simultaneous+2D+flow+velocity+and+gas+temperature+measurements+using+thermographic+phosphors%5BTitle%5D)+AND+%22Applied+Physics+B%22%5BJournal%5D)">Simultaneous 2D flow velocity and gas temperature measurements using thermographic phosphors.</a> Applied Physics B. 92 (1), 99-102 (2008).</li>
<li>Someya, S., Li, Y., Ishii, K., Okamoto, K. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Combined+two-dimensional+velocity+and+temperature+measurements+of+natural+convection+using+a+high-speed+camera+and+temperature-sensitive+particles%5BTitle%5D)+AND+%22Experiments+in+Fluids%22%5BJournal%5D)">Combined two-dimensional velocity and temperature measurements of natural convection using a high-speed camera and temperature-sensitive particles.</a> Experiments in Fluids. 50 (1), 65-73 (2010).</li>
<li>Peterson, B., Reuss, D. L., Sick, V. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((High-speed+imaging+analysis+of+misfires+in+a+spray-guided+direct+injection+engine%5BTitle%5D)+AND+%22Proceedings+of+the+Combustion+Institute%22%5BJournal%5D)">High-speed imaging analysis of misfires in a spray-guided direct injection engine.</a> Proceedings of the Combustion Institute. 33 (2), 3089-3096 (2011).</li>
</ol>

Autorzy nie mają nic do ujawnienia.

Podobnie jak w przypadku każdej techniki optycznego pomiaru przepływu, planowanie instalacji szybkiej prędkości obrazu cząstek (PIV) wymaga oceny ograniczeń i oceny najlepszych kompromisów dla danego zadania pomiarowego. Wybór powiększenia obrazu, liczby klatek na sekundę, właściwości arkusza laserowego i algorytmów analizy zależy od szczegółów badanego przepływu. W razie potrzeby należy przeprowadzić pomiary eksploracyjne w celu zidentyfikowania ustawień parametrów dla pomiarów o wysokiej wierności.
W tym artykule opisano ogólną procedurę i niektóre przykładowe wyniki dla szybkich PIV do badania warstwy granicznej przepływu wzdłuż płaskiej płyty. Sekwencja 500 obrazów została zarejestrowana z częstotliwością 5 kHz. Za pomocą mikroskopu dalekosiężnego uzyskano pole widzenia o wymiarach 2,4 x 1,8mm2 zlokalizowane na powierzchni płytki. Wysokiej jakości oświetlenie kropelek oleju z nasion uzyskano za pomocą wiązki z impulsowego lasera na ciele stałym pompowanego diodą, który został rozszerzony na arkusz świetlny za pomocą homogenizatora wiązki. Homogenizator wiązki zawiera układ mikrosoczewek składający się z małych cylindrycznych soczewek i dodatkowych, zintegrowanych teleskopów. Układ mikrosoczewek rozszerza wiązkę kołową w kierunku pionowym, dzieląc przychodzącą wiązkę na wiązki. Następnie następny teleskop nakłada na siebie belki, aby utworzyć arkusz światła o równomiernym rozkładzie natężenia światła w płaszczyźnie arkusza świetlnego prostopadłej do rozchodzenia się wiązki. Obrazy zostały przetworzone przy użyciu algorytmu korelacji krzyżowej PIV. Należy zauważyć, że homogenizowana wiązka jest pomocna, zwłaszcza podczas pracy w pobliżu powierzchni, ale nie jest kluczowa dla opisywanego tutaj zastosowania.
Metoda opisana w tej procedurze umożliwia nieinwazyjne badania przepływów o wysokiej rozdzielczości i dużej szybkości przy użyciu solidnych algorytmów korelacji. Kluczowymi zaletami tej techniki pomiarowej o wysokiej rozdzielczości i szybkości są wysoka rozdzielczość przestrzenna i czasowa oraz możliwość identyfikacji i śledzenia ewolucji struktur w przepływie. Korzystając z tych technik, Alharbi 6 i Jainski et al. 8 z nich wykazało zdolność do wizualizacji i śledzenia struktur wirowych w warstwie granicznej silnika spalinowego. Te kluczowe cechy umożliwiają prowadzenie badań nad strukturą i dynamiką przepływów o wysokich stanach nieustalonych. Co więcej, PIV może być rozszerzony poza dwuwymiarowe, dwuskładnikowe (2D-2C) pola prędkości (jak opisano tutaj) do rozdzielczości 3-składnikowej (3C) w płaszczyźnie (stereo-PIV) i w objętości (tomograficzne PIV, skanujące PIV, holograficzne PIV). Dodatkowo, PIV może być implementowany za pomocą innych technik, takich jak planarna fluorescencja indukowana laserem (PLIF), filtrowane rozpraszanie Rayleigha (FRS) i luminofory termograficzne w celu uzyskania jednoczesnych pomiarów 2D prędkości i innych skalarów (temperatura, stężenie gatunków, współczynniki równoważności) 11-14. Te optyczne, laserowe metody mogą być bezpośrednio stosowane do badania procesów wymiany masy i energii w wielu zastosowaniach, takich jak przepływy przyścienne w silniku spalinowym.

Wielowymiarowe pomiary prędkości i możliwość śledzenia pola przepływu w czasie dostarczają kluczowych informacji w wielu dziedzinach nauki, inżynierii i nauk o zdrowiu. Jedną z najczęściej stosowanych technik obrazowania przepływu jest prędkość obrazowania cząstek (PIV). Początkowo opracowana jako technika planarna, która mierzyła migawki dwóch składowych prędkości w płaszczyźnie, warianty PIV zostały opracowane w celu zapewnienia możliwości pomiaru trójskładnikowego i objętościowego. Wszystkie systemy PIV składają się z cząstek znacznikowych, jednego lub więcej źródeł światła oraz jednej lub więcej kamer. Cząstki stałe lub kropelki są powszechnie stosowane jako cząstki znacznikowe, ale pęcherzyki nieodłącznie związane z przepływem mogą być również używane jako cząstki znacznikowe. Następnie kamera (kamery) obrazuje rozproszone lub wyemitowane światło z cząstek znacznikowych po ich napromieniowaniu przez źródło (źródła) światła. Spośród szerokiej gamy wariantów 1,2 najczęstszą z nich jest uchwycenie dwóch składowych prędkości w płaszczyźnie z szybkością kilku klatek na sekundę. Ostatnio nowe oprzyrządowanie umożliwiło pomiary z dużą liczbą klatek na sekundę (> 1 kHz), które podążają za przepływem w turbulentnych skalach czasowych w zakresie kHz.
PIV określa pole prędkości, śledząc średni ruch grup cząstek z pary obrazów, które są oddzielone znanym opóźnieniem czasowym. Każdy obraz jest podzielony na siatkę regularnie rozmieszczonych okien przesłuchań. Najczęstszy rozmiar okna przesłuchania to 32 x 32 piksele. Algorytm oblicza funkcję korelacji krzyżowej dla wszystkich okien przesłuchań, co daje jeden wektor przemieszczenia na okno zapytania, a zatem tworzy regularną siatkę wektorów. Podzielenie pola wektora przemieszczenia przez opóźnienie czasowe określa następnie pole wektora prędkości.
Podczas planowania pomiarów PIV ważne jest, aby zdać sobie sprawę, że zazwyczaj wybór ustawień eksperymentalnych jest kompromisem między sprzecznymi wymaganiami. Innymi słowy, warunki eksperymentalne muszą być starannie zaplanowane, aby uchwycić aspekty przepływu, które są istotne dla danego badania. Książki Raffela i wsp.1 oraz Adriana i Westerweela 2 dostarczają doskonałych, dogłębnych dyskusji na temat tych ograniczeń. W tym miejscu zwracamy uwagę na kilka, które są najbardziej krytyczne w obecnym kontekście.
Wybór pola widzenia (FOV) ustawi punkt początkowy dla wyboru parametru tutaj. Liczba pikseli na chipie kamery określa następnie rozdzielczość przestrzenną i liczbę wektorów, które są uzyskiwane, przy założeniu, że ktoś zdecyduje się użyć rozmiaru okna przesłuchania 32 x 32 pikseli, często z 50% nałożeniem podczas procedury korelacji krzyżowej. Gęstość wysiewu wynosząca 8-10 cząstek na okno przesłuchania jest ogólnie pożądana, aby wspomóc funkcję korelacji krzyżowej. Istnieją jednak specjalne algorytmy, takie jak prędkość śledzenia cząstek (PTV) i uśrednione w czasie podejścia korelacji, które mogą być stosowane do rozwiązywania sytuacji o niskiej gęstości wysiewu (1-3 cząstki/okno przesłuchania), jak ma to miejsce w przypadku obrazowania w pobliżu powierzchni. Należy zauważyć, że gradienty prędkości w każdym oknie przesłuchania powinny być małe, aby uniknąć odchylenia w wynikowym wektorze reprezentatywnym dla tego okna.
Ustaloną zasadą jest, że przemieszczenia cząstek między pierwszą a drugą klatką nie powinny przekraczać 8 pikseli (1/4 rozmiaru okna przesłuchań), aby zmniejszyć liczbę strat parowania (utrata obrazów cząstek w oknie przesłuchania od pierwszej klatki do drugiej) dla korelacji. W związku z tym czas między dwoma kolejnymi impulsami laserowymi (dt) musi być odpowiednio dostosowany. Jednak zmniejszenie dt poniżej odpowiednika przemieszczeń 8 pikseli zmniejszy zakres dynamiki prędkości, ponieważ dolna granica rozdzielczości jest rzędu 0,1 przemieszczenia piksela.
Podobnie jak w przypadku 8-pikselowego przemieszczenia w płaszczyźnie obrazowania, cząstki o największej prędkości nie powinny przekraczać więcej niż 1/4 grubości arkusza światła, aby ponownie zmniejszyć liczbę strat parowania. Ponieważ opóźnienie czasowe między dwoma impulsami laserowymi jest wykorzystywane do zapewnienia najlepszych korelacji w płaszczyźnie arkusza świetlnego, grubość arkusza jest w tym kontekście zmienną. Chociaż jednorodność natężenia światła nie jest tak krytyczna, jak w przypadku pomiarów opartych na intensywności, takich jak planarne obrazowanie fluorescencyjne indukowane laserem 3, profil wiązki zbliżony do cylindra poprawia jakość PIV, zwłaszcza w przypadku obrazowania o wyższej rozdzielczości.
Ogólnie rzecz biorąc, niektóre założenia dotyczące charakteru badanego przepływu mogą być wykorzystane jako punkt wyjścia przy wyborze parametrów eksperymentalnych. Następnie mogą być konieczne eksperymenty eksploracyjne w celu doprecyzowania ustawień.
Tutaj opisujemy, jak skonfigurować eksperyment PIV, który umożliwia pomiary obrazowania z dużą liczbą klatek na sekundę dwóch składowych prędkości z rozdzielczością przestrzenną, która jest odpowiednia do rozdzielenia struktur warstwy granicznej. Osiąga się to za pomocą lasera na ciele stałego o wysokiej częstotliwości powtarzania TEM 00, pompowanego diodą TEM 00, mikroskopu dalekiego zasięgu i kamery CMOS o dużej liczbie klatek na sekundę. Uwzględniono również kilka szczegółów dotyczących obrazowania w pobliżu powierzchni.

<table><tbody><tr><td>Szybki laser Nd: YAG 532 nm </td><td>Quantronix</td><td>Model: Mikroskop</td> dalekiego zasięgu Hawk I <td></td></tr><tr><td>(QM-100)</td> <td>Questar</td> Model: Szybka kamera CMOS <td>QM-100</td><td></td></tr><tr><td>(Phantom v7.3)</td> <td>Model badań wizyjnych</td><td>: Atomizer Phantom v7.3</td><td></td></tr><tr><td>(TSI 9306)</td> Model <td>TSI</td><td>: 9306</td><td></td></tr><tr><td>Olej silikonowy </td><td>Dow Corning CST 510</td><td>CST 510 Płyn</td><td></td></tr><tr><td>Homogenizator wiązki </td><td>Fraunhofer</td><td></td><td>Część wykonana na zamówienie</td></tr><tr><td>45 & stopnie wysokiego odbicia (HR) Lustro obrotowe 532 nm </td><td>Laser Optik</td><td></td><td>Wielu dostawców</td></tr><tr><td>Przysłona</td><td></td><td></td><td>Wielu dostawców</td></tr><tr><td>Cel</td><td></td><td></td> kalibracji Oprogramowanie do <td></td></tr><tr><td>nagrywania i przetwarzania PIV </td> na zamówienie Oprogramowanie <td>LaVision</td><td>: Da Szybki kontroler Vis</td><td></td></tr><tr><td>(HSC)</td><td>LaVision</td><td></td><td></td></tr><tr><td>Szyna optyczna i nośniki </td><td></td><td></td><td>Wielu dostawców</td></tr><tr><td>Bloki i pułapki wiązki laserowe </td><td></td><td></td><td>Wielu dostawców</td></tr><tr><td>Mocowania do elementów </td><td></td><td></td> optycznych <td>Wielu dostawców</td></tr><tr><td>Etap tłumaczenia </td><td>Newport</td><td></td><td></td></tr><tr><td>Rurki metalowe do stworzenia przepływ strumieniowy</td><td>McMaster-Carr</td><td></td><td>Wielu dostawców</td></tr><tr><td>Kwadrat kombinowany i kwadrat</td><td></td><td></td> centrujący <td>Wielu dostawców</td></tr></tbody></table>

szybkie obrazowanie cząstek w pobliżu powierzchni

High-speed Particle Image Velocimetry Near Surfaces

Wielowymiarowe i przejściowe przepływy odgrywają kluczową rolę w wielu dziedzinach nauki, inżynierii i nauk o zdrowiu, ale często nie są dobrze rozumiane. Złożona natura tych przepływów może być badana za pomocą prędkości obrazowania cząstek (PIV), techniki obrazowania laserowego dla optycznie dostępnych przepływów. Chociaż istnieje wiele form PIV, które rozszerzają tę technikę poza pierwotne możliwości planarnego dwuskładnikowego pomiaru prędkości, podstawowy system PIV składa się ze źródła światła (lasera), kamery, cząstek znacznikowych i algorytmów analizy. Parametry obrazowania i nagrywania, źródło światła i algorytmy są dostosowywane w celu optymalizacji zapisu pod kątem przepływu zainteresowania i uzyskania prawidłowych danych dotyczących prędkości. 
Powszechne badania PIV mierzą dwuskładnikowe prędkości w samolocie z prędkością kilku klatek na sekundę. Jednak ostatnie osiągnięcia w oprzyrządowaniu ułatwiły pomiary z dużą liczbą klatek na sekundę (> 1 kHz), które są w stanie rozwiązać przepływy przejściowe z wysoką rozdzielczością czasową. W związku z tym pomiary z dużą liczbą klatek na sekundę umożliwiły badania nad ewolucją struktury i dynamiki przepływów o wysokich stanach nieustalonych. Badania te odgrywają kluczową rolę w zrozumieniu podstawowej fizyki złożonych przepływów.
Szczegółowy opis wykonywania szybkiego planarnego PIV o wysokiej rozdzielczości w celu zbadania przejściowego przepływu w pobliżu powierzchni płaskiej płyty jest przedstawiony tutaj. Dołączono szczegółowe informacje dotyczące dostosowywania ograniczeń parametrów, takich jak właściwości obrazu i zapisu, właściwości arkusza laserowego oraz algorytmy przetwarzania w celu dostosowania PIV do dowolnego przepływu zainteresowania.

Zdjęcie układu przedstawiono na Rysunku 1. Surowe obrazy cząstek w oknie analizy o rozmiarze 32 x 32 pikseli w pobliżu ścianki z dwóch kolejno zarejestrowanych obrazów przedstawiono na Rysunku 2. Cząstki na Rysunku 2a są przesunięte o 2-3 piksele w prawo na Rysunku 2b i spełniają „zasadę jednej czwartej”, która mówi, że przesunięcia cząstek w płaszczyźnie i poza nią nie powinny przekraczać ¼ rozmiaru okna analizy. Ponadto gęstość cząstek w oknie analizy powinna wynosić około 8-10 cząstek, ponieważ algorytmy korelacji PIV śledzą grupy cząstek. Jednak gęstość zasiać w badaniach PIV w pobliżu ścianki często wynosi od 1 do 3 cząstek. W związku z tym w badaniach z niższą gęstością zasiać należy stosować specjalne algorytmy, takie jak algorytmy wektoremetrii śledzenia cząstek (PTV), które śledzą poszczególne cząstki 1,2,4-6. Do rozwiązania problemów z niską gęstością zasiać można również zastosować podejście korelacji uśrednionej w czasie 7,8, co jednak zazwyczaj prowadzi do utraty rozdzielczości czasowej. Dodatkowo obrazowanie w pobliżu ścianek jest utrudnione przez jasne odbicia laserowe, które mogą negatywnie wpływać na korelacje PIV i generować błędne wektory. Jasne odbicia te ograniczają również położenie pierwszego poprawnego wektora prędkości w kierunku normalnym do ścianki. Przetwarzanie wstępne surowych obrazów cząstek jest niezbędne w celu zmniejszenia wpływu szumów tła pochodzących z takich źródeł jak odbicia laserowe. W tej demonstracji pierwszy poprawny wektor znajdował się w odległości 23 μm od ścianki.
Po przetworzeniu surowych obrazów cząstek za pomocą algorytmów korelacji PIV, należy ocenić jakość i poprawność wynikających z nich pól wektorów prędkości. W surowych polach wektorowych nie da się uniknąć wektorów błędnych, jednak posiadają one kilka charakterystycznych cech. Nieprawidłowe wektory często występują w pobliżu powierzchni, na krawędziach arkusza światła oraz na krawędziach przepływu. Ponadto wartość i kierunek nieprawidłowych wektorów znacząco różnią się od wektorów sąsiednich i nie mają uzasadnienia fizycznego. W przypadku tego przykładu przepływu w warstwie przyściennej, poprawne wektory prędkości powinny być skierowane od lewej do prawej, co wskazują przesunięcia cząstek na Ryc. 2. Dodatkowo prędkości powinny maleć w pobliżu ścianki ze względu na warunek braku poślizgu 9. Przedstawione na Ryc. 3 chwilowe pola prędkości spełniają oba te kryteria fizyczne. Inną przydatną miarą oceny poprawności wyników PIV jest określenie wyboru wektora (vector choice) dla każdego wektora w polu wektorów prędkości. Zasadniczo pole wektorowe powinno składać się w &gt;= 95% z wektorów pierwszego wyboru, tzn. takich, które nie wymagały postprocessingu, aby można było zastosować solidne algorytmy postprocessingu do wykrywania i zastępowania błędnych wektorów bez generowania znaczących artefaktów 2. Chwilowe pola wektorowe pokazane na Ryc. 3 składają się w całości z wektorów pierwszego wyboru.
Znaczenie pomiarów PIV o wysokiej prędkości, czyli kinematograficznych, staje się oczywiste po analizie sekwencji obrazów przepływu w czasie. Pola wektorowe prędkości chwilowej (Vi) oraz fluktuacji prędkości (V') na początku, w środku i na końcu sekwencji nagrań przedstawiono na Rysunku 3. Korzystając z dekompozycji Reynoldsa, Vi stanowi sumę uśrednionego pola prędkości (<img alt='Symbol prędkości „v” w równaniu fizycznym; ilustracja edukacyjna.' fo:src="/files/ftp_upload/50559/50559eq2highres.jpg" src="/files/ftp_upload/50559/50559eq2.jpg" data-altupdatedat="1787630408000" data-alt="Równanie">) oraz V' 10. W tym doświadczeniu <img alt='Symbol prędkości „v” w równaniu fizycznym; ilustracja edukacyjna.' fo:src="/files/ftp_upload/50559/50559eq2highres.jpg" src="/files/ftp_upload/50559/50559eq2.jpg" data-altupdatedat="1787630408000" data-alt="Równanie"> wyznaczono poprzez uśrednienie w czasie wszystkich obrazów w sekwencji. Chwilowe pola wektorowe w całej sekwencji nagrań są do siebie bardzo podobne i pokazują przepływ poruszający się od lewej do prawej strony. Wyniki te wskazują również, że przepływ odbywa się głównie w kierunku poziomym, ponieważ pozioma składowa prędkości (u) jest znacznie większa od pionowej składowej prędkości (v). Pola wektorowe fluktuacji wskazują również, że poziome fluktuacje prędkości (u') są większe od pionowych fluktuacji prędkości (V'). Jednakże fluktuacje wskazują także, że przepływ zwalnia, ponieważ u' zmienia swój kierunek w trakcie sekwencji nagrań.
Uśrednione w czasie oraz chwilowe profile u w kilku różnych momentach sekwencji pomiarowej przedstawiono na rysunku 4, co potwierdza, że prędkość przepływu maleje w czasie. Profile u wyznaczono poprzez uśrednienie czterech sąsiednich kolumn wektorów w celu zwiększenia istotności statystycznej wyników w pobliżu ścianki. Procedura ta była stosowana w poprzednich pracach 6,8. Słupki błędów wskazują dwukrotność odchylenia standardowego czterech sąsiednich kolumn wektorów. Największy słupek błędu występuje w pobliżu powierzchni płytki, co potwierdza trudność w stosowaniu algorytmów korelacji PIV w obszarach o niskiej gęstości cząstek śledzących. Kilka algorytmów analizy zostało zaprojektowanych w celu rozwiązania problemu niskiej gęstości cząstek śledzących, takich jak PTV 5,6 oraz metody korelacji uśrednionej w czasie 7,8.
<img alt="Układ optyczny lasera Nd:YAG z kamerą CMOS, mikroskopem, homogenizatorem wiązki i lustrem 532 nm." src="/files/ftp_upload/50559/50559fig1.jpg" data-altupdatedat="1787630408000" data-alt="Rysunek 1"> 
Rysunek 1. Montaż laboratoryjny.
<img alt="Obrazy absorpcji przejściowej przedstawiające dane rozdzielcze w czasie przy 0,2 ms i 0,4 ms; odległość w funkcji czasu." fo:content-width="5.5in" fo:src="/files/ftp_upload/50559/50559fig2highres.jpg" src="/files/ftp_upload/50559/50559fig2.jpg" data-altupdatedat="1787630408000" data-alt="Rycina 2"> 
Rycina 2. Obrazy cząstek w obszarze analizy o rozmiarze 32 x 32 pikseli w pobliżu ścianki przy a) t = 0,2 msec oraz b) t = 0,4 msec. Wymiary fizyczne okna analizy wynoszą 96 x 96 μm2.
<img alt="Zachowanie pola wektorów prędkości; wykres; dynamika płynów; przedstawia analizę zmienności czasowej i odległości powierzchniowej." src="/files/ftp_upload/50559/50559fig3.jpg" data-altupdatedat="1787630408000" data-alt="Rycina 3"> 
Rycina 3. Po lewej: chwilowe (Vi), a po prawej: fluktuacyjne (V') pola prędkości na początku, w środku i na końcu sekwencji rejestracji. Pola wektorowe składają się wyłącznie z wektorów pierwszego wyboru. Dla przejrzystości przedstawiono mniejszy podzbiór pól wektorowych. Pola Vi wskazują przepływ z lewej do prawej, podczas gdy V' mają kierunek przeciwny. Należy zauważyć, że dla przejrzystości wyświetlana jest tylko co czwarta kolumna wektorów w kierunku poziomym. Dodatkowo, skala prędkości między polami Vi a V' jest różna, co zaznaczono w lewym górnym rogu każdego obrazu.
<img alt="Wykres dynamiki płynów przedstawiający profil prędkości; pomiary uśrednione w czasie w porównaniu z chwilowymi." src="/files/ftp_upload/50559/50559fig4.jpg" data-altupdatedat="1787630408000" data-alt="Rysunek 4"> 
Rysunek 4. Profile prędkości poziomej (u) w różnych momentach podczas przepływu. Uśredniony w czasie profil u przedstawiono za pomocą kółek. Słupki błędów widoczne na profilu dla t = 0.1 msec są reprezentatywne dla słupków błędów we wszystkich pozostałych momentach. Historia czasowa profili u wykazuje spadek przepływu w czasie.

Procedura badania przepływów przejściowych w pobliżu granic za pomocą wysokiej rozdzielczości, szybkiej prędkości obrazowania cząstek (PIV) jest opisana tutaj. PIV to nieinwazyjna technika pomiarowa mająca zastosowanie do każdego optycznie dostępnego przepływu poprzez optymalizację kilku ograniczeń parametrów, takich jak właściwości obrazu i zapisu, właściwości arkusza laserowego i algorytmy analizy.

Szybkie obrazowanie cząstek w pobliżu powierzchni

Multi-dimensional and transient flows play a key role in many areas of science, engineering, and health sciences but are often not well understood. The ...

Ogólnym celem tego filmu jest zademonstrowanie sposobu konfiguracji i przeprowadzania eksperymentu z wykorzystaniem laboratoryjnej welocymetrii obrazowej cząstek (PIV). W pierwszej kolejności należy oszacować parametry rejestracji niezbędne do rozdzielenia prędkości i gradientów oczekiwanych w przepływie. Drugim krokiem jest dostosowanie gęstości rozsiewu, aby uzyskać od 8 do 10 cząstek na okno przesłuchowe.Następnie należy dostosować częstotliwość klatek kamery oraz czas między impulsami lasera. Aby uzyskać przesunięcie cząsteczek między parami obrazów mniejsze niż jedna czwarta rozmiaru okna interogacyjnego, ostatnim krokiem jest zapisanie danych i przetworzenie wyników za pomocą odpowiednich algorytmów przetwarzania. Ostatecznie technika PIV może zostać wykorzystana do badania struktury i dynamiki każdego przepływu z wystarczającym dostępem optycznym. Głównymi zaletami tej techniki w porównaniu z istniejącymi metodami pomiaru przepływu, takimi jak anemometria drutowa czy rurki Pitota, jest fakt, że jest ona nieinwazyjna i wielowymiarowa.Ta metoda może pomóc w odpowiedzi na kluczowe pytania dotyczące pola prędkości przepływu w wielu zastosowaniach, takich jak struktura i dynamika przejściowego przepływu warstwy przyściennej powstającej w silniku spalinowym. W eksperymentach tych wykorzystuje się lasery, dlatego przed przystąpieniem do pracy należy zapoznać się z zasadami bezpieczeństwa pracy z laserami i zapewnić sobie odpowiedni sprzęt ochronny. Każdy eksperymentator powinien nosić odpowiednie okulary ochronne.Zawieś zasłony ochronne przed laserem. Aby odizolować eksperyment laserowy, użyj lampy ostrzegawczej, która powiadomi innych o działaniu lasera. W pierwszym kroku dobierz odpowiedni obiektyw kamery – w tej konfiguracji należy zastosować obiektyw mikroskopowy do dużych odległości.Dodatkowo należy wyznaczyć częstotliwość klatek kamery oraz opóźnienie czasowe w oparciu o praktyczne wytyczne dla symetrii obrazu cząstek. W tym filmie badany będzie przepływ w pobliżu powierzchni płyty. Częstotliwość klatek i opóźnienie czasowe oparto na szacunkowych prędkościach przepływu w obszarze przyściennym.Należy upewnić się, że laser jest wypoziomowany, a optyka formująca wiązkę jest na swoim miejscu. W niniejszej demonstracji światło ulega rozproszeniu na wytworzonych kropelkach oleju silikonowego. Używając atomizera do oleju, należy uruchomić przepływ i dostosować gęstość osiadania za pomocą dysz atomizera oraz zaworów obejściowych w urządzeniu.Rozpocznij optymalizację, upewniając się, że szybki kontroler oraz oprogramowanie są odpowiednio skonfigurowane. Ustaw kamerę w tryb ciągłego nagrywania. Włącz laser oraz atomizer.Ustaw ostrość mikroskopu dalekiego zasięgu, aby uzyskać wyraźny obraz cząsteczek. Upewnij się również, że intensywność obrazów cząsteczek nie powoduje nasycenia kamery. Po uzyskaniu ostrych obrazów cząsteczek wyłącz tryb ciągłego pobierania danych.Po zakończeniu procesu zarejestruj kilkaset obrazów przepływu, a następnie sprawdź, czy obrazy są odpowiednio ustawione pod kątem ostrości, rozmieszczenia cząstek oraz czy dryf cząstek między obrazami wynosi około ośmiu pikseli. Poniżej znajduje się przykład rzadkiego zasiewu i dryfu cząstek wynoszącego około jednego piksela w oknie o rozmiarze 32 na 32 piksele, co służy do korekty gęstości zasiewu cząstek. Zwiększ lub zmniejsz liczbę dysz atomizera, aby dostosować wartość dryfu między obrazami.Zwiększyć lub zmniejszyć odstęp czasu między impulsami lasera. Niniejszy przykład jest bliski ideałowi, z ośmioma do 10 cząsteczek na okno analizy o rozmiarze 32 na 32 piksele, a przesunięcie cząsteczek wynosi mniej niż osiem pikseli. Jeśli śledzenie grup cząsteczek w serii obrazów jest utrudnione, tak jak w przypadku tych zdjęć, może to oznaczać zbyt duży ruch poza płaszczyzną obrazu.Skoryguj przesunięcie kamery względem punktu ogniskowania lub zwiększ jej odległość roboczą od arkusza światła, aby zrekompensować zbyt duży ruch poza płaszczyzną. Na koniec, po ponownym wprowadzeniu wszystkich korekt, zarejestruj kilkaset obrazów i sprawdź cząstkę. Dostosuj gęstość nasycenia i ostrość.Powtarzaj regulacje do momentu spełnienia kryteriów. Przed każdym uruchomieniem eksperymentu przeprowadź kalibrację intensywności kamery z założoną osłoną zespołu kamery, aby ustawić punkt odniesienia dla intensywności. Po zakończeniu kalibracji zdejmij osłonę.Ustaw laser na częstotliwość repetycji i prąd wyjściowy określone jako optymalne. Przed przełączeniem lasera w tryb zewnętrzny upewnij się, że laser otrzymuje ciągły sygnał wyzwalający zgodny z ustawioną częstotliwością. Włącz tryb rejestracji lasera.Sekwencja obrazów tła przedstawiających jedynie arkusz światła muskający powierzchnię płytki. Nie występuje istotne rozpraszanie tła, w związku z czym większość obrazów jest ciemna. Należy zapisać te obrazy do późniejszej obróbki.Następnie uruchom przepływ i poczekaj na jego stabilizację, ustawiając kamerę w tryb ciągłej akwizycji obrazów i upewniając się, że rejestrowane są ostre obrazy cząstek. Po wykonaniu tej czynności wyłącz tryb ciągłej akwizycji. Dla każdego pomiaru wprowadź żądaną liczbę obrazów do zarejestrowania.Następnie naciśnij przycisk nagrywania. Po zakończeniu nagrywania wyłącz przepływ i laser. Przejrzyj sekwencję obrazów, aby sprawdzić dryft cząstek, gęstość ich osadzenia oraz ostrość obrazu cząstek.Jeśli wynik jest satysfakcjonujący, zapisz nagranie. Wykonaj tyle serii, ile jest koniecznych, aby zagwarantować prawidłową konfigurację lub uzyskać odpowiednie dane statystyczne. Przy każdej zmianie zestawu kamery lub ostrości wykonaj zdjęcia, które pomogą określić położenie powierzchni w późniejszej analizie.W pierwszej kolejności należy zwiększyć czas ekspozycji kamery. Następnie należy umieścić wzorzec kalibracyjny w płaszczyźnie arkusza światła i upewnić się, że styka się on z płytką. Należy oświetlić wzorzec od tyłu za pomocą źródła światła, podczas gdy kamera stale rejestruje obrazy.Dostosuj cel tak, aby zarejestrowany obraz był ostry i niedodkształcony. Upewnij się, że punkt styku płytki z celem jest widoczny na obrazie. Jest to kluczowe dla określenia położenia płytki w zapisie obrazu. Zarejestruj 10 obrazów celu kalibracyjnego. Przetwarzanie danych odbywa się w oprogramowaniu.Najpierw uśrednij każdy zestaw 10 obrazów kalibracyjnych. W oknie dialogowym kalibracji wykorzystaj uzyskany obraz średni dla obrazów docelowych. Pozwoli to wyznaczyć rzeczywiste wymiary świata obrazów z akwizycji.Zastosuj każdy obraz kalibracyjny do odpowiadającego mu zestawu obrazów. Wyznacz położenie płytki na obrazach skalibrowanych. Następnie wyznacz średnią z obrazów tła.Porównaj wartości intensywności obrazu tła średniego z wartościami dla osiadłych cząsteczek, aby ustalić, czy odbicie lasera od powierzchni w istotny sposób wpływa na szum tła. W tym doświadczeniu odbicia nie miały istotnego wpływu na tło. Przetwórz wstępnie skalibrowane obrazy przepływu, stosując filtr górnoprzepustowy w celu usunięcia dużych fluktuacji intensywności w tle; w tym przypadku nie wpłynie to na sygnały z cząsteczek.Aby zdefiniować obszar w oprogramowaniu, w którym będzie przeprowadzana analiza asymetrii obrazu cząstek, należy zastosować maskę geometryczną. Następnie należy wybrać procedurę obliczania wektorów. Wybrana opcja oprogramowania wykonuje dwa wstępne przejścia z kwadratowymi oknami analizy o wymiarach 64 na 64 piksele i 50% zakładce.Następnie wykonuje się trzy przejścia, wykorzystując w oprogramowaniu kwadratowe okna przesłuchujące o rozmiarze 32 na 32 piksele z 50% zakładką. Należy zastosować podprogramy w celu usunięcia wektorów prędkości, które nie spełniają kryteriów ustanowionych w celu zapewnienia jakości wyników korelacji wzajemnej. Te przykładowe dane przedstawiają na czerwono pole wektorów prędkości w pobliżu powierzchni płytki jako funkcję odległości wzdłuż powierzchni oraz wysokości nad powierzchnią.Bezpośrednio po obliczeniach. Najpierw zastosuj podprogram w celu ponownego zdefiniowania obszaru analizy; podprogramy te najpierw określają stosunek szczytów, który porównuje najwyższe szczyty korelacji ze wspólną korelacją. Wektory tła ze stosunkiem szczytów bliskim jedności wskazują na fałszywy szczyt i są usuwane, co odzwierciedlają białe pola oraz brakujące strzałki.Następnie algorytmy odrzucają wektory, których składowa znajduje się poza zakresem jednego odchylenia standardowego (root mean square) od mediany dla danej składowej, obliczonej na podstawie podzbioru wektorów; algorytmy eliminują również wektory pozorne, które powstają w wyniku dużego nakładania się podczas obliczania wektora prędkości. W ostatnim kroku algorytmy wypełniają puste przestrzenie w polu, interpolując wartości pomiędzy sąsiednimi wektorami niezerowymi. Wektory te są przedstawione kolorem purpurowym.W ostatnim kroku oceniono jakość wyników; zmierzono przepływ w pobliżu powierzchni płaskiej płyty, a dane dotyczące pól prędkości chwilowej i prędkości fluktuacji przedstawiono na tym wykresie. Po lewej stronie pokazano podzbiór pola prędkości chwilowej w czasie T równym 0,1 ms, a po prawej pole prędkości fluktuacji. Wektory skali znajdują się w lewym górnym rogu każdego wykresu. W późniejszym czasie, w kroku w połowie sekwencji pomiarowej, pole prędkości chwilowej po lewej stronie pozostaje stabilne, ale pole fluktuacji uległo gwałtownej zmianie.Wartość wektorów uległa zmniejszeniu, a w kierunku pod koniec sekwencji nagrań występuje mniejsza jednorodność. Pole prędkości chwilowej nadal wskazuje ogólny kierunek przepływu z lewej do prawej strony. Pole fluktuacji zmieniło kierunek na przeciwny.Na tym rysunku przedstawiono profil prędkości poziomej w różnych odstępach czasu podczas przepływu. Słupki błędów widoczne na profilu dla T równym 0,1 ms są reprezentatywne dla słupków błędów w pozostałych punktach czasowych. Należy zauważyć, że historia czasowa wykazuje spadek przepływu w czasie.Profil uśredniony w czasie przedstawiono za pomocą kółek. Po obejrzeniu tego filmu powinni Państwo mieć dobrą wiedzę na temat tego, jak dostosować i zoptymalizować parametry symetrii obrazowania cząstek, aby przystosować procedurę do dowolnych badań przepływu. Należy pamiętać, że praca z laserami może być niezwykle niebezpieczna, dlatego podczas wykonywania tej procedury należy zawsze przestrzegać środków ostrożności, takich jak noszenie okularów ochronnych i używanie znaków ostrzegawczych w celu powiadomienia innych osób.

high-speed-particle-image-velocimetry-near-surfaces

Research

JoVE Journal

Engineering

Inżynieria

Przekładanie prędkości obrazu cząstek o dużej szybkości w pobliżu powierzchni

Niniejszy protokół zawiera instrukcje dotyczące obsługi autonomicznego urządzenia do pomiaru prędkości przepływu pod wodą (SCUVA), które zostało zaprojektowane do kwantyfikacji przepływów generowanych przez zwierzęta in situ. Ponadto protokół ten odnosi się do wyzwań wynikających z warunków terenowych, w tym do ruchów operatora, przewidywania pozycji zwierząt oraz orientacji urządzenia SCUVA.

Ilościowy pomiar przepływów in situ przy użyciu autonomicznego urządzenia do podwodnej welocymetrii (SCUVA)

Opisano system echoskopowej welocymetrii obrazowej cząstek (EPIV), który umożliwia pozyskiwanie dwuwymiarowych pól prędkości w płynach nieprzezroczystych optycznie lub przez geometrie nieprzezroczyste optycznie, oraz przedstawiono pomiary walidacyjne w przepływie w rurze.

Echo Particle Image Velocimetry

Przedstawiono technikę wykonywania ilościowego obrazowania trójwymiarowego (3D) dla różnych rodzajów przepływów płynów. Wykorzystując koncepcje z zakresu obrazowania pola światła (Light Field Imaging), rekonstruujemy objętości 3D z zestawów obrazów. Nasze wyniki 3D obejmują szeroki zakres danych, w tym pola prędkości oraz wielofazowe rozkłady wielkości pęcherzyków.

Wyznaczanie trójwymiarowych pól przepływu za pomocą wielokamerowego obrazowania pola światła

Prędkość obrazowania mikrocząstek (μPIV) służy do wizualizacji sparowanych obrazów mikrocząstek zasianych w przepływie krwi, które są skorelowane krzyżowo, aby uzyskać dokładny profil prędkości.  Na podstawie tych pomiarów można uzyskać szybkość ścinania, prędkość maksymalną, kształt profilu prędkości i natężenie przepływu, z których każdy ma zastosowania kliniczne.

Obrazowa prędkość mikrocząstek do pomiarów profilu prędkości mikroprzepływów krwi

Opisano dwa eksperymentalne urządzenia do badania uderzenia strumienia cieczy w poruszającą się powierzchnię o dużej prędkości: urządzenie z armatą powietrzną i urządzenie z wirującym dyskiem. Aparatura służy do określania optymalnych podejść do stosowania płynnego modyfikatora tarcia (LFM) na torach kolejowych w celu kontroli tarcia na szczycie szyny.

Wizualizacja zaparcia strumienia cieczy o dużej prędkości na ruchomej powierzchni

Trójwymiarowy system prędkości śledzenia cząstek (3D-PTV) oparty na szybkiej kamerze z czterostopniowym rozdzielaczem jest opisany tutaj. Technika ta jest stosowana do strumienia strumienia z okrągłej rury w pobliżu dziesięciu średnic w dół rzeki pod numerem Reynoldsa Re ≈ 7000.

Trójwymiarowa prędkość śledzenia cząstek do zastosowań turbulencyjnych: przypadek przepływu strumieniowego

Prezentowane są dwie różne metody charakteryzowania początkowego ruchu cząstek pojedynczego koralika jako funkcji geometrii złoża osadu od przepływu laminarnego do turbulentnego.

Wizualna charakterystyka początkowego ruchu cząstek w regularnych podłożach: od warunków laminarnych do turbulentnych

Technika opisana tutaj oferuje tanią i stosunkowo prostą metodę jednoczesnego pomiaru kinematyki cząstek i turbulencji w przepływach o niskich stężeniach cząstek. Turbulencje są mierzone za pomocą prędkości obrazu cząstek (PIV), a kinematyka cząstek jest obliczana na podstawie obrazów uzyskanych za pomocą szybkiej kamery w nakładającym się polu widzenia.

Jednoczesny pomiar turbulencji i kinematyki cząstek z wykorzystaniem technik obrazowania przepływowego

Tutaj prezentujemy uproszczone metody wytwarzania przezroczystych fantomów nerwowo-naczyniowych i charakteryzowania przepływu w nich. Zwracamy uwagę na kilka ważnych parametrów i pokazujemy ich związek z dokładnością pola.

Badania szybkościomierzy obrazu cząstek w mezoskali, przepływów nerwowo-naczyniowych in vitro

Obecny protokół opisuje pomiary prędkości obrazu cząstek (PIV) wykonywane w celu zbadania przepływu zatokowego przez układ in vitro przezcewnikowej zastawki aortalnej (TAV). Określane są również parametry hemodynamiczne na podstawie prędkości.

Badanie hemodynamiki za pomocą wizji obrazu cząstek za pomocą fantomu aorty

Poznaj metodę PIV do analizy przepływów płynów w pobliżu powierzchni przy użyciu obrazowania cząstek. Technika ta umożliwia nieinwazyjne pomiary wielowymiarowych pól prędkości z wysoką rozdzielczością czasową i przestrzenną.

Szybkie obrazowanie cząstek w pobliżu powierzchni

Przeglądaj więcej filmów