Matematická analýza a výpočtový algoritmus Rayleighovy-Plessetovy rovnice v okolí prudkého kolapsu bublinky.

Abstract

Rayleighova-Plessetova rovnice je běžně užívána k popisu dynamiky kulových bublin ve vodě. Avšak pozorování vysokorýchlostní kamerou ukazují, že bublina během prudkého kolapsu ztrácí svůj kulový tvar a často se rozpadá na mnoho menších bublin předtím, než se znovu zvětší nebo se úplně rozptýlí. Úkolem této práce je propojit mezeru mezi "sférickým" Rayleigh-Plessetovým přístupem a "nesférickým" chováním kolabující bubliny v podmínkách, kdy Rayleigh-Plessetova rovnice není už dále platná. Úkolu je dosaženo ve třech krocích. V prvním kroku jsme provedli analýzu tvarové stability sférického tvaru bubliny během kolapsu která vede k vývoji fyzikálního modelu obsahujícího ztrátu tvarové stability, transformaci povrchové energie a proces rozpadání bublin. Navíc odhaduje počet nových bublin a také podmínky tohoto rozpadu. Nakonec jsou teoretické výsledky včleněny do komplexního numerického kódu s důrazem na šíření zaokrouhlovacích chyb. Byly taktéž vypočteny numerické výsledky pro různé typické hydrodynamické kavitační situace. Výsledky této práce musí být dále potvrzeny experimenty.

The Rayleigh-Plesset equation is commonly used to describe the dynamics of spherical bubbles in water. However, as high-speed observations show, the bubble looses its spherical shape during the violent collapse and often ssions into a number of smaller bubbles before it rebounds again or dissolves completely. The objective of the work is to bridge the gap between the \spherical\ Rayleigh-Plesset approach and the \non-spherical\ behavior of the collapsing bubble in conditions when the Rayleigh-Plesset equation is no longer valid. The objective is achieved in 3 steps. First, the new criterion for the assessment of the reliability of the Rayleigh-Plesset equation is proposed. Second, the analysis of the shape stability of the spherical surface of the bubble during the collapse lead to the development of a physical model, which incorporates the loss of shape stability, surface energy dissipation and the ssion process. Moreover it estimates the number of bubble fragments as well as conditions when the ssion occurs. Finally, theoretical results are incorporated into the complex numerical code, with the focus on error estimation, specially propagation of rounding error. The sample numerical results for typical hydrodynamic cavitation situations were calculated. The results of this work have to be further con rmed...