Abstract
BACKGROUND: The human respiratory airway undergoes dramatic growth during infancy and childhood, which induces substantial variability in air flow pattern and particle deposition. However, deposition studies have typically focused on adult subjects, the results of which cannot be readily extrapolated to children. We developed models to quantify the growth of human nasal-laryngeal airways at early ages, and to evaluate the impact of that growth on breathing resistance and aerosol deposition.
METHODS: Four image-based nasal-laryngeal models were developed from 4 children, ages 10 days, 7 months, 3 years, and 5 years, and were compared to a nasal-laryngeal model of a 53-year-old adult. The airway dimensions were quantified in terms of different parameters (volume, cross-section area, and hydraulic diameter) and of different anatomies (nose, pharynx, and larynx). Breathing resistance and aerosol deposition were computed using a high-fidelity fluid-particle transport model, and were validated against the measurements made with the 3-dimensional models fabricated from the same airway computed tomography images.
RESULTS: Significant differences in nasal morphology were observed among the 5 subjects, in both morphology and dimension. The turbinate region appeared to experience the most noticeable growth during the first 5 years of life. The nasal airway volume ratios of the 10-day, 7-month, 3-year, and 5-year-old subjects were 6.4%, 18.8%, 24.2%, and 40.3% that of the adult, respectively. Remarkable inter-group variability was observed in air flow, pressure drop, deposition fraction, and particle accumulation. The computational fluid dynamics predicted pressure drops and deposition fractions were in close agreement with in vitro measurements.
CONCLUSIONS: Age effects are significant in both breathing resistance and micrometer particle deposition. The image/computational-fluid-dynamics coupled method provides an efficient and effective approach in understanding patient-specific air flows and particle deposition, which have important implications in pediatric inhalation drug delivery and respiratory disorder diagnosis.
- nasal morphology
- child-adult discrepancy
- infants
- breathing resistance
- aerosol deposition
- pediatric drug delivery
Introduction
Considerable structural changes occur in human respiratory tracts during early ages, which lead to vast variations in breathing disorders that are specific to children. Compared to adults, infants and children are more susceptible to respiratory distresses, due to their immature defense mechanism. Respiratory disease remains a leading cause of childhood morbidity in the United States and other developed countries, and is a leading cause of childhood deaths worldwide. Only by understanding the developmental respiratory system can we formulate effective protocols to treat an impaired physiology in pediatric patients.
A number of in vivo studies have been conducted to understand the development of respiratory morphology of early ages. Summaries of such studies can be found in recent reviews by Bassham et al1 and Ghuman et al.2 However, these studies mainly focused on the lower respiratory tract. Very few studies have considered the development of upper airway morphology.3 Clinical techniques for evaluating nasal pathology include visual analog scales,4 rhinostereometry,5 video-endoscopy,6 rhinomanometry,7 and acoustic rhinometry.8 A combination of the above techniques is usually employed in practice, based on the specific performance of each technique.9–12 Compared with clinical studies, computational fluid dynamics predictions have the advantage of providing detailed information on air flow and aerosol deposition, such as regional dosage, which are more relevant to health outcome than the average deposition. However, to our knowledge, very few computational fluid dynamics studies have been reported on the nasal morphology of children; one exception is Xi et al,13 who examined the upper airway structure and associated air flow dynamics in a 5-year-old child. The general neglect of child and infant airways in previous studies may largely be attributed to limited accessibility of pediatric medical images to computational fluid dynamics researchers, as well as the complexities involved in constructing physiologically realistic models of nasal passages.
The objective of this study was to characterize the growth of the upper airway with models created from images from 4 children, ages 10 days, 7 months, 3 years, and 5 years, using a coupled image/computational fluid dynamics approach. Our specific aims were: to develop anatomically accurate airway models of children, based on computed tomography and magnetic resonance imaging images; to quantify the airway dimensions; to numerically and experimentally determine the nasal-laryngeal resistance; and to numerically determine the inhaled particle deposition. The results of this study may improve our understanding of developmental respiratory morphology and its effects on children's health responses to environmental exposures and inhaled therapies.
QUICK LOOK
Current knowledge
Aerosol deposition in the respiratory tract is affected by various well studied anatomic, equipment, and breathing-pattern variables, but most of the data have been from adults.
What this paper contributes to our knowledge
Anatomic differences associated with age affect total and regional aerosol deposition. There were important differences in airway anatomy, air-flow dynamics, and aerosol deposition between subjects ages 10 days, 7 months, 3 years, 5 years, and 53 years.
Methods
Nasal-Laryngeal Airway Models
One example was used to illustrate the computer method to develop respiratory airway models based on computed tomography and magnetic resonance images. Figure 1 shows the procedures of translating a 2-dimensional magnetic resonance image (of a healthy 5-year-old boy, weight 21 kg, height 109 cm)13 into a 3-dimensional model. The image tracings (see Fig. 1A) contained 128 slices, each 1.5 mm apart, that spanned from the nostrils to the upper trachea. The multi-slice tracings were segmented in medical image processing (Mimics, Materialise, Leuven, Belgium), based on the contrast between soft tissue and intranasal air, to isolate the nasal-laryngeal airway, which was further converted into a set of cross-sectional contours that define the airway. Based on these contours, an internal nasal surface geometry was constructed in engineering-design software (ANSYS Workbench, CAE Associates, Middlebury, Connecticut). The airway model can be either printed with 3-dimensional prototyping techniques for experimental purposes (see Fig. 1B) or discretized into a computational mesh for numerical analysis (see Fig. 1C). Due to the high complexity of the model geometry, an unstructured tetrahedral mesh was created with high-resolution pentahedral elements in the near-wall region (see Fig. 1C). The respiratory geometry retained in this example extended from the nostrils to the upper trachea. In particular, anatomical details such as the epiglottal fold and laryngeal sinuses were retained (see Fig. 1A). The resulting model was intended to accurately represent the morphology of the upper airway, with only minor surface smoothing.
Study Design
Age effects on airway morphology, air flow, pressure drop, and inhaled particle deposition were assessed with nasal-laryngeal models created from computed tomography scans of 4 children (a 10-day-old female, a 7-month-old female, a 3-year-old female, and a 5-year-old male), seen at Arkansas Children Hospital, and one adult (a 53-year-old male). No health information was disclosed, and the use of the images was approved by the institutional review board of the University of Arkansas for Medical Science.
Steady inhalations were assumed for all simulations, with a wide spectrum of breathing conditions (1–45 L/min). Comparison of the 5 models under equivalent physical activities (quiet breathing) were conducted, based on published respiratory parameters.14–16 Table 1 shows the literature-based, quiet-breathing respiratory parameters. Interestingly, the tidal volumes in all 5 models intimately correlated with the nasal-laryngeal airway volumes. As shown in Tables 1 and 2, relative to the adult model, the tidal volume ratios for the 10-day-old, 7-month-old, 3-year-old, and 5-year-old were, respectively, 4.4%, 17.4%, 26.0%, and 35.4%, while the airway volume ratios were respectively 6.4%, 18.8%, 24.2%, and 40.3%. Based on the nostril areas in Table 3, the inlet velocities were similar between the 5 models under quiet breathing conditions. Initial particle velocities were assumed to be the same as the local fluid velocity. The airway surface was assumed to be smooth and rigid, with a no-slip (flow velocity at the wall [uwall] = 0) condition.
Pressure Measurement
We measured the pressure drops through the nasal-laryngeal models with a manometer (Magnehelic Gage, Dwyer Instrument, Michigan City, Indiana). The model nostrils were open to room air, and a vacuum was connected to the manometer at the model outlet. We tested constant inhalation flows between 1 and 48 L/min, by adjusting the valve on the vacuum line. Flow was measured at the model outlet, with a flow meter (4140, TSI, Shoreview, Minnesota) positioned upstream of the flow valve.
Numerical Method
The low Reynolds number k-ω model was selected to simulate the air flow dynamics, based on its ability to accurately predict pressure drop, velocity profiles, and shear stress for multi-regime flows.13,17–19 The low Reynolds number k-ω model applies to both laminar and turbulent flow regimes. Deposition fraction was calculated as the ratio of the amount of particles deposited on the airway surface to the amount of particles entering the nostrils. The transport and deposition of inhaled particles are simulated with a well tested discrete Lagrangian tracking model, enhanced with near-wall treatment. The inhaled particles were assumed to be dilute and to have no influence on the continuous phase (ie, one-way coupled particle motion). In our previous studies the Lagrangian tracking model, enhanced with user-defined routines, was shown to provide a close match to experimental deposition data in upper respiratory airways for both sub-micrometer18 and micrometer particles.17
To establish grid-independent results, convergence sensitivity analysis was conducted, following the methods of Xi et al.17 The final grids for reporting flow field consisted of 1.6–2.0 million cells, with a thin 5-layer pentahedral grid in the near-wall region and a cell height of 0.05 mm in the first layer (see Fig. 1C). The final number of particles tracked was 60,000; increasing the number of tracked particles did not alter the deposition fractions.
Results
Airway Morphological and Dimensional Variations
There were substantial morphological and dimensional differences between the 5 nasal-laryngeal airway models (Figs. 2 and 3). Younger subjects have smaller nostrils, a shorter turbinate region, a narrower nasopharynx, and a narrower pharynx-larynx. The nostril shape is more circular at birth, becomes more oval during infancy and childhood, and becomes wedge-shaped in adulthood.20 Of particular interest is the turbinate region, which seemed to be undeveloped in both the 10-day-old model and the 7-month-old model, and was much simpler in morphology than in our older children and adult models. The inferior meatus is missing in the 10-day-old. The inferior and middle meatuses apparently become both larger and more complex with age (see Fig. 2).
Figure 3 shows a quantitative comparison of the 5 models' airway dimensions, in terms of coronal cross-sectional area and perimeter, as a function of distance from the nose tip. Tables 2 and 3 compare the airway volumes and surface areas in the different anatomic sections. The nasal airway volumes of the 10-day-old, 7-month-old, 3-year-old, and 5-year-old were 6.4%, 18.8%, 24.2%, and 40.3% that of the adult, respectively. The nasal airway surface areas were 22.7%, 30.0%, 41.7%, and 65.7% that of the adult, respectively (see Table 3). There were 2 major disparities in the nasopharynx and nasal valve. First, compared to the adult, the children had a much smaller and narrower nasopharynx lumen. Figure 3A shows that the nasopharynx cross-sectional areas of the 10-day-old (31 mm2), the 7-month-old (78 mm2), the 3-year-old (29 mm2), and the 5-year-old (170 mm2) were only 5.8%, 15%, 5.4%, and 32% that of the adult (535 mm2), respectively. The nasopharynx hydraulic diameters (hydraulic diameter = (4 × cross-sectional area)/perimeter) of the 10-day-old (5.7 mm), 7-month-old (8.1 mm), 3-year-old (4.0 mm), and 5-year-old (12.2 mm) were about 1/4, 1/3, 1/6, and 1/2 that of the adult (24.0 mm), respectively. Specifically, the 3-year-old seemed to have an abnormal nasopharynx that is narrower than the other 3 young subjects. The second major disparity is in the nostril-to-valve distances, which were much shorter in the 2 infants (8.0 mm in the 10-day-old, and 11.2 mm in the 7-month-old, see Fig. 3A) than in the 2 older children (both ∼20 mm) and the adult (27.2 mm). The nasal valve cross-sectional areas were likewise smaller in the 2 youngest subjects (ie, ∼42 mm2 in the 10-day-old, and ∼56 mm2 in the 7-month-old) than in the 2 older children and the adult (see Fig. 3A).
Air Flow and Breathing Resistance
Figure 4A shows the nasal air flows during quiet breathing, as stream lines within the right nasal passage. High-velocity flows are observed in the middle portion of the nasal passage in all the models. The main flow changes direction dramatically from the nostrils to the nasopharynx, forming a nearly 180° curvature. However, this curvature is less severe in the 7-month-old model than in the other 3 models, presumably resulting from a more back-tilted head position during image acquisition. The air speed was higher in the nasopharynx of the 3-year-old, due to the severe nasopharynx constriction. No recirculation zone was observed in the nasopharynx of the 3 young subjects, due to a much smaller airway diameter in this region, which differs from the adult nasopharynx, in which flow recirculation is obvious (see Fig. 4A).
Air flow dynamics within the nasal passages are further visualized in Figure 4B as snapshots of particle locations at selected instants. Two thousand 1.0-μm particles were released at the right nostril, and particle positions were recorded at designated time points. Due to their small inertia, 1.0-μm particles are assumed to closely follow the air flow. Faster transport and deeper penetration of aerosols are apparent in the medial passages, while slow-moving particles are found near the airway walls. Because of the dramatic airway bend from the nostrils to the nasopharynx, a high-concentration of particles constantly adjust their directions, following the mean stream-line curvature of inhaled air flow. The seemingly random particle distributions in the pharynx indicate enhanced turbulent mixing in this region. Due to the smaller nasopharynx-larynx regions in children, particles are transported faster and reach the glottal aperture faster than in the adult. This difference was more pronounced in the 3-year-old model, which had an abnormally narrow naso- and velo-pharynx.
The comparison of measured and computational-fluid-dynamics-predicted pressure drop (ΔP) for the 5 models is shown in Figure 5, for the flow range 1–48 L/min. All the models are in close agreement. The in vitro experiments and the computational-fluid-dynamics-predictions were conducted in the same airway geometry models, so direct comparison is possible. In general, the pressure drop (ie, breathing resistance) decreased as age increased (see Fig. 5). For a given flow, the infant and children models had much higher breathing resistance than the adult model. However, it is also not surprising that the 3-year-old model had a much higher pressure drop than the 7-month-old model, considering the abnormally constricted nasopharynx in the 3-year-old model and the less severe airway curvature in the 7-month-old model. The flow-pressure relationships can be expressed as a power function (ΔP = a·Qb), which can be plotted as straight lines on a log-log scale with a slope of b, as shown in the logarithmic graph in Figure 5. The best-fitted coefficients a and b for the different model ages are listed in Table 4.
Particle Deposition
Particle deposition in the nasal-laryngeal airway of the 5-year-old is displayed in Figure 6A, for the particle range 2.5–40 μm, at a flow of 10 L/min. There was highly heterogeneous particle distribution for all particle sizes considered, and the variation of deposition location with particle size is clearly shown. Particles ≥ 40 μm are mostly filtered out in the nasal vestibule (gray points in Fig. 6A). Particles around 20 μm deposit mainly in the superior and middle nasal passages, due to their higher inertia. Particles around 10 μm tend to deposit further downward, mainly in the middle passage. Particles around 2.5 μm, which can closely follow the mainstream air flow, deposit mainly in the middle and lower passages.
The deposition fraction as a function of particle aerodynamic diameter is shown in Figure 6B, in comparison with experimental data from Zhou et al.13,21 The inhalation flows considered were 10 and 20 L/min, and the particle size range was 0.5–32 μm. The in vitro experiments of Zhou et al21 and the current computational-fluid-dynamics predictions are based on the same airway geometry model. Again, there is a close match between the predicted and the measured deposition, with slight underestimation by computational fluid dynamics, compared to the experimental measurements at 20 L/min. For both flows, particles larger than 10 μm are mostly filtered out in this pediatric nasal-laryngeal airway (ie, the cut size being 10 μm). For particles smaller than 10 μm the deposition rate exhibits a high sensitivity to the aerosol size and flow, which increases quickly as the aerosol size or flow increases.
The age-related effects on nasal deposition among the 5 models are shown in Figure 7 under equivalent physical activity (ie, sedentary) conditions. The particle sizes considered are 2.5 μm and 10 μm. Large variability is observed among the 5 models, both in deposition pattern (see Fig. 7A) and regional deposition fractions (see Fig. 7B). For small particles (2.5 μm), the deposition fraction is very low in the 3 regions (ie, anterior, middle, and posterior). One exception is in the 3-year-old model, in which elevated deposition is presumably attributable to the more constricted nasopharynx. For large particles (10 μm) nasal deposition increases by approximately one order of magnitude. Again, the highest deposition was in the 3-year-old model. Considering regional deposition, the allocation of deposition fractions in the 3 regions differs for 2.5 μm and 10 μm particles in each model. For 2.5 μm particles the anterior region (nasal vestibule and valve) generally receives more particles than the others. In contrast, the middle (turbinate) region receives more deposition for 10 μm particles (see Fig. 7B). A high deposition of 10 μm particles is found in the anterior region in the adult model, but not in the infant and child models. If considering the anterior and middle regions together (nasal passage), the adult model filters more 10 μm particles than the children models. This may result from the higher complexity of the adult nasal passages.
Discussion
There is enormous physiological growth during infancy and childhood. Substantial differences were noted in the nose-throat morphology among the 5 subjects of different ages. These differences manifest themselves not only in airway dimension but also in airway morphology. For example, the nasal-laryngeal airway volumes of the 10-day-old, 7-month-old, 3-year-old, and 5-year-old were 6.4%, 18.8%, 24.2%, and 40.3% that of the 53-year-old, respectively. The 4 young subjects have smaller nostrils, a shorter turbinate region, and a narrower nasopharynx. The results of this study suggest that the nasal valve and vestibule region mature between the ages 3 and 5. This is supported by the much shorter nostril-valve distance and much smaller cross-sectional area of the nasal valve in the 2 infants, compared to the 2 older children, as well as the proximity of the these 2 parameters between the adult and the 2 older children (see Fig. 3A). Specifically, the nostril-valve distances are 8.0–11.2 mm for the 2 infants, versus 19.0–20.8 mm for the 2 children and 27.2 mm for the adult. The nasal valve areas are 42–56 mm2 for the infants, versus 120–180 mm2 for the children and adult.
Our models suggest that the turbinate region experiences fast growth from birth to the age of 5, as indicated by the remarkable volume increase of this region: 1.57 cm3 in the 10-day-old, 2.83 cm3 in the 7-month-old, 5.83 cm3 in the 3-year-old, 11.03 cm3 in the 5-year-old, and 12.63 cm3 in the 53-year-old (see Table 2). However, a lack of similarity in shape between the 5-year-old and the adult may still indicate an undeveloped turbinate region at the age of 5. It is apparent that the nasopharynx grows the least, compared to other parts of the respiratory anatomy, during the first 5 years of life: the nasopharynx volume was 0.48 cm3 in the 10-day-old, 1.74 cm3 in the 7-month-old, 1.07 cm3 in the 3-year-old, 3.95 cm3 in the 5-year-old, and 16.33 cm3 in the 53-year-old (see Table 2). Further studies on more age groups are necessary to quantify the developmental respiratory morphology and its effects on breathing and aerosol filtration.
A power function (ΔP = a·Qb) has been used to correlate the breathing frequency and airway pressure drop in previous studies.22,23 The flow-pressure relation shows high sensitivity to age in this study, as indicated by the large variation of “a” for different age groups (5.86–0.21 from 10-D-F to 53-Y-M, Table 4). For an equivalent physical activity, the inhalation rates are different for the different ages, due to their particular respiratory variables such as frequency, inspiratory-expiratory ratio, and tidal volume. Based on the selected flows (3.8 L/min for the 10-day-old, 6.5 L/min for the 7-month-old, 9.5 L/min for the 3-year-old, 11.2 L/min for the 5-year-old, and 18 L/min for the 53-year-old), the nasal breathing resistance is highest for the 10-day-old model, and decreases with age in normal subjects. Furthermore, this decrease appears most dramatic in the first year (10-day-old vs 7-month-old model herein), and the rate of decrease becomes gradually smaller beyond this age. Interestingly, a similar observation in breathing frequencies was reported by Fleming et al,15 who studied 3,381 children and found a constant decline in breathing frequency from birth to adolescence, with the steepest fall in infants under 2 years of age.
The flow-pressure relationship also shows high sensitivity to airway abnormalities, as suggested by the abrupt increase in magnitude of slope b in the 3-year-old, who had a more constricted nasopharynx (see Table 4). In this case, the subject may expect increased respiratory effort when awake and obstructive sleep apnea symptoms during sleep. The results of this study indicate that computational modeling could be a supplemental tool to study breathing-related disorders. By providing detailed air flow information that is not readily measured by conventional diagnosis techniques,24 additional clues may be unveiled that are clinically relevant to breathing disorders such as snoring and sleep apnea.
Deposition patterns of inhaled aerosols are not uniform in human upper airways; there is heterogeneous deposition in the different anatomic regions.17,25,26 The regional deposition allocations we observed varied substantially among the different age models. Considering that tissues receiving higher depositions of toxicants are more vulnerable to injury, such results might help to elucidate the diverse etiology and symptoms of respiratory disorders in different subjects or age groups. The deposition in the 3 regions considered (vestibule-valve, turbinate, and nasopharynx) was remarkably different among the 5 models, indicating a different level of risks upon the region of interest, even when exposed to the same environment. The higher turbinate deposition in the 3-year-old model may suggest a higher chance of nasal inflammation and exacerbate existing symptoms of difficult breathing. On the other hand, outcomes from inhalation therapies require drugs to be delivered to the targeted tissues at a sufficient dose. The heterogeneity in regional deposition among the different ages implies that existing adult deposition results might not guarantee an accurate dosimetry planning for infants and children.
Limitations
Limitations of this study include the assumptions of steady flow, quiet breathing only, rigid airway walls, idealized particles, and a limited number of samples per age. Other studies have highlighted the physical importance of tidal breathing,27 effects of physical activities,27 airway wall motion,28 and nasal valve collapse.29 Environmental aerosols are mostly non-spherical,30 interacting among themselves,31 and undergo size changes due to hygroscopic effects32 or coagulation.33 Moreover, each model in this study was based on images of one subject, which does not account for inter-subject variability which can be substantial.34,35 Another limitation is the typical supine position of the subjects during data acquisition, which is different from sedentary breathing. Images acquired at the end of the inhalation may not reflect variations in airway geometry during a full breathing cycle. Therefore, future studies are needed that should be oriented toward improving physical realism and including a broader population group. Our knowledge of nasal deposition is currently lacking in subpopulations such as pediatrics, geriatrics, and patients with respiratory diseases. Due to physiological development, aging, or disease states, the airway morphology can be remarkably different from that of a healthy adult. Concentrating on these specific subpopulations will help to clarify inter-group and inter-individual variability and will allow for the design of more efficient pharmaceutical formulations and drug delivery protocols for different age groups. In addition, further methodology developments and deposition measurements are required to provide robust estimates of airway depositions of either airborne contaminants or inhaled pharmacologic particles.
Conclusions
Substantial variability exists in airway morphology, air flow dynamics, and aerosol deposition among subjects of different ages. Specific observations include:
The nasal airway volumes of the 10-day-old, 7-month-old, 3-year-old, and 5-year-old were 6.4%, 18.8%, 24.2%, and 40.3% that of the 53-year-old adult, respectively, and the nasal airway surface areas were 22.7%, 30.0%, 41.7%, and 65.7% that of the adult, respectively.
The airway pressure drop is sensitive to age and airway abnormalities. Flow-pressure correlations have been proposed for different age groups based on a wide range of breathing conditions.
Age effects are large in both total and regional aerosol deposition, and should be considered in future environmental health assessment and inhaled drug delivery.
Satisfactory agreements between computational-fluid-dynamics predictions and in vitro experiments were obtained in pressure drop and particle deposition, indicating that the image/computational-fluid-dynamics coupled method is a practical tool in diagnosing respiratory disorders and developing effective inhalation devices.
Footnotes
- Correspondence: Jinxiang Xi PhD, Department of Mechanical and Biomedical Engineering, Central Michigan University, 1200 South Franklin Street, Mount Pleasant MI 48858. E-mail: xi1j{at}cmich.edu.
Dr Berlinski has disclosed relationships with Johnson & Johnson, MPEX Pharmaceutical, Gilead, Philips, Genentech, Vertex, Abvie, and S&T Technologies. None of their products are discussed in this paper. The other authors have disclosed no conflicts of interest.
Dr Xi presented a version of this paper at the OPEN FORUM of the AARC Congress 2012, held November 10–13, 2012, in New Orleans, Louisiana; for this research he was awarded the 2012 Monaghan-Trudell Fellowship for Aerosol Technique Development from the American Respiratory Care Foundation.
- Copyright © 2014 by Daedalus Enterprises