Search NASA⌕ Search

SEARCH · Search NASA

Results for “Calculation”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Reactive Molecular Dynamics Simulations and Quantum Chemistry Calculations To Investigate Soot-Relevant Reaction Pathways for Hexylamine Isomers

Sooting tendencies of a series of nitrogen-containing hydrocarbons (NHCs) have been recently characterized experimentally using the yield sooting index (YSI) methodology. This work aims to identify soot-relevant reaction pathways for three selected C6H15N amines, namely, dipropylamine (DPA), diisopropylamine (DIPA), and 3,3-dimethylbutylamine (DMBA) using ReaxFF molecular dynamics (MD) simulations and quantum mechanical (QM) calculations and to interpret the experimentally observed trends. ReaxFF MD simulations are performed to determine the important intermediate species and radicals involved in the fuel decomposition and soot formation processes. QM calculations are employed to extensively search for chemical reactions involving these species and radicals based on the ReaxFF MD results and also to quantitatively characterize the potential energy surfaces. Specifically, ReaxFF simulations are carried out in the NVT ensemble at 1400, 1600, and 1800 K, where soot has been identified to form in the YSI experiment. These simulations account for the interactions among test fuel molecules and pre-existing radicals and intermediate species generated from rich methane combustion, using a recently proposed simulation framework. ReaxFF simulations predict that the reactivity of the amines decrease in the order DIPA > DPA > DMBA, independent of temperature. Both QM calculations and ReaxFF simulations predict that C2H4, C3H6, and C4H8 are the main nonaromatic soot precursors formed during the decomposition of DPA, DIPA, and DMBA, respectively, and the associated reaction pathways are identified for each amine. Both theoretical methods predict that sooting tendency increases in the order DPA, DIPA, and DMBA, consistent with the experimentally measured trend in YSI. This work demonstrates that sooting tendencies and soot-relevant reaction pathways of fuels with unknown chemical kinetics can be identified efficiently through combined ReaxFF and QM simulations. Overall, predictions from ReaxFF simulations and QM calculations are consistent, in terms of fuel reactivity, major intermediates, and major nonaromatic soot precursors.

fuel decomposition↗
DOE OSTI2015

Structures and magnetic properties of Co-Zr-B magnets studied by first-principles calculations

The structures and magnetic properties of Co-Zr-B alloys near the composition of Co{sub 5}Zr with B at. % ≤6% were studied using adaptive genetic algorithm and first-principles calculations. The energy and magnetic moment contour maps as a function of chemical composition were constructed for the Co-Zr-B magnet alloys through extensive structure searches and calculations. We found that Co-Zr-B system exhibits the same structure motif as the “Co{sub 11}Zr{sub 2}” polymorphs, and such motif plays a key role in achieving strong magnetic anisotropy. Boron atoms were found to be able to substitute cobalt atoms or occupy the “interruption” sites. First-principles calculations showed that the magnetocrystalline anisotropy energies of the boron-doped alloys are close to that of the high-temperature rhombohedral Co{sub 5}Zr phase and larger than that of the low-temperature Co{sub 5.25}Zr phase. Our calculations provide useful guidelines for further experimental optimization of the magnetic performances of these alloys.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Real-space representation of the quasiparticle self-consistent GW self-energy and its application to defect calculations

The quasiparticle self-consistent (QS) GW (G for Green's function, W for screened Coulomb interaction) approach incorporates the corrections of the quasiparticle energies from their Kohn-Sham density functional theory (DFT) eigenvalues by means of an energy-independent and Hermitian self-energy matrix usually given in the basis set of the DFT eigenstates. By expanding these into an atom-centered basis set (specifically here the linearized muffin-tin orbitals) a real space representation of the self-energy corrections becomes possible. In this work, We show that this representation is relatively short-ranged. This offers opportunities to construct the self-energy of a complex system from parts of the system by a cut-and-paste method. Specifically for a point defect, represented in a large supercell, the self-energy can be constructed from those of the host and a smaller defect-containing cell. The self-energy of the periodic host can be constructed simply from a GW calculation for the primitive cell. We show for the case of the As Ga in GaAs that the defect part can already be well represented by a minimal eight-atom cell and allows us to construct the self-energy for a 64-atom cell in good agreement with direct QSGW calculations for the large cell. Using this approach to an even larger 216-atom cell shows the defect band approaches an isolated defect level. The calculations also allow us to identify a second defect band which appears as a resonance near the conduction band minimum. The results on the extracted defect levels agree well with Green's function calculations for an isolated defect and with experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Calculating Time-Integral Quantities in Depletion Calculations

A method referred to as tally nuclides is presented for accurately and efficiently calculating the time-step averages and integrals of any quantities that are weighted sums of atomic densities with constant weights during the step. The method allows all such quantities to be calculated simultaneously as a part of a single depletion solution with existing depletion algorithms. Some examples of the results that can be extracted include step-average atomic densities and macroscopic reaction rates, the total number of fissions during the step, and the amount of energy released during the step. Furthermore, the method should be applicable with several depletion algorithms, and the integrals or averages should be calculated with an accuracy comparable to that reached by the selected algorithm for end-of-step atomic densities. The accuracy of the method is demonstrated in depletion calculations using the Chebyshev rational approximation method. Here, we demonstrate how the ability to calculate energy release in depletion calculations can be used to determine the accuracy of the normalization in a constant-power burnup calculation during the calculation without a need for a reference solution.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

Numerical calculation of protein-ligand binding rates through solution of the Smoluchowski equation using smooth particle hydrodynamics

We demonstrate the accuracy and effectiveness of a Lagrangian particle-based method, smoothed particle hydrodynamics (SPH), to study diffusion in biomolecular systems by numerically solving the time-dependent Smoluchowski equation for continuum diffusion. The numerical method is first verified in simple systems and then applied to the calculation of ligand binding to an acetylcholinesterase monomer. Unlike previous studies, a reactive Robin boundary condition (BC), rather than the absolute absorbing (Dirichlet) boundary condition, is considered on the reactive boundaries. This new boundary condition treatment allows for the analysis of enzymes with "imperfect" reaction rates. Rates for inhibitor binding to mAChE are calculated at various ionic strengths and compared with experiment and other numerical methods. We find that imposition of the Robin BC improves agreement between calculated and experimental reaction rates. Although this initial application focuses on a single monomer system, our new method provides a framework to explore broader applications of SPH in larger-scale biomolecular complexes by taking advantage of its Lagrangian particle-based nature.

Pan, Wenxiao↗
DOE OSTI2014

Low energy isomers of (H{sub 2}O){sub 25} from a hierarchical method based on Monte Carlo temperature basin paving and molecular tailoring approaches benchmarked by MP2 calculations

We report new global minimum candidate structures for the (H{sub 2}O){sub 25} cluster that are lower in energy than the ones reported previously and correspond to hydrogen bonded networks with 42 hydrogen bonds and an interior, fully coordinated water molecule. These were obtained as a result of a hierarchical approach based on initial Monte Carlo Temperature Basin Paving sampling of the cluster's Potential Energy Surface with the Effective Fragment Potential, subsequent geometry optimization using the Molecular Tailoring Approach with the fragments treated at the second order Møller-Plesset (MP2) perturbation (MTA-MP2) and final refinement of the entire cluster at the MP2 level of theory. The MTA-MP2 optimized cluster geometries, constructed from the fragments, were found to be within <0.5 kcal/mol from the minimum geometries obtained from the MP2 optimization of the entire (H{sub 2}O){sub 25} cluster. In addition, the grafting of the MTA-MP2 energies yields electronic energies that are within <0.3 kcal/mol from the MP2 energies of the entire cluster while preserving their energy rank order. Finally, the MTA-MP2 approach was found to reproduce the MP2 harmonic vibrational frequencies, constructed from the fragments, quite accurately when compared to the MP2 ones of the entire cluster in both the HOH bending and the OH stretching regions of the spectra.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Monte Carlo calculation and verification of the geometrical factors for the NPDGamma experiment

The NPDGamma experiment measures the parity-violating asymmetry in gamma-ray emission in the capture of polarized neutrons on liquid parahydrogen. The sensitivity to the asymmetry for each detector in the array is used as a parameter in the extraction of the physics asymmetry from the measured data. The detector array is approximately cylindrically symmetric around the target and a step-wise sinusoidal function has been used for the sensitivity in the previous iteration of the NPDGamma experiment, but deviations from cylindrical symmetry necessitate the use of a Monte Carlo model to determine corrections to the geometrical factors. For the calculations, source code modifications to MCNPX were done in order to calculate the sensitivity of each cesium iodide detector to the physics asymmetry. We describe the MCNPX model and results from calculations and how the results are validated through measurement of the parity violating asymmetry of gamma-rays from neutron capture on chlorine.

Grammer, K. B.↗
DOE OSTI2015

Adsorption and diffusion of Ru adatoms on Ru(0001)-supported graphene: Large-scale first-principles calculations

Large-scale first-principles density functional theory calculations are performed to investigate the adsorption and diffusion of Ru adatoms on monolayer graphene (G) supported on Ru(0001). The G sheet exhibits a periodic moiré-cell superstructure due to lattice mismatch. Within a moiré cell, there are three distinct regions: fcc, hcp, and mound, in which the C{sub 6}-ring center is above a fcc site, a hcp site, and a surface Ru atom of Ru(0001), respectively. The adsorption energy of a Ru adatom is evaluated at specific sites in these distinct regions. We find the strongest binding at an adsorption site above a C atom in the fcc region, next strongest in the hcp region, then the fcc-hcp boundary (ridge) between these regions, and the weakest binding in the mound region. Behavior is similar to that observed from small-unit-cell calculations of Habenicht et al. [Top. Catal. 57, 69 (2014)], which differ from previous large-scale calculations. We determine the minimum-energy path for local diffusion near the center of the fcc region and obtain a local diffusion barrier of ∼0.48 eV. We also estimate a significantly lower local diffusion barrier in the ridge region. These barriers and information on the adsorption energy variation facilitate development of a realistic model for the global potential energy surface for Ru adatoms. This in turn enables simulation studies elucidating diffusion-mediated directed-assembly of Ru nanoclusters during deposition of Ru on G/Ru(0001)

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

Thermoelectric properties of rocksalt ZnO from first-principles calculations

Zinc oxide (ZnO) undergoes a pressure-induced structural transition from its normal ambient-pressure wurtzite (WZ) phase to a rocksalt (RS) phase around 10 GPa. A recent experiment shows that the high-pressure RS ZnO phase can be recovered and stabilized at ambient conditions, which raises exciting prospects of expanding the range of properties of ZnO. For a fundamental understanding of the RS ZnO phase, we have performed first-principles calculations to determine its electronic, phonon, and thermodynamic properties at high (20 GPa) and ambient (0 GPa) pressure. Furthermore, we have calculated its electrical and thermal transport properties, which allow an evaluation of its thermoelectric figure of merit ZT at different temperature and doping levels. Our calculations show that the ambient-pressure RS ZnO phase can reach ZT values of 0.25 to 0.3 under both n-type and p-type doping in a large temperature range of 400 K to 800 K, which is considerably lower than the temperature range of 1400 K to 1600 K where WZ ZnO reaches similar ZT values. These results establish RS ZnO as a promising material for thermoelectric devices designed to operate at temperatures desirable for many heat recovery applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Coil current and vacuum magnetic flux calculation for axisymmetric equilibria

In fixed-boundary axisymmetric equilibrium calculations the plasma shape is assigned from input. In several circumstances, the plasma shape may not be known a priori, or one may desire to also compute the magnetic field in the volume surrounding the plasma through the calculation of a free-boundary equilibrium. This requires either the coil currents or the magnetic poloidal flux on a curve in the vacuum region to be assigned as input for the free-boundary equilibrium calculation. The FREE-FIX code presented in this article is a general tool for calculating coil currents being given a fixed-boundary calculation. Finally, a new formulation is presented, which considerably reduces the computational cost of the calculation. FREE-FIX performs well for different geometries and experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Evaluation of The NSUF Reactor Activation and Damage (RAD) Calculator: Assessing the Accuracy of the Activation Calculator

Nuclear Science User Facilities developed an activation calculator as part of the Combined Material Experiment Toolbox project. The activation calculator estimates the radionuclide concentration in an irradiated sample to indicate when and where the samples may best be examined to only help with scoping the experiment. To verify the accuracy of the calculator, the results from the specific gamma dose rate [mrem/hr/g] at various cooldown times were compared to equivalent results generated in ORIGEN. Most elements were tested as a sample composed solely of that element. These samples were exposed to thirteen different lengths of irradiation in all reactor positions included in the calculator. For elements of interest (i.e., ones currently included in the Nuclear Fuels and Materials Library), 83% of the test cases fell within the allowed tolerance for this scoping tool of ±50%. The overall error is logarithmic with irradiation length and is almost completely independent of reactor position. The primary source of error is missing nuclear data and reactions. NSUF believes that the RAD Calculator can be deployed to the user community because it is designed to provide estimates to aid researchers in planning experiments, and users are cautioned about the limitations of the calculator.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2017

Application of the new IAPWS Guideline on the fast and accurate calculation of steam and water properties with the Spline-Based Table Look-Up Method (SBTL) in RELAP-7

The numerical simulation of thermalhydraulic processes in nuclear power plants requires very accurate and extremely fast algorithms for calculating the thermophysical properties of water and steam. In order to provide such algorithms, the International Association for the Properties of Water and Steam (IAPWS) has adopted the new “IAPWS Guideline on the Fast Calculation of Steam and Water Properties with the Spline-Based Table Look-Up Method (SBTL)”. In this article, the SBTL method is applied to property functions of specific volume and specific internal energy (v,e) based on the scientific formulation IAPWS-95 and the latest IAPWS formulations for transport properties. From the newly generated SBTL functions, thermodynamic and transport properties as well as their derivatives and inverse functions are calculable in the fluid range of state for pressures up to 100 MPa and for temperatures up to 1273 K, including the metastable liquid and the metastable vapor regions. The SBTL functions reproduce the underlying formulations with an accuracy of 10–100 ppm and significantly reduced computing times. The SBTL method has been implemented into the nuclear reactor system safety analysis code RELAP-7 [2] to consider the real fluid behavior of water and steam in a novel 7-equation two-phase flow model.

42 ENGINEERING↗
DOE OSTI2021

Verification and Validation Tests of Gamma Library of MC2-3 for Coupled Neutron and Gamma Heating Calculation

For the accurate assessment of the heat generation rate in fast reactors, the gamma library of MC2 -3 and the MC2 -3 + GAMSOR procedure employing the coupled neutron and gamma heating calculation has been thoroughly verified against Monte Carlo results and validated using the ZPPR-15D gamma dose measurement data. NJOY outputs are post-processed in a consistent way with the NJOY procedure to avoid any missing data or double counting of data. Prompt heating for 379 out of 391 isotopes in the gamma library was verified against MCNP6.2 to within 1% relative error in total heating for most isotopes. For both a simple one-dimensional slab problem representing a sodium cooled fast reactor and the Experimental Breeder Reactor II (EBR-II) benchmark problem, the root-mean-square values of assembly power error were less than 0.5% for fuel assemblies, ~1% in blankets and ~1 to ~3% in reflectors compared to MCNP6.2 results. The most plausible cause for the 3% error in a reflector assembly is believed to be the error in the multigroup neutron cross sections for the reflector assembly. For validation, gamma doses measured with thermoluminiscent dosimeters (TLDs) in the ZPPR-15D experiment were calculated using GAMSOR. Due to the uncertainty in the TLD measurement with regards to the energy deposition of photons and neutrons, the validation data leads to a 12.7% uncertainty on the experimental measurement. With this uncertainty bound, the calculated doses all fell within one standard deviation of the measured value. Combined with the accurate calculation of reaction rate distributions and neutron spectrum measurements, these results indicate good agreement for the neutron and gamma heating calculations that were performed.

coupled neutron and gamma heating↗
DOE OSTI2014

CPMC-Lab: A Matlab package for Constrained Path Monte Carlo calculations

We describe CPMC-Lab, a Matlab program for the constrained-path and phaseless auxiliary-field Monte Carlo methods. These methods have allowed applications ranging from the study of strongly correlated models, such as the Hubbard model, to ab initio calculations in molecules and solids. The present package implements the full ground-state constrained-path Monte Carlo (CPMC) method in Matlab with a graphical interface, using the Hubbard model as an example. The package can perform calculations in finite supercells in any dimensions, under periodic or twist boundary conditions. Importance sampling and all other algorithmic details of a total energy calculation are included and illustrated. This open-source tool allows users to experiment with various model and run parameters and visualize the results. It provides a direct and interactive environment to learn the method and study the code with minimal overhead for setup. Furthermore, the package can be easily generalized for auxiliary-field quantum Monte Carlo (AFQMC) calculations in many other models for correlated electron systems, and can serve as a template for developing a production code for AFQMC total energy calculations in real materials. Several illustrative studies are carried out in one- and two-dimensional lattices on total energy, kinetic energy, potential energy, and charge- and spin-gaps.

Computer Science↗
DOE OSTI2025

GPR_calculator: An on-the-fly surrogate model to accelerate massive nudged elastic band calculations

We present GPR_calculator, a package based on Python and C++ programming languages to build an on-the-fly surrogate model using Gaussian Process Regression (GPR) to approximate computationally expensive electronic structure calculations. The key idea is to dynamically train a GPR model during the simulation that can accurately predict energies and forces with uncertainty quantification. When the uncertainty is high, the costly electronic structure calculation is performed to obtain the ground truth data, which is then used to update the GPR model. To illustrate the effectiveness of GPR_calculator, we demonstrate its application in Nudged Elastic Band (NEB) simulations of surface diffusion and reactions, achieving 3-10 times acceleration compared to pure ab initio calculations. The source code is available at https://github.com/MaterSim/GPR_calculator.

Gaussian process regression↗
DOE OSTI2020

Potassium iodide cluster based superhalogens and superalkalis: Theoretical calculations and experimental confirmation

We studied a series of potassium-iodide clusters with formulas of (KI)nI- and (KI)nK+ (n = 1-3) by quantum chemical calculations. The calculated vertical detachment energies (VDEs) of (KI)nI- clusters are all higher than those of halogen anions, and thus can be classified as superhalogen anions; while the calculated vertical electron affinities (VEAs) of (KI)nK+ clusters are all lower than those of alkali metal cations, and thus they can be recognized as superalkali cations. Our calculated VEA for (KI)K+, i.e., K2I+, agrees well with the previously measured value by thermal ionization mass spectrometry. We confirmed the theoretical predictions by measuring the negative ion photoelectron (NIPE) spectra of (KI)nI- anions, from which the measured VDEs are found to be in good agreement with the theoretically calculated values. In addition, we also explored the possibility of supersalts formed by these superhalogen anions and superalkali cations, for example, K3I3 can be regarded as a supersalt formed by K2I+ and KI2– subunits.

potassium-iodide clusters, superhalogen, superalka↗
DOE OSTI2018

Kinetics calculation of fast periodic pulsed reactors using MCNP6

Fast periodic pulsed reactor is a type of reactor in which the fission bursts are formed entirely with external reactivity modulation with a specified time periodicity. This type of reactors could generate much larger intensity of neutron beams for experimental use, compared with the steady state reactors. In the design of fast periodic pulsed reactors, the time dependent simulation of the power pulse is majorly based on a point kinetic model, which is known to have limitations. A more accurate calculation method is desired for the design analyses of fast periodic pulsed reactors. Monte Carlo computer code MCNP6 is used for this task due to its three dimensional transport capability with a continuous energy library. Some new routines were added to simulate the rotation of the movable reflector parts in the time dependent calculation. Fast periodic pulsed reactor IBR-2M was utilized to validate the new routines. This reactor is periodically in prompt supercritical state, which lasts for similar to 400 mu s, during the equilibrium state. This generates long neutron fission chains, which requires tremendously large amount of computation time during Monte Carlo simulations. Russian Roulette was applied for these very long neutron chains in MCNP6 calculation, combined with other approaches to improve the efficiency of the simulations. In the power pulse of the IBR-2M at equilibrium state, there is some discrepancy between the experimental measurements and the calculated results using the point kinetics model. MCNP6 results matches better the experimental measurements, which shows the merit of using MCNP6 calculation relative to the point kinetics model.

Zhong, Z.↗