A reduced space simultaneous approach for optimization of nonlinear index-1 differential algebraic equations.
Abstract not provided.
Parker, Robert↗
SEARCH · Engineering Papers
Results for “Nonlinear equations”
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.
A reduced-order model for nonlinear radiative transfer problems based on moment equations and POD-Petrov-Galerkin projection of the normalized Boltzmann transport equation
A data-driven projection-based reduced-order model (ROM) for nonlinear thermal radiative transfer (TRT) problems is presented. The TRT ROM is formulated by (i) a hierarchy of low-order quasidiffusion (aka variable Eddington factor) equations for moments of the radiation intensity and (ii) the normalized Boltzmann transport equation (BTE). The multilevel system of moment equations is derived by projection of the BTE onto a sequence of subspaces which represent elements of the phase space of the problem. Exact closure for the moment equations is provided by the Eddington tensor. A Petrov-Galerkin (PG) projection of the normalized BTE is formulated using a proper orthogonal decomposition (POD) basis representing the normalized radiation intensity over the whole phase space and time. The Eddington tensor linearly depends on the solution of the normalized BTE. By linear superposition of the POD basis functions, a low-rank expansion of the Eddington tensor is constructed with coefficients defined by the PG projected normalized BTE. The material energy balance (MEB) equation is coupled with the effective gray low-order equations which exist on the same dimensional scale as the MEB equation. The resulting TRT ROM is structure and asymptotic preserving. A detailed analysis of the ROM is performed on the classical Fleck-Cummings (F-C) TRT multigroup test problem in 2D geometry. Numerical results are presented to demonstrate the ROM's effectiveness in the simulation of radiation wave phenomena. Importantly, the ROM is shown to produce solutions with sufficiently high accuracy while using low-rank approximation of the normalized BTE solution. Essential physical characteristics of supersonic radiation wave are preserved in the ROM solutions.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Extended Lagrangian Born–Oppenheimer molecular dynamics using a Krylov subspace approximation
It is shown how the electronic equations of motion in extended Lagrangian Born–Oppenheimer molecular dynamics simulations can be integrated using low-rank approximations of the inverse Jacobian kernel. This kernel determines the metric tensor in the harmonic oscillator extension of the Lagrangian that drives the evolution of the electronic degrees of freedom. The proposed kernel approximation is derived from a pseudoinverse of a low-rank estimate of the Jacobian, which is expressed in terms of a generalized set of directional derivatives with directions that are given from a Krylov subspace approximation. The approach allows a tunable and adaptive approximation that can take advantage of efficient preconditioning techniques. The proposed kernel approximation for the integration of the electronic equations of motion makes it possible to apply extended Lagrangian first-principles molecular dynamics simulations to a broader range of problems, including reactive chemical systems with numerically sensitive and unsteady charge solutions. This can be achieved without requiring exact full calculations of the inverse Jacobian kernel in each time step or relying on iterative non-linear self-consistent field optimization of the electronic ground state prior to the force evaluations as in regular direct Born–Oppenheimer molecular dynamics. We note the low-rank approximation of the Jacobian is directly related to Broyden’s class of quasi-Newton algorithms and Jacobian-free Newton–Krylov methods and provides a complementary formulation for the solution of nonlinear systems of equations.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Growing neoclassical tearing modes seeded via transient-induced-multimode interactions
Nonlinear extended magnetohydrodynamics simulations demonstrating seeding of neoclassical tearing modes (NTMs) via magnetohydrodynamics-transient-induced multimode interactions are presented. Simulations of NTMs are enabled by two recent NIMROD code developments: the implementation of heuristic neoclassical stresses and the application of transient magnetic perturbations (MPs) at the boundary. NTMs are driven unstable by the inherently pressure driven kinetic bootstrap current, which arises due to collisional viscosity between passing and trapped electrons. These simulations use heuristic closures that model the neoclassical electron and ion stresses. NTM growth requires a seed island, which is generated by a transiently applied MP in simulations. The capability is demonstrated using kinetic-based reconstructions with flow of a DIII-D ITER Baseline Scenario discharge (La Haye et al., in Proceedings IAEA FEC, 2020). The applied MP seeds a 2/1 NTM that grows in two phases: a slow growth phase followed by a faster robust growth phase. Additionally, an evolving sequence of higher order core modes are excited at first. Power transfer analysis shows that nonlinear interactions between the core modes and the 2/1 helps drive the initial slow growth. Once the induced 2/1 magnetic island reaches a critical width, the NTM transitions to faster robust growth, which is well described by the nonlinear modified Rutherford equation. Furthermore, this work highlights the role of nonlinear mode coupling in seeding NTMs.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
The extended auxiliary equation mapping method to determine novel exact solitary wave solutions of the nonlinear fractional PDEs
Abstract In this paper, some new nonlinear fractional partial differential equations (PDEs) have been considered.Three models are including the space-time fractional-order Boussinesq equation, space-time (2 + 1)-dimensional breaking soliton equations, and space-time fractional-order SRLW equation describe the behavior of these equations in the diverse applications. Meanwhile, the fractional derivatives in the sense of β -derivative are defined. Some fractional PDEs will convert to the considered ordinary differential equations by the help of transformation of β -derivative. These equations are analyzed utilizing an integration scheme, namely, the extended auxiliary equation mapping method. The different kinds of traveling wave solutions, solitary, topological, dark soliton, periodic, kink, and rational, fall out as a by-product of this scheme. Finally, the existence of the solutions for the constraint conditions is also shown. The outcome indicates that some fractional PDEs are used as a growing finding in the engineering sciences, mathematical physics, and so forth.
Engineering↗
Physics-Informed Neural Network Solution of Point Kinetics Equations for a Nuclear Reactor Digital Twin
A digital twin (DT) for nuclear reactor monitoring can be implemented using either a differential equations-based physics model or a data-driven machine learning model. The challenge of a physics-model-based DT consists of achieving sufficient model fidelity to represent a complex experimental system, whereas the challenge of a data-driven DT consists of extensive training requirements and a potential lack of predictive ability. We investigate the performance of a hybrid approach, which is based on physics-informed neural networks (PINNs) that encode fundamental physical laws into the loss function of the neural network. We develop a PINN model to solve the point kinetic equations (PKEs), which are time-dependent, stiff, nonlinear, ordinary differential equations that constitute a nuclear reactor reduced-order model under the approximation of ignoring spatial dependence of the neutron flux. The PINN model solution of PKEs is developed to monitor the start-up transient of Purdue University Reactor Number One (PUR-1) using experimental parameters for the reactivity feedback schedule and the neutron source. The results demonstrate strong agreement between the PINN solution and finite difference numerical solution of PKEs. We investigate PINNs performance in both data interpolation and extrapolation. For the test cases considered, the extrapolation errors are comparable to those of interpolation predictions. Extrapolation accuracy decreases with increasing time interval.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
Shadow Lagrangian dynamics for superfluidity
Motivated by a similar approach for Born-Oppenheimer molecular dynamics, this paper proposes an extended "shadow" Lagrangian density for quantum states of superfluids. The extended Lagrangian contains an additional field variable that is forced to follow the wave function of the quantum state through a rapidly oscillating extended harmonic oscillator. By considering the adiabatic limit for large frequencies of the harmonic oscillator, we can derive the two equations of motions, a Schrödinger-type equation for the quantum state and a wave equation for the extended field variable. The equations are coupled in a nonlinear way, but each equation individually is linear with respect to the variable that it defines. The computational advantage of this new system is that it can be easily discretized using linear time stepping methods, where we propose to use a Crank-Nicolson-type approach for the Schrödinger equation and an extended leapfrog scheme for the wave equation. Furthermore, the difference between the quantum state and the extended field variable defines a consistency error that should go to zero if the frequency tends to infinity. By coupling the time-step size in our discretization to the frequency of the harmonic oscillator we can extract an easily computable consistency error indicator that can be used to estimate the numerical error without additional costs. The findings are illustrated in numerical experiments.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
PETN spark-gap detonators
The well-developed theory of Lorentz plasma that is dominated by electron–ion interactions is used to calculate the PETN arc characteristics. The spark-gap discharge current is a ramp with 10 to 25 ns rise time to peak and remaining constant subsequently. The approximate formulas for the arc channel conductivity, arc temperature, arc radius, and shock pressure from the arc are obtained from a system of nonlinear ordinary differential equations, which is the similarity solution of hydrodynamic equations similar to the Braginskii approximation. These arc parameters are given for the peak current ranging from 100 A to 1000 A and with different rise times. Representative cases are compared to the nonlinear ordinary differential equation code results. The shock pressures at the peak current are comparable to those from a typical commercial EBW bridgewire burst reported in the literature; the arc radius at the peak current is comparable to a typical bridgewire diameter of 0.0375 mm (e.g., RISI detonators, RP-1, and RP-80). The relevant Pop-Plot for low-density PETN is converted into an empirical detonation criterion, which is applicable to explosives subject to shocks of variable pressure. Finally, this criterion is then used to determine the detonation thresholds, which are comparable with test data obtained by Tucker, et al.
42 ENGINEERING↗
A Segregated Approach for Modeling the Electrochemistry in the 3-D Microstructure of Li-Ion Batteries and Its Acceleration Using Block Preconditioners
Abstract Battery performance is strongly correlated with electrode microstructure. Electrode materials for lithium-ion batteries have complex microstructure geometries that require millions of degrees of freedom to solve the electrochemical system at the microstructure scale. A fast-iterative solver with an appropriate preconditioner is then required to simulate large representative volume in a reasonable time. In this work, a finite element electrochemical model is developed to resolve the concentration and potential within the electrode active materials and the electrolyte domains at the microstructure scale, with an emphasis on numerical stability and scaling performances. The block Gauss-Seidel (BGS) numerical method is implemented because the system of equations within the electrodes is coupled only through the nonlinear Butler–Volmer equation, which governs the electrochemical reaction at the interface between the domains. The best solution strategy found in this work consists of splitting the system into two blocks—one for the concentration and one for the potential field—and then performing block generalized minimal residual preconditioned with algebraic multigrid, using the FEniCS and the Portable, Extensible Toolkit for Scientific Computation libraries. Significant improvements in terms of time to solution (six times faster) and memory usage (halving) are achieved compared with the MUltifrontal Massively Parallel sparse direct Solver. Additionally, BGS experiences decent strong parallel scaling within the electrode domains. Last, the system of equations is modified to specifically address numerical instability induced by electrolyte depletion, which is particularly valuable for simulating fast-charge scenarios relevant for automotive application.
25 ENERGY STORAGE↗
Fast Solution of Fully Implicit Runge--Kutta and Discontinuous Galerkin in Time for Numerical PDEs, Part II: Nonlinearities and DAEs
Fully implicit Runge--Kutta (IRK) methods have many desirable accuracy and stability properties as time integration schemes, but high-order IRK methods are not commonly used in practice with large-scale numerical PDEs because of the difficulty of solving the stage equations. This paper introduces a theoretical and algorithmic framework for solving the nonlinear equations that arise from IRK methods (and discontinuous Galerkin discretizations in time) applied to nonlinear numerical PDEs, including PDEs with algebraic constraints. Several new linearizations of the nonlinear IRK equations are developed, offering faster and more robust convergence than the often-considered simplified Newton, as well as an effective preconditioner for the true Jacobian if exact Newton iterations are desired. Inverting these linearizations requires solving a set of block 2 x 2 systems. Under quite general assumptions, it is proven that the preconditioned 2 x 2 operator's condition number is bounded by a small constant close to one, independent of the spatial discretization, spatial mesh, and time step, and with only weak dependence on the number of stages or integration accuracy. Moreover, the new method is built using the same preconditioners needed for backward Euler-type time stepping schemes, so can be readily added to existing codes. The new methods are applied to several challenging fluid flow problems, including the compressible Euler and Navier--Stokes equations, and the vorticity-streamfunction formulation of the incompressible Euler and Navier--Stokes equations. Up to 10th-order accuracy is demonstrated using Gauss IRK, while in all cases fourth-order Gauss IRK requires roughly half the number of preconditioner applications as required by standard Singly diagonally implicit Runge--Kutta methods.
97 MATHEMATICS AND COMPUTING↗
Non-intrusive data-driven model reduction for differential–algebraic equations derived from lifting transformations
In this paper we present a non-intrusive data-driven approach for model reduction of nonlinear systems. The approach considers the particular case of nonlinear partial differential equations (PDEs) that form systems of partial differential–algebraic equations (PDAEs) when lifted to polynomial form. Such systems arise, for example, when the governing equations include Arrhenius reaction terms (e.g., in reacting flow models) and thermodynamic terms (e.g., the Helmholtz free energy terms in a phase-field solidification model). Using the known structured form of the lifted algebraic equations, the approach computes the reduced operators for the algebraic equations explicitly, using straightforward linear algebra operations on the basis matrices. The reduced operators for the differential equations are inferred from lifted snapshot data using operator inference, which solves a linear least squares regression problem. The approach is illustrated for the nonlinear model of solidification of a pure material. The lifting transformations reformulate the solidification PDEs as a system of PDAEs that have cubic structure. The operators of the lifted system for this solidification example have affine dependence on key process parameters, permitting us to learn a parametric reduced model with operator inference. Numerical experiments show the effectiveness of the resulting reduced models in capturing key aspects of the solidification dynamics.
42 ENGINEERING↗
Modeling the AC Power Flow Equations with Optimally Compact Neural Networks: Application to Unit Commitment
Nonlinear power flow constraints render a variety of power system optimization problems computationally intractable. Emerging research shows, however, that the nonlinear AC power flow equations can be successfully modeled using neural networks. These neural networks can be exactly transformed into mixed integer linear programs and embedded inside challenging optimization problems, thus replacing nonlinearities that are intractable for many applications with tractable piecewise linear approximations. Such approaches, though, suffer from an explosion of the number of binary variables needed to represent the neural network. Accordingly, this paper develops a technique for training an "optimally compact'' neural network, i.e., one that can represent the power flow equations with a sufficiently high degree of accuracy while still maintaining a tractable number of binary variables. We demonstrate the use of this neural network as an approximator of the nonlinear power flow equations by embedding it in the AC unit commitment problem, transforming the problem from a mixed integer nonlinear program into a more manageable mixed integer linear program. We use the 14-, 57-, and 89-bus networks as test cases and compare the AC-feasibility of commitment decisions resulting from the neural network, DC, and linearized power flow approximations. Our results show that the neural network model outperforms both the DC and linearized power flow approximations when embedded in the unit commitment problem. The neural network formulation most often selects a feasible unit commitment schedule, and furthermore, it only s
AC power flow↗
Phenomena-based graph representations and applications to chemical process simulation
Rapid and robust simulation of chemical processes is critical to conduct process design, optimization, techno-economic analysis, and sustainability analysis. Yet, efficiently solving simulation models remains a challenge due to the highly coupled and nonlinear nature of the underlying algebraic equations that capture the physical phenomena taking place in the process (e.g., material and energy conservation, phase equilibrium, reactions). In this work, we show that graph-theoretic representations of the physical phenomena within unit operations can help navigate and decompose equations to systematically identify alternative approaches for fast and robust numerical solutions. Specifically, we present a graph-theoretic abstraction that captures the connectivity between the model variables/equations and use this abstraction to group variables/equations into fundamental phenomena. We show that phenomena-based decomposition of the underlying equations can help decouple nonlinearities and enforce material/energy conservation at the process level to accelerate convergence. The proposed decomposition approach differs from the more traditional sequential modular simulation approach, in which equations are grouped and decomposed by unit operations. We implemented the phenomena-based decomposition in BioSTEAM—an open-source process simulation platform in Python—and demonstrated that this approach can converge a variety of separation process models. Compared to sequential modular simulation, the phenomena-based approach can converge idealized systems faster, but it can be slower for (or even fail to converge) highly coupled and nonideal process systems.
Convergence↗
Nonlinear evolution of quadratic gravity in 3 + 1 dimensions
We present a numerically stable system of ( 3 + 1 ) evolution equations for the nonlinear gravitational dynamics of quadratic-curvature corrections to General Relativity (quadratic gravity). We also report on the numerical implementation of these evolution equations. We recover a well-known linear instability and gather evidence that—aside from said instability—quadratic gravity exhibits a physically stable Ricci-flat subsector. In particular, we demonstrate that Teukolsky-wave perturbations of a Schwarzschild black hole as well as a full binary inspiral (evolved up to merger) remain Ricci flat throughout evolution. Furthermore, this suggests that, at least in vacuum, classical quadratic gravity can mimic General Relativity, even in the fully nonlinear strong-gravity regime.
79 ASTRONOMY AND ASTROPHYSICS↗
Numerical and experimental evaluation of mechanical performance of the multifunctional energy storage composites
This work presents numerical simulation methods to model the mechanical behavior of the multifunctional energy storage composites (MESCs), which consist of a stack of multiple thin battery layers reinforced with through-the-hole polymer rivets and embedded inside carbon fiber composite laminates. MESC has been demonstrated through earlier experiments on its exceptional behavior as a structural element as well as a battery. However, the inherent complex infrastructure of the MESC design has created significant challenges in simulation and modeling. A novel homogenization technique was adopted to characterize the multi-layer properties of battery material using physics-based constitutive equations combined with nonlinear deformation theories to handle the interface between the battery layers. Second, mechanical damage and failure modes among battery materials, polymer reinforcements, and carbon fiber-polymer interfaces were characterized through appropriate models and experiments. The model of MESCs has been implemented in a commercial finite element code in ABAQUS. A comparison of structural response and failure modes from numerical simulations and experimental tests are presented. The results of the study showed that the predictions of elastic and damage responses of MESCs at various loading conditions agreed well with the experimental data. © 2021
Materials Science↗
A Realistic Theory of Quantum Measurement
Abstract We propose that the ontic understanding of quantum mechanics can be extended to a fully realistic theory that describes the evolution of the wavefunction at all times, including during a measurement. In such an approach the wave equation should reduce to the standard wave equation when there is no measurement, and describe state reduction when the system is measured. The general wave equation must be nonlinear and nonlocal, and we require it to be time-symmetric; consequently, this approach is not a new interpretation but a new theory. The wave equation is an integrodifferential equation (IDE). The time symmetry requirement leads to a retrocausal approach, in which the wave equation is solved subject to initial and final conditions to determine history at intermediate times. We propose that different outcomes from (apparently) identically prepared experiments may result from uncontrolled parameters; both the nonlocality and the retrocausality of the theory imply that Bell’s Theorem cannot rule out such “hidden variables.” Beginning with Hamilton’s principle, we demonstrate the construction of such a theory by replacing the action with a functional designed to give rise to a nonlinear, nonlocal IDE as the wave equation. This IDE reduces to the standard wave equation (a differential equation) in the absence of a measurement, but exhibits state reduction to a single eigenvalue when the system interacts with another system with the properties of a measurement apparatus. We demonstrate several desirable features of this theory; for other properties we indicate their plausibility and possible avenues to a proof.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Model predictive control of nonlinear processes using neural ordinary differential equation models
Not Available
Luo, Junwei (ORCID:000900072609629X)↗
A nonlinear moment model for radiative transfer equation in slab geometry
Not Available
Fan, Yuwei (ORCID:0000000195666996)↗