Projekt badawczy
Przetwarzanie aerozolu przez chmury - budowa wszechstronnego zaprogramowanego obiektowo narzędzia do symulacji numerycznych
prof. dr hab. Hanna Pawłowska | Kierownik projektu |
dr Sylwester Arabas | Współwykonawca |
dr Piotr Dziekan | Współwykonawca |
dr Anna Jaruga | Współwykonawca |
dr Maciej Waruszewski | Współwykonawca |
Cel badań
Celem projektu jest stworzenie modelu numerycznego do badania oddziaływań chmura-aerozol. Projekt jest odpowiedzią na palącą potrzebę znalezienia efektywnego narzędzia badawczego, gdyż istniejące obecnie, nawet najbardziej wyszukane ale tradycyjne metody opisu mikrofizyki, okazują się albo niewystarczające albo niepraktyczne.
Aerozole i chmury wpływają na bilans energetyczny Ziemi i odgrywają ważną, aczkolwiek ciągle słabo zrozumiałą rolę w prognozach zmian przyszłego klimatu (IPCC, 2007). Chmury wpływają na aerozol przez jego wymywanie oraz poprzez regenerację jąder kondensacji w wyniku odparowania kropelek chmurowych i mżawki. Zrozumienie w jaki sposób chmury wpływają na aerozol jest niezbędne ażeby wyczerpująco opisać wpływ aerozoli, chmur i opadu na radiację. Złożoność zachodzących oddziaływań sprawia, że jedynym sposobem poszukiwania odpowiedzi na pytanie jak naturalne i antropogeniczne aerozole wpływają na środowisko stają się symulacje numeryczne.
Metoda badawcza
W ramach projektu zostanie stworzony model numeryczny oparty na schemacie symulacji dużych wirów (ang. Large Eddy Simulation, LES) opisujący w wyczerpujący sposób procesy uczestniczące w przetwarzaniu aerozolu przez chmury. Model będzie stanowił rozszerzenie istniejącego pakietu oprogramowania (icicle) rozwijanego w Instytucie Geofizyki na Wydziale Fizyki Uniwersytetu Warszawskiego (IGFUW). Dotychczas icicle pozwalał na prowadzenie dwu-wymiarowych kinematycznych symulacji o zadanym przepływie, w którym mikrofizyka aerozoli i chmur jest oparta na śledzeniu cząstek. Świeżo stworzony moduł dynamiczny icicle będzie oparty na implementacji kluczowych komponentów NCAR-owskiego modelu chmurowego EULAG. Nowatorska natura projektu polega na połączeniu bardzo efektywnego opisu dynamiki (bazującego na algorytmie MPDATA modelu EULAG) z nowoczesnym opisem mikrofizyki aerozoli, kropel chmurowych i opadowych metodą Super Droplet. Warto podkreślić, że zarówno algorytm rozwiązywania nieelastycznych równań dynamiki oparty na MPDATA, jak i opis mikrofizyki bazujący na śledzeniu cząstek nie są powszechnie używane w obecnie rozwijanych modelach chmurowych, częściowo ze względu na duży stopień złożoności tych metod.
W ramach projektu planujemy również rozszerzyć mikrofizyczny moduł icicle oparty na metodzie Super Droplet o możliwość opisu reakcji chemicznych zachodzących w kroplach.
Projekt będzie realizowany w zgodzie z zasadami typowymi dla wolnego i otwartego oprogramowania oraz wykorzystując techniki programowania obiektowego.
Wpływ rezultatów
Otwartość i przejrzystość kodu modelu oraz pełna dokumentacja sprawią, że będzie on czytelny i możliwy do łatwego wykorzystania w środowisku naukowym, oferując niniejszym znakomite narzędzie do wykorzystania w badaniach.
Projekt ma zarówno cel naukowy jak edukacyjny. Stworzony model pozwoli na uzyskanie nowych wyników naukowych będących przedmiotem licznych publikacji w recenzowanych prestiżowych czasopismach. Powstały program komputerowy, wraz z kompletna dokumentacją techniczną, dostarczy narzędzia dla przyszłych magistrantów i doktorantów.
Współpraca międzynarodowa
Nowe moduły icicle będą budowane we współpracy z dr. Piotrem Smolarkiewiczem i dr. hab. Wojciechem Grabowskim z National Center for Atmospheric Research z Boulder, Colorado, USA. Stanowić będą nową oryginalną implementację wersji LES nieelastycznego numerycznego modelu chmurowego EULAG stworzonego i rozwijanego przez nich w NCAR. Partnerzy projektu są uznanymi w świecie naukowym ekspertami w dziedzinie zastosowań obliczeniowej dynamiki płynów w szerokiej gamie problemów atmosferycznych. Mają bogatą historię współpracy z IGFUW. Obecny projekt jest kolejnym krokiem w realizacji programu badawczego, który jest wzajemnie korzystny zarówno dla NCAR jak IGFUW i otwiera nowy rozdział współpracy, gdyż po raz pierwszy podejmuje wyzwanie budowy nowego narzędzia badawczego.
- Grabowski W. W., Dziekan P., and Pawlowska H., 2018, Lagrangian condensation microphysics with Twomey CCN activation, Geoscientific Model Development, vol. 11, 103-120, 10.5194/gmd-11-1-2018
- Jaruga A., and Pawlowska H., 2018, libcloudph++ 2.0: aqueous-phase chemistry extension of particle-based cloud microphysics scheme, Geoscientific Model Development, vol. 11, 3623-3645, 10.5194/gmd-11-3623-2018
- Waruszewski, M., Kühnlein, C., Pawlowska, H., Smolarkiewicz, P.K., 2018, MPDATA: Third-order accuracy for variable flows, Journal of Computational Physics, vol. 359, 361-379, 10.1016/j.jcp.2018.01.005
- Arabas S. and Shima S., 2017, On the CCN (de)activation nonlinearities, Nonlinear Processes in Geophysics, vol. 24, 535-542, 10.5194/npg-24-535-2017
- Dziekan P. and Pawlowska H., 2017, Stochastic coalescence in Lagrangian cloud microphysics, Atmospheric Chemistry and Physics, vol. 17 (2), 13509-13520, 10.5194/acp-17-13509-2017
- Arabas S., Jaruga A., Pawlowska H., and Grabowski W. W., 2015, libcloudph++ 1.0: a single-moment bulk, double-moment bulk, and particle-based warm-rain microphysics library in C++, Geoscientific Model Development, vol. 8, 1677-1707, doi: 10.5194/gmd-8-1677-2015
- Jarecka D., Jaruga A., Smolarkiewicz P.K., 2015, A spreading drop of shallow water, Journal of Computational Physics, vol. 289, 53–61, 10.1016/j.jcp.2015.02.003
- Jaruga A., Arabas S., Jarecka D., Pawlowska H., Smolarkiewicz P. K., Waruszewski M., 2015, libmpdata++ 1.0: a library of parallel MPDATA solvers for systems of generalised transport equations, Geoscientific Model Development, vol. 8, 1005-1032, doi:10.5194/gmd-8-1005-2015