Computational Analysis of Microbubble Flows in Bifurcating Airways: Role of Gravity, Inertia, and Surface Tension

Abstract

Although mechanical ventilation is a life-saving therapy for patients with severe lung disorders, the microbubble flows generated during ventilation generate hydrodynamic stresses, including pressure and shear stress gradients, which damage the pulmonary epithelium. In this study, we used computational fluid dynamics to investigate how gravity, inertia, and surface tension influence both microbubble flow patterns in bifurcating airways and the magnitude/distribution of hydrodynamic stresses on the airway wall. Direct interface tracking and finite element techniques were used to simulate bubble propagation in a two-dimensional (2D) liquid-filled bifurcating airway. Computational solutions of the full incompressible Navier-Stokes equation were used to investigate how inertia, gravity, and surface tension forces as characterized by the Reynolds (Re), Bond (Bo), and Capillary (Ca) numbers influence pressure and shear stress gradients at the airway wall. Gravity had a significant impact on flow patterns and hydrodynamic stress magnitudes where Bo > 1 led to dramatic changes in bubble shape and increased pressure and shear stress gradients in the upper daughter airway. Interestingly, increased pressure gradients near the bifurcation point (i.e., carina) were only elevated during asymmetric bubble splitting. Although changes in pressure gradient magnitudes were generally more sensitive to Ca, under large Re conditions, both Re and Ca significantly altered the pressure gradient magnitude. We conclude that inertia, gravity, and surface tension can all have a significant impact on microbubble flow patterns and hydrodynamic stresses in bifurcating airways.

Quantification of nasal airflow resistance in English bulldogs using computed tomography and computational fluid dynamics

“…For this analysis, the airways were filled with an incompressible, Newtonian fluid with constant density, and viscosity that mimicked the properties of air at environmental conditions. Iterations of the Navier‐Stokes equation for incompressible flow were used to solve the momentum balance and continuity equations governing airflow where ρ is the density, ∂ is a partial derivative, is the fluid velocity vector, p ' is the pressure, μ is the dynamic viscosity, ∇ is the spatial gradient operator, and g is the acceleration due to gravity. Boundary conditions of no flow on airway walls, zero reference pressure at the nares, and an airflow rate of 0.5 l/s at the most distal airway section near the nasopharynx were used to solve Equation in dimensional formats.…”

Multi-scale modeling of an upper respiratory airway: Effect of mucosal adhesion on Eustachian tube function in young children

“…mucolytics and surfactants. Therefore, future studies could use advanced fluid-structure interaction (FSI) and free-surface CFD techniques 34 to directly simulate the dynamics of air-liquid interface movement during ET opening. In addition, molecular scale interactions between muco-adhesive glycoproteins may involve non-linear forcedisplacement relationships and plastic deformation.…”

Multi‐scale finite element modeling of Eustachian tube function: influence of mucosal adhesion

“…In addition to increases in adhesion because of mucus glycoprotein interactions, modeled as adhesive springs in this study, other biophysical changes including changes in mucosal fluid viscosity and surface tension may also alter adhesion dynamics. Therefore, future studies could use sophisticated free‐surface CFD techniques to model air–liquid flows within the ET.…”