SEARCH · Search NASA

Results for “model potential calculations”

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.

SourceResearch connectionFrom yearThrough year

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

At least 181 records · Page 10

All sources NASA NTRS DOE OSTI

NASA NTRS2000

Simulated Active Region Emission and Dynamics: A STEREO Perspective

We present detailed three-dimensional simulations of active regions resulting from numerical models of sunspots coupled to coronal excitation and emission. The models incorporate a fully three-dimensional magnetoconvection calculation, described in a poster by Hurlburt and Rucklidge, potential field extrapolations from the sunspot model boundary conditions, steady-state coronal loops powered by the convective motions at the surface, EUV and X-ray instrument response functions, and a full voxel rendering. The result is a simulated dynamical active region in three dimensions which enables us to explore coronal heating and its relationship to the dynamics of the photosphere and convection zone. The 3D rendering of the resulting EUV emission allows us to investigate the expected coronal signatures of the EUVI instruments on board the twin STEREO spacecraft to be launched in 2004. The hybrid model developed here also provides a simulation testbed for the development of future STEREO image reconstruction tools: an integral component in the access to STEREO data by the solar physics community.

Alexander, David↗

DOE OSTI2022

Proteome-scale Structure Prediction Data - Pseudodesulfovibrio mercurii

The number of proteins predicted for Pseudodesulfovibrio mercurii is 3,446, each of which have five predicted structures from an AlphaFold run, as well as structural alignment results using the TMscore-based structural alignment method within the APoc program. Specifically, AlphaFold outputs the atoms and coordinates of the protein model in human-readable PDB files and quantitative prediction metrics in Python PICKLE files. The 5 models have been ranked based on the predicted TM-score (pTMS), a quantitative confidence metric output by AlphaFold that reports on protein model quality. The top ranked model has undergone an energy minimization calculation to relax and remove any potential clashes in the atomic coordinates. Structural alignment results are stored in two files for each protein; the top ranked model (as discussed above) is used for all alignment analyses. Both are compressed gzip files that, once unpacked, are human readable. The first file is the TMalign score results and contains the quantitative metrics for the top alignments between the predicted structure and experimental structures from the PDB70, a curated non-redundant database of about 80,000 experimental structures developed by the Soding lab. Each data point in this file is directly associated with one experimental structure; PDB ID and brief meta-data about the protein taken from the PDB70 file are reported alongside the quantitative metrics. The second results file contains the raw results associated with each alignment reported in the score results file. Specifically, the translation and rotation arrays for each alignment are provided so that the structural alignment can be recreated. Additionally, residue-level scores are reported to quantify the closeness of the aligned residues between the predicted and experimental models.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

Multiphysics Coupling Methods for Molten Salt Reactor Modeling and Simulation in VERA

New multiphysics coupling capabilities for molten salt reactor (MSR) analysis have been developed in the Virtual Environment for Reactor Applications (VERA). This development consisted of two main efforts. First, a generic species transport module was added in the CTF code, which is the thermal-hydraulics (TH) code for VERA. This module uses the velocity fields for which CTF solves during the TH calculation to transport species through the core and around the primary loop. Additionally, a gas sparging model has been added to CTF to model the movement of certain species, namely, fission products such as xenon gas, to transport between the molten salt and gas bubbles present in the salt. The second effort in this development was coupling this capability to VERA’s neutron transport code MPACT. This effort focused on coupling the detailed TH transport models in CTF to MPACT to account for feedback effects in the neutron transport calculations. Finally, the thermochemistry code Thermochimica has also been coupled to VERA. Thermochimica performs pointwise calculations for chemical potential and Gibbs free energy and determines what phases are produced by the temperature, pressure, and elemental concentrations at different locations in the primary loop. These capabilities are demonstrated using a model of the Molten Salt Reactor Experiment (MSRE). Furthermore, this reactor operated at Oak Ridge National Laboratory in the 1960s, providing sources of experimental data that were used to develop the model. Various combinations of species were modeled using VERA’s new multiphysics coupling capabilities. Species distributions and reactivity effects behaved as expected for the MSRE model, demonstrating that the coupling is behaving correctly and causing appropriate feedback. The results of these calculations show the potential for VERA to be used for a wide variety of MSR analyses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS2018

Modeling the Interactive Effects of Spruce Beetle Infestation and Climate on Subalpine Vegetation

In the subalpine zone of the Rocky Mountains, climate change is predicted to result in an increase in the frequency and severity of spruce beetle outbreaks. Climate change itself may affect vegetation, potentially leading to changes in species composition. The direct and indirect effects of climate and disturbances on forest composition, biomass, and dynamics open the possibility for non-linear ecosystem responses. Modeling studies allow for the study of the interaction of these effects and their impact on the forest system. University of Virginia Forest Model Enhanced (UVAFME), an individual-based gap model that simulates forest dynamics and characteristics, is updated with a spruce beetle subroutine that calculates the probability for beetle infestation and potential mortalityof each tree on a plot. The updated model is then run with multiple scenarios that combine beetle infestation with current or altered climate at sites across the southern Rocky Mountains. Results show that spruce beetle infestations acted to facilitate competition with invading lower-elevation species, resulting in an increase in the biomass of historically lower elevation species and a further decline in Engelmann spruce biomass than occurred with solely bark beetle disturbance or solely climate change. We also found an initial enhancing effect between spruce beetle infestation and climate change; however, by the end of 100 yr of climate change and potential beetle infestation, climate had a dampening effect on spruce beetle infestation, through loss of host trees. These results are an important step in understanding the possible futures for vegetation of the Rocky Mountains as well as for spruce forests across the western United States and Canada.

Foster, Adrianna C.↗

DOE OSTI2025

Microscopic optical potentials from a Green's function approach

Optical potentials are a standard tool in the study of nuclear reactions, as they describe the interaction between a target nucleus and a projectile. The use of phenomenological optical potentials built using experimental data on stable isotopes is widespread. Although successful in their dedicated domain, it is unclear whether these phenomenological potentials can provide reliable predictions for unstable isotopes. To address this problem, optical potentials based on microscopic nuclear structure input calculations prove to be crucial and are an important current line of research. In this work we present an explicit implementation of the Feshbach formalism for the systematic derivation of optical potentials using input from nuclear structure models. Numerical tools for the derivation of Green's functions associated with nonlocal potentials are presented. In conclusion, the new optical potential, based on the valence shell model, is applied to the calculations of 𝑛 + 24 Mg elastic scattering and yields a close agreement with the experimental data.

Direct reactions↗

NASA NTRS1986

Optical model calculations of heavy-ion target fragmentation

The fragmentation of target nuclei by relativistic protons and heavy ions is described within the context of a simple abrasion-ablation-final-state interaction model. Abrasion is described by a quantum mechanical formalism utilizing an optical model potential approximation. Nuclear charge distributions of the excited prefragments are calculated by both a hypergeometric distribution and a method based upon the zero-point oscillations of the giant dipole resonance. Excitation energies are estimated from the excess surface energy resulting from the abrasion process and the additional energy deposited by frictional spectator interactions of the abraded nucleons. The ablation probabilities are obtained from the EVA-3 computer program. Isotope production cross sections for the spallation of copper targets by relativistic protons and for the fragmenting of carbon targets by relativistic carbon, neon, and iron projectiles are calculated and compared with available experimental data.

Townsend, L. W.↗

DOE OSTI2020

Weak Transitions in Light Nuclei

Nuclei are used for high-precision tests of the Standard Model and for studies of physics beyond the Standard Model. Without a thorough understanding of nuclei, we will not be able to meaningfully interpret the growing body of experimental data nor will we be able to disentangle new physics signals from underlying nuclear effects. This calls for accurate calculations of nuclear structure and reactions. In this work, we focus on electroweak decays in nuclei with mass number A ≤ 10 and report on ab initio Quantum Monte Carlo calculations of reduced matrix elements entering beta decays and electron captures in nuclei with mass number A ≤ 10. The many-body wave functions are calculated using selected Norfolk two- and three-nucleon potential models and associated one- and two-body axial currents at tree-level obtained from a chiral effective field theory with pions, nucleons, and Δ. The agreement with the experimental data is satisfactory except for transitions in A = 8 nuclei. In this specific case, the theory significantly underpredicts the experimental data, which indicates the need of further improvements in the corresponding nuclear wave functions. In this study, emphasis is placed on the contributions of two-body axial currents that are carefully analyzed using two-body transition densities. This allow us to study the spatial distribution and short-range behavior of two-body dynamics. In particular, the transition densities when scaled to peak at 1.0 exhibit universal short-range behavior across the considered nuclei, while they differ in the long-range tails.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Adsorbate Free Energies from DFT-Derived Translational Energy Landscapes

Adsorption free energies are fundamental to surface chemistry and catalysis. Standard models combine some assumed analytical form of the translational potential energy surface, often parameterized against density functional theory (DFT) calculations, with an analytical expression for the resultant translational densities of states (DOS), free energy, and entropy. Here we compare the performance of such models against numerical evaluations of the DOS and thermodynamic functions derived from solutions to the translational Schrodinger equation. We compare results for a translational potential energy surface (PES) derived from nudged eleastic band calculations with those obtained from adsorbate rastering across a series of monatomic (O, S, C, N, H) and polyatomic (NH x ) adsorbates on (100) Pt and Au facets. We find that analytical models as commonly parameterized have mixed performance for describing the translational PES and that the consequences for computed free energies are modest but potentially significant in microkinetic models. Here, numerical solutions are possible for modest to no additional computational cost over analytical models and thus should be considered when reliable free energy estimates are needed or translational potential energy surfaces are available.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1971

Semiclassical vibration-rotation transition probabilities for motion in molecular state averaged potentials.

Collision-induced vibration-rotation transition probabilities are calculated from a semiclassical three-dimensional model, in which the collision trajectory is determined by the classical motion in the interaction potential that is averaged over the molecular rotational state, and compared with those for which the motion is governed by a spherically averaged potential. For molecules that are in highly excited rotational states, thus dominating the vibrational relaxation rate at high temperature, it is found that the transition probability for rotational state averaging is smaller than that for spherical averaging. For typical collisions, the transition cross section is decreased by a factor of about 1.5 to 2.

Stallcop, J. R.↗

DOE OSTI2024

Atomically Dispersed Zn/Co–N–C as ORR Electrocatalysts for Alkaline Fuel Cells

Hydrogen fuel cells have drawn increasing attention as one of the most promising next-generation power sources for future automotive transportation. Developing efficient, durable, and low-cost electrocatalysts, to accelerate the sluggish oxygen reduction reaction (ORR) kinetics, is urgently needed to advance fuel cell technologies. Herein, we report on metal–organic frameworks-derived nonprecious dual metal single-atom catalysts (SACs) (Zn/Co–N–C), consisting of Co–N 4 and Zn–N 4 local structures. These catalysts exhibited superior ORR activity with a half-wave potential (E 1/2 ) of 0.938 V versus RHE (reversible hydrogen electrode) and robust stability (ΔE 1/2 = –8.5 mV) after 50k electrochemical cycles. Moreover, this remarkable performance was validated under realistic fuel cell working conditions, achieving a record-high peak power density of ~1 W cm –2 among the reported SACs for alkaline fuel cells. Operando X-ray absorption spectroscopy was conducted to identify the active sites and reveal catalytic mechanistic insights. The results indicated that the Co atom in the Co–N 4 structure was the main catalytically active center, where one axial oxygenated species binds to form an O ads –Co–N 4 moiety during the ORR. In addition, theoretical studies, based on a potential-dependent microkinetic model and core-level shift calculations, showed good agreement with the experimental results and provided insights into the bonding of oxygen species on Co–N 4 centers during the ORR. Furthermore, this work provides a comprehensive mechanistic understanding of the active sites in the Zn/Co–N–C catalysts and will pave the way for the future design and advancement of high-performance single-site electrocatalysts for fuel cells and other energy applications.

25 ENERGY STORAGE↗

DOE OSTI2024

Understanding electronic structure tunability by metal dopants for promoting MgB 2 hydrogenation

Hydrogen is a promising energy carrier, but its onboard application is limited by the need for compact, low-pressure storage solutions. Solid-state complex metal hydride systems, such as MgB 2 /Mg(BH 4 ) 2 , offer high storage capacities but suffer from sluggish kinetics and poor reversibility. One avenue for improving reactivity is to introduce metal dopants to alter electronic and atomic properties, but the role of these chemical additives remains poorly understood, particularly for the hydrogenation reaction. In this work, we used density functional theory calculations on model MgB 2 systems to rationalize the potential role of metal dopants in destabilizing B–B bonding within the MgB 2 lattice. We carried out detailed electronic structure analyses for 28 different metal dopant adatoms to identify properties that contribute to a dopant’s efficacy. Based on the simulation results, we propose that an intermediate ionic and covalent character of the bonds between adatoms and B atoms is desirable for facilitating charge redistribution, disrupting the B–B bond network, and promoting H 2 dissociation and H atom chemisorption on MgB 2 .

08 HYDROGEN↗

DOE OSTI2024

Density functional theory of material design: fundamentals and applications—II

Abstract This is the second and the final part of the review on density functional theory (DFT), referred to as DFT-II. In the first review, DFT-I, we have discussed wavefunction-based methods, their complexity, and basics of density functional theory. In DFT-II, we focus on fundamentals of DFT and their implications for the betterment of the theory. We start our presentation with the exact DFT results followed by the concept of exchange-correlation (xc) or Fermi-Coulomb hole and its relationship with xc energy functional. We also provide the exact conditions for the xc-hole, xc-energy and xc-potential along with their physical interpretation. Next, we describe the extension of DFT for non-integer number of electrons, the piecewise linearity of total energy and discontinuity of chemical potential at integer particle numbers, and derivative discontinuity of the xc potential, which has consequences on fundamental gap of solids. After that, we present how one obtains more accurate xc energy functionals by going beyond the LDA. We discuss the gradient expansion approximation (GEA), generalized gradient approximation (GGA), and hybrid functional approaches to designing better xc energy functionals that give accurate total energies. However, these functionals fail to predict properties like the ionization potential and the band gap. Thus, we next describe different methods of modelling these potentials and results of their application for calculation of the band gaps of different solids to highlight accuracy of different xc potentials. Finally, we conclude with a glimpse on orbital-free density functional theory and the machine learning approach.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Plating and Stripping of Lithium Metal Stabilized by a Block Copolymer Electrolyte: Local Current Density Measurement and Modeling

Here, plating and stripping of lithium protrusions in lithium metal symmetric cells containing a solid block copolymer electrolyte was studied as a function of time in 3D using time-resolved X-ray tomography. These measurements enabled determination of the spatial variation in current densities at the plating and stripping electrodes. The initial interelectrode distance was 27 μm. Correlation functions were calculated to reveal the relationships between current densities at the two electrodes and local electrolyte thickness. Current densities at opposing electrode locations during protrusion growth is uncorrelated until the local interelectrode distance decreases to less than 6 μm, just before the cell shorts. Mass balance was used to determine the area from which lithium ions that form a protrusion were stripped. Computational modeling of the plating and stripping process reveals the interplay between electrochemical and mechanical driving forces and their effect on nonuniform current distribution. Model predictions were compared with experiments without resorting to any adjustable parameters. The computed correlation functions were in qualitative agreement with experiments. Finally, the model was used to calculate contour plots of electrochemical potential within the electrolyte, shedding light on how geometry, salt concentration, interelectrode distance, and mechanical stress influence local rates of electrochemical reaction.

36 MATERIALS SCIENCE↗

NASA NTRS1985

Spectral constraints on models of gas in clusters of galaxies

The HEAO 1A2 spectra of clusters of galaxies are used to determine the temperature profile which characterizes the X-ray emitting gas. Strong evidence of nonisothermality is found for the Coma, A85, and A1795 clusters. Properties of the cluster potential which binds the gas are calculated for a range of model parameters. The typical binding mass, if the gas is adiabatic, is 2-4E14 solar masses and is quite centrally concentrated. In addition, the Fe abundance in Coma is .26 + or - .06 solar, less than the typical value (.5) found for rich clusters. The results for the gas in Coma may imply a physical description of the cluster which is quite different from what was previously believed.

Henriksen, M. J.↗

NASA NTRS1985

Spectral constraints on models of gas in clusters of galaxies

Henriksen, M. J.↗

DOE OSTI2024

"Source Term Modeling for Advanced Gas Micro-Reactors"

Maintaining the safety of the public, environment, and operating personnel is the most important factor in designing, operating, maintaining, and decommissioning nuclear reactors. In recent years, there has been a growing interest in the development of micro-reactors employing TRi-structural ISOtropic (TRISO)-coated particle fuel. In gas reactors, TRISO fuel plays an important role in the safety case for high temperature reactors because of the fission product retention properties of the fuel. This ability enables the use of a functional containment strategy for the reactor where multiple barriers are used to prevent fission product release to the environment. Part of the safety analysis of these advanced reactors is the assessment of radionuclide releases under normal and accident conditions through the multiple credited safety barriers. Using conservative assumptions, a mechanistic analysis can be performed to quantify these releases that combines the probabilistic assessment of failure with analytic solutions to radionuclide transport equations. Source term modeling for TRISO fuel has been performed for previous reactor designs; however, these models are outdated, in many cases proprietary, and need updates to be applied to the current state of TRISO fuel technology and alternative gas reactor core configurations [1]. Currently, the only publicly available source term assessment for gas reactors is an expert-based Monte Carlo simulation based on the effectiveness of the fuel kernel, coating layers, and graphite block in a modular high temperature gas reactor [2]. Thus, there is a need to develop a simple, versatile, and mechanistic model of fission product release and transport in gas reactor cores that could be applied to a variety of reactors through user inputs and reactor-specific radionuclide inventories. The release is calculated by the diffusion of the key safety important fission products through the kernel, silicon carbide (SiC), graphite for both intact and defective TRISO particles based on fuel and graphite temperatures in the reactor under normal operation. These releases from the fuel enter the coolant where they can plate-out on cooler surfaces. A clean-up model is included for designs with a coolant purification system to remove fission gases. This initial distribution of fission products in the reactor serves as an initial condition for potential releases under postulated accident conditions. The model then can calculate the fission product release for any transient temperature profile and fission product releases can then be used to assess radiological dose to the workers and the public using conventional dose tools. Data on the diffusion of fission products is based on historic German TRISO experiments and the more current Department of Energy (DOE) Advanced Gas Reactor (AGR) TRISO fuel development program. The model is coded in python with inputs and outputs in excel spreadsheets, as well as python plotting utilities to aid in the interpretation of the results. References: [1] INL, NGNP Mechanistic Source Term White Paper, INL-10-17997, July 2010. [2] David A. Petti, Richard R. Hobbins, Peter Lowry, Hans Gougar, “Representative Source Terms and The Influence of Reactor Attributes on Functional Containment in Modular High Temperature Gas-cooled Reactors,” Nuclear Technology, Vol. 184, p. 181-197, Nov. 2013.

07 ISOTOPE AND RADIATION SOURCES↗

NASA NTRS1978

Ames interactive molecular model building system - A 3-D computer modelling system applied to the study of the origin of life

The investigation of specific interactions among biological molecules must take into consideration the stereochemistry of the structures. Thus, models of the molecules are essential for describing the spatial organization of potentially interacting groups, and estimations of conformation are required for a description of spatial organization. Both the function of visualizing molecules, and that of estimating conformation through calculations of energy, are part of the molecular modeling system described in the present paper. The potential uses of the system in investigating some aspects of the origin of life rest on the assumption that translation of conformation from genetic elements to catalytic elements would have been required for the development of the first replicating systems subject to the process of biological evolution.

Coeckelenbergh, Y.↗

NASA NTRS1976

A computer program to calculate the longitudinal aerodynamic characteristics of wing-flap configurations with externally blown flaps

A vortex lattice lifting-surface method is used to model the wing and multiple flaps. Each lifting surface may be of arbitrary planform having camber and twist, and the multiple-slotted trailing-edge flap system may consist of up to ten flaps with different spans and deflection angles. The engine wakes model consists of a series of closely spaced vortex rings with circular or elliptic cross sections. The rings are normal to a wake centerline which is free to move vertically and laterally to accommodate the local flow field beneath the wing and flaps. The two potential flow models are used in an iterative fashion to calculate the wing-flap loading distribution including the influence of the waves from up to two turbofan engines on the semispan. The method is limited to the condition where the flow and geometry of the configurations are symmetric about the vertical plane containing the wing root chord. The calculation procedure starts with arbitrarily positioned wake centerlines and the iterative calculation continues until the total configuration loading converges within a prescribed tolerance. Program results include total configuration forces and moments, individual lifting-surface load distributions, including pressure distributions, individual flap hinge moments, and flow field calculation at arbitrary field points.

Mendenhall, M. R.↗