The objective of this study is to systematically assess the influences of the larynopharyneal anatomical details on air flow and particle behaviors during exhalation by means of image-based modeling. A physiologically realistic nose-throat airway was developed with medical images. Individual airway anatomy such as uvula, pharynx, and larynx were then iso lated for examination by progressively simplifying this image-based model geometry. Low Reynolds number (LRN) k - w model and Langrangian tracking model were used to simulate the dynamics of airflow and particle transport for a wide range of exhalation conditions (4 - 45 L/min) and particle sizes (1 nm - 1 μm). Results showed that pharyngeal anatomical details exerted a significant impact on breathing resistance and particle profiles. Abrupt pressure drop resulting from the uvula-related airway obstruction was observed. Even though the total deposition rate in the nasal airway is largely unaf fected by the upstream effect, the local deposition patterns vary notably. Results of this study also indicate that the pressure drop appears to be an appropriate parameter to characterize the geometric variations for diffusive depositions. Inclusion of pressure drop (D0.5Q-0.62dp0.07) gives an improved correlation than using the conventional diffusion factor (D0.5Q﹣0.28).
A significant issue in evaluating nasal airflow and aerosol deposition during exhalation is determining the extent to which the larynx and pharynx affect their behaviors before entering the nasal cavity. The larynx consists of vocal folds that form an elliptical or triangular crosssectional area of flow constriction (i.e., glottal aperture). It is approximately 6 cm long and its cross-sectional area varies with mean flow rate [
The nasal deposition of submicrometer aerosols during exhalation has been considered in human volunteers in a limited number of studies [3-5]. A common disadvantage of such in vivo studies is the difficulty in determining local deposition values even though imaging methods now make it possible [
A number of numerical studies have considered the transport and deposition of fine and ultrafine particles in the nasal airways [10-15]. Similar CFD studies have also evaluated the transport and absorption of vapors in the nasal passages [16-18]. Comparisons of CFD results to experimental deposition data in the nasal cavity are often difficult due to differences in the geometric models that vary in their physical realisms. Untested modeling assumptions limit the accuracy of both in vitro and numerical models. During exhalation, the influence of the upstream laryngeal and pharyngeal anatomy on transport and deposition in the nasal cavity may be significant. Longest and Vinchurkar [
The validity of assuming simplified airway geometry to facilitate modeling is routinely adopted in previous studies but remains largely untested. The objective of this study is to systematically evaluate the effect of the laryngopharyngeal anatomical details on both airflow and aerosol depositions in the nasal airways of an adult. Starting from an image-based anatomically accurate airway model, we will progressively simplify this geometry, one anatomy at a time, to gain a more accurate understanding of the physiological roles of each anatomy in regulating airflows, breathing resistance, and particle filtering. The glottal aperture, pharynx, and uvula are especially of interest in this study.
In order to assess the impact of the laryngopharyngeal anatomy on airflow and particle deposition in the nasal airways, four computational models with varying geometric details have been considered. We started from a highly realistic nasal-laryngeal airway model that was recently developed based on MRI images of a 53-yearold male (weight 73 kg and height 173 cm) [
Bedsides the anatomically accurate nasal passages, three other key features characterize this image-based model, namely, a hanging uvula, a flat-plated pharynx, and a triangular-shaped glottal aperture. Movement of the uvula alters the airway morphology that connects the nasoand oropharynx. In this study, the uvula rests on the back of the throat and partially obstructs the inferior nasopharynx, resulting in two flow passages. This might be attributed to the supine position of the patient during image acquisition. Posterior to the uvula, the pharynx is featured by narrow and flat air channels, which converge into the wedge-shaped glottis in the larynx. In order to approximate the in vivo airway morphology that was captured in the images, anatomical details such as the epiglottis and the two pharyngeal sinuses on either side of the larynx were also retained, as shown in
As a first simplification, the hanging uvula was removed in Model 1, resulting in an un-obstructed nanopharynx in Model 2 (
ratio of 2.5 for healthy subjects. Thus, for Model 3, we assume that the pharynx is approximately an elliptical pipe with the major axis 2.5 times as large as the minor axis. The fourth geometry model (i.e., Model 4) excludes the triangular-shaped larynx and is composed of an elliptical pipe downstream of the nasal cavity (
During exhalation, airflow enters the upper trachea, travels through the larynx, pharynx, and nasal turbinates, and exits the nostrils. By retaining the laryngeal-pharyngeal region, more physiologically realistic inlet flow conditions are provided to the nasal cavity. As a result, the flow field and particle deposition characteristics in the nasal cavity considered in this study will better represent in vivo conditions compared with some previous studies that had excluded the larynx and pharynx. Therefore, inclusion of the laryngeal-pharyngeal region allows for the results of this computational study to be directly compared with previous in vivo nasal deposition data. In order to characterize deposition distributions, the nosethroat airway were divided into different sub-regions that include the vestibule, valve region (VR), turbinate region (TR), olfactory region (OR), nasopharynx (NP), pharynx, and larynx (
Steady exhalation was assumed for all simulations with uniform velocity profiles and particle distributions at the tracheal inlet (
The flow conditions considered in this study are assumed to be isothermal and incompressible. The mean inlet Reynolds number at the trachea varies from 368 to 3302. The maximum Reynolds number based on the hydraulic diameter of the glottal aperture is approximately 8037. Therefore, laminar, transitional, and fully turbulent conditions in the nasal-laryngeal model are expected. To resolve these multiple flow regimes, the low Reynolds number (LRN) k - w model was selected based on its ability to accurately predict pressure drop, velocity profiles, and shear stress for transitional and turbulent flows. Moreover, the LRN k - ω model was shown to provide an accurate solution for laminar flow as the turbulent viscosity approaches zero [
The transport and deposition of the submicrometer particles are simulated with a well-tested discrete Lagrangian tracking model, which is enhanced with userdefined functions (UDFs) accounting for the finite particle inertial effects that might be significant for submitcrometer particle depositions. The aerosols evaluated in this study had a tracheal Stokes number (Stk =) with a range of 1.68 ´ 10−8 to 1.0 ´ 10−3 and were assumed to be dilute and to not influence the continuous-phase, i.e., one-way coupled particle motion. In our previous studies, the UDF-enhanced Lagrangian model with near-wall interpolation algorithm [33,34] has been shown to provide close agreement with experimental deposition data in upper respiratory airways for both submicrometer [
To solve the governing mass and momentum conservation equations in each of the cases considered, the CFD package ANSYS Fluent was employed. User-supplied Fortran and C programs were implemented for the calculation of initial particle profiles, particle deposition factors, grid convergence, and deposition enhancement factors. For this study, a specific set of user-defined functions were applied that considered the Brownian force, anisotropic turbulence effect, and near-wall velocity interpolation. All transport equations were discretized to be at least second order accurate in space. A segregated implicit solver was applied to evaluate the resulting linear system of equations. This solver uses the Gauss-Seidel method in conjunction with an algebraic multigrid approach to improve the calculation performance on tetrahedral meshes. Convergence of the flow field solution was assumed when the global mass residual was reduced from its original value by five orders of magnitude and the residual-reduction-rates for both mass and momentum were sufficiently small.
The computational meshes of the four nasal-laryngeal airway models were generated with ANSYS IECM CFD (Ansys, Inc). Due to the high complexity of the model geometries, unstructured tetrahedral meshes were generated with high-resolution prismatic cells in the near-wall region (
For discrete Lagrangian tracking, the number of seeded particles required to produce count-independent depositions was considered. Particle count sensitivity testing was performed by incrementally releasing groups of 10,000 particles. The number of groups was increased until the deposition rate change was less than 1%. Due to the low deposition rates, more particles were required for fine aerosols to generate count-independent results compared with ultrafine aerosols. The final number of particles tracked for 1 - 40 nm and 100 - 1000 nm aerosols were 150,000 and 600,000, respectively.
Dimensions of the four models considered are shown in
There is less variation in the geometric parameters of Models 3 and 4. In the laryngeal-pharyngeal region, both models (i.e., 3 and 4) have smaller flow area (dh) than Models 1 and 2. The dip in cross-sectional area at a distance of 47 mm corresponds to the epiglottis that spreads above the glottal aperture at an angle of approximately 30˚ to the posterior wall of the pharynx (Fig ure 1(a), Slice C-C’). The larynx is located at a distance of 78 mm or so downstream of the nasopharynx.
The influence of the pharyngo-laryngeal geometries on exhalation breathing resistance is illustrated in
as a function of the axial distance from the tracheal inlet at an exhalation flow rate of 30 L/min for the four models considered. For Model 1, the maximum pressure drop is observed in the pharynx/nasopharynx region which is about 170 Pa in magnitude and is about 76% of the expiratory pressure drop between the upper trachea and nostrils. This abrupt pressure drop results from the airway obstruction as the hanging uvula partially blocks the pharyngeal passage. As discussed, Model 1 was based on images of a subject in a supine position during image acquisition. It is expected that changing from upright to supine position has the potential of causing pharyngeal obstruction as the uvula yields due to its own weight and that the degree of obstruction varies depending on the elasticity properties of the subject’s soft palate and uvula. Considering that this pressure drop is about three times that of the laryngeal region, which is the maximum of the other three models (Models 2, 3, 4), it is possible that the patient of Model 1 has an over-flexible uvula and hence a high degree of airway obstruction. Respiratory tissues in the pharyngeal region are quite collapsible. Severe pressure drop could possibly trigger or expedite the collapse of airway walls and induce the symptoms of hyponea (i.e., reduced airflow capacity) or apnea (i.e., complete airway obstruction). Meanwhile, flow instability may cause vibration of the pliable pendant uvula, resulting in snoring of the supine subject. It is therefore interesting to postulate that an over-flexible uvula may be one major reason for sleep disorders such as sleep apnea, whose constriction of the pharyngeal airway and the resultant elevated pressure drop induces further pharyngeal airway collapse.
For Models 2 and 3, which corresponds to the upright position, the glottis and nasal cavity each constitutes a breathing resistance of about 50 Pa. However, there is less pressure recovery in the pharynx of Model 2 comparing to that of Model 3. Considering the fourth model, which totally excludes the anatomical details of glottis and pharynx, there is negligible pressure drop within the larynx-pharynx region, and almost all the pressure drop comes from the convoluted nasal passages.