Search NASA⌕ Search

SEARCH · Search NASA

Results for “non-linear systems”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
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

Preserving nonlinear constraints in variational flow filtering data assimilation

Data assimilation aims to estimate the states of a dynamical system by optimally combining sparse and noisy observations of the physical system with uncertain forecasts produced by a computational model. The states of many dynamical systems of interest obey nonlinear physical constraints, and the corresponding dynamics is confined to a certain sub-manifold of the state space. Standard data assimilation techniques applied to such systems yield posterior states lying outside the manifold, violating the physical constraints. This work focuses on particle flow filters which use stochastic differential equations to evolve state samples from a prior distribution to samples from an observation-informed posterior distribution. The variational Fokker-Planck (VFP)—a generic particle flow filtering framework—is extended to incorporate non-linear, equality state constraints in the analysis. To this end, two algorithmic approaches that modify the VFP stochastic differential equation are discussed: (i) VFPSTAB, to inexactly preserve constraints with the addition of a stabilizing drift term, and (ii) VFPDAE, to exactly preserve constraints by treating the VFP dynamics as a stochastic differential-algebraic equation (SDAE). Additionally, an implicit-explicit time integrator is developed to evolve the VFPDAE dynamics. The strength of the proposed approach for constraint preservation in data assimilation is demonstrated on three test problems: the double pendulum, Korteweg-de-Vries, and the incompressible Navier-Stokes equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Complex Electrical Conductivity of a Single‐Fractured Rock: Fracture‐and‐Matrix Coupling Mechanism and Aperture Size Predictions

Fractured rocks play a crucial role in myriad natural and engineered systems, including Earth's critical zone, oil/gas/geothermal reservoirs, and geological CO 2 /H 2 /waste storage systems. While complex electrical conductivity is extensively used to estimate the pore and grain sizes of conventional porous rocks and soils, it is rarely used to predict the aperture size of fractured rocks and this remains poorly understood. Here, integrating theory, simulations, and experiments, we show that under external fields, fractured rocks follow the fracture-and-matrix coupling to make the bulk complex conductivity non-linear with respect to water conductivity. We find that the relaxation time and quadrature conductivity for porous media do not apply to fractured rocks, but, instead, reasonably accurate predictions of aperture size can be made based on the true formation factor. This study unravels the fundamental mechanism governing conduction and polarization of fractured rocks and paves the way for the non-invasive investigation of global fractured rocks.

58 GEOSCIENCES↗
DOE OSTI2026

Understanding spectral dependance of laser-induced damage precursors in dielectric materials (Abbreviated report_23-ERD-006)

High peak and average power laser systems are typically limited by the handling fluence of the optical components. In particular, the multilayer dielectric coatings are known to be much lower operational fluence than the more ideal bulk materials. In this project we proposed and succeeded in probing different established classes of damage prone precursors at different wavelengths to understand their fundamental laser damage response as a function of wavelength. This study helped to shed light on the fundamental physics of the non-linear precursors that govern laser damage phenomena for ns-regime pulsed laser damage. In this study, we utilized the onsite coating capabilities (VPL, IBS coating lab) to purposefully generate laser damage-prone precursors in hafnia-based coatings (both single and multi-layer coatings). Specifically, we engineered coatings with craze lines initiated by nodules, generated coatings with our xenon-based coating process to suppress nanobubble formation and generated hafnia coatings under controlled oxygen flow conditions to study hafnia sub-oxides and oxygen flow dependance. In these studies we found that craze lines are rife with precursors that are sensitive to ultra-violet light but not to infrared light; we found that the removal of nanobubbles helps with all wavelengths tested, but is most impactful for ultra-violet light; we also found that ultra-violet laser damage performance of hafnia is closely matched to oxygen flow rate, while the infrared performance may be slightly better at lower flow rates. During this LDRD we also successfully stood up a new laser damage testing capability, namely a wavelength agile damage test station to study the spectral response of known laser-induced damage precursors. This is a unique and important capability for Lawrence Livermore National Lab, allowing us to understand the spectral response of materials under high intensity irradiation and damage. This is a vital tool to understand non-linear optical response at wavelengths that we have previously been unable to test at.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Third integer resonant extraction transit time simulation studies

In this work, we present the investigation of transit time of particles in the non-linear third-integer resonant extraction process. Transit time is defined as the number of turns a particle takes to get extracted once it is in the unstable region in the phase space, i.e., outside the triangular separatrix in case of third-integer resonance. The study of transit time is important because transit time directly contributes to the beam response time during resonant extraction and thus knowing it apriori would be practically useful in designing of the extraction system. In this work, we shall investigate the analytical derivation of the transit time of particles (to the first order Kobayashi Hamiltonian) in different parts of the phase space distribution and compare against the analytical results. We also compare the simulation result of the transit time of particles (with higher statistics) for the static as well as dynamic extraction conditions cases, particularly in the context of resonant extraction parameters for Mu2e experiment at Fermilab.

Narayanan, Aakaash [Fermilab]↗
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 OSTI2026

Aerosol-deep convection interaction based on joint cell-thermal tracking in Large Eddy Simulations during the TRACER campaign

In cumulus clouds, aerosol concentrations control cloud droplet concentrations, modifying cloud radiative properties, precipitation processes, and cloud electrification. However, mechanisms of aerosol-deep convection interactions are not well understood due to complex cloud dynamics and microphysics. We investigate the interaction of aerosols with isolated deep convection using Large Eddy Simulations of two cases during the TRacking Aerosol Convection interactions ExpeRiment (TRACER) near Houston, Texas, using a joint cell-thermal tracking algorithm. Cumulus thermals are droplet generators, since supersaturation and droplet nucleation coincide with thermal centers, where the strongest updrafts occur. Primary ice crystal formation does not take place inside thermals, but at layers where previous thermals detrained moisture. As subsequent thermals containing supercooled droplets penetrate these layers, hail and graupel form at or near these thermals. Higher aerosol concentrations result in higher droplet concentrations that suppress drizzle, delay warm rain processes, and transport more moisture aloft. This increases snow and ice amount, as well as graupel and hail, leading to more lightning. Polluted thermals initiate at slightly higher altitudes, and are slightly larger and faster, suggesting a weak invigoration. We also find more thermals per cell, but fewer isolated cells, since convection is more aggregated and intense, especially near the end of the 24 h simulation. Non-linear mesoscale feedback likely triggered by temperature and moisture responses to aerosol-thermal interactions causes the aggregation. Time-lagged aerosol-reinitialization experiments show that the mesoscale response is the predominant forcing for the invigoration. These changes happen within one day, on a smaller scale than previously suggested.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Adaptive Control for Load-Following of Boiling Water Reactors Part I: Linear Systems and Fully-Observable Dynamics

Automation control is a key strategy to improve the economic competitiveness of nuclear power plants. Not only does it help reduce operational costs, but it also extends the value proposition of these plants to nontraditional markets, including unattended operations in remote villages and space. However, the dynamics of the operating environments of nuclear reactors are subject to changes over time, and there are no widely adopted methods to ensure that the automation strategy will remain effective over the extended durations required for these applications. Adaptive control is a discipline that offers the possibility to accommodate such changes online. However, it relies on mathematical assumptions that must be respected to ensure robustness and reliability. In this work, we derive an adaptive control formulation for linear systems in which all states are observable and apply it to an instance of load-follow operation for Boiling Water Reactors. We assumed uncertainty in two factors: the temperature coefficient and the control rod worth, both of which are affected over time by the evolution of the nuclear reactor core environment. With an arbitrary penalty factor of 5, we found that the mean absolute and integral time absolute errors can be reduced by more than 90%, underscoring the strength of adaptive control. To extend the application to more challenges, different uncertainties and load-follow trajectories, as well as new formulations that include non-linearity and partial observability, are currently being developed.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

High-frequency Data Integration for Landscape Model Calibration of Carbon Fluxes Across Diverse Tidal Marshes

Terrestrial Aquatic Interfaces (TAIs), and tidal wetlands in particular, store large amounts of carbon yet are not well represented in Earth System Models (ESMs). Predictions of carbon cycling and greenhouse gas (GHG) emissions in tidal wetlands are highly uncertain. Eddy covariance (EC) towers provide ecosystem-scale GHG flux data at a temporal resolution (every 30min) that is helpful for parameterizing and improving mechanistic realism in ESMs. We propose to use a network of eddy covariance towers and standardized ancillary data streams, along with mesocosm experiments and statistical analyses, across diverse tidal wetlands of North America to develop and improve biogeochemical modeling at the TAI. Our overarching objective is to improve understanding and process-based modeling of gross primary productivity (GPP) and CH4 emission responses, both non-linear and asynchronous, to stressors including plant inundation, disturbance, salinity and nitrogen loading.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

The Energy Exascale Earth System Model Version 3: 2. Overview of the Coupled System

The Energy Exascale Earth System Model version 3 (E3SMv3) represents the latest advancement in Earth system modeling developed by the U.S. Department of Energy (DOE). Building upon previous versions, E3SMv3 introduces significant updates across its coupled components to enhance capability and improve fidelity. The atmosphere component incorporates advancements in chemistry, aerosol-cloud interactions, convection, and microphysics. The ocean features a new time-stepping scheme and a higher-resolution unstructured mesh with sub-ice-shelf cavities, while the sea ice model integrates advanced snow and ice physics for more realistic cryospheric simulations. The land model introduces prognostic vegetation dynamics and a new sub-grid topographic treatment of solar radiation. A new tri-grid configuration harmonizes the horizontal grids of the land and river components for improved process coupling. It is enabled by a new non-linear remapping between the atmosphere and land. E3SMv3 underwent extensive testing through a comprehensive simulation campaign, including pre-industrial control, idealized CO 2 experiments, and historical simulations spanning 1850–2024. The model demonstrates significant improvements in simulating the evolution of the historical surface temperature, particularly addressing the “pothole cooling” bias in earlier versions. Reduced aerosol-related forcing contributes to more realistic radiative forcing and better alignment with the observational record. Ocean heat content (OHC) and sea ice trends are also improved as a result.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Evolution of the Antarctic Ice Sheet from 2000–2300 and beyond: model sensitivity and uncertainty analysis using MPAS-Albany Land Ice

We present a description of the Antarctic Ice Sheet model configuration submitted to the ISMIP6-Antarctica-2300 experiment using the MPAS-Albany Land Ice model, along with three new sets of simulations: (1) a set of extended simulations to 2500 for three forced experiments and to 2775 for the control experiment; (2) a sensitivity analysis of our model configuration to parameters controlling basal sliding and sub-shelf melt, and to model structural choices including the choice of the energy and stress balances; and (3) a 72-member ensemble run on graphics processing units (GPUs) and analysis of variance to determine the primary sources of uncertainty in our ice-sheet model projections. Our extended simulations predict rapid retreat beginning after 2300 for SSP1-2.6 forcing and after 2500 for present-day (control) forcing, primarily in the Amundsen Sea Embayment. We find that varying the sub-shelf melt parameter between the 5th to 95th percentile values for a mean-Antarctic calibration target results in an up to ∼ ± 40 % change in sea-level contribution relative to our baseline simulations that used the median value. Using a linear basal sliding law reduces sea-level contribution by 51 %–73 % relative to our baseline nonlinear sliding law with an exponent of 1/5. When using basal sliding law exponents of 1/3 and 1/10, the overall difference from our baseline simulations at 2300 is on the order of 10 %. The Amundsen Sea Embayment region displays a strongly non-linear dependence of mass loss on the sliding law exponent, with no discernible relationship between the sliding law exponent and the mass loss by 2300, while the sectors feeding the Ross and Filchner-Ronne ice shelves exhibit more mass loss with a more-plastic sliding law. Our model fidelity sensitivity experiments reveal a 9 %–31 % increase in sea-level contribution when using a depth-integrated stress balance approximation relative to our three-dimensional solver, while using a fixed-in-time temperature field increases sea-level contribution by 14 %–88 % relative to two thermomechanically coupled configurations. Our 72-member ensemble and analysis of variance show that the uncertainty in long-term projections is dominated by the choice of Earth system model forcing and the presence or absence of hydrofracture forcing, rather than uncertainty in sliding and sub-shelf melt parameters.

58 GEOSCIENCES↗
DOE OSTI2025

Variational AutoEncoders Reveal Intensifying GPP Extremes in Continental United States based on CESM2 Simulations

Climate extremes significantly impact terrestrial carbon cycle dynamics, necessitating robust methods for detecting and analyzing anomalous behavior in plant productivity. This study presents a novel application of variational autoencoders (VAE) for identifying extreme events in gross primary productivity (GPP) from Community Earth System Model version 2 simulations across four AR6 regions in the Continental United States. We compare VAE-based anomaly detection with traditional singular spectral analysis (SSA) methods across three time periods: 1850-80, 1950-80, and 2050-80 under SSP5-8.5 scenario. The VAE architecture employs three dense layers and a latent space with input sequence length of 12 months, training on normalized GPP time series to reconstruct the GPP and identify anomalies based on reconstruction errors. Extreme events are defined using 5th percentile thresholds applied to both VAE and SSA anomalies. Results demonstrate strong regional agreement between VAE and SSA methods in spatial patterns of extreme event frequencies, despite VAE consistently producing higher threshold values (179-756 GgC for VAE vs. 100-784 GgC for SSA across regions and periods). Both methods reveal increasing magnitudes and frequencies of negative carbon cycle extremes toward 2050-80, particularly in Western and Central North America. The VAE approach shows comparable performance to established SSA techniques while offering computational advantages and enhanced capability for capturing non-linear temporal dependencies in carbon cycle variability. This research demonstrates the potential of deep learning approaches for extremes detection and provides a foundation for improved understanding of future carbon cycle risks under future conditions.

Sharma, Bharat [ORNL] (ORCID:0000000266982487)↗
DOE OSTI2026

Robust optimization of flexible diafiltration systems for critical mineral separations

This paper provides major contributions in expanding the literature for membrane process design with critical mineral recovery applications and showcasing the importance of robust design techniques for reducing risks of underperformance in such systems. Here, a membrane process flowsheet featuring PrOMMiS membrane models for recovering lithium/cobalt from spent batteries is showcased, uncertainty in membrane sieving and localized fouling are considered, and robust designs are obtained using the PyROS toolset. This paper is intended for a general audience of researchers working in critical minerals, membranes, and optimization related areas.

36 MATERIALS SCIENCE↗
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 OSTI2026

Uncertainty quantification for Joule heating processes in fibrous pore-resolved media

Joule heating (JH) is an energy-efficient and sustainable technique for heating materials. Its application for industrial heating, particularly, has been gaining attention due to its potential for increasing the yield of various chemical products. The process involves the use of heating elements (materials that are highly conductive electrically and thermally) to heat up other materials or substances. These conductors, however, can exhbit varying degrees of uncertainty due to non-linearities in their temperature-dependent properties, which could result in variable material behavior. In this work, we carry out uncertainty quantification (UQ) at the pore scale to describe the uncertainty of such materials. In so doing, we applied the non-intrusive polynomial chaos expansion (PCE) technique to quantify the uncertainty within the system. The steady state Joule heating equation was solved numerically at the pore scale mimicking conditions within a heating chamber for propane dehydrogenation, and various electro-thermal profiles were obtained. We also examined the effect of the number of sampling points (20 – 100) and order of the PCE coefficients (2 – 5) on the accuracy of the temperature evaluations. The results were then benchmarked with the standard Monte Carlo (MC) method. The average temperature of the 4th-order global PCE showed good agreement with the MC results (which were positively skewed). Orders greater than 4 gave an underestimation of the temperatures while predictions for the peak temperature improved as the number of sampling points increased.

Fagbemi, Samuel [ORNL] (ORCID:0000000236995025)↗
DOE OSTI2026

Dynamically reconfigurable topological routing in nonlinear photonic systems

The propagation path of topologically protected states is bound to the interface between regions with different topology, and as such, the functionality of linear photonic devices leveraging these states is fixed during fabrication. Here, we propose a mechanism for dynamic control over a driven dissipative system’s local topology, yielding reconfigurable topological interfaces and thus tunable paths for protected routing. We illustrate our approach in non-resonantly pumped polariton lattices, where the nonlinear interaction between the polaritons and the exciton reservoir due to non-resonant pumping can yield picosecond-scale changes in the propagation paths of the chiral edge states. To analytically confirm the numerically observed topological dynamics, we generalize the spectral localizer framework to non-linear non-Hermitian Chern materials and apply this framework to a continuous model of the polariton system based on a driven-dissipative Gross-Pitaevskii equation. In doing so, we show that the local changes in the polariton lattice’s topology are captured by a local Chern marker. Looking forward, we anticipate such reconfigurable topological routing will enable the realization of novel classes of topological photonic devices.

Wong, Stephan [Sandia National Laboratories (SNL-N↗
DOE OSTI2025

Coupling to rotational manifolds to improve gas-phase pump–probe spectroscopic models

The physical picture of gas-phase optical transitions is normally presented as an isolated two-level system balanced by upward and downward processes. Isolated models assume a phenomenological treatment of collisional dephasing but do not strictly account for collisional population exchange with the rotational baths. While this assumption is valid under low-intensity conditions, where excitation is rate-limiting, isolated models can deviate from Beer’s Law at sufficient pressures and monochromatic intensities when both collisional broadening and power broadening are comparable to (or greater than) lifetime broadening, which are not uncommon conditions for cavity enhanced spectroscopies in the mid-IR spectral range. Although this problem has been addressed by rate-equation models for linear absorption measurements, a general treatment for multi-level quantum mechanical models suitable for non-linear absorption measurements (two-photon/two-color/pump–probe) is lacking. Isolated models require physical parameter inputs that disagree with expected values by at least an order of magnitude. These non-physical models undermine the ability to predict non-linear signal strengths under untested conditions and thereby limit the potential to optimize the sensitivity of non-linear spectroscopies and to expand their analytical applications (e.g., new analytes and/or buffer gases, changes in cavity free-spectral-range, changes in intracavity powers or wavelengths, and accurate investigation of physical phenomena). In this study, we derive bath-coupled models for gaseous pump–probe spectroscopy by application of the quantum Lindblad equation and detailed balance. Bath-coupled models are shown to fit data consistently across variations in intensity and agree with all physically expected values.

Cavity ring-down spectroscopy↗
DOE OSTI2025

Impact of high-temperature annealing on hafnia-silica composite coatings deposited via ion beam sputtering for high-peak power 1064 nm lasers

The maximum power handling fluence of high-peak and average power laser systems is often limited by the laser damage of the coatings on optical components. Furthermore, these multilayer dielectric coatings are limited in their maximum power handling due to laser-damage-prone defects in the lower optical bandgap, higher optical index material. Some of these defects can be mitigated by thermal annealing to high temperatures, which can greatly reduce the linear absorbative precursors. Typically, hafnia and silica are the materials of choice for high-peak and average power laser systems in the ultraviolet through infrared spectral range; however, hafnia crystalizes readily when annealed at high temperatures. In this study, we prepare composite HfO 2 -SiO 2 coatings by co-sputtering hafnia and silica in an ion beam sputtering system and compare them to pure hafnia-based coatings. We demonstrate that crystallinity in hafnia can be completely suppressed when it is mixed with silica, such as the composite coatings in this study. High reflectors were fabricated and annealed, demonstrating that the multilayer dielectric stacks can survive high-temperature annealing and exhibit an excellent linear absorption of 0.2 +/- 1 ppm at 1064 nm. Short- and long-pulse laser damage was explored, demonstrating the complex relationship between linear absorption and the non-linear absorption which drives pulsed laser damage. These results provide an excellent route to the creation of very low linear absorption optical coatings, which also utilize low scattering materials that are best suited for high-peak and average power applications.

Harthcock, Colin [Lawrence Livermore National Labo↗