Modeling Flow in a Compromised Pediatric Airway Breathing Air and Heliox:

Buckling critical pressures in collapsible tubes relevant for biomedical flows

The behaviour of collapsed or stenotic vessels in the human body can be studied by means of simplified geometries like a collapsible tube. The objective of this work is to determine the value of the buckling critical pressure of a collapsible tube by employing Landau's theory of phase transition. The methodology is based on the implementation of an experimentally validated 3D numerical model of a collapsible tube. The buckling critical pressure is estimated for different values of geometric parameters of the system by treating the relation between the intramural pressure and the area of the central cross-section as the order parameter function of the system. The results show the dependence of the buckling critical pressures on the geometric parameters of a collapsible tube. General non-dimensional equations for the buckling critical pressures are derived. The advantage of this method is that it does not require any geometric assumption, but it is solely based on the observation that the buckling of a collapsible tube can be treated as a second-order phase transition. The investigated geometric and elastic parameters are sensible for biomedical application, with particular interest to the study of the bronchial tree under pathophysiological conditions like asthma.

Predicting critical closing pressure in children with obstructive sleep apnea using fluid-structure interaction

Surgical treatment of obstructive sleep apnea (OSA) in children requires knowledge of upper airway dynamics, including the closing pressure (Pcrit), a measure of airway collapsibility. We applied a flow-structure interaction (FSI) computational model to estimate Pcrit in patient-specific upper airway models obtained from magnetic resonance imaging (MRI) scans. We sought to examine the agreement between measured and estimated Pcrit from FSI models in children with Down syndrome. We hypothesized that the estimated Pcrit would accurately reflect measured Pcrit during sleep and therefore reflect the severity of OSA as measured by the obstructive apnea-hypopnea index (AHI). All participants (n = 41) underwent polysomnography and sedated sleep MRI scans. We used Bland-Altman plots to examine the agreement between measured and estimated Pcrit. We determined associations between estimated Pcrit and OSA severity, as measured by AHI, using regression models. The agreement between passive and estimated Pcrit showed a fixed bias of -1.31 [confidence interval (CI) = -2.78, 0.15] and a nonsignificant proportional bias. A weaker agreement with active Pcrit was observed. A model including AHI, gender, an interaction term for AHI, and gender and neck circumference explained the largest variation (R2 = 0.61) in the relationship between AHI and estimated Pcrit (P < 0.0001). Overlap between the areas of the airway with the lowest stiffness, and areas of collapse on dynamic MRI, was 77.4 ± 30% for the nasopharyngeal region and 78.6 ± 33% for the retroglossal region. The agreement between measured and estimated Pcrit and the significant association with AHI supports the validity of Pcrit estimates from the FSI model.NEW & NOTEWORTHY We present a noninvasive method for estimating critical closing pressure (Pcrit) using fluid-structure interaction (FSI) simulations and magnetic resonance imaging (MRI) scans in patients with obstructive sleep apnea (OSA). We used patient-specific stiffness measures in our FSI model to account for any individual variability in the elasticity of soft tissues surrounding the upper airway. We validated this model by measuring the degree of agreement between measured and estimated Pcrit.

Effect of inhaled gas density on the pendelluft-induced lung injury

Helium, sulfur hexafluoride-oxygen, and air were modeled to examine the role of the gas density on the pendelluft-induced lung injury (PILI) under high frequency oscillatory ventilation (HFOV). Large eddy simulation coupled with physiological resistance-compliance boundary conditions was applied to capture pendelluft-induced gas entrapment and mechanical stresses in an image-based human lung model. The flow characteristics were strongly dependent on the inspired gas density. The flow partitioning, globally between the left and right lung and locally between adjacent units branches, was significantly affected by the density of inhaled gas and was more balanced when inspiring lighter gas. The incomplete loops of flow-volume and volume-pressure curves were significantly influenced by the variations of the flow redistribution, resistance, and turbulence associated with the pendelluft mechanism. Inhaling light gas reduced the entrapped gas volume and mechanical stress surrounding carina ridges signifying the important role of inhaled gas properties on PILI. In general, lung ventilation by HFOV with a gas mixture of large amounts of Helium is thought to mitigate ventilator complications.

The effects of curvature and constriction on airflow and energy loss in pathological tracheas

This paper considers factors that play a significant role in determining inspiratory pressure and energy losses in the human trachea. Previous characterisations of pathological geometry changes have focussed on relating airway constriction and subsequent pressure loss, however many pathologies that affect the trachea cause deviation, increased curvature, constriction or a combination of these. This study investigates the effects of these measures on tracheal flow mechanics, using the compressive goitre (a thyroid gland enlargement) as an example. Computational fluid dynamics simulations were performed in airways affected by goitres (with differing geometric consequences) and a normal geometry for comparison. Realistic airways, derived from medical images, were used because idealised geometries often oversimplify the complex anatomy of the larynx and its effects on the flow. Two mechanisms, distinct from stenosis, were found to strongly affect airflow energy dissipation in the pathological tracheas. The jet emanating from the glottis displayed different impingement and breakdown patterns in pathological geometries and increased loss was associated with curvature.

Planning human upper airway surgery using computational fluid dynamics

The study advances the idea of using computational fluid dynamics in the process of planning surgical treatment modalities for patients with obstructive airway disorders. It is hypothesized that the a priori knowledge of the functional outcome of surgical intervention on the flow and airway resistance can guide the surgeon in choosing an effective surgical strategy. Computed tomography images spanning the respiratory tract of an adult patient with a combined glottic and subglottic stenosis are used to reconstruct three-dimensional geometrical models of the airway. Computational fluid dynamics is used to obtain airway flow patterns during inspiration and expiration in these models. Numerical predictions about flow velocity, pressure distribution on the airway lumen, wall shear stress, and airway resistance are obtained so that the relevance of each individual stenotic level is quantified. Four different virtual surgeries in different combinations are assessed in order to remedy the constricted airway. The virtual surgery based airway models are evaluated by comparisons with the pre-treatment flow modeling results. The predicted numerical data revealed that the removal of the constriction at the level of the vocal folds will have the most significant effect on the airway resistance. The flow simulations offer a quantitative method of evaluating the airway resistance in patients with combined glottic and subglottic stenoses. Predictions of airway resistances and other numerical calculations from different virtual surgeries give additional inputs for the surgeon, in deciding the most appropriate surgery on a case-by-case basis.

Influence of subglottic stenosis on the flow-induced vibration of a computational vocal fold model

The effect of subglottic stenosis on vocal fold vibration is investigated. An idealized stenosis is defined, parameterized, and incorporated into a two-dimensional, fully-coupled finite element model of the vocal folds and laryngeal airway. Flow-induced responses of the vocal fold model to varying severities of stenosis are compared. The model vibration was not appreciably affected by stenosis severities of up to 60% occlusion. Model vibration was altered by stenosis severities of 90% or greater, evidenced by decreased superior model displacement, glottal width amplitude, and flow rate amplitude. Predictions of vibration frequency and maximum flow declination rate were also altered by high stenosis severities. The observed changes became more pronounced with increasing stenosis severity and inlet pressure, and the trends correlated well with flow resistance calculations. Flow visualization was used to characterize subglottal flow patterns in the space between the stenosis and the vocal folds. Underlying mechanisms for the observed changes, possible implications for human voice production, and suggestions for future work are discussed.

Patterns in pharyngeal airflow associated with sleep-disordered breathing

OBJECTIVE

To establish the feasibility of a noninvasive method to identify pharyngeal airflow characteristics in sleep-disordered breathing.

METHODS

Four patients with sleep-disordered breathing who underwent surgery or used positive airway pressure devices and four normal healthy controls were studied. Three-dimensional CT imaging and computational fluid dynamics modeling with standard steady-state numerical formulation were used to characterize pharyngeal airflow behavior in normals and pre-and post-treatment in patients. Dynamic flow simulations using an unsteady approach were performed in one patient.

RESULTS

The pre-treatment pharyngeal airway below the minimum cross-sectional area obstruction site showed airflow separation. This generated recirculation airflow regions and enhanced turbulence zones where vortices developed. This interaction induced large fluctuations in airflow variables and increased aerodynamic forces acting on the pharyngeal wall. At post-treatment, for the same volumetric flow rate, airflow field instabilities vanished and airflow characteristics improved. Mean maximum airflow velocity during inspiration reduced from 18.3±5.7 m/s pre-treatment to 6.3±4.5 m/s post-treatment (P=0.002), leading to a reduction in maximum wall shear stress from 4.8±1.7 Pa pre-treatment to 0.9±1.0 Pa post-treatment (P=0.01). The airway resistance improved from 4.3±1.4 Pa/L/min at pre-treatment to 0.7±0.7 Pa/L/min at post-treatment (P=0.004). Post-treatment airflow characteristics were not different from normal controls (all P ≥ 0.39).

CONCLUSION

This study demonstrates that pharyngeal airflow variables may be derived from CT imaging and computational fluid dynamics modeling, resulting in high quality visualizations of airflow characteristics of axial velocity, static pressure, and wall shear stress in sleep-disordered breathing.

Maxillary, mandibular, and chin advancement: treatment planning based on airway anatomy in obstructive sleep apnea

Surgical correction of obstructive sleep apnea (OSA) syndrome involves understanding a number of parameters, of which the 3-dimensional airway anatomy is important. Visualization of the upper airway based on cone beam computed tomography scans and automated computer analysis is an aid in understanding normal and abnormal airway conditions and their response to surgery. The goal of surgical treatment of OSA syndrome is to enlarge the velo-oropharyngeal airway by anterior/lateral displacement of the soft tissues and musculature by maxillary, mandibular, and possibly, genioglossus advancement. Knowledge of the specific airway obstruction and characteristics based on 3-dimensional studies permits a directed surgical treatment plan that can successfully address the area or areas of airway obstruction. The end occlusal result can be improved when orthodontic treatment is combined with the surgical plan. The individual with OSA, though, is more complicated than the usual orthognathic patient, and both the medical condition and treatment length need to be judiciously managed when OSA and associated conditions are present. The perioperative management of the patient with OSA is more complex and the margin for error is reduced, and this needs to be taken into consideration and the care altered as indicated.