SEARCH · Search NASA

Results for “nonlinearity”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 199 records · Page 11

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Nonlinear model predictive control for mode‐switching operation of reversible solid oxide cell systems

Abstract Solid oxide cells (SOCs) are a promising dual‐mode technology for the production of hydrogen through high‐temperature water electrolysis, and the generation of power through a fuel cell reaction that consumes hydrogen. Switching between these two modes as the price of electricity fluctuates requires reversible SOC operation and accurate tracking of hydrogen and power production set points. Moreover, a well‐functioning control system is important to avoid cell degradation during mode‐switching operation. In this article, we apply nonlinear model predictive control (NMPC) to an SOC module and supporting equipment and compare NMPC performance to classical proportional‐integral (PI) control strategies, while switching between the modes of hydrogen and power production. While both control methods provide similar performance across various metrics during mode switching, NMPC demonstrates a significant advantage in reducing cell thermal gradients and curvatures (mixed spatial‐temporal partial derivatives), thereby helping to mitigate long‐term degradation.

08 HYDROGEN↗

DOE OSTI2024

Numerical analysis of a time discretized method for nonlinear filtering problem with Lévy process observations

Abstract In this paper, we consider a nonlinear filtering model with observations driven by correlated Wiener processes and point processes. We first derive a Zakai equation whose solution is an unnormalized probability density function of the filter solution. Then, we apply a splitting-up technique to decompose the Zakai equation into three stochastic differential equations, based on which we construct a splitting-up approximate solution and prove its half-order convergence. Furthermore, we apply a finite difference method to construct a time semi-discrete approximate solution to the splitting-up system and prove its half-order convergence to the exact solution of the Zakai equation. Finally, we present some numerical experiments to demonstrate the theoretical analysis.

Mathematics↗

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 OSTI2025

Neural Active Manifolds: Nonlinear Dimensionality Reduction for Uncertainty Quantification

We present a new approach for nonlinear dimensionality reduction, specifically designed for computationally expensive mathematical models. We leverage autoencoders to discover a one-dimensional neural active manifold (NeurAM) capturing the model output variability, through the aid of a simultaneously learnt surrogate model with inputs on this manifold. Our method only relies on model evaluations and does not require the knowledge of gradients. The proposed dimensionality reduction framework can then be applied to assist outer loop many-query tasks in scientific computing, like sensitivity analysis and multifidelity uncertainty propagation. In particular, we prove, both theoretically under idealized conditions, and numerically in challenging test cases, how NeurAM can be used to obtain multifidelity sampling estimators with reduced variance by sampling the models on the discovered low-dimensional and shared manifold among models. Several numerical examples illustrate the main features of the proposed dimensionality reduction strategy and highlight its advantages with respect to existing approaches in the literature.

Autoencoders↗

DOE OSTI2026

Implementing a unified solver for nonlinearly constrained optimization

SQP and interior-point methods (also referred to as Lagrange-Newton methods) typically share key algorithmic components, such as strategies for computing descent directions and mechanisms that promote global convergence. Building on this insight, we introduce a unifying framework with eight building blocks that abstracts the workflows of Lagrange-Newton methods. We then present Uno, a modular C++ solver that implements our unifying framework and allows the automatic combination of a wide range of strategies with no programming effort from the user. Uno is meant to (1) organize mathematical optimization strategies into a coherent hierarchy; (2) offer a wide range of efficient and robust methods that can be compared for a given instance; (3) enable researchers to experiment with novel optimization strategies; and (4) reduce the cost of development and maintenance of multiple optimization solvers. Uno’s software design allows user to compose new customized solvers for emerging optimization areas such as robust optimization or optimization problems with complementarity constraints, while building on reliable nonlinear optimization techniques. We demonstrate that Uno is highly competitive against state-of-the-art solvers filterSQP, IPOPT, SNOPT, MINOS, LANCELOT, LOQO, and CONOPT on a subset of 429 small problems from the CUTE collection. Uno is available as open-source software under the MIT license at https://github.com/cvanaret/Uno and via its C, Julia, Python, Fortran, and AMPL interfaces.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Nonlinear behavior in high-frequency aggregate control of thermostatically controlled loads

Coordinated control of electric loads can provide valuable grid services, such as frequency regulation. However, due to the nonlinear characteristics of such load ensembles, it is important to systematically analyze their behavior and establish a thorough understanding of undesirable phenomena that can potentially arise. In this paper, we analyze the frequency response of an aggregate control scheme, with the goal of exploring controller performance limits. We show that rapid switching commands can induce oscillations in the power output due to the inherent lockout mechanism of the underlying devices. Here, we demonstrate that highly detailed aggregate models are required to capture such phenomena. Such models enable deeper understanding of the control boundaries and therefore play an important role in avoiding the introduction of undesirable effects on the grid.

24 POWER TRANSMISSION AND DISTRIBUTION↗

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

Asymptotic-preserving dynamical low-rank method for the stiff nonlinear Boltzmann equation

In kinetic theory, numerically solving the full Boltzmann equation is extremely expensive. This is because the Boltzmann collision operator involves a high-dimensional, nonlinear integral that must be evaluated at each spatial grid point and every time step. The challenge becomes even more pronounced in the fluid (strong collisionality) regime, where the collision operator exhibits strong stiffness, causing explicit time integrators to impose severe stability restrictions. In this paper, we propose addressing this problem through a dynamical low-rank (DLR) approximation. The resulting algorithm requires evaluating the Boltzmann collision operator only r 2 times, where r, the rank of the approximation, is much smaller than the number of spatial grid points. We propose a novel DLR integrator, called the XL integrator, which reduces the number of steps compared to the available alternatives (such as the projector splitting or basis update & Galerkin (BUG) integrator). For a class of problems including the Boltzmann collision operator which enjoys a separation property between physical and velocity space, we further propose a specialized version of the XL integrator, called the sXL integrator. This version requires solving only one differential equation to update the low-rank factors. Furthermore, the proposed low-rank schemes are asymptotic-preserving, meaning they can capture the asymptotic fluid limit in the case of strong collisionality. Our numerical experiments demonstrate the efficiency and accuracy of the proposed methods across a wide range of regimes, from non-stiff (kinetic) to stiff (fluid).

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Investigation of nonlinear mechanical behavior of two superfine-grained graphites with DIC assisted disc splitting test

Using standardized uniaxial tensile test specimens is not practical due to the limited volumes in irradiation capsules, molten salt degradation facilities, or oxidation apparatus. The ASTM International Standard Test Method for Tensile Strength Estimate by Disc Compression of Manufactured Graphite (ASTM D8289) was developed to provide a convenient way for estimating tensile strength. Unlike the traditional uniaxial tensile test in ASTM C749 (International Standard Test Method for Tensile Stress–Strain of Carbon and Graphite), which uses dog-bone shaped specimens larger than 12.95 mm × 120.65 mm, the ASTM D8289 standard uses smaller discs with diameters of 6–12.7 mm. A digital image correlation (DIC) system, which uses a full-field noncontact surface displacement measurement technique, was applied along with the ASTM D8289 disc splitting test on IG-110 samples and compared with previous measurements from Mersen 2114 samples. Results confirmed that the DIC technique can measure the surface displacement/strain on these small (Ø6 mm × 3 mm) graphite specimens with good repeatability. However, Mersen 2114 and IG-110 samples exhibited strain discrepancies when DIC measurements were compared with analytical and finite element simulation values. The loading history and strain results also indicated different mechanical behaviors between Mersen 2114 and IG-110, particularly the nonlinear behavior of the IG-110 samples. Good agreement was observed by comparing the splitting tensile strengths of two superfine-grained grades with results from other work. The specimen size effect is discussed when comparing the splitting tensile strength with corresponding uniaxial tensile strength of these two graphite grades.

Lin, Lianshan [ORNL] (ORCID:0000000203399219)↗

DOE OSTI2026

Physics-aware adaptive checkpointing with shadow systems for nonlinear PDE simulations

Large-scale simulations of nonlinear partial differential equations (PDEs) that exhibit strongly transient behavior and pattern-forming dynamics produce enormous amounts of data, which, even with modern storage systems, cannot be stored for later curation. Current I/O strategies either write dense time series of snapshots, which is often prohibitive in I/O and storage, or store a few checkpoints that enable restart but incur expensive recomputation cost and provide no control over post-restart error growth, especially when lossy compression is used. Moreover, most, if not all, existing strategies take no account of the actual physical state of the system. Here, we present a simple physics-aware I/O framework in which a low-cost shadow system adaptively triggers lossy checkpoints when the shadow system deviates from the fine-scale simulation. The shadow system can be a coarsened replica of the fine-scale simulation that evolves concurrently. This means that checkpoints are taken based on the physical state of the system: fewer checkpoints are triggered when the system is quiescent while more are taken when the system undergoes a rapid change. This type of behavior is observed in many systems such as Brusselator and FitzHugh–Nagumo. We illustrate that our framework maintains stable restarts, keeps fine-scale restart errors bounded by shadow errors, and reconstructs the time history with significantly lower error and storage than interpolating fixed-interval snapshots, with low-cost shadow replay and modest online synchronization overhead.

Gong, Qian [ORNL] (ORCID:0000000235704142)↗

DOE OSTI2024

Collocation methods for nonlinear differential equations on low-rank manifolds

We introduce new methods for integrating nonlinear differential equations on low-rank manifolds. These methods rely on interpolatory projections onto the tangent space, enabling low-rank time integration of vector fields that can be evaluated entry-wise. A key advantage of our approach is that it does not require the vector field to exhibit low-rank structure, thereby overcoming significant limitations of traditional dynamical low-rank methods based on orthogonal projection. To construct the interpolatory projectors, we develop a sparse tensor sampling algorithm based on the discrete empirical interpolation method (DEIM) that parameterizes tensor train manifolds and their tangent spaces with cross interpolation. Using these projectors, we propose two time integration schemes on low-rank tensor train manifolds. The first scheme integrates the solution at selected interpolation indices and constructs the solution with cross interpolation. The second scheme generalizes the well-known orthogonal projector-splitting integrator to interpolatory projectors. We demonstrate the proposed methods with applications to several tensor differential equations arising from the discretization of partial differential equations.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Implementing nonlinear optics from Off-Energy closed orbit at NSLS-II

To characterize the second-order (chromatic sextupole) magnet lattice with high precision, we implemented nonlinear optics from off-energy closed orbit (NOECO) tool based on the linear optics from closed orbit modulation (LOCOM) method, named LOCOM-NOECO. The preliminary numerical study indicates that 1–2% precision can be achieved for the calibration of chromatic sextupoles. Further, this accuracy could potentially help in resolving some long-standing challenges of NSLS-II (e.g., the discrepancy between the designed and measured dynamic apertures) if such high precision can be fulfilled. As an independent crosscheck, we also implemented NOECO based on the independent component analysis (ICA) method using turn-by-turn (TBT) BPM data, named ICA-NOECO. Both ICA-NOECO and LOCOM-NOECO have been successfully applied to identify the pre-dialed random errors of a chromatic sextupole family including five power supplies, and achieved the root mean square (RMS) residual error of 1% and peak error less than 2%. Moreover, to mitigate the chromatic sextupole error effect, we applied the correction and achieved significant improvements in the injection efficiency as well as the dynamic apertures.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2026

Joint Optimization of Multimodal Transit Frequency and Shared Autonomous Vehicle Fleet Size with Hybrid Metaheuristic and Nonlinear Programming

Shared autonomous vehicles (SAVs) bring competition to traditional transit services but redesigning multimodal transit network can utilize SAVs as feeders to enhance service efficiency and coverage. This paper presents an optimization framework for the joint multimodal transit frequency and SAV fleet size problem, a variant of the transit network frequency setting problem. The objective is to maximize total transit ridership (including SAV-fed trips and subtracting boarding rejections) across multiple time periods under budget constraints, considering endogenous mode choice (transit, point-to-point SAVs, driving) and route selection, while allowing for strategic route removal by setting frequencies to zero. Due to the problem’s non-linear, non-convex nature and the computational challenges of large-scale networks, we develop a hybrid solution approach that combines a metaheuristic approach (particle swarm optimization) with nonlinear programming for local solution refinement. To ensure computational tractability, the framework integrates analytical approximation models for SAV waiting times based on fleet utilization, multimodal network assignment for route choice, and multinomial logit mode choice behavior, bypassing the need for computationally intensive simulations within the main optimization loop. Applied to the Chicago metropolitan area’s multimodal network, our method illustrates a 33.3% increase in transit ridership through optimized transit route frequencies and SAV integration, particularly enhancing off-peak service accessibility and strategically reallocating resources.

Ng, Max↗

DOE OSTI2024

The Tribomechadynamics Research Challenge: Confronting blind predictions for the linear and nonlinear dynamics of a thin-walled jointed structure with measurement results

The present article summarizes the submissions to the Tribomechadynamics Research Challenge announced in 2021. The task was a blind prediction of the vibration behavior of a system comprising a thin plate clamped on two sides via bolted joints. Both geometric and frictional contact nonlinearities are expected to be relevant. Provided were the CAD models and technical drawings of all parts as well as assembly instructions. The main objective was to predict the frequency and damping ratio of the lowest-frequency mode as function of the amplitude. Many different prediction approaches were pursued, ranging from well-known methods to very recently developed ones. After the submission deadline, the system has been fabricated and tested. The aim of this article is to evaluate the current state of the art in modeling and vibration prediction, and to provide directions for future methodological advancements.

42 ENGINEERING↗

DOE OSTI2024

FeSi 4 P 4 and CoSi 3 P 3 : Hidden Gems of Ternary Tetrel Pnictides with Outstanding Nonlinear Optical Properties

Metal silicon phosphides have shown promise as nonlinear optical materials. To be practically useful and cheap, earth-abundant 3d transition metals are preferred over their scarcer and more expensive 4d and 5d counterparts. Here, we developed a synthetic method to produce polycrystalline bulk powders and millimeter-sized single crystals of ternary compounds FeSi 4 P 4 and CoSi 3 P 3 . Both studied compounds have noncentrosymmetric and chiral crystal structures with ordered Si/P arrangements as was confirmed by single-crystal X-ray diffraction and solid-state NMR. Despite the presence of the transition metal, FeSi 4 P 4 and CoSi 3 P 3 are semiconductors with direct band gaps of 1.3 and 1.6 eV, respectively, indicating low-spin d 6 electronic configuration for octahedral Fe 2+ and Co 3+ . Relative to reported sulfide materials, FeSi 4 P 4 and CoSi 3 P 3 small band gap semiconductors demonstrate an outstanding combination of second-harmonic generation (SHG) activity and laser damage threshold (LDT). Both studied materials are phase-matchable with a 2.09 μm laser and not only exhibit 2.5–3.0 times stronger SHG signal than that of the state-of-the-art AgGaS 2 standard but also demonstrate an LDT response of 2.3–2.5 times higher than that of AgGaS 2 (at 1.09 μm laser with a pulse width of 10 ns)-which is unprecedented for small band gap semiconductors.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Multiresolution Quantum Chemistry: Nonlinear Response Properties at the Basis Set Limit

We benchmark the accuracy of Dunning correlation-consistent Gaussian basis sets for computing frequencydependent second-order hyperpolarizabilities relevant to second-harmonic generation (SHG), using multiresolution analysis (MRA) as a reference. Basis set errors are analyzed using a unit-sphere representation of the effective hyperpolarizability vector, enabling direct assessment of directional error structure. We introduce a relative RMS total error metric that integrates directional deviations over the unit sphere and complement it with signed projection errors that distinguish over- and underestimation. Unsupervised clustering based on these signed directional metrics reveals four distinct convergence behaviors across a set of 68 molecules. Unitsphere visualizations of representative systems show that basis set errors are often highly anisotropic and localized along specific bond directions, even when global error measures appear small. Doubly augmented basis sets consistently outperform singly augmented ones, and core-polarization functions are required for uniform convergence in second-row systems. Overall, this work demonstrates that directional analysis combined with clustering provides a robust framework for understanding basis set convergence in nonlinear optical response properties.

Basis sets↗

DOE OSTI2025

Nonlinear Scaling of Water–Ion Interactions and Dynamics in Alkaline Solutions

Water-ion interactions govern many solvent properties critical to solution-phase chemistry and the behavior of liquid water. The water-ion interactions in alkaline conditions were probed using two-dimensional infrared spectroscopy (2D IR), small-angle x-ray scattering (SAXS), and nuclear magnetic resonance spectroscopy (NMR). Energy transfer between the donor molecule KSeCN, used as a 2D IR probe, and the acceptor molecule NaOD was used to track the average separation distance of ions in a D 2 O solution, while SAXS measurements showed effects in the bulk D 2 O solvent. Here, we observe consistent nonlinear scaling in the SeCN - and OD - average separation distance as a function of NaOD concentration while bulk solution D 2 O-to-D 2 O average separation distance remained highly linear. Ultrafast measurements of solution dynamics via 2D IR and polarization-selective pump-probe spectroscopy show consistent scaling in correlation times as a function of concentration. These results suggests that the SeCN - and OD - anions participate in a water-ion network that significantly reduces the degrees of freedom for the distribution of ions in solution below the standard stochastic for ion distribution.

2D IR Spectroscopy↗

DOE OSTI2026

Sedimentation and Nonlinear Trapping in Texas Reservoirs Identified Using Remote Sensing and Bathymetric Survey Records

Decreasing reservoir storage capacity due to sedimentation poses great challenges to aging U.S. reservoirs, as it reduces the efficacy and reliability of their socio‐economic services. However, systematic assessments of reservoir sedimentation rates and associated issues remain limited because of sparse and infrequent bathymetry survey data. In this study, we use remote sensing‐driven estimates of sediment concentrations to identify regions experiencing rapid reservoir capacity loss, as observed in repeated bathymetry surveys. Our analysis focuses on Texas, where one of the most reliable state‐level reservoir capacity loss data sets is available through a unique long‐term monitoring program by the Texas Water Development Board. We find that reservoirs with large storage capacities and high sedimentation rate are concentrated in Northeast Texas. We also show that reduced forest, increased barren land, and erosive soil properties are co‐varying with high reservoir sedimentation rates. Temporal changes in the longitudinal gradient of the remotely sensed sediment flux highlight the nonlinear nature of sediment trapping processes, and can be used to estimate the reservoir storage capacity loss over time. In combination with standard bathymetry surveys, our approach shows potential for more cost‐effective and frequent assessment of reservoir storage loss.

Lee, Jiyong [ORNL] (ORCID:0000000198957406)↗