Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “exponential integrators”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Navigating Integration: Key Challenges for Data Centers, Nuclear Stakeholders, and Utility Operators

The rapid expansion of data centers, driven by the exponential growth in data-processing and storage needs, presents significant challenges and opportunities for various stakeholders, including data center developers, nuclear energy providers, and utility companies. Data centers are projected to consume 6.7–12% of United States (U.S.) electricity by 2028, driven by artificial intelligence (AI) and cloud-computing demands. Nuclear energy offers reliability and dispatchable baseload power, but data centers need power now while nuclear still needs time to address siting, fast power ramping, and regulatory hurdles. Utilities must keep pace with the unprecedented acceleration of large load interconnection requests and urgently adapt to high-density loads while maintaining grid stability, reliability, and accelerating interconnection timelines. This report dives into these challenges and proposes key collaboration strategies to streamline data center integration that aligns with recent federal initiatives like America’s AI Action Plan and related executive orders that emphasize the importance of data center growth, nuclear energy expansion, and maintaining a competitive edge in the global AI race.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
DOE OSTI2020

Fast dynamic aperture optimization with forward-reversal integration

A fast dynamic aperture (DA) optimization method for storage rings has been developed through the use of reversal integration. Even if dynamical systems have an exact reversal symmetry, a numerical forward integration differs from its reversal. For a chaotic trajectory, cumulative round-off errors are scaled, which results in an exponential growth on the difference. The exponential effect is a generic chaos indicator which represents the sensitivity of the chaotic motion to its initial condition. The chaos indicator of the charged particle motion can be obtained by comparing the forward integrations of particle trajectories with corresponding reversals, a.k.a. “backward integrations.” The indicator is observable even through short-term particle tracking simulations. Therefore, adopting it as an objective function could speed up optimization. Finally, the DA of the National Synchrotron Light Source II storage ring, and another test diffraction-limited light source ring, were optimized using this method for the purpose of demonstration.

43 PARTICLE ACCELERATORS↗
NASA NTRS1987

Vibration suppression using a constrained rate-feedback Threshold control strategy

Quasi-closed form solutions are derived for the finite time, minimum force rate-feedback threshold controller to bring a system with or without known external disturbances back into an 'allowable' state manifold in finite time. The disturbances are assumed to be expandable in terms of Fourier series. The quasi-closed form solutions replace the solution of the two-point boundary value problem and definite integral constraints with the solution of algebraic equations and the calculation of matrix exponentials. Examples demonstrate the threshold control technique and compare the quasi-closed form solutions with MACSYMA generated exact solutions (for small system order) and with the numerical solution of the two-point boundary value problem.

Zimmerman, D. C.↗
NASA NTRS1991

Vibration suppression using a constrained rate-feedback-threshold control strategy

A finite time, minimum force rate-feedback-threshold controller is developed to bring a system with or without known external disturbances back into an 'allowable' state bound in finite time. The disturbances are assumed to be expandable in terms of Fourier series. The optimal control is defined by a two-point boundary value problem coupled to a set of definite integral constraints. Quasi-closed form solutions are derived which replace the solution of the two-point boundary value problem and definite integral constraints with the solution of algebraic equations and the calculation of matrix exponentials. Examples are provided which demonstrate the threshold control technique and compare the quasi-closed form solutions with numerical and MACSYMA generated exact solutions.

Zimmerman, D. C.↗
NASA NTRS1987

Investigation of the feasibility of an analytical method of accounting for the effects of atmospheric drag on satellite motion

An analytic technique for accounting for the joint effects of Earth oblateness and atmospheric drag on close-Earth satellites is investigated. The technique is analytic in the sense that explicit solutions to the Lagrange planetary equations are given; consequently, no numerical integrations are required in the solution process. The atmospheric density in the technique described is represented by a rotating spherical exponential model with superposed effects of the oblate atmosphere and the diurnal variations. A computer program implementing the process is discussed and sample output is compared with output from program NSEP (Numerical Satellite Ephemeris Program). NSEP uses a numerical integration technique to account for atmospheric drag effects.

Bozeman, Robert E.↗
NASA NTRS1972

Origin, evolution and present thermal state of the moon.

The relative absence of lunar volcanism in the last 3,000,000,000 years and the Apollo 15 heat flow measurement suggest that present-day temperatures in the moon are approximately steady-state to depths of about 100 km. An exponential distribution of heat sources with depth may then be scaled by equating the surface heat flow to the integrated heat production of this exterior shell. Presumed present-day interior temperatures, as well as the present-day surface heat flow of about 30 erg/sq cm-sec, may be obtained with an initial temperature roughly corresponding to the Apollo 11 basalt solidus and exponential scaling of heat sources. The concentrations of U for an originally homogeneous moon are estimated to be about 0.009 micrograms per gram, close to that measured for eucrites and inferred for primitive inclusions of the Allende meteorite. The estimated homogeneous concentrations of U, the chemistry of the lunar surface material and inferences to modest depth, and the short accretion time of the moon necessary to provide large-scale differentiation at 4.6 AE suggest that the moon had its origin in the rapid accretion of compounds first condensing from the protoplanetary nebula.

Hanks, T. C.↗
DOE OSTI2020

Fast Dynamic Aperture Optimization with Reversal Integration

A fast method for dynamic aperture (DA) optimization of storage rings has been developed through the use of reversal integration. Even if dynamical systems have an exact reversal symme- try, a numerical forward integration di ers from its reversal. For chaotic trajectories, cumulative round-o errors are scaled, which results in an exponential growth in the di erence. The exponential e ect, intrinsically associated with the Lyapunov exponent, is a generic indicator of chaos because it represents the sensitivity of chaotic motion to an initial condition. A chaos indicator of the charged particle motion is then obtained by comparing the forward integrations of particle trajectories with corresponding reversals, a.k.a. \backward integrations." The indicator was con rmed to be ob- servable through short-term particle tracking simulations. Therefore, adopting it as an objective function could speed up optimization. The DA of the National Synchrotron Light Source II storage ring, and another test di raction-limited light source ring, were optimized using this method for the purpose of demonstration.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Mocking the u-plane integral

Abstract The u -plane integral is the contribution of the Coulomb branch to correlation functions of $${\mathcal {N}}=2$$ N = 2 gauge theory on a compact four-manifold. We consider the u -plane integral for correlators of point and surface observables of topologically twisted theories with gauge group $$\mathrm{SU}(2)$$ SU ( 2 ) , for an arbitrary four-manifold with $$(b_1,b_2^+)=(0,1)$$ ( b 1 , b 2 + ) = ( 0 , 1 ) . The u -plane contribution equals the full correlator in the absence of Seiberg–Witten contributions at strong coupling, and coincides with the mathematically defined Donaldson invariants in such cases. We demonstrate that the u -plane correlators are efficiently determined using mock modular forms for point observables, and Appell–Lerch sums for surface observables. We use these results to discuss the asymptotic behavior of correlators as function of the number of observables. Our findings suggest that the vev of exponentiated point and surface observables is an entire function of the fugacities.

Mathematics↗
DOE OSTI2021

Complexity growth in integrable and chaotic models

We use the SYK family of models with N Majorana fermions to study the complexity of time evolution, formulated as the shortest geodesic length on the unitary group manifold between the identity and the time evolution operator, in free, integrable, and chaotic systems. Initially, the shortest geodesic follows the time evolution trajectory, and hence complexity grows linearly in time. We study how this linear growth is eventually truncated by the appearance and accumulation of conjugate points, which signal the presence of shorter geodesics intersecting the time evolution trajectory. By explicitly locating such “shortcuts” through analytical and numerical methods, we demonstrate that: (a) in the free theory, time evolution encounters conjugate points at a polynomial time; consequently complexity growth truncates at O($\sqrt{N}$), and we find an explicit operator which “fast-forwards” the free N-fermion time evolution with this complexity, (b) in a class of interacting integrable theories, the complexity is upper bounded by O(poly(N)), and (c) in chaotic theories, we argue that conjugate points do not occur until exponential times O(e N ), after which it becomes possible to find infinitesimally nearby geodesics which approximate the time evolution operator. Finally, we explore the notion of eigenstate complexity in free, integrable, and chaotic models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Exact general solutions for cosmological scalar field evolution in a background-dominated expansion

We derive exact general solutions (as opposed to attractor particular solutions) and corresponding first integrals for the evolution of a scalar field $\phi$ in a universe dominated by a background fluid with equation of state parameter wB. In addition to the previously examined linear [V($\phi$)=V 0 $\phi$] and quadratic [V($\phi$)=V 0 $\phi$ 2 ] potentials, we show that exact solutions exist for the power-law potential V($\phi$)=V 0 $\phi$ n with n=4(1+w B )/(1-w B )+2 and n=2(1+w B )/(1-w B ). These correspond to the potentials V($\phi$)=V 0 $\phi$ 6 and V($\phi$)=V 0 $\phi$ 2 for matter domination and V($\phi$)=V 0 $\phi$ 10 and V($\phi$)=V 0 $\phi$ 4 for radiation domination. The $\phi$ 6 and $\phi$ 10 potentials can yield either oscillatory or nonoscillatory evolution, and we use the first integrals to determine how the initial conditions map onto each form of evolution. The exponential potential yields an exact solution for a stiff/kination (w B =1) background. We use this exact solution to derive an analytic expression for the evolution of the equation of state parameter, $w_\phi$, for this case.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1973

Lp-stability (1 less than or equal to p less than or equal to infinity) of multivariable nonlinear time-varying feedback systems that are open-loop unstable

A class of multivariable, nonlinear time-varying feedback systems with an unstable convolution subsystem as feedforward and a time-varying nonlinear gain as feedback was considered. The impulse response of the convolution subsystem is the sum of a finite number of increasing exponentials multiplied by nonnegative powers of the time t, a term that is absolutely integrable and an infinite series of delayed impulses. The main result is a theorem. It essentially states that if the unstable convolution subsystem can be stabilized by a constant feedback gain F and if incremental gain of the difference between the nonlinear gain function and F is sufficiently small, then the nonlinear system is L(p)-stable for any p between one and infinity. Furthermore, the solutions of the nonlinear system depend continuously on the inputs in any L(p)-norm. The fixed point theorem is crucial in deriving the above theorem.

Callier, F. M.↗
NASA NTRS1977

Sampling errors in the measurement of rain and hail parameters

Attention is given to a general derivation of the fractional standard deviation (FSD) of any integrated property X such that X(D) = cD to the n. This work extends that of Joss and Waldvogel (1969). The equation is applicable to measuring integrated properties of cloud, rain or hail populations (such as water content, precipitation rate, kinetic energy, or radar reflectivity) which are subject to statistical sampling errors due to the Poisson distributed fluctuations of particles sampled in each particle size interval and the weighted sum of the associated variances in proportion to their contribution to the integral parameter to be measured. Universal curves are presented which are applicable to the exponential size distribution permitting FSD estimation of any parameters from n = 0 to n = 6. The equations and curves also permit corrections for finite upper limits in the size spectrum and a realistic fall speed law.

Gertzman, H. S.↗
DOE OSTI2020

Optimized Lie–Trotter–Suzuki decompositions for two and three non-commuting terms

Lie–Trotter–Suzuki decompositions are an efficient way to approximate operator exponentials exp ( t H ) when H is a sum of n (non-commuting) terms which, individually, can be exponentiated easily. They are employed in time-evolution algorithms for tensor network states, digital quantum simulation protocols, path integral methods like quantum Monte Carlo, and splitting methods for symplectic integrators in classical Hamiltonian systems. Here, we provide optimized decompositions up to order t 6 . The leading error term is expanded in nested commutators (Hall bases) and we minimize the 1-norm of the coefficients. For n = 2 terms, several of the optima we find are close to those in McLachlan (1995). Generally, our results substantially improve over unoptimized decompositions by Forest, Ruth, Yoshida, and Suzuki. We explain why these decompositions are sufficient to efficiently simulate any one- or two-dimensional lattice model with finite-range interactions. This follows by solving a partitioning problem for the interaction graph.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Grid Integration of Small-Scale Photovoltaic Systems in Secondary Distribution Network—A Review

The relative share of renewable energy, specifically the solar photovoltaic (PV), is increasing exponentially in the world electric energy sector. This is a cumulative result of reduction in the cost of solar panels, improvement in the panel efficiency, and advancement in the associated power electronics. Among different types of PV plants, installation of small-scale rooftop PV is growing rapidly due to direct end-user benefits and lucrative governmental incentives. There are various standards developed in regards to grid integration of PVs and other distributed generations (DGs). Different power converter topologies are developed to interface the PV panel with the utility grid. To keep up with the stringent regulations imposed by the standards, various control strategies and grid synchronization methods have been developed. This review article amalgamates and summarizes all of the aforementioned aspects of a grid-integrated PV system including various standards, power stage architectures, grid synchronization methods, operation under extreme events, and control methodologies, pertaining to small-scale PV plants. This article will help freshman researchers to gain some familiarity with the topic and introduce them to some of the key issues encountered in this field.

14 SOLAR ENERGY↗
NASA NTRS1993

Gaia and the colonization of Mars

The Gaia hypothesis states that the atmosphere, hydrosphere, surface sediments, and life on Earth behave dynamically as a single integrated physiological system. What has been traditionally viewed as the passive environment is a highly active, integral part of the gaian system. Aspects of the surface temperature and chemistry are regulated by the sum of life, the biota. Formulated first by James E. Lovelock, in the late 1960s, the Gaia hypothesis has been in the scientific literature for more than 25 years. Because of its properties of exponential growth and propagation, life is a powerful geologic force. A useful aspect of the Gaia idea is that it requires integration of scientific disciplines for the study of Earth. The recently touted Earth system science is broadly parallel with the gaian concept of the physiochemical regulation of Earth's surface. We discuss here, in a gaian context, the colonization of Mars by Earth organisms. Although colonizing Mars may be impossible, its accomplishment would be exactly equivalent to "the reproduction of Gaia by budding.".

NASA Discipline Exobiology↗
NASA NTRS1986

Design of multivariable controllers for robot manipulators

The paper presents a simple method for the design of linear multivariable controllers for multi-link robot manipulators. The control scheme consists of multivariable feedforward and feedback controllers. The feedforward controller is the minimal inverse of the linearized model of robot dynamics and contains only proportional-double-derivative (PD2) terms. This controller ensures that the manipulator joint angles track any reference trajectories. The feedback controller is of proportional-integral-derivative (PID) type and achieves pole placement. This controller reduces any initial tracking error to zero as desired and also ensures that robust steady-state tracking of step-plus-exponential trajectories is achieved by the joint angles. The two controllers are independent of each other and are designed separately based on the linearized robot model and then integrated in the overall control scheme. The proposed scheme is simple and can be implemented for real-time control of robot manipulators.

Seraji, H.↗
DOE OSTI2024

Gravitational index of the heterotic string

The fundamental heterotic string has a tower of BPS states whose supersymmetric index has an exponential growth in the charges. We construct the saddle-point of the gravitational path integral corresponding to this index. The saddle-point configuration is a supersymmetric rotating non-extremal Euclidean black hole. This configuration is singular in the two-derivative theory. We show that the addition of higher-derivative terms in four-dimensional N = 2 supergravity resolves the singularity. In doing so, we extend the recently-developed “new attractor mechanism” to include the effect of higher-derivative terms. Remarkably, the one-loop, four-derivative F-term contribution to the prepotential leads to a precise match of the gravitational and microscopic index. We also comment, using the effective theory near the horizon, on the possibility of a string-size near-extremal black hole. Our results clarify the meaning of different descriptions of this system in the literature. The thermal state transitions to a winding condensate and a gas of strings without ever reaching a small black hole, while the index is captured by the rotating Euclidean black hole solution and is constant and thus smoothly connected to the microscopic ensemble.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Quantifying Quantum Chaos through Microcanonical Distributions of Entanglement

A characteristic feature of “quantum chaotic” systems is that their eigenspectra and eigenstates display universal statistical properties described by random matrix theory (RMT). However, eigenstates of local systems also encode structure beyond RMT. To capture this feature, we introduce a framework that allows us to compare the properties of eigenstates in local systems with those of pure random states. In particular, our framework defines a notion of distance between quantum state ensembles that utilizes the Kullback-Leibler divergence to compare the microcanonical distribution of entanglement entropy (EE) of eigenstates with a reference RMT distribution generated by pure random states (with appropriate constraints). This notion gives rise to a quantitative metric for quantum chaos that not only accounts for averages of the distributions but also higher moments. The differences in moments are compared on a highly resolved scale set by the standard deviation of the RMT distribution, which is exponentially small in system size. As a result, the metric can distinguish between chaotic and integrable behaviors and, in addition, quantify and compare the of chaos (in terms of proximity to RMT behavior) between two systems that are assumed to be chaotic. We implement our framework in local, minimally structured, Floquet random circuits, as well as a canonical family of many-body Hamiltonians, the mixed-field Ising model (MFIM). Importantly, for Hamiltonian systems, we find that the reference random distribution must be appropriately constrained to incorporate the effect of energy conservation in order to describe the ensemble properties of midspectrum eigenstates. The metric captures deviations from RMT across all models and parameters, including those that have been previously identified as strongly chaotic, and for which other diagnostics of chaos such as level spacing statistics look strongly thermal. In Floquet circuits, the dominant source of deviations is the second moment of the distribution, and this persists for all system sizes. For the MFIM, we find significant variation of the KL divergence in parameter space. Notably, we find a small region where deviations from RMT are minimized, suggesting that “maximally chaotic” Hamiltonians may exist in fine-tuned pockets of parameter space. Published by the American Physical Society 2024

Physics↗