Development of a CFD boundary condition to model transient vapor absorption in the respiratory airways

A Hybrid Semi-Lagrangian Cut Cell Method for Advection-Diffusion Problems with Robin Boundary Conditions in Moving Domains

We present a new discretization approach to advection-diffusion problems with Robin boundary conditions on complex, time-dependent domains. The method is based on second order cut cell finite volume methods introduced by Bochkov et al. [8] to discretize the Laplace operator and Robin boundary condition. To overcome the small cell problem, we use a splitting scheme along with a semi-Lagrangian method to treat advection. We demonstrate second order accuracy in the L ¹, L ², and L ^∞ norms for both analytic test problems and numerical convergence studies. We also demonstrate the ability of the scheme to convert one chemical species to another across a moving boundary.

Secondary indoor air pollution and passive smoking associated with cannabis smoking using electric cigarette device-demonstrative in silico study

With electronic (e)-liquids containing cannabis components easily available, many anecdotal examples of cannabis vaping using electronic cigarette devices have been reported. For electronic cigarette cannabis vaping, there are potential risks of secondary indoor air pollution from vapers. However, quantitative and accurate prediction of the inhalation and dermal exposure of a passive smoker in the same room is difficult to achieve due to the ethical constraints on subject experiments. The numerical method, i.e., in silico method, is a powerful tool to complement these experiments with real humans. In this study, we adopted a computer-simulated person that has been validated from multiple perspectives for prediction accuracy. We then conducted an in silico study to elucidate secondary indoor air pollution and passive smoking associated with cannabis vaping using an electronic cigarette device in an indoor environment. The aerosols exhaled by a cannabis vaper were confirmed to be a secondary emission source in an indoor environment; non-smokers were exposed to these aerosols via respiratory and dermal pathways. Tetrahydrocannabinol was used as a model chemical compound for the exposure study. Its uptake by the non-smoker through inhalation and dermal exposure under a worst-case scenario was estimated to be 5.9% and 2.6% of the exhaled quantity from an e-cigarette cannabis user, respectively.

How can we best mitigate unintended passive smoking in an indoor environment? As marijuana tends to be legalized in more countries, there is an increasing number of cases of vaping cannabis using e-cigarette devices. E-cigarette vaping is presumed to cause relatively low levels of indoor air pollution due to the absence of a direct combustion process. In this study, we developed a numerical simulation model to quantitatively predict the impact of first- and second-hand cannabis vaping in an indoor environment. The study was conducted in response to vulnerable residents who are concerned regarding the deterioration of indoor air quality and informs policymakers of the potential risk of second-hand cannabis vaping exposure.

Influences of puff protocols and upper airway anatomy on cannabis pharmacokinetics: A CFPD-PK study

Predicting the optimal administration doses of the inhaled Δ9-tetrahydrocannabinol (THC), i.e., one of the major natural compounds in cannabis, is critical for maximizing the therapeutic outcomes and minimizing the toxic side effects. Thus, it is essential to developing an aerosol dosimetry model to simulate the transport, deposition, and translocation of inhaled THC aerosols from the human respiratory system to the systemic region. In this study, a computational fluid-particle dynamics (CFPD) plus pharmacokinetics (PK) model was developed and validated to quantify the localized vapor and particle uptake rates of THC and the resultant THC-plasma concentrations using two human upper airway geometries. In addition, two different puff protocols (4.0/10.0 s and 1.6/11.4 s as the inhalation/holding time ratios) were employed, associated with two different inhaled THC doses (2.0 mg and 8.82 mg, respectively). The computational results demonstrated that multiple parameters had noticeable influences on THC particle deposition and vapor absorption in the upper airways, as well as the resultant pharmacokinetic behaviors. These factors include anatomical features of the upper airway, puff flow rate, duration, and holding time. The results indicated that puff protocol with 4.0/10.0 s inhalation/holding time ratio would be recommended if the treatment needs THC delivery to the deeper lung. Furthermore, the inhaled THC dose had a dominant effect on the THC-plasma PK profiles, which could override the influences of anatomical variability and puff protocols. The developed CFPD-PK modeling framework has the potential to provide localized lung absorption data and PK profiles for in vitro-in vivo correlation, as well as supporting the development and assessment of drug products containing cannabis or cannabis-derived compounds.

A numerical investigation of the potential effects of e-cigarette smoking on local tissue dosimetry and the deterioration of indoor air quality

Electronic (e)-cigarette smoking is considered to be less harmful than traditional tobacco smoking because of the lack of a combustion process. However, e-cigarettes have the potential to release harmful chemicals depending on the constituents of the vapor. To date, there has been significant evidence on the adverse health effects of e-cigarette usage. However, what is less known are the impacts of the chemicals contained in exhaled air from an e-cigarette smoker on indoor air quality, the second-hand passive smoking of residents, and the toxicity of the exhaled air. In this study, we develop a comprehensive numerical model and computer-simulated person to investigate the potential effects of e-cigarette smoking on local tissue dosimetry and the deterioration of indoor air quality. We also conducted demonstrative numerical analyses for first-hand and second-hand e-cigarette smoking in an indoor environment. To investigate local tissue dosimetry, we used newly developed physiologically based pharmacokinetic/toxicokinetic models that reproduce inhalation exposure by way of the respiratory tract and dermal exposure through the human skin surface. These models were integrated into the computer-simulated person. Our numerical simulation results quantitatively demonstrated the potential impacts of e-cigarette smoking in enclosed spaces on indoor air quality.

Can the inhalation exposure of a specific worker in a cross-ventilated factory be evaluated by time- and spatial-averaged contaminant concentration?

Industry implies economic growth; however, outdoor and indoor air pollution generated by industrial activities represents a widespread problem for the environment and human beings. In terms of human health, indoor air quality assessment has become essential in a society where people spend most of their time in indoor dwellings, as in the case of industry workers. Because indoor air quality is strongly affected by the outdoor environment, especially under natural ventilation conditions (e.g., cross-ventilation), a comprehensive analysis that includes outdoor atmospheric-urban environment is needed to reproduce realistic scenarios. In this context, computational fluid dynamics (CFD) is a useful tool. To perform a precise analysis of the inhalation exposure of factory workers to potential gas-phase contaminants in the working environment (i.e., inhaled dose of contaminants and potential effects), the human body and respiratory tract need to be integrated in the analysis. Therefore, in this study, we performed an integrated occupational inhalation exposure/toxicology assessment in a factory building that applies a computer simulated person (CSP), a virtual human respiratory tract and integrated physiologically-based toxicokinetic (PBTK) model to predict tissue dosimetry distribution. Outdoor airflow variation was transported into the enclosure through an hourly change in wind pressure coefficient to calculate transient ventilation rate and indoor contaminant concentration between 08:00 and 17:00 h. Thereafter, the time-averaged contaminant concentration calculated at the nares of the human body was employed in a steady state calculation of airflow and contaminant distribution inside the virtual respiratory tract. Subsequently, we predicted adsorbed contaminant in the first layer of tissue of the human airways; highest adsorption took place in the nasal cavity. Finally, based on the results of the comprehensive coupled numerical analysis performed using the CFD-CSP-PBTK model, we quantitatively discussed differences between the inhalation exposure concentration and representative contaminant concentration in the factory space (e.g., time and volume-averaged concentration).

Bioaccessibility and bioavailability of environmental semi-volatile organic compounds via inhalation: A review of methods and models

Semi-volatile organic compounds (SVOCs) present in indoor environments are known to cause adverse health effects through multiple routes of exposure. To assess the aggregate exposure, the bioaccessibility and bioavailability of SVOCs need to be determined. In this review, we discussed measurements of the bioaccessibility and bioavailability of SVOCs after inhalation. Published literature related to this issue is available for 2,3,7,8-tetrachlorodibenzo-p-dioxin and a few polycyclic aromatic hydrocarbons, such as benzo[a]pyrene and phenanthrene. Then, we reviewed common modeling approaches for the characterization of the gas- and particle-phase partitioning of SVOCs during inhalation. The models are based on mass transfer mechanisms as well as the structure of the respiratory system, using common computational techniques, such as computational fluid dynamics. However, the existing models are restricted to special conditions and cannot predict SVOC bioaccessibility and bioavailability in the whole respiratory system. The present review notes two main challenges for the estimation of SVOC bioaccessibility and bioavailability via inhalation in humans. First, in vitro and in vivo methods need to be developed and validated for a wide range of SVOCs. The in vitro methods should be validated with in vivo tests to evaluate human exposures to SVOCs in airborne particles. Second, modeling approaches for SVOCs need to consider the whole respiratory system. Alterations of the respiratory cycle period and human biological variability may be considered in future studies.

Molecular Binding Contributes to Concentration Dependent Acrolein Deposition in Rat Upper Airways: CFD and Molecular Dynamics Analyses

Existing in vivo experiments show significantly decreased acrolein uptake in rats with increasing inhaled acrolein concentrations. Considering that high-polarity chemicals are prone to bond with each other, it is hypothesized that molecular binding between acrolein and water will contribute to the experimentally observed deposition decrease by decreasing the effective diffusivity. The objective of this study is to quantify the probability of molecular binding for acrolein, as well as its effects on acrolein deposition, using multiscale simulations. An image-based rat airway geometry was used to predict the transport and deposition of acrolein using the chemical species model. The low Reynolds number turbulence model was used to simulate the airflows. Molecular dynamic (MD) simulations were used to study the molecular binding of acrolein in different media and at different acrolein concentrations. MD results show that significant molecular binding can happen between acrolein and water molecules in human and rat airways. With 72 acrolein embedded in 800 water molecules, about 48% of acrolein compounds contain one hydrogen bond and 10% contain two hydrogen bonds, which agreed favorably with previous MD results. The percentage of hydrogen-bonded acrolein compounds is higher at higher acrolein concentrations or in a medium with higher polarity. Computational dosimetry results show that the size increase caused by the molecular binding reduces the effective diffusivity of acrolein and lowers the chemical deposition onto the airway surfaces. This result is consistent with the experimentally observed deposition decrease at higher concentrations. However, this size increase can only explain part of the concentration-dependent variation of the acrolein uptake and acts as a concurrent mechanism with the uptake-limiting tissue ration rate. Intermolecular interactions and associated variation in diffusivity should be considered in future dosimetry modeling of high-polarity chemicals such as acrolein.

Linking Suspension Nasal Spray Drug Deposition Patterns to Pharmacokinetic Profiles: A Proof-of-Concept Study Using Computational Fluid Dynamics

The objective of this study was to link regional nasal spray deposition patterns of suspension formulations, predicted with computational fluid dynamics, to in vivo human pharmacokinetic plasma concentration profiles. This is accomplished through the use of computational fluid dynamics simulations coupled with compartmental pharmacokinetic modeling. Results showed a rapid initial rise in plasma concentration that is due to the absorption of drug particles deposited in the nasal middle passages, followed by a slower increase in plasma concentration that is governed by the transport of drug particles from the nasal vestibule to the middle passages. Although drug deposition locations in the nasal cavity had a significant effect on the shape of the concentration profile, the absolute bioavailability remained constant provided that all the drug remained in the nose over the course of the simulation. Loss of drug through the nostrils even after long periods resulted in a significant decrease in bioavailability and increased variability. The results of this study quantify how differences in nasal drug deposition affect transient plasma concentrations and overall bioavailability. These findings are potentially useful for establishing bioequivalence for nasal spray devices and reducing the burden of in vitro testing, pharmacodynamics, and clinical studies.

Small Airway Absorption and Microdosimetry of Inhaled Corticosteroid Particles after Deposition

PURPOSE

To predict the cellular-level epithelial absorbed dose from deposited inhaled corticosteroid (ICS) particles in a model of an expanding and contracting small airway segment for different particle forms.

METHODS

A computational fluid dynamics (CFD)-based model of drug dissolution, absorption and clearance occurring in the surface liquid of a representative small airway generation (G13) was developed and used to evaluate epithelial dose for the same deposited drug mass of conventional microparticles, nanoaggregates and a true nanoaerosol. The ICS medications considered were budesonide (BD) and fluticasone propionate (FP). Within G13, total epithelial absorption efficiency (AE) and dose uniformity (microdosimetry) were evaluated.

RESULTS

Conventional microparticles resulted in very poor AE of FP (0.37%) and highly nonuniform epithelial absorption, such that <5% of cells received drug. Nanoaggregates improved AE of FP by a factor of 57-fold and improved dose delivery to reach approximately 40% of epithelial cells. True nanoaerosol resulted in near 100% AE for both drugs and more uniform drug delivery to all cells.

CONCLUSIONS

Current ICS therapies are absorbed by respiratory epithelial cells in a highly nonuniform manner that may partially explain poor clinical performance in the small airways. Both nanoaggregates and nanoaerosols can significantly improve ICS absorption efficiency and uniformity.

Absorption and Clearance of Pharmaceutical Aerosols in the Human Nose: Development of a CFD Model

PURPOSE

The objective of this study was to develop a computational fluid dynamics (CFD) model to predict the deposition, dissolution, clearance, and absorption of pharmaceutical particles in the human nasal cavity.

METHODS

A three-dimensional nasal cavity geometry was converted to a surface-based model, providing an anatomically-accurate domain for the simulations. Particle deposition data from a commercial nasal spray product was mapped onto the surface model, and a mucus velocity field was calculated and validated with in vivo nasal clearance rates. A submodel for the dissolution of deposited particles was developed and validated based on comparisons to existing in vitro data for multiple pharmaceutical products. A parametric study was then performed to assess sensitivity of epithelial drug uptake to model conditions and assumptions.

RESULTS

The particle displacement distance (depth) in the mucus layer had a modest effect on overall drug absorption, while the mucociliary clearance rate was found to be primarily responsible for drug uptake over the timescale of nasal clearance for the corticosteroid mometasone furoate (MF). The model revealed that drug deposition in the nasal vestibule (NV) could slowly be transported into the main passage (MP) and then absorbed through connection of the liquid layer in the NV and MP regions. As a result, high intersubject variability in cumulative uptake was predicted, depending on the length of time the NV dose was left undisturbed without blowing or wiping the nose.

CONCLUSIONS

This study has developed, for the first time, a complete CFD model of nasal aerosol delivery from the point of spray formation through absorption at the respiratory epithelial surface. For the development and assessment of nasal aerosol products, this CFD-based in silico model provides a new option to complement existing in vitro nasal cast studies of deposition and in vivo imaging experiments of clearance.

Absorption and Clearance of Pharmaceutical Aerosols in the Human Nose: Effects of Nasal Spray Suspension Particle Size and Properties

PURPOSE

The objective of this study was to use a recently developed nasal dissolution, absorption, and clearance (DAC) model to evaluate the extent to which suspended drug particle size influences nasal epithelial drug absorption for a spray product.

METHODS

Computational fluid dynamics (CFD) simulations of mucociliary clearance and drug dissolution were used to calculate total and microscale epithelial absorption of drug delivered with a nasal spray pump. Ranges of suspended particle sizes, drug solubilities, and partition coefficients were evaluated.

RESULTS

Considering mometasone furoate as an example, suspended drug particle sizes in the range of 1-5 μm did not affect the total nasal epithelial uptake. However, the microscale absorption of suspended drug particles with low solubilities was affected by particle size and this controlled the extent to which the drug penetrated into the distal nasal regions.

CONCLUSIONS

The nasal-DAC model was demonstrated to be a useful tool in determining the nasal exposure of spray formulations with different drug particle sizes and solubilities. Furthermore, the model illustrated a new strategy for topical nasal drug delivery in which drug particle size is selected to increase the region of epithelial surface exposure using mucociliary clearance while minimizing the drug dose exiting the nasopharynx.

Comparative Risks of Aldehyde Constituents in Cigarette Smoke Using Transient Computational Fluid Dynamics/Physiologically Based Pharmacokinetic Models of the Rat and Human Respiratory Tracts

Computational fluid dynamics (CFD) modeling is well suited for addressing species-specific anatomy and physiology in calculating respiratory tissue exposures to inhaled materials. In this study, we overcame prior CFD model limitations to demonstrate the importance of realistic, transient breathing patterns for predicting site-specific tissue dose. Specifically, extended airway CFD models of the rat and human were coupled with airway region-specific physiologically based pharmacokinetic (PBPK) tissue models to describe the kinetics of 3 reactive constituents of cigarette smoke: acrolein, acetaldehyde and formaldehyde. Simulations of aldehyde no-observed-adverse-effect levels for nasal toxicity in the rat were conducted until breath-by-breath tissue concentration profiles reached steady state. Human oral breathing simulations were conducted using representative aldehyde yields from cigarette smoke, measured puff ventilation profiles and numbers of cigarettes smoked per day. As with prior steady-state CFD/PBPK simulations, the anterior respiratory nasal epithelial tissues received the greatest initial uptake rates for each aldehyde in the rat. However, integrated time- and tissue depth-dependent area under the curve (AUC) concentrations were typically greater in the anterior dorsal olfactory epithelium using the more realistic transient breathing profiles. For human simulations, oral and laryngeal tissues received the highest local tissue dose with greater penetration to pulmonary tissues than predicted in the rat. Based upon lifetime average daily dose comparisons of tissue hot-spot AUCs (top 2.5% of surface area-normalized AUCs in each region) and numbers of cigarettes smoked/day, the order of concern for human exposures was acrolein > formaldehyde > acetaldehyde even though acetaldehyde yields were 10-fold greater than formaldehyde and acrolein.

A Bidirectional Coupling Procedure Applied to Multiscale Respiratory Modeling

In this study, we present a novel multiscale computational framework for efficiently linking multiple lower-dimensional models describing the distal lung mechanics to imaging-based 3D computational fluid dynamics (CFD) models of the upper pulmonary airways in order to incorporate physiologically appropriate outlet boundary conditions. The framework is an extension of the Modified Newton's Method with nonlinear Krylov accelerator developed by Carlson and Miller [1, 2, 3]. Our extensions include the retention of subspace information over multiple timesteps, and a special correction at the end of a timestep that allows for corrections to be accepted with verified low residual with as little as a single residual evaluation per timestep on average. In the case of a single residual evaluation per timestep, the method has zero additional computational cost compared to uncoupled or unidirectionally coupled simulations. We expect these enhancements to be generally applicable to other multiscale coupling applications where timestepping occurs. In addition we have developed a "pressure-drop" residual which allows for stable coupling of flows between a 3D incompressible CFD application and another (lower-dimensional) fluid system. We expect this residual to also be useful for coupling non-respiratory incompressible fluid applications, such as multiscale simulations involving blood flow. The lower-dimensional models that are considered in this study are sets of simple ordinary differential equations (ODEs) representing the compliant mechanics of symmetric human pulmonary airway trees. To validate the method, we compare the predictions of hybrid CFD-ODE models against an ODE-only model of pulmonary airflow in an idealized geometry. Subsequently, we couple multiple sets of ODEs describing the distal lung to an imaging-based human lung geometry. Boundary conditions in these models consist of atmospheric pressure at the mouth and intrapleural pressure applied to the multiple sets of ODEs. In both the simplified geometry and in the imaging-based geometry, the performance of the method was comparable to that of monolithic schemes, in most cases requiring only a single CFD evaluation per time step. Thus, this new accelerator allows us to begin combining pulmonary CFD models with lower-dimensional models of pulmonary mechanics with little computational overhead. Moreover, because the CFD and lower-dimensional models are totally separate, this framework affords great flexibility in terms of the type and breadth of the adopted lower-dimensional model, allowing the biomedical researcher to appropriately focus on model design. Research funded by the National Heart and Blood Institute Award 1RO1HL073598.

In situ casting and imaging of the rat airway tree for accurate 3D reconstruction

The use of anatomically accurate, animal-specific airway geometries is important for understanding and modeling the physiology of the respiratory system. One approach for acquiring detailed airway architecture is to create a bronchial cast of the conducting airways. However, typical casting procedures either do not faithfully preserve the in vivo branching angles or produce rigid casts that when removed for imaging are fragile and thus easily damaged. We address these problems by creating an in situ bronchial cast of the conducting airways in rats that can be subsequently imaged in situ using three-dimensional micro-CT imaging. We also demonstrate that deformations in airway branch angles resulting from the casting procedure are small, and that these angle deformations can be reversed through an interactive adjustment of the segmented cast geometry. Animal work was approved by the Institutional Animal Care and Use Committee of Pacific Northwest National Laboratory.

Application of a new dosimetry program TAOCS to assess transient vapour absorption in the upper airways

Most previous models of vapour absorption in the respiratory tract have assumed steady state flow fields and steady state diffusion into the airway walls. However, recent studies have shown that transient absorption flux into the walls of the upper airways can significantly influence predicted uptake or deposition values. The disadvantage of accounting for transient absorption into the airway walls is a more complex boundary condition and numerical model. The objective of this study was to evaluate the effects of both transient flow fields and transient mass absorption on the uptake of highly and moderately soluble compounds in an upper airway model. The geometry consisted of the mouth-throat region coupled with a multilayer wall model containing air, mucus, tissue, and blood phases. Based on previous studies, a boundary condition that represents transient absorption into the airway walls was applied. A new dosimetry program, named transient absorption of chemical species (TAOCS) 1.0, was developed and implemented to determine the coefficients needed for the transient boundary condition expression and to apply the boundary condition to the computational fluid dynamics (CFD) model. Both steady state and transient conditions were considered for the airflow field and wall absorption. The case of perfect wall absorption with a zero surface concentration was also considered. Results indicated that steady state airflow provided a reasonable approximation to transient airflow conditions in terms of total and local deposition (values within 10-30%). However, the simulation of transient wall absorption was critical unless the compound was highly soluble (with a mucus-air partition coefficient ≥320), in which case a perfect absorption boundary condition was accurate to within a relative difference of 50%. Still, the perfect absorption boundary condition did not accurately capture local deposition enhancement factor values. Based on these findings, implementation of the transient absorption boundary condition appears critical to predict local deposition characteristics for even highly soluble compounds. Use of the TAOCS program simplified the implementation of the complex transient absorption condition making the CFD simulation process more efficient and user-friendly.