Dennis O Frank-Ito1, Matthew Wofford2, Jeffry D Schroeter3, Julia S Kimbell2. 1. 1 Division of Otolaryngology, Head and Neck Surgery, Duke University Medical Center , Durham, North Carolina. 2. 2 Department of Otolaryngology, Head and Neck Surgery, University of North Carolina , Chapel Hill, North Carolina. 3. 3 Applied Research Associates, Inc. , Raleigh, North Carolina.
Abstract
BACKGROUND: There are methodological ambiguities in the literature on mesh refinement analysis for computational fluid dynamics (CFD) modeling of physiologically realistic airflow dynamics and particle transport in the human sinonasal cavity. To investigate grid independence in discretization of the (sino)nasal geometry, researchers have considered CFD variables such as pressure drop, velocity profile, wall shear, airflow, and particle deposition fractions. Standardization in nasal geometry is also lacking: unilateral or bilateral nasal cavities with and without paranasal sinuses have been used. These methodological variants have led to inconsistencies in establishing grid-independent mesh densities. The aim of this study is to provide important insight in the role of mesh refinement analysis on airflow and particle deposition in sinonasal airway modeling. METHODS: A three-dimensional reconstruction of the complete sinonasal cavity was created from computed tomography images of a subject who had functional endoscopic sinus surgery. To investigate airflow grid independence, nine different tetrahedral mesh densities were generated. For particle transport mesh refinement analysis, hybrid tetrahedral-prism elements with near-wall prisms ranging from 1 to 6 layers were implemented. Steady-state, laminar inspiratory airflow simulations under physiologic pressure-driven conditions and nebulized particle transport simulations were performed with particle sizes ranging from 1-20 μm. RESULTS: Mesh independence for sinonasal airflow was achieved with approximately 4 million unstructured tetrahedral elements. The hybrid mesh containing 4 million tetrahedral cells with three prism layers demonstrated asymptotic behavior for sinonasal particle deposition. Inclusion of boundary prism layers reduced deposition fractions relative to tetrahedral-only meshes. CONCLUSIONS: To ensure numerically accurate simulation results, mesh refinement analyses should be performed for both airflow and particle transport simulations. Tetrahedral-only meshes overpredict particle deposition and are less accurate than hybrid tetrahedral-prism meshes.
BACKGROUND: There are methodological ambiguities in the literature on mesh refinement analysis for computational fluid dynamics (CFD) modeling of physiologically realistic airflow dynamics and particle transport in the human sinonasal cavity. To investigate grid independence in discretization of the (sino)nasal geometry, researchers have considered CFD variables such as pressure drop, velocity profile, wall shear, airflow, and particle deposition fractions. Standardization in nasal geometry is also lacking: unilateral or bilateral nasal cavities with and without paranasal sinuses have been used. These methodological variants have led to inconsistencies in establishing grid-independent mesh densities. The aim of this study is to provide important insight in the role of mesh refinement analysis on airflow and particle deposition in sinonasal airway modeling. METHODS: A three-dimensional reconstruction of the complete sinonasal cavity was created from computed tomography images of a subject who had functional endoscopic sinus surgery. To investigate airflow grid independence, nine different tetrahedral mesh densities were generated. For particle transport mesh refinement analysis, hybrid tetrahedral-prism elements with near-wall prisms ranging from 1 to 6 layers were implemented. Steady-state, laminar inspiratory airflow simulations under physiologic pressure-driven conditions and nebulized particle transport simulations were performed with particle sizes ranging from 1-20 μm. RESULTS: Mesh independence for sinonasal airflow was achieved with approximately 4 million unstructured tetrahedral elements. The hybrid mesh containing 4 million tetrahedral cells with three prism layers demonstrated asymptotic behavior for sinonasal particle deposition. Inclusion of boundary prism layers reduced deposition fractions relative to tetrahedral-only meshes. CONCLUSIONS: To ensure numerically accurate simulation results, mesh refinement analyses should be performed for both airflow and particle transport simulations. Tetrahedral-only meshes overpredict particle deposition and are less accurate than hybrid tetrahedral-prism meshes.
Authors: Jarrod A Keeler; Aniruddha Patki; Charles R Woodard; Dennis O Frank-Ito Journal: J Aerosol Med Pulm Drug Deliv Date: 2015-08-13 Impact factor: 2.849
Authors: Masoud Gh Moghaddam; Guilherme J M Garcia; Dennis O Frank-Ito; Julia S Kimbell; John S Rhee Journal: Int J Comput Assist Radiol Surg Date: 2020-02-20 Impact factor: 2.924
Authors: Elizabeth L Perkins; Saikat Basu; Guilherme J M Garcia; Robert A Buckmire; Rupali N Shah; Julia S Kimbell Journal: Otolaryngol Head Neck Surg Date: 2017-11-21 Impact factor: 3.497
Authors: Karl Bass; Susan Boc; Michael Hindle; Kelley Dodson; Worth Longest Journal: J Aerosol Med Pulm Drug Deliv Date: 2018-12-15 Impact factor: 2.849
Authors: Bryan M Brandon; Wesley H Stepp; Saikat Basu; Julia S Kimbell; Brent A Senior; William W Shockley; J Madison Clark Journal: Laryngoscope Date: 2020-05-04 Impact factor: 3.325
Authors: Julia S Kimbell; Saikat Basu; Guilherme J M Garcia; Dennis O Frank-Ito; Frances Lazarow; Erica Su; Dimitry Protsenko; Zhongping Chen; John S Rhee; Brian J Wong Journal: Lasers Surg Med Date: 2018-07-26 Impact factor: 4.025
Authors: Dennis O Frank-Ito; Julia S Kimbell; Azadeh A T Borojeni; Guilherme J M Garcia; John S Rhee Journal: Clin Biomech (Bristol, Avon) Date: 2018-12-19 Impact factor: 2.063