Search NASA⌕ Search

SEARCH · Search NASA

Results for “local solution bounds”

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 73 records · Page 4

Smart-Divert Powered Descent Guidance to Avoid the Backshell Landing Dispersion Ellipse

A smart-divert capability has been added into the Powered Descent Guidance (PDG) software originally developed for Mars pinpoint and precision landing. The smart-divert algorithm accounts for the landing dispersions of the entry backshell, which separates from the lander vehicle at the end of the parachute descent phase and prior to powered descent. The smart-divert PDG algorithm utilizes the onboard fuel and vehicle thrust vectoring to mitigate landing error in an intelligent way: ensuring that the lander touches down with minimum- fuel usage at the minimum distance from the desired landing location that also avoids impact by the descending backshell. The smart-divert PDG software implements a computationally efficient, convex formulation of the powered-descent guidance problem to provide pinpoint or precision-landing guidance solutions that are fuel-optimal and satisfy physical thrust bound and pointing constraints, as well as position and speed constraints. The initial smart-divert implementation enforced a lateral-divert corridor parallel to the ground velocity vector; this was based on guidance requirements for MSL (Mars Science Laboratory) landings. This initial method was overly conservative since the divert corridor was infinite in the down-range direction despite the backshell landing inside a calculable dispersion ellipse. Basing the divert constraint instead on a local tangent to the backshell dispersion ellipse in the direction of the desired landing site provides a far less conservative constraint. The resulting enhanced smart-divert PDG algorithm avoids impact with the descending backshell and has reduced conservatism.

Carson, John M.↗

nPINNs: Nonlocal physics-informed neural networks for a parametrized nonlocal universal Laplacian operator. Algorithms and applications

Physics-informed neural networks (PINNs) are effective in solving inverse problems based on differential and integro-differential equations with sparse, noisy, unstructured, and multifidelity data. PINNs incorporate all available information, including governing equations (reflecting physical laws), initial-boundary conditions, and observations of quantities of interest, into a loss function to be minimized, thus recasting the original problem into an optimization problem. In this paper, we extend PINNs to parameter and function inference for integral equations such as nonlocal Poisson and nonlocal turbulence models, and we refer to them as nonlocal PINNs (nPINNs). The contribution of the paper is three-fold. First, we propose a unified nonlocal Laplace operator, which converges to the classical Laplacian as one of the operator parameters, the nonlocal interaction radius $\delta$ goes to zero, and to the fractional Laplacian as $\delta$ goes to infinity. This universal operator forms a super-set of classical Laplacian and fractional Laplacian operators and, thus, has the potential to fit a broad spectrum of data sets. We also provide theoretical convergence rates with respect to $\delta$ and verify them via numerical experiments. Second, we use nPINNs to estimate the two parameters, $\delta$ and $\alpha$, characterizing the kernel of the unified operator. The strong non-convexity of the loss function yielding multiple (good) local minima reveals the occurrence of the operator mimicking phenomenon, that is, different pairs of estimated parameters could produce multiple solutions of comparable accuracy. Third, we propose another nonlocal operator with spatially variable order $\alpha(y)$, which is more suitable for modeling turbulent Couette flow. Our results show that nPINNs can jointly infer this function as well as $\delta$. More importantly, these parameters exhibit a universal behavior with respect to the Reynolds number, a finding that contributes to our understanding of nonlocal interactions in wall-bounded turbulence.

97 MATHEMATICS AND COMPUTING↗

nPINNs: nonlocal Physics-Informed Neural Networks for a parametrized nonlocal universal Laplacian operator. Algorithms and Applications

Physics-informed neural networks (PINNs) are effective in solving inverse problems based on differential and integro-differential equations with sparse, noisy, unstructured, and multifidelity data. PINNs incorporate all available information, including governing equations (reflecting physical laws), initial-boundary conditions, and observations of quantities of interest, into a loss function to be minimized, thus recasting the original problem into an optimization problem. In this paper, we extend PINNs to parameter and function inference for integral equations such as nonlocal Poisson and nonlocal turbulence models, and we refer to them as nonlocal PINNs (nPINNs). The contribution of the paper is three-fold. First, we propose a unified nonlocal Laplace operator, which converges to the classical Laplacian as one of the operator parameters, the nonlocal interaction radius δ goes to zero, and to the fractional Laplacian as δ goes to infinity. This universal operator forms a super-set of classical Laplacian and fractional Laplacian operators and, thus, has the potential to fit a broad spectrum of data sets. We provide theoretical convergence rates with respect to δ and verify them via numerical experiments. Second, we use nPINNs to estimate the two parameters, δ and α, characterizing the kernel of the unified operator. The strong non-convexity of the loss function yielding multiple (good) local minima reveals the occurrence of the operator mimicking phenomenon, that is, different pairs of estimated parameters could produce multiple solutions of comparable accuracy. Third, we propose another nonlocal operator with spatially variable order α(γ), which is more suitable for modeling turbulent Couette flow. Our results show that nPINNs can jointly infer this function as well as δ. More importantly, these parameters exhibit a universal behavior with respect to the Reynolds number, a finding that contributes to our understanding of nonlocal interactions in wall-bounded turbulence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Nambu and Compositeness

The Nambu--Jona-Lasinio model is the simplest field theory of a composite scalar boson, consisting of a pair of chiral fermions. A bound state emerges from an assumed point-like 4-fermion interaction and is described by local effective field, $\Phi(x)$. We review this in the context of the renormalization group and show some its phenomenological successes, such as the QCD chiral dynamics and the prediction of stable heavy-light resonances that reveal the single light quark chiral dynamics. We also review some NJL--inspired composite Higgs models. We then describe a novel UV completion of the NJL model, to an extended interaction, where the solution is described by a bilocal effective, Yukawa field, $\Phi(x,y)$. I conclude with a personal recollection of Yoichiro Nambu.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Electromagnetic entrapment in gravity

Abstract We derive specific properties of electromagnetism when gravitational effects are not negligible and analyze their impact on new physics at the horizons of black holes. We show that a neutral configuration of charges in a region of high redshift, characterized by a largeg tt , produces a highly localized electromagnetic field that vanishes just beyond that region. This phenomenon implies the existence of extensive families of spacetime structures generated by electromagnetic degrees of freedom that are as compact as black holes. We construct neutral bound states of extremal black holes in four dimensions and in five dimensions, where one direction is compact. These geometries are indistinguishable from a neutral black hole, referred to asdistorted Schwarzschild, except in an infinitesimal region near its horizon where the entrapped electromagnetic structures start to manifest. The five-dimensional solutions satisfy various criteria for describing black hole microstructure: they increase in size with the Newton constant, are as compact as the Schwarzschild black hole, and have an entropy that scales likeM 2 .

Physics↗

On the accurate long-time solution of the wave equation in exterior domains: Asymptotic expansions and corrected boundary conditions

We consider the solution of scattering problems for the wave equation using approximate boundary conditions at artificial boundaries. These conditions are explicitly viewed as approximations to an exact boundary condition satisfied by the solution on the unbounded domain. We study the short and long term behavior of the error. It is provided that, in two space dimensions, no local in time, constant coefficient boundary operator can lead to accurate results uniformly in time for the class of problems we consider. A variable coefficient operator is developed which attains better accuracy (uniformly in time) than is possible with constant coefficient approximations. The theory is illustrated by numerical examples. We also analyze the proposed boundary conditions using energy methods, leading to asymptotically correct error bounds.

Hagstrom, Thomas↗

Bootstrapping boundaries and branes

The study of conformal boundary conditions for two-dimensional conformal field theories (CFTs) has a long history, ranging from the description of impurities in one-dimensional quantum chains to the formulation of D-branes in string theory. Nevertheless, the landscape of conformal boundaries is largely unknown, including in rational CFTs, where the local operator data is completely determined. We initiate a systematic bootstrap study of conformal boundaries in 2d CFTs by investigating the bootstrap equation that arises from the open-closed consistency condition of the annulus partition function with identical boundaries. We find that this deceivingly simple bootstrap equation, when combined with unitarity, leads to surprisingly strong constraints on admissible boundary states. In particular, we derive universal bounds on the tension (boundary entropy) of stable boundary conditions, which provide a rigorous diagnostic for potential D-brane decays. We also find unique solutions to the bootstrap problem of stable branes in a number of rational CFTs. Along the way, we observe a curious connection between the annulus bootstrap and the sphere packing problem, which is a natural extension of previous work on the modular bootstrap. We also derive bounds on the boundary entropy at large central charge. These potentially have implications for end-of-the-world branes in pure gravity on AdS 3 .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Motion Planning Algorithms for Safety and Quantum Computing Efficiency

Motion planning remains a fundamental problem in robotics. Sampling-based algorithms use randomization to allow efficient solutions to this complex problem. As mobile robots and autonomous vehicles become more prevalent in everyday life, motion planning must be applied to increasingly challenging scenarios. Safety has become a paramount concern in motion planning for ensuring robotic applications enrich human lives. To date, many motion planning techniques to increase safety in the face of uncertain and dynamic environments have been developed. This dissertation first addresses distributional safety of Rapidly-Exploring Random Trees (RRT) through our algorithm W-Safe RRT. To acknowledge distributional uncertainty and poor modeling, W-Safe RRT uses the Wasserstein metric to provide a probabilistic bound on the distributional distance between a robot and obstacles. Human-interpretable environmental agent classification allows online safety margin adaptation. We propose and analyze an integrating region method for online classification that increases actor labeling accuracy based on behavioral feature values when compared to state of the art methods. The method performs class assignments based on local maximum likelihood in a created behavioral feature-space, allowing a notion of classification uncertainty. Model-based methods with safety guarantees can quickly become computationally in tractable, especially with multiple agents, higher dimensions, and plentiful unknowns. Sampling based algorithms have been parallelized for computation with multi-core computers and GPUs. We consider the use of quantum algorithms and computers for sampling-based motion planning for the first time. Quantum computing performs operations on superpositions of states and can solve certain problems much more efficiently than classical computers, but introduces previously unseen challenges. With Quantum-RRT, we recast the motion planning problem into a database-search structure and use Quantum Amplitude Amplification to find reachable states in the database with a quadratic performance increase over classical methods. We address two error sources with this method: quantum measurement and quantum oracle errors. We then extend this method to Parallel Quantum-RRT, which uses a manager-worker architecture with multiple parallel quantum workers to increase database search efficiency. We compare algorithm architectures and characterize probabilities of multiple workers finding solutions. Lastly, we test in simulation the quantum algorithms against classical versions in a wide variety of scenarios, concluding that a similar parallelization improvement is to be found in the quantum case as was found in the parallelization of classical RRT.

97 MATHEMATICS AND COMPUTING↗

A method for bounding high-order finite element functions: Applications to mesh validity and bounds-preserving limiters

We introduce a novel method for bounding high-order multi-dimensional polynomials in finite element approximations. The method involves precomputing optimal piecewise-linear bounding boxes for polynomial basis functions, which can then be used to locally bound any combination of these basis functions. This approach can be applied to any element/basis type at any approximation order, can provide local (i.e., subcell) extremum bounds to a desired level of accuracy, and can be evaluated efficiently on-the-fly in simulations. Furthermore, we show that this approach generally yields more accurate bounds in comparison to traditional methods based on convex hull properties (e.g., Bernstein polynomials). Furthermore, the efficacy of this technique is shown in applications such as mesh validity checks and optimization for high-order curved meshes, where positivity of the element Jacobian determinant can be ensured throughout the entire element, and continuously bounds-preserving limiters for hyperbolic systems, which can enforce maximum principle bounds across the entire solution polynomial.

Bounding box↗

General field evaluation in high-order meshes on GPUs

Robust and scalable function evaluation at any arbitrary point in the finite/spectral element mesh is required for querying the partial differential equation solution at points of interest, comparison of solution between different meshes, and Lagrangian particle tracking. This is a challenging problem, particularly for high-order unstructured meshes partitioned in parallel with MPI, as it requires identifying the element that overlaps a given point and computing the corresponding reference space coordinates. Here, we present a robust and efficient technique for general field evaluation in large-scale high-order meshes with quadrilaterals and hexahedra. In the proposed method, a combination of globally partitioned and processor-local maps are used to first determine a list of candidate MPI ranks, and then locally candidate elements that could contain a given point. Next, element-wise bounding boxes further reduce the list of candidate elements. Finally, Newton’s method with trust region is used to determine the overlapping element and corresponding reference space coordinates. Since GPU-based architectures have become popular for accelerating computational analyses using meshes with tensor-product elements, specialized kernels have been developed to utilize the proposed methodology on GPUs. The method is also extended to enable general field evaluation on surface meshes. The paper concludes by demonstrating the use of the proposed method in various applications ranging from mesh-to-mesh transfer during r-adaptivity to Lagrangian particle tracking.

97 MATHEMATICS AND COMPUTING↗

Water Network Shape-Dependence of Local Interactions with the Microhydrated –NO 2 – and –CO 2 – Anionic Head Groups by Cold Ion Vibrational Spectroscopy

We report the structural evolutions of water networks and solvatochromic response of the CH 3 NO 2 – radical anion in the OH and CH stretching regions by analysis of the vibrational spectra displayed by cryogenically cooled CH 3 NO 2 – ·(H 2 O) n=1–6 clusters. The OH stretching bands evolve with a surprisingly large discontinuity at n = 6, which features the emergence of an intense, strongly red-shifted band along with a weaker feature that appears in the region assigned to a free OH fundamental. Very similar behavior is displayed by the perdeuterated carboxylate clusters, RCO 2 – ·(H 2 O) n=5–7 (R = CD 3 CD 2 ), indicating that this behavior is a general feature in the microhydration of the triatomic anionic domain and not associated with CH oscillators. Electronic structure calculations trace this behavior to the formation of a “book” isomer of the water hexamer that adopts a configuration in which one of the water molecules resides in an acceptor–acceptor–donor (AAD) (A = acceptor, D = donor) H-bonding site. Excitation of the bound OH in the AAD site explores the local network topology best suited to stabilize an incipient –XO 2 H–OH–(H 2 O) 2 intracluster proton-transfer reaction. These systems thus provide particularly clear examples where the network shape controls the potential energy landscape that governs water network-mediated, intracluster proton transfer. The CH stretching bands of the CH 3 NO 2 – ·(H 2 O) n=1–6 clusters also exhibit strong solvatochromic shifts, but in this case, they smoothly blue-shift with increasing hydration with no discontinuity at n = 6. Furthermore, this behavior is analyzed in the context of the solute-ion polarizability response and partial charge transfer to the water networks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Corrected calculation for the non-local solution to the g − 2 anomaly and novel results in non-local QED

We provide the corrected calculation of the (g − 2) μ in non-local QED previously done in the literature. In specific, we show the proper technique for calculating loops in non-local QED and use it to find the form factors F 1 (q 2 ) and F 2 (q 2 ) in non-local QED. We also utilize this technique to calculate some novel results in non-local QED, including calculating the correction to the photon self-energy, the modification to the classical Coulomb potential, the modification to the energy levels of the hydrogen atom, and the contribution to the Lamb shift. We also discuss charge dequantization through non-locality, and show that the experimental bounds on the electric charge on Dirac neutrinos, translate into strong flavor-dependent bounds on the scale on non-locality that range between 10 5 −10 10 TeV. We also discuss the inconsistencies of unrenormalized non-local Quantum Field Theories (QFTs) and the need for renormalizing them, even when they are free from UV divergences.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Long-Range Allosteric Communication Modulated by Active Site Mn(II) Coordination Drives Catalysis in Xanthobacter autotrophicus Acetone Carboxylase

Acetone carboxylase (AC) from Xanthobacter autotrophicus is a 360 KDa α2β2γ2 heterohexamer that catalyzes the ATP-dependent formation of phosphorylated acetone and bicarbonate intermediates that react at Mn(II) metal active sites to form acetoacetate. Structural models of X. autotrophicus AC (XaAC) with and without nucleotides reveal that the binding and phosphorylation of the two substrates occurs ~40 Å from the Mn(II) active sites where acetoacetate is formed. Based on the crystal structures, a significant conformational change was proposed to open and close a tunnel that facilitates the passage of reaction intermediates between the sites for nucleotide binding and phosphorylation of substrates and Mn(II) sites of acetoacetate formation. We have employed electron paramagnetic resonance (EPR), kinetic assays, and hydrogen/deuterium exchange mass spectrometry (HDX-MS) of poised ligand-bound states and site-specific amino acid variants to complete an in-depth analysis of Mn(II) coordination and allosteric communication throughout the catalytic cycle. In contrast with the established paradigms for carboxylation, our analyses of XaAC suggested a carboxylate shift that couples both local and long-range structural transitions. Shifts in the coordination mode of a single carboxylic acid residue (αE89) mediate both catalysis proximal to a Mn(II) center and communication with an ATP active site in a separate subunit of a 180 kDa α2β2γ2 complex at a distance of 40 Å. This work demonstrates the power of combining structural models from X-ray crystallography with solution-phase spectroscopy and biophysical techniques to elucidate functional aspects of a multi-subunit enzyme.

Biochemistry & Molecular Biology↗

NICER Discovers the Ultracompact Orbit of the Accreting Millisecond Pulsar IGR J17062-6143

We present results of recent Neutron Star Interior Composition Explorer (NICER) observations of the accreting millisecond X-ray pulsar (AMXP) IGR J17062-6143 that show that it resides in a circular, ultracompact binary with a 38-minute orbital period. NICER observed the source for 26 kiloseconds over a 5.3-day span in 2017 August, and again for 14 and 11 kiloseconds in 2017 October and November, respectively. A power spectral analysis of the August exposure confirms the previous detection of pulsations at 163.656 Hertz in Rossi X-ray Timing Explorer (RXTE) data, and reveals phase modulation due to orbital motion of the neutron star. A coherent search for the orbital solution using the Z squared method finds a best-fitting circular orbit with a period of 2278.21 seconds (37.97 minutes), a projected semimajor axis of 0.00390 lt-s (Localization Test Statistic), and a barycentric pulsar frequency of 163.6561105 Hertz. This is currently the shortest known orbital period for an AMXP. The mass function is 9.12 times 10 (sup minus 8) solar mass, presently the smallest known for a stellar binary. The minimum donor mass ranges from approximately 0.005 to 0.007 times the solar mass for a neutron star mass from 1.2 to 2 times the solar mass. Assuming mass transfer is driven by gravitational radiation, we find donor mass and binary inclination bounds of 0.0175-0.0155 times the solar mass and 19 degrees less than i less than 27.5 degrees, where the lower and upper bounds correspond to 1.4 and 2 times the solar mass neutron stars, respectively. Folding the data accounting for the orbital modulation reveals a sinusoidal profile with fractional amplitude 2.04 plus or minus 0.11 percent (0.3-3.2 kiloelectronvolts).

Strohmayer, T. E.↗

Thread Batching for High-performance Energy-efficient GPU Memory Design

Massive multi-threading in GPU imposes tremendous pressure on memory subsystems. Due to rapid growth in thread-level parallelism of GPU and slowly improved peak memory bandwidth, memory becomes a bottleneck of GPU’s performance and energy efficiency. In this article, we propose an integrated architectural scheme to optimize the memory accesses and therefore boost the performance and energy efficiency of GPU. First, we propose a thread batch enabled memory partitioning (TEMP) to improve GPU memory access parallelism. In particular, TEMP groups multiple thread blocks that share the same set of pages into a thread batch and applies a page coloring mechanism to bound each stream multiprocessor (SM) to the dedicated memory banks. After that, TEMP dispatches the thread batch to an SM to ensure high-parallel memory-access streaming from the different thread blocks. Second, a thread batch-aware scheduling (TBAS) scheme is introduced to improve the GPU memory access locality and to reduce the contention on memory controllers and interconnection networks. Experimental results show that the integration of TEMP and TBAS can achieve up to 10.3% performance improvement and 11.3% DRAM energy reduction across diverse GPU applications. We also evaluate the performance interference of the mixed CPU+GPU workloads when they are run on a heterogeneous system that employs our proposed schemes. Our results show that a simple solution can effectively ensure the efficient execution of both GPU and CPU applications.

Computer Science↗

Americium Oxalate: An Experimental and Computational Investigation of Metal–Ligand Bonding

A novel actinide-containing coordination polymer, [Am(C 2 O 4 )(H 2 O) 3 Cl] (Am-1), has been synthesized and structurally characterized. The crystallographic analysis reveals that the structure is two-dimensional and comprised of pseudo-dimeric Am 3+ nodes that are bridged by oxalate ligands to form sheets. Each metal center is nine-coordinate, forming a distorted capped square antiprism geometry with a C 1 symmetry, and features bound oxalate, aqua, and chloro ligands. The Am 3+ –ligand bonds were probed computationally using the quantum theory of atoms in molecules nd natural localized molecular orbital approaches to investigate the underlying mechanisms and hybrid atomic orbital contributions therein. The analyses indicate that the bonds within Am-1 are predominantly ionic and the 5f shell of the Am 3+ metal centers does not add a significant covalent contribution to the bonds. Further, our bonding assessment is supported by measurements on the optical properties of Am-1 using diffuse reflectance and photoluminescence spectroscopies. The position of the principal absorption band at 507 nm ( 5 L 6' ← 7 F 0' ) is notable because it is consistent with previously reported americium oxalate complexes in solution, indicating similarities in the electronic structure and ionic bonding. Compound Am-1 is an active phosphor, featuring strong bright-blue oxalate-based luminescence with no evidence of metal-centered emission.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Overlapping Schwarz Methods Are Not Anisotropy‐Robust Multigrid Smoothers

We analyze overlapping multiplicative Schwarz methods as smoothers in the geometric multigrid solution of two-dimensional anisotropic diffusion problems. For diffusion equations, it is well known that the smoothing properties of point-wise smoothers, such as Gauss Seidel, rapidly deteriorate as the strength of anisotropy increases. On the other hand, global smoothers based on line smoothing are known to generally provide good smoothing for diffusion problems, independent of the anisotropy strength. Here, a natural question is whether global methods are really necessary to achieve good smoothing in such problems, or whether it can be obtained with locally overlapping block smoothers using sufficiently large blocks and overlap. Through local Fourier analysis and careful numerical experimentation, we show that global methods are indeed necessary to achieve anisotropy-robust smoothing. Specifically, for any fixed block size bounded sufficiently far away from the global domain size, we find that the smoothing properties of overlapping multiplicative Schwarz rapidly deteriorate with increasing anisotropy, irrespective of the amount of overlap between blocks. Moreover, our results indicate that anisotropy-robust smoothing requires blocks of diameter 𝒪⁡(𝜖 −1/2 ) for anisotropy ratio 𝜖 ∈(0,1] .

97 MATHEMATICS AND COMPUTING↗

The Principles of Turbulent Heat Transfer

The literature on turbulent heat transfer has in the course of years attained a considerable volume. Since this very complicated problem has not as yet found a complete solution, further studies in this field may be expected. The heat engineer must therefore accomodate himself to a constantly increasing number of theories and formulas. Since the theories generally start from hypothetical assumptions, and since they contain true and false assertions, verified knowledge and pure suppositions often being intermingled in a manner difficult to tell them apart, the specialist had difficulty in forming a correct evaluation of the individual studies. The need therefore arises for a presentation of the problem of turbulent heat transfer which is not initially bound by hypothetical assumptions and in which uninvestigated can be clearly distinguished form each other. Such a presentation will be given in the present treatment. Brief remarks with regard to the development of the theory of local heat transfer are included.

Reichardt, H.↗