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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A Low-Rank QTT-based Finite Element Method for Elasticity Problems

We present an efficient and robust numerical algorithm for solving the linear elasticity problem that combines the Quantized Tensor Train format and a domain partitioning strategy. This approach makes it possible to solve the linear elasticity problem on a computational domain that is more general than a square. By integrating Z-ordering and subdomain concatenation, our method substantially decreases memory usage and achieves a notable reduction in rank compared to established Finite Element implementations like the FEniCS platform. This efficiency is maintained while still guaranteeing exponential convergence with respect to the number of degrees of freedom. This performance gain, however, requires a fundamental rethinking of how core finite element operations are implemented. This includes changes to mesh discretization, node and degree of freedom ordering, stiffness matrix and internal nodal force assembly, and the execution of algebraic matrix-vector operations. In this work, we discuss all these aspects in detail and assess the method’s performance in the numerical approximation of three representative test cases.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1986

The radial distribution of galactic gamma rays. III - The distribution of cosmic rays in the Galaxy and the CO-H2 calibration

The cosmic ray (CR) distribution in the Galaxy inside the solar circle is determined by using the kinematics of both H I and CO to ascertain the spatial distribution of the interstellar gas. The conversion factor between integrated CO line intensity and the molecular hydrogen column density is found, and possible large-scale variations of this relationship throughout the Galaxy are investigated. The galactocentric CR distributions are shown to be described satisfactorily by exponential distributions for radii larger than about three, with a radial scale length of 4-11 kpc for electrons and a scale length greater than 18 kpc for nuclei.

Bloemen, J. B. G. M.↗
DOE OSTI2024

Hidden quantum criticality and entanglement in quench dynamics

Entanglement exhibits universal behavior near the ground-state critical point where correlations are long ranged and the thermodynamic entropy is vanishing. On the other hand, a quantum quench imparts extensive energy and results in a build up of entropy, hence no critical behavior is expected at long times. In this work, we present a new paradigm in the quench dynamics of integrable spin chains which exhibit a ground-state order-disorder phase transition at a critical line. Specifically, we consider a quench along the critical line which displays a volume-law behavior of the entropy and exponentially decaying correlations; however, we show that quantum criticality is hidden in higher-order correlations and becomes manifest via measures such as the mutual information and logarithmic negativity. Furthermore, we showcase the scale invariance of the Rényi mutual information between disjoint regions as further evidence for genuine critical behavior. We attribute the emergent quantum criticality to the soft mode not getting excited in spite of the quench. Moreover, the results presented here are universal to models whose low-energy or long-wavelength dynamics are well described by a free-fermionic field theory. Our results are amenable to an experimental realization on different quantum simulator platforms, particularly the Rydberg simulators. Published by the American Physical Society 2024

Paul, Sanku↗
NASA NTRS2003

Flutter and Divergence Analysis using the Generalized Aeroelastic Analysis Method

The Generalized Aeroelastic Analysis Method (GAAM) is applied to the analysis of three well-studied checkcases: restrained and unrestrained airfoil models, and a wing model. An eigenvalue iteration procedure is used for converging upon roots of the complex stability matrix. For the airfoil models, exact root loci are given which clearly illustrate the nature of the flutter and divergence instabilities. The singularities involved are enumerated, including an additional pole at the origin for the unrestrained airfoil case and the emergence of an additional pole on the positive real axis at the divergence speed for the restrained airfoil case. Inconsistencies and differences among published aeroelastic root loci and the new, exact results are discussed and resolved. The generalization of a Doublet Lattice Method computer code is described and the code is applied to the calculation of root loci for the wing model for incompressible and for subsonic flow conditions. The error introduced in the reduction of the singular integral equation underlying the unsteady lifting surface theory to a linear algebraic equation is discussed. Acknowledging this inherent error, the solutions of the algebraic equation by GAAM are termed 'exact.' The singularities of the problem are discussed and exponential series approximations used in the evaluation of the kernel function shown to introduce a dense collection of poles and zeroes on the negative real axis. Again, inconsistencies and differences among published aeroelastic root loci and the new 'exact' results are discussed and resolved. In all cases, aeroelastic flutter and divergence speeds and frequencies are in good agreement with published results. The GAAM solution procedure allows complete control over Mach number, velocity, density, and complex frequency. Thus all points on the computed root loci can be matched-point, consistent solutions without recourse to complex mode tracking logic or dataset interpolation, as in the k and p-k solution methods.

Edwards, John W.↗
DOE OSTI2025

Open string amplitudes: singularities, asymptotics and new representations

Open string amplitudes at tree level have been studied for over fifty years. However, there is no known analytic form for general n-point amplitudes, and their conventional representation in terms of worldsheet integrals does not make many of their most basic physical properties manifest. Recently, a formulation of these amplitudes exposing the underlying “binary geometry” via the use of “u” variables has given us many insights into their basic features. In this paper, we initiate a systematic exploration of fundamental aspects of open string amplitudes from this new point of view. We begin by finding explicit expressions for the factorization of amplitudes at general massive levels, which are seen to be determined by products of lower-point massless amplitudes with shifted kinematics. We then study the asymptotic behavior when subsets of kinematic variables become large, delineating regimes with exponential (generalized hard scattering) and power-law (generalized Regge) behavior. We also give precise expressions for the asymptotics, which reveal another example of the recently observed property of factorization away from poles. We derive new recursion relations for the amplitude, which when repeatedly applied reduce to infinite series representations with a wider domain of convergence than the usual integral representations. For the five-point case, we present a new closed-form expression for the amplitude that for the first time gives its analytic continuation to all of kinematic space. We also discuss novel relations between amplitudes at different kinematic points following from the recently observed “split” factorizations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1980

Reformulation of Possio's kernel with application to unsteady wind tunnel interference

An efficient method for computing the Possio kernel has remained elusive up to the present time. In this paper the Possio is reformulated so that it can be computed accurately using existing high precision numerical quadrature techniques. Convergence to the correct values is demonstrated and optimization of the integration procedures is discussed. Since more general kernels such as those associated with unsteady flows in ventilated wind tunnels are analytic perturbations of the Possio free air kernel, a more accurate evaluation of their collocation matrices results with an exponential improvement in convergence. An application to predicting frequency response of an airfoil-trailing edge control system in a wind tunnel compared with that in free air is given showing strong interference effects.

Fromme, J. A.↗
NASA NTRS1993

Measurement of p-mode energy propagation in the quiet solar photosphere

We have measured and analyzed the p-mode oscillations in the profile of the Mg I 4571 A line in a quiet region near disk center. The oscillations are found to be mostly standing waves, in agreement with previous work. However, a small propagating component is measured, and we determine the direction, magnitude, and vertical variation of the energy propagation. The work integral indicates an upward energy flow of about 2 x 10 exp 7 ergs/sq cm/s at a height of 50 km above the base of the photosphere for waves with frequencies of 2-16 mHz. This energy flow decreases exponentially with height and drops below 10 exp 5 ergs/sq cm/s in the uppermost photosphere. The energy flow leaving the upper photosphere is at least an order of magnitude too small to constitute a significant source of heating for the chromosphere. However, the p-mode damping in the lower photosphere approaches levels large enough to account for the measured p-mode line widths. The relative amplitudes and phases of the thermodynamic quantities indicate that the p-mode are neither adiabatic nor isothermal in the photosphere.

Fontenla, J. M.↗
NASA NTRS1994

A spectral boundary integral equation method for the 2-D Helmholtz equation

In this paper, we present a new numerical formulation of solving the boundary integral equations reformulated from the Helmholtz equation. The boundaries of the problems are assumed to be smooth closed contours. The solution on the boundary is treated as a periodic function, which is in turn approximated by a truncated Fourier series. A Fourier collocation method is followed in which the boundary integral equation is transformed into a system of algebraic equations. It is shown that in order to achieve spectral accuracy for the numerical formulation, the nonsmoothness of the integral kernels, associated with the Helmholtz equation, must be carefully removed. The emphasis of the paper is on investigating the essential elements of removing the nonsmoothness of the integral kernels in the spectral implementation. The present method is robust for a general boundary contour. Aspects of efficient implementation of the method using FFT are also discussed. A numerical example of wave scattering is given in which the exponential accuracy of the present numerical method is demonstrated.

Hu, Fang Q.↗
NASA NTRS1975

Surface photometry of galaxies. I - The SB0 NGC 2950

Surface photometry of NGC 2950, a classic example of the SB0 type galaxy, is presented. The analysis shows that except for the rudimentary bar the galaxy closely resembles an S0 galaxy. The disk component has an exponential light curve with a central surface brightness of B(0) = 21.65 per square second of arc and a scale size of 4.3 kpc. The nucleus, bar, and disk components contain 29, 22, and 48 percent, respectively, of the light. Integral and differential colors are presented and discussed in terms of the three components mentioned above.

Crane, P.↗
NASA NTRS2021

Landsat Derived Bathymetry of Lakes on the Arctic Coastal Plan of Northern Alaska

The Pleistocene sand sea on the Arctic Coastal Plain (ACP) of northern Alaska is underlain by anancient sand dune field, a geological feature that affects regional lake characteristics. Many of these lakes, whichcover approximately 20 % of the Pleistocene sand sea, are relatively deep (up to 25 m). In addition to the nat-ural importance of ACP sand sea lakes for water storage, energy balance, and ecological habitat, the need forwinter water for industrial development and exploration activities makes lakes in this region a valuable resource.However, ACP sand sea lakes have received little prior study. Here, we collect in situ bathymetric data to test12 model variants for predicting sand sea lake depth based on analysis of Landsat-8 Operational Land Imager(OLI) images. Lake depth gradients were measured at 17 lakes in midsummer 2017 using a Humminbird 798ciHD SI Combo automatic sonar system. The field-measured data points were compared to red–green–blue (RGB)bands of a Landsat-8 OLI image acquired on 8 August 2016 to select and calibrate the most accurate spectral-depth model for each study lake and map bathymetry. Exponential functions using a simple band ratio (withbands selected based on lake turbidity and bed substrate) yielded the most successful model variants. For eachlake, the most accurate model explained 81.8 % of the variation in depth, on average. Modeled lake bathymetrieswere integrated with remotely sensed lake surface area to quantify lake water storage volumes, which rangedfrom 1.056×10−3to 57.416×10−3km3. Due to variations in depth maxima, substrate, and turbidity betweenlakes, a regional model is currently infeasible, rendering necessary the acquisition of additional in situ datawith which to develop a regional model solution. Estimating lake water volumes using remote sensing will fa-cilitate better management of expanding development activities and serve as a baseline by which to evaluatefuture responses to ongoing and rapid climate change in the Arctic. All sonar depth data and modeled lakebathymetry rasters can be freely accessed at https://doi.org/10.18739/A2SN01440 (Simpson and Arp, 2018) andhttps://doi.org/10.18739/A2HT2GC6G (Simpson, 2019), respectively.

Claire E Simpson↗
NASA NTRS2013

Measurements of Cyclotron Features and Pulse Periods in the High-Mass X-Ray Binaries 4U 1538-522 and 4U 1907+09 with the International Gamma-Ray Astrophysics Laboratory

We present a spectral and timing analysis of International Gamma-Ray Astrophysics Laboratory (INTEGRAL) observations of two high-mass X-ray binaries, 4U 1538−522 and 4U 1907+09. Our timing measurements for 4U 1538−522 find the pulse period to have exhibited a spin-up trend until approximately 2009, after which there is evidence for a torque reversal, with the source beginning to spin down to the most recently measured period of 525.407 plus or minus 0.001 seconds. The most recent INTEGRAL observations of 4U 1907+09 are not found to yield statistically significant pulse periods due to the significantly lower flux from the source compared with 4U 1538−522. A spectral model consisting of a power-law continuum with an exponential cutoff and modified by two cyclotron resonance scattering features is found to fit both sources well, with the cyclotron scattering features detected at approximately 22 and approximately 49 kiloelectronvolts for 4U 1538−522 and at approximately 18 and approximately 36 kiloelectronvolts for 4U 1907+09. The spectral parameters of 4U 1538−522 are generally not found to vary significantly with flux and there is little to no variation across the torque reversal. Examining our results in conjunction with previous work, we find no evidence for a correlation between cyclotron line energy and luminosity for 4U 1538−522. 4U 1907+09 shows evidence for a positive correlation between cyclotron line energy and luminosity, which would make it the fourth, and lowest luminosity, cyclotron line source to exhibit this relationship.

stars: magnetic field↗
NASA NTRS2005

The High Energy Spectrum of NGC 4151

We present first INTEGRAL observations of the type 1.5 Seyfert galaxy NGC 4151. Combining several INTEGRAL observations performed during 2003, totaling approximately 400 ksec of exposure time, allow us to study the spectrum in the 3 - 300 keV range. The measurements presented here reveal an overall spectrum from X-rays up to the soft gamma-rays that can be described by an absorbed (N(sub H) approximately equal to 5 x 10(exp 22) per square centimeter) and non-variable thermal component, plus a Fe Kalpha line, and an exponential cutoff occurs at 110 keV, consistent with earlier claims. The Galactic hydrogen column density in the line of sight is N(sub H), Gal approximately equal to 2.1 x 10 (exp 20) per square centimeter. The time resolved analysis shows little variation of the spectral parameters. The comparison with CGRO/OSSE data shows that the same spectral model can be applied over a time span of 15 years, while the flux varied by a factor of 2. Applying a Compton reflection component improves the model fit to the INTEGRAL data. Nonetheless the data available to date cannot significantly confirm or exclude the existence of reflection, nor is a high iron overabundance in the absorber, as had been previously suggested, clearly detectable.

Beckmann, V.↗
DOE OSTI2020

Distributed Resources for the Earth System Grid Federation (ESGF) Advanced Management (DREAM). Final Report

Distributed Resources for the Earth System Grid Federation (ESGF) Advanced Management (DREAM) is a proposed system that will enable data from an infinite number of diverse sources to be organized and accessed from anywhere using any handheld or other computer device. The approach offers a powerful roadmap for the creation and integration of a unified knowledge base of an entire ecosystem, including its many geophysical, geographical, social, political, agricultural, energy, transportation, and cyber aspects. The resulting aggregation of data has the potential to generate an informational universe of unprecedented size that has never before been possible due to the prohibitive costs, managerial complexity, and technical barriers associated with ever-changing exponential-growth data flows. We envision that DREAM will accelerate discovery by enabling climate researchers, among other types of researchers, to manage, analyze, and visualize data from earth-scale measurements and simulations. DREAM’s success will be built on proven components that leverage existing services and resources. A key building block for DREAM will be the ESGF, chaired by Dean N. Williams. Expanding on the existing ESGF, the project will ensure that the access, storage, movement, and analysis of the large quantities of data that are processed and produced by diverse science projects can be dynamically distributed with proper resource management. Much of the Office of Science data is currently generated by multiple stand-alone facilities. DREAM can collect data accumulated from these facilities and incorporate it into a fully integrated network accessible from anywhere in the world. The result is a completely new paradigm shift for data management, analysis, and visualization enabling researchers to: Manage their calculations, data, tools, and research results; Ensure that all data are sharable, reproducible and (re)usable—accompanied by appropriate metadata describing its provenance, syntax, and semantics at creation; Advance application performance by selectively adapting APIs and services in response to scientific requirements and architectural complexities; and Provide scalable interactive resource management—navigate data and metadata at multiple levels, provide architecture-aware data integration, analysis and visualization tools. We will engage closely with DOE, NASA, and NOAA science groups working at the leading edge of computing. These engagements—in domains such as biology, climate, and hydrology—will allow us to advance disciplinary science goals and inform our development of technologies that can accelerate discovery across DOE more broadly. We will advertise and promote our technologies via dedicated workshops, tutorials, and sessions at conferences, stand-alone events with broad inter-disciplinary invitation, and engagements with leadership facilities.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Lefschetz thimble quantum Monte Carlo for spin systems

Monte Carlo simulations are useful tools for modeling quantum systems, but in some cases they suffer from a sign problem, leading to an exponential slow down in their convergence to a value. While solving the sign problem is generically NP hard, many techniques exist for mitigating the sign problem in specific cases; in particular, the technique of deforming the Monte Carlo simulation's plane of integration onto Lefschetz thimbles (complex hypersurfaces of stationary phase) has seen significant success in the context of quantum field theories. We extend this methodology to spin systems by utilizing spin coherent state path integrals to reexpress the spin system's partition function in terms of continuous variables. Using some toy systems, we demonstrate its effectiveness at lessening the sign problem in this setting, despite the fact that the initial mapping to spin coherent states introduces its own sign problem. The standard formulation of the spin coherent path integral is known to make use of uncontrolled approximations; despite this, for large spins they are typically considered to yield accurate results, so it is somewhat surprising that our results show significant systematic errors. Furthermore, possibly of independent interest, our use of Lefschetz thimbles to overcome the intrinsic sign problem in spin coherent state path integral Monte Carlo enables a novel numerical demonstration of a breakdown in the spin coherent path integral.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Carbon Utilization: Electrochemical Approach Using Novel Catalyst and System Integration

The environmental footprint from power generation can be mitigated by capturing emitted carbon dioxide (CO2) and converting it to value-added fuel sources. A multitude of different fuels can be produced from CO2 depending on the chemical reaction catalysts. Typically, these pathways can be developed electrochemically, or through hydrogenation. In recent years the electrochemical reduction of CO2 has grown exponentially, as it has shown potential to decrease atmospheric pollution and provide an additional means of on-demand electricity storage. Formic acid (FA) in particular has industrial applications such as in fuel cells or as a hydrogen storage medium, and requires a lower thermodynamic energy input than any other fuel produced for electrochemical CO2 reduction. However, product selectivity remains a challenge for conventional electrocatalysts, making FA production efficiencies above 50% difficult. To address the challenge of catalytic selectivity, UK CAER employs enzymatic catalysts to facilitate the FA production and minimize unwanted byproducts. Our process couples the catalyst with a charge mediator which is necessary to shuttle electrons, and an electrolyte pH buffer to prevent the catalyst from degrading. Presently, production efficiencies have approached 60% at a rate of over 10 mM FA/hour. System design considerations (flow rate, bulk volume) will also be discussed.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2025

Genesis: A Compiler Framework for Hamiltonian Simulation on Hybrid CV-DV Quantum Computers

We introduce Genesis, the first compiler designed to support Hamiltonian Simulation on hybrid continuous-variable (CV) and discrete-variable (DV) quantum computing systems. Genesis is a two-level compilation system. At the first level, it decomposes an input Hamiltonian into basis gates using the native instruction set of the target hybrid CV-DV quantum computer. At the second level, it tackles the mapping and routing of qumodes/qubits to implement long-range interactions for the gates decomposed from the first level. Rather than a typical implementation that relies on SWAP primitives similar to qubit-based (or DV-only) systems, we propose an integrated design of connectivity-aware gate synthesis and beamsplitter SWAP insertion tailored for hybrid CV-DV systems. We also introduce an OpenQASM-like domain-specific language (DSL) named CVDV-QASM to represent Hamiltonian in terms of Pauli-exponentials and basic gate sequences from the hybrid CVDV gate set. Genesis has successfully compiled several important Hamiltonians, including the Bose-Hubbard model, Z2−Higgs model, Hubbard-Holstein model, Heisenberg model and Electron-vibration coupling Hamiltonians, which are critical in domains like quantum field theory, condensed matter physics, and quantum chemistry. Our implementation is available at Genesis-CVDV-Compiler https://github.com/ruadapt/Genesis-CVDV-Compiler

Chen, Henry↗
NASA NTRS2012

The PROPEL Electrodynamic Tether Mission and Connecting to the Ionosphere

The exponential increase of launch system size.and cost.with delta-V makes missions that require large total impulse cost prohibitive. Led by NASA's Marshall Space Flight Center, a team from government, industry, and academia has developed a flight demonstration mission concept of an integrated electrodynamic (ED) tethered satellite system called PROPEL: "Propulsion using Electrodynamics". The PROPEL Mission is focused on demonstrating a versatile configuration of an ED tether to overcome the limitations of the rocket equation, enable new classes of missions currently unaffordable or infeasible, and significantly advance the Technology Readiness Level (TRL) to an operational level. We are also focused on establishing a far deeper understanding of critical processes and technologies to be able to scale and improve tether systems in the future. Here, we provide an overview of the proposed PROPEL mission. One of the critical processes for efficient ED tether operation is the ability to inject current to and collect current from the ionosphere. Because the PROPEL mission is planned to have both boost and deboost capability using a single tether, the tether current must be capable of flowing in both directions and at levels well over 1 A. Given the greater mobility of electrons over that of ions, this generally requires that both ends of the ED tether system can both collect and emit electrons. For example, hollow cathode plasma contactors (HCPCs) generally are viewed as state-of-the-art and high TRL devices; however, for ED tether applications important questions remain of how efficiently they can operate as both electron collectors and emitters. Other technologies will be highlighted that are being investigated as possible alternatives to the HCPC such as Solex that generates a plasma cloud from a solid material (Teflon) and electron emission (only) technologies such as cold-cathode electron field emission or photo-electron beam generation (PEBG) techniques.

Gilchrist, Brian↗