Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “complex simulations”

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 1,711 records · Page 95
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Modeling receiver flux of commercial power tower concentrating solar power plants using ray tracing: a round-robin comparison of SolTrace, Solstice, and TieSOL

This study presents a multi-stage, cross-validation comparison of three software packages for Monte Carlo ray tracing (MCRT) applied to central tower concentrating solar power (CSP) systems. The three packages evaluated are: (1) SolTrace, an open-source tool developed by the National Renewable Energy Laboratory (NREL); (2) Solstice, an open-source program created by CNRS-PROMES and Meso-Star, with enhancements for CSP applications (called solsticepy) from the Australian National University; and (3) TieSOL, a commercial software developed by Tietronix. This investigation extends previous ray tracing comparisons by incorporating models of multi-facet heliostats within a commercial-scale solar field, taking into account zoned focal lengths and canting configurations. Receiver flux distributions were compared across the tools using a series of case studies, including single-heliostat scenarios, isolated blocking situations, and comprehensive full-field simulations. The case studies were designed to diagnose differences across the models at varying levels of complexity, and to identify and resolve discrepancies as additional parameters were introduced. Key factors examined in the analysis include sun positions, heliostat location, facet and canting focus, and aimpoint strategies. The comparison aims to improve the accuracy and reliability of these tools while providing benchmark cases for validating future optical modeling tools.

14 SOLAR ENERGY↗
DOE OSTI2024

Data-driven modeling of dislocation mobility from atomistics using physics-informed machine learning

Dislocation mobility, which dictates the response of dislocations to an applied stress, is a fundamental property of crystalline materials that governs the evolution of plastic deformation. Traditional approaches for deriving mobility laws rely on phenomenological models of the underlying physics, whose free parameters are in turn fitted to a small number of intuition-driven atomic scale simulations under varying conditions of temperature and stress. This tedious and time-consuming approach becomes particularly cumbersome for materials with complex dependencies on stress, temperature, and local environment, such as body-centered cubic crystals (BCC) metals and alloys. In this paper, we present a novel, uncertainty quantification-driven active learning paradigm for learning dislocation mobility laws from automated high-throughput large-scale molecular dynamics simulations, using Graph Neural Networks (GNN) with a physics-informed architecture. We demonstrate that this Physics-informed Graph Neural Network (PI-GNN) framework captures the underlying physics more accurately compared to existing phenomenological mobility laws in BCC metals.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Local molecular environment drives speciation and reactivity of ion complexes in concentrated salt solution

The speciation and reactivity of aqueous salt solutions are important for a wide variety of applications. However, application of this information is inhibited by broad disagreement about composition, mechanisms, and structure in concentrated solutions especially. Here, neutron diffraction with isotopic substitution measurements on aqueous zinc chloride solutions are used to calibrate molecular dynamics simulations that include effective electronic polarization. This allows us to probe the origin of speciation and reactivity of zinc chloride-water ion complexes, ZnCl x (H 2 O) y 2-x (x ≤ 4 and y ≤ 6), by comparing the reactivity of species in concentrated (4.5 m) and dilute (0.01 m) conditions. Within the concentrated solution, it is found that the extended solvation environment is dominated by solvent-separated ion complexes whose presence increases the free energy of activation for interconversion of species while simultaneously enhancing their thermodynamic stability. This concentration-dependent reactivity and stability suggests that other reactions, such as the nucleation of solid phases, will also be affected.

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

Performance-based simulants for Hanford radioactive waste treatment process testing

Radioactive wastes from nuclear fuel processing are stored in large underground storage tanks at the Hanford Nuclear Reservation. Treatment and remediation require that waste feed from the storage tanks be delivered to the Hanford Tank Waste Treatment and Immobilization Plant continually for the duration of the treatment mission. Furthermore, the complex physical and chemical processes required for this mission, the significant scale, and the hazardous nature of the waste necessitate the use of simulants in process testing. A new class of simulants, “performance-based simulants,” has been developed to match the process performance of the simulant to actual waste performance data.

Hanford tank waste treatment↗
DOE OSTI2025

Defect modeling in semiconductors: the role of first principles simulations and machine learning

Abstract Point defects in semiconductors dictate their electronic and optical properties. Vacancies, interstitials, substitutional defects, and defect complexes can form in the semiconductor lattice and significantly impact its performance in applications such as solar absorption, light emission, electronics, and catalysis. Understanding the nature and energetics of point defects is essential for the design and optimization of next-generation semiconductor technologies. Here, we provide a comprehensive overview of the current state of research on point defects in semiconductors, focusing on the application of density functional theory (DFT) and machine learning (ML) in accelerating the prediction and understanding of defect properties. DFT has been instrumental in accurately calculating defect formation energies, charge transition levels, and other defect-related properties such as carrier recombination rates and lifetimes, and ion migration barriers. ML techniques, particularly neural networks, have emerged as powerful tools for enabling rapid prediction of defect properties at DFT-accuracy in order to overcome the expense of using large supercells and advanced functionals. We begin this article with a discussion of different types of point defects and complexes, their impact on semiconductor properties, and the experimental and DFT approaches typically used for their characterization. Through multiple case studies, we explore how DFT has been successfully applied to understand defect behavior across a variety of semiconductors, and how ML approaches integrated with DFT can efficiently predict defect properties and facilitate the discovery of new materials with tailored defect behavior. Overall, the advent of ‘DFT+ML’ promises to drive advancements in semiconductor technology, catalysis, and renewable energy applications, paving the way for the development of high-performance semiconductors which are defect-tolerant or have desirable dopability.

Rahman, Md Habibur (ORCID:000000027705984X)↗
DOE OSTI2024

Gray two-moment neutrino transport: Comprehensive tests and improvements for supernova simulations

Aims. In this work we extended an energy-integrated neutrino transport method to facilitate efficient, yet precise, modeling of compact astrophysical objects. We particularly focus on core-collapse supernovae. Methods. We implemented a gray neutrino-transport framework from the literature into FLASH and performed a detailed evaluation of its accuracy in core-collapse supernova simulations. Based on comparisons with results from simulations using energy-dependent neutrino transport, we incorporated several improvements to the original scheme. Results. Our analysis shows that our gray neutrino transport method successfully reproduces key aspects from more complex energy-dependent transport across a variety of progenitors and equations of state. We find both qualitative and reasonable quantitative agreement with multi-group M1 transport simulations. However, the gray scheme tends to slightly favor shock revival. In terms of gravitational wave and neutrino signals, there is a good alignment with the energy-dependent transport, although we find 15–30% discrepancies in the average energy and luminosity of heavy-lepton neutrinos. Simulations using the gray transport are around four times faster than those using energy-dependent transport.

Astronomy & Astrophysics↗
DOE OSTI2025

Simulation-Based Inference for Neutrino Interaction Model Parameter Tuning

High-energy physics experiments studying neutrinos rely heavily on simulations of their interactions with atomic nuclei. Limitations in the theoretical understanding of these interactions typically necessitate ad hoc tuning of simulation model parameters to data. Traditional tuning methods for neutrino experiments have largely relied on simple algorithms for numerical optimization. While adequate for the modest goals of initial efforts, the complexity of future neutrino tuning campaigns is expected to increase substantially, and new approaches will be needed to make progress. In this paper, we examine the application of simulation-based inference (SBI) to the neutrino interaction model tuning for the first time. Using a previous tuning study performed by the MicroBooNE experiment as a test case, we find that our SBI algorithm can correctly infer the tuned parameter values when confronted with a mock data set generated according to the MicroBooNE procedure. This initial proof-of-principle illustrates a promising new technique for next-generation simulation tuning campaigns for the neutrino experimental community.

Tame-Narvaez, Karla Maria [Fermilab]↗
DOE OSTI2025

Simulation-based inference for neutrino interaction model parameter tuning

High-energy physics experiments studying neutrinos rely heavily on simulations of their interactions with atomic nuclei. Limitations in the theoretical understanding of these interactions typically necessitate ad hoc tuning of simulation model parameters to data. Traditional tuning methods for neutrino experiments have largely relied on simple algorithms for numerical optimization. While adequate for the modest goals of initial efforts, the complexity of future neutrino tuning campaigns is expected to increase substantially, and new approaches will be needed to make progress. In this paper, we examine the application of simulation-based inference (SBI) to the neutrino interaction model tuning for the first time. Using a previous tuning study performed by the MicroBooNE experiment as a test case, we find that our SBI algorithm can correctly infer the tuned parameter values when confronted with a mock data set generated according to the MicroBooNE procedure. This initial proof-of-principle illustrates a promising new technique for next-generation simulation tuning campaigns for the neutrino experimental community.

Tame-Narvaez, Karla [Fermilab] (ORCID:000000022249↗
DOE OSTI2020

Direct Numerical Simulation of the Flow through a Randomly Packed Pebble Bed

The proposition for molten salt and high-temperature gas-cooled reactors has increased the focus on the dynamics and physics in randomly packed pebble beds. Research is being conducted on the validity of these designs as a possible contestant for the fourth-generation nuclear power systems. A detailed understanding of the coolant flow behavior is required in order to ensure proper cooling of the reactor core during normal and accident conditions. In order to increase the understanding of the flow through these complex geometries and enhance the accuracy of lower-fidelity modeling, high-fidelity approaches such as direct numerical simulation (DNS) can be utilized. Nek5000, a spectral-element computational fluid dynamics (CFD) code, was used to develop DNS fluid flow data. The flow domain consisted of 147 pebbles enclosed by a bounding wall. In the work presented, the Reynolds numbers ranged from 430 to 1050 based on the pebble diameter and inlet velocity. Characteristics of the flow domain such as volume averaged porosity, axial porosity, and radial porosity were studied and compared with correlations available in the literature. Friction factors from the DNS results for all Reynolds numbers were compared with correlations in the literature. The first- and second-order statistics show good agreement with the available experimental data. Turbulence length scales were analyzed in the flow. Reynolds stress anisotropy was characterized by utilizing invariant analysis. Overall, the results of the analysis in this study provide deeper understanding of the flow behavior and the effect of the wall in packed beds.

Yildiz, Mustafa Alper↗
DOE OSTI2024

A Gaussian Process-Based extended Goldak heat source model for finite element simulation of laser powder bed fusion additive manufacturing process

In this study, laser powder bed fusion (L-PBF) additive manufacturing (AM) is a key enabling technology to manufacture highly complex and integrated metallic structures. In L-PBF AM process, the melting of the metal powders and the layers underneath can be governed by either “conduction mode” or “keyhole mode”, with the keyhole mode reportedly leading to porosity and decreased strength and ductility by many studies. In part scale simulations, finite element (FE) model is often used to study the temperature distribution during printing and to predict the residual stress, where a volumetric heat flux with a Gaussian or a double ellipsoidal (Goldak) distribution is often applied as the laser heat source. However, the above heat source models can only capture the melt pool shape in the conduction mode, and fail to capture the transition to keyhole melting mode when the process parameters change. To overcome this inaccuracy, an extended Goldak heat source model is proposed by introducing a laser penetration term as a function of laser parameters obtained from a Gaussian-Process (GP) model. The model is validated by “2D pad” AlSi10Mg L-PBF experiments under a wide range of laser power, scan speed, and laser focus offset, and the results show the model successfully captures the measured melt pool shape in all conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Application of NEAMS Multiphysics Framework for Species Tracking in Molten Salt Reactors

This report from Idaho National Laboratory (INL) summarizes the key modeling and simulation activities conducted under the Department of Energy (DOE) Molten Salt Reactor (MSR) Campaign during the Fiscal Year 2023 (FY23). The focus of the work was to leverage state-of-the-art modeling capabilities from the DOE Nuclear Energy Advanced Modeling and Simulation (NEAMS) codes to enable novel multiphysics and multiscale modeling and simulation of MSRs. Through collaboration with NEAMS code developers, advanced multiphysics analysis capabilities for MSR systems were demonstrated by coupling depletion, thermal-hydraulics, and thermochemistry into an innovative framework for chemical species transport in MSRs. As a result, the framework can track nuclides throughout their lifetimes in the core, from production (depletion) to advection throughout the salt volume (thermal-hydraulics) and off-gassing or precipitation outside of the salt (thermochemistry). This work supports the near-term deployment of MSRs by integrating the synergistic efforts between the DOE’s MSR Campaign and NEAMS program. The resulting framework will help better connect system design modelers with experimentalists to better understand and predict complex physical behaviors in MSRs. Researchers and MSR developers alike can now leverage these new modeling and simulation capabilities to perform novel analyses with applications including: • MSR dynamics during normal operational transients and accident scenarios • Off-gas system design and performance for fuel cycle and depletion analysis • Corrosion and active chemistry control for reactor component health and lifetime determination • Source term, decay heat and activity determination in accident scenarios • Special nuclear material accountancy and chemical forensic analysis for safeguards • Digital twin development of experiments and experimental reactor demonstrations • Measurement requirements for instrumentation and control design • Uncertainty and sensitivity analysis of missing data to inform future experimental data collection.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Performance simulation of the soft gamma-ray concentrator

The soft gamma-ray concentrator is a telescope mission concept utilizing a suitable arrangement of bent multilayer structures of alternating low- and high-density materials. This lens is able to channel gamma-ray photons via total external reflection and concentrate the incident radiation to a point. The channeling technique offers the potential for concentrating gamma rays with focal lengths <10 m and energies >100 keV, beyond the reach of current grazing-incidence hard x-ray mirrors. For the performance estimation of such an instrument, we have developed a flexible set of computer modeling tools to compute the optical properties of multilayer structures, predict the channeling efficiency for a given multilayer configuration, and aid in the optimization of potential gamma-ray concentrator-based telescope designs. This modeling includes the multilayer optical properties calculated by the IMD software, the ray tracing using an IDL code, and the focal plane detector simulation by MEGAlib. We illustrate the potential of this approach by presenting simulated astronomical observations from a balloon-borne platform. The final result, including simulated effective area, instrument sensitivity, and polarization performance, shows that the gamma-ray concentrator will provide greatly increased sensitivity for next-generation soft gamma-ray missions with modest cost and complexity.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Simulation of Reverse Osmosis Membrane Compaction Using Material Point Method

Reverse Osmosis (RO) is a promising technology to address the impending water-crisis in the upcoming decades. RO at high salinities and high pressures is challenging due to membrane compaction that changes its porosity and permeability. In this study, we present a simulation methodology for membrane structural mechanics to understand pore size distribution and permeability variations under high pressure. We use the material point method (MPM), that solves the solid mechanics equations in a Lagrangian framework. MPM provides many features that make it well-suited for simulating mesoporous membranes. The Lagrangian framework allows for large deformations, easy integration of constitutive models and direct import of complex geometries as particles. The spatial discretization in our MPM solver is achieved using linear or cubic-spline shape functions while the time integration is carried out using the explicit Euler method. A series of images containing detailed pore structures obtained from X-ray tomography is first converted to a collection of material points. Compressive loads are applied to the top layer of the membrane to simulate the application of pressure. The membrane deformation and pore size distribution before and after load application are reported and compared with the experimental measurements. The presentation discusses the numerical methods used, the performance of the solver on high-performance computing machines, and the results of membrane compact in detail.

compaction↗
DOE OSTI2021

Impact of Blockchain Delay on Grid-Tied Solar Inverter Performance: Preprint

This paper investigates the impact of the delay resulting from a blockchain, a promising security measure, for a hierarchical control system of inverters connected to the grid. The blockchain communication network is designed at the secondary control layer for resilience against cyberattacks. To represent the latency in the communication channel, a model is developed based on the complexity of the blockchain framework. Taking this model into account, this work evaluates the plant’s performance subject to communication delays, introduced by the blockchain, among the hierarchical control agents. In addition, this article considers an optimal model-based control strategy that performs the system’s internal control loop. The work shows that the blockchain’s delay size influences the convergence of the power supplied by the inverter to the reference at the point of common coupling. In the results section, real-time simulations on OPAL-RT are performed to test the resilience of two parallel inverters with increasing blockchain complexity.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2022

Formation of Magnetically Truncated Accretion Disks in 3D Radiation-transport Two-temperature GRMHD Simulations

Multiwavelength observations suggest that the accretion disk in the hard and intermediate states of X-ray binaries (XRBs) and active galactic nucleus transitions from a cold, thin disk at large distances into a hot, thick flow close to the black hole (BH). However, the formation, structure, and dynamics of such truncated disks are poorly constrained due to the complexity of the thermodynamic, magnetic, and radiative processes involved. We present the first radiation-transport two-temperature general relativistic magnetohydrodynamic (GRMHD) simulations of truncated disks radiating at ~35% of the Eddington luminosity with and without large-scale poloidal magnetic flux. We demonstrate that when a geometrically thin accretion disk is threaded by large-scale net poloidal magnetic flux, it self-consistently transitions at small radii into a two-phase medium of cold gas clumps floating through a hot, magnetically dominated corona. This transition occurs at a well-defined truncation radius determined by the distance out to which the disk is saturated with magnetic flux. The average ion and electron temperatures in the semiopaque corona reach, respectively, $T_i$ ≳ 10 10 K and $T_e$ ≳ 5 × 10 8 K. The system produces radiation, powerful collimated jets, and broader winds at the total energy efficiency exceeding 90%, the highest ever energy extraction efficiency from a spinning BH by a radiatively efficient flow in a GRMHD simulation. This is consistent with jetted ejections observed during XRB outbursts. The two-phase medium may naturally lead to broadened iron line emission observed in the hard state.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Multifidelity deep operator networks for data-driven and physics-informed problems

Operator learning for complex nonlinear systems is increasingly common in modeling multi-physics and multi-scale systems. However, training such high-dimensional operators requires a large amount of expensive, high-fidelity data, either from experiments or simulations. In this work, we present a composite Deep Operator Network (DeepONet) for learning using two datasets with different levels of fidelity to accurately learn complex operators when sufficient high-fidelity data is not available. Additionally, we demonstrate that the presence of low-fidelity data can improve the predictions of physics-informed learning with DeepONets. We demonstrate the new multi-fidelity training in diverse examples, including modeling of the ice-sheet dynamics of the Humboldt glacier, Greenland, using two different fidelity models and also using the same physical model at two different resolutions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

SAXS-guided unbiased coarse-grained Monte Carlo simulation for identification of self-assembly nanostructures and dimensions

Recent studies have shown that solvated amphiphiles can form nanostructured self-assemblies called dynamic binary complexes (DBCs) in the presence of ions. Since the nanostructures of DBCs are directly related to their viscoelastic properties, it is important to understand how the nanostructures change under different solution conditions. However, it is challenging to obtain a three-dimensional molecular description of these nanostructures by utilizing conventional experimental characterization techniques or thermodynamic models. To this end, we combined the structural data from small angle X-ray scattering (SAXS) experiments and thermodynamic knowledge from coarse-grained Monte Carlo (CGMC) simulations to identify the detailed three-dimensional nanostructure of DBCs. Specifically, unbiased CGMC simulations are performed with SAXS-guided initial conditions, which aids us to sample accurate nanostructures in a computationally efficient fashion. As a result, an elliptical bilayer nanostructure is obtained as the most probable nanostructure of DBCs whose dimensions are validated by scanning electron microscope (SEM) images. Then, utilizing the obtained molecular model of DBCs, here we could also explain the pH tunability of the system. Overall, our results from SAXS-guided unbiased CGMC simulations highlight that using potential energy combined with SAXS data, we can distinguish otherwise degenerate nanostructures resulting from the inherent ambiguity of SAXS patterns.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Efficient Simulation of Logical Magic State Preparation Protocols

Developing space- and time-efficient logical magic state preparation (MSP) protocols will likely be an essential step toward building a large-scale fault-tolerant quantum computer. Motivated by this need, we introduce a scalable method for simulating logical MSP protocols under the standard circuit-level noise model. When applied to protocols based on code-switching, magic state cultivation, and magic state distillation, our method yields a complexity polynomial in (i) the number of qubits and (ii) the nonstabilizerness, e.g., stabilizer rank or Pauli rank, of the target encoded magic state. The efficiency of our simulation method is rooted in a curious fact: every circuit-level Pauli error in these protocols propagates to a Clifford error at the end. This property is satisfied by a large family of protocols, including those that repeatedly measure a transversal Clifford that squares to a Pauli. We provide a proof-of-principle numerical simulation that prepares a magic state using such logical Clifford measurements. Our work enables practical simulation of logical MSP protocols without resorting to approximations or resource-intensive state-vector simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗