Large Eddy Simulation and Reynolds-Averaged Navier-Stokes modeling of flow in a realistic pharyngeal airway model: an investigation of obstructive sleep apnea.

Computational fluid dynamics techniques employing primarily steady Reynolds-Averaged Navier-Stokes (RANS) methodology have been recently used to characterize the transitional/turbulent flow field in human airways. The use of RANS implies that flow phenomena are averaged over time, the flow dynamics not being captured. Further, RANS uses two-equation turbulence models that are not adequate for predicting anisotropic flows, flows with high streamline curvature, or flows where separation occurs. A more accurate approach for such flow situations that occur in the human airway is Large Eddy Simulation (LES). The paper considers flow modeling in a pharyngeal airway model reconstructed from cross-sectional magnetic resonance scans of a patient with obstructive sleep apnea. The airway model is characterized by a maximum narrowing at the site of retropalatal pharynx. Two flow-modeling strategies are employed: steady RANS and the LES approach. In the RANS modeling framework both k-epsilon and k-omega turbulence models are used. The paper discusses the differences between the airflow characteristics obtained from the RANS and LES calculations. The largest discrepancies were found in the axial velocity distributions downstream of the minimum cross-sectional area. This region is characterized by flow separation and large radial velocity gradients across the developed shear layers. The largest difference in static pressure distributions on the airway walls was found between the LES and the k-epsilon data at the site of maximum narrowing in the retropalatal pharynx.

[1]  R. Rogers,et al.  Obstructive sleep apnea. , 1996, The New England journal of medicine.

[2]  Timothy J. Pedley,et al.  Pulmonary Fluid Dynamics , 1977 .

[3]  K. Lilly On the application of the eddy viscosity concept in the Inertial sub-range of turbulence , 1966 .

[4]  J. P. V. Doormaal,et al.  ENHANCEMENTS OF THE SIMPLE METHOD FOR PREDICTING INCOMPRESSIBLE FLUID FLOWS , 1984 .

[5]  S A Shea,et al.  Pharyngeal size and shape during wakefulness and sleep in patients with obstructive sleep apnoea. , 1989, The Quarterly journal of medicine.

[6]  Hukam Mongia,et al.  Recent Progress in Comprehensive Modeling of Gas Turbine Combustion , 2008 .

[7]  Michael Breuer,et al.  Improving prediction of aerosol deposition in an idealized mouth using large-Eddy simulation. , 2006, Journal of aerosol medicine : the official journal of the International Society for Aerosols in Medicine.

[8]  T. E. Tokyay,et al.  Validation of a Large-Eddy Simulation Model to Simulate Flow in Pump Intakes of Realistic Geometry , 2006 .

[9]  P. Nithiarasu,et al.  Laminar and turbulent flow calculations through a model human upper airway using unstructured meshes , 2006 .

[10]  Steven H Frankel,et al.  Modeling transition to turbulence in eccentric stenotic flows. , 2008, Journal of biomechanical engineering.

[11]  Thierry Poinsot,et al.  Comparison of LES, RANS and experiments in an aeronautical gas turbine combustion chamber , 2007 .

[12]  A. Pack,et al.  Upper airway size analysis by magnetic resonance imaging of children with obstructive sleep apnea syndrome. , 2003, American journal of respiratory and critical care medicine.

[13]  P. Moin,et al.  A numerical method for large-eddy simulation in complex geometries , 2004 .

[14]  D. Hudgel,et al.  Variable site of airway narrowing among obstructive sleep apnea patients. , 1986, Journal of applied physiology.

[15]  Soo-Jin Jeong,et al.  Numerical investigation on the flow characteristics and aerodynamic force of the upper airway of patient with obstructive sleep apnea using computational fluid dynamics. , 2007, Medical engineering & physics.

[16]  Clement Kleinstreuer,et al.  Laminar-to-turbulent fluid-particle flows in a human airway model , 2003 .

[17]  Jochen Fröhlich,et al.  A Comparative Study of the Turbulent Flow Over a Periodic Arrangement of Smoothly Contoured Hills , 2006 .

[18]  J. Fröhlich,et al.  Investigation of wall-function approximations and subgrid-scale models in large eddy simulation of separated flow in a channel with streamwise periodic constrictions , 2003 .

[19]  Clement Kleinstreuer,et al.  Micro-particle transport and deposition in a human oral airway model , 2002 .

[20]  M. Badr,et al.  Progressive retropalatal narrowing preceding obstructive apnea. , 1998, American journal of respiratory and critical care medicine.

[21]  Sang-Jin Sung,et al.  Customized three-dimensional computational fluid dynamics simulation of the upper airway of obstructive sleep apnea. , 2006, The Angle orthodontist.

[22]  P. G. Young,et al.  A computational fluid dynamics study of inspiratory flow in orotracheal geometries , 2007, Medical & Biological Engineering & Computing.

[23]  D. Wilcox Turbulence modeling for CFD , 1993 .

[24]  Jianren Fan,et al.  Large eddy simulation of inhaled particle deposition within the human upper respiratory tract , 2007 .

[25]  Laszlo Fuchs,et al.  Comparison of Some SGS Models for Turbulent Fluxes of Momentum and a Passive Scalar in Swirling Co-Annular Jets , 2004 .

[26]  A I Pack,et al.  Upper airway and soft tissue anatomy in normal subjects and patients with sleep-disordered breathing. Significance of the lateral pharyngeal walls. , 1995, American journal of respiratory and critical care medicine.

[27]  X Y Luo,et al.  LES modelling of flow in a simple airway model. , 2004, Medical engineering & physics.

[28]  J. Udupa,et al.  Computational fluid dynamics modeling of the upper airway of children with obstructive sleep apnea syndrome in steady flow. , 2006, Journal of biomechanics.

[29]  P. Parizel,et al.  Correlation between severity of sleep apnea and upper airway morphology based on advanced anatomical and functional imaging. , 2007, Journal of biomechanics.

[30]  Atul Malhotra,et al.  Obstructive sleep apnoea , 2002, The Lancet.