Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation and Modeling Tools”

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 685 records · Page 38

Using neutrons and x rays to measure plasma conditions in a solid sphere of deuterated polyethylene compressed to densities of 35 g/cc at temperatures of 2 keV and pressures of 40 Gbar

This paper describes an experiment that shock compresses the center of a solid deuterated polyethylene sphere, CD 2 , to densities of 35 g/cc and temperatures of 2 keV with corresponding pressure of 40 Gbar. The design employs a strong spherically converging shock launched through a solid ball of material using a Hohlraum radiation drive. As the shock coalesces at the center it produces a hot spot that we characterize by measuring the x-ray self-emission and 2.45 MeV neutrons emitted. Two-dimensional images and time-resolved measurements of the x rays emitted determine the size and time duration of the hot spot, leading to an estimated 2k eV electron temperature. The neutron time of flight spectrometer measures an average ion temperature of 1.06 ± 0.15 keV and neutron yield of 7.0 (±0.5) × 10 9 DD neutrons. Our new distribution function tool enables us to create a forward model of the experimental data based on 1D radiation-hydrodynamic simulations, leading to a better understanding of the plasma conditions that produce the measured neutrons and x rays. Furthermore, our simulations indicate that the x rays are produced in a short-lived hot-dense core over tens of picoseconds, whereas the neutron emission continues for about 200 ps, as the hot core starts to expand, thereby leading to a lower mean temperature of the plasma during neutron production. This finding is in agreement with the experimental data, and we therefore conclude that the forward-modeling is a useful tool for inferring the conditions of the hot spot in a laser-driven implosion during burn.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Report of the Atmosphere to Electrons Land-Based Mesoscale-to-Microscale Coupling Project (FY2020)

The overall goal of the Mesoscale-to-Microscale Coupling (MMC) project is to improve coupling between mesoscale and microscale simulations via improved guidance and new strategies for setting up simulations and for the development of new tools that can be used across the community. Including the mesoscale forcing is critical to modeling the full energy transfer across scales in the atmosphere. The project-specific objectives include: 1) To apply rigorous verification and validation (V&V) techniques to the new modeling tools that are developed as part of the project to ensure the accuracy of our codes and results and to develop estimates of the relative uncertainty. 2) To improve computational performance of the coupled MMC models through the development of methods that can be used to reduce turbulence spin-up time and hence the size of computational domains. 3) To improve representation of the surface layer in microscale models to enhance simulations of hub-height wind speed. 4) To develop guidance for the community describing the best ways to couple mesoscale and microscale models, including specific spatial scales at which the handoff to the microscale model should occur. 5) To prepare documentation and a suite of software tools that can be used across the community. And 6) To transition MMC research to the offshore environment. Major progress was made in each of these areas during FY20. The land-based portion of the project is reported herein, while the offshore portion will be reported separately. The team continued to advance the MMC tools and methodologies as well as to document their performance in journal papers and conference presentations, although several planned conferences were cancelled due to the COVID-19 pandemic.

54 ENVIRONMENTAL SCIENCES↗

Predicting fracture behavior in single crystal nickel using a coupled crystal plasticity phase field damage approach

Understanding the fracture behavior of single crystal metals is critical for predicting material performance under mechanical loading. Here, in this study, we investigate the fracture characteristics of single crystal nickel tensile bars using a crystal plasticity coupled phase field damage (CP-PFD) model. Experimental tensile tests were conducted on 15 specimens spanning five crystallographic orientations and three thickness variants per orientation. The results revealed two distinct fracture modes: brittle fractures with 45-degree angled surfaces and ductile fractures characterized by significant necking. The CP-PFD model successfully replicated these fracture behaviors, demonstrating strong agreement with experimental observations. The model effectively predicted the strain at which necking and fracture occurred, as well as the orientation-dependent fracture mechanisms. By comparing experimental and simulated fracture surfaces, we establish the CP-PFD model as a robust tool for predicting single crystal behavior and damage evolution. This work provides insight into the microstructural dependence of fracture behavior and establishes a predictive framework for modeling orientation-dependent damage evolution in single-crystal nickel.

Crystal plasticity↗

Large‐Eddy Simulations of a Convection Cloud Chamber: Sensitivity to Bin Microphysics and Advection

Abstract Bin microphysics schemes are useful tools for cloud simulations and are often considered to provide a benchmark for model intercomparison. However, they may experience issues with numerical diffusion, which are not well quantified, and the transport of hydrometeors depends on the choice of advection scheme, which can also change cloud simulation results. Here, an atmospheric large‐eddy simulation model is adapted to simulate a statistically steady‐state cloud in a convection cloud chamber under well‐constrained conditions. Two bin microphysics schemes, a spectral bin method and the method of moments, as well as several advection methods for the transport of the microphysical variables are employed for model intercomparison. Results show that different combinations of microphysics and advection schemes can lead to considerable differences in simulated cloud properties, such as cloud droplet number concentration. We find that simulations using the advection scheme that suffers more from numerical diffusion tends to have a smaller droplet number concentration and liquid water content, while simulation with the microphysics scheme that suffers more from numerical diffusion tends to have a broader size distribution and thus larger mean droplet sizes. Sensitivities of simulations to bin resolution, spatial resolution, and temporal resolution are also tested. We find that refining the microphysical bin resolution leads to a broader cloud droplet size distribution due to the advection of hydrometeors. Our results provide insight for using different advection and microphysics schemes in cloud chamber simulations, which might also help understand the uncertainties of the schemes used in atmospheric cloud simulations.

54 ENVIRONMENTAL SCIENCES↗

Surrogate modelling the Baryonic Universe – I. The colour of star formation

The spectral energy distribution of a galaxy emerges from the complex interplay of many physical ingredients, including its star formation history (SFH), metallicity evolution, and dust properties. Using GALAXPY , a new galaxy spectral prediction tool, and SFHs predicted by the empirical model UNIVERSEMACHINE and the cosmological hydrodynamical simulation IllustrisTNG, we isolate the influence of SFH on optical and near-infrared colours from 320 to 1080 Å at z = 0. By carrying out a principal component analysis, we demonstrate that physically motivated SFH variations modify galaxy colours along a single direction in colour space: the SFH-direction. We find that the projection of a galaxy’s present-day colours on to the SFH-direction is almost completely regulated by the fraction of stellar mass that the galaxy formed over the last billion years. Together with cosmic downsizing, this results in galaxies becoming redder as their host halo mass increases. We additionally study the change in galaxy colours due to variations in metallicity, dust attenuation, and nebular emission lines, finding that these properties vary broad-band colours along distinct directions in colour space relative to the SFH-direction. Finally, we show that the colours of low-redshift Sloan Digital Sky Survey galaxies span an ellipsoid with significant extent along two independent dimensions, and that the SFH-direction is well-aligned with the major axis of this ellipsoid. Our analysis supports the conclusion that variations in SFH are the dominant influence on present-day galaxy colours, and that the nature of this influence is strikingly simple.

79 ASTRONOMY AND ASTROPHYSICS↗

Accelerating Scientific Simulations with Bi-Fidelity Weighted Transfer Learning

High-fidelity modeling is an essential design tool for many engineering applications. However, for complex systems, computational cost can be a limiting factor. Analyzing parameter sensitivity, uncertainty quantification, and design optimization require many model evaluations. Surrogate models are often used to develop the relationship between model parameters and quantities of interest. However, in the case of complex systems, surrogate models require several degrees of freedom and, thus, a large number of data points to determine the correct dependencies. For many applications, this may be prohibitively expensive. The reduction of computational requirements can be achieved by leveraging low-fidelity models. Low-fidelity models represent the system at a coarser resolution with the advantage of computational efficiency. Therefore, a bi-fidelity modeling paradigm, which augments the accuracy of a low-fidelity model in a computationally efficient manner by invoking limited runs of a high-fidelity model, can be leveraged to sufficiently balance the accuracy and computational requirements. In this work, a bi-fidelity weighted transfer learning method using neural networks was applied to a computational fluid dynamics heat transfer modeling problem. The transfer learning advantage was investigated as a function of hyperparameters. Our main finding is that the use of a bi-fidelity modeling paradigm achieves accuracy close to that of a high-fidelity Gaussian process model while significantly reducing computational cost. The bi-fidelity model achieves comparable performance with 90 high-fidelity samples-that is, 60% less than the samples needed to achieve similar accuracy without the use of bi-fidelity modeling,

Borowiec, Katarzyna↗

Development of an Optimal Variable-Pitch Controller for Floating Axial-Flow Marine Hydrokinetic Turbines

This article discusses the development of an optimal variable-pitch controller for floating, axial-flow marine turbines. Recently, OpenFAST, an open-source wind turbine modeling tool, has been extended to model marine turbines. A controller is necessary to simulate marine turbines for different load cases using OpenFAST, which greatly impacts the performance of the energy system. Previous studies have designed controllers using a linearized model of the marine turbine, which can be time-consuming and require the expertise of a control engineer. In this study, we use an automated approach that uses generic models of the marine turbine to identify the controller gains, which can expedite the process of designing a controller. Using an optimizer to identify the control system parameters can additionally improve the controller's performance. The optimal controller tuned using such an approach results in a 20% reduction in the tower-base damage equivalent loading and better tracking of the rated generator speed and power.

axial flow↗

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↗

NEML2 - THE NEW ENGINEERING MATERIAL MODEL LIBRARY, VERSION 2

SF-22-157 NEML2 is an extension of an earlier open source code developed at ANL (SF-xxxx) called NEML. These libraries provide tools for constructing and using models of engineering materials, with a focus on high temperature materials for nuclear reactor applications. The software include tools for setting up material models, calibrating models, storing models for later use in other simulations, and running simple "material point" simulations of the material behavior. The main difference between NEML and NEML2 is the new software has been complete restructured to run on GPUs and other accelerators in order to improve the performance of methods used to train constitutive models against data and to improve the performance of the final models in production environments, like finite element software. While over time we expect NEML2 to become more-or-less feature complete with NEML we are filing a new software disclosure because (1) the library has been complete restructured and (2) we want to host NEML2 in a separate repository.

MESSNER, MARK↗

MOOSE Reactor Module Meshing Enhancements to Support Reactor Physics Analysis

The U.S. Department of Energy Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) program develops an integrated suite of advanced reactor physics tools built upon the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. Each code generally requires an input finite element mesh on which the physics solution is calculated, reported, and transferred to other physics codes. The meshing process is often burdensome for the complex geometries present in reactors due to lack of easy-to-use, open-source meshing tools. To address the bottleneck associated with meshing complex geometries found in nuclear reactors, NEAMS initiated the development of the MOOSE Reactor Module starting in FY21. The Reactor Module builds off the existing MOOSE Mesh System to include targeted meshing capabilities such as the ability to generate hexagonal pin cells, assemblies with ducts, rotating control drums, cores, peripheral zones around a core, as well as the automatic labeling (“reporting IDs”) of pin, assembly, and planar zones to simplify post-processing of results. As a Physics Module in MOOSE, the Reactor Module is open-source, available with any MOOSE installation, directly compatible with MOOSE-based tools, and can be invoked from MOOSEbased applications to generate meshes. Functionality from the Reactor Module has been applied to several advanced reactor concepts to demonstrate user workflow improvements and accuracy. The primary objective of the Reactor Module is to improve useability of MOOSEbased tools by streamlining mesh generation and output inspection processes. During FY22, the functionality of the Reactor Module (and accompanying Mesh System) has been expanded based on user needs. First, the Reactor Geometry Mesh Builder capability developed primarily in FY21 has been refactored and merged to the public MOOSE repository. This capability wraps underlying Reactor Module mesh generators into a “Pin – Assembly – Core” workflow appropriate for conventional Cartesian and hexagonal geometries, and notably assigns material IDs during mesh generation stage and generates only the minimal number of blocks needed in order to reduce computational burden. Biasing and boundary layer options have been added to the base mesh generators as required by thermal hydraulics solvers. The reporting ID functionality has been expanded to differentiate ring-wise and azimuthal sectors within a pin for use with depletion algorithms, and VectorPostProcessor and Reporter objects are now available to integrate solution variables across zones based on ID combinations. Functionality to trim hexagonal meshes along the center or periphery has been developed so users may leverage symmetry and reflective boundary conditions to reduce the mesh size. A flexible and powerful tool to fill the space between two sidesets has been introduced to the framework and can be used for transition layers such as stitching two assemblies together with different numbers of pins, or for complex geometries which do not follow conventional Cartesian/hexagonal patterns. Finally, additional verification problems were performed with NEAMS physics tools in complement with existing NEAMS work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Simulation Tools for Characterizing Stress Distribution in Laser Welded Dissimilar Joints

This project focuses on developing a thermo-metallurgical-mechanical modeling method to accurately predict the microstructural evolution and residual stress in laser welding between dissimilar metals, such as HSLA steel and high carbon equivalent (CE) gear steel. The method leverages a comprehensive material database to model the temperature and rate dependent phase transformations, along with their associated effects on material properties, such as thermal expansion and flow stress, throughout the welding process. A key innovation is the incorporation of phase transformation and phase-specific properties, which enhances the accuracy of residual stress predictions. The mixture material in the fusion zone due to the dissimilar metals will also be addressed in the numerical model. This is especially critical in scenarios involving phase transformations in the fusion zone and heat-affected zone (HAZ), where the phase changes can induce substantial residual stress variations. The material database has been generated using JMatPro. The modeling approach is implemented through a custom User Material (UMAT) subroutine, executed with the commercial finite element software Abaqus.

36 MATERIALS SCIENCE↗

Simulation Tools for Characterizing Stress Distribution in Laser Welded Dissimilar Joints

This project focuses on developing a thermo-metallurgical-mechanical modeling method to accurately predict the microstructural evolution and residual stress in laser welding between dissimilar metals, such as HSLA steel and high carbon equivalent (CE) gear steel. The method leverages a comprehensive material database to model the temperature and rate dependent phase transformations, along with their associated effects on material properties, such as thermal expansion and flow stress, throughout the welding process. A key innovation is the incorporation of phase transformation and phase-specific properties, which enhances the accuracy of residual stress predictions. The mixture material in the fusion zone due to the dissimilar metals will also be addressed in the numerical model. This is especially critical in scenarios involving phase transformations in the fusion zone and heat-affected zone (HAZ), where the phase changes can induce substantial residual stress variations. The material database has been generated using JMatPro. The modeling approach is implemented through a custom User Material (UMAT) subroutine, executed with the commercial finite element software Abaqus.

36 MATERIALS SCIENCE↗

Analysis Program (2023 Annual Progress Report)

This document summarizes the progress of VTO Analysis projects supported during the fiscal year 2023. The VTO Analysis Program supports mission-critical technological, economic, and interdisciplinary analyses to assist in prioritizing VTO technology investments and to inform research portfolio planning. These efforts provide essential vehicle and market data, modeling and simulation, and integrated and applied analyses, using the unique capabilities, analytical tools, and expertise resident in the DOE’s national laboratory system. VTO Analysis projects also demonstrate additional capabilities and expertise provided by research partnerships that may include academia, the private sector, and non-profit organizations.

33 ADVANCED PROPULSION SYSTEMS↗

Understanding and Predicting Pedestal Scenarios on NSTX-U (Final Technical Report)

This project investigated gyrokinetic instabilities in the NSTX pedestal, identified the major transport mechanisms (some of them novel) in the NSTX pedestal, and developed and validated a predictive modeling capability for pedestal transport in spherical tokamaks, with particular emphasis on NSTX discharges. The work combined first-principles gyrokinetic simulations, reduced transport models, integrated transport calculations with ASTRA, and exploratory machine learning tools. The central outcome is a practical modeling capability for pedestal temperature profiles based on reduced models informed and constrained by gyrokinetic physics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Hybrid Energy-Powered Electrochemical Direct Ocean Capture Model

Offshore synthetic fuel production and marine carbon dioxide removal can be enabled by direct ocean capture, which extracts carbon dioxide from the ocean that then can be used as a feedstock for fuel production or sequestered underground. To maximize carbon capture, plants require a variety of low-carbon energy sources to operate, such as variable renewable energy. However, the impacts of variable power on direct ocean capture have not yet been thoroughly investigated. To facilitate future deployments, a generalizable model for electrodialysis-based direct ocean capture plants is created to evaluate plant performance and electricity costs under intermittent power availability. This open-source Python-based model captures key aspects of the electrochemistry, ocean chemistry, post-processing, and operation scenarios under various conditions. To incorporate realistic energy supply dynamics and cost estimates, the model is coupled with the National Renewable Energy Laboratory’s H2Integrate tool, which simulates hybrid energy system performance profiles and costs. This integrated framework is designed to provide system-level insights while maintaining computational efficiency and flexibility for scenario exploration. Initial evaluations show similar results to those predicted by the industry, and demonstrate how a given plant could function with variable power in different deployment locations, such as with wind energy off the coast of Texas and with wind and wave energy off the coast of Oregon. The results suggest that electrochemical systems with greater tolerances for power variability and low minimum power requirements may offer operational advantages in variable-energy contexts. However, further research is needed to quantify these benefits and evaluate their implications across different deployment scenarios.

16 TIDAL AND WAVE POWER↗

NUMERICAL PREDICTIONS OF MEAN PERFORMANCE AND DYNAMIC BEHAVIOR OF A 10 MWe SCO2 COMPRESSOR WITH TEST DATA VALIDATION

High-fidelity aerodynamic analysis has been demonstrated for mean performance and unsteady dynamics in a sCO2 compressor designed by Hanwha Power Systems Americas for a 10 MWe Concentrating Solar Power plant. Simulations were performed with CRUNCH CFD® software tool that was matured to accurately model near critical real fluid effects in sCO2. Pre-test predictions for mean performance were validated with test data collected later. Performance predictions were accurate and captured sensitivity of the efficiency to inlet temperature of CO2 as well as steep drop-off at high flow rates due to condensation in the inlet throat. Detailed analysis was performed to understand the source of these performance losses at near critical conditions. Unsteady dynamic effects in the compressor at off-design conditions were also identified and quantified. In particular a system wide “condensation surge” condition was detected at high flow coefficients that results in large amplitude pulsations with accompanying mass flow fluctuations at low frequencies and has potential to cause damage in closed loop systems.

hosangadi, ashvin↗

Pathways to a Sustainable Aviation Ecosystem: Aviation Energy Research and Operation Simulator (AEROSim)

The National Renewable Energy Laboratory (NREL) has developed highly accessible tools for ground vehicle modeling, making it easy to quickly compare powertrains and estimate the impact of technology improvements on fuel/energy consumption, performance, cost, and battery life. With AEROSim, NREL advanced aircraft, enabling the same level of ease and versatility in comparing aircraft designs across a range of scenarios and conditions.

advanced air mobility↗

Development of an Optimal Variable-Pitch Controller for Floating Axial-Flow Marine Hydrokinetic Turbines: Preprint

This article discusses the development of an optimal variable-pitch controller for floating, axial-flow marine turbines. Recently, OpenFAST, an open-source wind turbine modeling tool, has been extended to model marine turbines. A controller is necessary to simulate marine turbines for different load cases using OpenFAST, which greatly impacts the performance of the energy system. Previous studies have designed controllers using a linearized model of the marine turbine, which can be timeconsuming and require the expertise of a control engineer. In this study, we use an automated approach that uses generic models of the marine turbine to identify the controller gains, which can expedite the process of designing a controller. Using an optimizer to identify the control system parameters can additionally improve the controller's performance. The optimal controller tuned using such an approach results in a 20% reduction in the towerbase damage equivalent loading and better tracking of the rated generator speed and power.

closed-loop control↗