Name: Ewolucja konstrukcji klatek schodowych w konwekcji dyfuzyjnej (Wideo) | JoVE Journal
Uploaded: 2018-09-05T15:30:00.000+00:00
Duration: 7 min 28 s
Description: Poznaj metodę symulacji ewolucji struktur klatek schodowych w konwekcji dyfuzyjnej w układzie laboratoryjnym z płytami miedzianymi i termistorami. To podejście umożliwia obserwację procesów niemożliwych do bezpośredniego zobaczenia w oceanie.

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

Ta praca była wspierana przez chińskie granty NSF (41706033, 91752108 i 41476167), granty Grangdong NSF (2017A030313242 i 2016A030311042) oraz granty LTO (LTOZZ1801).

1. Zbiornik roboczy
Uwaga: Eksperyment przeprowadzono w zbiorniku o przekroju prostokątnym. Zbiornik składa się z płyt górnej i dolnej oraz ścianki bocznej. Płyty górna i dolna wykonane są z miedzi z powierzchniami galwanicznymi. W płycie górnej znajduje się komora wodna. W płycie dolnej umieszczono elektryczną matę grzewczą. Ścianka boczna wykonana jest z przezroczystego pleksiglasu. Wymiary zbiornika to Lx = 257 mm (długość), Ly = 65 mm (szerokość) oraz Lz = 257 mm (wysokość). Grubość ścianki bocznej wynosi 9,5 mm.
<ol>
	<li>Starannie wyczyść miedziane płyty oraz pleksiglasowe ścianki boczne wodą destylowaną.</li>
	<li>Zmontuj zbiornik za pomocą śrub, aby zapewnić jego szczelność.</li>
	<li>Ustaw stalową ramę wsporczą (wysokość 150 mm) na stole optycznym i zamocuj zbiornik nad ramą, przekładając pomiędzy nimi płytę termoizolacyjną, która ogranicza ucieczkę ciepła ze zbiornika roboczego do stołu.</li>
	<li>W każdej płycie umieść trzy termistory (stabilność temperatury 0,01 °C) i podłącz je do cyfrowego multimetru. Należy pamiętać, że termistory te służą do monitorowania temperatury płyt górnej i dolnej.</li>
	<li>Umieść mikropróbnik przewodności i temperatury (MSCTI) wewnątrz zbiornika i podłącz go do wielofunkcyjnego systemu akwizycji danych (MDA). Zamocuj MSCTI do zmotoryzowanego precyzyjnego stołu przesuwnego (MPTS). 
	Uwaga: MSCTI może być przesuwany w górę i w dół za pomocą ruchu pionowego, co pozwala na uzyskanie profili temperatury i zasolenia czynnika roboczego. W tym przypadku MSCTI posiada stabilność temperatury 0,01 °C i stabilność zasolenia 1%. MPTS posiada dokładność pozycjonowania 0,05 mm.</li>
	<li>W odpowiednich programach sterujących cyfrowym multimetrem oraz wielofunkcyjnym systemem akwizycji danych ustaw parametry, takie jak częstotliwości próbkowania, kanały akwizycji danych i ścieżki zapisu. Tutaj ustaw częstotliwości próbkowania cyfrowego multimetru i systemu akwizycji danych odpowiednio na 1,0 i 128 Hz.</li>
	<li>W programowaniu sterującym MPTS ustaw parametry ruchu MSCTI, w tym pozycję początkową, pozycje najniższą i najwyższą, prędkość przesuwania oraz przyspieszenie. Tutaj ustaw prędkość przesuwania na 1 mm/s, przyspieszenie na 0,5 mm/s2, a pozycje najniższą i najwyższą na odpowiednio 20 i 220 mm powyżej płyty dolnej. Skutkuje to okresem pracy MPTS wynoszącym 404 s dla jednego pomiaru w cyklu góra-dół. Ustaw pozycję początkową MSCTI w pozycji najniższej.</li>
	<li>Utrzymuj temperaturę w pomieszczeniu na stałym poziomie około 24 °C za pomocą dwóch mocnych klimatyzatorów (moc robocza 30 W).</li>
</ol>
2. Aparatura optyczna
Uwaga: W trakcie eksperymentu ewolucja schodka DC będzie monitorowana przy użyciu techniki shadowgraph, realizowanej zgodnie z poniższymi procedurami
<ol>
	<li>Przymocuj arkusz papieru kalkowego (25,7 cm x 25,7 cm) do zewnętrznej ścianki zbiornika.</li>
	<li>W jako źródła światła zastosuj lampę LED z wąską wiązką. Umieść źródło światła w odległości około 5 m od przeciwległej strony zbiornika, aby uzyskać niemal kolimowane światło. Należy zauważyć, że podczas eksperymentu struktura warstwowa płynu DC jest widoczna na papierze kalkowym ze względu na zmianę gęstości płynu (odpowiadającą zmianie współczynnika załamania światła).</li>
	<li>Ustaw kamerę wysokich prędkości po tej samej stronie co papier kalkowy. Kamera powinna znajdować się w odległości około 1 m od zbiornika, aby można było zarejestrować struktury warstwowe w pełnej skali zbiornika.</li>
	<li>Ustaw częstotliwość próbkowania kamery. Należy pamiętać, że częstotliwość próbkowania powinna być odpowiednia do uchwycenia szczegółów ewolucji schodkowej. W tym przypadku częstotliwość próbkowania kamery wynosi 25 Hz.</li>
	<li>Włącz lampę oraz kamerę i lekko skoryguj ich położenie oraz odległość, aby upewnić się, że kamera rejestruje wyraźne obrazy.</li>
</ol>
3. Ciecz robocza
<ol>
	<li>Przygotować wodę słoną i świeżą w dwóch zbiornikach.
	<ol>
		<li>Połącz dwa identyczne zbiorniki prostokątne (zbiornik A i zbiornik B) za pomocą elastycznej rurki (o długości 10 cm, średnicy wewnętrznej 6 mm i średnicy zewnętrznej 10 mm), mocując ją do dna każdego ze zbiorników.</li>
		<li>Napełnij zbiornik A wodą zasoloną, której stężenie masowe soli (tj. zasolenie) wynosi w tym przykładzie 60 g/kg.</li>
		<li>Napełnić zbiornik B równą objętością odgazowanej świeżej wody i użyć elektrycznego mieszadła magnetycznego do ciągłego homogenizowania cieczy.</li>
		<li>Utrzymaj początkową temperaturę cieczy w obu zbiornikach na poziomie temperatury pokojowej (24 °C).</li>
	</ol>
	</li>
	<li>Ustanów liniową stratyfikację gęstości w zbiorniku roboczym.
	<ol>
		<li>Zastosować metodę dwuzbiornikową29 w celu ustanowienia początkowej liniowej stratyfikacji słonej wody w zbiorniku roboczym.</li>
		<li>Ustaw zbiornik A i B na tej samej wysokości, która powinna być o 30 cm wyższa niż zbiornik roboczy. Połącz zbiornik B i zbiornik roboczy za pomocą innego elastycznego wężyka (o długości 50 cm, średnicy wewnętrznej 2 mm i średnicy zewnętrznej 5 mm) w ich dolnych częściach. Dzięki różnicy ciśnienia cieczy w tych dwóch zbiornikach, ciecz ze zbiornika B może być powoli wstrzykiwana do zbiornika roboczego.</li>
		<li>Kontroluj prędkość przepływu za pomocą pompy perystaltycznej z prędkością 0,45 ml/s. Należy pamiętać, że całkowity czas napełniania zbiornika roboczego wodą wynosi około 3 h. Oblicz zasolenie na dnie zbiornika roboczego na podstawie29 
		<img alt="Równanie chemiczne: SB = SA(V/2V₀); wzór na stężenie; bilans reakcji." src="/files/ftp_upload/58316/58316eq1.jpg" data-altupdatedat="1788250015000" data-alt="Równanie 1">              (1) 
		gdzie SA, V i V0 odpowiadają odpowiednio zasoleniu zbiornika A, końcowej objętości cieczy w zbiorniku roboczym oraz początkowej objętości cieczy w zbiorniku A (lub B). Wykorzystując zasolenie na dnie SB a świeża woda na górze, częstotliwość wyporności początkowego uwarstwienia N0 jest 
		<img alt="Wzór na częstotliwość wyporności z członami grawitacyjnymi i gradientem gęstości w równaniu termodynamicznym." src="/files/ftp_upload/58316/58316eq2.jpg" data-altupdatedat="1788250015000" data-alt="Równanie 2">              (2) 
		gdzie g jest przyspieszeniem grawitacyjnym, ρ0 jest gęstością odniesienia, a β jest współczynnikiem skurczu zasolenia. Należy zwrócić uwagę na N0 w tym przykładzie wynosi 1,14 rad/s.</li>
	</ol>
	</li>
</ol>
4. Przeprowadzenie eksperymentu
<ol>
	<li>Zdefiniuj warunki brzegowe dla zbiornika roboczego.
	<ol>
		<li>Podłącz komorę wodną górnej płyty do chłodzenia cyrkulacyjnego za pomocą ośmiu równomiernie rozmieszczonych miękkich plastikowych wężyków (o długości 150 cm, średnicy wewnętrznej 10 mm i średnicy zewnętrznej 15 mm). Należy pamiętać, że temperatura górnej płyty zależy od temperatury urządzenia do chłodzenia cyrkulacyjnego. Ustaw temperaturę górnej płyty tak, aby była ona równa temperaturze pokojowej (24 °C).</li>
		<li>Podłącz elektryczną płytkę grzejną znajdującą się wewnątrz dolnej płyty do zasilacza prądu stałego. Należy pamiętać, że podczas tego eksperymentu do czynnika roboczego dostarczany jest stały strumień ciepła, który oblicza się jako 
		<img alt="Równanie mechaniki płynów, \( F_h = \frac{U^2}{RA} \), ilustracja wzoru." src="/files/ftp_upload/58316/58316eq3.jpg" data-altupdatedat="1788250015000" data-alt="Równanie 3">              (3) 
		gdzie U, R i A oznaczają odpowiednio przyłożone napięcie, opór elektryczny i powierzchnię czynną elektrycznej maty grzejnej. W niniejszym przykładzie opór i powierzchnia czynna wynoszą odpowiednio 4,12 Ω oraz 1,89×10-2 m2Ustaw dostarczone napięcie na 60 V, aby całkowity strumień ciepła Fh wynosi 4317 W/m2.</li>
	</ol>
	</li>
	<li>Włącz kamcorder, aby zarejestrować wzorzec przepływu.</li>
	<li>Włącz cyfrowy multimetr oraz wielofunkcyjny system akwizycji danych, aby monitorować temperaturę płyt górnej i dolnej oraz temperaturę i zasolenie cieczy za pomocą urządzenia MSCTI.</li>
	<li>Włącz MPTS, aby przesuwać MSCTI w górę i w dół w celu uzyskania odpowiednich profili temperatury i zasolenia czynnika roboczego.</li>
	<li>Włącz chłodzący cyrkulator oraz zasilacz prądu stałego, aby uzyskać górne i dolne warunki brzegowe cieczy roboczej. 
	Uwaga: Należy pamiętać, że w całym przebiegu eksperymentu wystąpi generowanie, rozwój, łączenie się i zanikanie schodków prądu stałego (DC), a cały proces potrwa około 5 godzin. Po zaniku wszystkich schodków DC należy kolejno wyłączyć zasilacz prądu stałego, chłodnicę cyrkulacyjną, MPTS, cyfrowy multimetr, wielofunkcyjny system akwizycji danych oraz kamerę wideo.</li>
</ol>
5. Przetwarzanie danych
<ol>
	<li>Obraz shadowgraficzny
	<ol>
		<li>Użyj programu Matlab, aby przekonwertować nagranie wideo z kamery na następujące po sobie obrazy do dalszej analizy. Przetwórz obrazy tak, aby uwydatnić wzorzec przepływu wewnątrz zbiornika. Przyjmij intensywność cyfrowego obrazu jako I(x, z), gdzie (x, z) oznacza współrzędne poziome i pionowe z początkiem układu w lewym dolnym rogu obrazu. Należy zauważyć, że I(x, z) zmienia się w zakresie (0, 1) przy 256 poziomach szarości. Znormalizuj każdy obraz za pomocą obrazu tła jako30 
		<img alt="Równanie do analizy natężenia światła, wykazujące stosunek natężeń; istotne dla eksperymentów optycznych." src="/files/ftp_upload/58316/58316eq4.jpg" data-altupdatedat="1788250015000" data-alt="Równanie 4">              (4) 
		gdzie <img alt="Równanie funkcji natężenia I(x, z); powiązane z analizą optyczną lub badaniem rozkładu przestrzennego." src="/files/ftp_upload/58316/58316eq5.jpg" data-altupdatedat="1788250015000" data-alt="Równanie 5"> jest to średnia intensywność obrazu z 10 zdjęć wykonanych przed zastosowaniem chłodzenia i ogrzewania, <img alt="Wzór równowagi statycznej, \(I(x, z)\); równanie do analizy mechanicznej; zastosowanie w badaniach naukowych." src="/files/ftp_upload/58316/58316eq6.jpg" data-altupdatedat="1788250015000" data-alt="Równanie 6"> oznacza natężenie ith obrazu. W ten sposób można usunąć stałe defekty obrazów. Aby zbadać ewolucję czasową wzorca DC, każdy obraz można przekształcić w pojedynczy pionowy profil fluktuacji intensywności, <img alt="Równanie równowagi statycznej δI(z) w formule matematycznej." src="/files/ftp_upload/58316/58316eq7.jpg" data-altupdatedat="1788250015000" data-alt="Równanie 7">, poprzez obliczenie fluktuacji intensywności obrazu (t.j., wartość skuteczna intensywności) wzdłuż kierunku poziomego <img alt="Równanie równowagi statycznej δI(z) w postaci wzoru matematycznego." src="/files/ftp_upload/58316/58316eq7.jpg" data-altupdatedat="1788250015000" data-alt="Równanie 7">Narysuj profile fluktuacji intensywności <img alt="Równanie równowagi statycznej δI(z) w formule matematycznej." src="/files/ftp_upload/58316/58316eq7.jpg" data-altupdatedat="1788250015000" data-alt="Równanie 7"> kolejnych obrazów wraz z upływem czasu, aby pokazać ewolucję schodków prądu stałego (DC).</li>
	</ol>
	</li>
	<li>Profile temperatury i zasolenia
	<ol>
		<li>Należy zauważyć, że w tym doświadczeniu pionowe profile temperatury i zasolenia czynnika roboczego są mierzone przez poruszający się w górę i w dół MSCTI. Oblicz wysokość czasową, h(t), urządzenia MSCTI, przyjmując średnią prędkość ruchu w, czas t oraz czas początkowy t0.0 (odpowiadająca najniższej pozycji), najniższa pozycja hL or najwyższa pozycja hH, jako 
		<img alt="Wzór analizy rozdzielczości czasowej: (t-t₀)/tₚ = n+δ; równanie dla zdarzeń sekwencyjnych, schemat." src="/files/ftp_upload/58316/58316eq8.jpg" data-altupdatedat="1788250015000" data-alt="Równanie 8">              (5) 
		gdzie <img alt="Równanie równowagi statycznej tₚ=(hₕ−hₗ)/w, wzór do obliczania stabilności statycznej." src="/files/ftp_upload/58316/58316eq9.jpg" data-altupdatedat="1788250015000" data-alt="Równanie 9"> jeśli okres przemieszczania MSCTI przebiega od najniższej (najwyższej) do najwyższej (najniższej) pozycji, a n i δ są odpowiednio częścią całkowitą i ułamkową. Następnie oblicz wysokość czasową h(t) jako 
		<img alt="Równania transformaty falkowej; formuła h(t); proces matematyczny; przetwarzanie sygnałów." src="/files/ftp_upload/58316/58316eq10.jpg" data-altupdatedat="1788250015000" data-alt="Równanie 10">              (6) 
		Zauważmy, że w równaniu (6), jeśli n jest liczbą parzystą, MSCTI przemieszcza się w górę; w przeciwnym razie MSCTI przemieszcza się w dół. Należy wykreślić szeregi czasowe temperatury T(t) i zasolenia S(t) w funkcji wysokości h(t), aby otrzymać pionowe profile temperatury i zasolenia.</li>
	</ol>
	</li>
</ol>

<ol><li>Turner, J. S. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aTurner%2C+JS+%22Buoyancy+Effects+in+Fluids%22">Buoyancy Effects in Fluids</a>. , Cambridge Univ. Press, N. Y. 367(1973).</li>
<li>Schmitt, R. W. <a target="_blank" href="https://www.annualreviews.org/doi/abs/10.1146/annurev.fl.26.010194.001351">Double diffusion in oceanography.</a> Annual Review of Fluid Mechanics. 26, 255-285 (1994).</li>
<li>Turner, J. S., Gustafson, L. B. <a target="_blank" href="https://www.sciencedirect.com/science/article/pii/0377027381900184">Fluid motions and compositional gradients produced by crystallization or melting at vertical boundaries.</a> Journal of Volcanology and Geothermal Research. 11, 9S125(1981).</li>
<li>Robb, L. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aRobb%2C+L+%22Introduction+to+Ore-forming+Processes%22">Introduction to Ore-forming Processes</a>. , Blackwell Publishing. 373(2004).</li>
<li>Chabrier, G., Baraffe, I. <a target="_blank" href="http://iopscience.iop.org/article/10.1086/518473/meta">Heat transport in giant (exo)planets: a new perspective.</a> The Astrophysical Journal Letters. 661, 81-84 (2007).</li>
<li>Langlois, W. E. <a target="_blank" href="https://www.annualreviews.org/doi/10.1146/annurev.fl.17.010185.001203">Buoyancy-driven flows in crystal-growth melts.</a> Annual Review of Fluid Mechanics. 17, 191(1985).</li>
<li>Chen, C. -F., Johnson, D. H. <a target="_blank" href="https://www.cambridge.org/core/journals/journal-of-fluid-mechanics/article/doublediffusive-convection-a-report-on-an-engineering-foundation-conference/90E5354AEB593E052485C6D0B4990CB5">Double-diffusive convection: A report on an engineering foundation conference.</a> Journal of Fluid Mechanics. 138, 405-416 (1984).</li>
<li>Griffiths, R. W. <a target="_blank" href="https://www.cambridge.org/core/journals/journal-of-fluid-mechanics/article/layered-double-diffusive-convection-in-porous-media/4D66778DB2E10FA89A691DC5A5E7125B">Layered double-diffusive convection in porous media.</a> Journal of Fluid Mechanics. 102, 221-248 (1981).</li>
<li>You, Y. Z. <a target="_blank" href="https://www.sciencedirect.com/science/article/pii/S0967063702000997">A global ocean climatological atlas of the Turner angle: implications for double-diffusion and water-mass structure.</a> Deep Sea Research Part I: Oceanographic Research Papers. 49, 2075-2093 (2002).</li>
<li>Kelley, D. E., Fernando, H. J. S., Gargett, A. E., Tanny, J., Ozsoy, E. <a target="_blank" href="https://www.sciencedirect.com/science/article/pii/S0079661103000260">The diffusive regime of double diffusive convection.</a> Progress in Oceanography. 56, 461-481 (2003).</li>
<li>Timmermans, M. L., Toole, J., Krishfield, R., Winsor, P. <a target="_blank" href="https://agupubs.onlinelibrary.wiley.com/doi/abs/10.1029/2008JC004829">Ice-Tethered Profiler observations of the double-diffusive staircase in the Canada Basin thermocline.</a> Journal of Geophysical Research: Oceans. 113, 1-10 (2008).</li>
<li>Zhou, S. Q., Lu, Y. Z. <a target="_blank" href="https://agupubs.onlinelibrary.wiley.com/doi/abs/10.1002/2013JC009141">Characterization of double diffusive convection steps and heat budget in the deep Arctic Ocean.</a> Journal of Geophysical Research: Oceans. 118 (12), 6672-6686 (2013).</li>
<li>Turner, J. S. <a target="_blank" href="https://journals.ametsoc.org/doi/10.1175/2009JPO4279.1">The melting of ice in the arctic ocean: The influence of double-diffusive transport of heat from below.</a> Journal of Physical Oceanography. 40, 249-256 (2010).</li>
<li>Neal, V. T., Neshyba, S., Denner, W. <a target="_blank" href="https://www.ncbi.nlm.nih.gov/pubmed/17796548">Thermal stratification in the Arctic Ocean.</a> Science. 166 (3903), 373-374 (1969).</li>
<li>Padman, L., Dillon, T. M. <a target="_blank" href="https://agupubs.onlinelibrary.wiley.com/doi/abs/10.1029/JC092iC10p10799">Vertical heat fluxes through the Beaufort Sea thermohaline staircase.</a> Journal of Geophysical Research: Oceans. 92 (C10), 10799-10806 (1987).</li>
<li>Sirevaag, A., Fer, I. <a target="_blank" href="https://agupubs.onlinelibrary.wiley.com/doi/full/10.1029/2012JC007910">Vertical heat transfer in the Arctic Ocean: The role of double-diffusive mixing.</a> Journal of Geophysical Research: Oceans. 117 (C7), (2012).</li>
<li>Guthrie, J. D., Fer, I., Morison, J. <a target="_blank" href="https://agupubs.onlinelibrary.wiley.com/doi/abs/10.1002/2015JC010884">Observational validation of the diffusive convection flux laws in the Amundsen Basin, Arctic Ocean.</a> Journal of Geophysical Research: Oceans. 120 (12), 7880-7896 (2015).</li>
<li>Bebieva, Y., Timmermans, M. L. <a target="_blank" href="https://agupubs.onlinelibrary.wiley.com/doi/abs/10.1002/2015JC011105">An examination of double-diffusive processes in a mesoscale eddy in the Arctic Ocean.</a> Journal of Geophysical Research: Oceans. 121 (1), 457-475 (2016).</li>
<li>Shibley, N. C., Timmermans, M. L., Carpenter, J. R., Toole, J. M. <a target="_blank" href="https://agupubs.onlinelibrary.wiley.com/doi/full/10.1002/2016JC012419">Spatial variability of the Arctic Ocean's double-diffusive staircase.</a> Journal of Geophysical Research: Oceans. 122 (2), 980-994 (2017).</li>
<li>Schmid, M., Busbridge, M. <a target="_blank" href="https://aslopubs.onlinelibrary.wiley.com/doi/abs/10.4319/lo.2010.55.1.0225">Double-diffusive convection in Lake Kivu.</a> Limnology and Oceanography. 55 (1), 225-238 (2010).</li>
<li>Sommer, T., et al. <a target="_blank" href="https://agupubs.onlinelibrary.wiley.com/doi/abs/10.1002/2013JC009166">Interface structure and flux laws in a natural double-diffusive layering.</a> Journal of Geophysical Research: Oceans. 118 (11), 6092-6106 (2013).</li>
<li>Carpenter, J. R., Sommer, T., Wüest, A. <a target="_blank" href="https://www.cambridge.org/core/journals/journal-of-fluid-mechanics/article/simulations-of-a-double-diffusive-interface-in-the-diffusive-convection-regime/63D2ECE2AA41439E01A01F9A0D76F2E2">Simulations of a double-diffusive interface in the diffusive convection regime.</a> Journal of Fluid Mechanics. 711, 411-436 (2012).</li>
<li>Flanagan, J. D., Lefler, A. S., Radko, T. <a target="_blank" href="https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1002/grl.50440">Heat transport through diffusive interfaces.</a> Geophysical Research Letters. 40 (10), 2466-2470 (2013).</li>
<li>Radko, T., Flanagan, J. D., Stellmach, S., Timmermans, M. L. <a target="_blank" href="https://journals.ametsoc.org/doi/10.1175/JPO-D-13-0156.1">Double-diffusive recipes. Part II: Layer-merging events.</a> Journal of Physical Oceanography. 44 (5), 1285-1305 (2014).</li>
<li>Scheifele, B., Pawlowicz, R., Sommer, T., Wüest, A. <a target="_blank" href="https://journals.ametsoc.org/doi/abs/10.1175/JPO-D-14-0070.1"> Double diffusion in saline Powell Lake, British Columbia.</a> Journal of Physical Oceanography. 44 (11), 2893-2908 (2014).</li>
<li>Guo, S. X., Zhou, S. Q., Qu, L., Lu, Y. Z. <a target="_blank" href="https://agupubs.onlinelibrary.wiley.com/doi/full/10.1002/2016JC012172">Laboratory experiments on diffusive convection layer thickness and its oceanographic implications.</a> Journal of Geophysical Research: Oceans. 121 (10), 7517-7529 (2016).</li>
<li>Guo, S. X., Cen, X. R., Zhou, S. Q. <a target="_blank" href="https://agupubs.onlinelibrary.wiley.com/doi/abs/10.1002/2017JC013545">New parametrization for heat transport through diffusive convection interface.</a> Journal of Geophysical Research: Oceans. 123 (2), 1327-1338 (2018).</li>
<li>Turner, J. S. <a target="_blank" href="https://www.sciencedirect.com/science/article/pii/0017931065900220">The coupled turbulent transports of salt and heat across a sharp density interface.</a> International Journal of Heat and Mass Transfer. 8 (5), 759-767 (1965).</li>
<li>Hill, D. F. <a target="_blank" href="https://link.springer.com/article/10.1007/s00348-001-0376-5">General density gradients in general domains: the &#34;two-tank&#34; method revisited.</a> Experiments in Fluids. 32 (4), 434-440 (2002).</li>
<li>Zhou, S. Q., Ahlers, G. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Spatiotemporal+chaos+in+electroconvection+of+a+homeotropically+aligned+nematic+liquid+crystal%5BTitle%5D)+AND+%22Physical+Review+E%22%5BJournal%5D)">Spatiotemporal chaos in electroconvection of a homeotropically aligned nematic liquid crystal.</a> Physical Review E. 74 (4), 046212(2006).</li>
</ol>

Autorzy nie mają nic do ujawnienia.

W niniejszej pracy opisano szczegółowy protokół eksperymentalny służący do symulacji termohalinicznych struktur schodkowych typu DC (double-diffusive convection) w zbiorniku prostokątnym. Początkowe liniowe rozwarstwienie gęstości czynnika roboczego uzyskano metodą dwóch zbiorników. Górna płyta była utrzymywana w stałej temperaturze, natomiast dolna przy stałym strumieniu ciepła. Cały proces ewolucji schodków DC, obejmujący ich powstawanie, rozwój, łączenie się i zanikanie, został zwizualizowany za pomocą techniki shadow...

Konwekcja podwójnej dyfuzji (DDC) jest jednym z najważniejszych procesów mieszania pionowego. Występuje, gdy pionowy rozkład gęstości warstwowego słupa wody jest kontrolowany przez dwa lub więcej gradientów składowych skalarnych o przeciwnych kierunkach, gdzie składniki mają wyraźnie różne dyfuzyjności molekularne1. Powszechnie występuje w oceanografii2, atmosphere3, geology4, astrophysics5, material science6, metalurgii7, oraz inżynierii architektonicznej8. DDC jest obecny w prawie połowie globalnego oceanu i ma istotny wpływ na wieloskalowe procesy oceaniczne, a nawet zmiany klimatyczne9.
Istnieją dwa podstawowe tryby dla DDC: palec solny (SF) i konwekcja dyfuzyjna (DC). SF występuje, gdy ciepła, słona masa wody nakłada się na chłodniejszą, świeższą wodę w środowisku warstwowym. Kiedy ciepła i słona woda znajdzie się poniżej zimnej i słodkiej wody, utworzy się prąd stały. Niezwykłą cechą DC jest to, że pionowe profile temperatury, zasolenia i gęstości przypominają schody, składają się z naprzemiennych jednorodnych warstw konwekcyjnych i cienkich, silnie uwarstwionych granic. DC występuje głównie w oceanach na dużych szerokościach geograficznych i niektórych wewnętrznych słonych jeziorach, takich jak Ocean Arktyczny i Antarktyczny, Morze Ochockie, Morze Czerwone i Afrykańskie Jezioro Kiwu10. Na Oceanie Arktycznym w górnych i głębokich oceanach znajdują się ogólnobasenowe i trwałe schody DC11,12. Ma on istotny wpływ na mieszanie się diapiknów w górnym oceanie i może znacząco wpływać na topnienie lodu, co ostatnio budzi coraz większe zainteresowanie w społeczności oceanograficznej13.
Struktura schodów DC została po raz pierwszy odkryta na Oceanie Arktycznym w 1969 roku14. Następnie Padman & Dillon15, Timmermans et al. 11, Sirevaag & Fer16, Zhou & Lu12, Guthrie et al. 17, Bebieva & Timmermans18, oraz Shibley et al. 19 zmierzył klatki schodowe DC w różnych basenach Oceanu Arktycznego, w tym pionowe i poziome skale warstwy konwekcyjnej i interfejsu, głębokość i całkowitą grubość schodów, pionowy transfer ciepła, procesy DC w mezoskalowych wirach oraz czasowe i przestrzenne zmiany konstrukcji schodów. Schmid i wsp. 20 oraz Sommer et al.21 obserwował klatki schodowe DC za pomocą profilera mikrostruktury w jeziorze Kivu. Przedstawili oni główne cechy struktury i strumienie ciepła prądu stałego oraz porównali zmierzone strumienie ciepła z istniejącym wzorem parametrycznym. Wraz z poprawą szybkości przetwarzania komputerowego, symulacje numeryczne prądu stałego zostały ostatnio wykonane, na przykład w celu zbadania struktury i niestabilności interfejsu, wymiany ciepła przez interfejs, zdarzenia łączenia warstw i tak dalej22,23,24.
Obserwacja terenowa znacznie poszerzyła zrozumienie oceanicznego DC dla oceanografów, ale pomiar jest mocno ograniczony przez nieokreślone środowiska i instrumenty przepływu oceanicznego. Na przykład, interfejs DC ma bardzo małą skalę pionową, cieńszą niż 0,1 m w niektórych jeziorach i oceanach25 i potrzebne są specjalne instrumenty o wysokiej rozdzielczości. Eksperyment laboratoryjny pokazuje swoje unikalne zalety w badaniu podstawowych dynamicznych i termodynamicznych praw prądu stałego. Dzięki eksperymentowi laboratoryjnemu można obserwować ewolucję schodów DC, mierzyć temperaturę i zasolenie oraz proponować pewne parametryzacje dla zastosowań oceanicznych26,27. Ponadto w eksperymencie laboratoryjnym kontrolowane parametry i warunki można łatwo dostosować w razie potrzeby. Na przykład Turner po raz pierwszy przeprowadził symulację schodów DC w laboratorium w 1965 roku i zaproponował parametryzację wymiany ciepła przez interfejs dyfuzyjny, która była często aktualizowana i szeroko stosowana w obserwacjach oceanicznych in situ28.
W tym artykule opisano szczegółowy protokół eksperymentalny, który symuluje proces ewolucji schodów DC, w tym powstawanie, rozwój i zanikanie w warstwowej solanki ogrzewanej od dołu. Temperatura i zasolenie są mierzone za pomocą przyrządu w skali mikro, a schody DC są monitorowane techniką wykresu cienia. Szczegółowo opisano konfigurację eksperymentalną, proces ewolucji, analizę danych i dyskusję na temat wyników. Zmieniając warunki początkowe i brzegowe, obecna konfiguracja i metoda eksperymentalna mogą być wykorzystane do symulacji innych zjawisk oceanicznych, takich jak oceaniczna konwekcja pozioma, głębokowodne erupcje hydrotermalne, pogłębianie powierzchni warstwy mieszanej, wpływ podmorskiej geotermii na cyrkulację oceaniczną i tak dalej.

<table><tbody><tr>Zbiornik <td></td><td></td><td>prostokątny</td><td>Część wykonana na zamówienie</td></tr><tr><td>Pleksiglas</td><td></td><td></td><td>Część wykonana na zamówienie</td></tr><tr><td>Elektryczna poduszka</td><td></td><td></td> grzewcza <td>Część wykonana na zamówienie</td></tr><tr><td>Woda</td><td></td><td></td> destylowana <td>Wielu dostawców</td></tr><tr><td>Stół optyczny</td><td>Liansheng Inc.</td><td>Termiostory MRT-P / B</td><td></td></tr><tr><td></td><td>Część wykonana na zamówienie</td>Multimetr cyfrowy Keithley Inc Model 2700 Przyrząd do pomiaru przewodności i temperatury w mikroskali (MSCTI) PME. Inc.<td></td><td></td></tr><tr><td></td><td></td><td></td><td></td></tr><tr><td></td><td></td><td>Model 125</td><td></td></tr><tr><td>Wielofunkcyjna akwizycja danych (MDA)</td><td>MCC. Inc.</td><td>USB-2048</td><td></td></tr><tr><td>Zmotoryzowany etap precyzyjnego tłumaczenia (MPTS)</td><td>Thorlabs Inc.</td><td>LTS300</td><td></td></tr><tr><td>Kalka</td><td></td><td></td> kreślarska <td>Wielu dostawców Lampa</td></tr><tr><td></td><td></td><td></td> LED <td>Wielu dostawców</td></tr><tr><td>Kamera</td><td>Sony Inc.</td><td>XDR-XR550</td><td></td></tr><tr><td>Odgazowana słodka woda</td><td></td><td></td><td>Część wykonana na zamówienie</td></tr><tr><td>Woda</td><td></td><td></td> solankowa <td>Część wykonana na zamówienie</td></tr><tr><td>Elastyczna rura</td><td></td><td></td><td>Wielu dostawców</td></tr><tr><td>Elektryczne mieszadło magnetyczne&nbsp;</td><td>Meiyingpu Inc.</td><td>MYP2011-100</td><td></td></tr><tr><td>Pompa perystaltyczna</td><td>Zhisun Inc.</td><td>DDBT-201</td><td></td></tr><tr><td>Chłodnicza pompa cyrkulacyjna</td><td>Polyscience Inc.</td><td>Model 9702</td><td></td></tr><tr><td>Plastikowa miękka rurka</td><td></td><td></td><td>Wielu dostawców</td></tr><tr><td>Zasilacz prądu stałego</td><td>GE Inc.</td><td>GPS-3030</td><td></td></tr><tr><td>Matlab</td><td>MathWorks Inc.</td><td>Czynnik chłodniczy R2012a</td><td></td></tr></tbody></table>

ewolucja konstrukcji klatek schodowych w konwekcji dyfuzyjnej

Evolution of Staircase Structures in Diffusive Convection

Konwekcja dyfuzyjna (DC) występuje, gdy pionowa gęstość warstwowa jest kontrolowana przez dwa przeciwstawne gradienty skalarne, które mają wyraźnie różne dyfuzyjności molekularne, a większe i mniejsze gradienty skalarne dyfuzyjności mają odpowiednio ujemny i dodatni wpływ na rozkład gęstości. DC występuje w wielu procesach naturalnych i zastosowaniach inżynieryjnych, na przykład w oceanografii, astrofizyce i metalurgii. W oceanach jedną z najbardziej niezwykłych cech prądu stałego jest to, że pionowe profile temperatury i zasolenia są strukturą przypominającą schody, składającą się z następujących po sobie stopni z grubymi jednorodnymi warstwami konwekcji oraz stosunkowo cienkimi i wysokogradientowymi interfejsami. Schody DC zostały zaobserwowane w wielu oceanach, zwłaszcza w Oceanie Arktycznym i Antarktyce, i odgrywają ważną rolę w cyrkulacji oceanicznej i zmianach klimatycznych. Na Oceanie Arktycznym w górnych i głębokich oceanach znajdują się ogólnobasenowe i trwałe schody prądu stałego. Proces prądu stałego ma istotny wpływ na mieszanie się diapiknów w górnym oceanie i może znacząco wpływać na topnienie lodu powierzchniowego. W porównaniu z ograniczeniami obserwacji terenowych, eksperyment laboratoryjny pokazuje swoją unikalną zaletę w skutecznym badaniu procesów dynamicznych i termodynamicznych w prądzie stałym, ponieważ warunki brzegowe i kontrolowane parametry można ściśle regulować. W tym miejscu opisano szczegółowy protokół symulujący proces ewolucji struktury schodów prądu stałego, w tym jej powstawanie, rozwój i zanikanie, w prostokątnym zbiorniku wypełnionym warstwową wodą solankową. Szczegółowo opisano konfigurację eksperymentalną, proces ewolucji, analizę danych i dyskusję na temat wyników.

Rycina 1 przedstawia schemat układu eksperymentalnego. Jego komponenty zostały opisane w protokole. Główne części pokazano na Rycinnie 1a, a szczegółowy widok zbiornika roboczego na Rycinnie 1b. Rycina 2 przedstawia zmiany temperatury na płytce dolnej (Tb, krzywa czerwona) i górnej (Tt, krzywa czarna). Wskazano, że początkowo temperatury obu płytek są niemal identyczne z temperatur...

Konwekcja dyfuzyjna (DC) występuje powszechnie w naturalnych procesach i zastosowaniach inżynieryjnych, charakteryzując się szeregiem schodów z jednorodnymi warstwami konwekcji i warstwowymi interfejsami. Opisano procedurę eksperymentalną mającą na celu symulację procesu ewolucji struktury schodów prądu stałego, w tym powstawania, rozwoju i zaniku, w prostokątnym zbiorniku.

Ewolucja konstrukcji klatek schodowych w konwekcji dyfuzyjnej

Diffusive convection (DC) occurs when the vertical stratified density is controlled by two opposing scalar gradients that have distinctly different ...

Metoda ta może być wykorzystana do zbadania kluczowych zagadnień w dynamicznych procesach konstrukcji schodów dyfuzyjnych konwekcyjnych, takich jak ich powstawanie, rozwój i zanikanie. Przy obecnym układzie eksperymentalnym możemy wyraźnie zobaczyć ewolucję konwekcji dyfuzyjnej, która jest bardzo trudna, jeśli nie niemożliwa, do zaobserwowania bezpośrednio w oceanie. Zbierz elementy działającego zbiornika.Należą do nich miedziane płyty górne i dolne oraz akrylowy zespół ścian bocznych. Górna płyta jest galwanizowana i zamyka komorę wodną dostępną za pomocą rurek. Posiada również miejsca na termistory.Dolna płyta jest również galwanizowana i posiada szczeliny na termistory. Zawiera poduszkę grzewczą. Użyj wody destylowanej, aby dokładnie wyczyścić miedziane płyty i akrylowe ścianki boczne.Po zakończeniu zmontuj zbiornik za pomocą, aby upewnić się, że jest wodoszczelny. Ten ukończony zbiornik ma 257 milimetrów długości i wysokości oraz 65 milimetrów szerokości. Na stole optycznym ustaw ramę nośną ze stali nierdzewnej.Na górze ramy umieść płytę izolacyjną. Gdy wszystko jest na miejscu, umieść zmontowany zbiornik na płycie. Włóż termistory do górnej i dolnej płyty i podłącz je do systemu akwizycji danych.Termistor monitoruje temperaturę płytki, w której jest umieszczony. Następnie przesuń precyzyjny pionowy stolik translacyjny na miejsce w celu zamontowania sondy. Przez górną płytę umieść w zbiorniku czujnik przyrządu do pomiaru przewodności i temperatury w skali mikro.Przymocuj instrument do etapu tłumaczenia precyzyjnego. Teraz skonfiguruj oprogramowanie do akwizycji danych z czujników i odczytów termistorów. Wróć do zbiornika, aby wyregulować położenie czujnika przewodności i temperatury.Ustaw początkową pozycję czujnika w jego najniższym punkcie, tutaj 20 milimetrów nad dnem zbiornika. Następnie ustaw parametry ruchu dla etapu translacji dla eksperymentu. Użyj techniki wykresu cienia, aby monitorować eksperyment.W tym celu przymocuj kawałek kalki technicznej na zewnątrz po jednej stronie zbiornika. Po przeciwnej stronie zbiornika, w odległości około pięciu metrów, umieść źródło światła. Aby wytworzyć prawie skolimowaną wiązkę, użyj wąskiej wiązki LED jako źródła światła.Umieść szybką kamerę około metra przed kalką techniczną na ścieżce wiązki. Gdy lampa i kamera są włączone, dostosuj ich położenie, aby uzyskać wyraźny obraz. Uzyskaj dwa identyczne prostokątne zbiorniki na płyn roboczy.Zbiorniki są tej samej wielkości co zbiornik roboczy. Połącz je u podstawy za pomocą zaciśniętej elastycznej rurki o długości 10 centymetrów. Umieść dwa zbiorniki na tej samej wysokości.Zadbaj o to, aby w jednym zbiorniku znajdowało się mieszadło elektryczne. Następnie zaciśnij elastyczną 50-centymetrową rurkę w pompie perystaltycznej. Użyj rurki, aby połączyć zbiornik z mieszadłem ze zbiornikiem roboczym.Napełnij zbiornik, który nie ma mieszadła, roztworem soli o stężeniu 60 gramów na kilogram. Napełnij drugi zbiornik równą objętością odgazowanej świeżej wody. Po napełnieniu zbiorników odblokuj rurkę łączącą.Ciągłą homogenizację wody zarobowej za pomocą mieszadła elektrycznego. Kontroluj natężenie przepływu w zbiorniku roboczym za pomocą pompy perystaltycznej. Zbiornik roboczy potrzebuje około trzech godzin, aby się napełnić i być gotowym do eksperymentu.Do eksperymentu użyj cyrkulatora z chłodzeniem, aby ustawić warunki brzegowe górnej płyty. Na górnej płycie podłącz komorę wodną do cyrkulatora chłodniczego za pomocą ośmiu miękkich plastikowych rurek. W pompie cyrkulacyjnej ustaw temperaturę górnej płyty na temperaturę pokojową.Teraz podłącz poduszkę grzewczą dolnej płyty do źródła zasilania prądem stałym i ustaw jej poziom mocy. Włącz kamerę, aby nagrać wzorzec przepływu. Rozpocznij monitorowanie danych dotyczących temperatury i zasolenia.Następnie rozpocznij ruch pionowy czujnika. Na koniec włącz cyrkulator chłodniczy i zasilacz prądu stałego, aby uzyskać warunki brzegowe płynu roboczego. Jest to przykład obrazu cienia wykonanego, gdy górna płyta ma temperaturę pokojową, a dolna płyta jest podgrzewana.Na obrazie znajdują się trzy warstwy konwekcyjne, w których gęstość płynu jest jednorodna. Istnieją również trzy warstwy interfejsu, w których występują gradienty o dużej gęstości. Ten profil fluktuacji intensywności ma trzy szczyty, które odpowiadają interfejsom.Oto ewolucja profilu fluktuacji intensywności w czasie, która ujawnia powstawanie, rozwój i zanikanie warstw związanych z dyfuzyjnymi schodami konwekcyjnymi. Profil temperatury i profil zasolenia dostarczają innych poglądów na ewolucję systemu w czasie. W tych profilach każda linia ciągła reprezentuje dane zebrane w ciągu 404 sekund.Każdy zestaw danych w profilu temperatury jest przesunięty o 1,5 stopnia Celsjusza w stosunku do poprzedniego. W przypadku profilu zasolenia każdy zestaw danych jest przesunięty o trzy gramy na kilogram w stosunku do poprzedniego zestawu. Metoda ta może być również stosowana do symulacji innych zjawisk w oceanach, takich jak globalna cyrkulacja oceaniczna, pogłębianie warstw mieszanych i erupcja pióropuszy hydrotermalnych.Jeśli chcesz wypróbować tę procedurę, pamiętaj, aby najpierw ustawić czujniki przewodności i temperatury w najniższej pozycji stolika translacyjnego, aby zapobiec rozbiciu się czujnika o dolną płytę. Dzięki przedstawionym wynikom możemy dalej rozwijać nową parametryzację grubości warstwy, strumienia ciepła i dodawać dyfuzyjność w konwekcji dyfuzyjnej, co utorowałoby drogę społeczności oceanografii fizycznej do badania zjawisk w zastosowaniu oceanicznym.

evolution-of-staircase-structures-in-diffusive-convection

Research

JoVE Journal

Environment

Środowisko

Implanty stentów w zwężonych krzywiznach tętnic są podatne na awarie "typu IV" polegające na całkowitym złamaniu poprzecznym stentów i liniowym przemieszczeniu złamanych części. Przedstawiamy protokół wykrywania struktur przepływu wtórnego (wirowego) w modelu tętnicy zakrzywionej, po klinicznie istotnych awariach stentów "typu IV".

Eksperymentalne badanie wtórnych struktur przepływu za modelowym uszkodzeniem stentu typu IV w odcinku testowym tętnicy zakrzywionej pod kątem 180°

Tutaj prezentujemy protokół, który demonstruje powstawanie lodu, gdy woda jest wprowadzana do zimnej kąpieli solanki, jako wtórny czynnik chłodniczy, w zakresie temperatur znacznie poniżej temperatury zamarzania wody. Może być stosowany jako alternatywny sposób produkcji lodu dla przemysłu.

Wytwarzanie lodu a zjawiska wymiany ciepła i masy spowodowane wprowadzaniem wody do zimnej kąpieli solankowej

Przedstawiono protokół badania dyfuzji pasywnych znaczników w przepływie napędzanym ciśnieniem laminarnym. Procedura ma zastosowanie do różnych geometrii rur kapilarnych.

Dyfuzja pasywnych znaczników w laminarnym przepływie ścinającym

Prosty protokół jest dostępny do wytwarzania struktur hemiwickich o różnych rozmiarach, kształtach i materiałach. Protokół wykorzystuje kombinację fizycznego tłoczenia, formowania PDMS i cienkowarstwowych modyfikacji powierzchni za pomocą popularnych technik osadzania materiałów.

Skalowalne drukowanie pieczątek i produkcja powierzchni hemiwickingowych

Celem tego protokołu jest tworzenie symulowanych kominów hydrotermalnych za pomocą eksperymentów z wstrzykiwaniem w ogrodzie chemicznym i wprowadzenie gradientu termicznego przez nieorganiczną membranę osadu, za pomocą kondensatora do druku 3D, który można odtworzyć w celach edukacyjnych.

Symulacja kominów hydrotermalnych wczesnej Ziemi w środowisku gradientu termicznego

Ten artykuł opisuje, jak można stworzyć trzy różne reżimy wirów wodnych w hiperbolicznym lejku Schaubergera, ich najważniejsze cechy oraz jak można obliczyć powiązane parametry, takie jak szybkość transferu tlenu.

Przygotowanie hiperbolicznych wirów wodnych o swobodnej powierzchni

Spaliny z elektrowni są tanim źródłem CO2 dla wzrostu glonów. Zbudowaliśmy prototypowe systemy "od gazów spalinowych do hodowli glonów" i opisaliśmy, jak zwiększyć skalę procesu hodowli glonów. Wykazaliśmy zastosowanie modelu bioreakcji z przeniesieniem masy do symulacji i zaprojektowania optymalnej pracy gazów spalinowych do wzrostu Chlorella sp. w fotobioreaktorach glonów.

Optymalizacja ustawień spalin w celu promowania wzrostu mikroalg w fotobioreaktorach za pomocą symulacji komputerowych

Konsorcja mikrobiologiczne w ulach trzmieli wzbogacają i zachowują pyłek dla larw pszczół. Wykorzystując sekwencjonowanie nowej generacji, wraz z eksperymentami laboratoryjnymi i terenowymi, manuskrypt ten opisuje protokoły stosowane do testowania hipotezy, że pozostałości fungicydów zmieniają mikrobiom pyłku i demografię kolonii, ostatecznie prowadząc do utraty kolonii.

Techniki empiryczne, metagenomiczne i obliczeniowe rzucają światło na mechanizmy, za pomocą których fungicydy zagrażają zdrowiu pszczół

Obecność toksyn cyjanobakteryjnych w zbiornikach słodkiej wody przeznaczonej do spożycia przez ludzi jest głównym problemem dla organów zarządzających gospodarką wodną. Aby ocenić ryzyko zanieczyszczenia wody, w artykule opisano protokół wykrywania w terenie szczepów sinic w próbkach ciekłych i stałych przy użyciu chipa mikromacierzy przeciwciał.

Eksperymentalny protokół wykrywania sinic w próbkach ciekłych i stałych za pomocą chipa mikromacierzy przeciwciał

Ta praca opisuje nowy protokół oceny ekotoksyczności zanieczyszczeń, w tym pojawiających się zanieczyszczeń, takich jak nanomateriały, przy użyciu morskiej bakterii Vibrio anguillarum. Metoda ta pozwala na określenie LC50, czyli śmiertelności, stężenia, które powoduje 50% spadek zdolności hodowli bakterii, po 6 godzinach ekspozycji.

Metoda ekotoksykologiczna z bakteriami morskimi Vibrio anguillarum do oceny ostrej toksyczności zanieczyszczeń środowiska

Poznaj metodę symulacji ewolucji struktur klatek schodowych w konwekcji dyfuzyjnej w układzie laboratoryjnym z płytami miedzianymi i termistorami. To podejście umożliwia obserwację procesów niemożliwych do bezpośredniego zobaczenia w oceanie.

Ewolucja konstrukcji klatek schodowych w konwekcji dyfuzyjnej

Przeglądaj więcej filmów