Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulation package”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

FINETUNA: fine-tuning accelerated molecular simulations

Abstract Progress towards the energy breakthroughs needed to combat climate change can be significantly accelerated through the efficient simulation of atomistic systems. However, simulation techniques based on first principles, such as density functional theory (DFT), are limited in their practical use due to their high computational expense. Machine learning approaches have the potential to approximate DFT in a computationally efficient manner, which could dramatically increase the impact of computational simulations on real-world problems. However, they are limited by their accuracy and the cost of generating labeled data. Here, we present an online active learning framework for accelerating the simulation of atomic systems efficiently and accurately by incorporating prior physical information learned by large-scale pre-trained graph neural network models from the Open Catalyst Project. Accelerating these simulations enables useful data to be generated more cheaply, allowing better models to be trained and more atomistic systems to be screened. We also present a method of comparing local optimization techniques on the basis of both their speed and accuracy. Experiments on 30 benchmark adsorbate-catalyst systems show that our method of transfer learning to incorporate prior information from pre-trained models accelerates simulations by reducing the number of DFT calculations by 91%, while meeting an accuracy threshold of 0.02 eV 93% of the time. Finally, we demonstrate a technique for leveraging the interactive functionality built in to Vienna ab initio Simulation Package (VASP) to efficiently compute single point calculations within our online active learning framework without the significant startup costs. This allows VASP to work in tandem with our framework while requiring 75% fewer self-consistent cycles than conventional single point calculations. The online active learning implementation, and examples using the VASP interactive code, are available in the open source FINETUNA package on Github.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Explanation of the seasonal variation of cosmic multiple muon events observed with the NOvA Near Detector

The flux of cosmic ray muons at the Earth’s surface exhibits seasonal variations due to changes in the temperature of the atmosphere affecting the production and decay of mesons in the upper atmosphere. Using 1473 live days of data collected by the NuMI Off-axis 𝜈 𝑒 Appearance (NOvA) Near Detector during 2018–2022, we studied the seasonal pattern in the multiple-muon event rate. The data confirm an anticorrelation between the multiple-muon event rate and effective atmospheric temperature, consistent across all the years of data. Previous analyses from MINOS and NOvA saw a similar anticorrelation but did not include an explanation. We find that this anticorrelation is driven by altitude–geometry effects as the average muon production height changes with the season. This has been studied with a CORSIKA cosmic ray simulation package by varying atmospheric parameters, and provides an explanation to a longstanding discrepancy between the seasonal phases of single and multiple-muon events.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Neural conditional reweighting

There is a growing use of neural network classifiers as unbinned, high-dimensional (and variable-dimensional) reweighting functions. To date, the focus has been on marginal reweighting, where a subset of features are used for reweighting while all other features are integrated over. There are some situations, though, where it is preferable to condition on auxiliary features instead of marginalizing over them. Here, we introduce neural conditional reweighting, which extends neural marginal reweighting to the conditional case. This approach is particularly relevant in high-energy physics experiments for reweighting detector effects conditioned on particle-level truth information. Furthermore we leverage a custom loss function that not only allows us to achieve neural conditional reweighting through a single training procedure, but also yields sensible interpolation even in the presence of phase space holes. As a specific example, we apply neural conditional reweighting to the energy response of high-energy jets, which could be used to improve the modeling of physics objects in parametrized fast simulation packages.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

GIWAXS-SIIRkit: scattering intensity, indexing and refraction calculation toolkit for grazing-incidence wide-angle X-ray scattering of organic materials

Grazing-incidence wide-angle X-ray scattering (GIWAXS) has become an increasingly popular technique for quantitative structural characterization and comparison of thin films. For this purpose, accurate intensity normalization and peak position determination are crucial. At present, few tools exist to estimate the uncertainties of these measurements. Here, a simulation package is introduced called GIWAXS-SIIRkit , where SIIR stands for scattering intensity, indexing and refraction. Furthermore, the package contains several tools that are freely available for download and can be executed in MATLAB. The package includes three functionalities: estimation of the relative scattering intensity and the corresponding uncertainty based on experimental setup and sample dimensions; extraction and indexing of peak positions to approximate the crystal structure of organic materials starting from calibrated GIWAXS patterns; and analysis of the effects of refraction on peak positions. Additionally, each tool is based on a graphical user interface and designed to have a short learning curve. A user guide is provided with detailed usage instruction, tips for adding functionality and customization, and exemplary files.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Large-Scale Inverter Integration in Bulk Power Grids Using Heterogeneous Grid-Forming Control Strategies

The conventional generation within the grid is gradually being replaced by inverter based resources. However, the integration of inverters into realistic transmission grid models has not been thoroughly explored. To address this gap, this work develops an automation framework aimed at facilitating the integration of a large number of inverter-based resources into a large-scale grid. Two types of grid-forming controlled inverters—namely droop-controlled inverters and virtual synchronous machine-controlled inverters—are developed for use in positive-sequence simulation packages. Various penetration levels ranging from 23% to 100% of grid-forming inverters are tested in a grid setup in the US Western Interconnection. The dynamic performance of inverters under different control schemes is compared in response to various disturbances.

Lyu, Xue [BATTELLE (PACIFIC NW LAB)]↗
DOE OSTI2026

GPU-Accelerated Analytic Simulation of Sparse Ionization Signal Formation in Pixelated Projection Detector

This paper presents a GPU-accelerated simulation package, TRED, for next-generation neutrino detectors with pixelated charge readout, leveraging community-driven software ecosystems to ensure adaptability and extensibility. We introduce two generic contributions: (i) an effective-charge representation based on Gaussian quadrature rules, in which the linear- interpolation factors for the field response inside each voxel are absorbed into the effective charge, and (ii) a sparse, block- binned tensor representation that enables efficient FFT-based computation of induced signals on readout electrodes for sparsely activated detector volumes. The former captures structure inside a voxel without dense sampling, while the latter achieves low memory usage and scalable runtime, as demonstrated in bench- mark studies. The underlying data representation is applicable to large-scale detectors and to other computational problems involving sparse activity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Multi-Agent Simulation Based Framework for Power Restoration Time Estimation at Distribution Level

The growing frequency of power outages has prompted increased interest in developing a more resilient power grid that can quickly recover from weather-related damage. At the distribution level, power restoration is a complex, multi-stage process involving multiple response entities. Providing utility stakeholders, government regulators, and the public with information about outage duration and estimated time to restoration is crucial. The research employs a multi-agent simulation approach, which allows for the simulation of decision-making behaviors among different entities and the incorporation of various uncertainties. Specifically, the study uses the open-source simulation package Mesa-Geo in conjunction with the Python language and constructs a road network using the open-source network extension pgRouting for routing queries. The research design includes several experiments focused on Florida as a case study, comparing repair crew sizes, power outage numbers, and road damage scenarios. The findings could offer valuable managerial guidance on resource allocation in the restoration process.

Chen, Yang↗
DOE OSTI2022

Low energy neutrino detection with a compact water-based liquid scintillator detector

In this study, the conceptual design and physics simulations of a near-field Water-based Liquid Scintillator (WbLS) detector placed 100 m from the Akkuyu Nuclear Power Plant (ANPP), currently under construction and aiming at being Turkey’s first nuclear power plant, is presented. The ANPP is an excellent opportunity for neutrino studies and the development of an R &D program for neutrino detectors in Turkey. The Reactor Neutrino Experiments of Turkey (RNET) program includes a compact detector with a 2.5-ton volume of WbLS and a ~ 30% photo-coverage, and the program is planned to be expanded with a medium-size 30-ton detector that will be an international testbed for WbLS and new detector technologies through low energy neutrino studies. In the following, the focus will be on the smaller ~2.5 ton detector, instrumented with 8-in. high quantum efficiency PMTs and two layers of cosmic veto paddles, covering all sides of the detector, to track and veto cosmic particles. Inverse Beta Decay (IBD) events from electron anti-neutrinos generated in the reactor core are simulated using the RAT-PAC simulation package and several liquids with different percentages of Liquid Scintillator (LS) and Gadolinium (Gd) are investigated.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

RadSim: Detector Response

This package simulates the expected output gamma spectra from a gamma detector for given configuration.

Hangal, DhanushA [Lawrence Livermore National Labo↗
DOE OSTI2023

NbTaV_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution ternary alloys niobium-tantalum-vanadium (Nb-Ta-V). The first-principle code that has been used to run the calculations is the closed-source Vienna Ab-Initio Simulation Package (VASP). The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 8. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 10,500 randomized atomic structures over 105 chemical compositions. The calculations have been collected on NERSC-Perlmutter using the VASP 6.3.2. The VASP calculations for every atomic structure have been performed in 2 main steps: 1. Starting from an ideal body-centered-cubic (BCC) structure, geometry optimization with low precision has been executed to perform a preliminary optimization of the atomic structure. The output for this calculations is available in the files 0.CONTCAR, 0.OUTCAR, rlx1.out. 2. Using the atomic structure resulting from the preliminary geometry optimization, a second geometry optimization has been performed using normal precision. The output for this calculations is available in the files CONTCAR, OUTCAR, rlx2.out, vaspout.h5, and vasprun.xml. Every chemical composition sampled across the composition range in the dataset has its own directory. The convention used to name the directories for ternary alloys is AXBYCZ, where A, B, and C refer to the constituents, and X, Y, and Z are positive integers that represent the number of atoms for each constituent and their values still sum up to 128. Each atomic structure associated with a specific chemical composition has its own sub-directory within the directory of the corresponding chemical composition. The sub-directories for each atomic structure for each chemical composition are named 'case-*', where * is a positive integer that spans all the values from 1 through 100, extremes included. The files contained in each sub-directory 'case-*' for each atomic structure are as follows: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. 0.POSCAR: input file that defines the atomic structure of a system 3. 0.CONTCAR: output file that provides the atomic positions and cell parameters after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 4. 0.OUTCAR: output file that contains detailed information about the progress of a calculation after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 5. rlx1.out: file with diagnostic information about the execution of the first geometry optimization with precision variable set to PREC=Low in the INCAR file 6. POSCAR: input file that defines the atomic structure of a system after the first geometry optimization has been run at low precision. This represents the input for the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 7. CONTCAR: output file that provides the atomic positions and cell parameters after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 8. OUTCARL: output file that contains detailed information about the progress of a calculation after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 9. rlx2.out: file with diagnostic information about the execution of the second geometry optimization with precision variable set to PREC=Normal in the INCAR file 10. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 11. vasprun.xml: contains similar information to OUTCAR, but in an xml format. This research is sponsored by the Artificial Intelligence Initiative as part of the Laboratory Directed Research and Development (LDRD) Program of Oak Ridge National Laboratory, managed by UT-Battelle, LLC, for the US Department of Energy under contract DE-AC05-00OR22725. This work used resources of the Oak Ridge Leadership Computing Facility, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC05-00OR22725, under Directorate Discretionary awards MAT025 (Materials Science) and LRN026 (Machine Learning), and INCITE award MAT201. This work also used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231, under award ERCAP0025216.

36 MATERIALS SCIENCE↗
DOE OSTI2024

VZr_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys vanadium-zirconium (V-Zr). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach \cite{PAW}. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements V and Zr. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

TaTi_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys tantalum-titanium (Ta-Ti). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Hf and Nb. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on Air Force HPC11 cluster using the VASP 6.5.1. Additional methodology and file structure information is available in the dataset README.txt file.

36 MATERIALS SCIENCE↗
DOE OSTI2025

HfNb_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys hafnium-niobium (Hf-Nb). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Hf and Nb. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect BCC lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,029 randomized atomic structures over 31 chemical compositions. Further methodological and structural information is contained in the dataset README.txt file.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Thermoelasticity Model for Lithium Fluoride

Lithium fluoride (LiF) is a typical ionic solid, and at the same time a typical polymorphic (multi-phase) material. Its phase diagram contains two solid phases, B1 (rock salt) and B2 (CsCl), which are typical of ionic solids. The phase diagram of LiF was suggested by Smirnov. We have carried out extensive study of the properties of LiF using one of the most recent versions of VASP (Vienna Ab initio Simulation Package) which allows one to carry out calculations on ionic solids; more detail will be presented elsewhere.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Advanced CO 2 Capture Solvent Systems for Dynamic Power Generation

RTI International, in collaboration with Pacific Northwest National Laboratory (PNNL), Carbon Capture Simulation for Industry Impact (CCSI 2 ), Electricity Power Research Institute (EPRI), and West Virginia University (WVU), successfully completed a joint research effort in developing a cost-effective, resilient, load-following advanced CO 2 capture technology for natural gas power plants. The project’s objective was to develop a CO 2 capture process that maximizes the net present value (NPV) of the electricity sale by minimizing the levelized cost of electricity (LCOE) under dynamic plant loads and high renewable penetration environments. The two key innovations developed in this project were the use of (i) advanced water-lean solvents (WLSs) and (ii) process intensification equipment, such as a rotating packed bed (RPB) absorber and dual-stage flash regeneration. The process’s low CO 2 capture cost is realized through WLSs’ low energy required for solvent regeneration, which lowers the operating cost while RPBs intensify the absorption process and reduces the power plant capital cost. A suite of advanced computational and simulation packages was implemented to guide the process design, validate the dynamic response of the capture plant, evaluate system-wide performance, and maximize the power plant’s profit. The project also engaged with power producers and other stakeholders to ensure its technical relevance and techno-economic viability. The development of this highly disruptive CO 2 capture technology could accelerate the industry adoption and thereby lower the greenhouse gas emissions of the U.S. power sector. Deployment of this technology can increase the reliability and decrease the cost of electricity generation in the U.S. by enabling the use of low-carbon fossil fuels to balance fluctuations of renewable energy availability.

03 NATURAL GAS↗
DOE OSTI2024

Imaging Enhanced Simulation Validation for Directed Energy Deposition Additive Manufacturing (Final CRADA Report – NFE-22-09107)

Directed Energy Deposition (DED) is a welding-based metal Additive Manufacturing (AM) process that relies on the programmed rastering of an electric arc or laser induced weld pool to construct a component in a layerwise fashion. The induced complex thermal field and uneven thermal contraction depends on the printed geometry and scan pattern, and as such, accumulated residual stresses and distortions are complex and difficult to predict. Several prior works have resulted in tools to combat this issue; ANSYS has developed a thermoplastic simulation package targeting DED AM, and ORNL has developed an in-situ imaging sensor package ‘Stereo Correlated Optical and Pyrometric System’ (SCOPS) that can spatially monitor temperature and full field strain. In this CRADA, these tools are compared in order to validate the results and complimentarily address the weaknesses in each other.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Research and Development Studies for Reactor Neutrino Experiments in Turkey (RNET)

The program of the Reactor Neutrino Experiments of Turkey includes a small portable Water-based liquid scintillator detector to detect neutrinos from the Akkuyu nuclear power plant, planned to begin operating in 2023. The small near-field detector will weigh about 2-3 tons and will be placed less than 100 meters from the reactor cores. The Reactor Neutrino Experiments of Turkey program also includes a medium-size 30-ton Water-based liquid scintillator detector, which will be placed 1-2 km away from the reactor cores and will be used as a far detector. Both detectors and their responses to neutrino interactions were simulated using a GEANT4-based RAT-PAC simulation package. Here, we will share the technical and physical details of both detectors, and discuss the ongoing R&D effort for neutrino studies in Turkey.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

A Physical Model for Understanding the Activation of MoS 2 Basal-plane Sulfur Atoms for the Hydrogen Evolution Reaction

All the DFT calculations are done with the Vienna Ab Initio Simulation Package (VASP) using the projector augmented wave method. The Bayesian error estimation exchange-correlation functionals (BEEF) with van der Waals interactions are employed. A plane-wave cutoff energy of 400 eV is used together with PAW-PBE potentials where semi core p states are treated as valence. All the calculations allow for spin-polarization. The structures are relaxed until the force is converged to < 0.01 eV/Å. The lattice parameter of MoS 2 unit cell, optimized with this functional, is 3.19 Å. A (4×4) supercell is used to model all the transition metal doped MoS 2 systems studied here, including those with S vacancies. For calculations in the initial dopant structure exploration, the Brillouin zone is sampled with a 3x3x1 Monkhorst-Pack k-point mesh. A 6×6×1 Monkhorst-Pack k-point mesh is used for the H binding energy and density of states calculations. In all calculations, the vacuum layer is set as 15 Å to eliminate periodic interaction perpendicular to the basal plane.

77 NANOSCIENCE AND NANOTECHNOLOGY↗