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.

SourceResearch connectionFrom yearThrough year

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

At least 865 records · Page 48

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Challenges and Solutions for Real-Time Phasor Modeling of Large-scale Distribution Network with High PV Penetration

The conversion process of a practical large-scale feeder data from a quasi-static time series (OpenDSS) model into a real-time phasor model (ePHASORSIM in Opal-RT) is discussed in the paper. The process is implemented using an open source Python software. Previous reported implementations for the conversion process lead to several errors when applied to a larger-scale system such as the one considered here. Hence in this work, we describe the common issues in this conversion and propose a customized solution to enhance the efficiency of the conversion and reduce the complexity in the process. A quantitative validation of the enhanced conversion process is presented in this work using an actual high PV penetration feeder model that consists of 2230 buses, and using actual load and PV profile data. After a detailed analysis, this customized conversion software will be made available as an open source tool and is expected to be helpful for researchers who want to pursue a similar conversion. Solutions to various observed issues such as identifying the lines due to islanded network, representation of full impedance model of transformer/lines as sequential models, complexity in the representation of single phase buses/lines as three phase buses/lines to make it compatible with the simulator platform are discussed. Comparison of power flow and time series simulation results obtained from both OpenDSS and ePHASORsim models show very low errors, validating the accuracy of the proposed conversion process.

14 SOLAR ENERGY↗

DOE OSTI2023

Exploring structural transitions at grain boundaries in Nb using a generalized embedded atom interatomic potential

The advancement in experimental techniques, like the atom probe tomography and high resolution electron microscopy, is fueling interest in studying structural transformations of grain boundaries in metal and alloys to uncover correlations between mechanical properties and solute or impurity segregation to grain boundaries. Atomistic modeling is an important tool that can pinpoint the intricate dynamics of grain boundary phase transitions, but the lack of accurate interatomic potentials needed to simulate the complex dynamics of grain boundary structural transitions and identify different metastable phases has been the bottleneck. To this end, we use niobium as a model body centered cubic (BCC) metal and develop an interatomic potential to study grain boundary phase transitions. The potential for Nb is based on a generalization of the embedded atomic method potential and has sufficient flexibility to learn complex energy landscapes using a small set of training structures. We systematically test and validate the using data from ab initio density functional theory calculations and experiments. Using this potential, we calculate energies of multiple symmetric-tilt grain boundaries spanning a wide range of misorientation angles. Additionally, we explore different metastable structures of the Σ 27(552) [$1\overline{1}0$] grain boundary and use molecular dynamic simulations to study the coexistence of metastable phases and grain boundary transitions at finite temperature.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Enabling kilometer-scale E3SM land model simulation over North America: A new integrated framework solution

This study introduces a novel framework designed to enhance the performance, scalability, and portability of the kilometer-scale E3SM Land Model (km-ELM) within the E3SM modeling infrastructure. By seamlessly integrating cutting-edge data tools, we address existing challenges such as slow performance, limited scalability, and difficulties in software integration in current data-driven ELM simulation over large geographic areas. Our innovative approach leverages the KiloCraft data toolkit to generate unified inputs for simulations ranging from a single-cite case, to a 72,083-cell regional case to a continental configuration encompassing 21.6 million land grid cells at a 1 km × 1 km resolution. We conduct extensive strong- and weak-scaling experiments on three state-of-the-art supercomputers, utilizing up to 100,800 CPU cores across 2400 compute nodes to evaluate end-to-end metrics including wall-clock time, simulation-years-per-day (SYPD), initialization costs, and I/O throughput. Our results reveal the land (LND) component’s efficient scaling, demonstrating near-ideal weak scaling and strong-scaling parallel efficiencies reaching up to 87% at 50,400 cores. We confirm portability and reproducibility through bitwise-equivalent outputs across different machines using identical inputs over supported machines. Notably, at extreme scales, we identify I/O as a critical bottleneck and that leads to effective solution with the SCORPIO/ADIOS stack. Collectively, these findings validate the deployment of km-ELM at a continental scale with high parallel efficiency and provide essential guidance on configuration, decomposition, and I/O settings for optimized kilometer-scale land simulations in E3SM. This work emphasizes the innovative design and practical solutions that enhance the operational capabilities of km-ELM, focusing on software performance and scalability while leaving detailed scientific evaluations of simulated land processes for future investigations.

E3SM land model (ELM), km-ELM, scalability, perfor↗

DOE OSTI2021

WaterTAP3 (The Water Technoeconomic Assessment Pipe-Parity Platform)

The Water Technoeconomic Assessment Pipe-Parity Platform (WaterTAP3) was developed under the National Alliance for Water Innovation (NAWI) to facilitate consistent technoeconomic assessments of desalination treatment trains. The WaterTAP3 is an analytically robust modeling tool that can be used to evaluate water technology cost, energy, environmental, and resiliency tradeoffs across different water sources, sectors, and scales. The model simulates steady-state water treatment train performance and costs including flow and constituent mass balance across unit processes, based on source water conditions, configurations of treatment technologies, and system-level techno-economic assumptions. Users can build a new treatment train by connecting any number of unit processes, specific for their context and system, or selecting a train from the treatment train library. The model contains various technical and cost parameter options for a range of treatment processes and a library of influent water quality characteristics for a variety of source waters and case studies. Users can customize water quality parameters to evaluate the technology performance in their context. The model can be set up for different assessment needs including simulation, optimization, and uncertainty and sensitivity analyses. The results from WaterTAP3 can help identify trade-offs among the different system performance metrics, with insight on how particular technologies or systems promote pipe-parity. The flexibility and comprehensive scope of the tool makes it a promising solution to industry-wide water technoeconomic evaluations, leading to more informed water investment decisions and technologies. As a user-friendly, open-source platform, WaterTAP3 can be used by industry, academia, policymakers, planners, and those with or without extensive analytical experience.

Miara, Ariel↗

DOE OSTI2021

Water Technoeconomic Assessment Pipe-Parity Platform (WaterTAP3)

Miara, Ariel↗

DOE OSTI2025

Multiscale Thermal-hydraulic analysis of the MARVEL micro-reactor using coupled MOOSE Subchannel (SCM) and SAM

MARVEL is a natural-convection-cooled sodium-potassium microreactor that is anticipated to generate 85 kilowatts of thermal energy. It will operate within Idaho National Laboratory’s Transient Reactor Test Facility and is being developed by the DOE Microreactor Program. MARVEL will be used to test microreactor applications, evaluate systems for remote monitoring, and develop autonomous control technologies. A thermal-hydraulic computational model of this facility is a valuable tool to study important transients and calculate the safety limits of the micro-reactor design. For this purpose, the authors have chosen to use a multiscale coupled simulation: SCM for modeling the reactor core and SAM for the reactor’s primary cooling system. SCM is MOOSE physics module for subchannel analysis, which was designed to model single-phase flows through liquid-metal cooled, wire-wrapped fuel pin sub-assemblies, ordered in a triangular lattice. The SCM code was modified to be able to model MARVEL’s unique geometry. SAM is a systems analysis module based on the MOOSE framework. It aims to provide fast-running, whole-plant transient analyses capability with improved-fidelity for various advanced reactor types. The coupling between the two SCM and SAM for MARVEL modeling is done implementing a domain over-lapping approach. The resulting coupled simulation can model transients such as reactor startup/shutdown and provide an intermediate fidelity picture of the temperature field and other variables, in the core. Results for the steady-state simulations are presented in the article as well as flow blockage transient.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2021

Scalable Programming Workflows for Validation of Quantum Computers

Hybrid quantum-classical workflows have become standard methods for executing variational algorithms and other quantum simulation techniques, which are key applications for noisy intermediate scale quantum (NISQ) computers. Validating these simulations is an important task which helps gauge the progress of quantum computer development, and classical simulation can serve as a tool to this end. Both exact and more scalable approximate methods with quantifiable error bounds can be used in validation tasks where the applicable metrics include the distance from a calculable ground truth, the quality of an error model fit to data, etc. Here we present a library extension that includes methods for validation of quantum simulations based on scalable hybrid workflows executable on high performance computers. We provide examples that use approximate methods based on tensor networks and stabilizer simulators to bound the error of quantum simulations on NISQ hardware.

Nguyen, Thien↗

DOE OSTI2024

A simple model for short-range ordering kinetics in multi-principal element alloys

Short-range ordering (SRO) in multi-principal element alloys influences material properties such as strength and corrosion. While some degree of SRO is expected at equilibrium, predicting the kinetics of its formation is challenging. We present a simplified isothermal concentration-wave (CW) model to estimate an effective relaxation time of SRO formation. Estimates from the CW model agree to within a factor of five with relaxation times obtained from kinetic Monte Carlo (kMC) simulations when above the highest ordering instability temperature. Further, the advantage of the CW model is that it only requires mobility and thermodynamic parameters, which are readily obtained from alloy mobility databases and Metropolis Monte Carlo simulations, respectively. The simple parameterization of the CW model and its analytical nature makes it an attractive tool for the design of processing conditions to promote or suppress SRO in multicomponent alloys.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Coarse-Grained Simulations for the Characterization and Optimization of Hybrid Protein–DNA Nanostructures

We present here the combination of experimental and computational modeling tools for the design and characterization of protein–DNA hybrid nanostructures. Our work incorporates several features in the design of these nanostructures: (1) modeling of the protein–DNA linker identity and length; (2) optimizing the design of protein–DNA cages to account for mechanical stresses; (3) probing the incorporation efficiency of protein–DNA conjugates into DNA nanostructures. The modeling tools were experimentally validated using structural characterization methods like cryo-TEM and AFM. Our method can be used for fitting low-resolution electron density maps when structural insights cannot be deciphered from experiments, as well as enable in-silico validation of nanostructured systems before their experimental realization. Finally, these tools will facilitate the design of complex hybrid protein–DNA nanostructures that seamlessly integrate the two different biomolecules.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

In situ synchrotron diffraction and modeling of non-equilibrium solidification of a MnFeCoNiCu alloy

The solidification mechanism and segregation behavior of laser-melted Mn 35 Fe 5 Co 20 Ni 20 Cu 20 was firstly investigated via in situ synchrotron x-ray diffraction at millisecond temporal resolution. The transient composition evolution of the random solid solution during sequential solidification of dendritic and interdendritic regions complicates the analysis of synchrotron diffraction data via any single conventional tool, such as Rietveld refinement. Therefore, a novel approach combining a hard-sphere approximation model, thermodynamic simulation, thermal expansion measurement and microstructural characterization was developed to assist in a fundamental understanding of the evolution of local composition, lattice parameter, and dendrite volume fraction corresponding to the diffraction data. This methodology yields self-consistent results across different methods. Via this approach, four distinct stages were identified, including: (I) FCC dendrite solidification, (II) solidification of FCC interdendritic region, (III) solid-state interdiffusion and (IV) final cooling with marginal diffusion. It was found out that in Stage I, Cu and Mn were rejected into liquid as Mn 35 Fe 5 Co 20 Ni 20 Cu 20 solidified dendritically. During Stage II, the lattice parameter disparity between dendrite and interdendritic region escalated as Cu and Mn continued segregating into the interdendritic region. After complete solidification, during Stage III, the lattice parameter disparity gradually decreases, demonstrating a degree of composition homogenization. The volume fraction of dendrites slightly grew from 58.3 to 65.5%, based on the evolving composition profile across a dendrite/interdendritic interface in diffusion calculations. Postmortem metallography further confirmed that dendrites have a volume fraction of 64.7 ± 5.3% in the final microstructure.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Exploring scenarios for enhanced fuel compression and performance on the National Ignition Facility with machine-learning-aided design techniques

Recent fusion experiments on the National Ignition Facility (NIF) have achieved ignition, producing multi-MJ fusion yields for input laser energies of roughly 2 MJ [Abu-Shawareb et al., Phys. Rev. Lett. 132, 065102 (2024)]. Building on the success of the target designs that have achieved ignition, we explore new implosion scenarios predicted to generate significantly more compression of the dense DT ice layer and correspondingly higher yields while preserving many of the key physics characteristics of present-day ignition designs. Our main result is a novel 3-shock implosion scheme that effectively minimizes the shock-induced entropy in the dense, accelerating DT shell and maximizes the resulting fuel compression subject to a fixed leading shock strength consistent with present-day ignition experiments, which is necessary to melt the crystalline high-density carbon ablator. Compared to the first NIF experiment to fulfill Lawson's ignition criterion, shot N210808 [Abu-Shawareb et al., Phys. Rev. Lett. 129, 075001 (2022)], our design exhibits a 40% increase in simulated peak areal density (ρR) and a 5× increase in 1D fusion yield using a 4% lighter ablator and identical DT payloads. We also present a complete integrated 2D hohlraum design and laser pulse specifications capable of generating the desired 3-shock drive and maintaining control of the low-mode capsule implosion symmetry, where the increase in simulated 2D yield relative to N210808 is > 10×. This new implosion regime was discovered with help from a machine-learning-enabled capsule design optimization framework. We outline the workflow this automated tool uses to identify improved design candidates by running several rounds of capsule simulations, constructing a surrogate model mapping input variations to key physics output quantities, and querying the resulting statistical model to propose adjustments to the x-ray drive and capsule to reach a set of physics objectives prescribed by the designer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

ROM-Based Surrogate Systems Modeling of EBR-II

We report the System Analysis Module (SAM), developed and maintained by Argonne National Laboratory, is designed to provide whole-plant transient safety analysis capabilities for a number of advanced non-light water reactors, including sodium-cooled fast reactor (SFR), lead-cooled fast reactor (LFR), and molten salt reactor (MSR)/fluoride-salt-cooled high-temperature reactor (FHR) designs. SAM is primarily constructed as a systems-level analysis tool, with the potential to incorporate reduced order models from three-dimensional computational fluid dynamics (CFD) simulations to improve characterization of complex, multidimensional physics. It is recognized that the computational expense associated with CFD can be intractable for various engineering analyses, such as uncertainty quantification, inference, and design optimization. This paper explores the reducibility of a SAM model using recent advances in randomized linear algebra techniques, which attempt to find recurring patterns in the various realizations generated by a model after randomly perturbing all its input parameters. The reduction is described in terms of fewer degrees of freedom (DOFs), referred to as the active DOFs, for the model variables such as input model parameters and model responses. The results indicate that there is significant room for additional reduction that may be leveraged for additional computational gains when employing SAM for engineering-intensive analyses that require repeated model executions. Different from physics-based reduction approaches, the proposed approach allows one to estimate upper bounds on the reduction errors, which are rigorously developed in this work. Finally, different methods for surrogate model construction, such as regression and neural network-based training, are employed to correlate the input and output active DOFs, which are related back to the original variables using matrix-based linear transformations.

42 ENGINEERING↗

DOE OSTI2023

Stochastic and mixed density functional theory within the projector augmented wave formalism for simulation of warm dense matter

Stochastic density functional theory (DFT) and mixed stochastic-deterministic DFT are burgeoning approaches for the calculation of the equation of state and transport properties in materials under extreme conditions. In the intermediate warm dense matter regime, a state between correlated condensed matter and kinetic plasma, electrons can range from being highly localized around nuclei to delocalized over the whole simulation cell. The plane-wave basis pseudopotential approach is thus the typical tool of choice for modeling such systems at the DFT level. Unfortunately, stochastic DFT methods scale as the square of the maximum plane-wave energy in this basis. To reduce the effect of this scaling and improve the overall description of the electrons within the pseudopotential approximation, we present stochastic and mixed DFT approaches developed and implemented within the projector augmented wave formalism. In conclusion, we compare results between the different DFT approaches for both single-point and molecular dynamics trajectories and present calculations of self-diffusion coefficients of solid density carbon from 1 to 50 eV.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

NMACFoam (Numerical-model for Membrane-spacer Assemblies and Configurations in OpenFOAM) [SWR-23-57]

NMACFoam is a computational fluid dynamics tool for simulating flow and solids transport in membrane-spacer channels pertaining to reverse osmosis or nanofiltration applications. It is developed on the OpenFOAM framework and allows for simulating spacers with complex geometries. It includes porous wall boundary conditions for membranes and can accurately capture concentration polarization effects. One can perform large calculations at module-scale domain sizes using high performance computing resources with NMACFoam, which scales to 1000s of processors for problems with 1-10 million cells.

Sitaraman, Hariswaran↗

DOE OSTI2026

Ambient and Initial Temperature Effects on Energy Consumption Rate Modeled in FASTSim

Ambient and initial temperatures significantly impact the energy consumption rate (ECR) of battery electric vehicles (BEVs) due to auxiliary loads and the temperature dependence of battery efficiency. This study introduces a streamlined, physics-based thermal modeling approach within the FASTSim tool that bridges the gap between oversimplified constant-load models and computationally expensive high-fidelity simulations. By employing a lumped thermal mass framework, the model captures fundamental energy balances and critical non-linear energy penalties while maintaining the computational efficiency required for expansive sensitivity studies. The simulations evaluated a compact BEV hatchback with a resistive heater over city (UDDS) and highway (HWFET) test cycles. Compared to a 22 degrees Celsius initial and ambient temperature baseline, a -7 degrees Celsius initial/ambient temperature resulted in a 221% increase in the ECR for the city cycle and a 100% increase for the highway cycle. Conversely, a 45 degrees Celsius initial / 40 degrees Celsius ambient temperature resulted in a 40% increase for UDDS and an 18% increase for HWFET. These results demonstrate that while cold conditions impose the most severe energy penalties due to resistive heating, the impact is consistently more pronounced in city driving where auxiliary loads represent a larger proportion of total energy. This lightweight yet robust framework enables researchers to rapidly quantify BEV thermal sensitivity across diverse climates without the need for high-overhead simulation environments.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2024

Thermal-Mechanic Modeling of Fusion Components Using The MOOSE Framework

With fusion energy rapidly developing, there is a great need for an open source tool that can rapidly model different fusion tokamak designs. The open-source codes Multiphysics Object Oriented Simulation Environment (MOOSE), Tritium Migration Analysis Program Version 8 (TMAP8) , and Fusion ENergy Integrated Multiphysi-X (FENIX) are being used to demonstrate thermal-mechanic modeling of fusion components. In this study, an analysis is done on the fusion nuclear science facility (FNSF) blanket, and ITER toroidal magnet and diverter mono block. This study will lay the groundwork for fusion component analysis in the MOOSE framework and FENIX program.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Tiered Energy Use in Buildings (TEB)

This tool creates a set of building energy simulations and organizes them according to a tier-based model representing circuits that will be on or off depending on a microgrid controller (not included in this software). This Python based code uses the open-source RCBuildingSimulator that is under an amended MIT license. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525. SAND2021-2528 O

Villa, Daniel↗

DOE OSTI2021

Development of MCNP Training Modules for Safeguards Practitioners [Abstract]

The Monte-Carlo N-Particle (MCNP) software developed at LANL is the most widely used neutron transport code in the world. It is an essential tool for a variety of applications including detector development and design, nuclear fuel burnup simulation, criticality safety, and nondestructive assay system optimization. For this reason, it is indispensable within the safeguards and materials control & accountability (MC&A) communities. Multiple MCNP training courses have been created and taught over the last several decades by the MCNP development team at LANL, however there are no existing courses that cover specialized topics considered fundamental to NDA and safeguards models. To fill this gap, the MCNP team and Safeguards Science and Technology group at LANL have co-created a set of training modules customized to meet the specialized needs of the safeguards and MC&A communities. The basic modules cover concepts such as NDA system optimization, He-specific and other capture tallies, and tools for improved theoretical understanding. An advanced module was also created to cover topics including variance reduction for active interrogation simulations, use of the LANL MCNPTools post-processor, PTRAC (particle tracking) and list-mode data simulations, and fuel burnup simulations. The training modules teach to the latest and most state-of-the-art MCNP features and tools released by the development team at LANL and are intended to be taught jointly by the developers and safeguards experts. Ultimately, we hope that creation of these modules will serve to capture and convey the safeguards modeling and MCNP expertise at LANL, and that we will be able to share the modules more broadly with the MC&A and safeguards communities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗