Development of an effective two-equation turbulence modeling approach for simulating aerosol deposition across a range of turbulence levels

Pharmaceutical aerosol systems present a significant challenge to computational fluid dynamics (CFD) modeling based on the need to capture multiple levels of turbulence, frequent transition between laminar and turbulent flows, anisotropic turbulent particle dispersion, and near-wall particle transport phenomena often within geometrically complex systems over multiple time scales. Two-equation turbulence models, such as the k - ω family of approximations, offer a computationally efficient solution approach, but are known to require the use of near-wall (NW) corrections and eddy interaction model (EIM) modifications for accurate predictions of aerosol deposition. The objective of this study was to develop an efficient and effective two-equation turbulence modeling approach that enables accurate predictions of pharmaceutical aerosol deposition across a range of turbulence levels. Key systems considered were the traditional aerosol deposition benchmark cases of a 90-degree bend (R e = 6,000 ) and a vertical straight section of pipe (R e = 10,000 ), as well as a highly complex case of direct-to-infant (D2I) nose-to-lung pharmaceutical aerosol delivery from an air-jet dry powder inhaler (DPI) including a patient interface and infant nasal geometry through mid-trachea (500 < R e < 7,000 ). Of the k - ω family of models, the low Reynolds number (LRN) shear stress transport (SST) approach was determined to provide the best agreement with experimental aerosol deposition data in the D2I system, based on an improved simulation of turbulent jet flow that frequently occurs in DPIs. Considering NW corrections, a new correlation was developed to quantitatively predict best regional values of the y + l i m i t , within which anisotropic NW turbulence is approximated. Considering EIM modifications, a previously described drift correction approach was implemented in pharmaceutical aerosol simulations for the first time. Considering all model corrections and modifications applied to the D2I system, regional relative errors in deposition fractions between CFD predictions and new experimental data were improved from 19-207% (no modifications) to 2-15% (all modifications) with a notable decrease in computational time (up to ∼15%). In conclusion, the highly efficient two-equation k - ω models with physically realistic corrections and modifications provided a viable, efficient and accurate approach to simulate the transport and deposition of pharmaceutical aerosols in complex airway systems that include laminar, turbulent and transitional flows.

Establishing quantitative relationships between changes in nasal spray in vitro metrics and drug delivery to the posterior nasal region

Nasal sprays are typically characterized using in vitro spray metrics such as spray cone angle and droplet size distribution. It is currently not clear how these in vitro metrics correlate with regional nasal deposition, and these relationships could help explain the impact of product differences. In this study, the effects of changes in spray cone angle, spray velocity, spray ovality and droplet size distribution on regional nasal deposition were analyzed using a validated computational fluid dynamics model in recently developed adult characteristic nasal airway anatomies. The impact of the spray on the surrounding air phase was included. Results indicated that changes in spray cone angle largely influenced the nasal posterior deposition (PD) of the drug. Changes in the plume ovality and characteristic droplet size moderately influenced PD, but the results were dependent on the insertion conditions and nasal geometry. Changes in spray velocity and uniformity constant of the droplet size distribution had only minimal influence on PD. The rank order of metrics having the greatest to least impact on PD was cone angle ≫ plume ovality ≫ characteristic droplet size ≫ velocity ≫ size distribution uniformity constant. Overall, results from this study established quantitative relationships for predicting expected changes in PD.

Importance of Spray–Wall Interaction and Post-Deposition Liquid Motion in the Transport and Delivery of Pharmaceutical Nasal Sprays

Nasal sprays, which produce relatively large pharmaceutical droplets and have high momentum, are primarily used to deliver locally acting drugs to the nasal mucosa. Depending on spray pump administration conditions and insertion angles, nasal sprays may interact with the nasal surface in ways that creates complex droplet–wall interactions followed by significant liquid motion after initial wall contact. Additionally, liquid motion can occur after deposition as the spray liquid moves in bulk along the nasal surface. It is difficult or impossible to capture these conditions with commonly used computational fluid dynamics (CFD) models of spray droplet transport that typically employ a deposit-on-touch boundary condition. Hence, an updated CFD framework with a new spray–wall interaction (SWI) model in tandem with a post-deposition liquid motion (PDLM) model was developed and applied to evaluate nasal spray delivery for Flonase and Flonase Sensimist products. For both nasal spray products, CFD revealed significant effects of the spray momentum on surface liquid motion, as well as motion of the surface film due to airflow generated shear stress and gravity. With Flonase, these factors substantially influenced the final resting place of the liquid. For Flonase Sensimist, anterior and posterior liquid movements were approximately balanced over time. As a result, comparisons with concurrent in vitro experimental results were substantially improved for Flonase compared with the traditional deposit-on-touch boundary condition. The new SWI-PDLM model highlights the dynamicenvironment that occurs when a nasal spray interacts with a nasal wall surface and can be used to better understand the delivery of current nasal spray products as well as to develop new nasal drug delivery strategies with improved regional targeting.

Characterizing the Effects of Nasal Prong Interfaces on Aerosol Deposition in a Preterm Infant Nasal Model

The objective of this study was to characterize the effects of multiple nasal prong interface configurations on nasal depositional loss of pharmaceutical aerosols in a preterm infant nose-throat (NT) airway model. Benchmark in vitro experiments were performed in which a spray-dried powder formulation was delivered to a new preterm NT model with a positive-pressure infant air-jet dry powder inhaler using single- and dual-prong interfaces. These results were used to develop and validate a computational fluid dynamics (CFD) model of aerosol transport and deposition in the NT geometry. The validated CFD model was then used to explore the NT depositional characteristic of multiple prong types and configurations. The CFD model highlighted a turbulent jet effect emanating from the prong(s). Analysis of NT aerosol deposition efficiency curves for a characteristic particle size and delivery flowrate (3 µm and 1.4 L/min (LPM)) revealed little difference in NT aerosol deposition fraction (DF) across the prong insertion depths of 2-5 mm (DF = 16-24%) with the exception of a single prong with 5-mm insertion (DF = 36%). Dual prongs provided a modest reduction in deposition vs. a single aerosol delivery prong at the same flow for insertion depths < 5 mm. The presence of the prongs increased nasal depositional loss by absolute differences in the range of 20-70% compared with existing correlations for ambient aerosols. In conclusion, the use of nasal prongs was shown to have a significant impact on infant NT aerosol depositional loss prompting the need for prong design alterations to improve lung delivery efficiency.

In vitro evaluation of regional nasal drug delivery using multiple anatomical nasal replicas of adult human subjects and two nasal sprays

Quantifying drug delivery to the site of action using locally-acting nasal suspension sprays is a challenging but important step toward understanding bioequivalence (BE) between test and reference products. The main objective of this study was to investigate the in vitro deposition pattern of two common but different locally-acting nasal suspension sprays using multiple nasal cavities. Twenty anatomically accurate nasal replicas were developed from high-resolution sinonasal computed tomography scans of adults with healthy nasal airways. The airways were segmented into two regions of anterior and posterior to the internal nasal valve. Both sides of the septum were considered separately; hence, 40 nasal cavities were studied. The positioning of the spray nozzle in all 40 cavities was characterized by the head angle, coronal angle, and the insertion depth. Despite using a controlled protocol to minimize the anterior losses, a wide range of variability in posterior drug delivery was observed. The observed intersubject variability using this in vitro method may have important implications for understanding BE of locally-acting nasal suspension sprays.

Mice-to-men comparison of inhaled drug-aerosol deposition and clearance

Part of the effective prediction of the pharmacokinetics of drugs (or toxic particles) requires extrapolation of experimental data sets from animal studies to humans. As the respiratory tracts of rodents and humans are anatomically very different, there is a need to study airflow and drug-aerosol deposition patterns in lung airways of these laboratory animals and compare them to those of human lungs. As a first step, interspecies computational comparison modeling of inhaled nano-to-micron size drugs (50 nm < d<15μm) was performed using mouse and human upper airway models under realistic breathing conditions. Critical species-specific differences in lung physiology of the upper airways and subsequently in local drug deposition were simulated and analyzed. In addition, a hybrid modeling methodology, combining Computational Fluid-Particle Dynamics (CF-PD) simulations with deterministic lung deposition models, was developed and predicted total and regional drug-aerosol depositions in lung airways of both mouse and man were compared, accounting for the geometric, kinematic and dynamic differences. Interestingly, our results indicate that the total particle deposition fractions, especially for submicron particles, are comparable in rodent and human respiratory models for corresponding breathing conditions. However, care must be taken when extrapolating a given dosage as considerable differences were noted in the regional particle deposition pattern. Combined with the deposition model, the particle retention and clearance kinetics of deposited nanoparticles indicates that the clearance rate from the mouse lung is higher than that in the human lung. In summary, the presented computer simulation models provide detailed fluid-particle dynamics results for upper lung airways of representative human and mouse models with a comparative analysis of particle lung deposition data, including a novel mice-to-men correlation as well as a particle-clearance analysis both useful for pharmacokinetic and toxicokinetic studies.

Nanomedicine for Treatment of Acute Lung Injury and Acute Respiratory Distress Syndrome

Acute lung injury and acute respiratory distress syndrome (ARDS) represent a heterogenous group of lung disease in critically ill patients that continues to have high mortality. Despite the increased understanding of the molecular pathogenesis of ARDS, specific targeted treatments for ARDS have yet to be developed. ARDS represents an unmet medical need with an urgency to develop effective pharmacotherapies. Multiple promising targets have been identified that could lead to the development of potential therapies for ARDS; however, they have been limited because of difficulty with the mode of delivery, especially in critically ill patients. Nanobiotechnology is the basis of innovative techniques to deliver drugs targeted to the site of inflamed organs, such as the lungs. Nanoscale drug delivery systems have the ability to improve the pharmacokinetics and pharmacodynamics of agents, allowing an increase in the biodistribution of therapeutic agents to target organs and resulting in improved efficacy with reduction in drug toxicity. Although attractive, delivering nanomedicine to lungs can be challenging as it requires sophisticated systems. Here we review the potential of novel nanomedicine approaches that may prove to be therapeutically beneficial for the treatment of this devastating condition.

Computationally efficient analysis of particle transport and deposition in a human whole-lung-airway model. Part II: Dry powder inhaler application

Pulmonary drug delivery is becoming a favored route for administering drugs to treat both lung and systemic diseases. Examples of lung diseases include asthma, cystic fibrosis and chronic obstructive pulmonary disease (COPD) as well as respiratory distress syndrome (ARDS) and pulmonary fibrosis. Special respiratory drugs are administered to the lungs, using an appropriate inhaler device. Next to the pressurized metered-dose inhaler (pMDI), the dry powder inhaler (DPI) is a frequently used device because of the good drug stability and a minimal need for patient coordination. Specific DPI-designs and operations greatly affect drug-aerosol formation and hence local lung deposition. Simulating the fluid-particle dynamics after use of a DPI allows for the assessment of drug-aerosol deposition and can also assist in improving the device configuration and operation. In Part I of this study a first-generation whole lung-airway model (WLAM) was introduced and discussed to analyze particle transport and deposition in a human respiratory tract model. In the present Part II the drug-aerosols are assumed to be injected into the lung airways from a DPI mouth-piece, forming the mouth-inlet. The total as well as regional particle depositions in the WLAM, as inhaled from a DPI, were successfully compared with experimental data sets reported in the open literature. The validated modeling methodology was then employed to study the delivery of curcumin aerosols into lung airways using a commercial DPI. Curcumin has been implicated to possess high therapeutic potential as an antioxidant, anti-inflammatory and anti-cancer agent. However, efficacy of curcumin treatment is limited because of the low bioavailability of curcumin when ingested. Hence, alternative drug administration techniques, e.g., using inhalable curcumin-aerosols, are under investigation. Based on the present results, it can be concluded that use of a DPI leads to low lung deposition efficiencies because large amounts of drugs are deposited in the oral cavity. Hence, the output of a modified DPI has been evaluated to achieve improved drug delivery, especially needed when targeting the smaller lung airways. This study is the first to utilize CF-PD methodology to simulate drug-aerosol transport and deposition under actual breathing conditions in a whole lung model, using a commercial dry-powder inhaler for realistic inlet conditions.

Computationally efficient analysis of particle transport and deposition in a human whole-lung-airway model. Part I: Theory and model validation

Computational predictions of aerosol transport and deposition in the human respiratory tract can assist in evaluating detrimental or therapeutic health effects when inhaling toxic particles or administering drugs. However, the sheer complexity of the human lung, featuring a total of 16 million tubular airways, prohibits detailed computer simulations of the fluid-particle dynamics for the entire respiratory system. Thus, in order to obtain useful and efficient particle deposition results, an alternative modeling approach is necessary where the whole-lung geometry is approximated and physiological boundary conditions are implemented to simulate breathing. In Part I, the present new whole-lung-airway model (WLAM) represents the actual lung geometry via a basic 3-D mouth-to-trachea configuration while all subsequent airways are lumped together, i.e., reduced to an exponentially expanding 1-D conduit. The diameter for each generation of the 1-D extension can be obtained on a subject-specific basis from the calculated total volume which represents each generation of the individual. The alveolar volume was added based on the approximate number of alveoli per generation. A wall-displacement boundary condition was applied at the bottom surface of the first-generation WLAM, so that any breathing pattern due to the negative alveolar pressure can be reproduced. Specifically, different inhalation/exhalation scenarios (rest, exercise, etc.) were implemented by controlling the wall/mesh displacements to simulate realistic breathing cycles in the WLAM. Total and regional particle deposition results agree with experimental lung deposition results. The outcomes provide critical insight to and quantitative results of aerosol deposition in human whole-lung airways with modest computational resources. Hence, the WLAM can be used in analyzing human exposure to toxic particulate matter or it can assist in estimating pharmacological effects of administered drug-aerosols. As a practical WLAM application, the transport and deposition of asthma drugs from a commercial dry-powder inhaler is discussed in Part II.

Nanoparticle mass transfer from lung airways to systemic regions--Part I: Whole-lung aerosol dynamics

This is a two-part paper describing inhaled nanoparticle (NP) transport and deposition in a model of a human respiratory tract (Part I) as well as NP-mass transfer across barriers into systemic regions (Part II). Specifically, combining high-resolution computer simulation results of inhaled NP deposition in the human airways (Part I) with a multicompartmental model for NP-mass transfer (Part II) allows for the prediction of temporal NP accumulation in the blood and lymphatic systems as well as in organs. An understanding of nanoparticle transport and deposition in human respiratory airways is of great importance, as exposure to nanomaterial has been found to cause serious lung diseases, while the use of nanodrugs may have superior therapeutic effects. In Part I, the fluid-particle dynamics of a dilute NP suspension was simulated for the entire respiratory tract, assuming steady inhalation and planar airways. Thus, a realistic airway configuration was considered from nose/mouth to generation 3, and then an idealized triple-bifurcation unit was repeated in series and parallel to cover the remaining generations. Using the current model, the deposition of NPs in distinct regions of the lung, namely extrathoracic, bronchial, bronchiolar, and alveolar, was calculated. The region-specific NP-deposition results for the human lung model were used in Part II to determine the multicompartmental model parameters from experimental retention and clearance data in human lungs. The quantitative, experimentally validated results are useful in diverse fields, such as toxicology for exposure-risk analysis of ubiquitous nanomaterial as well as in pharmacology for nanodrug development and targeting.

Nanoparticle mass transfer from lung airways to systemic regions--Part II: Multi-compartmental modeling

This is the second article of a two-part paper, combining high-resolution computer simulation results of inhaled nanoparticle deposition in a human airway model (Kolanjiyil and Kleinstreuer, 2013, "Nanoparticle Mass Transfer From Lung Airways to Systemic Regions--Part I: Whole-Lung Aerosol Dynamics," ASME J. Biomech. Eng., 135(12), p. 121003) with a new multicompartmental model for insoluble nanoparticle barrier mass transfer into systemic regions. Specifically, it allows for the prediction of temporal nanoparticle accumulation in the blood and lymphatic systems and in organs. The multicompartmental model parameters were determined from experimental retention and clearance data in rat lungs and then the validated model was applied to humans based on pharmacokinetic cross-species extrapolation. This hybrid simulator is a computationally efficient tool to predict the nanoparticle kinetics in the human body. The study provides critical insight into nanomaterial deposition and distribution from the lungs to systemic regions. The quantitative results are useful in diverse fields such as toxicology for exposure-risk analysis of ubiquitous nanomaterial and pharmacology for nanodrug development and targeting.