Search NASA⌕ Search

SEARCH · Search NASA

Results for “Normal forms”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS2007

Broadband Phase Retrieval for Image-Based Wavefront Sensing

A focus-diverse phase-retrieval algorithm has been shown to perform adequately for the purpose of image-based wavefront sensing when (1) broadband light (typically spanning the visible spectrum) is used in forming the images by use of an optical system under test and (2) the assumption of monochromaticity is applied to the broadband image data. Heretofore, it had been assumed that in order to obtain adequate performance, it is necessary to use narrowband or monochromatic light. Some background information, including definitions of terms and a brief description of pertinent aspects of image-based phase retrieval, is prerequisite to a meaningful summary of the present development. Phase retrieval is a general term used in optics to denote estimation of optical imperfections or aberrations of an optical system under test. The term image-based wavefront sensing refers to a general class of algorithms that recover optical phase information, and phase-retrieval algorithms constitute a subset of this class. In phase retrieval, one utilizes the measured response of the optical system under test to produce a phase estimate. The optical response of the system is defined as the image of a point-source object, which could be a star or a laboratory point source. The phase-retrieval problem is characterized as image-based in the sense that a charge-coupled-device camera, preferably of scientific imaging quality, is used to collect image data where the optical system would normally form an image. In a variant of phase retrieval, denoted phase-diverse phase retrieval [which can include focus-diverse phase retrieval (in which various defocus planes are used)], an additional known aberration (or an equivalent diversity function) is superimposed as an aid in estimating unknown aberrations by use of an image-based wavefront-sensing algorithm. Image-based phase-retrieval differs from such other wavefront-sensing methods, such as interferometry, shearing interferometry, curvature wavefront sensing, and Shack-Hartmann sensing, all of which entail disadvantages in comparison with image-based methods. The main disadvantages of these non-image based methods are complexity of test equipment and the need for a wavefront reference.

Dean, Bruce H.↗
DOE OSTI2023

Incremental Interval Assignment by Integer Linear Algebra with Improvements

Interval Assignment (IA) is the problem of selecting the number of mesh edges (intervals) for each curve for conforming quad and hex meshing. The intervals x is fundamentally integer-valued. Many other approaches perform numerical optimization then convert a floating-point solution into an integer solution, which is slow and error prone. We avoid such steps: we start integer, and stay integer. Incremental Interval Assignment (IIA) uses integer linear algebra (Hermite normal form) to find an initial solution to the meshing constraints, satisfying the integer matrix equation Solving for reduced row echelon form provides integer vectors spanning the nullspace of A. Here we add vectors from the nullspace to improve the initial solution, maintaining Ax = b Heuristics find good integer linear combinations of nullspace vectors that provide strict improvement towards variable bounds or goals. IIA always produces an integer solution if one exists. In practice we usually achieve solutions close to the user goals, but there is no guarantee that the solution is optimal, nor even satisfies variable bounds, e.g. has positive intervals. We describe several algorithmic changes since first publication that tend to improve the final solution. The software is freely available.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Mushroom body input connections form independently of sensory activity in Drosophila melanogaster

Associative brain centers, such as the insect mushroom body, need to represent sensory information in an efficient manner. In Drosophila melanogaster, the Kenyon cells of the mushroom body integrate inputs from a random set of olfactory projection neurons, but some projection neurons—namely those activated by a few ethologically meaningful odors—connect to Kenyon cells more frequently than others. This biased and random connectivity pattern is conceivably advantageous, as it enables the mushroom body to represent a large number of odors as unique activity patterns while prioritizing the representation of a few specific odors. How this connectivity pattern is established remains largely unknown. Here, we test whether the mechanisms patterning the connections between Kenyon cells and projection neurons depend on sensory activity or whether they are hardwired. We mapped a large number of mushroom body input connections in partially anosmic flies—flies lacking the obligate odorant co-receptor Orco—and in wild-type flies. Statistical analyses of these datasets reveal that the random and biased connectivity pattern observed between Kenyon cells and projection neurons forms normally in the absence of most olfactory sensory activity. Finally, this finding supports the idea that even comparatively subtle, population-level patterns of neuronal connectivity can be encoded by fixed genetic programs and are likely to be the result of evolved prioritization of ecologically and ethologically salient stimuli.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Bicrystallography-informed Frenkel–Kontorova model for interlayer dislocations in strained 2D heterostructures

In recent years, van der Waals (vdW) heterostructures and homostructures, which consist of stacks of two-dimensional (2D) materials, have risen to prominence due to their association with exotic quantum phenomena originating from correlated electronic states harbored by them. Atomistic scale relaxation effects play an extremely important role in the electronic scale quantum physics of these systems, providing means of manipulation of these materials and allowing them to be tailored for emergent technologies. We investigate such structural relaxation effects in this work using atomistic and mesoscale models, within the context of twisted bilayer graphene — a well-known heterostructure system that features moiré patterns arising from the lattices of the two graphene layers. For small twist angles, atomic relaxation effects in this system are associated with the natural emergence of interface dislocations or strain solitons, which result from the cyclic nature of the generalized stacking fault energy (GSFE), that measures the interface energy based on the relative movement of the two layers. Here, in this work, we first demonstrate using atomistic simulations that atomic reconstruction in bilayer graphene under a large twist also results from interface dislocations, although the Burgers vectors of such dislocations are considerably smaller than those observed in small-twist systems. To reveal the translational invariance of the heterointerface responsible for the formation of such dislocations, we derive the translational symmetry of the GSFE of a 2D heterostructure using the notions of coincident site lattices (CSLs) and displacement shift complete lattices (DSCLs). The workhorse for this exercise is a recently developed Smith normal form bicrystallography framework. Next, we construct a bicrystallography-informed and frame-invariant Frenkel–Kontorova model, which can predict the formation of strain solitons in arbitrary 2D heterostructures, and apply it to study a heterostrained, large-twist bilayer graphene system. Our mesoscale model is found to produce results consistent with atomistic simulations. We anticipate that the model will be invaluable in predicting structural relaxation and for providing insights into various heterostructure systems, especially in cases where the fundamental unit cell is large and therefore, atomistic simulations are computationally expensive.

2D heterostructures↗
DOE OSTI2024

New Calabi–Yau manifolds from genetic algorithms

Calabi–Yau manifolds can be obtained as hypersurfaces in toric varieties built from reflexive polytopes. We generate reflexive polytopes in various dimensions using a genetic algorithm. As a proof of principle, we demonstrate that our algorithm reproduces the full set of reflexive polytopes in two and three dimensions, and in four dimensions with a small number of vertices and points. Motivated by this result, we construct five-dimensional reflexive polytopes with the lowest number of vertices and points. By calculating the normal form of the polytopes, we establish that many of these are not in existing datasets and therefore give rise to new Calabi–Yau four-folds. In some instances, the Hodge numbers we compute are new as well.

Astronomy & Astrophysics↗
DOE OSTI2024

Stellarator equilibrium axis-expansion to all orders in distance from the axis for arbitrary plasma beta

A systematic theory of the asymptotic expansion of the magnetohydrostatics (MHS) equilibrium in the distance from the magnetic axis is developed to include arbitrary smooth currents near the magnetic axis. Compared with the vacuum and the force-free system, an additional magnetic differential equation must be solved to obtain the pressure-driven currents. It is shown that there exist variables in which the rest of the MHS system closely mimics the vacuum system. Thus, a unified treatment of MHS fields is possible. The mathematical structure of the near-axis expansions to arbitrary order is examined carefully to show that the double-periodicity of physical quantities in a toroidal domain can be satisfied order by order. The essential role played by the leading-order Birkhoff–Gustavson normal form in solving the magnetic differential equations is highlighted. Several explicit examples of vacuum, force-free and MHS equilibrium in different geometries are presented.

Physics↗
DOE OSTI2022

Learning emergent partial differential equations in a learned emergent space

We propose an approach to learn effective evolution equations for large systems of interacting agents. This is demonstrated on two examples, a well-studied system of coupled normal form oscillators and a biologically motivated example of coupled Hodgkin-Huxley-like neurons. For such types of systems there is no obvious space coordinate in which to learn effective evolution laws in the form of partial differential equations. In our approach, we accomplish this by learning embedding coordinates from the time series data of the system using manifold learning as a first step. In these emergent coordinates, we then show how one can learn effective partial differential equations, using neural networks, that do not only reproduce the dynamics of the oscillator ensemble, but also capture the collective bifurcations when system parameters vary. The proposed approach thus integrates the automatic, data-driven extraction of emergent space coordinates parametrizing the agent dynamics, with machine-learning assisted identification of an emergent PDE description of the dynamics in this parametrization.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Task-oriented machine learning surrogates for tipping points of agent-based models

We present a machine learning framework bridging manifold learning, neural networks, Gaussian processes, and Equation-Free multiscale approach, for the construction of different types of effective reduced order models from detailed agent-based simulators and the systematic multiscale numerical analysis of their emergent dynamics. The specific tasks of interest here include the detection of tipping points, and the uncertainty quantification of rare events near them. Our illustrative examples are an event-driven, stochastic financial market model describing the mimetic behavior of traders, and a compartmental stochastic epidemic model on an Erdös-Rényi network. We contrast the pros and cons of the different types of surrogate models and the effort involved in learning them. Importantly, the proposed framework reveals that, around the tipping points, the emergent dynamics of both benchmark examples can be effectively described by a one-dimensional stochastic differential equation, thus revealing the intrinsic dimensionality of the normal form of the specific type of the tipping point. This allows a significant reduction in the computational cost of the tasks of interest.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Multiscale analysis of large twist ferroelectricity and swirling dislocations in bilayer hexagonal boron nitride

With its atomically thin structure and intrinsic ferroelectric properties, heterodeformed bilayer hexagonal boron nitride (hBN) has gained prominence in next-generation non-volatile memory applications. However, studies to date have focused almost exclusively on small-twist bilayer hBN, leaving the question of whether ferroelectricity can persist under small heterostrain and large heterodeformation entirely unexplored. In this work, we establish the crystallographic origin of ferroelectricity in bilayer hBN configurations heterodeformed relative to high-symmetry configurations such as AA-stacking and 21.786789° twisted configurations (Σ7), using Smith normal form bicrystallography. We then demonstrate out-of-plane ferroelectricity in bilayer hBN across configurations vicinal to both the AA and Σ7 stackings. Atomistic simulations reveal that AA-vicinal systems support ferroelectricity under both small twist and small strain, with polarization switching in the latter governed by the deformation of swirling dislocations rather than the straight interface dislocations seen in the former. For Σ7-vicinal systems, where existing interatomic potentials underperform particularly under extreme out-of-plane compression, we develop a density-functional-theory-informed continuum framework—the bicrystallography-informed frame-invariant multiscale (BFIM) model, which captures out-of-plane ferroelectricity in heterodeformed configurations vicinal to Σ7 stacking. Interface dislocations in these large heterodeformed bilayer configurations exhibit markedly smaller Burgers vectors compared to interface dislocations in small-twist and small-strain bilayer hBN. The BFIM model reproduces experimental results and provides a powerful, computationally efficient framework for predicting ferroelectricity in large-unit-cell heterostructures where atomistic simulations are prohibitively expensive.

Ahmed, Md Tusher [Univ. of Illinois at Urbana-Cham↗
DOE OSTI2022

Solving quasi-free and quadratic Lindblad master equations for open fermionic and bosonic systems

The dynamics of Markovian open quantum systems are described by Lindblad master equations. For fermionic and bosonic systems that are quasi-free, i.e. with Hamiltonians that are quadratic in the ladder operators and Lindblad operators that are linear in the ladder operators, we derive the equation of motion for the covariance matrix. This determines the evolution of Gaussian initial states and the steady states, which are also Gaussian. Using ladder super-operators (a.k.a. third quantization), here we show how the Liouvillian can be transformed to a many-body Jordan normal form which also reveals the full many-body spectrum. Extending previous work by Prosen and Seligman, we treat fermionic and bosonic systems on equal footing with Majorana operators, shorten and complete some derivations, also address the odd-parity sector for fermions, give a criterion for the existence of bosonic steady states, cover non-diagonalizable Liouvillians also for bosons, and include quadratic systems. In extension of the quasi-free open systems, quadratic open systems comprise additional Hermitian Lindblad operators that are quadratic in the ladder operators. While Gaussian states may then evolve into non-Gaussian states, the Liouvillian can still be transformed to a useful block-triangular form, and the equations of motion for k-point Green’s functions form a closed hierarchy. Based on this formalism, results on criticality and dissipative phase transitions in such models are discussed in a companion paper.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Dark and bright modes, and their coherent control in dipolar metasurface bilayers

Several plasmonic nanoparticles supporting dipolar resonances can couple to form normal modes. Here, we develop an analytical model to explain the formation of nonradiative “dark” and radiative “bright” modes through radiative coupling in bilayers consisting of dipolar nanoantenna arrays that are separated by a subwavelength distance. We also include near-field contributions in our model and show that the absorption and reflectance spectra obtained from our model agree reasonably well with the respective finite-difference time-domain simulation results for both perfectly aligned and misaligned bilayers. Here, the ability to vary the reflection and absorption spectra of these bilayers by changing the material and geometrical parameters has potential applications in the design of efficient spectral filters. We also show that we can selectively excite these modes by adjusting the phase between two counterpropagating normally incident fields, which has applications in all-optical modulators and switches based on purely linear interferometric effects.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

An Analytical Solution for the Initiation and Early Progression of Fretting Wear in Spherical Contacts

Abstract This article derives analytical solutions to calculate the wear volume at the initiation of fretting motion and its early progression over the first few oscillation cycles. The Archard-based model considers a deformable hemisphere that is contact with a deformable flat block. The material pairs investigated are special alloys, the Inconel 617/Incoloy 800H, and Inconel 617/Inconel 617. The analytical study begins with a unidirectional frictional sliding contact, where the local interfacial sliding distance and the nominal sliding distance at the initiation of gross slip are derived. The obtained analytical expressions for unidirectional sliding are then used to derive the corresponding wear volume for the initiation and early progression of gross slip and the wear volume for a general fretting cycle under elastic conditions. These analytical derivations are all verified by the finite element analysis (FEA). The FEA method and the analytical solutions render virtually identical results for both similar and dissimilar material pairs. The effects of plasticity on the wear volume under elastic–plastic conditions are also investigated. It is found that the fretting wear volumes obtained from the FEA simulations, which include plasticity, are close to those obtained from the analytical expressions for purely elastic regimes. All the results are presented in normalized forms, which can be easily generalized and applied to three-dimensional fretting wear of other material pairs.

Engineering↗
DOE OSTI2024

QED: A Powerful Query Equivalence Decider for SQL

Checking query equivalence is of great significance in database systems. Prior work in automated query equivalence checking sets the first steps in formally modeling and reasoning about query optimization rules, but only supports a limited number of query features. In this paper, we present Qed, a new framework for query equivalence checking based on bag semantics. Qed uses a new formalism called Q-expressions that models queries using different normal forms for efficient equivalence checking, and models features such as integrity constraints and NULLs in a principled way unlike prior work. Our formalism also allows us to define a new query fragment that encompasses many real-world queries with a complete equivalence checking algorithm, assuming a complete first-order theory solver. Empirically, Qed can verify 299 out of 444 query pairs extracted from the Calcite framework and 979 out of 1287 query pairs extracted from CockroachDB, which is more than 2× the number of cases proven by prior state-of-the-art solver.

Computer Science↗
DOE OSTI2025

SDA: a symbolic differential algebra package in C++

Truncated Power Series Algebra (TPSA), or Differential Algebra (DA), is a well-established tool in accelerator physics, commonly used for generating high-order maps of dynamic systems, as well as in symplectic tracking, normal form analysis, verified integration, optimization, and fast multipole methods. This package is the first to perform symbolic DA computations, enabling traceability of initial condition contributions and runtime reduction for repeated DA calculations, potentially expanding DA’s applications.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1970

On a theorem of K T Chen

Proving normal form of mappings of real line into itself by contracting mapping principle

Braun, M.↗
NASA NTRS1971

Radiofrequency recombination lines from the interstellar medium

Observations of recombination lines form normal H II regions, extended H II regions, nonthermal sources, and the H I medium are discussed. Detection of recombination lines from elements other than hydrogen may provide a means of identifying fossil Stromgren spheres at high temperature.

Dupree, A. K.↗