Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial constraints”

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 649 records · Page 36

Numerical Investigation of Film Coefficient Approximation for Chemically Reacting Boundary-Layer Flows

Aerothermal analysis of spacecraft planetary entry is heavily dependent on heritage engineering models. The film coefficient heat transfer model examined in this paper estimates the convective heating to the vehicle for a laminar, dissociated, chemically reacting boundary layer for an Earth atmosphere. This model requires information about the vehicle and flowfield for a given trajectory point and estimates a proportional relationship between enthalpy potential and convective heat flux. In practice it is the aerothermal engineer who must decide which assumptions are appropriate for his/her application. This work looks at numerous CFD simulations for an arbitrary, axisymmetric flight vehicle to analyze the relative importance of both the mass and energy constraints imposed at the wall boundary, as well as the effect of various diffusion models. Within the subset of tested energy boundary conditions, it is found that the most desirable energy boundary condition is the radiative equilibrium boundary condition, which permits conservative estimates of convective heat flux, but also generates flowfield-dependent spatial thermal distributions along the surface. Other key findings are presented in an effort to make the film coefficient engineering model readily available to design engineers across industry.

Film coefficient↗

New Constraints on the Volatile Deposit in Mercury’s North Polar Crater, Prokofiev

There is strong evidence from multi-wavelength, multi-instrument observations and thermal models for the presence of ice and other volatiles in Mercury’s north polar permanently shadowed regions (PSRs)[1-6]. Ground-based radar measurements identified bright regions in some polar craters interpreted to be water ice on or near the surface[1]. The MESSENGER spacecraft found that some radar-bright deposits appeared bright in the optical and near-infrared while others appeared dark[4-6]. The bright regions were interpreted to be exposed surface ice while the dark regions were hypothesized to be complex organic volatile lag deposits overlaying buried ice[3]. One example of the former is in the 112 km diameter crater, Prokofiev, for which the Mercury Dual Imaging System (MDIS) acquired broadband images while in Sun shadow but illuminated by scattered light from the crater walls[5]. These images revealed a relatively high-albedo region spatially correlated with the PSR. Subsequent work[7] found a small offset between the boundary of the radar-bright region and that of the PSR predicted by an illumination model based on the polar digital elevation model (DEM)from the Mercury Laser Altimeter (MLA). Such an offset, if true, could have implications for the ice properties and its delivery and destruction mechanisms. However, the low resolution (1 km/pix) of the polar MLA DEM precluded a definitive conclusion. In the present work, we build new higher-resolution topographic, illumination, and thermal models of Prokofiev to place stronger constraints on its PSR and volatile deposits. We also study the MLA reflectance data in this crater to quantitatively estimate the ice abundance.

Michael K Barker↗

Space-time modeling using environmental constraints in a mobile robot system

Grid-based models of a robot's local environment have been used by many researchers building mobile robot control systems. The attraction of grid-based models is their clear parallel between the internal model and the external world. However, the discrete nature of such representations does not match well with the continuous nature of actions and usually serves to limit the abilities of the robot. This work describes a spatial modeling system that extracts information from a grid-based representation to form a symbolic representation of the robot's local environment. The approach makes a separation between the representation provided by the sensing system and the representation used by the action system. Separation allows asynchronous operation between sensing and action in a mobile robot, as well as the generation of a more continuous representation upon which to base actions.

Slack, Marc G.↗

A Decade of Global Volcanic SO2 Emissions Measured from Space

The global flux of sulfur dioxide (SO2) emitted by passive volcanic degassing is a key parameter that constrains the fluxes of other volcanic gases (including carbon dioxide, CO2) and toxic trace metals (e.g., mercury). It is also a required input for atmospheric chemistry and climate models, since it impacts the tropospheric burden of sulfate aerosol, a major climate-forcing species. Despite its significance, an inventory of passive volcanic degassing is very difficult to produce, due largely to the patchy spatial and temporal coverage of ground-based SO2 measurements. We report here the first volcanic SO2 emissions inventory derived from global, coincident satellite measurements, made by the Ozone Monitoring Instrument (OMI) on NASA's Aura satellite in 2005-2015. The OMI measurements permit estimation of SO2 emissions from over 90 volcanoes, including new constraints on fluxes from Indonesia, Papua New Guinea, the Aleutian Islands, the Kuril Islands and Kamchatka. On average over the past decade, the volcanic SO2 sources consistently detected from space have discharged a total of approximately 63 kt/day SO2 during passive degassing, or approximately 23 +/- 2 Tg/yr. We find that approximately 30% of the sources show significant decadal trends in SO2 emissions, with positive trends observed at multiple volcanoes in some regions including Vanuatu, southern Japan, Peru and Chile.

Carn, S. A.↗

Improved constraints on hematite refractive index for estimating climatic effects of dust aerosols

Abstract Uncertainty in desert dust composition poses a big challenge to understanding Earth’s climate across different epochs. Of particular concern is hematite, an iron-oxide mineral dominating the solar absorption by dust particles, for which current estimates of absorption capacity vary by over two orders of magnitude. Here, we show that laboratory measurements of dust composition, absorption, and scattering provide valuable constraints on the absorption potential of hematite, substantially narrowing its range of plausible values. The success of this constraint is supported by results from an atmospheric transport model compared with station-based measurements. Additionally, we identify substantial bias in simulating hematite abundance in dust aerosols with current soil mineralogy descriptions, underscoring the necessity for improved data sources. Encouragingly, the next-generation imaging spectroscopy remote sensing data hold promise for capturing the spatial variability of hematite. These insights have implications for enhancing dust modeling, thus contributing to efforts in climate change mitigation and adaptation.

Environmental Sciences & Ecology↗

Multiscale and multidimensional modeling of particle acceleration and transport in solar flares

Multi-messenger, multi-viewpoint, and time-resolved observations of solar flares are now providing unprecedented constraints on particle acceleration sites, energy conversion, and energy transport. The interpretation of current observations, including microwave imaging spectroscopy from EOVSA, hard x-ray (HXR) imaging from Solar Orbiter/STIX, gamma-ray diagnostics from Fermi, and in situ measurements from Parker Solar Probe and Solar Orbiter, collectively demands modeling frameworks that go beyond traditional spatially unresolved, one-zone models or single-mechanism descriptions. This review surveys multiscale and multidimensional modeling approaches, including kinetic, magnetohydrodynamic (MHD), and macroscopic particle models, that are being developed to meet the need. Kinetic simulations reveal that three-dimensional (3D) effects, including field-line chaos and self-generated turbulence, are essential for sustained power-law particle acceleration. MHD simulations now capture flux-rope eruptions, plasmoid-unstable current sheets, and turbulent flare regions in realistic magnetic topologies. Macroscopic models coupling MHD with energetic-particle models produce spatially resolved electron distributions and synthetic HXR and microwave emissions for direct comparison with observations. Despite these advances, outstanding challenges remain in bridging kinetic and global scales, improving MHD simulations and macroscopic particle models, and achieving quantitative model-observation closure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Unified Model of Cosmic-Ray Propagation and Radio Extreme Scattering Events from Intermittent Interstellar Structures

Intermittent magnetic structures are a plausible candidate for explaining cosmic-ray (CR) diffusion rates derived from observed CR energy spectra. Independently, studies of extreme scattering events (ESEs) of radio quasars and pulsar scintillation have hinted that very straight, large aspect ratio magnetic current sheets may be responsible for the localized large scattering of radio waves. The required shortest axis of the typical structures producing ESEs is of the same scale (∼au) as the gyroradii of ∼GeV CRs. In this Letter, we propose that the same magnetic/density sheets can produce large scattering of both CRs and radio waves. We demonstrate that the geometry and volume-filling factor of the sheets derived from quasar ESEs can explain the observed mean free path of GeV CRs without introducing free parameters. The model places constraints on the sheet geometry, such as straightness and large aspect ratio, and assumes that the statistics of the sheets are similar throughout the Galactic volume. We therefore discuss observational tests of the sheet model, which includes observations of echoes in pulsars and fast radio bursts, gravitationally lensed quasars, the distribution of ESE durations, and spatial correlations between ESE events and rotation measure fluctuations. Such tests will be enabled by upcoming wide-field radio instruments, including the Canadian Hydrogen Observatory and Radio-transient Detector and Deep Synoptic Array 2000 Antennas.

Kempski, Philipp [Princeton Univ., NJ (United Stat↗

System implications of large radiometric array antennas

Current radiometric earth and atmospheric sensing systems in the centimeter wavelength range generally employ a directive antenna connected through a single terminal pair to a Dicke receiver. It is shown that this approach does not lend itself to systems with greatly increased spatial resolution. Signal to noise considerations relating to antenna efficiency force the introduction of active elements at the subarray level; thus, if Dicke switching is to be used, it must be distributed throughout the system. Some possible approaches are suggested. The introduction of active elements at the subarray level is found to ease the design constraints on time delay elements, necessary for bandwidth, and on multiple beam generation, required in order to achieve sufficient integration time with high resolution.

Levis, C. A.↗

Optical-system design for next-generation pushbroom sensors

Next-generation pushbroom sensors for earth observation require high-performance optics that provide high spatial resolution over wide fields of view. Specifically, blur diameters on the order of 10 to 15 microns are needed over 5 to 15 deg fields. In addition to this fundamental level of optical performance, other characteristics, such as spatial coregistration of spectral bands, flat focal plane, telecentricity, and workable pupil location are significant instrument design considerations. The detector-assembly design, optical line-of-sight pointing method and sensor packaging all hinge on these secondary attributes. Moreover, the need for broad spectral coverage, ranging from 0.4 to 12.5 microns, places an additional constraint on optical design. This paper presents alternative design forms that are candidates for wide-field pushbroom sensors, and discusses the instrument-design tradeoffs that are linked to the selection of these alternate optical approaches.

Mika, A. M.↗

Imaging the sun in hard X-rays - Spatial and rotating modulation collimators

Several approaches to imaging hard X-rays emitted from solar flares have been proposed or are planned for the nineties including the spatial modulation collimator (SMC) and the rotating modulation collimator (RMC). A survey of current solar flare theoretical literature indicates the desirability of spatial resolutions down to 1 arcsecond, field of views greater than the full solar disk (i.e., 32 arcminutes), and temporal resolutions down to 1 second. Although the sun typically provides relatively high flux levels, the requirement for 1 second temporal resolution raises the question as to the viability of Fourier telescopes subject to the aforementioned constraints. A basic photon counting, Monte Carlo 'end-to-end' model telescope was employed using the Astronomical Image Processing System (AIPS) for image reconstruction. The resulting solar flare hard X-ray images compared against typical observations indicated that both telescopes show promise for the future.

Campbell, Jonathan W.↗

Comet Shoemaker-Levy 9: An Active Comet

The important elements of the debate over the activity versus dormancy of comet Shoemaker Levy 9 (S-L 9) are reviewed. It is argued that the circularity of the isophotes in the inner comae of S-L 9 as well as the spatial dependencies of the comae brightness profiles are indicators of sustained dust production by S-L 9. It is also shown that the westward tail orientations, which were formerly interpreted as a sign of the comet's dormancy, are not a good indicator of either activity or dormancy. Rather, the tail orientations simply place constraints on the dust production rate for grains smaller than approx. equals 5(micron). All the available evidence points to S-L 9 as having been an active, dust-producing comet. Synthetic images of an active comet are fitted to Hubble Space Telescope images of the S-L 9 fragment K, and its grain size and outflow velocity distributions are extracted. These findings show that the appearance of the dust coma was dominated by large grains having radii between approx. equals 30 (micron) and approx. equals 3 mm, produced at a rate of M approx. equals 22 kg/ s, and ejected at outflow velocities of approx. equals 0.5 m/ s. Only upper limits on the production rates of smaller grains are obtained. The nucleus of fragment K was not observed directly but its size is restricted to lie within a rather narrow interval 0.4 less than or equal to Rf less than or equal to 1.2 km.

Rettig, Terrence W.↗

Super Smooth Optics for Extra-Solar Planet Detection

The goal of imaging planets around the nearby stars has important scientific significance but requires the use of advanced methods of controlling diffracted and scattered light. Over the last three years we have undertaken a study of coronagraphic methods of controlling diffracted light and of figuring hyper-contrast optics. Progress in these two general areas have led to a proposed space-based, 1.9 meter diameter coronagraphic telescope designed specifically for very high performance in the imaging of faint objects near bright sources. This instrument, called the Circumstellar Imaging Telescope (CIT), relies on a new high efficiency coronagraph design and the careful control of scattered light by extremely smooth optics. The high efficiency coronagraph uses focal plane apodization in order to concentrate diffracted light more efficiently in the pupil. This allows convenient removal of the diffracted light by masking off parts of the telescope pupil while not sacrificing the center of the field. Reductions of diffracted light by factors exceeding 1000 are not only possible but are required in order to detect extra-solar planets. Laboratory experiments with this new design have confirmed the theoretical diffraction reductions to the limits of the optics used (factors of about 300) . The extremely high efficiency of this coronagraph puts strong constraints on the narrow angle scattered light due to figure errors in the telescope mirror. Since planets orbiting nearby stars are expected at angular distances of about 1 arcsecond, it is in this small angular range in which scattering must be controlled. The figure errors responsible for scattering in this range come from mid-spatial frequencies corresponding to correlation lengths of about 10 cm on the primary mirror. A primary mirror about 15 times smoother than the Hubble Space Telescope mirror is required for the CIT. Laboratory experiments indicate that small test mirrors can be fabricated with existing technology which come within a factor of two of this requirement.

Terrile, Richard J.↗

Initial OSSE Sky Survey Observation: Deep Survey in the Virgo Region

The OSSE Virgo survey covered about 1% of the sky, using eight weeks of observing time on the Compton Gamma Ray Observatory, and was performed in conjunction with the EGRET and COMPTEL instruments on CGRO. A fundamental aspect of the survey was the development of data reduction techniques for the scanning data from the OSSE detectors. A direct algebraic approach for spatially deconvolving the data was developed. For the survey work, no background observations isolated from the survey region were acquired, unlike normal discrete source observations. The survey, reduction method covers the mapped region with a grid of pixels, and estimates the flux from each pixel by least-squares agreement with the data, assuming positivity constraints. This technique has led to the detection of compact sources at hard X-ray energies within the Virgo region, and the OSSE survey studies have been extended to the south galactic pole and north Ecliptic pole, with similar analysis techniques.

Cameron, Robert↗

Kinematic principles of primate rotational vestibulo-ocular reflex. II. Gravity-dependent modulation of primary eye position

The kinematic constraints of three-dimensional eye positions were investigated in rhesus monkeys during passive head and body rotations relative to gravity. We studied fast and slow phase components of the vestibulo-ocular reflex (VOR) elicited by constant-velocity yaw rotations and sinusoidal oscillations about an earth-horizontal axis. We found that the spatial orientation of both fast and slow phase eye positions could be described locally by a planar surface with torsional variation of <2.0 +/- 0.4 degrees (displacement planes) that systematically rotated and/or shifted relative to Listing's plane. In supine/prone positions, displacement planes pitched forward/backward; in left/right ear-down positions, displacement planes were parallel shifted along the positive/negative torsional axis. Dynamically changing primary eye positions were computed from displacement planes. Torsional and vertical components of primary eye position modulated as a sinusoidal function of head orientation in space. The torsional component was maximal in ear-down positions and approximately zero in supine/prone orientations. The opposite was observed for the vertical component. Modulation of the horizontal component of primary eye position exhibited a more complex dependence. In contrast to the torsional component, which was relatively independent of rotational speed, modulation of the vertical and horizontal components of primary position depended strongly on the speed of head rotation (i.e., on the frequency of oscillation of the gravity vector component): the faster the head rotated relative to gravity, the larger was the modulation. Corresponding results were obtained when a model based on a sinusoidal dependence of instantaneous displacement planes (and primary eye position) on head orientation relative to gravity was fitted to VOR fast phase positions. When VOR fast phase positions were expressed relative to primary eye position estimated from the model fits, they were confined approximately to a single plane with a small torsional standard deviation ( approximately 1.4-2.6 degrees). This reduced torsional variation was in contrast to the large torsional spread (well >10-15 degrees ) of fast phase positions when expressed relative to Listing's plane. We conclude that primary eye position depends dynamically on head orientation relative to space rather than being fixed to the head. It defines a gravity-dependent coordinate system relative to which the torsional variability of eye positions is minimized even when the head is moved passively and vestibulo-ocular reflexes are evoked. In this general sense, Listing's law is preserved with respect to an otolith-controlled reference system that is defined dynamically by gravity.

NASA Discipline Neuroscience↗

GBaTSv2: a revised synthesis of the likely basal thermal state of the Greenland Ice Sheet

The basal thermal state (frozen or thawed) of the Greenland Ice Sheet is under-constrained due to few direct measurements, yet knowledge of this state is becoming increasingly important to interpret modern changes in ice flow. The first synthesis of this state relied on inferences from widespread airborne and satellite observations and numerical models, for which most of the underlying datasets have since been updated. Further, new and independent constraints on the basal thermal state have been developed from analysis of basal and englacial reflections observed by airborne radar sounding. Here we synthesize constraints on the Greenland Ice Sheet's basal thermal state from boreholes, thermomechanical ice-flow models that participated in the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6; Coupled Model Intercomparison Project Phase 6), IceBridge BedMachine Greenland v4 bed topography, Making Earth Science Data Records for Use in Research Environments (MEaSUREs) Multi-Year Greenland Ice Sheet Velocity Mosaic v1 and multiple inferences of a thawed bed from airborne radar sounding. Most constraints can only identify where the bed is likely thawed rather than where it is frozen. This revised synthesis of the Greenland likely Basal Thermal State version 2 (GBaTSv2) indicates that 33 % of the ice sheet's bed is likely thawed, 40 % is likely frozen and the remainder (28 %) is too uncertain to specify. The spatial pattern of GBaTSv2 is broadly similar to the previous synthesis, including a scalloped frozen core and thawed outlet-glacier systems. Although the likely basal thermal state of nearly half (46 %) of the ice sheet changed designation, the assigned state changed from likely frozen to likely thawed (or vice versa) for less than 6 % of the ice sheet. This revised synthesis suggests that more of northern Greenland is likely thawed at its bed and conversely that more of southern Greenland is likely frozen, both of which influence interpretation of the ice sheet's present subglacial hydrology and models of its future evolution. The GBaTSv2 dataset, including both code that performed the analysis and the resulting datasets, is freely available at https://doi.org/10.5281/zenodo.6759384 (MacGregor, 2022).

Joseph A. MacGregor↗

Dynamic Learning of Correlation Potentials for a Time-Dependent Kohn-Sham System

We develop methods to learn the correlation potential for a time-dependent Kohn-Sham (TDKS) system in one spatial dimension. We start from a low-dimensional two-electron system for which we can numerically solve the time-dependent Schr¨odinger equation; this yields electron densities suitable for training models of the correlation potential. We frame the learning problem as one of optimizing a least-squares objective subject to the constraint that the dynamics obey the TDKS equation. Applying adjoints, we develop efficient methods to compute gradients and thereby learn models of the correlation potential. Our results show that it is possible to learn values of the correlation potential such that the resulting electron densities match ground truth densities. We also show how to learn correlation potential functionals with memory, demonstrating one such model that yields reasonable results for trajectories outside the training set.

97 MATHEMATICS AND COMPUTING↗

Constraints on the Recent Rate of Lunar Regolith Accumulation from Diviner Observations

Many large craters on the lunar nearside show radar CPR signatures consistent with the presence of blocky ejecta blankets, to distances pre dicted to be covered by continuous ejecta. However, most of these sur faces show limited enhancements in both derived rock abundance and rock-free regolith temperatures calculated from Diviner nighttime infrar ed observations. This indicates that the surface rocks are covered by a layer of thermally insulating regolith material. By matching the results of one-dimensional thermal models to Diviner nighttime temperat ures, we have constrained the thermophysical properties of the upper regolith, and the thickness of regolith overlying proximal ejecta. We find that for all of the regions surveyed (all in the nearside highla nds), the nighttime cooling curves are best fit by a density profile that varies exponentially with depth, consistent with a linear mixture of rocks and regolith fines, with increasing rock content with depth . Our results show significant spatial variations in the density e-folding depth, H, among young crater ejecta regions, indicating differen ces in the thickness of accumulated regolith. However, away from youn g craters, the average regional "equilibrium" value of H (Heq) is remarkably consistent, and is on the order of 5 cm. As expected, near-rim ejecta associated with young craters show lower values of H, indicating a high rock content in the shallow subsurface; for older craters, the average value of H approaches the regional value of Heq. Calculat ed H values for young craters show a clear correlation with published ages, providing the first observational constraint on the recent rate of lunar regolith accumulation. In addition, this result may help to resolve the apparent discrepancy between ages calculated from small crater counts on melt ponds versus counts on continuous ejecta (e.g., King crater; Ashley et al., 2011, LPSC 42, abstract 2437). This method could, in principle, be extended to other airless bodies (e.g., aste roids), which would in turn constrain the recent impactor flux.

Ghent, R. R.↗

A solution framework for linear PDE-constrained mixed-integer problems

Abstract We present a general numerical solution method for control problems with state variables defined by a linear PDE over a finite set of binary or continuous control variables. We show empirically that a naive approach that applies a numerical discretization scheme to the PDEs to derive constraints for a mixed-integer linear program (MILP) leads to systems that are too large to be solved with state-of-the-art solvers for MILPs, especially if we desire an accurate approximation of the state variables. Our framework comprises two techniques to mitigate the rise of computation times with increasing discretization level: First, the linear system is solved for a basis of the control space in a preprocessing step. Second, certain constraints are just imposed on demand via the IBM ILOG CPLEX feature of a lazy constraint callback. These techniques are compared with an approach where the relations obtained by the discretization of the continuous constraints are directly included in the MILP. We demonstrate our approach on two examples: modeling of the spread of wildfire and the mitigation of water contamination. In both examples the computational results demonstrate that the solution time is significantly reduced by our methods. In particular, the dependence of the computation time on the size of the spatial discretization of the PDE is significantly reduced.

97 MATHEMATICS AND COMPUTING↗