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

Summary of Structural Material Modeling Development for the NEAMS Program in Fiscal Year 2020

This report summarizes work performed during Fiscal Year (FY) 2020 at Idaho National Laboratory (INL) for the U.S. Department of Energy?s Nuclear Engineering Advanced Modeling and Simulation (NEAMS) program for the Structural Materials and Chemistry Technical Area in the work package entitled "MS- 20IN050104 - Structural Materials - INL." The Structural Materials and Chemistry Technical Area is a relatively new component of the NEAMS program, and is currently focusing on developing simulation capabilities to support the deployment of nuclear energy in the areas of molten salt reactor chemistry, light water reactor (LWR) structural material degradation, and structural material behavior for advanced reactor applications. INL performed work for to advance capabilities for simulation of structural material behavior in both LWR and advanced reactor applications in the work described here.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

HELIUM LEAK TEST MODELING OF A SPENT NUCLEAR FUEL CANISTER

The U.S. Department of Energy (DOE) is considering the development of one or more federal consolidated interim storage facilities (CISFs) to be used to store commercial spent nuclear fuel (SNF) at locations in the U.S. One of the first technical challenges of a CISF is performing an inspection of SNF canisters upon their receipt to confirm they can be placed into the CISF’s licensed storage configuration. The canister receipt inspection is critical to CISF site operations. The test is conceived as being a helium (He) leak check, intended to confirm that the confinement boundary of a SNF canister is intact. SNF canisters are filled with He when they are sealed, so detection of a He leak indicates that a through-wall flaw has occurred in the canister confinement boundary. Other measurements are planned to occur upon canister receipt in addition to the He leak check such as krypton-85 measurements, which would indicate confinement breaches of one or more fuel rods in addition to a breach of the SNF canister. However, the He leak check has been identified as one of such high importance and has such significant technical challenges that a full-scale demonstration is needed to confirm the He leak test’s viability and to assist in planning relative to its operational requirements. A modeling methodology for simulating the He detection test was developed to help inform the test plan and the design of the test vessels. To develop the modeling methodology a detailed computational fluid dynamics (CFD) benchmark model was constructed to compare against leak rate test data from a transportation package for radioactive material. This report is focused on modeling efforts to simulate the benchmark leak test.

Suffield, Sarah R.↗
DOE OSTI2020

Electron transport in gaseous detectors with a Python-based Monte Carlo simulation code

Understanding electron drift and diffusion in gases and gas mixtures is a topic of central importance for the development of modern particle detection instrumentation. The industry-standard MagBoltz code has become an invaluable tool during its 20 years of development, providing capability to solve for electron transport (‘swarm’) properties based on a growing encyclopedia of built-in collision cross sections. We have made a refactorization of this code from FORTRAN into Cython, and studied a range of gas mixtures of interest in high energy and nuclear physics. The results from the new open source PyBoltz package match the outputs from the original MagBoltz code, with comparable simulation speed. An extension to the capabilities of the original code is demonstrated, in implementation of a new Modified Effective Range Theory interface. We hope that the versatility afforded by the new Python code-base will encourage continued use and development of the MagBoltz tools by the particle physics community.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

BEEAM (Building Electrical Efficiency Analysis Model) [SWR-20-107]

Modern high-performance buildings exhibit an increasing number of building loads that use direct current (DC) electricity internally, rather than alternating current (AC), due to the advent of low-cost computing and advanced power electronics. Powering DC devices from the AC grid requires AC/DC power conversion, which introduces energy losses and reduces efficiency. As DC loads proliferate, the cumulative wasted energy associated with hundreds of millions of AC/DC converters has become one of the broadest energy savings opportunities in buildings. DC power distribution systems have been proposed as an elegant and transformative solution to the problem of DC devices. In a DC distribution system, the building's wires carry DC electricity, rather than AC; and a few centralized, highly efficient AC/DC converters replace the many smaller, less efficient converters that serve individual DC loads. Unfortunately, the trade-offs associated with DC distribution systems are not well understood. Reported energy savings associated with DC power distribution differ widely and have not been well validated. The Building Electrical Efficiency Analysis Model (BEEAM) is a Modelica library that simulates the efficiency of building electrical distribution systems using harmonic power flow. BEEAM can model a wide variety of building distribution topologies, including three-phase AC, single-phase AC, unipolar DC, bipolar DC, and hybrid networks under both balanced and unbalanced load conditions. BEEAM accurately models power electronic converter losses, provides granular estimates of losses throughout the distribution system, and properly models efficiency at part load conditions. Users can package BEEAM within a functional mockup unit (FMU), enabling co-simulation with other modeling platform, such as EnergyPlus. In summary, BEEAM provides a tool suite for fair and accurate comparison of the efficiency of building electrical distribution systems, including AC, DC, and hybrid systems.

Frank, Stephen↗
DOE OSTI2020

Composite Analysis Solid Waste Release Data Reduction of Activity Flux from Waste Sites to the Vadose Zone for Baseline Assessment

This environmental calculation file (ECF) describes the data reduction methodology for the solid waste release datasets generated by the Composite Analysis Solid Waste Release (CASWR) model (CP-62766, Model Package Report: Composite Analysis Solid Waste Release Model [CASWR Model]) and documents the resulting reduced datasets for the updated Composite Analysis (CA) baseline assessment. This ECF is in support of the updated Hanford Site CA. The purpose of the updated Hanford Site CA is to provide an estimate of the cumulative radiological impacts from active and planned low-level radioactive waste disposals and other potentially interacting radioactive waste disposal sources that will remain following Hanford Site closure. The vadose zone modeling activities associated with the CA includes solid waste forms as contaminant sources. To support the modeling of solid waste transport from the vadose zone to the groundwater at the Hanford Site, the CASWR model was designed to generate deterministic radionuclide release rates for the Hanford Site Central Plateau solid waste disposal sites. The complete CASWR model-generated dataset for the CA Baseline Assessment consists of 2,378,831 sets of time and release rates (1,073,869 sets for the 200 East Area and 1,304,962 sets for the 200 West Area), representing a course of the simulation from year 1953 to year 12,069, 16 radiological contaminants of potential concern (COPCs), and 237 solid waste sources as documented in ECFHANFORD- 19-0112, Solid Waste Release Calculations for the Composite Analysis Baseline Assessment. For several of the STOMP vadose zone models, the number of release rate data pairs required for release at all waste sites would exceed the limitations imposed by the STOMP modelling software package used to simulate contaminant transport through the vadose zone. As a result, a data reduction of the CASWR dataset is required prior to its use as a source input for the CA vadose zone modeling activities.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Simulating Superconducting Devices in Novel Materials Using the G4CMP Simulation Toolkit

Understanding the kinematics of phonons and charge propagation in superconducting devices is crucial for conducting low-threshold dark matter searches and minimizing correlated errors in superconducting qubits. For nearly a decade, the Geant4 Condensed Matter Physics (G4CMP) package, originally developed for the Cryogenic Dark Matter Search (CDMS) experiment, has been limited to simulating charge and phonon transport in silicon and germanium materials. In this work, we have expanded the capabilities of G4CMP to include novel substrate materials such as sapphire (Al2O3), Gallium Arsenide (GaAs), Lithium Fluoride (LiF), Calcium Tungstate (CaWO4), and Calcium Fluoride (CaF2). We demonstrate the use of this toolkit in generating phonon transport properties for these materials and compare the results with experimentally determined values where available. Additionally, we explore the performance of superconducting devices utilizing these materials.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

A coupled fluid-mechanical workflow to simulate the directed energy deposition additive manufacturing process

Simulation of additive manufacturing processes can provide essential insight into material behavior, residual stress, and ultimately, the performance of additively manufactured parts. In this work, we describe a new simulation based workflow utilizing both solid mechanics and fluid mechanics based formulations within the finite element software package SIERRA (Sierra Solid Mechanics Team in Sierra/Solid Mechanics 4.52 User’s Guide SAND2019-2715. Technical report, Sandia National Laboratories, 2011) to enable integrated simulations of directed energy deposition (DED) additive manufacturing processes. In this methodology, a high-fidelity fluid mechanics based model of additive manufacturing is employed as the first step in a simulation workflow. This fluid model uses a level set field to track the location of the boundary between the solid material and background gas and precisely predicts temperatures and material deposition shapes from additive manufacturing process parameters. Next, the resulting deposition shape and temperature field from the fluid model are then mapped into a solid mechanics formulation to provide a more accurate surface topology for radiation and convection boundary conditions and a prescribed temperature field. Solid mechanics simulations are then conducted to predict the evolution of material stresses and microstructure within a part. By combining thermal history and deposition shape from fluid mechanics with residual stress and material property evolutions from solid mechanics, additional fidelity and precision are incorporated into additive manufacturing process simulations providing new insight into complex DED builds.

36 MATERIALS SCIENCE↗
DOE OSTI2021

AENET–LAMMPS and AENET–TINKER : Interfaces for accurate and efficient molecular dynamics simulations with machine learning potentials

Machine-learning potentials (MLPs) trained on data from quantum-mechanics based first-principles methods can approach the accuracy of the reference method at a fraction of the computational cost. To facilitate efficient MLP-based molecular dynamics and Monte Carlo simulations, an integration of the MLPs with sampling software is needed. Here, we develop two interfaces that link the atomic energy network (ænet) MLP package with the popular sampling packages TINKER and LAMMPS. The three packages, ænet, TINKER, and LAMMPS, are free and open-source software that enable, in combination, accurate simulations of large and complex systems with low computational cost that scales linearly with the number of atoms. Scaling tests show that the parallel efficiency of the ænet–TINKER interface is nearly optimal but is limited to shared-memory systems. The ænet–LAMMPS interface achieves excellent parallel efficiency on highly parallel distributed memory systems and benefits from the highly optimized neighbor list implemented in LAMMPS. We demonstrate the utility of the two MLP interfaces for two relevant example applications: the investigation of diffusion phenomena in liquid water and the equilibration of nanostructured amorphous battery materials.

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

Benchmark Calculation for Turkey Point Unit 3 Cycles 1-3 Using the SCALE 6.3/Polaris–PARCS v3.4.2 Code Package

Benchmark calculations were performed for Turkey Point Unit 3 cycles 1–3 to validate the SCALE 6.3/Polaris–PARCS v3.4.2 code with the ENDF/B–VII.1 56–group library by comparing the simulated results with the measured data. The benchmark results will be used in evaluating the SCALE/Polaris–PARCS code package’s uncertainties for pressurized water reactor physics analysis. That future analysis will include key nuclear parameters such as reactivity, control bank worth, temperature coefficients, and pin and assembly power peaking factors. The present document details plant and fuel design specifications and input data for SCALE/Polaris, GenPMAXS, and PARCS. Additional details are provided with respect to the input and output files produced for the benchmark calculations. The benchmark results are summarized such that they can be used in evaluating uncertainties with other benchmark results for key nuclear parameters.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Role of E × B Drift in Divertor Detachment Control via Boron Powder Injection on EAST

Here, the effects of B powder injection on plasma detachment about EAST discharge were studied by using SOLPS-ITER code package with the effects of E × B drifts considered. The simulation results show that plasma detachment occurs at the inner target in favourable toroidal magnetic field (B t ) direction at a relatively low B powder flow rate, one order of magnitude lower than that at the outer target. In a similar scenario with unfavourable B t , it is found that the detachment thresholds of B flow rate for both the inner and outer targets are close and of the same order as that for the outer target with favourable Bt. In favourable B t direction at B powder flow rate of 1.2 × 10 21 atoms/s, a localized, broadened high-density region is formed near the inner target benefitted by the injection location and the E × B drift, and a radiation-intensified zone, mostly contributed by B 1+ and B 2+ , occurs there. The E × B drift facilitates plasma detachment at the inner target and simultaneously amplifies the in–out divertor asymmetry. In addition, the simulation results with three different injection locations show that the injection from outer strike point leads to the lowest Z eff inside the separatrix and has an intermediate flow rate for detachment at the outer target, comparing with the X-point and upstream locations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Accurate determination of solvation free energies of neutral organic compounds from first principles

The main goal of molecular simulation is to accurately predict experimental observables of molecular systems. Another long-standing goal is to devise models for arbitrary neutral organic molecules with little or no reliance on experimental data. While separately these goals have been met to various degrees, for an arbitrary system of molecules they have not been achieved simultaneously. For biophysical ensembles that exist at room temperature and pressure, and where the entropic contributions are on par with interaction strengths, it is the free energies that are both most important and most difficult to predict. We compute the free energies of solvation for a diverse set of neutral organic compounds using a polarizable force field fitted entirely to ab initio calculations. The mean absolute errors (MAE) of hydration, cyclohexane solvation, and corresponding partition coefficients are 0.2 kcal/mol, 0.3 kcal/mol and 0.22 log units, i.e. within chemical accuracy. The model (ARROW FF) is multipolar, polarizable, and its accompanying simulation stack includes nuclear quantum effects (NQE). The simulation tools’ computational efficiency is on a par with current state-of-the-art packages. The construction of a wide-coverage molecular modelling toolset from first principles, together with its excellent predictive ability in the liquid phase is a major advance in biomolecular simulation.

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

AtomAI framework for deep learning analysis of image and spectroscopy data in electron and scanning probe microscopy

Over the past several decades, electron and scanning probe microscopes have become critical components of condensed matter physics, materials science and chemistry research. At the same time, the infrastructure for establishing a connection between microscopy observations and materials behaviour over a broader parameter space is lacking. In this work, we introduce AtomAI, an open-source software package bridging instrument-specific Python libraries, deep learning and simulation tools into a single ecosystem. AtomAI allows direct applications of deep neural networks for atomic and mesoscopic image segmentation converting image and spectroscopy data into class-based local descriptors for downstream tasks such as statistical and graph analysis. For atomically resolved imaging data, the output is types and positions of atomic species, with an option for subsequent refinement. AtomAI further allows the implementation of a broad range of image and spectrum analysis functions, including invariant variational autoencoders for disentangling structural factors of variation and im2spec type of encoder–decoder models for mapping structure–property relationships. Finally, our framework allows seamless connection to the first principles modelling with a Python interface on the inferred atomic positions.

36 MATERIALS SCIENCE↗
DOE OSTI2023

WindSE: Wind Systems Engineering [SWR-17-03]

WindSE is a python package that uses a FEniCS backend to perform wind farm simulations and optimization. Documentation can be found at: https://windse.readthedocs.io/en/latest/

King, Ryan↗
DOE OSTI2025

Multiscale Machine-Learned Modeling Infrastructure

The Multiscale Machine-Learned Modeling Infrastructure (MuMMI) is a multiscale workflow management infrastructure that can concurrently orchestrate thousands of molecular dynamics (MD) simulations operating at different time and/or length scales, spanning nanoseconds to seconds and nanometers to micrometers. MuMMI uses machine learning (backed by biology experiments) to guide a massive ensemble of MD simulations that capture biologically relevant time and length scales with unprecedented resolution. MuMMI supports multiple MD codes such as GROMACS and ddcMD and can be fully deployed using the HPC package manager Spack. MuMMI has been used in many publications to run hundreds of thousands simulations, leading to significant biology breakthroughs.

Di Natale, Francesco [Lawrence Livermore National ↗
DOE OSTI2022

MCNP6.3 Unstructured Mesh Verification: GodivR and CANDU Models

A geometric cell of the Monte Carlo N-Particle (MCNP)1 transport code is traditionally created by using Boolean operators on defined surfaces. This constructive solid geometry (CSG) capability has been available in the MCNP code since its beginning. However, a CSG model approach is limited when it comes to constructing a representative geometry for a complex model in its ability to capture a correct model representation. Starting with the version 6.0, the MCNP code has the ability of embedding an unstructured mesh (UM) model into a CSG cell to create a hybrid geometry [1]. The MCNP UM feature provides the flexibility of defining very complex geometries because computer aided design (CAD) and mesh generation software packages can be utilized to construct UM models for MCNP simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Data-driven time-scale separation of ODE right-hand sides using dynamic mode decomposition and time delay embedding

Multi-physics simulation often involves multiple different scales. The ARKODE ODE solver package in the SUNDIALS library addresses multi-scale problems with a multi-rate time-integrator that can work with a right-hand side that has fast scale and slow scale components. In this report, we use dynamic mode decomposition and time delay embedding to extract the fast and slow components of the right-hand sides of a simple ODE from data. We then use the extracted components to solve the ODE with ARKODE. Finally, to move towards a real-world use case, we attempt to extract fast and slow scale dynamics from synthetic seismic modeling data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Status of PyGriffin Development for Integration into the NEAMS Workbench

The integration of Griffin into Workbench, by way of the PyGriffin code package, was initiated by the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program in FY-2022. PyGriffin was developed as a Python wrapper to Griffin to streamline the complex workflow involving mesh, cross section, and Griffin input generation, code execution, and results postprocessing, by improving user experience with high fidelity neutronics analysis through the Workbench GUI interface. PyGriffin can be used as a standalone application or through the PyARC code suite. FY-2023 saw the continued development of PyGriffin to expand its capabilities including advanced post-processing of simulation results and initial development of the Monte Carlo (MC) cross section generation workflow with the Shift MC code. PyGriffin was also approved for open-source software status in FY-2023, initiating plans to move the code repository into an open location. In addition to PyGriffin development, there were many improvements and developments made in PyARC leading to several software releases (latest version is v2.3.0), but the focus of this report is PyGriffin.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗