Senthil Kabilan

Comparative computational modeling of airflows and vapor dosimetry in the respiratory tracts of rat, monkey, and human

Computational fluid dynamics (CFD) models are useful for predicting site-specific dosimetry of airborne materials in the respiratory tract and elucidating the importance of species differences in anatomy, physiology, and breathing patterns. We improved the imaging and model development methods to the point where CFD models for the rat, monkey, and human now encompass airways from the nose or mouth to the lung. A total of 1272, 2172, and 135 pulmonary airways representing 17±7, 19±9, or 9±2 airway generations were included in the rat, monkey and human models, respectively. A CFD/physiologically based pharmacokinetic model previously developed for acrolein was adapted for these anatomically correct extended airway models. Model parameters were obtained from the literature or measured directly. Airflow and acrolein uptake patterns were determined under steady-state inhalation conditions to provide direct comparisons with prior data and nasal-only simulations. Results confirmed that regional uptake was sensitive to airway geometry, airflow rates, acrolein concentrations, air:tissue partition coefficients, tissue thickness, and the maximum rate of metabolism. Nasal extraction efficiencies were predicted to be greatest in the rat, followed by the monkey, and then the human. For both nasal and oral breathing modes in humans, higher uptake rates were predicted for lower tracheobronchial tissues than either the rat or monkey. These extended airway models provide a unique foundation for comparing material transport and site-specific tissue uptake across a significantly greater range of conducting airways in the rat, monkey, and human than prior CFD models.

Application of Magnetic Resonance (MR) Imaging for the Development and Validation of Computational Fluid Dynamic (CFD) Models of the Rat Respiratory System

Computational fluid dynamic (CFD) models of the respiratory system provide a quantitative basis for extrapolating the localized dose of inhaled materials and improving human health risk assessments based upon inhalation studies conducted in animals. Nevertheless, model development and validation have historically been tedious and time-consuming tasks. In recognition of this, we previously reported on the use of proton (1H) magnetic resonance (MR) imaging for visualizing nasal-sinus passages in the rat, and for speeding computational mesh generation. Here, the generation and refinement of meshes for rat nasal airways are described in more detail and simulated airflows are presented. To extend the CFD models to the complete respiratory tract, three-dimensional (3D) 1H MR imaging of rat pulmonary casts was also utilized to construct pulmonary airway meshes using procedures developed for the nasal airways. Furthermore, the feasibility of validating CFD predictions with MR was tested by imaging hyperpolarized 3He gas at physiological flow rates in a straight pipe with a diameter comparable to the rat trachea. Results from these diverse studies highlight the potential utility of MR imaging not only for speeding CFD development but also possibly for model validation.

Phase-Contrast MRI and CFD Modeling of Apparent 3He Gas Flow in Rat Pulmonary Airways

Phase-contrast (PC) magnetic resonance imaging (MRI) with hyperpolarized ³He is potentially useful for developing and testing patient-specific models of pulmonary airflow. One challenge, however, is that PC-MRI provides apparent values of local ³He velocity that not only depend on actual airflow but also on gas diffusion. This not only blurs laminar flow patterns in narrow airways but also introduces anomalous airflow structure that reflects gas-wall interactions. Here, both effects are predicted in a live rat using computational fluid dynamics (CFD), and for the first time, simulated patterns of apparent ³He gas velocity are compared with in vivo PC-MRI. Results show (1) that correlations (R²) between measured and simulated airflow patterns increase from 0.23 to 0.79 simply by accounting for apparent ³He transport, and (2) that remaining differences are mainly due to uncertain airway segmentation and partial volume effects stemming from relatively coarse MRI resolution. Higher-fidelity testing of pulmonary airflow predictions should therefore be possible with future imaging improvements.

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.

Stimuli-responsive/rheoreversible hydraulic fracturing fluids as a greener alternative to support geothermal and fossil energy production

Cost-effective yet safe creation of high-permeability reservoirs within deep bedrock is the primary challenge for the viability of enhanced geothermal systems (EGS) and unconventional oil/gas recovery. Although fracturing fluids are commonly used for oil/gas, standard fracturing methods are not developed or proven for EGS temperatures and pressures. Furthermore, the environmental impacts of currently used fracturing methods are only recently being determined. Widespread concerns about the environmental contamination have resulted in a number of regulations for fracturing fluids advocating for greener fracturing processes. To enable EGS feasibility and lessen environmental impact of reservoir stimulation, an environmentally benign, CO2-activated, rheoreversible fracturing fluid that enhances permeability through fracturing due to in situ volume expansion and gel formation is investigated herein. The chemical mechanism, stability, phase-change behavior, and rheology for a novel polyallylamine (PAA)-CO2 fracturing fluid was characterized at EGS temperatures and pressures. Hydrogel is formed upon reaction with CO2, and this process is reversible (via CO2 depressurization or solubilizing with a diluted acid) allowing potential removal from the formation and recycling, decreasing environmental impact. Rock obtained from the Coso geothermal field was fractured in laboratory-scale experiments under various EGS temperatures and pressures at significantly (at least an order of magnitude) lower effective stress than standard fracturing fluids, and the fractures were characterized with imaging, permeability measurement, and flow modeling. Although additional work is required to further understand the fluid properties, potential and limitations, this novel fracturing fluid and process represent a potential alternative to conventional fracturing fluids to vastly reduce water usage and the environmental impact of fracturing practices and effectively make EGS production and unconventional oil/gas exploitation cost-effective and cleaner.

Magnetic resonance imaging and computational fluid dynamics (CFD) simulations of rabbit nasal airflows for the development of hybrid CFD/PBPK models

The percentages of total airflows over the nasal respiratory and olfactory epithelium of female rabbits were calculated from computational fluid dynamics (CFD) simulations of steady-state inhalation. These airflow calculations, along with nasal airway geometry determinations, are critical parameters for hybrid CFD/physiologically based pharmacokinetic models that describe the nasal dosimetry of water-soluble or reactive gases and vapors in rabbits. CFD simulations were based upon three-dimensional computational meshes derived from magnetic resonance images of three adult female New Zealand White (NZW) rabbits. In the anterior portion of the nose, the maxillary turbinates of rabbits are considerably more complex than comparable regions in rats, mice, monkeys, or humans. This leads to a greater surface area to volume ratio in this region and thus the potential for increased extraction of water soluble or reactive gases and vapors in the anterior portion of the nose compared to many other species. Although there was considerable interanimal variability in the fine structures of the nasal turbinates and airflows in the anterior portions of the nose, there was remarkable consistency between rabbits in the percentage of total inspired airflows that reached the ethmoid turbinate region (~50%) that is presumably lined with olfactory epithelium. These latter results (airflows reaching the ethmoid turbinate region) were higher than previous published estimates for the male F344 rat (19%) and human (7%). These differences in regional airflows can have significant implications in interspecies extrapolations of nasal dosimetry.

Respiratory tract lung geometry and dosimetry model for male Sprague-Dawley rats.

Abstract While inhalation toxicological studies of various compounds have been conducted using a number of different strains of rats, mechanistic dosimetry models have only had tracheobronchial (TB) structural data for Long-Evans rats, detailed morphometric data on the alveolar region of Sprague-Dawley rats and limited alveolar data on other strains. Based upon CT imaging data for two male Sprague-Dawley rats, a 15-generation, symmetric typical path model was developed for the TB region. Literature data for the alveolar region of Sprague-Dawley rats were analyzed to develop an eight-generation model, and the two regions were joined to provide a complete lower respiratory tract model for Sprague-Dawley rats. The resulting lung model was used to examine particle deposition in Sprague-Dawley rats and to compare these results with predicted deposition in Long-Evans rats. Relationships of various physiologic variables and lung volumes were either developed in this study or extracted from the literature to provide the necessary input data for examining particle deposition. While the lengths, diameters and branching angles of the TB airways differed between the two Sprague-Dawley rats, the predicted deposition patterns in the three major respiratory tract regions were very similar. Between Sprague-Dawley and Long-Evans rats, significant differences in TB and alveolar predicted deposition fractions were observed over a wide range of particle sizes, with TB deposition fractions being up to 3- to 4-fold greater in Sprague-Dawley rats and alveolar deposition being significantly greater in Long-Evans rats. Thus, strain-specific lung geometry models should be used for particle deposition calculations and interspecies dose comparisons.

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

ABSTRACT 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.

CFD modeling and image analysis of exhaled aerosols due to a growing bronchial tumor: Towards non-invasive diagnosis and treatment of respiratory obstructive diseases

Diagnosis and prognosis of tumorigenesis are generally performed with CT, PET, or biopsy. Such methods are accurate, but have the limitations of high cost and posing additional health risks to patients. In this study, we introduce an alternative computer aided diagnostic tool that can locate malignant sites caused by tumorigenesis in a non-invasive and low-cost way. Our hypothesis is that exhaled aerosol distribution is unique to lung structure and is sensitive to airway structure variations. With appropriate approaches, it is possible to locate the disease site, determine the disease severity, and subsequently formulate a targeted drug delivery plan to treat the disease. This study numerically evaluated the feasibility of the proposed breath test in an image-based lung model with varying pathological stages of a bronchial squamous tumor. Large eddy simulations and a Lagrangian tracking approach were used to model respiratory airflows and aerosol dynamics. Respirations of tracer aerosols of 1 µm at a flow rate of 20 L/min were simulated, with the distributions of exhaled aerosols recorded on a filter at the mouth exit. Aerosol patterns were quantified with multiple analytical techniques such as concentration disparity, spatial scanning and fractal analysis. We demonstrated that a growing bronchial tumor induced notable variations in both the airflow and exhaled aerosol distribution. These variations became more apparent with increasing tumor severity. The exhaled aerosols exhibited distinctive pattern parameters such as spatial probability, fractal dimension, and multifractal spectrum. Results of this study show that morphometric measures of the exhaled aerosol pattern can be used to detect and monitor the pathological states of respiratory diseases in the upper airway. The proposed breath test also has the potential to locate the site of the disease, which is critical in developing a personalized, site-specific drug delivery protocol.

Anatomical Details of the Rabbit Nasal Passages and Their Implications in Breathing, Air Conditioning, and Olfaction.

The rabbit is commonly used as a laboratory animal for inhalation toxicology tests and detail knowledge of the rabbit airway morphometry is needed for outcome analysis or theoretical modeling. The objective of this study is to quantify the morphometric dimension of the nasal airway of a New Zealand white rabbit and to relate the morphology and functions through analytical and computational methods. Images of high‐resolution MRI scans of the rabbit were processed to measure the axial distribution of the cross‐sectional areas, perimeter, and complexity level. The lateral recess, which has functions other than respiration or olfaction, was isolated from the nasal airway and its dimension was quantified separately. A low Reynolds number turbulence model was implemented to simulate the airflow, heat transfer, vapor transport, and wall shear stress. Results of this study provide detailed morphological information of the rabbit that can be used in the studies of olfaction, inhalation toxicology, drug delivery, and physiology‐based pharmacokinetics modeling. For the first time, we reported a spiral nasal vestibule that splits into three paths leading to the dorsal meatus, maxilloturbinate, and ventral meatus, respectively. Both non‐dimensional functional analysis and CFD simulations suggested that the airflow in the rabbit nose is laminar and the unsteady effect is only significantly during sniffing. Due to the large surface‐to‐volume ratio, the maxilloturbinate is highly effective in warming and moistening the inhaled air to body conditions. The unique anatomical structure and respiratory airflow pattern may have important implications for designing new odorant detectors or electronic noses. Anat Rec, 299:853–868, 2016.