dizertační práce

Rozvoj lattice Boltzmannovy metody pro řešení transportu a depozice vláknových částic

Text práce 11.6 MB Teze 7.95 MB

Autor práce: Ing. František Prinz

Ak. rok: 2025/2026

Vedoucí: prof. Ing. Miroslav Jícha, CSc.

Oponenti: doc. Ing. Radek Fučík, Ph.D., prof. Ing. Jan Vimmr, Ph.D.

Abstrakt:

Transport a depozice aerosolů v lidském dýchacím traktu jsou aktivní oblastí výzkumu.
Inhalace aerosolů je jednou z metod dodávání léčiv do lidského těla. Na druhé straně však
vdechovaný vzduch může obsahovat také škodlivé částice, jako jsou pevné prachové částice.
Kromě sférických aerosolů jsou zkoumány také vláknité částice. Díky svému nepravidelnému
tvaru mohou pronikat hlouběji do dýchacího systému a mohou být účinnější pro lékařské
účely. Nicméně jejich složitý tvar představuje výzvu pro simulace, protože odporová síla
závisí na aktuální orientaci částice vůči proudění. Navíc v oblasti bifurkací mají vláknité
částice tendenci měnit orientaci vůči proudění. Metoda Euler-Lagrange Euler-Rotation
(ELER) počítá pohyb částic pomocí Lagrangeovským přístupem a rotaci částic pomocí
Eulerových kinematických rovnic. Tyto rovnice jsou výpočetně náročnější než standardní
Lagrangeovský přístup metodou efektivního průměru, ale poskytují přesnější výsledky.
Hlavním cílem této práce bylo implementovat metodu ELER do Lattice Boltzmannov metod
(LBM), což je mezoskopický přístup modelující proudění kapalin na základě statistického
chování částic tekutiny. Tento integrovaný rámec byl následně aplikován na simulaci
transportu vláknitých částic v realistickém anatomickém modelu lidských dýchacích cest,
sahajícím až do sedmé generace větvení. Výsledky potvrdily, že použitá metoda ELER
poskytuje vyšší přesnost při predikci depozičních frakcí ve srovnání s přístupem založeným
na efektivním průměru, jak bylo ověřeno na základě experimentálních dat. Nicméně podrobná
analýza rotačního pohybu vláken pomocí vysokorychlostní kamery odhalila stále
přetrvávající nesrovnalosti. Jako jejich možný zdroj byly identifikovány asymetrické deformace
vláken, vzniklé během výrobního procesu, které byly pozorovány v experimentech.
Pro zohlednění těchto reálných tvarových změn byl navržen a implementován nový přístup
pro částice s nehomogenním rozložením hmotnosti. Zatímco bylo prokázáno, že tyto nehomogenity
významně ovlivňují rotační pohyb ve srovnání s homogenními částicemi, jejich
dopad na trajektorii a depozici částic byl méně výrazný. V důsledku toho bylo prokázáno,
že metoda ELER ve spojení s LBM je přesným nástrojem pro predikci klíčových depozičních
charakteristik v biomedicíně a dosahuje přesnějších výsledků než metoda efektivního
průměru. Určité nesrovnalosti v rotačním pohybu však zůstávají předmětem budoucího
výzkumu.

Klíčová slova:

Lattice Boltzmannova metoda, vlákna, Euler-Lagrange Euler-Rotation, depozice.

Termín obhajoby

21.11.2025

Výsledek obhajoby

obhájeno (práce byla úspěšně obhájena)

znamkaPznamka

Průběh obhajoby

Práce nabízí nový podrobný přístup k řešení pohybu vláknových částic v proudu vzduchu v dílčí části dýchacího traktu. Student zodpověděl položené dotazy a komentoval dílčí nesrovnalosti práce. Práce má přínos v oblasti lékařství, konkrétně distribuce léčiv v oblasti plic.

Jazyk práce

čeština

Fakulta

Ústav

Studijní program

Inženýrská mechanika (D-IME-P)

Složení komise

prof. Ing. Jiří Pospíšil, Ph.D. (předseda)
prof. Ing. Jan Vimmr, Ph.D. (člen)
doc. Ing. Jaroslav Katolický, Ph.D. (člen)
doc. Ing. Bc. Jan Fišer, Ph.D. (člen)
doc. Ing. Radek Fučík, Ph.D. (člen)
prof. Dr. Mgr. Ivan Cimrák (člen)
dr inż. Robert Straka (člen)

Ing. František Prinz nastoupil doktorandské studium v akademickém roce 2019/2020. Současně v tomto školním roce nastoupil na studium 5. ročníku Matematického inženýrství, které ukončil s výborným hodnocením. Již na počátku roku 2019 projevil zájem spolupracovat na vývoji metody lattice Boltzmann, kterou pak začal rozvíjet v rámci diplomní práce na Matematickém inženýrství a současně v rámci doktorského studia. Od roku 2019 začala úzká spolupráce s týmem profesora M. Krauseho z Karlsruhe university (KIT), a též na vývoji jejich open source kódu OpenLB.

Od počátku svého studie, jak magisterského na Matematickém inženýrství, tak doktorského na Energetickém ústavu, organizoval s německými partnery pravidelné online meetingy, v kterých si vyjasňoval nejasnosti ohledně LB. Postupně se zúčastnil několika mezinárodních seminářů LB.

Rozvoj metody LB zaměřil na dvoufázové diskrétní soustavy „kapka-plyn“, či „pevné částice-plyn“. Po rozsáhlých úvodních studiích problematiky LB se zaměřil na zcela novou oblast, tj. na soustavu „vláknové částice-plyn“ s aplikací na transport a depozici těchto částic v dýchacím systému člověka. Toto téma je zcela nové a jeho řešení představuje významný posun v aplikaci lattice Boltzmann.

Od roku 2019 se zapojil do řešení grantového projektu GAČR prof. Jíchy „Influence of infants and children lungs development on flow field and aerosol deposition–computational modelling and experimental validation“ a stal se také součástí týmu projektu VUT KInG „Droplet dispersion in low turbulent airflow and droplet impact on the wire mesh“. V letech 20222 až 2024 byl spoluřešitel Polsko-českého projektu Lead Agency „Holistic approach to Rotating Packed Bed (RPB) absorption process with the use of 3D CFD, visual studies, and mass transfer experiments“. Od roku 2025 se zapojil do řešení mezinárodního Česko-Polsko-Německého projektu programu Weave a současně projektu Theta II, kde se zabývá řešením dvoufázového proudění s varem na rotujícím disku. V projektech je jeho úkolem rozvoj a aplikace metod počítačových simulací pro vícefázové soustavy s použitím jak klasického CFD, tak tam, kde to složitost modelu umožní, i metodou lattice Boltzmann.

Závěrem chci zdůraznit, že úkolem doktoranda nebyl vývoj zcela nového přístupu k LB metodě, ale zaměřit se na možnou praktickou aplikaci metody s dílčími vývojovými prvky ve spolupráci s německým partnerem.

Známka navržená vedoucím: A

viz. posudek v PDF
Soubor vložený oponentem Velikost
Posudek oponenta [.pdf] 97,90 kB

viz. posudek v PDF
Soubor vložený oponentem Velikost
Posudek oponenta [.pdf] 232,84 kB

Odpovědnost: Mgr. et Mgr. Hana Odstrčilová