Flow transport and gas mixing during invasive high frequency oscillatory ventilation

The effect of budesonide delivered by high-frequency oscillatory ventilation on acute inflammatory response in severe lung injury in adult rabbits

The inflammation present in acute respiratory distress syndrome (ARDS) and thereby associated injury to the alveolar-capillary membrane and pulmonary surfactant can potentiate respiratory failure. Even considering the high mortality rate of severe ARDS, glucocorticoids appear to be a reasonable treatment option along with an appropriate route of delivery to the distal lung. This study aimed to investigate the effect of budesonide therapy delivered intratracheally by high-frequency oscillatory ventilation (HFOV) on lung function and inflammation in severe ARDS. Adult New Zealand rabbits with respiratory failure (P/F<13.3 kPa) induced by intratracheal instillation of hydrochloric acid (HCl, 3 ml/kg, pH 1.5) followed by high tidal ventilation (VT 20 ml/kg) to mimic ventilator-induced lung injury (VILI) were treated with intratracheal bolus of budesonide (0.25 mg/kg, Pulmicort) delivered by HFOV (frequency 8 Hz, MAP 1 kPa, deltaP 0.9 kPa). Saline instead of HCl without VILI with HFOV delivered air bolus instead of therapy served as healthy control. All animals were subjected to lung-protective ventilation for 4 h, and respiratory parameters were monitored regularly. Postmortem, lung injury, wet-to-dry weight ratio, leukocyte shifts, and levels of cytokines in plasma and lung were evaluated. Budesonide therapy improved the lung function (P/F ratio, oxygenation index, and compliance), decreased the cytokine levels, reduced lung edema and neutrophils influx into the lung, and improved lung architecture in interstitial congestion, hyaline membrane, and atelectasis formation compared to untreated animals. This study indicates that HFOV delivered budesonide effectively ameliorated respiratory function, and attenuated acid-induced lung injury in a rabbit model of severe ARDS.

Endotracheal tube, by the venturi effect, reduces the efficacy of increasing inlet pressure in improving pendelluft

In mechanically ventilated severe acute respiratory distress syndrome patients, spontaneous inspiratory effort generates more negative pressure in the dorsal lung than in the ventral lung. The airflow caused by this pressure difference is called pendelluft, which is a possible mechanisms of patient self-inflicted lung injury. This study aimed to use computer simulation to understand how the endotracheal tube and insufficient ventilatory support contribute to pendelluft. We established two models. In the invasive model, an endotracheal tube was connected to the tracheobronchial tree with 34 outlets grouped into six locations: the right and left upper, lower, and middle lobes. In the non-invasive model, the upper airway, including the glottis, was connected to the tracheobronchial tree. To recreate the inspiratory effort of acute respiratory distress syndrome patients, the lower lobe pressure was set at -13 cmH2O, while the upper and middle lobe pressure was set at -6.4 cmH2O. The inlet pressure was set from 10 to 30 cmH2O to recreate ventilatory support. Using the finite volume method, the total flow rates through each model and toward each lobe were calculated. The invasive model had half the total flow rate of the non-invasive model (1.92 L/s versus 3.73 L/s under 10 cmH2O, respectively). More pendelluft (gas flow into the model from the outlets) was observed in the invasive model than in the non-invasive model. The inlet pressure increase from 10 to 30 cmH2O decreased pendelluft by 11% and 29% in the invasive and non-invasive models, respectively. In the invasive model, a faster jet flowed from the tip of the endotracheal tube toward the lower lobes, consequently entraining gas from the upper and middle lobes. Increasing ventilatory support intensifies the jet from the endotracheal tube, causing a venturi effect at the bifurcation in the tracheobronchial tree. Clinically acceptable ventilatory support cannot completely prevent pendelluft.

Ventilation‐induced epithelial injury drives biological onset of lung trauma in vitro and is mitigated with prophylactic anti‐inflammatory therapeutics

Mortality rates among patients suffering from acute respiratory failure remain perplexingly high despite the maintenance of blood oxygen homeostasis during ventilatory support. The biotrauma hypothesis advocates that mechanical forces from invasive ventilation trigger immunological mediators that spread systemically. Yet, how these forces elicit an immune response remains unclear. Here, a biomimetic in vitro three‐dimensional (3D) upper airways model allows to recapitulate lung injury and immune responses induced during invasive mechanical ventilation in neonates. Under such ventilatory support, flow‐induced stresses injure the bronchial epithelium of the intubated airways model and directly modulate epithelial cell inflammatory cytokine secretion associated with pulmonary injury. Fluorescence microscopy and biochemical analyses reveal site‐specific susceptibility to epithelial erosion in airways from jet‐flow impaction and are linked to increases in cell apoptosis and modulated secretions of cytokines IL‐6, ‐8, and ‐10. In an effort to mitigate the onset of biotrauma, prophylactic pharmacological treatment with Montelukast, a leukotriene receptor antagonist, reduces apoptosis and pro‐inflammatory signaling during invasive ventilation of the in vitro model. This 3D airway platform points to a previously overlooked origin of lung injury and showcases translational opportunities in preclinical pulmonary research toward protective therapies and improved protocols for patient care.

Independent Lung Ventilation-Experimental Studies on a 3D Printed Respiratory Tract Model

Independent lung ventilation (ILV) is a life-saving procedure in unilateral pulmonary pathologies. ILV is underused in clinical practice, mostly due to the technically demanding placement of a double lumen endotracheal tube (ETT). Moreover, the determination of ventilation parameters for each lung in vivo is limited. In recent years, the development of 3D printing techniques enabled the production of highly accurate physical models of anatomical structures used for in vitro research, considering the high risk of in vivo studies. The purpose of this study was to assess the influence of double-lumen ETT on the gas transport and mixing in the anatomically accurate 3D-printed model of the bronchial tree, with lung lobes of different compliances, using various ventilation modes. The bronchial tree was obtained from Respiratory Drug Delivery (RDD Online, Richmond, VA, USA), processed and printed by a dual extruder FFF 3D printer. The test system was also composed of left side double-lumen endotracheal tube, Siemens Test Lung 190 and anesthetic breathing bag (as lobes). Pressure and flow measurements were taken at the outlets of the secondary bronchus. The measured resistance increased six times in the presence of double-lumen ETT. Differences between the flow distribution to the less and more compliant lobe were more significant for the airways with double-lumen ETT. The ability to predict the actual flow distribution in model airways is necessary to conduct effective ILV in clinical conditions.

The effect of the volumetric flow rate and endotracheal tube diameter on the pressure distribution in human airways

The comprehension of the fluid flow in the upper airways is of paramount importance when treating patients under clinical conditions that demand mechanical ventilation. Barotrauma and overdistension are related to undesirable pressures and might be responsible for morbidity and mortality. In the current work we use computational fluid dynamics to investigate the pressure field in the upper respiratory airways. We performed a set of simulations varying the volumetric flow rate of mechanical ventilators and we have shown that the pressure profile can be calculated by means of the volumetric flow rate in accordance with a mathematical expression given by P_av=aV˙², where P_av is the average pressure at selected sections of the upper airways and V˙ is the volumetric flow rate. Numerical findings provide evidence that the constant a varies with the location of the plane in the upper airways. We also show that some particular diameters of endotracheal tubes (ETT) must be used with care for a given range of volumetric flow rates. Overall, we document an important relationship among pressure, volumetric flow rate and selected internal diameters from ETT.

Effect of Non-Newtonian Dynamics on the Clearance of Mucus From Bifurcating Lung Airway Models

Mucus hypersecretion is a common pathophysiological manifestation of several obstructive airway diseases in which the mucociliary clearance is impaired, and the airflow generated by a cough or a forced expiratory maneuver called the huff is primarily responsible for clearing mucus. This airflow driven clearance of mucus is a complex process that is affected by the mucus rheology, airflow rate, airway geometry, and gravity. This study examines the role of mucus rheology in the transport and distribution of mucus in idealized 3D airway geometries. The complex air-mucus interface was tracked by the volume-of-fluid (VOF) model, and the turbulence in the core airflow was modeled using the k-ω shear stress transport (SST) model. Mucus was modeled as a shear-thinning liquid by using a power-law model. The computational model was validated using in vitro experimental data available in the literature. Gravity-dominated eccentric core-annular flow was observed with the core biased toward the outer wall in the inclined daughter branches of the bifurcation models, which transitions into concentric core-annular flow in the trachea. The increase in tangential shear at the interface due to the secondary flow structures developed in the flow divider location resulted in a region of enhanced mucus clearance with reduced mucus layer thickness. Secondary flow developed due to the curvature in the airway geometry resulted in a local redistribution of mucus that reduced the eccentricity. The accumulation of mucus around the carinal ridges and the regions with reduced clearance are sites with the potential for microbial growth.

The Effectiveness of An Averaged Airway Model in Predicting the Airflow and Particle Transport Through the Airway

Background: In this study, we proposed an averaged airway model design based on four healthy subjects and numerically evaluated its effectiveness for predicting the airflow and particle transport through an airway. Methods: Direct-averaged models of the conducting airways of four subjects were restored by averaging the three-dimensional (3D) skeletons of four healthy airways, which were calculated using an inverse 3D thinning algorithm. We simulated the airflow and particle transport in the individual and the averaged airway models using computational fluid dynamics. Results: The bifurcation geometry differs even among healthy subjects, but the averaged model retains the typical geometrical characteristics of the airways. The Reynolds number of the averaged model varied within the range found in the individual subject models, and the averaged model had similar inspiratory flow characteristics as the individual subject models. The deposition fractions at almost all individual lobes ranged within the variation observed in the subjects, however, the deposition fraction was higher in only one lobe. The deposition distribution at the main bifurcation point differed among the healthy subjects, but the characteristics of the averaged model fell within the variation observed in the individual subject models. On the contrary, the deposition fraction of the averaged model was higher than that of the average of the individual subject models and deviated from the range observed in the subject models. Conclusion: These results indicate that the direct-averaged model may be useful for predicting the individual airflow and particle transport on a macroscopic scale.

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.

Fluid flow and particle transport in mechanically ventilated airways. Part I. Fluid flow structures

A large eddy simulation-based computational study of fluid flow and particle transport in upper tracheobronchial airways is carried out to investigate the effect of ventilation parameters on pulmonary fluid flow. Respiratory waveforms commonly used by commercial mechanical ventilators are used to study the effect of ventilation parameters and ventilation circuit on pulmonary fluid dynamics. A companion paper (Alzahrany et al. in Med Biol Eng Comput, 2014) reports our findings on the effect of the ventilation parameters and circuit on particle transport and aerosolized drug delivery. The endotracheal tube (ETT) was found to be an important geometric feature and resulted in a fluid jet that caused an increase in turbulence and created a recirculation zone with high wall shear stress in the main bronchi. Stronger turbulence was found in lower airways than would be found under normal breathing conditions due to the presence of the jet caused by the ETT. The pressure-controlled sinusoidal waveform induced the lowest wall shear stress on the airways wall.

Fluid flow and particle transport in mechanically ventilated airways. Part II: particle transport

The flow mechanisms that play a role on aerosol deposition were identified and presented in a companion paper (Timothy et al. in Med Biol Eng Comput. doi: 10.1007/s11517-015-1407-3 , 2015). In the current paper, the effects of invasive conventional mechanical ventilation waveforms and endotracheal tube (ETT) on the aerosol transport were investigated. In addition to the enhanced deposition seen at the carinas of the airway bifurcations, enhanced deposition was also seen in the right main bronchus due to impaction and turbulent dispersion resulting from the fluid structures created by jet caused by the ETT. The orientation of the ETT toward right bronchus resulted in a substantial deposition inside right lung compared to left lung. The deposition inside right lung was ~12-fold higher than left lung for all considered cases, except for the case of using pressure-controlled sinusoidal waveform where a reduction of this ratio by ~50 % was found. The total deposition during pressure constant, volume ramp, and ascending ramp waveforms was similar and ~1.44 times higher than deposition fraction when using pressure sinusoidal waveform. Varying respiratory waveform demonstrated a significant role on the deposition enhancement factors and give evidence of drug aerosol concentrations in key deposition sites, which may be significant for drugs with negative side effects in high concentrations. These observations are thought to be important for ventilation treatment strategy.

A biomechanical model of pendelluft induced lung injury

Lung ventilation using high frequency oscillatory techniques have been documented to attain adequate gas exchange through various gas transport mechanisms. Among them, the pendelluft flow is considered one of the most crucial mechanisms. In this work, we computationally investigate the induction of abnormal mechanical stresses and a regionally trapped volume of gas due to pendelluft flow. Large eddy simulation was used to model the turbulence in an upper tracheobronchial lung geometry that was derived from CT scans. The pendelluft flow was captured by modeling physiological boundary conditions at the truncated level of the lung model that is sensitive to the coupled resistance and compliance of individual patients. The flow-volume and volume-pressure loops are characterized by irregular shapes and suggest abnormal regional lung ventilation. Incomplete loops were observed indicating gas trapping in these regions signifying a potential for local injury due to incomplete ventilation from a residual volume build-up at the end of the expiration phase. In addition, the gas exchange between units was observed to create a velocity gradient causing a region of high wall shear stress surrounding the carina ridges. The recurrence of the pendelluft flow could cause a rupture to the lung epithelium layer. The trapped gas and wall shear stress were observed to amplify with increasing compliance asymmetry and ventilator operating frequency. In general, despite the significant contribution of the pendelluft flow to the gas exchange augmentation there exists significant risks of localized lung injury, phenomena we describe as pendelluft induced lung injury or PILI.

Effect of mechanical ventilation waveforms on airway wall shear

Better understanding of airway wall shear stress/strain rate is very important in order to prevent inflammation in patients undergoing mechanical ventilation due to respiratory problems in intensive-care medicine. The objective of this study was to investigate the role of mechanical ventilation waveforms on airway wall shear/strain rate using computational fluid dynamics analysis. Six different waveforms were considered to investigate the airway wall shear stress (WSS) from fluid dynamics analysis for the airway geometry of two-to-three generations. The simulation results showed that Original with Sine Inhale Waveform (OSIW) produced the highest WSS value and the Near True Sine Waveform produced the lowest WSS value. Also, the Original with Sine Inhale Waveform and the Short Sine Inhale with Long Sine Exhale Waveform (SSILSEW) produced a higher shear strain rate in comparison to the Original Waveform (OW). These results, combined with optimization, suggest that it is possible to develop a set of mechanical ventilation waveform strategies to avoid inflammation in the lung.

The role of coupled resistance-compliance in upper tracheobronchial airways under high frequency oscillatory ventilation

A large eddy simulation (LES) based computational fluid dynamics (CFD) study was conducted to investigate lung lobar ventilation and gas exchange under high frequency oscillatory ventilation conditions. Time-dependent pressure coupled with the airways resistance and compliance (R&C) were imposed as boundary conditions (BCs) in the upper tracheobronchial tree of patient-specific lung geometry. The flow distribution in the left and right lungs demonstrated significant variations compared to the case in which traditional BCs based on mass flow rate fractions was used and is in agreement with the in vivo data available in the literature. The gas transport due to the pendelluft mechanism was captured in the different lung lobes and units. The computed pendelluft elapsed time was consistent with available physiological data. In contrast to in vivo studies, our simulations were able to predict the volume associated with the pendelluft elapsed time at different frequencies. Significant differences in coaxial counter flow and flow structures were observed between different BCs. The consistency of the results with the physiological in vivo data indicates that computations with coupled R&C BCs provide a suitable alternative tool for understanding the gas transport, diagnosing lung pathway disease severity, and optimizing ventilation management techniques.