Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear analysis”

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

Accelerating uncertainty quantification in incremental dynamic analysis using dimension reduction-based surrogate modeling

We propose a surrogate modeling framework based on dimension reduction to facilitate the quantification of seismic risk of structural systems in performance-based earthquake engineering. The framework adopts incremental dynamic analysis (IDA) for addressing hazard variability, and promotes significant computational efficiency improvement for propagating epistemic uncertainties associated with the structural models. It utilizes both linear and nonlinear dimension reduction approaches, equipped with inverse mappings, to learn a functional between the input parameter space (e.g., the epistemic uncertainties of the structure) to the high-dimensional output space created through the IDA implementation across different ground motions and seismic intensity levels. Polynomial chaos expansion is adopted as the surrogate model to learn this functional in the reduced space. A nine-story steel moment-resisting frame with uncertain structural properties is used as a testbed. Furthermore, we select the seismic fragility curves as a measure of the structure’s seismic performance, since it provides an estimate of the probability of entering specified damage states for given levels of ground shaking.

42 ENGINEERING↗
DOE OSTI2021

Sensitive parameter identification and uncertainty quantification for the stability of pipeline conveying fluid

In this study, several uncertainty quantification and sensitivity analysis methods are used to determine the most sensitive geometric and material input parameters of a cantilevered pipeline conveying fluid when uncertainty is introduced to the system at the onset of instability. The full nonlinear equations of motion are modeled using the extended Hamilton’s principle and then discretized using Galerkin’s method. A parametric study is first performed, and the Morris elementary effects are calculated to obtain a preliminary understanding of how the onset speed changes when each parameter is introduced to a ± 5% uncertainty. Then, four different input uncertainty distributions, mainly, uniform and Gaussian distribution, are chosen to investigate how input distributions affect uncertainty in the output. A convergence analysis is used to determine the number of samples needed to maintain simulation accuracy while saving the most computational time. Then, Monte Carlo simulations are run, and the output distributions for each input distribution at ± 1%, ± 3% and ± 5% input uncertainty range are found and discussed. Additionally, the Pearson correlation coefficients are evaluated for different uncertainty ranges. A final Monte Carlo study is performed in which single parameters are held constant while all others still have uncertainty. Overall, the flow speed at the onset of instability is the most sensitive to changes in the outer diameter of the pipe.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Gravity Well Commercial Economics Assessment: Potential Revenue and Cost: Cooperative Research and Development (Final Report)

In the Gravity Well Revenue Study, we evaluate the potential revenue from energy storage using historical energy-only electricity prices, forward-looking projections of hourly electricity prices, and actual reported revenue. This analysis examines the impact of storage duration and round-trip efficiency, as well as the location of the storage, on storage revenue within the current and projected U.S. power system. We also investigated the impact of round-trip efficiency on storage revenue. We found that the relationship between storage revenue and round-trip efficiency is nonlinear. The value of improved round-trip efficiency declines as round-trip efficiency increases. In the Gravity Well Future Cost Study, we applied learning curves to predict the future cost trajectory of Gravity Wells (GrWs). Two types of analysis were implemented. The first was a bottom-up analysis that used historical learning rates for cost components, such as motors and gearboxes, and cost categories (e.g., engineering and design, etc.) to determine the learning-by-doing based single-factor learning curve. The single factor learning curve expresses the relationship between the cost of GrW and the number of units deployed (or the cumulative capacity). In the second analysis, we predicted future GrW costs via a top-down approach. This approach accounts for historical cost trends in other renewable energy and storage technologies, which have similarities with GrWs. Using a multifactor learning curve that accounts for both intrinsic (cumulative capacity) and extrinsic (the elasticity in the price of steel) factors, we estimated the future cost of GrWs.

25 ENERGY STORAGE↗
DOE OSTI2024

Nonlinear dynamics of multiple Alfvén modes driven by trapped energetic ions in tokamaks: A triplet paradigm

Non-linear dynamics of multiple infernal Alfvén eigenmodes—a subset of global Alfvén eigenmodes in tokamak plasmas with extended low-shear central core [Marchenko et al., Phys. Plasmas 16, 092502 (2009)]—is studied. The analysis is carried out for a mode triplet with toroidal mode-numbers n=1, 2, 3. It was assumed that the n = 1 mode was linearly unstable due to precession resonance with trapped fast ions, whereas the other modes were linearly damped. The modes were coupled due to a non-linearity in a bounce-averaged drift kinetic equation for the distribution function of fast ions. Nonlinear equations for the mode amplitudes and phases are derived and solved numerically. It is found that the temporal evolution of the amplitudes and the phase (responsible for the frequency chirping) of the modes exhibit Hopf bifurcations to stable limit cycles. This can explain a synchronous cyclic destabilization of multiple modes in Alfvén avalanches (sudden growth of amplitudes of the mode cluster with different n and approximately equal frequency spacing) in NSTX and bursting modes in MAST—events, which resulted in enhanced loss of fast ions.

Marchenko, V. S. (ORCID:0000000261480533)↗
DOE OSTI2025

Towards Floquet Chern insulators of light

Topological photonics explores photonic systems that exhibit robustness against defects and disorder, enabled by protection from underlying topological phases. These phases are typically realized in linear optical systems and characterized by their intrinsic photonic band structures. Here we experimentally study Floquet Chern insulators in periodically driven nonlinear photonic crystals, where the topological phase is controlled by the polarization and the frequency of the driving field. Our transient sum-frequency generation measurements reveal strong hybridization of the Floquet photonic bands. The measured spectrum remains gapless under a linearly polarized drive but becomes gapped under a circularly polarized drive. Theoretical analysis confirms that the Floquet gap is topological, characterized by a non-zero Chern number—a consequence of time-reversal symmetry breaking induced by the circularly polarized driving field. Furthermore, this work offers opportunities to explore the role of classical optical nonlinearity in topological phases and their applications in nonlinear optoelectronics.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Autoencoder Based Analysis of RF Parameters in the Fermilab Low Energy Linac

Machine learning (ML) has the potential for significant impact on the modeling, operation, and control of particle accelerators due to its ability to model nonlinear behavior, interpolate on complicated surfaces, and adapt to system changes over time. Anomaly detection in particular has been highlighted as an area where ML can significantly impact the operation of accelerators. These algorithms work by identifying subtle behaviors of key variables prior to negative events. Efforts to apply ML to anomaly detection have largely focused on subsystems such as RF cavities, superconducting magnets, and losses in rings. However, dedicated efforts to understand how to apply ML for anomaly detection in linear accelerators have been limited. In this paper the use of autoencoders is explored to identify anomalous behavior in measured data from the Fermilab low-energy linear accelerator.

Edelen, Jonathan P. (ORCID:0000000215180652)↗
DOE OSTI2026

Quantitative Analysis and Prediction of Thermal Runaway Metrics of High-Nickel Oxide Cathodes by Machine Learning Models

The pursuit of higher energy density in lithium-ion batteries has made high-nickel (Ni) layered oxides leading cathode candidates for next-generation electric vehicles. However, their poor thermal stability, particularly at Ni contents ≥ 90%, increases the risk of cathode-initiated thermal runaway. Furthermore, we present a data-driven framework combining linear and nonlinear machine learning models to predict key thermal runaway descriptors from a high-throughput differential scanning calorimetry database. With cathode composition and state of charge (SOC) as input features, the ensemble model accurately predicts peak temperature, heat release, and peak heat flow. SHAP analysis identifies Ni content and SOC as the dominant factors controlling thermal runaway temperature, while SOC primarily governs heat release and peak heat flow. Al, Mg, and Mn improve thermal stability by strengthening metal–oxygen bonding and delaying structural transformation, whereas B mainly reduces heat release through surface passivation. Validation with a new cathode composition confirms accurate prediction of SOC-dependent thermal runaway behavior and critical SOC.

25 ENERGY STORAGE↗
DOE OSTI2021

On-chip non-Hermitian optical parametric amplifiers with a large bandwidth

Recently, our groups have introduced the notion of optical parametric amplification based on non-Hermitian phase matching wherein the incorporation of loss can lead to gain in this nonlinear optical process. Previous simulation results using second-order nonlinear optical coupled-mode theory have demonstrated the potential of this technique as an alternative to the stringent phase-matching condition, which is often difficult to achieve in semiconductor platforms. Here we fortify this notion for the case of third-order nonlinearity by considering parametric amplification in silicon nanowires and illustrate the feasibility of these devices by employing rigorous finite-difference time-domain analysis using realistic materials and geometric parameters. Particularly, we demonstrate that by systematic control of the optical loss of the idler in a four-wave mixing process, we can achieve efficient unidirectional energy conversion from the pump to the signal component even when the typical phase-matching condition is violated. Importantly, our simulations show that a signal gain of ∼ <#comment/> 9 d B for a waveguide length of a few millimeters is possible over a large bandwidth of several hundreds of nanometers ( ∼ <#comment/> 600 n m ). This bandwidth is nearly 2 orders of magnitude larger than what can be achieved in the conventional silicon-photonics-based four-wave mixing process.

Ahmed, Asif (ORCID:0000000152234392)↗
DOE OSTI2023

Aeroelastic Analysis Using Deforming Cartesian Grids

Ongoing work in air-vehicle design illustrates the potential of advanced concepts to provide significant improvements in efficiency; but with their incorporation of lightweight flexible structures, such configurations may require active control systems to ensure reliability and safety. However, many contemporary analysis methods are inefficient for aeroelastic analysis and design of such configurations. This paper describes the development of a new approach that automates the geometry setup, mesh generation, and assembly of fluid–structural coupling interfaces to enable efficient aeroelastic and aeroservoelastic analysis of advanced concepts. The core elements for this approach are a cut-cell Cartesian grid-based computational fluid dynamics solver, a nonlinear beam element structural model, a conservative fluid–structural interface treatment, and the formulation and implementation of a new deforming grid capability within the cut-cell Cartesian grid solver. In this paper, emphasis is on this latter component with detailed description given of the mesh motion strategy, evaluation of fluxes and structural loads at the surface, and computation of geometrical properties such as cell volume, directed face areas, centroids, and motion-induced fluxes for deforming Cartesian grids required to advance the flow states. Aeroelastic simulations exercising the capability show favorable agreement with data and predictions in the literature for subsonic and supersonic applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Comprehensive assessment of deep reinforcement learning approaches for economic dispatch in nuclear-driven microgrids

As the electrical grid integrates more variable renewable energy sources such as wind and solar, the demand for distributed and flexible systems to address this increased variability becomes critical. Nuclear-driven microgrids provide a promising solution by offering stable generation to complement intermittent renewables, ensuring grid reliability and operating efficiency. This paper proposes a recurrent deep reinforcement learning framework for optimal economic dispatch in a nuclear-powered microgrid integrating renewable energy sources, small modular reactors, battery storage systems, and balance-of-plant dynamics. A three-agent control architecture is developed, where demand and renewable energy agents act as forecasters, and a reinforcement learning-based dispatch agent performs real-time energy allocation. A nonlinear programming formulation is first used to generate an optimal baseline for benchmarking. The proposed dispatch controller, based on Proximal Policy Optimization enhanced with Long Short-Term Memory networks, exploits temporal correlations in system dynamics by taking advantage of the time series used as inputs to improve policy robustness under uncertainty. Comparative analysis against established deep reinforcement learning methods, including Proximal Policy Optimization with a feedforward architecture, Soft Actor-Critic, and Twin Delayed Deep Deterministic Policy Gradient, demonstrates superior performance. Numerical results indicate that the proposed controller achieves a 0.39% cost reduction relative to the nonlinear programming benchmark and outperforms other learning-based methods by generating additional revenue of up to 0.35%. All reinforcement learning controllers compute dispatch actions in less than 0.3 s, resulting in a computational speedup of more than three orders of magnitude over the nonlinear programming baseline. The findings of this paper highlight their applicability for real-time operation and control in nuclear-integrated microgrids under volatile operating conditions.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Spatio-spectral optical fission in time-varying subwavelength layers

Transparent conducting oxides are highly doped semiconductors that exhibit favourable optical features compared with metals, including reduced material losses, tuneable electronic and optical properties, and enhanced damage thresholds. Recently, the photonic community has renewed its attention towards these materials, recognizing their remarkable nonlinear optical properties in the near-infrared spectrum. The exceptionally large and ultrafast change in the refractive index, which can be optically induced in these compounds, extends beyond the boundaries of conventional perturbative analysis and makes this class of materials the closest approximation to a time-varying system. Here we report the spatio-spectral fission of an ultrafast pulse trespassing a thin film of aluminium zinc oxide with a non-stationary refractive index. By applying phase conservation to this time-varying layer, our model can account for both space and time refraction and explain, in quantitative terms, the spatial separation of both spectrum and energy. Our findings represent an example of extreme nonlinear phenomena on subwavelength propagation distances, which provides new insights into transparent conducting oxides’ transient optical properties. This can be critical for the ongoing research on photonic time crystals, on-chip generation of non-classical states of light, integrated optical neural networks, ultrafast beam steering and frequency-division multiplexing.

integrated optics↗
DOE OSTI2023

Ideal MHD induced temperature flattening in spherical tokamaks

This paper extends the analysis first presented in Jardin et al. [Phys. Rev. Lett. 128, 245001 (2022)] to more thoroughly examine the stability of spherical torus equilibrium to ideal magnetohydrodynamic (MHD) infernal modes and their nonlinear consequences. Here we demonstrate that in a 3D resistive magnetohydrodynamic (MHD) simulation of a NSTX discharge, anomalous transport can occur due to these instabilities. Here we generate a family of equilibrium of differing β and use this to show that these instabilities could explain the experimentally observed flattening of the electron temperature profile at modest β. The modes studied in this paper are found to occur with poloidal mode number m and toroidal mode number n when the ratio m/ n is in the range of 1.2–1.5, when the central safety factor is in this range or slightly lower, and when the central region has very low magnetic shear. Our analysis gives some insight as to why the unstable linear growth rates are oscillatory functions of the toroidal mode number n. We present a simulation of an initially stable configuration that passes through a stability boundary at a critical β as it is heated. We also show that a particular NSTX discharge is unstable to these modes over a timescale of several hundred ms. We conclude that these modes must be taken into account when performing predictive modeling. An appendix shows that similar modes can be found in R / a = 4 tokamaks for certain q-profiles and β values.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Sharing is caring: An extensive analysis of parameter-based transfer learning for the prediction of building thermal dynamics

In recent years deep neural networks have been proposed as a lightweight data-driven model to capture high-dimensional, nonlinear physical processes to predict building thermal responses. However, the need of a large amount of data for the training process of deep neural networks clashes with the potential limited data availability in most existing or new buildings. Transfer learning aims to enhance the performance of a target learner exploiting knowledge from related and similar environments. This study conducted a suite of experiments that leveraged 250 data-driven models based on a synthetic dataset of a building archetype to study the influence of data availability, energy efficiency level, occupancy and climate for the transfer process of thermal dynamics. The performance of the transfer learning process was compared against a classical machine learning approach. Here, the results suggest that building thermal dynamics can be effectively transferred under the same climatic conditions, increasing performance when dealing with different occupancy schedules, efficiency levels and low data availability. Furthermore, the paper compares the performance of both transfer learning and machine learning approaches in an online fashion, to support the implementation in real-world deployment.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Probability of Initiation in Neutron Transport

We discuss the numerical solution of the nonlinear integro-differential equation for the probability of a divergent neutron chain in a stationary system (i.e., the probability of initiation (POI)). We follow the development described in Bell’s classic paper on the stochastic theory of neutron transport. As noted by Bell, the linearized form of this equation resembles the linear adjoint neutron transport equation. A matrix formalism for the discretized steady state (or forward) neutron equation in slab geometry is first developed and is then used to derive the discrete adjoint equation. A main advantage of this discrete development is that the resulting discrete adjoint equation does not depend upon how the multigroup cross sections for the forward problem are obtained. That is, we derive the discrete adjoint directly from the discrete forward equations rather than discretizing directly the adjoint equation. This also guarantees that the discrete adjoint operator is consistent with the inner product used to define the adjoint operator. We discuss three approaches for the numerical solution of the POI equations, and present numerical results on several test problems. The three solution methods are a simple fixed-point iteration, a second approach that is akin to a nonlinear Power iteration, and a third approach which uses a Newton-Krylov nonlinear solver. We also give sufficient conditions to guarantee the existence and uniqueness of nontrivial solutions to our discrete POI equations when the discrete system is supercritical, and that only the trivial solution exists when the discrete system is subcritical. Our approach is modeled after the analysis presented for the continuous POI equations by Mokhtar-Kharroubi and Jarmouni-Idrissi, and by Pazy and Rabinowitz.

42 ENGINEERING↗
DOE OSTI2024

Nonlinear elasticity with the Shifted Boundary Method

Here, we propose a new unfitted/immersed computational framework for nonlinear solid mechanics, which bypasses the complexities associated with the generation of CAD representations and subsequent body-fitted meshing. This approach allows to speed up the cycle of design and analysis in complex geometry and requires relatively simple computer graphics representations of the surface geometries to be simulated, such as the Standard Tessellation Language (STL format). Complex data structures and integration on cut elements are avoided by means of an approximate boundary representation and a modification (shifting) of the boundary conditions to maintain optimal accuracy. An extensive set of computational experiments in two and three dimensions is included.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Data-Driven Modeling of a High Capacity Cryogenic System for Control Optimization

The Cryogenic Moderator System (CMS) is responsible for maintaining a steady flow of cold neutrons for numerous physics experiments at the Spallation Neutron Source (SNS) in Oak Ridge National Laboratory (ORNL). Sudden losses in beam power, known as beam trips, cause a major disturbance to the CMS due to large step changes in cooling demands. Ongoing efforts on upgrading the neutron beam power from 1.4 to 2.0MW are expected to generate larger transients that can further strain the CMS subsystems if they are not properly controlled. To manage such disturbances, four flow valves and one electric heater are adjusted by five decentralized proportional-integral-derivative (PID) controllers. However, the original PID gains were calibrated empirically based only on tracking performance and not based on disturbance rejection. To address this issue without compromising current CMS operations, a control-oriented model was developed to recalibrate the PID controllers offline. The zero-dimensional (0-D) model was based on simple physics-based principles and data-driven system identification techniques. The CMS was broken into several subsystems for analysis, each of which corresponds to a parametric model tied to the thermodynamic states of the working fluid. The model parameters were identified using the nonlinear least squares method where the residuals were calculated from available sensor data. Simulation results show that the proposed model can capture the dynamics of the CMS at steady state and during beam trips.

Maldonado Puente, Bryan↗
DOE OSTI2025

A generalizable machine learning-assisted fast Fourier transform algorithm to simulate the large strain phenomena in polycrystalline materials

Machine learning methods have shown initial promise in constitutive modeling for single crystals or homogenized polycrystals, delivering notable computational efficiency. However, existing machine learning-based constitutive models often lack generalizability, limiting their application across diverse boundary value problems. This study introduces a thermodynamics-informed artificial neural network model to accelerate rate-tangent crystal plasticity fast Fourier transform simulations for cross-scale deformation behaviors of polycrystals under complex loading. Our model integrates microstructural variability and local interactions effectively. To address local effects in each grain, we employ K-means clustering to group Gauss points within the microstructure into clusters assumed to be in similar mechanical states. This approach, based on self-clustering analysis, extends model scope from macroscopic stress response to the granular level, capturing mechanical responses and orientation evolution across grains. This reduces the number of nonlinear problems to solve, with cluster responses propagated throughout each group. The thermodynamics-based artificial neural network-extracted features are further processed using local material state clusters to account for history-dependent deformation and evolving microstructures. Additionally, representative volume element simulations with rate-tangent crystal plasticity fast Fourier transform provide reliable datasets for model training. The proposed model demonstrates high efficiency, accuracy, self-consistency, and enhanced generalizability in predicting strain–stress responses and orientation evolution at both individual grain and aggregate scales under complex loading conditions, such as biaxial tension and arbitrary loading scenarios.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Total kinetic energy release in the fast-neutron-induced fission of Np 237

The total kinetic energy (TKE) release in the fast neutron-induced fission of 237Np was measured for neutron energies from En = 2.6–100 MeV at the Los Alamos Neutron Science Center at the Weapons Neutron Research facility. The post-TKE release decreases nonlinearly with increasing incident neutron energy and can be represented as TKE(MeV) = (174.38 ± 0.72) – (5.11 ± 0.5821) log 10 E n for En > 1 MeV. Analysis of the fragment mass distributions indicates that the decrease in TKE with increasing E n is a consequence of the fading out of shell effects at high excitation energies (resulting in an increase of symmetric fission) and the decrease of the total kinetic energy associated with asymmetric fission with increasing E n .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗