SEARCH · Search NASA

Results for “asymptotic preserving”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

66 records · Page 4

All sources NASA NTRS DOE OSTI

NASA NTRS1994

Blocking reduction of Landsat Thematic Mapper JPEG browse images using optimal PSNR estimated spectra adaptive postfiltering

Two representative sample images of Band 4 of the Landsat Thematic Mapper are compressed with the JPEG algorithm at 8:1, 16:1 and 24:1 Compression Ratios for experimental browsing purposes. We then apply the Optimal PSNR Estimated Spectra Adaptive Postfiltering (ESAP) algorithm to reduce the DCT blocking distortion. ESAP reduces the blocking distortion while preserving most of the image's edge information by adaptively postfiltering the decoded image using the block's spectral information already obtainable from each block's DCT coefficients. The algorithm iteratively applied a one dimensional log-sigmoid weighting function to the separable interpolated local block estimated spectra of the decoded image until it converges to the optimal PSNR with respect to the original using a 2-D steepest ascent search. Convergence is obtained in a few iterations for integer parameters. The optimal logsig parameters are transmitted to the decoder as a negligible byte of overhead data. A unique maxima is guaranteed due to the 2-D asymptotic exponential overshoot shape of the surface generated by the algorithm. ESAP is based on a DFT analysis of the DCT basis functions. It is implemented with pixel-by-pixel spatially adaptive separable FIR postfilters. PSNR objective improvements between 0.4 to 0.8 dB are shown together with their corresponding optimal PSNR adaptive postfiltered images.

Linares, Irving↗

NASA NTRS1995

Solutions of the benchmark problems by the dispersion-relation-preserving scheme

The 7-point stencil Dispersion-Relation-Preserving scheme of Tam and Webb is used to solve all the six categories of the CAA benchmark problems. The purpose is to show that the scheme is capable of solving linear, as well as nonlinear aeroacoustics problems accurately. Nonlinearities, inevitably, lead to the generation of spurious short wave length numerical waves. Often, these spurious waves would overwhelm the entire numerical solution. In this work, the spurious waves are removed by the addition of artificial selective damping terms to the discretized equations. Category 3 problems are for testing radiation and outflow boundary conditions. In solving these problems, the radiation and outflow boundary conditions of Tam and Webb are used. These conditions are derived from the asymptotic solutions of the linearized Euler equations. Category 4 problems involved solid walls. Here, the wall boundary conditions for high-order schemes of Tam and Dong are employed. These conditions require the use of one ghost value per boundary point per physical boundary condition. In the second problem of this category, the governing equations, when written in cylindrical coordinates, are singular along the axis of the radial coordinate. The proper boundary conditions at the axis are derived by applying the limiting process of r approaches 0 to the governing equations. The Category 5 problem deals with the numerical noise issue. In the present approach, the time-independent mean flow solution is computed first. Once the residual drops to the machine noise level, the incident sound wave is turned on gradually. The solution is marched in time until a time-periodic state is reached. No exact solution is known for the Category 6 problem. Because of this, the problem is formulated in two totally different ways, first as a scattering problem then as a direct simulation problem. There is good agreement between the two numerical solutions. This offers confidence in the computed results. Both formulations are solved as initial value problems. As such, no Kutta condition is required at the trailing edge of the airfoil.

Tam, Christopher K. W.↗

DOE OSTI2026

Geometric decoherence time in Lindbladian dynamics

The onset of decoherence in open many-body systems lacks a dynamical timescale grounded in the loss of bipartite entanglement. Here, we introduce the geometric decoherence time, defined as the earliest moment the monotone relation between logarithmic negativity and Rényi-$\frac{1}{2}$ entropy—exactly equal across any bipartition for pure states—breaks down under open-system evolution, signaling entropy growth without accompanying entanglement growth. We establish this criterion in both single-particle Gaussian dynamics and many-body Lindbladian evolution. We show that quantum mutual information provides a complementary long-time diagnostic: Its asymptotic vanishing is equivalent to factorization of the steady state across the bipartition, a condition strictly stronger than separability, and whenever a product steady state is approached exponentially in trace norm, negativity and mutual information share the same decay rate. In the presence of a strong symmetry, this tracking can fail—residual classical correlations can survive after entanglement has vanished. In the Kitaev chain with balanced gain and loss, we derive a closed-form solution and show that the topological phase sustains longer coherence times than the trivial phase at identical dissipation, with a local minimum at the chiral-symmetric point. In the interacting XXZ chain, exact many-body evolution shows that local 𝑍 dephasing preserves residual classical correlations, whereas gain and loss restore the mutual-information tracking of negativity. Furthermore, our results establish the geometric decoherence time as a dynamical scale tracking the onset of decoherence.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2021

Theory of Trotter Error with Commutator Scaling

The Lie-Trotter formula, together with its higher-order generalizations, provides a simple approach to decomposing the exponential of a sum of operators. Despite significant effort, the error scaling of such product formulas remains poorly understood. We develop a theory of Trotter error that overcomes the limitations of truncating the Baker-Campbell-Hausdorff expansion. Our analysis directly exploits the commutativity of operator summands, producing tighter error bounds for both real- and imaginary-time evolutions. Whereas previous work achieves similar goals for systems with geometric locality or Lie-algebraic structure, our approach holds in general. We give a host of improved algorithms for digital quantum simulation and quantum Monte Carlo methods, nearly matching or even outperforming the best previous results. Our applications include: (i) a simulation of second-quantized plane-wave electronic structure, nearly matching the interaction-picture algorithm of Low and Wiebe; (ii) a simulation of $k$-local Hamiltonians almost with induced one-norm scaling, faster than the qubitization algorithm of Low and Chuang; (iii) a simulation of rapidly decaying power-law interactions, outperforming the Lieb-Robinson-based approach of Tran et al.; (iv) a hybrid simulation of clustered Hamiltonians, dramatically improving the result of Peng, Harrow, Ozols, and Wu; and (v) quantum Monte Carlo simulations of the transverse field Ising model and quantum ferromagnets, tightening previous analyses of Bravyi and Gosset. We obtain further speedups using the fact that product formulas can preserve the locality of the simulated system. Specifically, we show that local observables can be simulated with complexity independent of the system size for power-law interacting systems, which implies a Lieb-Robinson bound nearly matching a recent result of Tran et al. Our analysis reproduces known tight bounds for first- and second-order formulas. We further investigate the tightness of our bounds for higher-order formulas. For quantum simulation of a one-dimensional Heisenberg model with an even-odd ordering of terms, our result overestimates the complexity by only a factor of $5$. Our bound is also close to tight for power-law interactions and other orderings of terms. This suggests that our theory can accurately characterize Trotter error in terms of both the asymptotic scaling and the constant prefactor.

quantum computing, numerical analysis↗

NASA NTRS1997

Low-Dispersion Scheme for Nonlinear Acoustic Waves in Nonuniform Flow

The linear dispersion-relation-preserving scheme and its boundary conditions have been extended to the nonlinear Euler equations. This allowed computing, a nonuniform flowfield and a nonlinear acoustic wave propagation in such a medium, by the same scheme. By casting all the equations, boundary conditions, and the solution scheme in generalized curvilinear coordinates, the solutions were made possible for non-Cartesian domains and, for the better deployment of the grid points, nonuniform grid step sizes could be used. It has been tested for a number of simple initial-value and periodic-source problems. A simple demonstration of the difference between a linear and nonlinear propagation was conducted. The wall boundary condition, derived from the momentum equations and implemented through a pressure at a ghost point, and the radiation boundary condition, derived from the asymptotic solution to the Euler equations, have proven to be effective for the nonlinear equations and nonuniform flows. The nonreflective characteristic boundary conditions also have shown success but limited to the nonlinear waves in no mean flow, and failed for nonlinear waves in nonuniform flow.

Baysal, Oktay↗

DOE OSTI2022

Sharpening the Distance Conjecture in diverse dimensions

The Distance Conjecture holds that any infinite-distance limit in the scalar field moduli space of a consistent theory of quantum gravity must be accompanied by a tower of light particles whose masses scale exponentially with proper field distance ||$\phi$|| as m ~ exp(–λ||$\phi$||), where λ is order-one in Planck units. While the evidence for this conjecture is formidable, there is at present no consensus on which values of λ are allowed. In this paper, we propose a sharp lower bound for the lightest tower in a given infinite-distance limit in d dimensions: λ ≥ $1/\sqrt{d-2}$. In support of this proposal, we show that (1) it is exactly preserved under dimensional reduction, (2) it is saturated in many examples of string/M-theory compactifications, including maximal supergravity in d = 4 – 10 dimensions, and (3) it is saturated in many examples of minimal supergravity in d = 4 – 10 dimensions, assuming appropriate versions of the Weak Gravity Conjecture. We argue that towers with λ < $1/\sqrt{d-2}$ discussed previously in the literature are always accompanied by even lighter towers with λ ≥ $1/\sqrt{d-2}$, thereby satisfying our proposed bound. We discuss connections with and implications for the Emergent String Conjecture, the Scalar Weak Gravity Conjecture, the Repulsive Force Conjecture, large-field inflation, and scalar field potentials in quantum gravity. In particular, we argue that if our proposed bound applies beyond massless moduli spaces to scalar fields with potentials, then accelerated cosmological expansion cannot occur in asymptotic regimes of scalar field space in quantum gravity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Dispersion and the speed-limited particle-in-cell algorithm

This paper discusses temporally continuous and discrete forms of the speed-limited particle-in-cell (SLPIC) method first treated by Werner et al. [Phys. Plasmas 25, 123512 (2018)]. The dispersion relation for a 1D1V electrostatic plasma whose fast particles are speed-limited is derived and analyzed. By examining the normal modes of this dispersion relation, we show that the imposed speed-limiting substantially reduces the frequency of fast electron plasma oscillations while preserving the correct physics of lower-frequency plasma dynamics (e.g. ion acoustic wave dispersion and damping). We then demonstrate how the timestep constraints of conventional electrostatic particle-in-cell methods are relaxed by the speed-limiting approach, thus enabling larger timesteps and faster simulations. Here, these results indicate that the SLPIC method is a fast, accurate, and powerful technique for modeling plasmas wherein electron kinetic behavior is nontrivial (such that a fluid/Boltzmann representation for electrons is inadequate) but evolution is on ion timescales.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1993

Diffusion of Na and K in the uppermost regolith of Mercury

We examine the conditions under which ions impacting the Hermean surface can act as a regional source of enhanced atmospheric column through ion implantation and subsequent release at the surface. At most latitudes, energetic (a few keV) Na or K ions which impact the nightside surface, are released quickly (well before noon) upon warming of the surface. We show that the relative sunrise/sunset difference produced by ion implantation is ((n)(sub SR) - (n)(sub SS)/(n)(sub Ave) = f, where f is the fraction of the photo-ions recycled, (n)(sub SR) is the average zenith column above the sunrise portion of the illuminated hemisphere, (n)(sub SS) is the average zenith column in the sunset portion, and (n)(sub Ave) is the average column over the sunlit disk. Thus, to produce a large sunrise/sunset difference via ion implantation and subsequent release requires efficient (close to total) recycling. We show that the most extensive set of available data reduced to Na column abundance does not show any sunrise enhancement. We argue that the K data do not permit an unambiguous interpretation in favor of sunrise/sunset differences. We further find that if an efficient surface loss process for the alkali is not operating after sunrise, the initial Na and K distribution will relax into the bulk of the solid. Preserved abundance gradients of Na and K in lunar glasses that the Arrhenius coefficients for impact glasses are likely more modest than those for laboratory glasses derived from rock by nonimpact processes, but the effect of the more modest diffusion rates is only to delay the efficient loss of Na by a few Earth days. We argue that implantation can lead to observable regional increases in the observed Na or K column densities only if it occurs at very high latitudes, where diffusion is slow. It is typically lower energy ions which impact at high latitude and these are both more numerous than the high energy ions and possess smaller average penetration depths; thus there are additional reasons to favor a high-latitude locus for any possible prompt return related increases in zenith column. We find that we cannot rule out sputtering as a source process, as the low-energy ions are efficient sputterers. Finally, we argue that the observed Na/K ratio in the atmopshere may be the expression of their different asymptotic rates of loss from the interiors of the regolith grains. If so, this allows us to fix the importance of sputtering and photon stimulated desorption relative to impact vaporization.

Killen, R. M.↗

DOE OSTI2022

Ladder symmetries of black holes. Implications for love numbers and no-hair theorems

It is well known that asymptotically flat black holes in general relativity have a vanishing static, conservative tidal response. Here, we show that this is a result of linearly realized symmetries governing static (spin 0,1,2) perturbations around black holes. The symmetries have a geometric origin: in the scalar case, they arise from the (E)AdS isometries of a dimensionally reduced black hole spacetime. Underlying the symmetries is a ladder structure which can be used to construct the full tower of solutions, and derive their general properties: (1) solutions that decay with radius spontaneously break the symmetries, and must diverge at the horizon; (2) solutions regular at the horizon respect the symmetries, and take the form of a finite polynomial that grows with radius. Taken together, these two properties imply that static response coefficients — and in particular Love numbers — vanish. Moreover, property (1) is consistent with the absence of black holes with linear (perturbative) hair. We also discuss the manifestation of these symmetries in the effective point particle description of a black hole, showing explicitly that for scalar probes the worldline couplings associated with a non-trivial tidal response and scalar hair must vanish in order for the symmetries to be preserved.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Distributed State Estimation Over Time-Varying Graphs: Exploiting the Age-of-Information

Here, we study the problem of designing a distributed observer for an LTI system over a time-varying communication graph. The limited existing work on this topic imposes various restrictions either on the observation model or on the sequence of communication graphs. In contrast, we propose a single-time-scale distributed observer that works under mild assumptions. Specifically, our communication model only requires strong-connectivity to be preserved over non-overlapping, contiguous intervals that are even allowed to grow unbounded over time. We show that under suitable conditions that bound the growth of such intervals, joint observability is sufficient to track the state of any discrete-time LTI system exponentially fast, at any desired rate. We also develop a variant of our algorithm that is provably robust to worst-case adversarial attacks, provided the sequence of graphs is sufficiently connected over time. The key to our approach is the notion of a "freshness-index" that keeps track of the age-of-information being diffused across the network. Such indices enable nodes to reject stale estimates of the state, and, in turn, contribute to stability of the error dynamics.

42 ENGINEERING↗

DOE OSTI2020

Overset meshes for incompressible flows: On preserving accuracy of underlying discretizations

This study on overset meshes for incompressible-flow simulations is motivated by accurate prediction of wind farm aerodynamics involving large motions and deformations of components with complex geometry. Specifically, using first-order hyperbolic and elliptic equation proxies for the incompressible Navier-Stokes (NS) equations, we investigate the influence of information exchange between overset meshes on numerical performance where the underlying discretization is second-order accurate. The first aspect of information exchange surrounds interpolation of solution where we examine Lagrange and point-cloud-based interpolation for creating constraint equations between overset meshes. To maintain overall second-order accuracy, higher-order interpolation is required for elliptic problems, but linear interpolation is sufficient for hyperbolic problems in first-order form. Higher-order point-cloud-based interpolation provides a pathway to maintaining accuracy in unstructured meshes, but at higher complexity. The second aspect of information exchange focuses on comparing the approaches of overset single system (OSS) and overset Additive Schwarz (OAS) for coupling the linear systems of the overlapping meshes. While the former involves a single linear system, in the latter the discrete linear systems are solved separately, and solving the global system is accomplished through outer iterations and sequential information exchange in a Jacobi fashion. For the test cases studied, accuracy for hyperbolic problems is maintained by performing two outer iterations, whereas many outer iterations are required for elliptic systems. The order-of-accuracy studies explored here are critical for verifying the overset-mesh coupling algorithms used in engineering simulations. Accuracy of these simulations themselves is, however, quantified using engineering quantities of interest such as drag, power, etc. Consequently, we conclude with numerical experiments using NS equations for incompressible flows where we show that linear interpolation and few outer iterations are sufficient for achieving asymptotic convergence of engineering quantities of interest.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2020

An Evolutionary System of Mineralogy. Part I: Stellar Mineralogy (>13 to 4.6 Ga)

Minerals preserve records of the physical, chemical, and biological histories of their origins and subsequent alteration, and thus provide a vivid narrative of the evolution of Earth and other worlds through billions of years of cosmic history. Mineral properties, including trace and minor elements, ratios of isotopes, solid and fluid inclusions, external morphologies, and other idiosyncratic attributes, represent information that points to specific modes of formation and subsequent environmental histories—information essential to understanding the co-evolving geosphere and biosphere. This perspective suggests an opportunity to amplify the existing system of mineral classification, by which minerals are defined solely on idealized end-member chemical compositions and crystal structures. Here we present the first in a series of contributions to explore a complementary evolutionary system of mineralogy—a classification scheme that links mineral species to their paragenetic modes. The earliest stage of mineral evolution commenced with the appearance of the first crystals in the universe at >13 Ga and continues today in the expanding, cooling atmospheres of countless evolved stars, which host the high-temperature (T > 1000 K), low-pressure (P < 10-2 atm) condensation of refractory minerals and amorphous phases. Most stardust is thought to originate in three distinct processes in carbon- and/or oxygen-rich mineral-forming stars: (1) condensation in the cooling, expanding atmospheres of asymptotic giant branch stars; (2) during the catastrophic explosions of supernovae, most commonly core collapse (Type II) supernovae; and (3) classical novae explosions, the consequence of runaway fusion reactions at the surface of a binary white dwarf star. Each stellar environment imparts distinctive isotopic and trace element signatures to the micro- and nanoscale stardust grains that are recovered from meteorites and micrometeorites collected on Earth’s surface, by atmospheric sampling, and from asteroids and comets. Although our understanding of the diverse mineral-forming environments of stars is as yet incomplete, we present a preliminary catalog of 41 distinct natural kinds of stellar minerals, representing 22 official International Mineralogical Association (IMA) mineral species, as well as 2 as yet unapproved crystalline phases and 3 kinds of non-crystalline condensed phases not codified by the IMA.

Robert M. Hazen↗