Search NASA⌕ Search

SEARCH · Search NASA

Results for “Approximate inverses”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Effective Model for Fractional Topological Corner Modes in Quasicrystals

High-order topological insulators (HOTIs), as generalized from topological crystalline insulators, are characterized with lower-dimensional metallic boundary states protected by spatial symmetries of a crystal, whose theoretical framework based on band inversion at special k points cannot be readily extended to quasicrystals because quasicrystals contain rotational symmetries that are not compatible with crystals, and momentum is no longer a good quantum number. In this work, we develop a low-energy effective model underlying HOTI states in 2D quasicrystals for all possible rotational symmetries. By implementing a novel Fourier transform developed recently for quasicrystals and approximating the long-wavelength behavior by their large-scale average, we construct an effective k · p Hamiltonian to capture the band inversion at the center of a pseudo-Brillouin zone. We show that an in-plane Zeeman field can induce mass kinks at the intersection of adjacent edges of a 2D quasicrystal topological insulators and generate corner modes (CMs) with fractional charge, protected by rotational symmetries. Our model predictions are confirmed by numerical tight-binding calculations. Furthermore, when the quasicrystal is proximitized by an s-wave superconductor, Majorana CMs can also be created by tuning the field strength and chemical potential. Our work affords a generic approach to studying the low-energy physics of quasicrystals, in association with topological excitations and fractional statistics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Reconnection from a turbulence perspective

The spectral properties associated with laminar, anti-parallel reconnection are examined using a 2.5D kinetic particle in cell simulation. Both the reconnection rate and the energy spectrum exhibit three distinct phases: an initiation phase where the reconnection rate grows, a quasi-steady phase, and a declining phase where both the reconnection rate and the energy spectrum decrease. During the steady phase, the energy spectrum exhibits approximately a double power law behavior, with a slope near −5/3 at wave numbers smaller than the inverse ion inertial length and a slope steeper than −8/3 for larger wave numbers up to the inverse electron inertial length. This behavior is consistent with a Kolmogorov energy cascade and implies that laminar reconnection may fundamentally be an energy cascade process. Consistent with this idea is the fact that the reconnection rate exhibits a rough correlation with the energy spectrum at wave numbers near the inverse ion inertial length. The 2D spectrum is strongly anisotropic with most energy associated with the wave vector direction normal to the current sheet. Reconnection acts to isotropize the energy spectrum, reducing the Shebalin angle from an initial value of 70° to about 48° (nearly isotropic) by the end of the simulation. The distribution of energy over length scales is further analyzed by dividing the domain into spatial subregions and employing structure functions.

Adhikari, S. (ORCID:0000000321607066)↗
DOE OSTI2025

Learning local and semi-local density functionals from exact exchange-correlation potentials and energies

Finding accurate exchange-correlation (XC) functionals remains the defining challenge in density functional theory (DFT). Despite 40 years of active development, attaining general purpose chemical accuracy is still elusive with existing functionals. We present a data-driven pathway to learn the XC functional by using the exact density, XC energy, and XC potential. While the exact densities are obtained from accurate configuration interaction (CI), the exact XC energies and XC potentials are obtained via inverse DFT calculations on the CI densities. We demonstrate how simple neural network (NN)–based local density approximation (LDA) and generalized gradient approximation (GGA), trained on just five atoms and two molecules, provide remarkable improvement in total energies and densities. Particularly, the NN-based GGA functional attains similar accuracy as the higher rung SCAN meta-GGA on various thermochemistry datasets. These results underscore the promise of using the XC potential in modeling XC functionals and can pave the way for systematic learning of increasingly accurate XC functionals.

Science & Technology - Other Topics↗
NASA NTRS1988

Active galactic nuclei. II - The acceleration of relativistic particles in a cluster of accreting black holes

An accreting cluster of black holes in an active galactic nucleus is a natural site for a system of shock structures with a hierarchy of sizes, corresponding to the distribution of masses in the cluster. Accreted gas containing some magnetic fields and supersonically falling onto the core forms shocks on the outside of each hole and these shocks are capable of accelerating relativistic particles. The energies reached in a single shock are size rather than acceleration time limited and are proportional to the mass of the hole with a proportionality constant being a function of the position of the hole within a cluster and the model of the cluster and the shock formation. These energies are adequate to explain the observed properties of synchrotron and inverse-Compton radiation from these objects. The resulting energy spectrum of particles in the cluster in 'zeroth' approximation has the form of a doubly broken power law with indices of two and three on both extremes of the energy domain respectively, bridged by an index of about 2.5.

Pacholczyk, A. G.↗
NASA NTRS1993

The dependence of the dissipation on rotation

A form of the evolution equation for the dissipation showing the effects of rotation is derived directly from the Navier-Stokes equations. Rotation affects the dissipation if the flow is inhomogeneous or compressible. This is due to the fact that pressure scales with inverse Rossby number and certain pressure correlations usually neglected have to be retained when Re Ro is approximately equal to 1. The term reflecting the effects of rotation on the dissipation is a volume integral of the two-point velocity correlation similar in structure to the 'rapid-pressure' correlation appearing in the second-moment equations. Some of its properties are indicated.

Ristorcelli, J. R.↗
NASA NTRS1999

Efficient Implementations of the Quadrature-Free Discontinuous Galerkin Method

The efficiency of the quadrature-free form of the dis- continuous Galerkin method in two dimensions, and briefly in three dimensions, is examined. Most of the work for constant-coefficient, linear problems involves the volume and edge integrations, and the transformation of information from the volume to the edges. These operations can be viewed as matrix-vector multiplications. Many of the matrices are sparse as a result of symmetry, and blocking and specialized multiplication routines are used to account for the sparsity. By optimizing these operations, a 35% reduction in total CPU time is achieved. For nonlinear problems, the calculation of the flux becomes dominant because of the cost associated with polynomial products and inversion. This component of the work can be reduced by up to 75% when the products are approximated by truncating terms. Because the cost is high for nonlinear problems on general elements, it is suggested that simplified physics and the most efficient element types be used over most of the domain.

Lockard, David P.↗
NASA NTRS1986

Do large rate coefficients for ion-polar neutral reactions have a serious effect on chemical models of dense clouds?

In order to incorporate large ion-polar neutral rate coefficients into existing gas phase reaction networks, it is necessary to utilize simplified theoretical treatments because of the significant number of rate coefficients needed. The authors have used two simple theoretical treatments: the locked dipole approach of Moran and Hamill for linear polar neutrals and the trajectory scaling approach of Su and Chesnavich for nonlinear polar neutrals. The former approach is suitable for linear species because in the interstellar medium these are rotationally relaxed to a large extent and the incoming charged reactants can lock their dipoles into the lowest energy configuration. The latter approach is a better approximation for nonlinear neutral species, in which rotational relaxation is normally less severe and the incoming charged reactants are not as effective at locking the dipoles. The treatments are in reasonable agreement with more detailed long range theories and predict an inverse square root dependence on kinetic temperature for the rate coefficient. Compared with the locked dipole method, the trajectory scaling approach results in rate coefficients smaller by a factor of approximately 2.5.

Herbst, E.↗
NASA NTRS1983

A time-dependent approach for calculating steady inverse boundary-layer flows with separation

An unsteady inverse boundary-layer method is developed which can be used to calculate steady flows with separation. Two versions of Keller's box method with the Mechul function formulation developed by Cebeci (1976) are employed, depending on the complexity of the flow. The regular box is employed in regions of positive streamwise velocity component u, whereas the zigzag box is employed in regions where u becomes negative (t greater than 0). The regular box with the FLARE approximation is employed when t = 0 and u becomes negative in some region across the layer. Results of calculations show that the use of a time-dependent inverse boundary-layer method in which time is used as an iteration parameter provides a good approach in improving the accuracy of the solutions obtained from the FLARE approximation.

Cebeci, T.↗
DOE OSTI2022

Nuclear cusps and singularities in the nonadditive kinetic potential bifunctional from analytical inversion

The non-additive kinetic potential $v^{\text{NAD}}$ is a key quantity in density-functional theory (DFT) embedding methods, such as frozen density embedding theory and partition DFT. $v^{\text{NAD}}$ is a bi-functional of electron densities $\rho_{\rm B}$ and $\rho_{\rm tot} = \rho_{\rm A} + \rho_{\rm B}$. It can be evaluated using approximate kinetic-energy functionals, but accurate approximations are challenging. The behavior of $v^{\text{NAD}}$ in the vicinity of the nuclei has long been questioned, and singularities were seen in some approximate calculations. In this article, the existence of singularities in $v^{\text{NAD}}$ is analyzed analytically for various choices of $\rho_{\rm B}$ and $\rho_{\rm tot}$, using the nuclear cusp conditions for the density and Kohn-Sham potential. It is shown that no singularities arise from smoothly partitioned ground-state Kohn-Sham densities. We confirm this result by numerical calculations on diatomic test systems HeHe, HeLi$^+$, and H$_2$, using analytical inversion to obtain a numerically exact $v^{\rm NAD}$ for the local density approximation. We examine features of $v^{\rm NAD}$ which can be used for development and testing of approximations to $v^{\rm NAD}[\rho_{\rm B},\rho_{\rm tot}]$ and kinetic-energy functionals.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1973

Effects of gravity on ontogeny in animals.

Inversion or clinostat rotation of frogs' eggs demonstrated a critical period of enhanced gravity-sensitivity prior to the first cleavage. These results were corroborated by centrifugation studies which localized the period of maximal sensitivity at approximately 20 minutes post-fertilization. Eggs of various invertebrates suspended in aqueous solution proved capable of normal development following brief ultracentrifugation. Among fly larvae, grasshopper nymphs, turtles, mice, rats, and chickens, growth rates were inversely related to G-force, and maximal chronic acceleration tolerated was inversely related to body size. Body composition data demonstrated the importance of separately evaluating fat and the fat-free portion of the body in studies of the effect of acceleration on growth in homeotherms. The attempt to evaluate the effect of weightlessness on the development of frogs' eggs in Biosatellite 2 was inconclusive for technical reasons.

Pitts, G. C.↗
DOE OSTI2026

Temperature Inversions below 1 km from a V-Band Scanning Radiometer at the North Slope of Alaska

A single-channel (56.7 GHz) scanning radiometer was deployed in August 2022 at the Atmospheric Radiation Measurement (ARM) North Slope of Alaska site near Utqiaġvik. The radiometer is designed to provide temperature profiles between 0 and 1 km every 5 min. Averaging kernels show that this single-channel radiometer, taking observations at 10 discrete elevation angles, yields approximately the same information as a seven-channel V-band radiometer scanning three elevation angles. The instrument is able to reproduce the occurrence of temperature inversions between the surface and 1 km and their strength showing a correlation of 0.85, bias of −0.6 K, and slope of 1.04 with respect to radiosondes. Uncertainty in the inversion base height varies from 50 m near the surface to ∼300 m above 0.4 km when compared with radiosondes. Conversely, the inversion top height is overestimated and has higher uncertainty due to the degrading effects of the averaging kernels on the vertical resolution of the retrievals. Thanks to the high temporal resolution of the retrievals, the diurnal cycle of boundary layer temperature was evaluated showing that the radiometer can capture some aspects of the boundary layer thermal structure. The present analysis provides an overview of the capabilities of this simple observing configuration for selected applications.

Atmospheric profilers↗
DOE OSTI2023

Uncertainty quantification for neutrino opacities in core-collapse supernovae and neutron star mergers

We perform an extensive study of the correlations between the neutrino-nucleon inverse mean free paths (IMFPs) and the underlying equation of states (EoSs). Strong interaction uncertainties in the neutrino mean free path are investigated in different density regimes. The nucleon effective mass, the nucleon chemical potentials, and the residual interactions in the medium play an important role in determining neutrino-nucleon interactions in a density-dependent manner. Here we study how the above quantities are constrained by an EoS consistent with (i) nuclear mass measurements, (ii) proton-proton scattering phase shifts, and (iii) neutron star observations. We then study the uncertainties of both the charged current and the neutral current neutrino-nucleon inverse mean free paths due to the variation of these quantities, using the Hartree-Fock+random phase approximation method. Finally, we calculate the Pearson correlation coefficients between (i) the EoS-based quantities and the EoS-based quantities; (ii) the EoS-based quantities and the IMFPs; (iii) the IMFPs and the IMFPs. We find a strong impact of residual interactions on neutrino opacity in the spin and spin-isospin channels, which are not well constrained by current nuclear modelings.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Control Design for a Marine Hydrokinetic Cycloturbine Vehicle

Abstract Marine Hydrokinetic (MHK) cycloturbines generate sustainable power by exploiting tidal currents. By powering the turbines and using pitching foils for control, a vehicle comprised of MHK cycloturbines also has the ability to station keep and maneuver. The vehicle consists of four counter-rotating cycloturbines, with hydrofoils oriented perpendicular to the flow in a paddlewheel configuration. Lift and drag generated from these foils sum together to produce thrust. An experimentally tuned simulation model that solves the six-degrees-of-freedom rigid body equations of motion for the MHK vehicle subject to hydrodynamic, hydrostatic, and propulsive forces is used to aid the design of vehicle controllers. Global feedback controllers are initially designed by applying classical control methods to an approximate linear model of the system dynamics. A higher performing nonlinear controller is designed using the nonlinear dynamic inversion (NDI) method. NDI accounts for the nonlinearities of the MHK system and therefore is suitable for a wide range of operating conditions. The response of the classical and NDI controllers to speed, depth, roll, pitch, and yaw commands are evaluated and compared in simulation. The classical controller outperforms the NDI controller for small amplitude maneuvers, although the degradation with NDI is minor. However, in the nonlinear operating regime the NDI controller outperforms the classical controller and the classical controller exhibits instability.

Automation & Control Systems↗
DOE OSTI2023

An Analytic Formula for Entraining CAPE in Midlatitude Storm Environments

Abstract This article introduces an analytic formula for entraining convective available potential energy (ECAPE) with an entrainment rate that is determined directly from an environmental sounding, rather than prescribed by the formula user. Entrainment is connected to the background environment using an eddy diffusivity approximation for lateral mixing, updraft geometry assumptions, and mass continuity. These approximations result in a direct correspondence between the storm-relative flow and the updraft radius and an inverse scaling between the updraft radius squared and entrainment rate. The aforementioned concepts, combined with the assumption of adiabatic conservation of moist static energy, yield an explicit analytic equation for ECAPE that depends entirely on state variables in an atmospheric profile and a few constant parameters with values that are established in past literature. Using a simplified Bernoulli-like equation, the ECAPE formula is modified to account for updraft enhancement via kinetic energy extracted from the cloud’s background environment. CAPE and ECAPE can be viewed as predictors of the maximum vertical velocity w max in an updraft. Hence, these formulas are evaluated using w max from past numerical modeling studies. Both of the new formulas improve predictions of w max substantially over commonly used diagnostic parameters, including undiluted CAPE and ECAPE with a constant prescribed entrainment rate. The formula that incorporates environmental kinetic energy contribution to the updraft correctly predicts instances of exceedance of by w max , and provides a conceptual explanation for why such exceedance is rare among past simulations. These formulas are potentially useful in nowcasting and forecasting thunderstorms and as thunderstorm proxies in climate change studies. Significance Statement Substantial mixing occurs between the upward-moving air currents in thunderstorms (updrafts) and the surrounding comparatively dry environmental air, through a process called entrainment. Entrainment controls thunderstorm intensity via its diluting effect on the buoyancy of air within updrafts. A challenge to representing entrainment in forecasting and predictions of the intensity of updrafts in future climates is to determine how much entrainment will occur in a given thunderstorm environment without a computationally expensive high-resolution simulation. To address this gap, this article derives a new formula that computes entrainment from the properties of a single environmental profile. This formula is shown to predict updraft vertical velocity more accurately than past diagnostics, and can be used in forecasting and climate prediction to improve predictions of thunderstorm behavior and impacts.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1984

Discrete Random Media Techniques for Microwave Modeling of Vegetated Terrain

Microwave remote sensing of agricultural crops and forested regions is studied. Long term goals of the research involve modeling vegetation so that radar signatures can be used to infer the parameters which characterize the vegetation and underlying ground. Vegetation is modeled by discrete scatterers viz, leaves, stems, branches and trunks. These are replaced by glossy dielectric discs and cylinders. Rough surfaces are represented by their mean and spectral characteristics. Average scattered power is then calculated by employing discrete random media methodology such as the distorted Born approximation or transport theory. Both coherent and incoherent multiple scattering techniques are explored. Once direct methods are developed, inversion techniques can be investigated.

Lang, R. H.↗
NASA NTRS1989

VLBI using a telescope in Earth orbit. II - Brightness temperatures exceeding the inverse Compton limit

VLBI observations at 2.3 GHz were conducted using an antenna of the Tracking and Data Relay Satellite System (TDRSS), in geosynchronous orbit, and two ground stations. Twenty-three of 24 sources were detected, on baselines as long as 2.15 earth diameters. These baseline lengths gave the interferometer much better sensitivity to high brightness temperatures than any earth-based observations. Brightness temperatures of 1-4 times the 10 to the 12th K inverse Compton limit were measured for 10 sources, suggesting bulk relativistic motion in these sources. Coherence values of approximately 85 percent for integration times of 360 s were obtained.

Linfield, R. P.↗
NASA NTRS1992

Recursive Inversion Of Externally Defined Linear Systems

Technical memorandum discusses mathematical technique described in "Recursive Inversion by Finite-Impulse-Response Filters" (ARC-12247). Technique is recursive algorithm yielding finite-impulse-response approximation of unknown single-input/single-output, causal, time-invariant, linear, real system, response of which is sequence of impulses. Useful in such diverse applications as medical diagnoses, identification of military targets, geophysical exploration, and nondestructive testing.

Bach, Ralph E., Jr.↗
NASA NTRS2013

Algorithm Development for Predicting Biodiversity Based on Phytoplankton Absorption

Ocean color remote sensing has provided the scientific community with unprecedented global coverage of chlorophyll a, an indicator of phytoplankton biomass. Together, satellite-derived chlorophyll a and knowledge of Phytoplankton Functional Types (PFTs) will improve our limited understanding of marine ecosystem responses to physiochemical climate drivers involved in carbon cycle dynamics and linkages. Using cruise data from the Gulf of Maine and the Middle Atlantic Bight (N 269 pairs of HPLC and phytoplankton absorption samples), two modeling approaches were utilized to predict phytoplankton absorption and pigments. Algorithm I predicts the chlorophyll-specific absorption coefficient (a*(sub ph) (m(exp 2) mg chl a(exp - 1)) using inputs of temperature, light, and chlorophyll a. Modeled r(exp 2) values (400700 nm) ranged from 0.79 to 0.99 when compared to in situ observations with approximately 25 percent lower r(exp 2) values in the UV region. Algorithm II-a utilizes matrix inversion analysis to predict a(sub ph)(m(exp -1), 400700 nm) and r(exp 2) values ranged from 0.89 to 0.99. The prediction of phytoplankton pigments with Algorithm II-b produced r(exp 2) values that ranged from 0.40 to 0.93. When used in combination, Algorithm I, and Algorithm II-a are able to use satellite products of SST, PAR, and chlorophyll a (Algorithm I) to predict pigment concentrations and ratios to describe the phytoplankton community. The results of this study demonstrate that the spatial variation in modeled pigment ratios differ significantly from the 10-year SeaWiFS average chlorophyll a data set. Contiguous observations of chlorophyll a and phytoplankton biodiversity will elucidate ecosystem responses with unprecedented complexity.

phytoplankton↗