Search NASA⌕ Search

SEARCH · Search NASA

Results for “numerical validation”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Quantifying and modeling of stress-driven short-circuits in lithium-ion batteries in electrified vehicles

Despite the huge expansion of electric vehicle sales in the market, customers are discouraged by the possible catastrophic consequences brought by the safety issues of lithium-ion batteries, such as internal short circuits, especially in crash scenarios. Herein, we reveal the quantitative relationship between the deformations of the battery and the internal short circuit. By in situ and ex situ observations and post mortem characterization of the component materials, we quantify the stress-driven internal short circuit and failure behavior of the component material. With the aid of a validated numerical computational model as well as the in situ characterization of the global-field temperature, we successfully identify the minor and major short circuits of the cells upon various mechanical abusive loadings. Finally, we establish the internal short circuit criteria for typical formats of batteries. Furthermore, this discovery also provides a fundamental understanding of both internal and external stress-driven short circuits in a much broader context.

25 ENERGY STORAGE↗
DOE OSTI2025

On the discretization error of the discrete generalized quantum master equation

The transfer tensor method (TTM) [Cerrillo and Cao, Phys. Rev. Lett. 112 , 110401 (2014)] can be considered a discrete-time formulation of the Nakajima–Zwanzig quantum master equation (NZ-QME) for modeling non-Markovian quantum dynamics. A recent paper [Makri, J. Chem. Theory Comput. 21 , 5037 (2025)] raised concerns regarding the consistency of the TTM discretization, particularly a spurious term at the initial time t = 0. Here, this work presents a detailed analysis of the discretization structure of the TTM, clarifying the origin of the initial-time correction and establishing a consistent relationship between the TTM discrete-time memory kernel K N and the continuous-time NZ-QME kernel $\mathscr{K}$( N Δ t ). This relationship is validated numerically using the spin-boson model, demonstrating convergence of reconstructed memory kernels and accurate dynamical evolution as Δ t → 0. While the TTM provides a consistent discretization, we note that alternative schemes are also viable, such as the midpoint derivative/midpoint integral scheme proposed in Makri’s work. The relative performance of various schemes for either computing accurate $\mathscr{K}$( N Δ t ) from exact dynamics or obtaining accurate dynamics from exact $\mathscr{K}$( N Δ t ) warrants further investigation.

Density-matrix↗
DOE OSTI2021

Steady-State Mechanical Analysis for Target Assembly in the Material Plasma Exposure eXperiment Facility

The Material Plasma Exposure eXperiment (MPEX) steady-state linear plasma facility is currently under design at Oak Ridge National Laboratory to expose target specimens to fusion divertor regimes. The neutron-irradiated target is actively cooled and remote handled in the MPEX facility for conducting plasma-material–interaction (PMI) experiments. Here, the steady-state stresses in the target and target assembly system are investigated using two-dimensional (2-D) and three-dimensional (3-D) models to provide expected stresses/strains under the heat loads to which various system components would be exposed during MPEX operation. The calculated temperatures from the 2-D axisymmetric mechanical model were found to be in excellent agreement with those from the full 3-D thermohydraulic model, providing a strong model validation. Numerical simulation results for the steady-state mechanical model indicate nonuniform distributions for the temperature, stress, and deformation within the critical components. For the initial design, the deformation results indicate possible gap openings between contacting surfaces below the plasma-facing materials. To reduce the possibility of interfacial gap opening, the target assembly was slightly changed and evaluated using the 2-D stress model. Numerical simulation results indicate that the interfacial gap openings can be minimized without drastically changing the entire target assembly. The stress-strain conditions for the target will be further used to assess the appropriate operation during MPEX experiments and gain insight into materials science phenomena during PMI.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Sweeping Horndeski canvas: new growth-rate parameterization for modified-gravity theories

We propose and numerically validate a new fitting formula that is sufficiently accurate to model the growth of structure in Horndeski theories of modified gravity for upcoming Stage IV and V large-scale structure surveys. Based on an analysis of more than 18,000 Horndeski models and adopting the popular parameterization of the growth rate f(z) = Ω M (z) γ , we generalize the constant growth index γ to a two-parameter redshift-dependent quantity, γ(z), that more accurately fits these models. We demonstrate that the functional form γ(z) = γ 0 + γ 0 z 2 /(1 + z) improves the median χ2 of the fit to viable Horndeski models by a factor of ~40 relative to that of a constant γ, and is sufficient to obtain unbiased results even for precise measurements expected in Stage IV and V surveys. Finally, we constrain the parameters of the new fitting formula using current cosmological data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Designing stellarators using perpendicular permanent magnets

We have developed a fast method to design perpendicular permanent magnets for simplifying stellarator coils based on existing codes. Coil complexity is one of the main challenges for stellarators. To date, only electromagnetic coils have been used to generate 3D fields for stellarators. Permanent magnets provide an alternative way to produce the desired magnetic field for optimized stellarators. In this paper, we revisit the concept of representing surface current using magnetic dipoles and carry out numerical validations. A surface magnetization is proven to be equivalent to the surface current that can be linearly solved by existing coil design codes. An incremental multi-layer method has been developed to obtain a practical solution that is attainable with present permanent magnets. With this method, we can reproduce a half-Tesla NCSX configuration using specially designed neodymium magnets together with simple planar coils. It shows that stellarator coils could be substantially simplified by adopting permanent magnets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Modeling of frequency-sweeping Alfvén modes in the TJ-II stellarator

Alfvénic activity has been observed in the TJ-II stellarator which resembles the frequency sweeping demonstrated by Alfvén cascade modes in tokamaks. A numerical validation study was conducted using a reduced magnetohydrodynamic (MHD) model to show that such modes could only have been observed in discharges where the rotational transform profile was non-monotonic. During experiments, coil current was varied which resulted in shifting of the minimum value of the rotational transform profile. To mimic this effect, we study the Alfvénic activity predicted by the reduced MHD model for a set of input rotational transform profiles with varying minima. A mode is found whose toroidal and poloidal mode numbers match those predicted in experiments which sweeps downward/upward in frequency as the minimum value of the rotational transform profile is increased/decreased. The results serve as a demonstration of the validity and utility of MHD spectroscopy.

Alfvén cascades↗
DOE OSTI2024

An international benchmark for wind plant wakes from the American WAKE ExperimeNt (AWAKEN)

This article introduces the first benchmark study within the International Energy Agency Wind Task 57 framework, focusing on wind plant wakes. Leveraging data from the American WAKE ExperimeNt (AWAKEN), the benchmark aims to assess the accuracy of simulation tools in modeling wind plant wakes and their impact on the downstream flow under diverse inflow conditions. The AWAKEN field campaign, conducted in Oklahoma from 2022 to 2024, provides unprecedented observations of wind plant-atmosphere interactions, thus offering a large dataset to validate numerical models of different complexity. The benchmark will include three phases—code calibration, blind comparison, and iteration—allowing participants to refine their numerical models based on the feedback from the benchmark team. This article describes the benchmark case study selected from observations providing details on atmospheric conditions, wake evidence, and wind turbine operation. The benchmark’s structure and timeline, along with the expected publication of results, are discussed as well. This collaborative effort aims to enhance the accuracy of wind plant wake simulations, thus contributing to the improvement of wind energy production estimates.

17 WIND ENERGY↗
DOE OSTI2026

Vapor-cavity-QED system for quantum computation and communication

In this work, we propose performing key operations in quantum computation and communication using room-temperature atoms moving across a grid of high-quality-factor, small-mode-volume cavities. These cavities enable high-cooperativity interactions with single atoms to be achieved with a characteristic timescale much shorter than the atomic transit time, allowing multiple coherent operations to take place. We study scenarios where we can drive a Raman transition to generate photons with specific temporal shapes and to absorb, and hence detect, single photons. The strong atom-cavity interaction can also be used to implement the atom-photon controlled-phase gate, which can then be used to construct photon-photon gates, create photonic cluster states, and perform nondemolition detection of single photons. We provide numerics validating our methods and discuss the implications of our results for several applications.

Austin, Sharoon [National Institute of Standards a↗
DOE OSTI2020

Framework for simulating gauge theories with dipolar spin systems

Gauge theories appear broadly in physics, ranging from the standard model of particle physics to long-wavelength descriptions of topological systems in condensed matter. However, systems with sign problems are largely inaccessible to classical computations and also beyond the current limitations of digital quantum hardware. In this work, we develop an analog approach to simulating gauge theories with an experimental setup that employs dipolar spins (molecules or Rydberg atoms). We consider molecules fixed in space and interacting through dipole-dipole interactions, avoiding the need for itinerant degrees of freedom. Each molecule represents either a site or gauge degree of freedom, and Gauss's law is preserved by a direct and programmatic tuning of positions and internal state energies. This approach can be regarded as a form of analog systems programming and charts a path forward for near-term quantum simulation. As a first step, we numerically validate this scheme in a small-system study of U(1) quantum link models in (1+1) dimensions with link spin S = 1/2 and S =1 and illustrate how dynamical phenomena such as string inversion and string breaking could be observed in near-term experiments. Our work brings together methods from atomic and molecular physics, condensed matter physics, high-energy physics, and quantum information science for the study of nonperturbative processes in gauge theories.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Factorization Machine Learning for Disaggregation of Transmission Load Profiles with High Penetration of Behind-the-Meter Solar

The ever-growing high penetration of ubiquitously distributed energy resources, especially behind-the-meter solar (BTM) generations, has significant impacts on nodal load (i.e., net injection) profiles and consequently caused imperative operational challenges to system operators such as regional transmission organizations (RTOs). Illustrated by real-world nodal data and examples at PJM Interconnection, this paper first discusses the application and necessity of effectively extracting daily nodal load profiles in a non-intrusive manner. More importantly, a novel bi-level architecture, including Factorization Machines (FM) learning procedure has been proposed to effectively disaggregate not only one node but every node in an RTO service territory. Specifically, FM leaning is adopted to capture the interconnections between related features to better utilize the correlation between buses in the same region and between a single bus and the zonal load. The proposed bi-level technique is numerically validated using real-world, minute-level, normalized, and anonymized nodal data at PJM service territory.

behind the meter solar, load disaggregation, load ↗
DOE OSTI2021

On Distributed Model-Free Reinforcement Learning Control with Stability Guarantee

Distributed learning can enable scalable and effective decision making in numerous complex cyber-physical systems such as smart transportation, robotics swarm, power systems, etc. However, the stability of the system is usually not guaranteed in most existing learning paradigms; and this limitation can hinder the wide deployment of machine learning in decision making of safety-critical systems. This paper presents a stability guaranteed distributed reinforcement learning (SGDRL) framework for interconnected linear subsystems, without knowing the subsystem models. While the learning process requires data from a peer-to-peer (p2p) communication architecture, the control implementation of each subsystem is only based on its local state. The stability of the interconnected subsystems will be ensured by a diagonally dominant eigenvalue condition, which will then be used in a model-free RL algorithm to learn the feedback control gains. The RL algorithm structure follows an off-policy iterative framework, with interleaved policy evaluation and policy update steps. We numerically validate our theoretical results by performing simulations on four interconnected sub-systems.

Mukherjee, Sayak↗
DOE OSTI2024

A Transactive Approach for Service Restoration Utilizing Customer Load Flexibility and Grid-Edge Resources

This paper develops a transactive energy system model to restore electricity to customers in an isolated distribution system after an outage. The model is developed to engage a variety of customer types - prosumers, flexible loads, critical/noncritical customers, and distributed generators - as active participants in the restoration process. Unlike many existing transactive approaches, the proposed model is developed for service restoration and accounts for various customer types and their autonomy and privacy through an iterative approach to determine the optimal market price, while maintaining systemlevel power flow and voltage constraints. The advantages of the proposed approach are numerically validated on a modified IEEE 123-bus test system.

distributed energy resources↗
DOE OSTI2020

Distributed Conditions for Small-signal Stability of Power Grids and Local Control Design

Operating modern power grids with stability guarantees is markedly important. Typical methods for analyzing and certifying power grid stability are largely centralized relying on the ability of the system operator to gather network-wide information and accurately compute the system's eigenvalues. These methods are oftentimes not privacy-preserving and computationally burdensome. They are therefore, not well-suited to modern power grids where small-signal stability has to be evaluated timely, efficiently and in a privacy-preserving fashion. Herein, we introduce a distributed methodology for certifying small-signal stability of power grids and designing the local controllers. First, we analytically derive distributed conditions for network-wide stability that bus agents can inspect using local information. By leveraging these conditions, we then introduce a distributed control design algorithm (DCDA) that can guide the local control design so that stability of the interconnected system is guaranteed. The agents that adopt the proposed distributed algorithm are responsible for tuning their local controllers, producing their local control commands and ensuring that their local stability condition is met. The system operator is only responsible for verifying network-wide stability upon receiving affirmative responses from all agents and, announcing, that the overall system is stable. The proposed DCDA algorithm is numerically validated via simulations using the IEEE 39-bus system.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Implementation of a Model Predictive Control Strategy to Regulate Temperature Inside Plug-Flow Solar Reactor With Countercurrent Flow

Abstract Solar-driven thermochemical energy storage systems are proven to be promising energy carriers (solar fuels) to utilize solar energy by using reactive solid-state pellets. However, the production of solar fuel requires a quasi-steady-state process temperature, which represents the main challenge due to the transient nature of solar power. In this work, an adaptive model predictive controller (MPC) is presented to regulate the temperature inside a tubular solar reactor to produce solid-state solar fuel for long-term thermal storage systems. The solar reactor system consists of a vertical tube heated circumferentially over a segment of its length by concentrated solar power, and the reactive pellets (MgMn2O4) are fed from the top end and flow downwards through the heated tube. A countercurrent flowing gas supplied from the lower end interacts with flowing pellets to reduce it thermochemically at a temperature range of 1000—1500 °C. A low-order physical model was developed to simulate the dynamics of the solar reactor including the reaction kinetics, and the proposed model was validated numerically by using a 7-kW electric furnace. The numerical model then was utilized to design the MPC controller, where the control system consists of an MPC code linked to an adaptive system identification code that updates system parameters online to ensure system robustness against external disturbances (sudden change in the flow inside the reactor), model mismatches, and uncertainty. The MPC controller parameters are tuned to enhance the system performance with minimum steady-state error and overshoot. The controller is tested to track different temperature ranges between 500 °C and 1400 °C with different particles/gas mass flowrates and ramping temperature profiles. Results show that the MPC controller successfully regulated the reactor temperature within ± 1 °C of its setpoint and maintained robust performance with minimum input effort when subjected to sudden changes in the amount of flowing media and the presence of chemical reaction.

Engineering↗
DOE OSTI2023

An Experimental Investigation of Heat Transfer for Supercritical Carbon Dioxide Cooling in a Staggered Pin Fin Array

Pin fin geometry in the trailing edge as turbulators in tur-bine blades has been a subject thoroughly researched for the past several decades. To improve the efficiency of turbomachinery, higher turbine inlet temperatures are necessary, introducing the need for internal cooling channel features, such as pin fin arrays. These solutions are necessary to avoid material failure at higher operating conditions by enhancing the heat transfer of the blade. Over the past decade, a shift toward applying supercritical carbon dioxide or sCO2 cycles has shown increased potential compared to the industry standard air cycles. Promising higher cycle ef¬ficiency, smaller footprints, and zero greenhouse gas emissions, sCO2 has become an attractive alternative to well-established air cycles. This alternative cycle introduces the need to re-visit the heat transfer capabilities of internal cooling geometries within this novel operating fluid environment. The first objective of this paper is to experimentally demonstrate the heat transfer charac-teristics of sCO2 as the secondary internal cooling fluid environ¬ment through testing a pin fin array in a closed loop at 200 bar and 400 Celsius inlet conditions, which stands well within the supercritical region of sCO2. Secondly, the experimental process aims to validate numerical simulation methodology for studying sCO2 environments for internal cooling of turbine blades in a matching pin fin geometry.

Wardell, Ryan J↗
DOE OSTI2023

Superdiffusion from Nonabelian Symmetries in Nearly Integrable Systems

The Heisenberg spin chain is a canonical integrable model. As such, it features stable ballistically propagating quasiparticles, but spin transport is subballistic at any nonzero temperature: An initially localized spin fluctuation spreads in time t to a width t 2/3 . This exponent as well as the functional form of the dynamical spin correlation function suggest that spin transport is in the Kardar–Parisi–Zhang (KPZ) universality class. However, the full counting statistics of magnetization is manifestly incompatible with KPZ scaling. A simple two-mode hydrodynamic description, derivable from microscopic principles, captures both the KPZ scaling of the correlation function and the coarse features of the full counting statistics, but remains to be numerically validated. These results generalize to any integrable spin chain invariant under a continuous nonabelian symmetry and are surprisingly robust against moderately strong integrability-breaking perturbations that respect the nonabelian symmetry.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Development of a three-dimensional APOLLO3 neutrons deterministic scheme for the CABRI reactor

CABRI is an experimental reactor to study the fuel behavior during reactivity injection transients. These transients being highly multiphysics, the development of suitable modeling and simulation tools to simulate them is important for the optimization of the tests and the control of the experimental conditions. This paper focuses on the development of an APOLLO3 deterministic core calculation tool dedicated to the CABRI transient analysis. It represents the first stage of the incremental process for the implementation of a multiphysics time-dependent modeling of the CABRI transient. The neutron calculation scheme is based on a classical two-step approach. The first step consists of a 281-energy group calculation flux with the TDT-MOC (Method Of Characteristics) solver for cross-section space and energy (23 groups) collapsing for the CABRI different assembly clusters. The bias on a 2D core neutron calculation due to the self-shielding calculation and collapsing on a restricted pattern are investigated thanks to a comparison with a direct full 2D calculation on a quarter of core. The second step relies on a pin-resolved transport 3D transport core calculation with the SN solver MINARET. A progressive numerical validation process is followed to quantify the calculation biases on reactivity and reaction rates at each step using reference calculations with the stochastic code TRIPOLI4. The next development stage toward a multiphysics scheme will be the implementation of the 3D-kinetics equation resolution and the coupling with a core thermal-hydraulics model. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

A Bayesian approach to nanoparticle characterization

We introduce and numerically validate a Bayesian method for plasmonic nanometrology. Applicable to any system described by a scattering cross section, our approach quantifies uncertainty automatically and enables model comparison through Bayes factors.

Lukens, Joseph↗