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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Hilbert space fragmentation in a 2D quantum spin system with subsystem symmetries

We consider a 2D quantum spin model with ring-exchange interaction that has subsystem symmetries associated to conserved magnetization along rows and columns of a square lattice, which implies the conservation of the global dipole moment. The model is not integrable, but violates the eigenstate thermalization hypothesis through an extensive Hilbert space fragmentation, including an exponential number of inert subsectors with trivial dynamics, arising from kinetic constraints. While subsystem symmetries are quite restrictive for the dynamics, we show that they alone cannot account for such a number of inert states, even with infinite-range interactions. We present a procedure for constructing shielding structures that can separate and disentangle dynamically active regions from each other. Notably, subsystem symmetries allow the thickness of the shields to be dependent only on the interaction range rather than on the size of the active regions, unlike in the case of generic dipole-conserving systems.

Khudorozhkov, Alexey↗
DOE OSTI2020

Innovations in Renewable Energy Technologies, Systems, and Energy Analysis

Costs for renewable energy technologies have declined rapidly in the past decade and their use for residential, commercial, and utility scale electricity has therefore grown exponentially. As the energy mix diversifies, energy analysis becomes increasing important, requiring renewable energy to integrate with other reduced emission energy sources, such as small modular nuclear reactors and carbon capture and utilization. Energy analysis can also explore the expansion of clean energy into new sectors, such as industrial processes and electrification. The U.S. National Renewable Energy Laboratory (NREL) and the Joint Institute for Strategic Energy Analysis (JISEA), which is a partnership of NREL and universities, study the integration and application of clean energy technologies into existing and new sectors. This presentation includes an overview of NREL and JISEA and summarizes innovations in renewable energy technologies and energy systems analysis.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2020

DEEP Solar: Data DrivEn Modeling and Analytics for Enhanced System Layer ImPlementation

Realizing the SETO 2030 mission of reducing solar energy costs to 3-5 c/kWh will require innovative enabling research on effective, cost-efficient integration of local PV within distribution systems. However, the intermittent and variable nature of PVs compels operators to impose conservative hosting capacity constraints. Given the extremely high variability of (intermittent and unpredictable) solar energy generation, relaxing the capacity constraints (which are currently around 15%) and achieving 100% or greater integration of renewables will require a fundamental transformation of the power grid via the utilization of exponentially larger amounts of AMI enabled fine-grained data. To address the challenges in increasing the penetration of renewable energy based DERs, this project envisions an Enhanced System Layer (ESL) at the distribution network level that is reliable, cost-effective and scalable to millions of Distributed Energy Resources (DERs)/devices. This includes developing: 1) Transformative and highly scalable machine learning based predictive analytics tools that plug into distribution system planning and provide real-time situational awareness at the distribution level for short and long-term operational planning. The tools will be built using novel data-driven energy models of millions of active nodes with AMI, 2) Adaptive stochastic analysis and optimization algorithms for real-time grid operations, 3) Dynamic Scenario Analysis using parallel Cloudenabled implementations with < 1 minute computational cycle times.

14 SOLAR ENERGY↗
DOE OSTI2026

Quantum mechanics and observers for gravity in a closed universe

Recent arguments based on the quantum extremal surface formula or the gravitational path integral have given fairly compelling evidence that the Hilbert space of quantum gravity in a closed universe is one-dimensional and real. How can this be consistent with the complexity of our own experiences? In this paper we propose that the experiences of any observer Ob in a closed universe can be approximately described by a quantum mechanical theory with a Hilbert space whose dimension is roughly e S Ob , where S Ob is the number of degrees of freedom of Ob. Moreover we argue that the errors in this description are exponentially small in S Ob . We give evidence for this proposal by incorporating it into the gravitational path integral and the coding interpretation of holography in simple models and seeing that it works, and we explain how similar effects arise in black hole physics in appropriate circumstances.

AdS-CFT Correspondence↗
DOE OSTI2020

Subleading power resummation of rapidity logarithms: the energy-energy correlator in $ \mathcal{N} $ = 4 SYM

We derive and solve renormalization group equations that allow for the resummation of subleading power rapidity logarithms. Our equations involve operator mixing into a new class of operators, which we term the “rapidity identity operators”, that will generically appear at subleading power in problems involving both rapidity and virtuality scales. To illustrate our formalism, we analytically solve these equations to resum the power suppressed logarithms appearing in the back-to-back (double light cone) limit of the Energy-Energy Correlator (EEC) in \( \mathcal{N} \) = 4 super-Yang-Mills. These logarithms can also be extracted to \( \mathcal{O}\left({\alpha}_s^3\right) \) from a recent perturbative calculation, and we find perfect agreement to this order. Instead of the standard Sudakov exponential, our resummed result for the subleading power logarithms is expressed in terms of Dawson’s integral, with an argument related to the cusp anomalous dimension. We call this functional form “Dawson’s Sudakov”. Our formalism is widely applicable for the resummation of subleading power rapidity logarithms in other more phenomenologically relevant observables, such as the EEC in QCD, the p T spectrum for color singlet boson production at hadron colliders, and the resummation of power suppressed logarithms in the Regge limit.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Nonperturbative quantum gravity in a closed Lorentzian universe

We study how meaningful physical predictions can arise in nonperturbative quantum gravity in a closed Lorentzian universe. In such settings, recent developments suggest that the quantum gravitational Hilbert space is one-dimensional and real for each α-sector, as induced by spacetime wormholes. This appears to obstruct the conventional quantum-mechanical prescription of assigning probabilities via projection onto a basis of states. While previous approaches have introduced external observers or augmented the theory to resolve this issue, we argue that quantum gravity itself contains all the necessary ingredients to make physical predictions. We demonstrate that the emergence of classical observables and probabilistic outcomes can be understood as a consequence of partial observability: physical observers access only a subsystem of the universe. Tracing out the inaccessible degrees of freedom yields reduced density matrices that encode classical information, with uncertainties exponentially suppressed by the environment’s entropy. We develop this perspective using both the Lorentzian path integral and operator formalisms and support it with a simple microscopic model. Our results show that quantum gravity in a closed universe naturally gives rise to meaningful, robust predictions without recourse to external constructs.

AdS-CFT Correspondence↗
NASA NTRS1987

Balloon-borne observations of the development and vertical structure of the Antarctic ozone hole in 1986

The vertical distribution of ozone measured at McMurdo Station, Antarctica using balloon-borne sensors on 33 occasions during November 6, 1986 - August 25, 1986 is described. These observations suggest a highly structured cavity confined to the 12-20 km altitude region. In the 17-19 km altitude range, the ozone volume mixing ratio declined from about 2 ppm at the end of August to about 0.5 ppm by mid-October. The average decay in this region can be described as exponential with a half life of about 25 days. While total ozone, as obtained from profile integration, declined only about 35 percent, the integrated ozone between 14 and 18 km declined more than 70 percent. Vertical ozone profiles in the vortex revealed unusual structure with major features from 1 to 5 km thick which had suffered ozone depletions as great as 90 percent.

Hofmann, D. J.↗
NASA NTRS1991

Application of Krylov exponential propagation to fluid dynamics equations

This paper presents an application of matrix exponentiation via Krylov subspace projection, to the solution of fluid dynamics problems. The main idea is to approximate the operation exp(A)v by means of a projection-like process onto a Krylov subspace. This results in a computation of an exponential matrix vector product similar to the one above but of a much smaller size. Time integration schemes can then be devised to exploit this basic computational kernel. The motivation of this approach is to provide time-integration schemes that are essentially of an explicit nature but which have good stability properties.

Saad, Y.↗
NASA NTRS1991

Application of Krylov exponential propagation to fluid dynamics equations

An application of matrix exponentiation via Krylov subspace projection to the solution of fluid dynamics problems is presented. The main idea is to approximate the operation exp(A)v by means of a projection-like process onto a krylov subspace. This results in a computation of an exponential matrix vector product similar to the one above but of a much smaller size. Time integration schemes can then be devised to exploit this basic computational kernel. The motivation of this approach is to provide time-integration schemes that are essentially of an explicit nature but which have good stability properties.

Saad, Youcef↗
NASA NTRS1993

AGB stars as a source of short-lived radioactive nuclei in the solar nebula

The purpose is to estimate the possible contribution of some short-lived nuclei to the early solar nebula from asymptotic giant branch (AGB) sources. Low mass (1 to 3 solar mass) AGB stars appear to provide a site for synthesis of the main s process component for solar system material with an exponential distribution of neutron irradiations varies as exp(-tau/tau(sub 0)) (where tau is the time integrated neutron flux with a mean neutron exposure tau(sub 0)) for solar abundances with tau(sub 0) = 0.28 mb(sup -1). Previous workers estimated the synthesis of key short-lived nuclei which might be produced in AGB stars. While these calculations exhibit the basic characteristics of nuclei production by neutron exposure, there is need for a self-consistent calculation that follows AGB evolution and takes into account the net production from a star and dilution with the cloud medium. Many of the general approaches and the conclusions arrived at were presented earlier by Cameron. The production of nuclei for a star of 1.5 solar mass during the thermal pulsing of the AGB phase was evaluated. Calculations were done for a series of thermal pulses with tau(sub 0) = 0.12 and 0.28 mb(sup -1). These pulses involve s nucleosynthesis in the burning shell at the base of the He zone followed by the ignition of the H burning shell at the top of the He zone. After about 10-15 cycles the abundances of the various nuclei in the He zone become constant. Computations of the abundances of all nuclei in the He zone were made following Gallino. The mass of the solar nebula was considered to consist of some initial material of approximately solar composition plus some contributions from AGB stars. The ratios of the masses required from the AGB He burning zone to the ISM necessary to produce the observed value of Pd-107/Pd-108 in the early solar system were calculated and this dilution factor was applied to all other relevant nuclei.

Wasserburg, G. J.↗
DOE OSTI2020

Rossi-alpha Uncertainty Quantification by Analytic, Bootstrap, and Sample Methods to Inform Fitting Best Practices

The prompt neutron period (the negative reciprocal of the prompt neutron decay constant) can be estimated using the Rossi-alpha technique that is predicated on fitting Rossi alpha histograms and of interest in nuclear criticality safety and nonproliferation. The histograms are traditionally fit with a one-exponential model; however, recent work has proposed a two-exponential model to account for reflector induced phenomenon. Until recently, the uncertainty quantification for either model was inadequate (inaccurate and demanded large measurement times). Measurement uncertainty quantification by sample and analytic methods was developed and validated in Ref. The purpose of this transaction is to (i) validate a new bootstrap method by comparing bin-by-bin error bar estimates and (ii) demonstrate how to choose bin widths and reset times to optimize precision and accuracy.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Spectral form factor of a quantum spin glass

It is widely expected that systems which fully thermalize are chaotic in the sense of exhibiting random-matrix statistics of their energy level spacings, whereas integrable systems exhibit Poissonian statistics. In this paper, we investigate a third class: spin glasses. These systems are partially chaotic but do not achieve full thermalization due to large free energy barriers. We examine the level spacing statistics of a canonical infinite-range quantum spin glass, the quantum p-spherical model, using an analytic path integral approach. We find statistics consistent with a direct sum of independent random matrices, and show that the number of such matrices is equal to the number of distinct metastable configurations — the exponential of the spin glass “complexity” as obtained from the quantum Thouless-Anderson-Palmer equations. We also consider the statistical properties of the complexity itself and identify a set of contributions to the path integral which suggest a Poissonian distribution for the number of metastable configurations. Our results show that level spacing statistics can probe the ergodicity-breaking in quantum spin glasses and provide a way to generalize the notion of spin glass complexity beyond models with a semi-classical limit.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Code Generators for Floating-Point Unit Design in Integrated Circuits (OpenFloat) v1.0

This IP provides a comprehensive set of code generators for various floating-point units (FPUs) essential for integrated circuit design and integration, targeting a broad spectrum of applications, including machine learning and scientific computing. The suite includes FP adders, multipliers, subtractors, dividers, reciprocals, exponentials, square roots, trigonometric functions (sine, cosine, arctangent), and more. It supports customizable hardware design parameters, such as precision (16, 32, 64, and 128 bits) and pipeline depths, offering users enhanced flexibility and productivity. The generated code is in an industry-standard hardware description language, ensuring compatibility with standard design flows, including simulation, verification, synthesis, and implementation on both field-programmable gate arrays (FPGAs) and application-specific integrated circuits (ASICs).

Shalf, JohnM. [Lawrence Berkeley National Laborato↗
DOE OSTI2024

Simple and General Unitarity Conserving Numerical Real-Time Propagators of the Time-Dependent Schrödinger Equation Based on Magnus Expansion

Magnus expansion (ME) provides a general way to expand the real-time propagator of a time-dependent Hamiltonian within the exponential such that the unitarity is satisfied at any order. Here, we use this property and explicit integration of Lagrange interpolation formulas for the time-dependent Hamiltonian within each time interval and derive approximations that preserve unitarity for the differential time evolution operators of general time-dependent Hamiltonians. The resulting second-order approximation is the same as using the average of Hamiltonians for two end points of time. We identify three fourth-order approximations involving commutators of Hamiltonians at different times and also derive a sixth-order expression. A test of these approximations along with other available expressions for a two-state time-dependent Hamiltonian with sinusoidal time dependences provides information on the relative performance of these approximations and suggests that the derived expressions can serve as useful numerical tools for time evolution in time-resolved spectroscopy, quantum control, quantum sensing, real-time ab initio quantum dynamics, and open system quantum dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1972

Origin, evolution and present thermal state of the moon

The relative absence of lunar volcanism in the last 3 b.y. and the Apollo 15 heat flow measurement suggest that present-day temperatures in the moon are approximately steady state to depths of 100 km. An exponential distribution of heat sources with depth is scaled by equating the surface heat flow to the integrated heat production of this exterior shell. Presumed present-day interior temperatures and the present-day surface heat flow of 30 ergs/cm2-sec are obtained. 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. The present thermal state of the moon may involve at least some partial melting through all the lunar interior deeper than 200 km. Such a thermal configuration is inconsistent neither with temperatures inferred from electrical conductivity studies nor with the nonhydrostatic shape of the moon.

Hanks, T. C.↗
NASA NTRS1973

Hot H2 and interstellar shocks

The rotational excitation of H2 molecules by interstellar shock waves is calculated. Level populations integrated through a shock in a cloud with a density of 10 per cu cm have nearly exponential distributions, with effective rotational temperatures T sub r = 15 V to the 1.5 power K, where V is the shock velocity in km/sec. Ortho-para conversion is effective for high shock velocities, but not for low ones, so that the degree of alternation of intensities of rotational lines in the Lyman absorption band should be correlated with V, and therefore with T sub 2. The results are compared with observations of H2 rotational levels made by the Copernicus satellite.

Aannestad, P. A.↗
NASA NTRS1985

Tracking and disturbance rejection of MIMO nonlinear systems with PI controller

The tracking and disturbance rejection of a class of MIMO nonlinear systems with linear proportional plus integral (PI) compensator is studied. Roughly speaking, it is shown that if the given nonlinear plant is exponentially stable and has a strictly increasing dc steady-state I/O map, then a simple PI compensator can be used to yield a stable unity-feedback closed-loop system which asymptotically tracks reference inputs that tend to constant vectors and asymptotically rejects disturbances that tend to constant vectors.

Desoer, C. A.↗
NASA NTRS1985

Tracking and disturbance rejection of MIMO nonlinear systems with PI controller

The tracking and disturbance rejection of a class of MIMO nonlinear systems with a linear proportional plus integral (PI) compensator is studied. Roughly speaking, it is shown that if the given nonlinear plant is exponentially stable and has a strictly increasing dc steady-state I/O map, then a simple PI compensator can be used to yield a stable unity-feedback closed-loop system which asymptotically tracks reference inputs that tend to constant vectors and asymptotically rejects disturbances that tend to constant vectors.

Desoer, C. A.↗