Andrew P. Kuprat

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.

Modeling biofilms with dual extracellular electron transfer mechanisms

Electrochemically active biofilms have a unique form of respiration in which they utilize solid external materials as terminal electron acceptors for their metabolism. Currently, two primary mechanisms have been identified for long-range extracellular electron transfer (EET): a diffusion- and a conduction-based mechanism. Evidence in the literature suggests that some biofilms, particularly Shewanella oneidensis, produce the requisite components for both mechanisms. In this study, a generic model is presented that incorporates the diffusion- and the conduction-based mechanisms and allows electrochemically active biofilms to utilize both simultaneously. The model was applied to S. oneidensis and Geobacter sulfurreducens biofilms using experimentally generated data found in the literature. Our simulation results show that (1) biofilms having both mechanisms available, especially if they can interact, may have a metabolic advantage over biofilms that can use only a single mechanism; (2) the thickness of G. sulfurreducens biofilms is likely not limited by conductivity; (3) accurate intrabiofilm diffusion coefficient values are critical for current generation predictions; and (4) the local biofilm potential and redox potential are two distinct parameters and cannot be assumed to have identical values. Finally, we determined that simulated cyclic and squarewave voltammetry based on our model are currently not capable of determining the specific percentages of extracellular electron transfer mechanisms in a biofilm. The developed model will be a critical tool for designing experiments to explain EET mechanisms.

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.

Modeling microstructure evolution in three dimensions with Grain3D and LaGriT

This paper will describe modeling microstructure evolution using a combination of our gradient-weighted moving finite elements code, Grain3D and our 3-D unstructured grid generation and optimization code, LaGriT. Grain boundaries evolve by mean curvature motion, and Grain3D allows for the incorporation of grain boundary orientation dependence modeled as anisotropic mobility and energy. We also describe the process of generating an initial computational grid from images obtained from electron backscatter diffraction. We present the grid optimization operations developed to respond to changes in the physical topology such as the collapse of grains and to maintain uniform computational grid quality. For 3-D columnar microstructures, validation of the method is demonstrated by comparison with experiments. For large systems of fully 3-D microstructures, simulations compare favorably to the parabolic law of normal grain growth.

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.

An automated self-similarity analysis of the pulmonary tree of the Sprague-Dawley rat.

We present the results of an automated analysis of the morphometry of the pulmonary airway trees of the Sprague–Dawley rat. Our work is motivated by a need to inform lower‐dimensional mathematical models to prescribe realistic boundary conditions for multiscale hybrid models of rat lung mechanics. Silicone casts were made from three age‐matched, male Sprague–Dawley rats, immersed in a gel containing a contrast agent and subsequently imaged with magnetic resonance (MR). From a segmentation of this data, we extracted a connected graph, representing the airway centerline. Segment statistics (lengths and diameters) were derived from this graph. To validate this MR imaging/digital analysis method, airway segment measurements were compared with nearly 1,000 measurements collected by hand using an optical microscope from one of the rat lung casts. To evaluate the reproducibility of the MR imaging/digital analysis method, two lung casts were each imaged three times with randomized orientations in the MR bore. Diameters and lengths of randomly selected airways were compared among each of the repeated imaging datasets to estimate the variability. Finally, we analyzed the morphometry of the airway tree by assembling individual airway segments into structures that span multiple generations, which we call branches. We show that branches not segments are the fundamental repeating unit in the rat lung and develop simple mathematical relationships describing these structures for the entire lung. Our analysis shows that airway diameters and lengths have both a deterministic and stochastic character. Anat Rec, 2008.

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.

Visualization of high resolution spatial mass spectrometric data during acquisition

Mass Spectrometric Imaging (MSI) allows the generation of 2D ion density maps that help visualize molecules present in sections of tissues and cells. The combination of spatial resolution and mass resolution results in very large and complex data sets. New capabilities are necessary for efficient analysis and interpretation of this data. This work details the development and application of the capability to process, visualize, query, and analyze spatial mass spectrometry data. Applications include the generation of 2D maps for selected spectra, the manipulation of the heat maps, and the identification of spectral peaks. Heat maps are generated by projecting the sum of intensity vs. time spectra of each pixel for selected m/z value or range. These capabilities take the form of a new interactive software toolkit, MSI QuickView. This software approach is a significant advance over the previous state-of-the art methods that required the conversion of the RAW data using one software, manual assembly of the data, and visualization in another software.

Linking Experimental Characterization and Computational Modeling of Grain Growth in Al-Foil

Experimental results on grain boundary properties and grain growth obtained using the Electron Backscattered Diffraction (EBSD) technique are compared with the Finite Element simulation results of an Al-foil with a columnar grain structure. The starting microstructure and grain boundary properties are implemented as an input for the three-dimensional grain growth simulation. In the computational model, minimization of the interface energy is the driving force for the grain boundary motion. The computed evolved microstructure is compared with the final experimental microstructure, after annealing at 550°C. Good agreement is observed between the experimentally obtained microstructure and the simulated microstructure. The constitutive description of the grain boundary properties was based on a 1-parameter characterization of the variation in mobility with misorientation angle.