Search NASA⌕ Search

SEARCH · Search NASA

Results for “Non linear”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Measure this, not that: Optimizing the cost and model-based information content of measurements

Model-based design of experiments (MBDoE) is a powerful framework for selecting and calibrating science-based mathematical models from data. Here, this work extends popular MBDoE workflows by proposing a convex mixed integer (non)linear programming (MINLP) to optimize the selection of measurements. The solver MindtPy is modified to support calculating the D-optimality objective and its gradient via an external package, scipy, using the grey-box module in Pyomo. The new approach is demonstrated in two case studies: estimating highly correlated kinetics from a batch reactor and estimating transport parameters in a large-scale rotary packed bed for CO 2 capture. Both case studies show how examining the Pareto optimal trade-offs between information content measured by A- and D-optimality versus measurement budget offers practical guidance for selecting measurements for scientific experiments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A high-order explicit Runge-Kutta approximation technique for the shallow water equations

Here, we introduce a high-order space–time approximation of the Shallow Water Equations with sources that is invariant-domain preserving (IDP), well-balanced with respect to rest states, and employs a novel explicit Runge–Kutta (ERK) introduced in Ern and Guermond (SIAM J. Sci. Comput. 44(5), A3366–A3392, 2022) for systems of non-linear conservation equations. The resulting method is then numerically illustrated through verification and validation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Novel artificial neural network model for instantaneous power losses and operational efficiency mapping of MW-scale vanadium redox flow battery for improved technoeconomic analysis

A novel data-driven, machine-learning-based method for modeling the instantaneous power losses of a distribution-sited 2 MW/8MWh vanadium redox flow battery (VRFB), a grid-scale electrochemical storage technology, is introduced and compared against benchmark empirical modeling approaches, including symmetric and asymmetric models, as well as a recent convex hull modeling approach. The novel loss modeling method introduces several advantages over the benchmark models and over simplistic efficiency estimates, the most significant of which is that the model can accurately reflect the stepwise and non-linear parasitic losses associated with the duty cycles of mechanical auxiliary systems like pump motor drives and blower fans. Residuals of the models are compared; the proposed data driven model features significantly improved accuracy over the benchmark models. The model's coefficient of determination is also improved relative to that of the benchmark models. Furthermore, a novel method for visualization of operational efficiency of the grid-scale storage technology is introduced. To demonstrate the benefits of the novel data-driven method for modeling the VRFB, the benchmark models and the proposed models are embedded into an Open DSS distribution network model to study two applications of the grid-scale electrical storage system: load leveling for grid support and energy arbitrage. This article demonstrates that the accuracy of the instantaneous power loss model significantly impacts the understanding of the state of charge of the VRFB. In turn, the accuracy of the efficiency modeling of the VRFB impacts the understanding of the potential economic value and technical benefits to the distribution network operators. In conclusion, the presented power loss modeling approach is, therefore, highly relevant for utility-stakeholders, battery asset owners, system engineers, system designers, and financial planners interested in evaluating or optimizing the operation of grid-scale VRFBs.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Numerical simulation of compressible fluid-dynamics in the chamber of inertial fusion energy systems

Here, this paper aims to establish new and innovative modeling capabilities for analyzing chambers in Inertial Fusion Energy (IFE) systems. IFE is emerging as a promising method to achieve fusion power production, but several challenges must be overcome to develop an IFE pilot plant or deploy commercial IFE systems. These challenges are both theoretical and technical, encompassing a deeper understanding of the underlying physical phenomena and the development of new technologies and materials. One of the needs is to develop mathematical models to describe IFE systems and numerical tools to simulate them. This paper contributes to this endeavor by presenting a new OpenFOAM solver for IFE systems, focusing on gas dynamics in their chambers. The analysis and development of chamber designs will play a significant role in the transition from single-shot experiments to high-repetition rates, as there is a need to protect the chamber walls from the intense radiation fields produced by fusion reactions. A promising design option, normally referred to as thick wall chamber design, consists in using lithium or molten salt jet arrays within the chamber. A critical phenomenon is the venting of high-pressure gases from the center to the external part of the chamber, passing through the blanket array. This process involves the propagation and attenuation of strong pressure waves, requiring suitable modeling approaches for compressible fluid-dynamics. The solver proposed in this work implements a multi-material hydrodynamics model tailored to accurately describe the non-linear propagation of pressure waves while avoiding numerical oscillation issues typical of high-velocity compressible simulation. This solver is verified against numerical test cases, validated against experimental data, and applied to the analysis of the High-Yield Lithium-Injection Fusion-Energy (HYLIFE-I) concept. The relevance of this paper is threefold. Firstly, it contributes to developing and testing modeling approaches for compressible fluid-dynamics phenomena, with specific focus on the new and unexplored topic of IFE thick-liquid-wall blanket modeling. Secondly, it marks one of the first applications of the OpenFOAM library in the research field of IFE systems. Finally, the investigated problem is of practical interest for IFE developers, as it provides useful indications about relevant phenomena in pressure wave propagation in the chamber of these systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

An in-situ view cell system for investigating swelling behavior of elastomers upon high-pressure hydrogen exposure

The transition to hydrogen as a clean and efficient energy carrier is impeded by challenges in the compatibility of hydrogen with materials used within hydrogen infrastructure. Elastomers, crucial in sealing components, often exhibit premature failures in high-pressure hydrogen environments due to excessive swelling. This study employs an innovative in-situ view cell system to assess the swelling behavior of hydrogenated nitrile butadiene rubber (HNBR) under various hydrogen conditions. The system, designed to withstand pressures up to 96.5 MPa, incorporates Digital Image Correlation (DIC) for strain measurements and volume estimation. Results reveal non-linear volume increases during depressurization, challenging conventional assumptions. Furthermore, investigations into peak hydrogen pressures and pressure-holding scenarios during decompression highlight complex swelling trends. The introduction of a novel computer vision (CV) method enhances precision in volume estimation, overcoming DIC limitations. The study provides insights into mitigating elastomer swelling, crucial for developing robust materials to support future hydrogen-driven energy systems.

Elastomer↗
DOE OSTI2025

Droplet breakup by multimodal nonlinear Rayleigh Taylor instability

A droplet impacted by a shock wave will undergo a process of fragmentation due to the development of interfacial hydrodynamic instabilities. The interface experiences variable acceleration and shear that result in the development of both inertial (Rayleigh–Taylor) and shear (Kelvin–Helmholtz) instabilities. These perturbations grow in time and drive the fragmentation and breakup of the deformed droplet. Experiments are performed on nominally 0.86 mm water droplet subjected to a Mach 7.6 detonation wave, resulting in a high Weber number (~36,000) breakup event. Perturbation growth is measured from a series of high-speed (> MHz) shadowgraph images. It is proposed that, given the size of the large-scale perturbations observed in experiments, these instabilities are growing in the non-linear regime and can be described by bubble-merger models for nonlinear mixing. Calculations are performed for the growth rates and size of these instabilities using deformation and external flow models to establish the time-dependent boundary conditions. The concurrence of the measured perturbation widths and the predictions of the simple model lend credence to the theory. This novel approach serves to open a new avenue in the characterization of droplet breakup via hydrodynamic instabilities.

Droplet breakup↗
DOE OSTI2025

High-pressure Cr 3+ luminescence and Raman spectroscopy of a natural MgAl 2 O 4 spinel to ∼60 GPa

Cr 3+ luminescence and Raman spectroscopy of a natural spinel were collected up to ∼60 GPa at ambient temperature using a Ne or a 4:1 methanol-ethanol mixture as the pressure transmitting media. Here, in the Ne experiments above ∼40 GPa, and in the 4:1 methanol–ethanol experiments above ∼20 GPa, all luminescence emission lines shift non-linearly, and by 59 GPa the spectra are broadened and are significantly changed. The Raman modes shift smoothly to higher frequencies with pressure, with a slight change in slope above 40 GPa, with no kinks or discontinuities. By ∼60 GPa, only the most intense Eg mode remains resolvable, appearing broadened and reduced in intensity. The change in slope in the Raman modes and luminescence bands may be due to a change in the compression mechanism of spinel. When spinel was decompressed from 40 GPa all changes were fully reversible; however, when the samples were decompressed from 54 to 60 GPa in Ne and 59 GPa in 4:1 methanol-ethanol the spectra were irreversibly changed. The change in the luminescence spectra and the Raman spectra are likely due to a pressure induced disorder of the anions or partial amorphization between 40 and 60 GPa.

Cr3+↗
DOE OSTI2026

Multi-material additive manufacturing of aluminum 6061-T6 alloy with stainless steel 304: Suppression of intermetallic compounds and interface growth mechanism

Fabricating integrated aluminum-steel structures is difficult because of metallurgical incompatibilities that promote intermetallic compound (IMC) formation at their interfaces. This work introduces a practical route for additively manufacturing 6061-T6 aluminum alloy (AA6061-T6) directly onto 304 stainless steel (SS304) using a high shear strain rate-assisted interfacial bonding mechanism. Through the friction extrusion deposition-based additive manufacturing method, both single- and multi-layer deposits of AA6061-T6 were successfully fabricated on SS304, yielding a robust hybrid multi-material structure. Comprehensive analyses of deposition quality, interface porosity, and bonding performance showed that the interface achieved a tensile strength exceeding 154 MPa under quasi-static loading. Microstructural observations revealed that the deposited aluminum layers experienced severe plastic deformation, resulting in pronounced grain refinement. Importantly, the interfacial zone was found to be free of brittle IMCs, instead containing a thin amorphous layer that evolved through a non-linear reaction-zone growth law. This solid-state additive strategy establishes a promising pathway for lightweight structural systems, nuclear energy component cladding, and multifunctional engineering components, redefining how dissimilar metals can be integrated for advanced applications.

Aluminum alloy↗
DOE OSTI2025

Direct sensitivity analysis on the parameterization of crystal plasticity models

Various methods for calibrating crystal plasticity finite element (CPFE) models lead to non-unique input parameter values, which subsequently introduce uncertainty in the predicted mechanical response. Sensitivity analysis (SA) conducted on crystal plasticity models is used to identify how variability in these parameters contribute to output uncertainty. Traditional SA on CPFE parameters uses simplified surrogate models to save computational time. However, the accuracy of the surrogate models depends on the quantity of training data used, and any modeling error can propagate into the SA results, potentially affecting their reliability. In this work, the elementary effects test (EET) method, a global SA technique using direct CPFE simulations was employed, and the results obtained were compared with the First Order Second Moment (FOSM) method. ExaConstit, an open-source GPU-enabled CPFE code, was used to perform the simulations and direct SA. The EET method was accurately able to capture the non-linear effects of all the input parameters on the output and is a valuable approach for reliably attributing parameter sensitivities in CPFE models. Based on the results, efficient strategies to perform future parameter calibration and SA are discussed. Additionally, the SA trends observed in different single crystal orientations closely mirrored the activity of the slip systems.

Elementary Effects Test↗
DOE OSTI2025

Sensitivity study of hydrogen Soret transport in yttrium Hydride-Based nuclear fuel

Yttrium hydride is an excellent solid neutron moderator material for high temperature nuclear reactor applications due to its high hydrogen density and exceptional hydride stability at high temperatures. Despite these attractive characteristics, the details of how hydrogen behaves within yttrium hydride while temperature gradients exist are still not well understood. The evolution of the hydrogen composition profile resulting from a temperature gradient requires knowledge of hydrogen’s heat of transport, a critical parameter that has not yet been measured for this material. In this work, we perform hydride redistribution, hydrogen dissociation, and hydrogen leakage calculations while varying the Soret heat of transport of hydrogen in yttrium hydride to elucidate the sensitivity of hydride stability under temperature gradients to this parameter. This study analyzes hydride stability of a hypothetical uranium-yttrium hydride nuclear fuel design during operation of a high temperature liquid metal-cooled nuclear reactor. Assuming U-YH x could be fabricated in a physically stabilized manner, this fuel system can likely maintain hydride stability while operating at very high power densities and temperatures. We find that even though the hydrogen dissociation pressure in the gas gap does vary by several percent as the heat of transport temperature parameter is varied, the hydrogen content in the U-YH x fuel meat is relatively insensitive to this parameter over the course of a high burnup fuel cycle; this is due to yttrium hydride’s excellent hydrogen retention under the high temperature conditions considered here. Here, this suggests that hydride stability analyses are insensitive to the value of the Soret heat of transport in U-YH x under steady state liquid metal-cooled reactor conditions. However, the susceptibility to internal gas overpressurization-induced stress-rupture of the cladding during a high temperature transient is more sensitive to this parameter due to the non-linear dependence of hydrogen gas dissociation pressure vs. composition and temperature.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

A fast and robust computational modeling approach for density and shape predictions in powder metallurgy hot isostatic pressing

Powder metallurgy hot isostatic pressing (PM-HIP) is an advanced manufacturing process that produces near-net-shape parts with high material utilization and uniform microstructures. PM-HIP is frequently used for producing small-scale parts with complicated geometries and is potentially economical for producing large-scale parts. However, excessive post-HIP shape distortions can reduce its effectiveness and economic advantage, especially for larger parts. A PM-HIP computational model can predict and help mitigate these distortions. However, due to complex deformation mechanisms and thermo-mechanical coupling present in PM-HIP processes, these non-linear computational models sometimes become numerically unstable. The numerical instabilities in these models can lead to very slow convergence or no convergence at all, which often translates to slow and unreliable models. These limitations are more pronounced in large models with complicated geometries. Hence, in this work, an alternative modeling approach is presented that improves numerical stability and computational performance. The presented approach achieves these improvements through approximating the fully coupled thermo-mechanical PM-HIP model as a decoupled model and adding inertial damping to the model’s mechanical part. In conclusion, a comparison with the fully coupled model indicated a slight dip in prediction accuracy (<5% error) but significant improvements in numerical stability (>20 times larger time step size) and computational performance (5-10 times speed-up with less computational resource usage) when using the presented approach.

Hot isostatic pressing↗
DOE OSTI2026

Oil price states and drivers: An analysis of the second-month spot-futures price differential

Oil remains a dominant component of global energy use, and its price, characterized by frequent changes and an ever-present potential for large swings, continues to be a focus of industry participants, policymakers and analysts attention. Here, this study examines the behavior of future spot oil prices using a continuous-time hidden Markov model (HMM) and daily price data spanning years 2007 to 2024. We identify six states in the second-month WTI spot-futures price differential and assess the roles of eleven futures price, quantity, financial market, and geopolitical/volatility variables in each state. The model yields several insights into the workings of the oil market and the relative roles of these drivers. We find support for several theoretical and empirical findings in the oil market literature, including the role of inventory, volatility/risk, and market responses to contango/backwardation in futures markets. A novel finding is that “normal contango” conditions represent a significant portion of next-day states in our in-sample data. Under the most volatile normal contango state, many of the oil market drivers differ markedly in coefficient signs and magnitudes from those in other states. The resulting model also performed well out-of-sample and would, in addition to understanding the impact of market drivers, be useful for short-term forecasting. Overall, the findings highlight the highly non-linear, regime-dependent interactions of the oil price and its drivers, emphasizing the importance of detailed information to market stakeholders.

Oladosu, Gbadebo A. [Oak Ridge National Laboratory↗
DOE OSTI2026

Extreme-scale EV charging infrastructure planning for last-mile delivery using high-performance parallel computing

Here, this paper addresses stochastic charger location and allocation problems under queue congestion for last-mile delivery using electric vehicles (EVs). The objective is to decide where to open charging stations and how many chargers of each type to install, subject to budgetary and waiting-time constraints. We formulate the problem as a mixed-integer non-linear program, where each station-charger pair is modeled as a multiserver queue with stochastic arrivals and service times to capture the notion of waiting in fleet operations. The model is extremely large, with billions of variables and constraints for a typical metropolitan area; even loading the model in solver memory is difficult, let alone solving it. To address this challenge, we develop a Lagrangian-based dual decomposition framework that decomposes the problem by station and leverages parallelization on high-performance computing systems, where the subproblems are solved by using a cutting plane method and their solutions are collected at the master level. We also develop a three-step rounding heuristic to transform the fractional subproblem solutions into feasible integral solutions. Computational experiments on data from the Chicago metropolitan area with hundreds of thousands of households and thousands of candidate stations show that our approach produces high-quality solutions in cases where existing exact methods cannot even load the model in memory. We also analyze various policy scenarios, demonstrating that combining existing depots with newly built stations under multiagency collaboration substantially reduces costs and congestion. These findings offer a scalable and efficient framework for developing sustainable large-scale EV charging networks.

Capacity allocation↗
DOE OSTI2025

Detecting Magnon-phonon coupling in yttrium iron garnet with variable temperature STEM-EELS

Magnons, quanta of spin wave excitations in magnetically ordered materials, have been identified as candidates for several potentially transformative technologies in recent years. Macroscopic techniques, such as neutron scattering or Raman spectroscopy, can be used to identify and analyze magnons, but provide relatively delocalized information about the sample. Understanding how the bonding and local structure of a material interacts with, and influences, the magnon population in a material is a crucial step toward the ability to produce any real-world application utilizing magnons. By leveraging the combined spatial resolution of scanning transmission electron microscopy (STEM) and the energy resolution of monochromated electron energy-loss spectroscopy (EELS) nanoscale analysis of magnons can be performed. While the weak interaction of magnons with the electron beam makes magnon EELS challenging on reasonable timescales, magnon-phonon coupling can be leveraged to understand magnons through their effect on the more easily measured phonons. Here, we examine yttrium iron garnet (YIG) flakes, and demonstrate non-linear, temperature-dependent shifts in the phonon frequencies, consistent with previously described magnon-phonon coupling effects. The ability to measure the temperature-dependence of vibrational frequencies with high precision in individual nanoscale flakes, demonstrates the ability to study magnon-phonon coupling in the STEM with unprecedented spatial resolution.

Reifsnyder, Alexander [ORNL]↗
DOE OSTI2025

A Bayesian framework to investigate radiation reaction in strong fields

Recent experiments aiming to measure phenomena predicted by strong-field quantum electrodynamics (SFQED) have done so by colliding relativistic electron beams and high-power lasers. In such experiments, measurements of collision parameters are not always feasible. However, precise knowledge of these parameters is required to accurately test SFQED. Here, we present a novel Bayesian inference procedure that infers collision parameters that could not be measured on-shot. This procedure is applicable to all-optical non-linear Compton scattering experiments investigating radiation reaction. The framework allows multiple diagnostics to be combined self-consistently and facilitates the inclusion of known information pertaining to the collision parameters. Using this Bayesian analysis, the relative validity of the classical, quantum-continuous and quantum-stochastic models of radiation reaction was compared for several test cases, which demonstrates the accuracy and model selection capability of the framework and highlight its robustness if the experimental values of fixed parameters differ from their values in the models.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Antielectrophoretic Response-Driven Bending–Tilting Deformation of Cationic Polyelectrolyte Brushes Drives Nonlinear Electroosmotic Transport in Brush-Grafted Nanochannels

In this paper, we use all-atom molecular dynamics (MD) simulations to describe a non-linearly enhanced electroosmotic (EOS) flow, where, in a nanochannel grafted with cationic PMETAC ([Poly(2-(Methacryloyloxy)Ethyl) Trimethylammonium Chloride]) brushes, a two-fold increase in the electric field strength leads to a several-fold (more than two-fold) increase in the EOS flow strength and volume flow rate. The electric field enforces the PMETAC brushes to undergo a bending-tilting driven deformation with a significant portion of the brush layer becoming parallel to the grafting surface. In response, a substantial fraction of the counterions leave the brush layer (hence become more mobile), but instead of going into the bulk, accumulate at the brush-bulk interface, i.e., stay in proximity of the brush segments aligned parallel to the grafting surface. This creates an interesting situation, where the counterions are not completely within the brush layer, yet they fully screen the brush charges. Such “freer” conditions enable the counterions to achieve very high velocity, thereby ensuring that the water solvating the counterions themselves move very fast triggering the significantly augmented EOS transport. Probing deeper we can identify that the bending-tilting driven brush deformation, enforcing the brushes to align parallel to the substrate, results from the anti-electrophoretic behavior of the brushes, where despite being positively charged, the brushes move against the electric field direction. Such an anti-electrophoretic behavior of the PE brushes, which has not been reported before, can be associated with the very fast velocities of the negatively charged counterions and the electrostatic and hydrodynamic coupling of the counterions with the positive functional groups of the brushes. Here, we anticipate that the findings of this paper will shed light on strategies for nanochannel flow, the anti-electrophoretic response of charged polymer chains, and the significance of capturing the detailed chemical architecture of polyelectrolytes in nanoscale science and engineering.

36 MATERIALS SCIENCE↗
DOE OSTI2024

On the Interplay Between Distributed Bulk Plasticity and Local Fault Slip in Evolving Fault Zone Complexity

Abstract We numerically investigate the role of plastic strain accumulation on the mechanical response of a planar strike‐slip fault. Our models show that fault‐zone strength significantly impacts the ensuing sequence of earthquakes. Weaker fault zones accumulating more plastic strain promote more complexity in the seismicity pattern through aperiodic earthquake occurrences and intermittent episodes of rupture and arrest. However, if the fault zone strength is high enough, the overall earthquake sequence is characterized by periodic fault‐spanning events. We find that both the fault normal stress and the fault geometric profile evolve throughout the earthquake sequence, suggesting a self‐roughening mechanism. Despite the significant impact of plasticity on the fault response, the width of the plastically deforming region in the fault zone is small compared to the fault length. Our results suggest a rich behavior in dynamically evolving fault zones and support the need for further high‐resolution studies of the highly non‐linear near‐fault region.

Geology↗
DOE OSTI2024

Shallow Soil Response to a Buried Chemical Explosion With Geophones and Distributed Acoustic Sensing

Abstract Shallow sediments can respond non‐linearly to large dynamic strains and undergo a subsequent healing phase as the material gradually recovers following the passing of seismic waves. This study focuses on the physical changes in the subsurface caused by the shaking from a buried chemical explosion detonated in a borehole in Nevada, USA, as a part of the Source Physics Experiment Phase II. The explosion damaged the shallow subsurface and modified the frequency content recorded by 491 geophones and 2240 Distributed Acoustic Sensing (DAS) channels within 2.5 km from surface ground zero. We observe a gradual shift of resonance frequencies in the 10–25 Hz frequency band in the hours following the explosion and develop a method to characterize the related logarithm‐type healing process of the shallow (i.e., upper ∼25 m) subsurface. We find that stronger levels of ground motion increase the relative degree of damage and duration of the subsurface healing; with the spall region exhibiting the largest degree of damage and longest healing recovery time. We observe coherent spatial patterns of damage with the region located to the southeast of the explosion exhibiting more damage than the southwest region. This study demonstrates that both DAS and co‐located geophones capture similar temporal changes associated with the physical processes occurring in the subsurface, with the high‐density sampling of DAS measurements enabling a new capability to monitor the fine‐scale changes of the Earth's shallow subsurface following the detonation of a buried explosion.

Viens, Loïc↗