Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computer models”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9

Computer Modeling of Transport of Oxidizing Species in Grain Boundaries during Zirconium Corrosion

Zirconium (Zr) based alloys are widely used as the cladding materials in light-water reactors. The water-side corrosion of these alloys degrades their structural integrity and poses serious safety concerns. During the Zr corrosion process, a thin Zr oxide (ZrO2) layer forms on the alloy surface and serves as a barrier layer for further corrosion. The majority of the oxide has the monoclinic phase. At the transition region between the oxide and the metal, the oxide contains a thin layer of stabilized tetragonal phase. It is found that the texture of the tetragonal layer determines the protectiveness of the oxide for corrosion. The transport of oxidizing species, such as anion defects, cation defects, and electron through the tetragonal oxide layer could be the rate limiting step of the corrosion. The defect diffusion can be affected by the growing stresses and microstructures such as grain boundaries and dislocations. In this work molecular dynamics simulations are used to investigate the anion and cation diffusion in bulk and at grain boundaries in tetragonal ZrO2. The results show that defect diffusion at grain boundaries is complex and the behavior strongly depends on the grain boundary type. For most of the grain boundaries studied the defect diffusion are much slower than in the bulk, implying that grain boundaries may not be fast defect transport paths during corrosion. The connection between the modeling results and published experimental work will also be discussed. This work is funded by the Laboratory Directed Research and Development (LDRD) program at Idaho National Laboratory.

36 MATERIALS SCIENCE↗

High Performance Computational Models of Synthetic Jet Actuators for Step Change Improvements in Freight Efficiency

The project presents results of a computational fluid dynamics (CFD) model developed by the Oak Ridge National Laboratory team for a synthetic jet produced by an actuator. Numerical results were compared against experimental data provided by Actasys, Inc. Three CFD packages were used for this work that are Star-CCM+, Nek5000 and OpenFOAM. A sensitivity analysis was performed with Star-CCM+ and provided more insights on how to optimize the design of the actuator. Nek5000 was used to obtain a high-resolution solution of the synthetic jet. In a exploratory stage, a CFD model of the actuator was also developed in OpenFOAM. All CFD models were able to capture the main flow features.

42 ENGINEERING↗

Combining computational modeling and experimental library screening to affinity-mature VEEV-neutralizing antibody F5

Engineered monoclonal antibodies have proven to be highly effective therapeutics in recent viral outbreaks. However, despite technical advancements, an ability to rapidly adapt or increase antibody affinity and by extension, therapeutic efficacy, has yet to be fully realized. We endeavored to stand-up such a pipeline using molecular modeling combined with experimental library screening to increase the affinity of F5, a monoclonal antibody with potent neutralizing activity against Venezuelan Equine Encephalitis Virus (VEEV), to recombinant VEEV (IAB) E1E2 antigen. We modeled the F5/E1E2 binding interface and generated predictions for mutations to improve binding using a Rosetta-based approach and dTERMen, an informatics approach. The modeling was complicated by the fact that a high-resolution structure of F5 is not available and the H3 loop of F5 exceeds the length for which current modeling approaches can determine a unique structure. A subset of the predicted mutations from both methods were incorporated into a phage display library of scFvs. This library and a library generated by error-prone PCR were screened for binding affinity to the recombinant antigen. Results from the screens identified favorable mutations which were incorporated into 12 human-IgG1 variants. The best variant, containing eight mutations, improved KD from 0.63 nM (parental) to 0.01 nM. While this did not improve neutralization or therapeutic potency of F5 against IAB, it did increase cross-reactivity to other closely related VEEV epizootic and enzootic strains, demonstrating the potential of this method to rapidly adapt existing therapeutics to emerging viral strains.

affinity-maturation↗

DFO Computational Modeling Project - Catalytic Upgrading of Bio-based Furfural to 1,5-Pentanediol: A New Renewable Monomer for the Coatings Industry (Final Report)

This report summarizes the results of a collaborative effort between Oak Ridge National Laboratory (ORNL) and Pyran™ Inc. to utilize modeling capabilities developed by the CCPC (Consortium for Computational Physics and Chemistry) to assist Pyran in scaling up its proprietary process for thermocatalytic conversion of furfural to 1,5 pentanediol (PDO). Pyran’s bio-based PDO is a direct replacement for petroleum-based PDO and 1,6 hexanediol (HDO) currently used in polymers, additives, coatings, adhesives and sealants. The current market is in excess of $\$$1B/yr. According to Pyran, their production process results in 95+% reduction in fossil CO 2 at lower production costs compared to the existing petroleum-based routes for PDO & HDO. Other chemicals based on intermediates from this process have current markets in excess of $\$$10B/yr.

36 MATERIALS SCIENCE↗

Comparison of realistic and idealized breathing patterns in computational models of airflow and vapor dosimetry in the rodent upper respiratory tract

Abstract Context: Computational fluid dynamics (CFD) simulations of airflows coupled with physiologically based pharmacokinetic (PBPK) modeling of respiratory tissue doses of airborne materials have traditionally used either steady-state inhalation or a sinusoidal approximation of the breathing cycle for airflow simulations despite their differences from normal breathing patterns. Objective: Evaluate the impact of realistic breathing patterns, including sniffing, on predicted nasal tissue concentrations of a reactive vapor that targets the nose in rats as a case study. Materials and methods: Whole-body plethysmography measurements from a free-breathing rat were used to produce profiles of normal breathing, sniffing and combinations of both as flow inputs to CFD/PBPK simulations of acetaldehyde exposure. Results: For the normal measured ventilation profile, modest reductions in time- and tissue depth-dependent areas under the curve (AUC) acetaldehyde concentrations were predicted in the wet squamous, respiratory and transitional epithelium along the main airflow path, while corresponding increases were predicted in the olfactory epithelium, especially the most distal regions of the ethmoid turbinates, versus the idealized profile. The higher amplitude/frequency sniffing profile produced greater AUC increases over the idealized profile in the olfactory epithelium, especially in the posterior region. Conclusions: The differences in tissue AUCs at known lesion-forming regions for acetaldehyde between normal and idealized profiles were minimal, suggesting that sinusoidal profiles may be used for this chemical and exposure concentration. However, depending upon the chemical, exposure system and concentration and the time spent sniffing, the use of realistic breathing profiles, including sniffing, could become an important modulator for local tissue dose predictions.

dosimetry↗

Preliminary Phase Field Computational Model Development

This interim report presents progress towards the development of meso-scale models of magnetic behavior that incorporate microstructural information. Modeling magnetic signatures in irradiated materials with complex microstructures (such as structural steels) is a significant challenge. The complexity is addressed incrementally, using the monocrystalline Fe (i.e., ferrite) film as model systems to develop and validate initial models, followed by polycrystalline Fe films, and by more complicated and representative alloys. In addition, the modeling incrementally addresses inclusion of other major phases (e.g., martensite, austenite), minor magnetic phases (e.g., carbides, FeCr precipitates), and minor nonmagnetic phases (e.g., Cu precipitates, voids). The focus of the magnetic modeling is on phase-field models. The models are based on the numerical solution to the Landau-Lifshitz-Gilbert equation. From the computational standpoint, phase-field modeling allows the simulation of large enough systems that relevant defect structures and their effects on functional properties like magnetism can be simulated. To date, two phase-field models have been generated in support of this work. First, a bulk iron model with periodic boundary conditions was generated as a proof-of-concept to investigate major loop effects of single versus polycrystalline bulk iron and effects of single non-magnetic defects. More recently, to support the experimental program herein using iron thin films, a new model was generated that uses finite boundary conditions representing surfaces and edges. This model has provided key insights into the domain structures observed in magnetic force microscopy (MFM) measurements. Simulation results for single crystal thin-film iron indicate the feasibility of the model for determining magnetic domain wall thickness and mobility in an externally applied field. Because the phase-field model dimensions are limited relative to the size of most specimens used in experiments, special experimental methods were devised to create similar boundary conditions in the iron films. Preliminary MFM studies conducted on single and polycrystalline iron films with small sub-areas created with focused ion beam have correlated quite well qualitatively with phase-field simulations. However, phase-field model dimensions are still small relative to experiments thus far. We are in the process of increasing the size of the models and decreasing specimen size so both have identical dimensions. Ongoing research is focused on validation of the phase-field model. Validation is being accomplished through comparison with experimentally obtained MFM images (in progress), and planned measurements of major hysteresis loops and first order reversal curves. Extrapolation of simulation sizes to represent a more stochastic bulk-like system will require sampling of various simulations (i.e., with single non-magnetic defect, single magnetic defect, single grain boundary, single dislocation, etc.) with distributions of input parameters. These outputs can then be compared to laboratory magnetic measurements and ultimately to simulate magnetic Barkhausen noise signals.

36 MATERIALS SCIENCE↗

GEO2D - Two-Dimensional Computer Model of a Ground Source Heat Pump System

This file contains a zipped file that contains many files required to run GEO2D. GEO2D is a computer code for simulating ground source heat pump (GSHP) systems in two-dimensions. GEO2D performs a detailed finite difference simulation of the heat transfer occurring within the working fluid, the tube wall, the grout, and the ground. Both horizontal and vertical wells can be simulated with this program, but it should be noted that the vertical wall is modeled as a single tube. This program also models the heat pump in conjunction with the heat transfer occurring. GEO2D simulates the heat pump and ground loop as a system. Many results are produced by GEO2D as a function of time and position, such as heat transfer rates, temperatures and heat pump performance. On top of this information from an economic comparison between the geothermal system simulated and a comparable air heat pump systems or a comparable gas, oil or propane heating systems with a vapor compression air conditioner. The version of GEO2D in the attached file has been coupled to the DOE heating and cooling load software called ENERGYPLUS. This is a great convenience for the user because heating and cooling loads are an input to GEO2D. GEO2D is a user friendly program that uses a graphical user interface for inputs and outputs. These make entering data simple and they produce many plotted results that are easy to understand. In order to run GEO2D access to MATLAB is required. If this program is not available on your computer you can download the program MCRInstaller.exe, the 64 bit version, from the MATLAB website or from the Geothermal Data Repository. This is a free download which will enable you to run GEO2D.

15 GEOTHERMAL ENERGY↗

Noise analysis of genome-scale protein synthesis using a discrete computational model of translation

Noise in genetic networks has been the subject of extensive experimental and computational studies. However, very few of these studies have considered noise properties using mechanistic models that account for the discrete movement of ribosomes and RNA polymerases along their corresponding templates (messenger RNA (mRNA) and DNA). The large size of these systems, which scales with the number of genes, mRNA copies, codons per mRNA, and ribosomes, is responsible for some of the challenges. Additionally, one should be able to describe the dynamics of ribosome exchange between the free ribosome pool and those bound to mRNAs, as well as how mRNA species compete for ribosomes. We developed an efficient algorithm for stochastic simulations that addresses these issues and used it to study the contribution and trade-offs of noise to translation properties (rates, time delays, and rate-limiting steps). The algorithm scales linearly with the number of mRNA copies, which allowed us to study the importance of genome-scale competition between mRNAs for the same ribosomes. We determined that noise is minimized under conditions maximizing the specific synthesis rate. Moreover, sensitivity analysis of the stochastic system revealed the importance of the elongation rate in the resultant noise, whereas the translation initiation rate constant was more closely related to the average protein synthesis rate. We observed significant differences between our results and the noise properties of the most commonly used translation models. Overall, our studies demonstrate that the use of full mechanistic models is essential for the study of noise in translation and transcription.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A computational model to generate simulated three-dimensional breast masses

Purpose: To develop algorithms for creating realistic three-dimensional (3D) simulated breast masses and embedding them within actual clinical mammograms. The proposed techniques yield high-resolution simulated breast masses having randomized shapes, with user-defined mass type, size, location, and shape characteristics. Methods: The authors describe a method of producing 3D digital simulations of breast masses and a technique for embedding these simulated masses within actual digitized mammograms. Simulated 3D breast masses were generated by using a modified stochastic Gaussian random sphere model to generate a central tumor mass, and an iterative fractal branching algorithm to add complex spicule structures. The simulated masses were embedded within actual digitized mammograms. The authors evaluated the realism of the resulting hybrid phantoms by generating corresponding left- and right-breast image pairs, consisting of one breast image containing a real mass, and the opposite breast image of the same patient containing a similar simulated mass. The authors then used computer-aided diagnosis (CAD) methods and expert radiologist readers to determine whether significant differences can be observed between the real and hybrid images. Results: The authors found no statistically significant difference between the CAD features obtained from the real and simulated images of masses with either spiculated or nonspiculated margins. Likewise, the authors found that expert human readers performed very poorly in discriminating their hybrid images from real mammograms. Conclusions: The authors’ proposed method permits the realistic simulation of 3D breast masses having user-defined characteristics, enabling the creation of a large set of hybrid breast images containing a well-characterized mass, embedded within real breast background. The computational nature of the model makes it suitable for detectability studies, evaluation of computer aided diagnosis algorithms, and teaching purposes.

60 APPLIED LIFE SCIENCES↗

Progress on Computational Modeling of Water-Based NSTF (FY2021)

This report documents the FY21 progress and achievements made in the computational analyses of the water-based NSTF. Both system-level and high-fidelity Computational Fluid Dynamics (CFD) analyses were performed to gain a complete understanding of the complex flow and heat transfer phenomena in natural convection systems. The progress on the waterbased NSTF experimental testing is summarized in a companion report (ANL-ART-230). As a continuation of progress from previous years, in FY21 the RELAP5 model of the NSTF was first updated to include a representative heated cavity that incorporated convection and conduction means of heat transfer. Two different cavity models were developed, one that models the radiation heat exchange among all surfaces in one enclosure (Cavity Model 1), and another that employs multiple radiation enclosures with each considering the surfaces at one axial level (Cavity Model 2). This updated RELAP5 deck was benchmarked with single-phase test data to tune form loss coefficients of the elbows and tees, along with heat transfer coefficients within the cavity and off the external insulation panels. This tuned model was then found to accurately predict fluid temperatures within 2.5% and both system and riser flow rates within 6% during single-phase, steady-state simulations. For transient simulations, there is a slight over prediction of the flow rate, fluid temperatures, and heater temperatures during the heat-up period, which indicate that the transient heat loss may be underestimated. Following, the capabilities of RELAP5-3D were examined, which feature a conduction enclosure model that is not present in the current version of RELAP5-MOD3.3 used for all previous NSTF simulations. As an initial comparison of the two codes, single-phase simulations were run that used similar input decks. The results indicate that the single-phase flow predicted by the two codes was very similar, however RELAP5-3D appears to predict a higher system flow during the two-phase transient when compared to the flow rates predicted by RELAP5-MOD3.3. Furthermore, RELAP5-3D run did not capture the density wave oscillations or geysering phenomena, which are large oscillation flow instabilities observed in the experimental facility. To further the prediction capabilities for capturing two-phase instabilities, Options 55, 57, 58, 61, which are implemented by default in RELAP5-MOD3.3 but not in RELAP5-3D, were investigated. When enabling Option 61, RELAP5-3D agrees well with MOD3.3 in predicting the total system flow, while other options have negligible effect on the predicted total system flow by RELAP5-3D. However, none of these options enables RELAP5-3D to produce the same flow instabilities. For future steps, efforts are currently underway to investigate RELAP5-3D’s different behavior in predicting low-pressure two-phase instabilities. STAR-CCM+ was used to perform CFD analysis of the water-based test facility and provide simulation results for comparisons to experimental data. Efforts first focused on describing the characteristics of flow within the heated cavity. A study was performed on whether possible leakage flow between the two sides of the cavity could substantially impact and explain the results. An additional second assumption was examined which defined adiabatic side walls, so physical insulation structures were added to the models to allow conduction across this outer insulation. The results suggest that potential porosity of the insulation could be a major factor in the temperature distribution within the cavity and should be considered in any studies where it is desirable to know the cavity wall temperatures to higher accuracy.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

High Fidelity Computational Model for Fluidized Bed Experiments

MFiX was integrated with preconditioners and linear solver packages in Trilinos via MFiX, Fortran, C, and C++ wrappers. The MFiX wrapper interpreted the structure of the matrix and vector from MFiX while the Fortran wrapper transfered this information to a low-level language, C. The C-wrapper transformed the memory semantic between Fortran and C++ language. The C++ wrapper passed the matrix and vector to a Trilinos package and got the solution. This solution was transferred to MFiX via C++, C, Fortran, MFiX wrapper. The framework was used to integrate first, second, and third generation linear solvers in Trilinos with MFiX.

20 FOSSIL-FUELED POWER PLANTS↗

Computational modeling of charge hopping dynamics along a disordered one-dimensional wire with energy gradients in quantum environments

This computational study investigates the effects of energy gradients on charge hopping dynamics along a one-dimensional chain of discrete sites coupled to quantum bath, which is modeled at the level of Pauli master equation (PME). This study also assesses the performance of different approximations for the hopping rates. Three different methods for solving the PME, a fourth order Runge–Kutta method, numerical diagonalization of the rate matrix followed by analytic propagation, and kinetic Monte Carlo simulation method, are tested and confirmed to produce virtually identical values of time dependent mean square displacement, diffusion constant, and mobility. Five different rate expressions, exact numerical evaluation of Fermi’s Golden Rule (FGR) rate, stationary phase interpolation (SPI) approximation, semiclassical approximation, classical Marcus rate, and Miller–Abrahams rate, are tested to help understand the effects of approximations in representing quantum environments in the presence of energy gradients. The results based on direct numerical evaluation of FGR rate exhibit transition from diffusive to non-diffusive behavior with the increase in the gradient and show that the charge transport in the quantum bath is more sensitive to the magnitude of the gradient and the disorder than in the classical bath. Among all the four approximations for the hopping rates, the SPI approximation is confirmed to work best overall. A comparison of two different methods to calculate the mobility identifies drift motion of the population distribution as the major source of non-diffusive behavior and provides more reliable information on the contribution of quantum bath.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Web-Based Computational Model Builder for Nuclear Energy Advanced Modeling & Simulation (Phase I Final Report)

Today, the importance of nuclear power in the US is a geopolitical matter as much as it is an economic one. This is due to the notion that increasing the use of nuclear power reduces the US's reliance on imported oil and gas. In order to design future systems and continue operational improvements at existing US plants, advanced modeling and simulation of nuclear power reactors is crucial. A variety of challenges exist to using advanced modeling and simulation in practice, especially for the nuclear energy community and other organizations with fewer computational resources and limited expertise. There is an increasing need of access to highly capable compute, data, and software resources from network systems, which are becoming more administratively restricted from access through installation of traditional client-server applications. The nuclear energy community is not capable of maintaining exponential growth of computing capability and data storage coupled with the necessary end-to-end network bandwidth improvements. Applications that provide higher fidelity with increased physics and decreased assumptions, by their very nature, require more and more computational resources. We believe Web-based, Cloud-enabled, end-to-end simulation software for nuclear energy advanced modeling and simulation to be the solution to these challenges.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Discovering Active Subspaces for High-Dimensional Computer Models

Dimension reduction techniques have long been an important topic in statistics, and active subspaces (AS) have received much attention this past decade in the computer experiments literature. The most common approach towards estimating the AS is to use Monte Carlo with numerical gradient evaluation. While sensible in some settings, this approach has obvious drawbacks. Recent research has demonstrated that active subspace calculations can be obtained in closed form, conditional on a Gaussian process (GP) surrogate, which can be limiting in high-dimensional settings for computational reasons. In this paper, we produce the relevant calculations for a more general case when the model of interest is a linear combination of tensor products. These general equations can be applied to the GP, recovering previous results as a special case, or applied to the models constructed by other regression techniques including multivariate adaptive regression splines (MARS). Furthermore, using a MARS surrogate has many advantages including improved scaling, better estimation of active subspaces in high dimensions and the ability to handle a large number of prior distributions in closed form. In one real-world example, we obtain the active subspace of a radiation-transport code with 240 inputs and 9,372 model runs in under half an hour.

97 MATHEMATICS AND COMPUTING↗

Localizing Clinical Patterns of Blast Traumatic Brain Injury Through Computational Modeling and Simulation

Blast traumatic brain injury is ubiquitous in modern military conflict with significant morbidity and mortality. Yet the mechanism by which blast overpressure waves cause specific intracranial injury in humans remains unclear. Reviewing of both the clinical experience of neurointensivists and neurosurgeons who treated service members exposed to blast have revealed a pattern of injury to cerebral blood vessels, manifested as subarachnoid hemorrhage, pseudoaneurysm, and early diffuse cerebral edema. Additionally, a seminal neuropathologic case series of victims of blast traumatic brain injury (TBI) showed unique astroglial scarring patterns at the following tissue interfaces: subpial glial plate, perivascular, periventricular, and cerebral gray-white interface. The uniting feature of both the clinical and neuropathologic findings in blast TBI is the co-location of injury to material interfaces, be it solid-fluid or solid-solid interface. This motivates the hypothesis that blast TBI is an injury at the intracranial mechanical interfaces. In order to investigate the intracranial interface dynamics, we performed a novel set of computational simulations using a model human head simplified but containing models of gyri, sulci, cerebrospinal fluid (CSF), ventricles, and vasculature with high spatial resolution of the mechanical interfaces. Simulations were performed within a hybrid Eulerian—Lagrangian simulation suite (CTH coupled via Zapotec to Sierra Mechanics). Because of the large computational meshes, simulations required high performance computing resources. Twenty simulations were performed across multiple exposure scenarios—overpressures of 150, 250, and 500 kPa with 1 ms overpressure durations—for multiple blast exposures (front blast, side blast, and wall blast) across large variations in material model parameters (brain shear properties, skull elastic moduli). All simulations predict fluid cavitation within CSF (where intracerebral vasculature reside) with cavitation occurring deep and diffusely into cerebral sulci. These cavitation events are adjacent to high interface strain rates at the subpial glial plate. Larger overpressure simulations (250 and 500kPa) demonstrated intraventricular cavitation—also associated with adjacent high periventricular strain rates. Additionally, models of embedded intraparenchymal vascular structures—with diameters as small as 0.6 mm—predicted intravascular cavitation with adjacent high perivascular strain rates. The co-location of local maxima of strain rates near several of the regions that appear to be preferentially damaged in blast TBI (vascular structures, subpial glial plate, perivascular regions, and periventricular regions) suggest that intracranial interface dynamics may be important in understanding how blast overpressures leads to intracranial injury.

59 BASIC BIOLOGICAL SCIENCES↗

Advanced High-Performance Computational Modeling of the Seismic Response of High-Hazard and/or Nuclear Facilities and Critical Infrastructure at the NNSS

New methods for predicting the amplitude and variability of ground shaking from earthquakes (and explosions) are needed for seismic hazard analysis for buildings, nuclear power plants, and critical infrastructure at the NNSS. We are comparing existing 1-D and new 3-D geophysical methods for estimating the shear-wave velocity structure in the upper 30 meters of the ground surface (Vs30), which plays a major role in ground motion amplification and seismic response of buildings. We evaluate the performance of these methodologies at the U1a Complex at the NNSS and develop simple 1-D and high-resolution 3-D Vs30 models. We then emplace these high-resolution models into a background seismic velocity model. We will collaborate with Lawrence Livermore National Laboratory (LLNL) to conduct numerical modeling of the ground shaking at the NNSS using their high-performance computing technology and state-of-the-art ground motion simulation methodology. The primary work that was completed in FY 2019 was to acquire the seismic systems and familiarize staff at the NNSS with their use. We also worked on developing a collection plan with the Device Assembly Facility (DAF) at the NNSS, but due to time constraints and other ongoing projects at the DAF, we had to use U1a as a backup. We were able to coordinate the seismic survey, and we will complete the collection of seismic data in FY 2020. Additionally, during FY 2019, we completed the geologic framework model (GFM) for the U1a Complex and modeled the Yucca fault. LLNL worked with us through FY 2019 to prepare the data files for their modeling software and tested the software for reliability. In FY 2020 we will develop the end-to-end capability so that any facility could easily be modeled and the expected shaking from a local earthquake understood. The work in FY 2020 will include building fault models from the GFM and finalizing the velocity model analysis. The final simulations will be run for multiple rupture models, and final assessments will demonstrate the seismic hazard at the U1a Complex.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗