Computational fluid dynamics endpoints for assessment of adenotonsillectomy outcome in obese children with obstructive sleep apnea syndrome

Airway stability in sleep apnea: Assessing continuous positive airway pressure efficiency

Obstructive Sleep Apnea Syndrome (OSAS) disrupts millions of lives with its burden of airway obstruction during sleep. Continuous Positive Airway Pressure (CPAP) therapy has been scrutinized for its biomechanical impact on the respiratory tract. This study leverages computational fluid dynamics to investigate CPAP's effects at 9 cm H2O (882.6 Pa) on the computed-tomography-based nasal-to-14-generation full respiratory tract model compared to ambient conditions, focusing on static pressure, airflow velocity, and shear stress. Our findings reveal that CPAP significantly increases static pressure, enhancing airway patency without adverse changes in airflow velocity or harmful shear stress on lung tissue, challenging prior concerns about its safety. Notably, the larynx experiences the highest shear stress due to its narrow anatomy, yet CPAP therapy overall supports airway walls against collapse. This investigation highlights CPAP's critical role in OSAS treatment, offering reassurance about its safety and efficacy. By clarifying CPAP therapy's physiological impacts, our study contributes vital insights for optimizing OSAS management strategies, affirming CPAP's benefit in maintaining open airways with minimal tissue strain.

Analysis of the aerodynamic characteristics of the upper airway in obstructive sleep apnea patients

Background/purpose

This study was designed to explore upper airway aerodynamic characteristics in individuals diagnosed with obstructive sleep apnea (OSA) and to evaluate correlations between these characteristics and other anatomical upper airway findings in these patients.

Materials and methods

This was a retrospective study of 40 OSA patients (22 male, 18 female) who were stratified into groups with mild, moderate, and severe disease based upon overnight polysomnographic (PSG) recording results. Newtom5G cone-beam CT scans (CBCT) were conducted for all patients, and the resultant images were used to reconstruct three-dimensional images of the upper airways which were used to calculate aerodynamic characteristics. Differences in these characteristics between groups were evaluated with one-way ANOVAs, while relationships between anatomical and aerodynamic characteristics were assessed through Pearson correlation analyses.

Results

The aerodynamic of the upper airway has typical characteristic in severe group. There was a significant negative correlation in severe group between resistance during inspiration (Rin) and volume (V) (r = -0.693, P = 0.013), minimum axial area (MMA) (r = -0.685, P = 0.014), and lateral dimension (LAT) (r = -0.724, P = 0.008), resistance during expiration (Rex) and LAT (r = -0.923, P < 0.001).

Conclusion

This study showed that airway resistance during inspiration and expiration are most closely associated with upper airway collapse in OSA patients, with repetitive collapse occurring during both of these breathing processes. LAT may be an important anatomical factor associated with OSA pathogenesis and treatment.

Effect of different degrees of adenoid hypertrophy on pediatric upper airway aerodynamics: a computational fluid dynamics study

To improve the diagnostic accuracy of adenoid hypertrophy (AH) in children and prevent further complications in time, it is important to study and quantify the effects of different degrees of AH on pediatric upper airway (UA) aerodynamics. In this study, based on computed tomography (CT) scans of a child with AH, UA models with different degrees of obstruction (adenoidal-nasopharyngeal (AN) ratio of 0.9, 0.8, 0.7, and 0.6) and no obstruction (AN ratio of 0.5) were constructed through virtual surgery to quantitatively analyze the aerodynamic characteristics of UA with different degrees of obstruction in terms of the peak velocity, pressure drop (△P), and maximum wall shear stress (WSS). We found that two obvious whirlpools are formed in the anterior upper part of the pediatric nasal cavity and in the oropharynx, which is caused by the sudden increase in the nasal cross-section area, resulting in local flow separation and counterflow. In addition, when the AN ratio was ≥ 0.7, the airflow velocity peaked at the protruding area in the nasopharynx, with an increase 1.1-2.7 times greater than that in the nasal valve area; the △P in the nasopharynx was significantly increased, with an increase 1.1-6.8 times greater than that in the nasal cavity; and the maximum WSS of the posterior wall of the nasopharynx was 1.1-4.4 times larger than that of the nasal cavity. The results showed that the size of the adenoid plays an important role in the patency of the pediatric UA.

Effects of vertical control on anatomic and aerodynamic characteristics of the oropharyngeal airway during premolar extraction treatment of Class II hyperdivergent nonsevere crowding malocclusion

INTRODUCTION

This study aimed to analyze the effects of premolar extraction treatment with vertical control on changes in the anatomy and aerodynamics of the oropharynx in Class II hyperdivergent malocclusion with nonsevere crowding.

METHODS

Thirty-nine patients with Class II hyperdivergent malocclusion were enrolled consecutively. All the participants underwent 4 premolar extractions. The high-pull J-hook and mini-implants were used to provide vertical control. Cone-beam computed tomography was performed before and after treatment. The participants were divided into a decreased lower vertical facial height group (n = 23) and an increased lower vertical facial height group (n = 16) on the basis of superimposition. The aerodynamic characteristics, including airway resistance (inspiration, R_in; expiration, R_ex) and maximum velocity (inspiration, Vmax_in; expiration, Vmax_ex) at inspiration and expiration, were calculated using computational fluid dynamics. Anatomic characteristics, including volume and cross-sectional area (CSA_min), were measured using the Dolphin Imaging software (Dolphin Imaging and Management Solutions, Chatsworth, Calif).

RESULTS

After treatment, the median volume and CSA_min increased by 2357 mm³ and 43 mm², respectively, and median R_in and Vmax_ex decreased by 0.15 Pa/L/min and 0.24 m×s^-1, respectively, in decreased lower vertical facial height group. In contrast, the median CSA_min decreased by 9.5 mm² in the increased lower vertical facial height group. All the changes were statistically significant (all P <0.05). Significant differences in volume, CSA_min, R_in, and Vmax_ex were observed between the 2 groups.

CONCLUSIONS

Vertical control might improve the anatomic and aerodynamic characteristics of the oropharyngeal airway during premolar extraction treatment of Class II hyperdivergent malocclusion with nonsevere crowding.

The influence of mouth opening on pharyngeal pressure loss and its underlying mechanism: A computational fluid dynamic analysis

Objective: During inspiration, mechanical energy generated from respiratory muscle produces a negative pressure gradient to fulfill enough pulmonary ventilation. The pressure loss, a surrogate for energy loss, is considered as the portion of negative pressure without converting into the kinetic energy of airflow. Mouth opening (MO) during sleep is a common symptom in patients with obstructive sleep apnoea-hypopnea syndrome (OSAHS). This study aimed to evaluate the effects of mouth opening on pharyngeal pressure loss using computational fluid dynamics (CFD) simulation. Methods: A total of four subjects who were morphologically distinct in the pharyngeal characteristics based on Friedman tongue position (FTP) grades were selected. Upper airway computed tomography (CT) scan was performed under two conditions: Mouth closing (MC) and mouth opening, in order to reconstruct the upper airway models. computational fluid dynamics was used to simulate the flow on the two different occasions: Mouth closing and mouth opening. Results: The pharyngeal jet was the typical aerodynamic feature and its formation and development were different from mouth closing to mouth opening in subjects with different Friedman tongue position grades. For FTP I with mouth closing, a pharyngeal jet gradually formed with proximity to the velopharyngeal minimum area plane (plane_Amin). Downstream the plane_Amin, the jet impingement on the pharyngeal wall resulted in the frictional loss associated with wall shear stress (WSS). A rapid luminal expansion led to flow separation and large recirculation region, corresponding to the interior flow loss. They all contributed to the pharyngeal total pressure loss. While for FTP I with mouth opening, the improved velopharyngeal constriction led to smoother flow and a lower total pressure loss. For FTP IV, the narrower the plane_Amin after mouth opening, the stronger the jet formation and its impingement on the pharyngeal wall, predicting a higher frictional loss resulted from higher WSS. Besides, a longer length of the mouth opening-associated constant constrictive segment was another important morphological factor promoting frictional loss. Conclusion: For certain OSAHS patients with higher Friedman tongue position grade, mouth opening-related stronger jet formation, more jet breakdown and stronger jet flow separation might contribute to the increased pharyngeal pressure loss. It might require compensation from more inspiratory negative static pressure that would potentially increase the severity of OSAHS.

Postoperative functional evaluation of obstructive sleep apnea syndrome by computational fluid dynamics

The aim of this study is to evaluate the efficacy of uvulopalatopharyngoplasty and the corresponding postoperative airflow. Eleven patients diagnosed with obstructive sleep apnea syndrome who complained of snoring and apnea were enrolled in this study. Computational fluid dynamics (CFD) was implemented. CFD could be accomplished in nine cases. Airflow analysis was not possible in cases with a high respiratory event index (REI) score. Before surgery, stenosis was identified in the oropharynx and epiglottic area. And the airflow velocity and pressure were found to have significantly decreased in the oropharynx postoperatively, while in the epiglottic area, those data had increased postoperatively in some cases. The velocity and pressure of the oropharynx are related to REI score. From the CFD analysis, airflow analysis is important for evaluating the apnea state. It is suggested that the postoperative function can now be predicted preoperatively.

Fluid dynamics of the upper airway in pediatric patients with severe laryngomalacia

Laryngomalacia is the top cause of pediatric laryngeal wheeze. We used computational fluid dynamics to study the inspiratory airflow dynamics in severe pediatric laryngomalacia. Computed tomography was performed on the upper airways of two infants, one with severe laryngomalacia and one with normal airway, and 3D models were reconstructed. ANSYS CFD-POST software was used to simulate airflow in these models to compare the volumetric flow rate, flow velocity, pressure, wall shear, and vortex. The volume flow rate in the laryngomalacia model was significantly reduced compared with the control model. Under inspiratory pressures, the peak flow velocity, pressure, and shear force in the control model appeared at the soft palate stenosis, while that in the laryngomalacia model appeared at the supraglottis stenosis. In both models, the maximum flow velocity and shear force increased with decreasing inspiratory pressure, while the minimum pressure decreased with decreasing inspiratory pressure. In the control model, the airflow vortex appeared anteriorly below the posterior section of the soft palate. In the laryngomalacia model, the vortex appeared anteriorly below the posterior section of the soft palate and anteriorly below the vocal folds. Our methodology provides a new mechanistic understanding of pediatric laryngomalacia.

Upper Airway Flow Dynamics in Obstructive Sleep Apnea Patients with Various Apnea-Hypopnea Index

BACKGROUND AND AIM

This study evaluates the upper airway flow characteristics, anatomical features and analyzes their correlations with AHI in patients with varied degrees of OSA severity seeking for discernments of the underlying pathophysiological profile.

MATERIALS AND METHODS

Patient-specific computational fluid dynamics models were reconstructed from high-resolution cone-beam computed tomography images for 4 OSA patients classified as minimal, mild, moderate, and severe according to AHI.

RESULTS

The parameters, minimal cross-sectional area (MCA), and the pharyngeal airway volume did not show clear correlations with the OSA severity defined according to AHI. No correlations were found between the classically defined resistance of the airway in terms of pressure drop and AHI. The flow analysis further showed that the fluid mechanisms likely to cause airway collapse are associated with the degree of narrowing in the pharyngeal airway rather than AHI. Results also suggested that some patients classified as severe OSA according to the AHI can show less susceptibility to airway collapse than patients with relatively lower AHI values and vice versa.

CONCLUSIONS

The relative contribution of anatomical and non-anatomical causes to the OSA severity can significantly vary between patients. AHI alone is inadequate to be used as a marker of the pathophysiological profile of OSA. Combining airflow analysis with AHI in diagnosing OSA severity may provide additional details about the underlying pathophysiology, subsequently improving the individualized clinical outcomes.

Upper airway flow characteristics of childhood obstructive sleep apnea-hypopnea syndrome

Revealing the structural morphology and inner flow field of the upper airway is important for understanding obstructive sleep apnea-hypopnea syndrome (OSAHS) incidence phenomena and pathological diagnosis in children. However, prior work on this topic has been focused on adults and the findings cannot be directly extrapolated to children because of different inducing factors. Therefore, this paper employs a simulation method to investigate upper airway flow characteristics of childhood OSAHS. It is found that the Reynold number changes highly throughout the whole upper airway, and the laminar assumption is no longer suitable for low Reynold number flow, which is much unlike classic fluid mechanics. Turbulent models of Standard k-ω and Spalart-Allmaras were developed prior to suggestion. The simulation is validated by experiments with an error of approximately 20%. Additionally, carried out in this analysis is the influence of adenoidal hypertrophy with different narrow levels. The cross-sectional area, flow velocity, pressure drop and volume rate will change greatly when the narrow level is above 64% of the upper airway, which can be a quantitative explanation for medical intervention if adenoid hypertrophy blocks 2/3 of the upper airway in the common clinical judgment of otorhinolaryngology. It is expected that this paper can be a meaningful instruction on OSAHS surgery plan making as well as recovery evaluation postoperatively.

A preoperative predictive study of advantages of airway changes after maxillomandibular advancement surgery using computational fluid dynamics analysis

The purpose of this study was to develop a simulation approach for predicting maxillomandibular advancement-induced airway changes using computational fluid dynamics. Eight patients with jaw deformities who underwent maxillomandibular advancement and genioglossus advancement surgery were included in this study. Computed tomography scans and rhinomanometric readings were performed both preoperatively and postoperatively. Computational fluid dynamics models were created, and airflow simulations were performed using computational fluid dynamics software; the preferable number of computational mesh points was at least 10 million cells. The results for the right and left nares, including simulation and postoperative measurements, were qualitatively consistent, and surgery reduced airflow pressure loss. Geometry prediction simulation results were qualitatively consistent with the postoperative stereolithography data and postoperative simulation results. Simulations were performed with either the right or left naris blocked, and the predicted values were similar to those found clinically. In addition, geometry prediction simulation results were qualitatively consistent with the postoperative stereolithography data and postoperative simulation results. These findings suggest that geometry prediction simulation facilitates the preoperative prediction of the postoperative structural outcome.

A New Method for Anterior Boundary Demarcation of the Nasopharynx in Three-Dimensional Analysis

ABSTRACT

Three-dimensional (3D) measurements of the upper airway have been extensively applied and researched, but the division of the airway is carried out in various ways, especially when demarcating the anterior boundary of the nasopharynx. The present study was to propose a new method based on the anatomical definition for the anterior boundary demarcation of the nasopharynx used in three-dimensional analysis. Twenty computed tomography scans (age 9.5 ± 2.5 years, 11 males, and 9 females) of head and neck were randomly selected and transferred to Materialism's interactive medical image control system 19.0 for segmentation of the nasopharynx. Precise localization of the reference points that determining the nasopharyngeal anterior and inferior boundaries was managed by recording their coordinates in the interface of the software. The area of the anterior and the inferior boundaries, and the volume of the nasopharynx were measured and repeated with a 2-week interval for the consistency test. Both the interobserver reliability as well as the intra-observer reliability were very high (intraclass correlation coefficients, 0.985-0.997). Paired t test showed no significant difference between the first and the second examinations. This new simple method proposed for demarcation of the nasopharyngeal anterior boundary was based on the innate anatomical boundary, which was statistically reliable, technically convenient, and clinically reasonable.

Effect of sleep on upper airway dynamics in obese adolescents with obstructive sleep apnea syndrome

STUDY OBJECTIVES

The biomechanical basis of obstructive sleep apnea syndrome (OSAS) may influence upper airway dynamics. In this study, we investigate dynamic changes during respiration in wakefulness and sleep in obese adolescents with and without OSAS.

METHODS

Respiratory-gated dynamic magnetic resonance imaging (MRI) at the retropalatal and retroglossal regions was performed with simultaneous measurement of SpO2 and nasal-oral mask airflow and pressure. Airway cross-sectional area (CSA) was determined using AMIRA. Percent change in CSA was calculated from five continuous tidal breaths in states of wakefulness and sleep. Mixed effects models were used to evaluate interactions between group (OSAS/control), site (retropalatal/retroglossal), and stage (wake/sleep).

RESULTS

We studied 24 children with OSAS (mean age 15.49 ± 2.00 years, mean apnea-hypopnea index [AHI] 16.53 ± 8.72 events/h) and 19 controls (mean age 14.86 ± 1.75 years, mean AHI 2.12 ± 1.69 events/h). Groups were similar in age, sex, height, weight, and BMI Z-score. Participants with OSAS had a 48.17% greater increase in percent change of airway CSA during sleep than controls (p < 0.0001), while there was no difference between groups during wakefulness (p = 0.6589). Additionally, participants with OSAS had a 48.80% increase in percent change of airway CSA during sleep as compared with wakefulness (p < 0.0001), whereas no such relationship was observed in controls (p = 0.5513).

CONCLUSIONS

This study demonstrates significant effects of sleep on upper airway dynamics in obese children with OSAS. Dynamic MRI with physiological data can potentially provide further insight into the biomechanical basis of OSAS and assist in more effective management.

Open oral cavity has little effects on upper airway aerodynamics in children with obstructive sleep apnea syndrome: A computational fluid dynamics study based on patient-specific models

Obstructive sleep apnea syndrome (OSAS) is a common disorder with recurrent pharyngeal airway collapse and sleep disruption. Recently, great progress has been made in investigating the physical mechanism of OSAS development and treatment using computational fluid dynamics (CFD). However, previous studies always neglected the oral cavity artificially in the patient's upper airway CFD model, but did not give any specific explanation. The oral cavity effect on the OSAS upper airway flow is still a matter of unclear. This paper reconstructed the patient-specific upper airway models based on the cone beam computed tomography images of ten children subjects (seven boys and three girls) and used CFD to simulate both the steady and unsteady expiration and inspiration states in the upper airway model with or without the oral cavity. A series of pressure measurement experiments based on the in vitro 1:1 scaled airway model were performed to validate the reliability of the present CFD methods. Finally, the CFD results indicate that the open oral cavity is almost a region of flow stasis with constant pressure, and both the upper airway aerodynamics with and without the oral cavity have the similar trends, with the maximum average relative difference less than 6%. The present study shows that the open oral cavity causes very little impacts on the upper airway flow of the children patients with OSAS using the nasal respiration only, and confirms the reasonability of ignoring the oral cavity for CFD simulation.

A comparative analysis of aerodynamic and anatomic characteristics of upper airway before and after mini-implant-assisted rapid maxillary expansion

INTRODUCTION

The objective of this research was to observe changes in aerodynamics and anatomic characteristics of the upper airway after mini-implants assisted rapid maxillary expansion and to evaluate the correlation between the 2 changes of the upper airway in young adults.

METHODS

Thirty consecutive patients (mean age, 23.82 ± 3.90 years; median, 24.5 years; 9 males, 21 females) were involved. Cone-beam computed tomography was taken before activation and over 3 months. Three-dimensional models of the upper airway were reconstructed on the basis of cone-beam computed tomography. The anatomic characteristics of the upper airway, including volume, area, transverse, and sagittal diameter, were measured. The aerodynamic characteristics of the upper airway were calculated on the basis of 3-dimensional models using computational fluid dynamics. The correlation between the changes in aerodynamics and anatomic characteristics of the upper airway was explored.

RESULTS

The enlargements of the volume of the total pharynx, nasopharynx, and oropharynx were found (9.99%, 20.7%, and 8.84%, respectively). The minimum cross-sectional area increased significantly (13.6%). The airway resistance (R) and maximum velocity (V_max) decreased significantly in both the inspiration and expiration phase (inspiration: R, -26.8%, V_max, -15.7%; expiration: R, -24.7%, V_max, -16.5%). The minimum wall shear stress reduced significantly only in the inspiration phase (-26.3%). The correlations between decreased R and increased volume and minimum cross-sectional area were observed.

CONCLUSIONS

Mini-implants assisted rapid maxillary expansion is an effective device for improving anatomic characteristics represented by the total volume of the upper airway and minimum cross-sectional area, which contributed to the respiratory function depending on the favorable changes of aerodynamic characteristics including resistance, velocity, and minimum wall shear stress.

Morphologic and Aerodynamic Changes of Upper Airway in Pediatric Hemifacial Microsomia Patients Undergoing Distraction Osteogenesis

Current studies on hemifacial microsomia (HFM) patients undergoing unilateral distraction osteogenesis (DO) mainly confined to description of facial morphology. This study was to investigate the effect of unilateral DO on upper airway in pediatric HFM patients using computational fluid dynamics. The investigators implemented the present retrospective study composed of 20 patients (age 9.5 ± 2.5 years, 11 males, 9 females) with moderate HFM performed unilateral DO on the ramus of mandible. Computational fluid dynamics models of the upper airway were obtained based on the computed tomography data sets which were taken before and 6 months after distraction. Morphologic and aerodynamic parameters were respectively computed and compared at peak inspiration. Paired t-tests were used to compare the differences between the before and after parameters with the significance set at P < 0.05. The authors observed that after DO, oropharynx and nasopharynx were the major expanded regions, average pressure increased in all cross-sections, average velocity decreased significantly except in the superior border of the epiglottis plane, and the airflow resistance significantly reduced in the whole upper airway (P < 0.05). The results indicate that unilateral DO may expand the constricted oropharynx and nasopharynx, obtain a relatively symmetrical airway shape and decrease the airflow resistance which consequently reduce the workload necessary for breathing and facilitate inspiration.

Computational fluid dynamics modelling of human upper airway: A review

BACKGROUND AND OBJECTIVE

Human upper airway (HUA) has been widely investigated by many researchers covering various aspects, such as the effects of geometrical parameters on the pressure, velocity and airflow characteristics. Clinically significant obstruction can develop anywhere throughout the upper airway, leading to asphyxia and death; this is where recognition and treatment are essential and lifesaving. The availability of advanced computer, either hardware or software, and rapid development in numerical method have encouraged researchers to simulate the airflow characteristics and properties of HUA by using various patient conditions at different ranges of geometry and operating conditions. Computational fluid dynamics (CFD) has emerged as an efficient alternative tool to understand the airflow of HUA and in preparing patients to undergo surgery. The main objective of this article is to review the literature that deals with the CFD approach and modeling in analyzing HUA.

METHODS

This review article discusses the experimental and computational methods in the study of HUA. The discussion includes computational fluid dynamics approach and steps involved in the modeling used to investigate the flow characteristics of HUA. From inception to May 2020, databases of PubMed, Embase, Scopus, the Cochrane Library, BioMed Central, and Web of Science have been utilized to conduct a thorough investigation of the literature. There had been no language restrictions in publication and study design of the database searches. A total of 117 articles relevant to the topic under investigation were thoroughly and critically reviewed to give a clear information about the subject. The article summarizes the review in the form of method of studying the HUA, CFD approach in HUA, and the application of CFD for predicting HUA obstacle, including the type of CFD commercial software are used in this research area.

RESULTS

This review found that the human upper airway was well studied through the application of computational fluid dynamics, which had considerably enhanced the understanding of flow in HUA. In addition, it assisted in making strategic and reasonable decision regarding the adoption of treatment methods in clinical settings. The literature suggests that most studies were related to HUA simulation that considerably focused on the aspects of fluid dynamics. However, there is a literature gap in obtaining information on the effects of fluid-structure interaction (FSI). The application of FSI in HUA is still limited in the literature; as such, this could be a potential area for future researchers. Furthermore, majority of researchers present the findings of their work through the mechanism of airflow, such as that of velocity, pressure, and shear stress. This includes the use of Navier-Stokes equation via CFD to help visualize the actual mechanism of the airflow. The above-mentioned technique expresses the turbulent kinetic energy (TKE) in its result to demonstrate the real mechanism of the airflow. Apart from that, key result such as wall shear stress (WSS) can be revealed via turbulent kinetic energy (TKE) and turbulent energy dissipation (TED), where it can be suggestive of wall injury and collapsibility tissue to the HUA.

Primary site identification in children with obstructive sleep apnea by computational fluid dynamics analysis of the upper airway

STUDY OBJECTIVES

Obstructive sleep apnea (OSA) is a respiratory disorder caused by the obstruction of the upper airway during sleep. The identification of the primary site of OSA is essential to determine treatment strategy. This study aimed to establish computational fluid dynamics (CFD) analysis for determining the clinical severity of OSA and the primary site of OSA.

METHODS

Twenty children (mean age, 6 years) were divided into OSA and control groups according to their apnea-hypopnea index. Three-dimensional airways were constructed from computed tomography data. The pharyngeal airway morphology and the pressure and velocity of the upper airway were evaluated using CFD analysis.

RESULTS

The maximum velocity and negative pressure of the upper airway in the OSA group were significantly correlated with the severity of OSA (r_s = .741, P < .001; r_s = -.653, P = .002). A velocity higher than 12 m/s indicated the primary site of OSA. In addition, we found that the primary site of OSA is not necessarily the same as the collapsible conduit site.

CONCLUSIONS

CFD analysis allows both the evaluation of the disease severity of OSA and the identification of the primary site of OSA in children. The primary site of OSA is not necessarily the same as the collapsible conduit site; therefore, CFD analysis can be used to identify the appropriate intervention for treating OSA.

Numerical simulation of nasal airflows and thermal air modification in newborns

Warming, filtering, and humidification of inspired air are major functions of the upper airway, which can be negatively altered by local disorders or surgical interventions. These functions have not been described in neonates because of ethical and technical problems difficult to solve. Numerical simulations can get around these limitations. The objective of this study was to analyze physiological nasal airflow and thermal distribution using computational fluid dynamics (CFD) techniques in neonates. CT imaging of neonates was collected from the Pediatric Radiology Department of our center. CFD has been used to simulate nasal airflow numerically, with ambient air set at 19 °C, following the recommendations for a neonate's bedroom. Thermal distribution within the nasal cavity was analyzed and coupled with airflow patterns over complete respiratory cycles. Sixteen patients have been included in the study. During inspiration, important air warming is noticed in the first centimeter of the nasal cavity (+ 8 °C at the anterior end of the inferior turbinate). During the expiration phase, the temperature decreases slightly (- 3 °C) between the pharynx and the nostrils. A model with asymmetric nasal fossae showed that nasal obstruction leads to decreased airflow and abnormally high temperatures in the obstructed side (+ 2 °C at the nasal valve, + 4 °C at the choana). According to our results, the nasal valve area is of crucial importance in air warming in neonates, when ambient air is 19 °C, since about 70% of air warming is performed in this area. When needed, surgical interventions should respect the anatomy of this zone and restore normal airflows and warming. Graphical abstract .

Evaluation of human obstructive sleep apnea using computational fluid dynamics

Obstructive sleep apnea (OSA) severity might be correlated to the flow characteristics of the upper airways. We aimed to investigate the severity of OSA based on 3D models constructed from CT scans coupled with computational fluid dynamics (CFD) simulations. The CT scans of seven adult patients diagnosed with OSA were used to reconstruct the 3D models of the upper airways and CFD modeling and analyses were performed. Results from the fluid simulations were compared with the apnea-hypopnea index. Here we show a correlation between a CFD-based parameter, the adjusted pressure coefficient (Cp*), and the respective apnea-hypopnea index (Pearson's r = 0.91, p = 0.004), which suggests that the anatomical-based model coupled with CFD could provide functional and localized information for different regions of the upper airways.

From CPAP to tailored therapy for obstructive sleep Apnoea

Obstructive Sleep Apnoea (OSA) is a common sleep disorder that is associated with daytime symptoms and a range of comorbidity and mortality. Continuous Positive Airway Pressure (CPAP) therapy is highly efficacious at preventing OSA when in use and has long been the standard treatment for newly diagnosed patients. However, CPAP therapy has well recognised limitations in real world effectiveness due to issues with patient acceptance and suboptimal usage. There is a clear need to enhance OSA treatment strategies and options. Although there are a range of alternative treatments (e.g. weight loss, oral appliances, positional devices, surgery, and emerging therapies such as sedatives and oxygen), generally there are individual differences in efficacy and often OSA will not be completely eliminated. There is increasing recognition that OSA is a heterogeneous disorder in terms of risk factors, clinical presentation, pathophysiology and comorbidity. Better characterisation of OSA heterogeneity will enable tailored approaches to therapy to ensure treatment effectiveness. Tools to elucidate individual anatomical and pathophysiological phenotypes in clinical practice are receiving attention. Additionally, recognising patient preferences, treatment enhancement strategies and broader assessment of treatment effectiveness are part of tailoring therapy at the individual level. This review provides a narrative of current treatment approaches and limitations and the future potential for individual tailoring to enhance treatment effectiveness.

A novel method to generate dynamic boundary conditions for airway CFD by mapping upper airway movement with non-rigid registration of dynamic and static MRI

Computational fluid dynamics (CFD) simulations of airflow in the human airways have the potential to provide a great deal of information that can aid clinicians in case management and surgical decision making, such as airway resistance, energy expenditure, airflow distribution, heat and moisture transfer, and particle deposition, as well as the change in each of these due to surgical interventions. However, the clinical relevance of CFD simulations has been limited to date, as previous models either did not incorporate neuromuscular motion or any motion at all. Many common airway pathologies, such as obstructive sleep apnea (OSA) and tracheomalacia, involve large movements of the structures surrounding the airway, such as the tongue and soft palate. Airway wall motion may be due to many factors including neuromuscular motion, internal aerodynamic forces, and external forces such as gravity. Therefore, to realistically model these airway diseases, a method is required to derive the airway wall motion, whatever the cause, and apply it as a boundary condition to CFD simulations. This paper presents and validates a novel method of capturing in vivo motion of airway walls from magnetic resonance images with high spatiotemporal resolution, through a novel combination of non-rigid image, surface, and surface-normal-vector registration. Coupled with image-synchronous pneumotachography, this technique provides the necessary boundary conditions for dynamic CFD simulations of breathing, allowing the effect of the airway's complex motion to be calculated for the first time, in both normal subjects and those with conditions such as OSA.

Effect of Nasal Valve Shape on Downstream Volume, Airflow, and Pressure Drop: Importance of the Nasal Valve Revisited

OBJECTIVES

The relative importance of the nasal valve relative to the remainder of the nasal airway remains unknown. The goal of this article was to objectively measure the shape of the nasal inlet and its effect on downstream airflow and nasal cavity volume using a physical model and a physiologic flow model.

METHODS

A patient who had isolated nasal valve surgery and had pre- and postoperative computed tomography scans available for analysis was studied. Nasal inlet shape measurements, computational fluid dynamics, and nasal volume analysis were performed using the computed tomography data. In addition, a physical model was used to determine the effect of nasal obstruction on downstream soft tissue.

RESULTS

The postoperative shape of the nasal inlet was improved in terms of length and degree of tortuosity. Whereas the operated-on region at the nasal inlet showed an only 25% increase in cross-sectional area postoperatively, downstream nonoperated sites in the nasal cavity revealed increases in area ranging from 33% to 51%. Computational fluid dynamics analysis showed that airway resistance decreased by 42%, and pressure drop was reduced by 43%. Intraluminal mucosal expansion was found with nasal obstruction in the physical model.

CONCLUSION

By decreasing the degree of bending and length at the nasal valve, inspiratory downstream nonoperated sites of the nasal cavity showed improvement in volume and airflow, suggesting that the nasal valve could play an oversized role in modulating the aerodynamics of the airway. This was confirmed with the physical model of nasal obstruction on downstream mucosa.

Upper airway reconstruction using long-range optical coherence tomography: Effects of airway curvature on airflow resistance

OBJECTIVES

Adenotonsillectomy (AT) is commonly used to treat upper airway obstruction in children, but selection of patients who will benefit most from AT is challenging. The need for diagnostic evaluation tools without sedation, radiation, or high costs has motivated the development of long-range optical coherence tomography (LR-OCT), providing real-time cross-sectional airway imaging during endoscopy. Since the endoscope channel location is not tracked in conventional LR-OCT, airway curvature must be estimated and may affect predicted airway resistance. The study objective was to assess effects of three realistic airway curvatures on predicted airway resistance using computational fluid dynamics (CFD) in LR-OCT reconstructions of the upper airways of pediatric patients, before and after AT.

METHODS

Eight subjects (five males, three females, aged 4-9 years) were imaged using LR-OCT before and after AT during sedated endoscopy. Three-dimensional (3D) airway reconstructions included three airway curvatures. Steady-state, inspiratory airflow simulations were conducted under laminar conditions, along with turbulent simulations for one subject using the k-ω turbulence model. Airway resistance (pressure drop/flow) was compared using two-tailed Wilcoxon signed rank tests.

RESULTS

Regardless of the airway curvatures, CFD findings corroborate a surgical end-goal with computed post-operative airway resistance significantly less than pre-operative (P < 0.01). The individual resistances did not vary significantly for different airway curvatures (P > 0.25). Resistances computed using turbulent simulations differed from laminar results by less than ∼5%.

CONCLUSIONS

The results suggest that reconstruction of the upper airways from LR-OCT imaging data may not need to account for airway curvature to be predictive of surgical effects on airway resistance. Lasers Surg. Med. 51:150-160, 2019. © 2018 Wiley Periodicals, Inc.

Modeling congenital nasal pyriform aperture stenosis using computational fluid dynamics

OBJECTIVES

Congenital nasal pyriform aperture stenosis (CNPAS) is a rare cause of airway obstruction in the neonate. Computational airway modeling has not been done in neonates and young infants to understand the impact of stenosis on functional nasal airflow. In this study, we 1) applied computational fluid dynamics (CFD) model to the airway of a neonate with CNPAS and 2) compare airflow dynamics of a normal and CNPAS airway.

METHODS

Three-dimensional models of the nasal airway of a normal neonate and a neonate with CNPAS were created using computed tomography scans of the facial bones. Measured anatomic parameters included volume, surface area, and cross-sectional area. CFD simulation was then performed. Simulated flow parameters included pressure, average velocity, and resistance.

RESULTS

The neonate with CNPAS had a lesser volume (2.74 cm³ vs. 4.50 cm³) and surface area (18.8 cm² vs. 45.5 cm²) than the normal airway. The CNPAS airway had a lesser bilateral cross-sectional area and average cross-sectional velocity throughout the length of the model. While there is a large pressure drop in the normal airway immediately after the entry point, the pressure drop in the CNPAS airway occurs more posteriorly. The total nasal resistance was approximately eight-fold greater in the CNPAS airway than the normal.

CONCLUSIONS

CFD analysis can be performed on airways of neonates with nasal obstruction, such as in CNPAS. A CFD model may help characterize severity of airway obstruction as it can predict the three-dimensional pattern of airflow. Determining the role of CFD in clinical management of CNPAS requires further investigation.

Biomechanics of the soft-palate in sleep apnea patients with polycystic ovarian syndrome

Highly compliant tissue supporting the pharynx and low muscle tone enhance the possibility of upper airway occlusion in children with obstructive sleep apnea (OSA). The present study describes subject-specific computational modeling of flow-induced velopharyngeal narrowing in a female child with polycystic ovarian syndrome (PCOS) with OSA and a non-OSA control. Anatomically accurate three-dimensional geometries of the upper airway and soft-palate were reconstructed for both subjects using magnetic resonance (MR) images. A fluid-structure interaction (FSI) shape registration analysis was performed using subject-specific values of flow rate to iteratively compute the biomechanical properties of the soft-palate. The optimized shear modulus for the control was 38 percent higher than the corresponding value for the OSA patient. The proposed computational FSI model was then employed for planning surgical treatment for the apneic subject. A virtual surgery comprising of a combined adenoidectomy, palatoplasty and genioglossus advancement was performed to estimate the resulting post-operative patterns of airflow and tissue displacement. Maximum flow velocity and velopharyngeal resistance decreased by 80 percent and 66 percent respectively following surgery. Post-operative flow-induced forces on the anterior and posterior faces of the soft-palate were equilibrated and the resulting magnitude of tissue displacement was 63 percent lower compared to the pre-operative case. Results from this pilot study indicate that FSI computational modeling can be employed to characterize the mechanical properties of pharyngeal tissue and evaluate the effectiveness of various upper airway surgeries prior to their application.

The Role of Functional Respiratory Imaging in Treatment Selection of Children With Obstructive Sleep Apnea and Down Syndrome

STUDY OBJECTIVES

The complexity of the pathogenesis of obstructive sleep apnea (OSA) in children with Down syndrome (DS) is illustrated by a prevalence of residual OSA after adenotonsillectomy. The aim of this study was to investigate whether upper airway imaging combined with computation fluid dynamics could characterize treatment outcome after adenotonsillectomy in these children.

METHODS

Children with DS and OSA were prospectively included. All children underwent an evaluation of the upper airway and an ultra-low dose computed tomography scan of the upper airway before adenotonsillectomy. The upper airway tract was extracted from the scan and combined with computational fluid dynamics. Results were evaluated using control polysomnography after adenotonsillectomy.

RESULTS

Thirty-three children were included: 18 boys, age 4.3 ± 2.3 years, median body mass index z-score 0.6 (-2.9 to 3.0), and median obstructive apnea-hypopnea index was 15.7 (3-70) events/h. The minimal upper airway cross-sectional area was significantly smaller in children with more severe OSA (P = .03). Nineteen children underwent a second polysomnography after adenotonsillectomy. Seventy-nine percent had persistent OSA (obstructive apneahypopnea index > 2 events/h). A greater than 50% decrease in obstructive apnea-hypopnea index was observed in 79% and these children had a significantly higher volume of the regions below the tonsils.

CONCLUSIONS

This is the first study to characterize treatment outcome in children with DS and OSA using computed tomography upper airway imaging. At baseline, children with more severe OSA had a smaller upper airway. Children with a less favorable response to adenotonsillectomy had a smaller volume of regions below the tonsils, which could be due to enlargement of the lingual tonsils, glossoptosis, or macroglossia.

COMMENTARY

A commentary on this article appears in this issue on page 501.

Analyses of aerodynamic characteristics of the oropharynx applying CBCT: obstructive sleep apnea patients versus control subjects

OBJECTIVES

To determine the most relevant aerodynamic characteristic of the oropharynx related to the collapse of the upper airway in obstructive sleep apnea (OSA) patients; and to determine the correlation between the most relevant aerodynamic characteristic(s) of the oropharynx and anatomical characteristics of the oropharynx in OSA patients.

METHODS

31 mild to moderate OSA patients (mean ± SD age = 43.5 ± 9.7 years) and 13 control subjects (mean ± SD age = 48.5 ± 16.2 years) were included in this prospective study. The diagnosis of OSA patients was based on an overnight polysomnographic recording. To exclude the presence of OSA in the control subjects, they were asked to fill out a validated questionnaire to determine the risk of OSA. NewTom5G cone beam CT (CBCT) scans were obtained from both OSA patients and control subjects. Computational models of the oropharynx were reconstructed based on CBCT images. The aerodynamic characteristics of the oropharynx were calculated based on these computational models. Pearson correlation analysis was used to analyse the correlation between the most relevant aerodynamic characteristic(s) and anatomical characteristics of the oropharynx in OSA patients.

RESULTS

Compared with controls, the airway resistance during expiration (R_ex) of the OSA patients was significantly higher (p = 0.04). There was a significant negative correlation between R_ex and the minimum cross-sectional area (CSA_min) of the oropharynx (r = -0.41, p = 0.02), and between R_ex and the volume of the oropharynx (r = -0.48, p = 0.01) in OSA patients. After excluding an outlier, there is only significant correlation between R_ex and the CSA_min of the oropharynx (r = -0.45, p = 0.01).

CONCLUSIONS

Within the limitations of this study, we concluded that the most relevant aerodynamic characteristic of the oropharynx in the collapse of the upper airway in OSA patients is R_ex. Therefore, the repetitive collapse of the upper airway in OSA patients may be explained by a high R_ex, which is related to the CSA_min of the oropharynx.

Predicting the effect of treatment in paediatric OSA by clinical examination and functional respiratory imaging

OBJECTIVE

The aim of this study was to investigate whether functional respiratory imaging (FRI) or clinical examination could predict treatment outcome for obstructive sleep apnea (OSA) in normal-weight, non-syndromic children.

METHODS

Normal weight children diagnosed with OSA by polysomnography were prospectively included. All children got a thorough evaluation and an ultra-low dose computed tomography scan of the upper airway (UA). A 3-D reconstruction was built combined with computational fluid dynamics for FRI. Decisions on the need and type of surgery were based upon findings during drug-induced sleep endoscopy. A second polysomnography was performed 3-12 months after surgery.

RESULTS

Ninety-one children were included: 62 boys, 5.0 ± 2.7 years, and BMI z-score of -0.1 ± 1.2. Children with more severe OSA had a smaller volume of the overlap region between the adenoids and tonsils. Nineteen out of 60 patients had persistent OSA (oAHI >2/h). A lower conductance in the UA and a higher tonsil score predicted successful treatment.

CONCLUSIONS

A less constricted airway, as characterized by both FRI and a lower tonsil score, was associated with a less favorable response to (adeno) tonsillectomy. Further studies after treatment using FRI and DISE are warranted to further characterize the UA of these subjects.

Efficacy of minimally invasive tonsil surgery for treatment of obstructive sleep apnea-hypopnea syndrome in children

This study aimed to investigate the efficacy of minimally invasive tonsil surgery for the treatment of obstructive sleep apnea-hypopnea syndrome (OSAHS) in children. Tonsil ablation or turbinate reduction was performed on 49 pediatric patients with OSAHS by minimally invasive tonsil surgery. In order to evaluate the efficacy of surgery, a comparison was conducted between pre-operation and post-operation data in terms of the symptoms, signs and polysomnography test. Total effectiveness rate of the surgery was 83.7%. Subgroup analysis was also performed based on the severity of their conditions: mild, moderate, and severe groups had an effectiveness rate of 90.0, 88.9, and 66.7%, respectively (Hc=6.665, P<0.05). Postoperatively, the apnea-hypopnea index, the minimum oxygen saturation (SaO2), and corresponding symptoms improved compared to pre-operation conditions (P<0.05). Minimally invasive tonsil surgery was a safe and effective method for treating OSAHS in children.

Computational fluid dynamics study of the pharyngeal airway space before and after mandibular setback surgery in patients with mandibular prognathism

The purpose of this study was to investigate the relationship between the pressure drop in the pharyngeal airway space (ΔP_PAS) and the minimum cross-sectional area (minCSA) of the pharyngeal airway before and after mandibular setback surgery using computational fluid dynamics, in order to prevent iatrogenic obstructive sleep apnoea. Eleven patients with mandibular prognathism underwent bilateral sagittal split osteotomy for mandibular setback. Three-dimensional models of the upper airway were reconstructed from preoperative and postoperative computed tomography images, and simulations were performed using computational fluid dynamics. ΔP_PAS and the minCSA of the pharyngeal airway were calculated, and the relationship between them was evaluated by non-linear regression analysis. In all cases, the minCSA was found at the level of the velopharynx. After surgery, ΔP_PAS increased significantly and the minCSA decreased significantly. The non-linear regression equation expressing the relationship between these variables was ΔP_PAS=3.73×minCSA^-2.06. When the minCSA was <1cm², ΔP_PAS increased greatly. The results of this study suggest that surgeons should consider bimaxillary orthognathic surgery rather than mandibular setback surgery to prevent the development of iatrogenic obstructive sleep apnoea when correcting a skeletal class III malocclusion.

Anatomically correct visualization of the human upper airway using a high-speed long range optical coherence tomography system with an integrated positioning sensor

The upper airway is a complex tissue structure that is prone to collapse. Current methods for studying airway obstruction are inadequate in safety, cost, or availability, such as CT or MRI, or only provide localized qualitative information such as flexible endoscopy. Long range optical coherence tomography (OCT) has been used to visualize the human airway in vivo, however the limited imaging range has prevented full delineation of the various shapes and sizes of the lumen. We present a new long range OCT system that integrates high speed imaging with a real-time position tracker to allow for the acquisition of an accurate 3D anatomical structure in vivo. The new system can achieve an imaging range of 30 mm at a frame rate of 200 Hz. The system is capable of generating a rapid and complete visualization and quantification of the airway, which can then be used in computational simulations to determine obstruction sites.

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.

Computational fluid dynamics upper airway effective compliance, critical closing pressure, and obstructive sleep apnea severity in obese adolescent girls

Obstructive sleep apnea syndrome (OSAS) is associated with anatomical abnormalities restricting upper airway size and functional factors decreasing pharyngeal dilator activity in sleep. In this study we hypothesized that OSAS is also associated with altered pharyngeal mechanical compliance during wakefulness. Five OSAS and six control obese girls between 14 and 18 years of age were studied. All underwent polysomnography, critical closing pressure (P_crit) studies, and dynamic MRI of the upper airway during awake tidal breathing. Effective airway compliance was defined as the slope of cross-sectional area vs. average pressure between maximum inspiration and maximum expiration along the pharyngeal airway. Pharyngeal pressure fields were calculated by using image-based computational fluid dynamics and nasal resistance. Spearman correlations were calculated to test associations between apnea-hypopnea index (AHI), P_crit, and airway compliance. Effective compliances in the nasopharynx (C_NP) and velopharynx (C_VP) were lower and negative in OSAS compared with controls: -4.4 vs. 1.9 (mm²/cmH₂O, P = 0.012) and -2.1 vs. 3.9 (mm²/cmH₂O, P = 0.021), respectively, suggesting a strong phasic pharyngeal dilator activity during inspiration in OSAS compared with controls. For all subjects, C_NP and AHI correlated negatively (r_S = -0.69, P = 0.02), and passive P_crit correlated with C_NP (r_S = -0.76, P = 0.006) and with AHI (r_S = 0.86, P = 0.0006). Pharyngeal mechanics obtained during wakefulness could be used to characterize subjects with OSAS. Moreover, negative effective compliance during wakefulness and its correlation to AHI and P_crit suggest that phasic dilator activity of the upper pharynx compensates for negative pressure loads in these subjects.

Is the difference in the volume of the pharyngeal space, as measured by acoustic pharyngometry, before and after tonsillectomy proportional to the volume of the excised tonsils?

OBJECTIVE

Adenotonsillectomy is recognized as an effective therapy for snoring and sleep disorders in children. It is important to understand whether adenotonsillectomy significantly increases the volume of the pharyngeal space. The goal of this study was to evaluate the change in oropharyngeal volume after adenotonsillectomy and the correlation of this change with the objective volume of the tonsils and body mass index.

METHODS

We included 27 subjects (14 males) with snoring caused by tonsil and adenoid hypertrophy. The mean age of the subjects was 7.92 (±2.52) years. Children with craniofacial malformations or neuromuscular diseases or syndromes were excluded. The parents/caregivers answered an adapted questionnaire regarding sleep-disordered breathing. All patients were subjected to weight and height measurements and body mass index was calculated. The subjects underwent pharyngometry before and after adenotonsillectomy and the volume of both excised tonsils together was measured in cm3 in the operating room.

RESULTS

Pharyngometric analysis showed that the mean pharyngeal volume was 28.63 (±5.57) cm3 before surgery and 31.23 (±6.76) cm3 after surgery; the volume of the oropharynx was significantly increased post-surgery (p=0.015, Wilcoxon test). No correlation was found between the objective tonsil volume and the post-surgical volume increase (p=0.6885). There was a fair correlation between the oropharyngeal volume and body mass index (p=0.0224).

CONCLUSION

Adenotonsillectomy increases the volume of the pharyngeal space, but this increase does not correlate with the objective tonsil size. Furthermore, greater BMI was associated with a smaller increase in the pharyngeal volume. Oropharyngeal structures and craniofacial morphology may also play a role in the increase in oropharyngeal volume.

Bone Age and Serum Osteocalcin Levels in Children With Obstructive Sleep Apnea Hypopnea Syndrome Before and After Adenotonsillectomy

Our aim was to study the changes in bone age and serum osteocalcin levels before and after adenotonsillectomy (AT) in children with obstructive sleep apnea hypopnea syndrome (OSAHS). A total of 58 OSAHS children (37 males and 21 females) with the mean age of 6.68 ± 1.11 years were enrolled and assessed by x-ray-based bone age estimation and enzyme-linked immunosorbent assay-based measurement of serum osteocalcin levels, before surgery and 6 months after AT. SPSS 19.0 software was used for statistical analysis. Our results revealed that bone age and serum osteocalcin levels in OSAHS patients were significantly lower than normal controls before AT (P < 0.05). Within 6 months after surgery, the bone age and the serum osteocalcin levels were significantly elevated in OSAHS patients (P < 0.05), compared with those before surgery. Serum osteocalcin levels and bone age are negatively correlated with apnea-hypopnea index, oxygen desaturation index, the percentage of the total recorded time spent below 90% oxygen saturation, and Epworth sleepiness scale scores (all P < 0.05). Our findings suggested that bone age and serum osteocalcin levels may be correlated with the development of OSAHS in children. AT may improve bone age and serum osteocalcin levels in OSAHS children.

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.