Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Recursion”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Orthogonal Polynomials Defined by Self-Similar Measures with Overlaps

Here, we study orthogonal polynomials with respect to self-similar measures, focusing on the class of infinite Bernoulli convolutions, which are defined by iterated function systems with overlaps, especially those defined by the Pisot, Garsia, and Salem numbers. By using an algorithm of Mantica, we obtain graphs of the coefficients of the 3-term recursion relation defining the orthogonal polynomials. We use these graphs to predict whether the singular infinite Bernoulli convolutions belong to the Nevai class. Based on our numerical results, we conjecture that all infinite Bernoulli convolutions with contraction ratios greater than or equal to 1/2 belong to Nevai’s class, regardless of the probability weights assigned to the self-similar measures.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

From binaries to multiples. II. Hierarchical multiplicity of F and G dwarfs

Statistics of hierarchical multiplicity among solar-type dwarfs are studied using the distance-limited sample of 4847 targets presented in the accompanying Paper I. Known facts about binaries (multiplicity fraction 0.46, lognormal period distribution with median period 100 yr and logarithmic dispersion 2.4, and nearly uniform mass-ratio distribution independent of the period) are confirmed with a high statistical significance. The fraction of hierarchies with three or more components is 0.13 ± 0.01, and the fractions of targets with n = 1, 2, 3, ... components are 54:33:8:4:1. Subsystems in the secondary components are almost as frequent as in the primary components, but in half of such cases both inner pairs are present. The high frequency of those 2+2 hierarchies (4%) suggests that both inner pairs were formed by a common process. The statistics of hierarchies can be reproduced by simulations, assuming that the field is a mixture coming from binary-rich and binary-poor environments. Periods of the outer and inner binaries are selected recursively from the same lognormal distribution, subject to the stability constraint and accounting for the correlation between inner subsystems. The simulator can be used to evaluate the frequency of multiple systems with specified parameters. However, it does not reproduce the observed excess of inner periods shorter than 10 days, caused by tidal evolution.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Prolate spheroidal harmonic expansion of gravitational field

As a modification of the oblate spheroidal case, a recursive method is developed to compute the point value and a few low-order derivatives of the prolate spheroidal harmonics of the second kind, Q{sub nm} (y), namely the unnormalized associated Legendre function (ALF) of the second kind with its argument in the domain, 1 < y < ∞. They are required in evaluating the prolate spheroidal harmonic expansion of the gravitational field in addition to the point value and the low-order derivatives of P-bar {sub nm}(t), the 4π fully normalized ALF of the first kind with its argument in the domain, |t| ≤ 1. The new method will be useful in the gravitational field computation of elongated celestial objects.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2011

PHOTOIONIZATIONAL PLASMAS. I. THEORY

In this paper, an attempt is made to define the subject of photoionizational plasmas through rigorous mathematical formalism. The central results of this paper are the following. (1) A set of recursive equations is introduced for the computation of the charge state distributions in photoionizational plasmas. (2) Quantitative validity limits are given for both the collisional and the photoionizational domains. (3) A parameter that determines the charge state distribution in the photoionizational regime is introduced. (4) We provide detailed discussion about the most important components of the emission spectrum in the different equilibrium domains of the emitting plasma. Some of these components have not been reported so far as included in the analysis of such plasmas.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

PHYSICAL PARAMETERS OF STANDARD AND BLOWOUT JETS

The X-ray Telescope on board the Hinode mission revealed the occurrence, in polar coronal holes, of much more numerous jets than previously indicated by the Yohkoh/Soft X-ray Telescope. These plasma ejections can be of two types, depending on whether they fit the standard reconnection scenario for coronal jets or if they include a blowout-like eruption. In this work, we analyze two jets, one standard and one blowout, that have been observed by the Hinode and STEREO experiments. We aim to infer differences in the physical parameters that correspond to the different morphologies of the events. To this end, we adopt spectroscopic techniques and determine the profiles of the plasma temperature, density, and outflow speed versus time and position along the jets. The blowout jet has a higher outflow speed, a marginally higher temperature, and is rooted in a stronger magnetic field region than the standard event. Our data provide evidence for recursively occurring reconnection episodes within both the standard and the blowout jet, pointing either to bursty reconnection or to reconnection occurring at different locations over the jet lifetimes. We make a crude estimate of the energy budget of the two jets and show how energy is partitioned among different forms. Also, we show that the magnetic energy that feeds the blowout jet is a factor of 10 higher than the magnetic energy that fuels the standard event.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

A classical and a relativistic law of motion for spherical supernovae

In this paper we derive some first order differential equations which model the classical and the relativistic thin layer approximations. The circumstellar medium is assumed to follow a density profile of the Plummer type, the Lane-Emden (n = 5) type, or a power law. The first order differential equations are solved analytically, numerically, by a series expansion, or by recursion. The initial conditions are chosen in order to model the temporal evolution of SN 1993J over 10 yr and a smaller chi-squared is obtained for the Plummer case with η = 6. The stellar mass ejected by the SN progenitor prior to the explosion, expressed in solar mass, is identified with the total mass associated with the selected density profile and varies from 0.217 to 0.402 when the central number density is 10{sup 7} particles per cubic centimeter. The FWHM of the three density profiles, which can be identified with the size of the pre-SN 1993J envelope, varies from 0.0071 pc to 0.0092 pc.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

MAGNETOHYDRODYNAMIC MODELING OF SOLAR SYSTEM PROCESSES ON GEODESIC GRIDS

This report describes a new magnetohydrodynamic numerical model based on a hexagonal spherical geodesic grid. The model is designed to simulate astrophysical flows of partially ionized plasmas around a central compact object, such as a star or a planet with a magnetic field. The geodesic grid, produced by a recursive subdivision of a base platonic solid (an icosahedron), is free from control volume singularities inherent in spherical polar grids. Multiple populations of plasma and neutral particles, coupled via charge-exchange interactions, can be simulated simultaneously with this model. Our numerical scheme uses piecewise linear reconstruction on a surface of a sphere in a local two-dimensional 'Cartesian' frame. The code employs Haarten-Lax-van-Leer-type approximate Riemann solvers and includes facilities to control the divergence of the magnetic field and maintain pressure positivity. Several test solutions are discussed, including a problem of an interaction between the solar wind and the local interstellar medium, and a simulation of Earth's magnetosphere.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Linear-algebraic bath transformation for simulating complex open quantum systems

In studying open quantum systems, the environment is often approximated as a collection of non-interacting harmonic oscillators, a configuration also known as the star-bath model. It is also well known that the star-bath can be transformed into a nearest-neighbor interacting chain of oscillators. The chain-bath model has been widely used in renormalization group approaches. The transformation can be obtained by recursion relations or orthogonal polynomials. Based on a simple linear algebraic approach, we propose a bath partition strategy to reduce the system-bath coupling strength. As a result, the non-interacting star-bath is transformed into a set of weakly coupled multiple parallel chains. Furthermore, the transformed bath model allows complex problems to be practically implemented on quantum simulators, and it can also be employed in various numerical simulations of open quantum dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2012

Self-consistent initial conditions for primordial black hole formation

For an arbitrarily strong, spherically symmetric super-horizon curvature perturbation, we present analytic solutions of the Einstein equations in terms of the asymptotic expansion over the ratio of the Hubble radius to the length-scale of the curvature perturbation to set initial conditions for numerical computations of primordial black hole formation. To obtain this solution we develop a recursive method of quasi-linearization which reduces the problem to a system of coupled ordinary differential equations for the Nth order terms in the asymptotic expansion with sources consisting of a non-linear combination of the lower order terms.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Excursions through KK modes

In this article we study Kaluza-Klein (KK) dimensional reduction of massive Abelian gauge theories with charged matter fields on a circle. Since local gauge transformations change position dependence of the charged fields, the decomposition of the charged matter fields into KK modes is gauge dependent. While whole KK mass spectrum is independent of the gauge choice, the mode number depends on the gauge. The masses of the KK modes also depend on the field value of the zero-mode of the extra dimensional component of the gauge field. In particular, one of the KK modes in the KK tower of each massless 5D charged field becomes massless at particular values of the extra-dimensional component of the gauge field. When the extra-dimensional component of the gauge field is identified with the inflaton, this structure leads to recursive cosmological particle productions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

The matrix product approximation for the dynamic cavity method

Stochastic dynamics of classical degrees of freedom, defined on vertices of locally tree-like graphs, can be studied in the framework of the dynamic cavity method which is exact for tree graphs. Such models correspond for example to spin-glass systems, Boolean networks, neural networks, and other technical, biological, and social networks. The central objects in the cavity method are edge messages—conditional probabilities of two vertex variable trajectories. In this paper, we discuss a rather pedagogical derivation for the dynamic cavity method, give a detailed account of the novel matrix product edge message (MPEM) algorithm for the solution of the dynamic cavity equation as introduced, and present optimizations and extensions. Matrix product approximations of the edge messages are constructed recursively in an iteration over time. Computation costs and precision can be tuned by controlling the matrix dimensions of the MPEM in truncations. Without truncations, the dynamics is exact. Data for Glauber–Ising dynamics shows a linear growth of computation costs in time. In contrast to Monte Carlo simulations, the approach has a much better error scaling. Hence, it gives for example access to low probability events and decaying observables like temporal correlations. Here, we discuss optimized truncation schemes and an extension that allows to capture models which have a continuum time limit.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

A canonical Hamiltonian for open quantum systems

Abstract If an open quantum system is initially uncorrelated from its environment, then its dynamics can be written in terms of a Lindblad-form master equation. The master equation is divided into a unitary piece, represented by an effective Hamiltonian, and a dissipative piece, represented by a hermiticity-preserving superoperator; however, the division of open system dynamics into unitary and dissipative pieces is non-unique. For finite-dimensional quantum systems, we resolve this non-uniqueness by specifying a norm on the space of dissipative superoperators and defining the canonical Hamiltonian to be the one whose dissipator is minimal. We show that the canonical Hamiltonian thus defined is equivalent to the Hamiltonian initially defined by Lindblad, and that it is uniquely specified by requiring the dissipator’s jump operators to be traceless, extending a uniqueness result known previously in the special case of Markovian master equations. For a system weakly coupled to its environment, we give a recursive formula for computing the canonical effective Hamiltonian to arbitrary orders in perturbation theory, which we can think of as a perturbative scheme for renormalizing the system’s bare Hamiltonian.

Physics↗
DOE OSTI2026

Analog and symbolic computation through the Koopman framework

We develop a Koopman operator framework for studying the computational structure of dynamical systems. Specifically, we show that the resolvent of the Koopman operator provides a natural abstraction of halting, yielding a ‘Koopman halting problem’ that is recursively enumerable in general. For symbolic systems, such as those defined on Cantor space, this operator formulation captures reachability between clopen sets, while for equicontinuous systems we prove that the Koopman halting problem is decidable. Our framework demonstrates that absorbing (halting) states in coarse-grained finite automata correspond to Koopman eigenfunctions with eigenvalue one, while cycles in the transition graph impose spectral constraints associated with periodic dynamics. These results provide a unifying perspective on computation in symbolic and analog systems, showing how computational universality is reflected in operator spectra, invariant subspaces, and algebraic structures. Beyond symbolic dynamics, this operator-theoretic lens opens pathways to analyze the computational properties of a broader class of dynamical systems, including polynomial and analog models, and suggests that computational hardness may admit dynamical signatures in terms of Koopman spectral structure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Strategies for connecting whole-building LCA to the low-carbon design process

Abstract Decarbonization is essential to meeting urgent climate goals. With the building sector in the United States accounting for 35% of total U.S. carbon emissions, reducing environmental impacts within the built environment is critical. Whole-building life cycle analysis (WBLCA) quantifies the impacts of a building throughout its life cycle. Despite being a powerful tool, WBLCA is not standard practice in the integrated design process. When WBLCA is used, it is typically either speculative and based on early design information or conducted only after design completion as an accounting measure, with virtually no opportunity to impact the actual design. This work proposes a workflow for fully incorporating WBLCA into the building design process in an iterative, recursive manner, where design decisions impact the WBLCA, which in turn informs future design decisions. We use the example of a negative-operational carbon modular building seeking negative upfront embodied carbon using bio-based materials for carbon sequestration as a case study for demonstrating the utility of the framework. Key contributions of this work include a framework of computational processes for conducting iterative WBLCA, using a combination of an existing building WBLCA tool (Tally) within the building information modeling superstructure (Revit) and a custom script (in R) for materials, life cycle stages, and workflows not available in the WBLCA tool. Additionally, we provide strategies for harmonizing the environmental impacts of novel materials or processes from various life cycle inventory sources with materials or processes in existing building WBLCA tool repositories. These strategies are useful for those involved in building design with an interest in reducing their environmental impact. For example, this framework would be useful for researchers who are conducting WBLCAs on projects that include new or unusual materials and for design teams who want to integrate WBLCA more fully into their design process in order to ensure the building materials are consciously chosen to advance climate goals, while still ensuring best performance by traditional measures.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Linear complexity

We present factorization and solution phases for a new linear complexity direct solver designed for concurrent batch operations on fine-grained parallel architectures, for matrices amenable to hierarchical representation. We focus on the strong-admissibility-based $\mathscr{H}^{2}$ format, where strong recursive skeletonization factorization compresses remote interactions. We build upon previous implementations of $\mathscr{H}^{2}$ matrix construction for efficient factorization and solution algorithm design, which are illustrated graphically in stepwise detail. The algorithms are ‘blackbox’ in the sense that the only inputs are the matrix and right-hand side, without analytical or geometrical information about the origin of the system. We demonstrate linear complexity scaling in both time and memory on four representative families of dense matrices up to one million in size. Parallel scaling up to 16 threads is enabled by a multi-level matrix graph coloring and avoidance of dynamic memory allocations thanks to prefix-sum memory management. An experimental backward error analysis is included. We break down the timings of different phases, identify phases that are memory-bandwidth limited, and discuss alternatives for phases that may be sensitive to the trend to employ lower precisions for performance.

Boukaram, Wajih↗
DOE OSTI2018

On the Lebesgue constant of weighted Leja points for Lagrange interpolation on unbounded domains

Here, this work focuses on weighted Lagrange interpolation on an unbounded domain and analyzes the Lebesgue constant for a sequence of weighted Leja points. The standard Leja points are a nested sequence of points defined on a compact subset of the real line and can be extended to unbounded domains with the introduction of a weight function w:R → [0,1]. Due to a simple recursive formulation in one dimension, such abscissas provide a foundation for high-dimensional approximation methods such as sparse grid collocation, deterministic least squares and compressed sensing. Just as in the unweighted case of interpolation on a compact domain, we use results from potential theory to prove that the Lebesgue constant for the Leja points grows subexponentially with the number of interpolation nodes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

On decoupling the integrals of cosmological perturbation theory

ABSTRACT Perturbation theory (PT) is often used to model statistical observables capturing the translation and rotation-invariant information in cosmological density fields. PT produces higher order corrections by integration over linear statistics of the density fields weighted by kernels resulting from recursive solution of the fluid equations. These integrals quickly become high dimensional and naively require increasing computational resources the higher the order of the corrections. Here, we show how to decouple the integrands that often produce this issue, enabling PT corrections to be computed as a sum of products of independent 1D integrals. Our approach is related to a commonly used method for calculating multiloop Feynman integrals in quantum field theory, the Gegenbauer Polynomial x-Space Technique. We explicitly reduce the three terms entering the 2-loop power spectrum, formally requiring 9D integrations, to sums over successive 1D radial integrals. These 1D integrals can further be performed as convolutions, rendering the scaling of this method Nglog Ng with Ng the number of grid points used for each fast Fourier transform. This method should be highly enabling for upcoming large-scale structure redshift surveys where model predictions at an enormous number of cosmological parameter combinations will be required by Monte Carlo Markov Chain searches for the best-fitting values.

Slepian, Zachary↗
DOE OSTI2020

Combining the synergistic control capabilities of modeling and experiments: Illustration of finding a minimum-time quantum objective

A common way to manipulate a quantum system, for example spins or artificial atoms, is to use properly tailored control pulses. In order to accomplish quantum information tasks before coherence is lost, it is crucial to implement the control in the shortest possible time. Here we report the near time-optimal preparation of a Bell state with fidelity higher than 99% in an NMR experiment, which is feasible by combining the synergistic capabilities of modeling and experiments operating in tandem. In this work, the pulses preparing the Bell state are found by experiments that are recursively assisted with a gradient-based optimization algorithm working with a model. Thus, we exploit the interplay between model-based numerical optimal design and experimental-based learning control. Utilizing the balanced synergism between the dual approaches, as dictated by the case specific capabilities of each approach, should have broad applications for accelerating the search for optimal quantum controls.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗