Multiscale Computational Models for Respiratory Aerosol Dynamics with Medical Applications
View this Special IssueResearch Article  Open Access
Sanghun Choi, Shinjiro Miyawaki, ChingLong Lin, "A Feasible Computational Fluid Dynamics Study for Relationships of Structural and Functional Alterations with Particle Depositions in Severe Asthmatic Lungs", Computational and Mathematical Methods in Medicine, vol. 2018, Article ID 6564854, 12 pages, 2018. https://doi.org/10.1155/2018/6564854
A Feasible Computational Fluid Dynamics Study for Relationships of Structural and Functional Alterations with Particle Depositions in Severe Asthmatic Lungs
Abstract
This study aims to investigate the effect of altered structures and functions in severe asthma on particle deposition by using computational fluid dynamics (CFD) models. Airway geometrical models of two healthy subjects and two severe asthmatics were reconstructed from computed tomography (CT) images. Subjectspecific flow boundary conditions were obtained by image registration to account for regional functional alterations of severe asthmatics. A large eddy simulation (LES) model for transitional and turbulent flows was applied to simulate airflows, and particle transport simulations were then performed for 2.5, 5, and 10 μm particles using CFDpredicted flow fields. Compared to the healthy subjects, the severe asthmatics had a smaller airvolume change in the lower lobes and a larger airvolume change in the upper lobes. Both severe asthmatics had smaller airway circularity (Cr), but one of them had a significant reduction of hydraulic diameter (). In severe asthmatics, the larger airvolume change in the upper lobes resulted in more particles in the upper lobes, especially for the small 2.5 μm particles. The structural alterations measured by Cr and were associated with a higher particle deposition. was found to be the most important metric which affects the specific location of particle deposition. This study demonstrates the relationship of CTbased structural and functional alterations in severe asthma with flow and particle dynamics.
1. Introduction
Asthma is pathologically characterized by combined phenotypes of airflow obstruction, bronchial hyperresponsiveness, and airway inflammation [1]. However, how structural and functional alterations of asthma affect flow structure and particle deposition is yet to be investigated. In imaging studies of asthma, ventilation defects and airway structural changes have been investigated by using magnetic resonance image (MRI), positron emission tomography (PET), and singlephoton emission computed tomography (SPECT) [2–4]. In addition, via quantitative computed tomography (QCT) imaging, several studies [5–8] have demonstrated significant alterations such as reduced airway diameter as well as increased wall thickness and air trapping. Although QCT can provide structural measurements of the airways up to segmental airways (~2 mm), quantification of local functional variables is still limited.
Image registration technique has been utilized to provide functional information by matching images at different inflation levels [9]. The registration derivedvariables were validated by comparing ventilation maps from different imaging modalities [10]. Furthermore, this technique has shown strengths when characterizing functional alterations of diseased lungs [11, 12]. For instance, a study used the technique to differentiate airway vs. parenchymal phenotypes in a chronic obstructive pulmonary disease (COPD) [11]. We have recently demonstrated that volume changes of severe asthmatic lungs are preferentially smaller near basal regions and the smaller volume changes are compensated with airvolume change in apical regions [12]. We also demonstrated that severe asthmatics were characterized by the reduced airway diameters and noncircular shapes, unlike healthy and nonsevere asthmatics [13].
Computational fluid dynamics (CFD) technique has been used to analyze flow characteristics and particle depositions in the human lungs [14]. With regard to CFD simulations of asthmatic lungs, a study [15] showed alterations of particle depositions with an asthmatic subject before and after asthma attack, and another study [16] correlated forced expiratory volume in one second (FEV_{1}) with CFDbased resistance before and after bronchodilator. However, both of them imposed parabolic velocity profiles to the trachea regions and uniform pressure boundary conditions to QCTresolved ending branches, under the Reynoldsaveraged NavierStokes equations (RANS) model or laminar assumption. Heenan et al. [17] and Jayaraju et al. [18] have compared CFD simulations using RANS model with experimental results, concluding that ones with RANS provided less accurate prediction in both flow structure and particle deposition. On the other hand, Longest et al. [19] and Tian et al. [20] have demonstrated that lowReynolds number  approximation with nearwall correction improves the prediction of particle deposition. With five asthmatics, Vinchurkar et al. [21] have also demonstrated a good agreement of aerosol delivery between imagingbased CFD method and in vivo experimental data. Meanwhile, Backer et al. [22] have demonstrated that lobar distributions of airvolume change between SPECT and QCT are consistent, so they have emphasized the importance of subjectspecific boundary conditions in CFD simulations.
The main objective of this study is to investigate how functional and structural alterations of severe asthmatic lungs affect flow structures and particle depositions. This study further evaluates clinical applications of CFD in localizing hot spots. According to our existing populationbased analysis [8, 12, 13], severe asthmatics were characterized by the shifted lobar airvolume change from lower lobes to upper lobes, decreased airway circularity (Cr), and reduced hydraulic diameter (). To reproduce regional airvolume change, we employed an image registration technique for a subjectspecific boundary condition. Next, to reflect subjectspecific airway features such as Cr and , we employed a surface fitting method from 1D skeleton to 3D CTresolved airways [23]. We then performed computational fluidparticle dynamics simulations to obtain local particle distribution and deposition. To achieve main objective, three major components in this study are (1) reproducing regional airvolume change, (2) quantifying subjectspecific airway features, and (3) quantifying local particle distribution and deposition.
2. Methods
2.1. Human Subjects
The imaging study and protocols for acquiring QCT images at both total lung capacity (TLC) and functional residual capacity (FRC) were approved by Institutional Review Board (IRB) of University of Pittsburgh as a part of multicenter Severe Asthma Research Program (SARP) [24]. Four human subjects were employed in this study, among whom two subjects were healthy and two subjects were severe asthmatic (Table 1). The CT images were taken via GE helical VCT64 slice scanner with the slice thickness of 0.625 mm. Major criteria used to define severe asthma included treatments with oral corticosteroids and highdose inhaled corticosteroids, besides several minor criteria [24]. At least, one major and two minor criteria were required to be classified as severe asthmatics. The CT scans were acquired in a supine position, and 3D airways, 1D skeletons, lungs, and lobes were segmented using Apollo software (VIDA Diagnostics, Coralville, Iowa).
 
All subjects are Caucasians and nonsmokers. PFT and CT measurements were obtained in upright and supine positions, respectively. CT scans were performed after bronchodilator. 
2.2. Flow Simulation and Boundary Condition
Filtered continuity and incompressible NavierStokes equations readwhere , , , , and are filtered velocity vector, fluid density, filtered pressure, fluid kinematic viscosity, and subgridscale turbulent eddy viscosity [25], respectively. A large eddy simulation (LES) model was adopted to resolve laminartransitionalturbulent flows in the central airways. The properties of and were set to 1.12 kg/m^{3} and 1.64 × 10^{5 }m^{2}/s, respectively. A characteristic Galerkin finite element method [26] was employed to discretize (1a) and (1b). A moderate steadyinspiratory flowrate of 3.27 × 10^{4 }m^{3}/s (≈20 liters/min) was imposed as the inlet condition, being equivalent to the peak flowrate of a sinusoidal waveform with a tidal volume of 500 mL and a period of 4.8 s. Reynolds numbers (Re) in trachea range from 1300 to 1700 in the given flowrate and individual tracheal sizes of two healthy and two severe asthmatics.
Figure 1(a) shows a flow chart from the image segmentation and registration to a realistic CFD simulation. First, an image registration technique [9] was employed to estimate local airvolume change between TLC and FRC at lung parenchyma (top row), and onedimensional (1D) centerline (CL) tree structures obtained by a volume filling method [27] bridged threedimensional (3D) resolved ending branches and lung parenchyma (middle row). Airvolume change measured at parenchyma was summed at CTresolved ending branch to estimate regional flowrate ratios. With the flowrate ratios, parabolic profiles of velocity were essentially imposed at the CT ending branch. Lastly, a surface fitting method [23] together with Gmesh [28] was used to fit an initial CLbased 3D surface geometry to CTsegmented airway surface and construct a CLCTbased airway geometrical model (bottom row). The resulting model captured subjectspecific airway features such as Cr and . At this step, we created a synthetic turbulence [29] right above the glottal constricted regions, because the SARP study [24] did not gather oropharynx scans. The turbulent intensity and largest eddy size were set to 0.29 and 8 mm, respectively, to mimic the turbulent flows found in the CFD simulations with oropharynx [30]. The number of elements ranges from 3.9 to 5.0 million for four cases. Time steps of 5.0 × 10^{−6} s (HS 1, HS 2, and SA 1) and 3.0 × 10^{−6} s (SA 2) were chosen to satisfy Courant–Friedrichs–Lewy (CFL) number less than 1.
(a)
(b)
2.3. Structural Metrics Associated with Aerosol Delivery
Cr and could be directly associated with airflows and particle delivery. In this study, we measured the structural metrics at 31 segmental airways (Figure 1(b)). Cr was computed to assess the degree of elliptical shape of an airway cross section as follows:where and are average diameter and perimeter of an airway. was computed as follows: , where is the crosssectional area of the airway. Next, in order to assess a degree of airway narrowing, was computed as follows:
Note that Cr decreases as airway lumen becomes elliptic, and decreases as airway narrowing is dominant. Cr and were decreased in severe asthmatics compared with healthy subjects, in a populationbased comparison [13]. In this study, Cr decreases in both severe asthmatics (SA 1 and SA 2), but only decreases in a severe asthmatic subject (SA 2). In Results, the effects of Cr and on particle delivery are investigated.
2.4. Particle Simulation
To compare the characteristics of particle distribution and deposition such as pharmaceutical aerosols between healthy subjects and severe asthmatics, particle transport analysis was conducted using LESpredicted airflow fields. Quasisteady airflow fields for 2.4 s were collected, e.g., 400 datasets × 0.006 s, after initial 2.4 s airflow fields. For particle simulation, Lagrangian particle tracking algorithm [31] was adopted to track particle trajectories as follows:where , , and are particle velocity, particle density, and gravitational acceleration, respectively. In the equation, Stokes number (Stk) is defined as follows:where Q, d, , , and α are the flowrate of the branch, the diameter of particles, fluid dynamic viscosity, the Cunningham slip correction factor, and the particleparticle interaction factor, respectively. The detailed descriptions of particle transport simulation could be found in [30]. At the beginning of particle simulation, the spherical particles were uniformly distributed in a cylinder with a radius of 10 mm and a depth of 4 mm at the trachea inlet. They were then released 9 different times to obtain an ensemble average of particle deposition. The number of particles was set to 10,000, and three different spherical particle sizes of 2.5, 5, and 10 μm were chosen in this study. The aerosol size was chosen with a general distribution of aerosols using dry powder inhaler (DPI) and soft mist inhaler (SMI) [32]. This is to investigate the relationship between aerosol size and airway structure. The particle density and mean free path were given as 1000 kg/m^{3} and 68 nm, respectively. “Particle distribution”, “deposition”, and “advection” by lobe are defined as “particles entering each lobe”, “those deposited in 3D segments of each lobe”, and “those exiting 3D ending branches of each lobe without being deposited in 3D segments”, respectively.
3. Results
3.1. Pulmonary Function Test (PFT)
Table 1 shows demographic as well as PFT and CTbased measurements. Bronchodilator was performed for both healthy subjects and severe asthmatics to obtain maximal lung functions, and CT scans were acquired after bronchodilator because the aim of SARP study was to assess lung function of stable asthma [33]. First, the baseline and maximal FEV_{1} and FVC %predicted values in two healthy subjects were close to normal ranges (~100% predicted). In contrast, baseline FEV_{1} % predicted values of SA 1 (34%) and SA 2 (40%) were much smaller than normal ranges (≥80%). The maximal FEV_{1}% predicted value (~46%) of SA 2 was still lower even with bronchodilator. In the same subject, the maximal FVC % predicted value was very close to normal ranges, leading to a significantly reduced FEV_{1}/FVC (38%). Thus, the baseline FEV_{1}% predicted values of both SA 1 and SA 2 were small values, but the lung function of SA 1 (38%) is relatively reversible than that of SA 2 (6%). According to baseline and maximal PFTs, SA 2 is likely to have significant airway narrowing despite bronchodilator, whereas airways of SA 1 might have dilated with the aid of bronchodilator.
3.2. CTBased Structural and Functional Characteristics
We compared structural quantities of Cr and in two healthy subjects and two severe asthmatics (Figure 2). Figure 2(a) shows that Cr of two severe asthmatics was reduced compared to those of two healthy subjects. For example, Cr of Trachea, RMB, TriRUL, RLL7, RB6, RB9+10, and LMB in two severe asthmatics were deviated from those of the same airway segments in healthy subjects. Both SA 1 and SA 2 had the prominently reduced Cr, but only SA 2 demonstrated significantly smaller compared with HS 1, HS 2, and SA 1 in CTresolved airways (Figure 2(b)), which was also reflected in PFT measurements at postbronchodilator. The effect of reduced Cr in RMB and TriRUL is assessed between “HS 1, HS 2” and “SA 1, SA 2”. Next, the effect of constricted branch on particle deposition is mainly assessed with SA 2. In addition, SA 1 and SA 2 were characterized by reduced airvolume change in lower lungs along with elevated airvolume change in upper lungs, as compared with the healthy subjects, being quantified by the airvolume change in upper lobes to airvolume change in middle and lower lobes, i.e., (Figure 3).
(a)
(b)
3.3. Regional Volume Change vs. Particle Distribution
Given the regional airvolume change, particle distribution was evaluated with three different particle sizes (Figure 3). With 2.5 μm particles, the (the particle distribution ratio of upper lobes to middle and lower lobes) was close to the ratio of airvolume changes (), because small particles are likely to follow airflow streamlines. Next, decreased as particle size increases due in part to the inertial effect that particles can easily move along with the flow to the lower lobes. Overall, SA 1 and SA 2 had an increased particle distribution to upper lobes compared with HS 1 and HS 2, being consistent with the alteration of airvolume change distribution in severe asthmatics. Furthermore, the particle distribution to upper lobes in severe asthmatics became more evident with decreasing particle size.
3.4. Circularity and Particle Deposition
Around RMB and TriRUL regions (Figure 1(b)), we compared particle depositions of two healthy subjects without elliptic shapes vs. two severe asthmatics with elliptic shapes, because SA 1 and SA 2 exhibit significantly decreased Cr in both RMB and TriRUL (Table 2). The two selected severe asthmatics (83° and 88°) have similar bifurcation angles as two healthy subjects (90° and 91°) at RMB, expecting similar particle depositions given the same range of Stk numbers. In Table 2, of TriRUL in SA 1 is the largest among four subjects, but was slightly smaller than those of both HS 1 and HS 2 due to the reduced Cr. Figure 4(a) shows the increased particle deposition of SA 1 and SA 2 in TriRUL, as compared with HS 1 and HS 2 in the same Stk range of the parent branch (RMB). Note that the comparisons based on the same Stk numbers were made to control the effects of flowrate and airway narrowing. Figure 5 also displays the distributions of particle deposition in TriRUL. Thus, an increase of particle deposition in both severe asthmatics was possibly due to the decreased airway Cr.

(a)
(b)
(c)
(a)
(b)
3.5. Constriction vs. Particle Deposition
Particle deposition efficiency of the branches based on Stk numbers, daughter branches of RB9+10 and LB10, were plotted in Figures 4(b) and 4(c), respectively. Figure 6(a) (around RB9+10) and Figure 6(b) (around LB10) show particle depositions of SA 2 with narrowed and elliptic airways. Compared with HS 1, HS 2, and SA 1, Stk is increased for SA 2 due to airway narrowing, because Stk is proportional to with a given flowrate (see (5)). In the case of 5 μm particles, particle deposition efficiency reached 60% and 90% for RB9+10 and LB10, respectively. For 10 μm particles, particle deposition reached 100%; thus large particles could not ventilate beyond these airways. These results imply that particle deposition would be very sensitive to the reduced airway diameters due to airway constriction (see Figures 6(a) and 6(b)). We further investigated flow structure, wall shear stress, and pressure drop around LB10 (Figure 6(c)). In this airway, constrictioninduced jetflow, higher wall shear stress and large pressure drop were simultaneously observed.
(a)
(b)
(c)
4. Discussion
Our previous studies [8, 12] have demonstrated that airvolume change is smaller in lower lobes and larger in the upper lobes in severe asthmatics, as compared with healthy subjects. In addition, severe asthmatics were structurally characterized by smaller Cr and . In this study, we focused on the effect of structural and functional alterations in severe asthmatics on particle depositions under a moderate inspiratory flow condition (~20 liters/min). The two severe asthmatics had altered lung function characterized by larger than the two healthy subjects. Both SA 1 and SA 2 had small % predicted values in the baseline lung function and showed smaller Cr, whereas only SA 2 subject had small % predicted in the maximal lung function and showed smaller in CTresolved regions. This observation implied that bronchodilator may help recover narrowed airways, but it may not recover noncircular airway shape. This emphasizes the importance of investigating airway shape besides airway diameter in particle delivery. In this study, SA 1 and SA 2 were used to evaluate the effect of Cr on particle deposition, and SA 2 was used to evaluate the effect of on particle deposition.
First, investigating characteristics of particle deposition in severe asthmatics is important for inhaled pharmaceutical drugs [21, 34], airborne bacteria, or air pollutants. Lobar distribution of small particles was consistent with airvolume distribution (Figure 3). With increasing particle size, the delivery of particles to lower lobes may increase because of an increased inertial effect of large particles. When targeting subjectspecific treatments of aerosols, understanding the mechanism of particle distribution would be critical. For example, if lower lobar bronchi need be targeted for some specific subjects, relatively larger aerosol should be treated. However, one should also consider sizes of targeted branches, because larger aerosols could not be delivered to small airways with higher generation.
We previously found that in severe asthmatics Cr of RMB is significantly smaller [13], so we investigated the effect of smaller Cr on particle deposition. With the same Stk and similar bifurcation angle of RMB, particle depositions of SA 1 and SA 2 were greater in TriRUL than healthy subjects, which could be attributed to smaller Cr. Airway structures with smaller Cr might increase the sensitivity of bacterial inflammation due to particle deposition, which would potentially increase airway wall thickness. While both SA 1 and SA 2 demonstrated smaller Cr, only SA 2 had constricted airways, especially in lower lobes. Particle deposition efficiency is the function of based on Stk numbers, given flowrate. In other words, airway constriction could be the most important structural characteristic in association with increased particle deposition efficiency. Thus, if airways are chronically constricted in a subject, the subject may have an increased exposure of pollutants, so further reducing airway diameters along with airway inflammation. Mechanisms among airway constriction, wall shear stress, pressure drop, and particle deposition may be strongly coupled as described in Figure 6(c). First, constrictioninduced high velocity creates high wall shear stress in association with high velocity gradient on the wall. A large pressure drop is required to overcome large wall shear stress. In addition, the constrictioninduced high velocity might lead to an increase of particle deposition due to particle impaction in both constricted parent and daughter branches (Figures 6(a) and 6(b)). As a result, both high velocity and small may contribute to be a large Stk (see (5)), leading to an increase of particle deposition.
In this study, we only investigated two healthy subjects and two severe asthmatics, which could prevent statistical reliability, so a CFDbased population study is necessary. However, performing many CFD simulations to obtain normal distribution is computationally expensive. We have recently performed CT imagingbased clustering analysis and found clinically meaningful subgroups [35]. Since the clustering membership employed airway structure and lung function, each cluster has similar airway structure and lung function. We believe that such an approach using clustering analysis possibly reduces the number of samples by transitioning an interest from intersubject study to intercluster study. In this study, we imposed the same flowrate on trachea inlet, due to a lack of subjectspecific flowrate information. If the subjectspecific flowrate is measured, we would further improve the comparison between healthy subjects and severe asthmatics. Furthermore, we have recently investigated airway resistance with symmetric branching angles up to the 20th generation for inspiration and expiration [36]. Then an empirical airway resistance model was introduced to estimate pressure drop due to kinetic energy and viscous dissipation. However, we have not discussed airway resistance in this study, because the number of airways was limited to investigate airway resistance. In the future, airway resistance study of severe asthmatics with realistic multiscale airways would be desirable.
In conclusion, we applied a highfidelity CFD model together with CT imagebased airway models to study particle deposition in both healthy and severe asthmatic lungs. With the aid of image registration technique, the subjectspecific physiologically realisticflow boundary condition is derived based on airvolume difference between two CT lung images of the same human subject. As expected from flowrate distribution, particle distribution to upper lobes was larger in severe asthmatics, relative to healthy subjects. This phenomenon was more prominent when using smaller particles. In both healthy subjects and severe asthmatics, with increasing particle size, particles are distributed more toward lower lobar regions due to inertial effects. Alterations of airway Cr and were found to be associated with particle deposition. With the same Stk, reduced Cr increases particle deposition. On the other hand, reduced significantly elevates Stk, resulting in greater deposition efficiency. The constricted airways contribute to high wall shear stress, elevated pressure drop, and significantly increased particle deposition. It is concluded that , Cr, and shall be carefully considered to target subjectspecific aerosols.
Disclosure
Preliminary results of the current manuscript were presented as an abstract in the following link: http://meetings.aps.org/Meeting/DFD14/Session/E15.1.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The authors thank Drs. Eric A Hoffman, Sally E Wenzel, and Jiwoong Choi for assisting in imaging data acquisition and analysis. Computational time was supported by XSEDE and by the National Institute of Supercomputing and Network/Korea Institute of Science and Technology Information with supercomputing resources including technical support (KSC2017S10021). This study was supported in part by NIH Grants U01 HL114494, R01 HL094315, and S10 RR022421 and in part by Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (NRF2017R1D1A1B03034157).
References
 W. W. Busse and R. F. Lemanske Jr., “Asthma,” The New England Journal of Medicine, vol. 344, no. 5, pp. 350–362, 2001. View at: Publisher Site  Google Scholar
 R. S. Harris, H. FujiiRios, T. Winkler, G. Musch, M. F. Vidal Melo, and J. G. Venegas, “Ventilation Defect Formation in Healthy and Asthma Subjects Is Determined by Lung Inflation,” PLoS ONE, vol. 7, no. 12, Article ID e53216, 2012. View at: Publisher Site  Google Scholar
 G. G. King, B. Harris, and S. Mahadev, “V/Q SPECT: Utility for investigation of pulmonary physiology,” Seminars in Nuclear Medicine, vol. 40, no. 6, pp. 467–473, 2010. View at: Publisher Site  Google Scholar
 J. G. Venegas, T. Schroeder, S. Harris, R. T. Winkler, and M. F. Vidal Melo, “The distribution of ventilation during bronchoconstriction is patchy and bimodal: A PET imaging study,” Respiratory Physiology & Neurobiology, vol. 148, no. 12, pp. 57–64, 2005. View at: Publisher Site  Google Scholar
 M. Montaudon, M. Lederlin, S. Reich et al., “Bronchial measurements in patients with asthma: Comparison of quantitative thinsection CT findings with those in healthy subjects and correlation with pathologic findings,” Radiology, vol. 253, no. 3, pp. 844–853, 2009. View at: Publisher Site  Google Scholar
 A. Busacker, J. D. Newell Jr., T. Keefe et al., “A multivariate analysis of risk factors for the airtrapping asthmatic phenotype as measured by quantitative CT analysis,” CHEST, vol. 135, no. 1, pp. 48–56, 2009. View at: Publisher Site  Google Scholar
 R. S. Aysola, E. A. Hoffman, D. Gierada et al., “Airway remodeling measured by multidetector CT is increased in severe asthma and correlates with pathology,” CHEST, vol. 134, no. 6, pp. 1183–1191, 2008. View at: Publisher Site  Google Scholar
 S. Choi, E. A. Hoffman, S. E. Wenzel, M. Castro, and C.L. Lin, “Improved CTbased estimate of pulmonary gas trapping accounting for scanner and lungvolume variations in a multicenter asthmatic study,” Journal of Applied Physiology, vol. 117, no. 6, pp. 593–603, 2014. View at: Publisher Site  Google Scholar
 Y. Yin, E. A. Hoffman, and C.L. Lin, “Mass preserving nonrigid registration of CT lung images using cubic Bspline,” Medical Physics, vol. 36, no. 9, pp. 4213–4222, 2009. View at: Publisher Site  Google Scholar
 J. M. Reinhardt, K. Ding, K. Cao, G. E. Christensen, E. A. Hoffman, and S. V. Bodas, “Registrationbased estimates of local lung tissue expansion compared to xenon CT measures of specific ventilation,” Medical Image Analysis, vol. 12, no. 6, pp. 752–763, 2008. View at: Publisher Site  Google Scholar
 C. J. Galbán, M. K. Han, J. L. Boes et al., “Computed tomographybased biomarker provides unique signature for diagnosis of COPD phenotypes and disease progression,” Nature Medicine, vol. 18, no. 11, pp. 1711–1715, 2012. View at: Publisher Site  Google Scholar
 S. Choi, E. A. Hoffman, S. E. Wenzel et al., “Registrationbased assessment of regional lung function via volumetric CT images of normal subjects vs. severe asthmatics,” Journal of Applied Physiology, vol. 115, no. 5, pp. 730–742, 2013. View at: Publisher Site  Google Scholar
 S. Choi, E. A. Hoffman, S. E. Wenzel et al., “Quantitative assessment of multiscale structural and functional alterations in asthmatic populations,” Journal of Applied Physiology, vol. 118, no. 10, pp. 1286–1298, 2015. View at: Publisher Site  Google Scholar
 H. Zhang and G. Papadakis, “Computational analysis of flow structure and particle deposition in a single asthmatic human airway bifurcation,” Journal of Biomechanics, vol. 43, no. 13, pp. 2453–2459, 2010. View at: Publisher Site  Google Scholar
 K. Inthavong, J. Tu, Y. Ye, S. Ding, A. Subic, and F. Thien, “Effects of airway obstruction induced by asthma attack on particle deposition,” Journal of Aerosol Science, vol. 41, no. 6, pp. 587–601, 2010. View at: Publisher Site  Google Scholar
 J. W. De Backer, W. G. Vos, A. Devolder et al., “Computational fluid dynamics can detect changes in airway resistance in asthmatics after acute bronchodilation,” Journal of Biomechanics, vol. 41, no. 1, pp. 106–113, 2008. View at: Publisher Site  Google Scholar
 A. F. Heenan, E. Matida, A. Pollard, and W. H. Finlay, “Experimental measurements and computational modeling of the flow field in an idealized human oropharynx,” Experiments in Fluids, vol. 35, no. 1, pp. 70–84, 2003. View at: Publisher Site  Google Scholar
 S. T. Jayaraju, M. Brouns, C. Lacor, B. Belkassem, and S. Verbanck, “Large eddy and detached eddy simulations of fluid flow and particle deposition in a human mouththroat,” Journal of Aerosol Science, vol. 39, no. 10, pp. 862–875, 2008. View at: Publisher Site  Google Scholar
 P. W. Longest, M. Hindle, S. D. Choudhuri, and P. R. Byron, “Numerical simulations of capillary aerosol generation: CFD model development and comparisons with experimental data,” Aerosol Science and Technology, vol. 41, no. 10, pp. 952–973, 2007. View at: Publisher Site  Google Scholar
 G. Tian, P. W. Longest, G. Su, R. L. Walenga, and M. Hindle, “Development of a stochastic individual path (SIP) model for predicting the tracheobronchial deposition of pharmaceutical aerosols: Effects of transient inhalation and sampling the airways,” Journal of Aerosol Science, vol. 42, no. 11, pp. 781–799, 2011. View at: Publisher Site  Google Scholar
 S. Vinchurkar, L. De Backer, W. Vos, C. Van Holsbeke, J. De Backer, and W. De Backer, “A case series on lung deposition analysis of inhaled medication using functional imaging based computational fluid dynamics in asthmatic patients: Effect of upper airway morphology and comparison with in vivo data,” Inhalation Toxicology, vol. 24, no. 2, pp. 81–88, 2012. View at: Publisher Site  Google Scholar
 J. W. De Backer, W. G. Vos, S. C. Vinchurkar et al., “Validation of computational fluid dynamics in CTbased airway models with SPECT/CT,” Radiology, vol. 257, no. 3, pp. 854–862, 2010. View at: Publisher Site  Google Scholar
 S. Miyawaki, M. H. Tawhai, E. A. Hoffman, S. E. Wenzel, and C.L. Lin, “Automatic construction of subjectspecific human airway geometry including trifurcations based on a CTsegmented airway skeleton and surface,” Biomechanics and Modeling in Mechanobiology, vol. 16, no. 2, pp. 583–596, 2017. View at: Publisher Site  Google Scholar
 S. E. Wenzel and W. W. Busse, “Severe asthma: Lessons from the Severe Asthma Research Program,” The Journal of Allergy and Clinical Immunology, vol. 119, no. 1, pp. 14–21, 2007. View at: Publisher Site  Google Scholar
 A. W. Vreman, “An eddyviscosity subgridscale model for turbulent shear flow: algebraic theory and applications,” Physics of Fluids, vol. 16, no. 10, pp. 3670–3681, 2004. View at: Publisher Site  Google Scholar
 C.L. Lin, H. Lee, T. Lee, and L. J. Weber, “A level set characteristic Galerkin finite element method for free surface flows,” International Journal for Numerical Methods in Fluids, vol. 49, no. 5, pp. 521–547, 2005. View at: Publisher Site  Google Scholar
 M. Howatson Tawhai, A. J. Pullan, and P. J. Hunter, “Generation of an Anatomically Based ThreeDimensional Model of the Conducting Airways,” Annals of Biomedical Engineering, vol. 28, no. 7, pp. 793–802, 2000. View at: Publisher Site  Google Scholar
 C. Geuzaine and J. F. Remacle, “Gmsh: A 3D finite element mesh generator with builtin pre and postprocessing facilities,” International Journal for Numerical Methods in Engineering, vol. 79, no. 11, pp. 1309–1331, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 S. Miyawaki, E. A. Hoffman, and C.L. Lin, “Numerical simulations of aerosol delivery to the human lung with an idealized laryngeal model, imagebased airway model, and automatic meshing algorithm,” Computers & Fluids, vol. 148, pp. 1–9, 2017. View at: Publisher Site  Google Scholar  MathSciNet
 S. Miyawaki, M. H. Tawhai, E. A. Hoffman, and C.L. Lin, “Effect of carrier gas properties on aerosol distribution in a CTbased human airway numerical model,” Annals of Biomedical Engineering, vol. 40, no. 7, pp. 1495–1507, 2012. View at: Publisher Site  Google Scholar
 M. R. Maxey and J. J. Riley, “Equation of motion for a small rigid sphere in a nonuniform flow,” Physics of Fluids, vol. 26, no. 4, pp. 883–889, 1983. View at: Publisher Site  Google Scholar
 P. W. Longest, G. Tian, R. L. Walenga, and M. Hindle, “Comparing MDI and DPI aerosol deposition using in vitro experiments and a new stochastic individual path (SIP) model of the conducting airways,” Pharmaceutical Research, vol. 29, no. 6, pp. 1670–1688, 2012. View at: Publisher Site  Google Scholar
 R. L. Sorkness, E. R. Bleecker, W. W. Busse et al., “Lung function in adults with stable but severe asthma: Air trapping and incomplete reversal of obstruction with bronchodilation,” Journal of Applied Physiology, vol. 104, no. 2, pp. 394–403, 2008. View at: Publisher Site  Google Scholar
 R. L. Walenga and P. W. Longest, “Current Inhalers Deliver Very Small Doses to the Lower Tracheobronchial Airways: Assessment of Healthy and Constricted Lungs,” Journal of Pharmaceutical Sciences, vol. 105, no. 1, pp. 147–159, 2016. View at: Publisher Site  Google Scholar
 S. Choi, E. A. Hoffman, S. E. Wenzel et al., “Quantitative computed tomographic imaging–based clustering differentiates asthmatic subgroups with distinctive clinical phenotypes,” The Journal of Allergy and Clinical Immunology, vol. 140, no. 3, pp. 690–700.e8, 2017. View at: Publisher Site  Google Scholar
 S. Choi, J. Choi, and C.L. Lin, “Contributions of Kinetic Energy and Viscous Dissipation to Airway Resistance in Pulmonary Inspiratory and Expiratory Airflows in Successive Symmetric Airway Models with Various Bifurcation Angles,” Journal of Biomechanical Engineering, vol. 140, no. 1, 2018. View at: Google Scholar
Copyright
Copyright © 2018 Sanghun Choi et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.