Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “horizon matching”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Love symmetry in higher-dimensional rotating black hole spacetimes

We develop a method for constructing a 1-parameter family of globally-defined Love symmetry generators in rotating black hole spacetimes of general dimension. The key ingredient is to focus on the vicinity of the (physical) outer horizon, matching only the radial derivative and the outer horizon pole pieces of the Klein-Gordon operator in the black hole spacetime to the SL(2, ℝ) Casimir operator. After revisiting the 4D Kerr and 5D Myers-Perry cases, the procedure is illustrated on generalized Lense-Thirring spacetimes which describe a wide variety of slowly rotating black hole metrics in any number of dimensions. Such spacetimes are known to admit an extended tower of Killing tensor and Killing vector symmetries and, as demonstrated in this paper, allow for separability of the massive scalar wave equation in Myers-Perry-like coordinates. Interestingly, separability also occurs in the horizon-penetrating Painlevé–Gullstrand coordinates associated with the freely infalling observer who registers flat space around her all the way to singularity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2016

Horizon Based Orientation Estimation for Planetary Surface Navigation

Planetary rovers navigate in extreme environments for which a Global Positioning System (GPS) is unavailable, maps are restricted to relatively low resolution provided by orbital imagery, and compass information is often lacking due to weak or not existent magnetic fields. However, an accurate rover localization is particularly important to achieve the mission success by reaching the science targets, avoiding negative obstacles visible only in orbital maps, and maintaining good communication connections with ground. This paper describes a horizon solution for precise rover orientation estimation. The detected horizon in imagery provided by the on board navigation cameras is matched with the horizon rendered over the existing terrain model. The set of rotation parameters (roll, pitch yaw) that minimize the cost function between the two horizon curves corresponds to the rover estimated pose.

localization↗
DOE OSTI2025

Revenue-Maximizing Shared Parking and Electric Vehicle Charging Management in Multi-Unit Dwellings

In urban areas, searching for parking and electric vehicle (EV) charging can result in cruising, congestion, and environmental externalities. Recognizing the business opportunity of offering private parking and charging infrastructure access within multi-unit dwellings (MUDs) during daytime, we model a shared parking and EV charging management system. We maximize the revenue of MUD charging hubs in mixed land use, catering to public demand. Our approach accounts for the objectives of the two stakeholders involved: a demand model is fitted on the choices of EV charging users, and the supply model optimizes the allocation of parking and charging requests in an MUD parking lot. A binary integer linear programming model for the allocation of parking and charging spaces with a rolling horizon is integrated with matching rules that handle both parking and charging requests. In our numerical experiments in a neighborhood of Chicago, Illinois, we estimate the performance of the MUD parking and charging system with metrics that include revenue, number of matchings, and utilization rates. At any given time, MUDs with lower prices attract more charging requests, particularly those of longer duration, resulting in higher revenue and greater charging utilization. Dynamic pricing facilitates a more equitable distribution of requests; as MUD parking lots reach capacity and their fees increase, other MUDs become more competitive, attracting additional requests. Comparing our method against first-come-first-served and optimal-solution benchmarks, we demonstrate our model’s effectiveness in dynamically managing mixed parking and charging demand in MUD charging hubs.

electric vehicle, multi-unit dwelling, charging in↗
NASA NTRS2010

What Can Geolocated Sferics Tell Us About Terrestrial Gamma-Ray Flashes?

The Fermi Gamma-ray Burst Monitor (GBM) has been detecting TGFs with increasing sensitivity over the past two years, owing to changes in flight software that have lowered its threshold for triggering and, recently, allowed a search for TGFs weaker than those which would cause an onboard trigger. Associations between TGFs detected in the first 18 months of operation and sferics detected using the World Wide Lightning Location Network (WWLLN) show that TGF peaks and lightning discharges are simultaneous to within tens of microseconds, and that GBM triggered on TGFs that occurred up to a distance of 300 km from the sub-spacecraft position. In the work presented here, we look for associations between TGFs detected by the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) and WWLLN sferics over the same 18 months, and we compare the match rate and detection horizon of the two instruments. We also look for associations between WWLLN sferics and more recent GBM TGFs, both triggered events and weaker TGFs uncovered in our untriggered search. We discuss whether in this new mode, GBM is detecting TGFs that are more distant from the sub-spacecraft point than 300 km, or whether the weaker TGFs are instead indicative of a luminosity distribution, either because the weaker ones originate deeper in the atmosphere or because they are intrinsically dimmer.

Connaughton, V.↗
NASA NTRS2021

Flight Demonstration of the Tailored Arrival Manager

A flight demonstration of arrival traffic management automation was conducted in partnership between NASA, FAA and Boeing as an element of the latter’s ecoDemonstrator2020 flight program. For the demonstration, NASA’s prototype Tailored Arrival Manager(TAM) was used to compute efficient trajectory-based solutions to traffic management problems representing those encountered during time-based metering operations today. TAM solutions involving route modifications were uplinked to a Boeing 787 airplane using Controller Pilot Data Link Communications and seamlessly loaded into the airplane’s Flight Management System (FMS). The paper describes the concept and technology behind TAM along with the data exchanges and procedures involved with the demonstration. All TAM solutions were delivered in a timely manner and successfully integrated with the airplane’s FMS, thereby demonstrating the basic feasibility of trajectory-based arrival management using currently available data communications and avionics equipage. Although TAM solutions were not executed during this initial demonstration, a limited study of trajectory prediction accuracy was possible given that the airplane flew uninterrupted, automated descents with known route and speed intent. Analysis revealed that TAM meter-fix arrival time predictions were accurate to within ± 30 seconds for time horizons of 30 minutes or less, which matched closely with FMS predictions acquired through real-time data exchange. Top-of-descent predictions over similar time horizons were found accurate to within ± 5 nautical miles.

air traffic management↗
NASA NTRSDate not supplied

Flight Demonstration of the Tailored Arrival Manager

A flight demonstration of arrival traffic management automation was conducted in partnership between NASA, FAA and Boeing as an element of the latter’s ecoDemonstrator2020 flight program. For the demonstration, NASA’s prototype Tailored Arrival Manager(TAM) was used to compute efficient trajectory-based solutions to traffic management problems representing those encountered during time-based metering operations today. TAM solutions involving route modifications were uplinked to a Boeing 787 airplane using Controller Pilot Data Link Communications and seamlessly loaded into the airplane’s Flight Management System (FMS). The paper describes the concept and technology behind TAM along with the data exchanges and procedures involved with the demonstration. All TAM solutions were delivered in a timely manner and successfully integrated with the airplane’s FMS, thereby demonstrating the basic feasibility of trajectory-based arrival management using currently available data communications and avionics equipage. Although TAM solutions were not executed during this initial demonstration, a limited study of trajectory prediction accuracy was possible given that the airplane flew uninterrupted, automated descents with known route and speed intent. Analysis revealed that TAM meter-fix arrival time predictions were accurate to within ± 30 seconds for time horizons of 30 minutes or less, which matched closely with FMS predictions acquired through real-time data exchange. Top-of-descent predictions over similar time horizons were found accurate to within ± 5 nautical miles.

Air traffic managment↗
DOE OSTI2025

Demographic Microsimulator for Integrated Urban Systems: Adapting Panel Survey of Income Dynamics to Capture the Continuum of Life

Agent-based models (ABMs) in transportation modeling simulate activity and travel decisions at the disaggregate level of households and individuals. To do this, ABMs require detailed and realistic information on agents’ socioeconomic and demographic characteristics. Various synthetic population generators have been proposed to address this need. However, most of those currently in practice are cross-sectional in nature and do not account for the dynamics within households and individuals as they progress through life events over time. This is a major shortcoming, as literature has shown that transportation decisions are affected by the transition between and co-occurrence of life cycle events. While some demographic evolution simulators have been proposed to address this issue, they are developed using cross-sectional data and capture only a small set of life cycle events and their interdependence. Addressing these drawbacks, we propose a demographic microsimulator (DEMOS) that captures the “continuum of life” by considering a range of household- and individual-level life cycle events. DEMOS is developed using the Panel Survey of Income Dynamics, one of the world’s longest-running longitudinal surveys. The DEMOS submodels consider key life cycle events that are influenced by agents’ demographic variables. DEMOS is applied to evolve the population of the San Francisco Bay Area over a 9-year horizon. Results demonstrate how DEMOS generates life trajectories and how DEMOS outputs match the observed demographic trends. DEMOS is expected to enable longitudinal analysis in the context of ABMs and expand ABMs analyses relating to dynamic processes such as household-level vehicle transactions.

Demographic evolution↗
DOE OSTI2025

Jipole: A Differentiable ipole-based Code for Radiative Transfer in Curved Spacetimes

Recent imaging of supermassive black holes by the Event Horizon Telescope has relied on exhaustive parameter-space searches, matching observations to large, precomputed libraries of theoretical models. As observational data become increasingly precise, the limitations of this computationally expensive approach grow more acute, creating a pressing need for more efficient methods. In this work, we present Jipole, an automatically differentiable (AD), ipole-based code for radiative transfer in curved spacetimes, designed to compute image gradients with respect to underlying model parameters. These gradients quantify how parameter changes—such as the black hole’s spin or the observer’s inclination—affect the image, enabling more efficient parameter estimation and reducing the number of required images. We validate Jipole against ipole in two analytical tests and then compare pixelwise intensity derivatives from AD with those from finite-difference methods. We then demonstrate the utility of these gradients by performing parameter recovery for an analytical model in three increasingly complex cases for the injected image: ideal, blurred, and blurred with added noise. In most cases, high-accuracy fits are obtained in only a few optimization steps, failing only in cases with extremely low signal-to-noise ratios. These results highlight the potential of AD-based methods to accelerate robust, high-fidelity model-data comparisons in current and future black hole imaging efforts.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2002

Flexible Wing Base Micro Aerial Vehicles: Towards Flight Autonomy: Vision-Based Horizon Detection for Micro Air Vehicles

Recently substantial progress has been made towards design building and testifying remotely piloted Micro Air Vehicles (MAVs). This progress in overcoming the aerodynamic obstacles to flight at very small scales has, unfortunately, not been matched by similar progress in autonomous MAV flight. Thus, we propose a robust, vision-based horizon detection algorithm as the first step towards autonomous MAVs. In this paper, we first motivate the use of computer vision for the horizon detection task by examining the flight of birds (biological MAVs) and considering other practical factors. We then describe our vision-based horizon detection algorithm, which has been demonstrated at 30 Hz with over 99.9% correct horizon identification, over terrain that includes roads, buildings large and small, meadows, wooded areas, and a lake. We conclude with some sample horizon detection results and preview a companion paper, where the work discussed here forms the core of a complete autonomous flight stability system.

Nechyba, Michael C.↗
NASA NTRS1992

The psychophysical function for perceived gravitational-inertial force does not depend on the orientation of the otolith organs

It has generally been believed that the perceived intensity of a gravitational-inertial force depends on both the magnitude and orientation of the force with respect to the otolith organs, as does the elevator illusion. In this study, we examine the perceived intensity of Gz force and the elevator illusion as a function of the applied force and the orientation of the subject's head. Methods: Each of 7 male subjects was seated upright in a swinging chair mounted in the Ames 20-G Human Centrifuge while he set a visual target to his apparent horizon and judged the perceived intensity of Gz forces by cross-modal matches on a hand dynamometer. Plateau Gz levels were 1.00 1.25, 1.50, 2.00, 2.25, and 2.50; a 30 second ramp to plateau was used in all cases, and the duration of exposure at each plateau was 120 seconds. All measures were obtained both with the subject's head erect and pitched forward 30 degrees. Results: Although the elevator illusion changed with head orientation (F(6,60) = 7.56; p less than 0.001) the perceived intensity of Gz was essentially the same for both orientations of the head (F (6,60) = 0.61; p greater than 50). Conclusions: The results of this experiment suggest that the perceived intensity of gravitational-inertial force does not depend on otolith mechanisms in the same way as does the elevator illusion and that somesthetic, tactile, and other proprioceptive inputs are important for the psychophysical function.

Cohen, M.↗
DOE OSTI2025

Coupling Approaches with Non-matching Grids for Classical Linear Elasticity and Bond-based Peridynamic Models in 1D

Local-nonlocal coupling approaches provide a means to combine the computational efficiency of local models and the accuracy of nonlocal models. To facilitate the coupling of the two models, non-matching grids are often desirable as nonlocal grids usually require a finer resolution than local grids. In that case, it is often convenient to resort to interpolation operators so that models can exchange information in the overlap regions when nodes from the two grids do not coincide. This paper studies three existing coupling approaches, namely 1) a method that enforces matching displacements in an overlap region, 2) a variant that enforces a constraint on the stresses instead, and 3) a method that considers a variable horizon in the vicinity of the interfaces. Further, the effect of the interpolation order and of the grid ratio on the performance of the three coupling methods with non-matching grids is carefully studied on one-dimensional examples using polynomial manufactured solutions. The numerical results show that the degree of the interpolants should be chosen with care to avoid introducing additional modeling errors, or simply minimize these errors, in the coupling approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Gravitational index of the heterotic string

The fundamental heterotic string has a tower of BPS states whose supersymmetric index has an exponential growth in the charges. We construct the saddle-point of the gravitational path integral corresponding to this index. The saddle-point configuration is a supersymmetric rotating non-extremal Euclidean black hole. This configuration is singular in the two-derivative theory. We show that the addition of higher-derivative terms in four-dimensional N = 2 supergravity resolves the singularity. In doing so, we extend the recently-developed “new attractor mechanism” to include the effect of higher-derivative terms. Remarkably, the one-loop, four-derivative F-term contribution to the prepotential leads to a precise match of the gravitational and microscopic index. We also comment, using the effective theory near the horizon, on the possibility of a string-size near-extremal black hole. Our results clarify the meaning of different descriptions of this system in the literature. The thermal state transitions to a winding condensate and a gas of strings without ever reaching a small black hole, while the index is captured by the rotating Euclidean black hole solution and is constant and thus smoothly connected to the microscopic ensemble.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Quantum fate of the Choptuik naked singularity

Classical critical collapse provides a dynamical route from smooth initial data to a naked singularity, representing a sharper violation of predictability than ordinary black hole singularities. We argue that this distinction is erased by quantum backreaction. Building on the semiclassical interior analysis, where quantum self-energy of the collapsing matter generates a universal growing mode and a finite mass gap, we study the exterior naked singularity region that determines global visibility in the Einstein-scalar system. We analyze controlled exterior models in both $2+1$ and $3+1$ dimensions. In the former, smooth matching and physical boundary conditions analytically select a vacuum polarization state, whose backreaction cloaks the classically naked region by a quantum trapped branch. In the latter, numerical horizon tracing shows that near a quantum-shifted threshold the exterior develops finite-mass marginally trapped surfaces rather than a zero-mass naked endpoint. These results suggest a global quantum picture in which the Choptuik naked singularity shares the fate of an ordinary black hole singularity: quantum effects push the putative Cauchy horizon behind a quantum-generated horizon, thereby reducing the loss of predictability to the standard black hole evaporation problem.

FOS: Physical sciences↗
NASA NTRS2018

Remote Sensing of the Ocean Surface Refractive Index via Short-Wave Infrared Polarimetry

The linear polarization of sunlight reflected by ocean surfaces in the shortwave infrared (SWIR), at geometries where specular reflection dominates the signal, is a direct function of the refractive index of the surface microlayer(SML). This simple physical concept is at the base of a novel technique presented in this study. We invert observations obtained by the airborne Research Scanning Polarimeter (RSP) in the sun glint region, where each pixel's radiance is dominated by the signal originating from the wave slopes oriented precisely to cause specular reflection. The SWIR wavelength ensures minimization of aerosol interference when radiance travels through the atmosphere; strong absorption by the water body then limits the penetration depth to the first micrometer or so, effectively probing the SML. The resulting Degree of Linear Polarization (DoLP) is then governed by the refractive index via the Fresnel law for the specfic pixel geometry, independently of the wind speed. The selected dataset concerns several field deployments from both low- and high-altitude aircraft, including total reflectance measurements with the sole purpose of accounting for the residual aerosol effect.Stable retrievals from transects above pure seawater yielded values of refractive index that match the values published in the literature within an accuracy of 5 x 10 [superscript 4]. Flying over the oil spill caused by the explosion of the Deepwater Horizon platform, detected variations were found compatible with the presence of an oil slick. The robustness of the results, guaranteedby the high RSP polarimetric accuracy (less than or equal to 0:2%), opens the possibility for remote-sensing detection of other entities that similarly affect the refractive index including whitecaps, microplastics, biological gels, seaweed and grass mats.

Surface Microlayer↗
DOE OSTI2025

Enhancing Autonomous Control of Microreactors Using Multi-Agent Reinforcement Learning

In order for microreactors to be economically competitive, operation costs will need to be minimized through some degree of autonomous control. Previous work has demonstrated the effectiveness of reinforcement learning (RL) for load-following control in a drum-controlled microreactor. This study extends that work by exploring the potential of RL to independently control each of the reactor’s drums. We compare a single-agent RL approach with a multi-agent RL (MARL) framework, testing them for generalization across different load-following power profiles and control timescales, and for robustness in cases of randomly disabled control drums. Since the point kinetics simulation environment used in this study cannot resolve spatial effects, we assume that in the absence of spatially localized disturbances, optimal drum movements should be symmetrical. We demonstrate that single-agent RL is able to achieve accurate performance only when symmetric actions are ignored; otherwise, it fails to train a useful controller. Meanwhile, the MARL framework performs symmetric actions by design and trains a robust, accurate agent, as evidenced by mean absolute errors in power matching of 0.41% for the training power profile, 0.68% for a profile with half the drums disabled, and 0.21% for a profile on a realistic load-following time horizon.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRSDate not supplied

An Updated Secondary Lunar Meteoroid Ejecta Model for Engineering Design

Introduction: The surface of the Moon is constantly being bombarded by a flux of meteoroids of various sizes. Impacts due to these meteoroids produce secondary ejecta material at much lower speeds but with a total mass larger than the original impactor. Details about the secondary ejecta are important for planning missions on the lunar surface. In this work, an updated ejecta model is presented called the Meteoroid Model of Secondary Ejecta (MeMoSeE), to re-place the Apollo-era ejecta model, NASA SP-8013 [1], in the SLS-SPEC 159 Design Specification for Natural Environments (DSNE) [2]. The model pro-duces secondary ejecta flux environments for a user-specified location on the lunar surface, and sorts the incoming secondary flux by angular direction and speed. Methods: MeMoSeE is separated into three parts: the inputs, the conversion step, and the integration of fluxes. Inputs to the model include the primary meteoroid fluxes and the primary near-Earth object (NEO) fluxes. Meteoroid fluxes, both asteroidal and cometary, are calculated using the Meteoroid Engineering Model (MEM3) [3] for different locations on the Moon. For each surface location, an ephemeris is generated using the JPL HORIZONS System [4] that feeds into MEM3. The NEO fluxes are approximated by the high-density population of MEM3 (i.e., only the directionality), where the speed distribution of the NEO fluxes is renormalized to match observations [5]. The regolith properties are used as defined in the DSNE [2]. The conversion step utilizes scaling laws given by Housen & Holsapple 2011 [6] to convert the primary impactor flux to the total mass of secondary ejecta. The ejecta distribution, at the point-of-impact (POI), is separated into a zenith angle and azimuthal angle distribution. The zenith angle distribution follows a beta distribution where the peak depends on the impact altitude angle and the impact azimuth [7]. We employ an ejecta azimuth distribution that is based on Rival & Mandeville 1999 [8] which focuses ejecta in the downstream direction for more oblique impacts. Finally, during the integration step, we sum secondary ejecta number fluxes at a particular region-of-interest (ROI) that originated from many POI locations over the entire surface of the Moon. We keep track of both altitude and azimuth angle bins as well as a range of speed bins, following the igloo gridding as done in MEM3 [9]. The ejecta particle size distribution and density is assumed to be the same as the lunar regolith [10, 2]. Results: The primary fluxes are computed for one Metonic cycle (19 years) for various locations over the lunar surface with a fixed orientation. Both the angular and speed distributions of the primary fluxes are dependent on the latitude and longitude. The over-all primary fluxes show a roughly 13% increase from the eastern limb to the western limb. In general, the speed distribution of the secondary fluxes span from a user-defined minimum speed to the escape speed of the Moon (2.38 km/s), roughly following a power-law relation [6]. Different parts of the speed distribution come from different primary impact locations on the Moon. The secondary ejecta is dominated by the slowest speeds, where these particles originate nearby the ROI. For speeds around 71% of the escape speed, the secondary ejecta originates from locations near the antipodal point. On the other hand, for secondary ejecta speeds that exceed roughly 90% the escape speed, the ejecta particles come from all over the lunar surface to the ROI. Comparing the secondary ejecta fluxes from MeMoSeE with NASA SP-8013 [1], there is a reduction by about 2-3 orders of magnitude for secondary ejecta particles greater than 1 µg. These estimates agree with recent findings from Bjorkman & Christiansen 2019 [11]. The secondary ejecta fluxes are also compared with the primary fluxes, where the ejecta fluxes are roughly an order of magnitude greater than the primary fluxes. References: [1] Cour-Palais, B. G., (1969) NASA SP-8013. [2] NASA SLS-SPEC-159 Rev. H (2020). [3] Moorhead, A. V., et al. (2019) JS&R, 1-17. [4] Giorgini, J. D., (2015) IAUGA, 29, 2256293. [5] Moorhead, A. V., (2020) Memo OSMA/MEO/Lunar-001. [6] Housen, K. R., and Holsapple, K. A., (2011) Icarus, 211(1), 856-875. [7] Gault, D. E., and Wedekind, J. A., (1978) L&PSCP, 9, 3843-3875. [8] Rival, M., and Mandeville, J., (1999) Space Debris, 1(1), 45-57. [9] Moorhead, A. V., (2019) MEM3 User Guide. [10] Carrier, W. D., (1973) The Moon, 6(3-4), 250-263. [11] Bjorkman, M. D. & Christiansen, E. L., (2019) ODC

impacts↗
NASA NTRS1977

Ultraviolet spectroscopy of the zodiacal light at 20-deg elongation

The zodiacal light at 20-deg elongation and 10-deg inclination was observed by rocket ultraviolet spectrometers at 10-15-A resolution in the spectral range 1200-3200 A during an experiment designed to observe comet Kohoutek (1973 XII). The data were obtained above 180 km when scattered horizon light in the startracker caused a loss of tracking on the comet. Airglow emission due to NO and O(+), identified spectroscopically and by its variation with altitude, is significant between 1900 and 2500 A. Longward of 2600 A, the spectrum matches that of the sun, and the derived value of the color ratio, relative to the visible, is 0.90 + or - 0.20. At 1600 A, an upper limit on the zodiacal-light emission of 0.07 R per A or 7 hundred-millionths erg/s per sq cm/sterad per A is obtained.

Feldman, P. D.↗