SEARCH · Engineering Papers

Results for “Nonlinear problems”

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 451 records · Page 25

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Fourier-MIONet: Fourier-enhanced multiple-input neural operators for multiphase modeling of geological carbon sequestration

Geologic carbon sequestration (GCS) is a safety-critical technology that aims to reduce the amount of carbon dioxide in the atmosphere, which also places high demands on reliability. Multiphase flow in porous media is essential to understand CO 2 migration and pressure fields in the subsurface associated with GCS. However, numerical simulation for such problems in 4D is computationally challenging and expensive, due to the multiphysics and multiscale nature of the highly nonlinear governing partial differential equations (PDEs). It prevents us from considering multiple subsurface scenarios and conducting real-time optimization. Here, we develop a Fourier-enhanced multiple-input neural operator (Fourier-MIONet) to learn the solution operator of the problem of multiphase flow in porous media. Fourier-MIONet utilizes the recently developed framework of the multiple-input deep neural operators (MIONet) and incorporates the Fourier neural operator (FNO) in the network architecture. Once Fourier-MIONet is trained, it can predict the evolution of saturation and pressure of the multiphase flow under various reservoir conditions, such as permeability and porosity heterogeneity, anisotropy, injection configurations, and multiphase flow properties. Compared to the enhanced FNO (U-FNO), the proposed Fourier-MIONet has 90% fewer unknown parameters, and it can be trained in significantly less time (about 3.5 times faster) with much lower CPU memory (<15%) and GPU memory (<35%) requirements, to achieve similar prediction accuracy. In addition to the lower computational cost, Fourier-MIONet can be trained with only 6 snapshots of time to predict the PDE solutions for 30 years. Furthermore, we observed that Fourier-MIONet can maintain good accuracy when predicting out-of-distribution (OOD) data. The excellent generalizability of Fourier-MIONet is enabled by its adherence to the physical principle that the solution to a PDE is continuous over time. Furthermore, the developed Fourier-MIONet makes it possible to solve the long-time evolution of geological carbon sequestration in a large-scale three-dimensional space accurately and efficiently.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Nonlinear control of the minimum safety factor in tokamaks by optimal allocation of spatially moving electron cyclotron current drive

The minimum value of the safety factor profile is related to the magnetohydrodynamic (MHD) stability of the plasma confined in a tokamak. Therefore, active control of the minimum safety factor may mitigate MHD instabilities that can degrade or even terminate plasma confinement. Typically, in most tokamak scenarios, the minimum safety factor evolves spatially with time, i.e., the location at which the safety factor achieves the minimum value changes with time. In addition to the inherent nonlinearities in the minimum safety factor evolution, its spatial variation makes the control design challenging. In particular, complexity in control design may arise from the need for time-dependent nonlinear models that account for spatial variation of the minimum safety factor. Furthermore, the minimum safety factor may drift to locations where the actuator authority is low. The problem of minimum safety factor control with target location tracking and moving electron cyclotron current drive (ECCD) is addressed in this work. A nonlinear time-dependent model that incorporates the spatial variation of the minimum safety factor is presented. A nonlinear controller based on optimal feedback linearization is developed to track a target minimum safety factor. The proposed controller treats the ECCD position as a controllable variable. In other words, the controller prescribes the ECCD position (in addition to the non-inductive powers) in real time based on an optimal criterion that is defined a priori. This work also presents the steps necessary to integrate the minimum safety factor controller with a total energy controller to achieve multiple control objectives simultaneously. In conclusion, the proposed integrated control algorithm is tested using nonlinear simulations in the Control Oriented Transport SIMulator (COTSIM) for a DIII-D tokamak scenario.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Mixed-strategy approach to band-edge analysis and modeling in semiconductors

In semiconductor physics, the Urbach rule describes an exponential dependence of the UV-Vis absorption coefficient on photon energy at the band edge. The rule is so ubiquitously applied that a single common mechanism, variability in the local potential, is generally considered to be the cause of exponential band edges. Variability in the local potential may be caused by alloy inhomogeneity, structural defects, impurities, thermal excitation, or any other mechanism that interrupts ideal crystallinity and creates disorder. A crucial point of the Urbach rule is its development for band edges which approach linearity when plotted on a logarithmic scale. The rule was not developed in the context of nonlinear band edges with complex fine structure. Even so, this model has been applied extensively to band edges with such features. Using differential analysis and deconvolution, we develop an analytical technique capable of separating the components responsible for the fine structure observed. This approach reduces challenging signal-artifact convolution problems to a simple linear combination. Furthermore, we develop a model that can be applied to the nonlinear band-edge signals extracted. In this work, we show that in the limiting case of linearity, our model is consistent with and produces the same results as the linear Urbach model. Our model is capable of bypassing “tailing” as an indirect measure of the spread in the local potential and gives a direct measure of this quantity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Implicit shock tracking for unsteady flows by the method of lines

Here, a recently developed high-order implicit shock tracking (HOIST) framework for resolving discontinuous solutions of inviscid, steady conservation laws is extended to the unsteady case. Central to the framework is an optimization problem which simultaneously computes a discontinuity-aligned mesh and the corresponding high-order approximation to the flow, which provides nonlinear stabilization and a high-order approximation to the solution. This work extends the implicit shock tracking framework to the case of unsteady conservation laws using a method of lines discretization via a diagonally implicit Runge-Kutta method by “solving a steady problem at each timestep”. We formulate and solve an optimization problem that produces a feature-aligned mesh and solution at each Runge-Kutta stage of each timestep, and advance this solution in time by standard Runge-Kutta update formulas. A Rankine-Hugoniot based prediction of the shock location together with a high-order, untangling mesh smoothing procedure provides a high-quality initial guess for the optimization problem at each time, which results in rapid convergence of the sequential quadratic programing (SQP) optimization solver. This method is shown to deliver highly accurate solutions on coarse, high-order discretizations without nonlinear stabilization and recover the design accuracy of the Runge-Kutta scheme. We demonstrate this framework on a series of inviscid, unsteady conservation laws in both one- and two- dimensions. We also verify that our method is able to recover the design order of accuracy of our time integrator in the presence of a strong discontinuity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Analytical Large-Signal Modeling of Inverter-Based Microgrids With Koopman Operator Theory for Autonomous Control

The microgrid (MG) plays a crucial role in the energy transition, but its nonlinearity presents a significant challenge for large-signal power systems studies in the electromagnetic transient (EMT) time scale. In this paper, we develop a large-signal linear MG model that considers the detailed dynamics of the primary and zero-control levels based on the Koopman operator (KO) theory. Firstly, a set of observable functions is carefully designed to capture the nonlinear dynamics of the MG. The corresponding linear KO is then analytically derived based on these observables, resulting in the linear representation of the original nonlinear MG with observables as the new coordinate. The influence of external input on the system dynamics is also considered during the derivation, enabling control of the MG. We solve the voltage control problem using the traditional linear quadratic integrator (LQI) method to demonstrate that textbook linear control techniques can accurately control the original nonlinear MG via the developed KO-linearized MG model. Our proposed KO linearization method is generic and can be easily extended for different control objectives and MG structures using our analytical derivation procedure. Here, we validate the effectiveness of our methodology through various case studies.

25 ENERGY STORAGE↗

DOE OSTI2020

Relative importance of nonlinear electron-phonon coupling and vertex corrections in the Holstein model

Determining the range of validity of Migdal’s approximation for electron-phonon (e-ph) coupled systems is a long-standing problem. Many attempts to answer this question employ the Holstein Hamiltonian, where the electron density couples linearly to local lattice displacements. When these displacements are large, however, nonlinear corrections to the interaction must also be included, which can significantly alter the physical picture obtained from this model. Using determinant quantum Monte Carlo and the self-consistent Migdal approximation, we compared superconducting and charge-density-wave correlations in the Holstein model with and without second-order nonlinear interactions. We find a disagreement between the two cases, even for relatively small values of the e-ph coupling strength, and, importantly, that this can occur in the same parameter regions where Migdal’s approximation holds. Our results demonstrate that questions regarding the validity of Migdal’s approximation go hand in hand with questions of the validity of a linear e-ph interaction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Integrability technique for fluid flow induced deformation of a boundary hair

The deformation of a dense carpet of hair due to Stokes flow in a channel can be described by a nonlinear integrodifferential equation for the shape of a single hair, which possesses several solutions for a given choice of parameters. Although it was posed in a previous study and it bears a resemblance to the pendulum problem from mechanics, this equation has not been analytically solved until now. Despite the presence on an integral with a nonlinear functional dependence on the dependent variable, the system is integrable. We compare the analytically obtained solution to a finite-difference numerical approach, identify the physically realizable solution branch, and briefly study the solution structure through a conserved energylike quantity. Time-dependent fluid-structure interactions are a rich and complex subject to investigate, and we argue that the solution discussed herein can be used as a basis for understanding these systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Bayesian sequential optimal experimental design for nonlinear models using policy gradient reinforcement learning

We present a mathematical framework and computational methods for optimally designing a finite sequence of experiments. This sequential optimal experimental design (sOED) problem is formulated as a finite-horizon partially observable Markov decision process (POMDP) under a Bayesian setting and with information-theoretic utilities. The formulation is general and may accommodate continuous random variables, non-Gaussian posteriors, and nonlinear forward models. The sOED design policy incorporates elements of feedback and lookahead simultaneously, and we show it to generalize the commonly-used batch and greedy design strategies. We solve for the sOED policy using the policy gradient (PG) method from reinforcement learning, and provide a derivation for the PG expression in the sOED context. Adopting an actor-critic approach, the policy and value functions are parameterized using deep neural networks and improved via PG estimates produced from simulated episodes of designs and observations. The new PG-sOED algorithm is first validated on a linear-Gaussian benchmark, and then compared against other design baselines on a sensor movement problem for contaminant source inversion in a convection-diffusion field. As a result, we provide explanation for the policy behaviors using knowledge of the underlying physical process.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Measurement of the Alfvén Wave Parametric Decay Instability Growth Rate

Alfvén waves, a fundamental mode of magnetized plasmas, are ubiquitous in space and laboratory plasmas. The nonlinear behavior of these modes is thought to play a key role in important problems in space plasma, such as the heating of the solar corona and solar wind turbulence. In particular, theoretical predictions show that these Alfvén waves may be unstable to various parametric instabilities, but space observations of these processes are limited. We demonstrate the first measurement of the Alfvén wave parametric decay instability (PDI) growth rate. Experiments are conducted on the Large Plasma Device at UCLA in which a high amplitude 𝛿⁢𝐵/𝐵 0 ∼ 0.7% pump Alfvén wave is launched from one end of the device and a smaller seed Alfvén wave is launched from the other side. When the frequency of the seed wave is chosen to match the backward wave expected from PDI, damping of the seed wave is reduced. We compare this reduction in damping to the theoretically expected PDI growth rate while accounting for acoustic mode damping. Results show agreement between measurements and theoretical predictions. As a result, this not only provides critical validation for PDI theories and simulations that could help interpret future space observations but also suggests a new way of studying similar nonlinear wave phenomena.

Alfvén waves↗

DOE OSTI2024

Spatiotemporal dispersion compensation for a 200-THz noncollinear optical parametric amplifier

A noncollinear optical parametric amplifier (NOPA) can produce few-cycle femtosecond laser pulses that are ideally suited for time-resolved optical spectroscopy measurements. However, the nonlinear-optical process giving rise to ultrabroadband pulses is susceptible to spatiotemporal dispersion problems. Here, we detail refinements, including chirped-pulse amplification (CPA) and pulse-front matching (PFM), that minimize spatiotemporal dispersion and thereby improve the properties of ultrabroadband pulses produced by a NOPA. Here, the description includes a rationale behind the choices of optical and optomechanical components, as well as assessment protocols. We demonstrate these techniques using a 1 kHz, second-harmonic Ti:sapphire pump configuration, which produces ~5-fs duration pulses that span from about 500 to 800 nm with a bandwidth of about 200 THz. To demonstrate the utility of the CPA-PFM-NOPA, we measure vibrational quantum beats in the transient–absorption spectrum of methylene blue, a dye molecule that serves as a reference standard.

47 OTHER INSTRUMENTATION↗

DOE OSTI2020

Mean-field control for efficient mixing of energy loads

In this work, we pose an engineering challenge of controlling an ensemble of energy devices via coordinated, implementation-light, and randomized on/off switching as a problem in nonequilibrium statistical mechanics. We show that mean-field control with nonlinear feedback on the cumulative consumption, assumed available to the aggregator via direct physical measurements of the energy flow, allows the ensemble to recover from its use in the demand response regime, i.e., transition to a statistical steady state, significantly faster than in the case of the fixed feedback. Moreover when the nonlinearity is sufficiently strong, one observes the phenomenon of “super-relaxation,” where the total instantaneous energy consumption of the ensemble transitions to the steady state much faster than the underlying probability distribution of the devices over their state space, while also leaving almost no devices outside of the comfort zone.

25 ENERGY STORAGE↗

DOE OSTI2021

Stability and Conservation properties of Hermite-based approximations of the Vlasov-Poisson System

Spectral approximation based on Hermite-Fourier expansion of the Vlasov-Poisson model for a collisionless plasma in the electro-static limit is provided by adding high-order artificial collision operators of Lenard-Bernstein type. These differential operators are suitably designed in order to preserve the physically-meaningful invariants (number of particles, momentum, energy). In view of time-discretization, stability results in appropriate norms are presented. In this study, necessary conditions link the magnitude of the artificial collision term, the number of spectral modes of the discretization, as well as the time-step. The analysis, carried out in full for the Hermite discretization of a simple linear problem in one-dimension, is then partly extended to cover the complete nonlinear Vlasov-Poisson model.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Nonlinear Model Predictive Control Based on Real-Time Iteration Scheme for Wave Energy Converters Using WEC-Sim

One of several challenges that wave energy technologies face is their inability to generate electricity cost-competitively with other grid-scale energy generation sources. Several studies have identified two approaches to lower the levelised cost of electricity: reduce the cost over the device's lifetime or increase its overall electrical energy production. Several advanced control strategies have been developed to address the latter. However, only a few take into account the overall efficiency of the power take-off (PTO) system, and none of them solve the optimisation problem that arises at each sampling time on real-time. In this paper, a detailed Nonlinear model predictive control (NMPC) approach based on the real-time iteration (RTI) scheme is presented, and the controller performance is evaluated using a time-domain hydrodynamics model (WEC-Sim). The proposed control law incorporates the PTO system's efficiency in a control law to maximise the energy extracted. The study also revealed that RTI-NMPC clearly outperforms a simple resistive controller.

model predictive control↗

DOE OSTI2020

Rapid data-driven model reduction of nonlinear dynamical systems including chemical reaction networks using ℓ 1 -regularization

We develop a new data-driven paradigm for efficient model reduction of a broad class of nonlinear dynamical systems. Our model reduction method directly enables the interpretation of key components of the dynamical system, unlike traditional projection-based model reduction methods that focus on reducing computational complexity more than interpretability. Our method is not application specific and is simple to implement on nonlinear dynamical systems arising from a variety of different fields. It requires minimal parameterization using a single parameter to trade-off between model complexity and estimation error. We use a data-driven paradigm to formulate model reduction as an efficient convex optimization problem that scales polynomially in the original size of the complex system, enabling systems with as many as thousands of components to be reduced in a matter of minutes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Nonlinear Model Predictive Control Based on Real-Time Iteration Scheme for Wave Energy Converters Using WEC-Sim: Preprint

model predictive control↗

DOE OSTI2022

Characterization of nonlinear ultrasonic waves behavior while interacting with poor interlayer bonds in large-scale additive manufactured materials

Over the past decades, researchers have developed several nonlinear ultrasonic techniques for quality control of materials commonly used in different applications. Owing to the superior sensitivity of nonlinear ultrasound waves to small defects such as micro-cracks, their applicability in different nondestructive testing (NDT) problems has been investigated in numerous studies. These studies utilize frequency domain analysis to detect the generation of higher harmonics because of the formation of defects in the inspected medium. Frequency domain analysis based on the Fourier transform is a significant approach used in linear systems; however, it may not perform adequately on nonlinear systems. Hence, studies on nonlinear dynamics and physics consider analyzing systems' behavior in the phase-space domain. In contrast to the frequency domain analysis, which can result in information loss, analysis in the phase-space domain retains all the information regarding a system's states. Here, we investigate the nonlinearities induced by poor interlayer bonding in polymer-based additive manufactured parts in the phase-space domain. It is convenient to characterize the nonlinearity in the phase-space domain because it provides a geometrical representation of a system's states. Two types of low quality interlayer bond are considered. The first type is simulated artificially while the second type is manufactured by reducing the bond quality during the printing process. The analysis verified that the received ultrasonic signals exhibit classical nonlinear behavior in the phase-space domain while interacting with simulated poor interlayer bonds. In addition, the results showed that the behavior of ultrasonic waves is amplitude-dependent and evolves into models that have not been previously reported. Furthermore, Largest Lyapunov Exponent (LLE) is used to quantify the behavior of nonlinear ultrasonic waves while interacting with poor interlayer bonds. Using LLE, it was observed that the divergence rate of the phase-space trajectories depends on the amplitude of the excitation. This observation quantitatively proves that nonlinear behavior of ultrasound while interacting with poor interlayer bonds can be amplitude-dependent. The results of both simulated and inherent poor interlayer bond cases showed that LLE can be used as a reliable quantitative damage-sensitive feature to detect and potentially characterize weak bonds, which are difficult to detect using conventional approaches. Additionally, the reported results in the phase-space domain provide a basis for proposing a new mathematical model for ultrasonic waves interacting with poor interlayer bonds.

36 MATERIALS SCIENCE↗

DOE OSTI2020

An Improved Genetic Algorithm approach to the Unit Commitment/Economic Dispatch problem

The deployment of new technologies, the importance of accurately modeling the dynamics of the generating units and the introduction of new policies are making the solution of the Unit Commitment/Economic Dispatch problem more and more complicated.In the present scenario, traditionally followed scheduling criteria might not lead to the optimal fleet configuration any more. In addition, most of the widely used techniques have limited capabilities at modeling the nonlinear dynamics of committed power plants. When realistic power systems comprising of several tens of generating units are modeled, the resulting optimization problem turns to be computationally intensive for the current computing capabilities. In this paper, an improved version of a GA-based optimization algorithm is presented. A detailed methodology aimed at obtaining a more efficient version of the GA, and a more detailed and accurate description of the flexible operation flexibility of the power plants is described.

genetic algorithm↗

DOE OSTI2024

Learning Constrained Parametric Differentiable Predictive Control Policies With Guarantees

We present differentiable predictive control (DPC), a method for offline learning of constrained neural control policies for nonlinear dynamical systems with performance guarantees. We show that the sensitivities of the parametric optimal control problem can be used to obtain direct policy gradients. Specifically, we employ automatic differentiation (AD) to efficiently compute the sensitivities of the model predictive control (MPC) objective function and constraints penalties. To guarantee safety upon deployment, we derive probabilistic guarantees on closed-loop stability and constraint satisfaction based on indicator functions and Hoeffding’s inequality. We empirically demonstrate that the proposed method can learn neural control policies for various parametric optimal control tasks. In particular, we show that the proposed DPC method can stabilize systems with unstable dynamics, track time-varying references, and satisfy nonlinear state and input constraints. Our DPC method has practical time savings compared to alternative approaches for fast and memory-efficient controller design. Specifically, DPC does not depend on a supervisory controller as opposed to approximate MPC based on imitation learning. We demonstrate that, without losing performance, DPC is scalable with greatly reduced demands on memory and computation compared to implicit and explicit MPC while being more sample efficient than model-free reinforcement learning (RL) algorithms.

97 MATHEMATICS AND COMPUTING↗