Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “histograms”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
NASA NTRS2016

Full Field X-Ray Fluorescence Imaging Using Micro Pore Optics for Planetary Surface Exploration

Many planetary surface processes leave evidence as small features in the sub-millimetre scale. Current planetary X-ray fluorescence spectrometers lack the spatial resolution to analyse such small features as they only provide global analyses of areas greater than 100 mm(exp 2). A micro-XRF spectrometer will be deployed on the NASA Mars 2020 rover to analyse spots as small as 120m. When using its line-scanning capacity combined to perpendicular scanning by the rover arm, elemental maps can be generated. We present a new instrument that provides full-field XRF imaging, alleviating the need for precise positioning and scanning mechanisms. The Mapping X-ray Fluorescence Spectrometer - "Map-X" - will allow elemental imaging with approximately 100μm spatial resolution and simultaneously provide elemental chemistry at the scale where many relict physical, chemical and biological features can be imaged in ancient rocks. The arm-mounted Map-X instrument is placed directly on the surface of an object and held in a fixed position during measurements. A 25x25 mm(exp 2) surface area is uniformly illuminated with X-rays or alpha-particles and gamma-rays. A novel Micro Pore Optic focusses a fraction of the emitted X-ray fluorescence onto a CCD operated at a few frames per second. On board processing allows measuring the energy and coordinates of each X-ray photon collected. Large sets of frames are reduced into 2d histograms used to compute higher level data products such as elemental maps and XRF spectra from selected regions of interest. XRF spectra are processed on the ground to further determine quantitative elemental compositions. The instrument development will be presented with an emphasis on the characterization and modelling of the X-ray focussing Micro Pore Optic. An outlook on possible alternative XRF imaging applications will be discussed.

micro-XRF spectrometer↗
NASA NTRS2016

Post-Test Inspection of NASA's Evolutionary Xenon Thruster Long-Duration Test Hardware: Discharge Chamber

The NEXT Long-Duration Test is part of a comprehensive thruster service life assessment intended to demonstrate overall throughput capability, validate service life models, quantify wear rates as a function of time and operating condition, and identify any unknown life-limiting mechanisms. The test was voluntarily terminated in February 2014 after demonstrating 51,184 hours of high-voltage operation, 918 kg of propellant throughput, and 35.5 MN-s of total impulse. The post-test inspection of the thruster hardware began shortly afterwards with a combination of non-destructive and destructive analysis techniques, and is presently nearing completion. This paper presents relevant results of the post-test inspection for the discharge chamber as well as other miscellaneous components such as the high-voltage propellant isolators and electrical cabling. Comparison of magnetic field measurements taken during pretest and post-test inspections indicate that the field strength did not degrade, consistent with performance data obtained during the test. Inspection of discharge chamber mesh samples show a deposition coating primarily composed of grid material that is approximately 15 micrometers in thickness. This thickness is well within the retention capability of the mesh and is therefore not expected to present any issues. Approximately 3.1 grams of deposition flakes were found at the bottom of the discharge chamber, composed primarily of grid material and carbon. Calculated size histograms of these flakes indicate that 99% have a maximum dimension of 200 micrometers or smaller, which is significantly less than the ion optics grid gap. Larger flakes that are capable of causing a grid-to-grid short will be analyzed to determine if their formation will occur in flight or is a facility effect. The high-voltage propellant isolators as well as numerous other electrical insulators were inspected and no evidence of arcing or any other issues were found.

Electric Propulsion↗
NASA NTRS2016

Survey of the Frequency Dependent Latitudinal Distribution of the Fast Magnetosonic Wave Mode from Van Allen Probes Electric and Magnetic Field Instrument and Integrated Science Waveform Receiver Plasma Wave Analysis

We present a statistical survey of the latitudinal structure of the fast magnetosonic wave mode detected by the Van Allen Probes spanning the time interval of 21 September 2012 to 1 August 2014. We show that statistically, the latitudinal occurrence of the wave frequency (f) normalized by the local proton cyclotron frequency (f(sub cP)) has a distinct funnel-shaped appearance in latitude about the magnetic equator similar to that found in case studies. By comparing the observed E/B ratios with the model E/B ratio, using the observed plasma density and background magnetic field magnitude as input to the model E/B ratio, we show that this mode is consistent with the extra-ordinary (whistler) mode at wave normal angles (theta(sub k)) near 90 deg. Performing polarization analysis on synthetic waveforms composed from a superposition of extra-ordinary mode plane waves with theta(sub k) randomly chosen between 87 and 90 deg, we show that the uncertainty in the derived wave normal is substantially broadened, with a tail extending down to theta(sub k) of 60 deg, suggesting that another approach is necessary to estimate the true distribution of theta(sub k). We find that the histograms of the synthetically derived ellipticities and theta(sub k) are consistent with the observations of ellipticities and theta(sub k) derived using polarization analysis.We make estimates of the median equatorial theta(sub k) by comparing observed and model ray tracing frequency-dependent probability occurrence with latitude and give preliminary frequency dependent estimates of the equatorial theta(sub k) distribution around noon and 4 R(sub E), with the median of approximately 4 to 7 deg from 90 deg at f/f(sub cP) = 2 and dropping to approximately 0.5 deg from 90 deg at f/f(sub cP) = 30. The occurrence of waves in this mode peaks around noon near the equator at all radial distances, and we find that the overall intensity of these waves increases with AE*, similar to findings of other studies.

Boardsen, Scott A.↗
NASA NTRS2016

Prospects of the ICESat-2 Laser Altimetry Mission for Savanna Ecosystem Structural Studies Based on Airborne Simulation Data

The next planned spaceborne lidar mission is the Ice, Cloud and land Elevation Satellite 2 (ICESat-2), which will use the Advanced Topographic Laser Altimeter System (ATLAS) sensor, a photon counting technique. To pre-validate the capability of this mission for studying three dimensional vegetation structure in savannas, we assessed the potential of the measurement approach to estimate canopy height in an oak savanna landscape. We used data from the Multiple Altimeter Beam Experimental Lidar (MABEL), an airborne photon counting lidar sensor developed by NASA's Goddard Space Flight Center. ATLAS-like data was generated using the MATLAS simulator, which adjusts MABEL data's detected number of signal and noise photons to that expected from the ATLAS instrument. Transects flown over the Tejon ranch conservancy in Kern County, California, USA were used for this work. For each transect we chose to use data from the near infrared channel that had the highest number of photons. We segmented each transect into 50 m, 25 m and 14 m long blocks and aggregated the photons in each block into a histogram based on their elevation values. We then used an automated algorithm to identify cut off points where the cumulative density of photons from the highest elevation indicates the presence of the canopy top and likewise where such cumulative density from the lowest elevation indicates the mean terrain elevation. MABEL derived height metrics were moderately correlated to discrete return lidar (DRL) derived height metrics r(sub 2) and RMSE values ranging from 0.60 to 0.73 and 2.9 m to 4.4 m respectively) but MATLAS simulation resulted in more modest correlations with DRL indices r(sub 2) ranging from 0.5 to 0.64 and RMSE from 3.6 m to 4.6 m). Simulations also indicated that the expected number of signal photons from ATLAS will be substantially lower, a situation that reduces canopy height estimation precision especially in areas of low density vegetation cover. On the basis of the simulated data, there is reason to believe that the ability of ICESat-2 to estimate height in savannas will be comparable to the original ICESat mission although the respective sensors have different measurement principles.

Gwenzi, David↗
NASA NTRS2016

Regime-Based Evaluation of Cloudiness in CMIP5 Models

The concept of Cloud Regimes (CRs) is used to develop a framework for evaluating the cloudiness of 12 fifth Coupled Model Intercomparison Project (CMIP5) models. Reference CRs come from existing global International Satellite Cloud Climatology Project (ISCCP) weather states. The evaluation is made possible by the implementation in several CMIP5 models of the ISCCP simulator generating for each gridcell daily joint histograms of cloud optical thickness and cloud top pressure. Model performance is assessed with several metrics such as CR global cloud fraction (CF), CR relative frequency of occurrence (RFO), their product (long-term average total cloud amount [TCA]), cross-correlations of CR RFO maps, and a metric of resemblance between model and ISCCP CRs. In terms of CR global RFO, arguably the most fundamental metric, the models perform unsatisfactorily overall, except for CRs representing thick storm clouds. Because model CR CF is internally constrained by our method, RFO discrepancies yield also substantial TCA errors. Our findings support previous studies showing that CMIP5 models underestimate cloudiness. The multi-model mean performs well in matching observed RFO maps for many CRs, but is not the best for this or other metrics. When overall performance across all CRs is assessed, some models, despite their shortcomings, apparently outperform Moderate Resolution Imaging Spectroradiometer (MODIS) cloud observations evaluated against ISCCP as if they were another model output. Lastly, cloud simulation performance is contrasted with each model's equilibrium climate sensitivity (ECS) in order to gain insight on whether good cloud simulation pairs with particular values of this parameter.

cloud climatology↗
NASA NTRS2017

The Resource Prospector Neutron Spectrometer System: RP's Bloodhound

The primary goal of the Resource Prospector Neutron Spectrometer System (NSS) is to locate and characterize hydrogen-bearing volatile deposits, especially subsurface ice, that may exist at the lunar poles. A key objective is to detect water-equivalent hydrogen concentrations of 0.5 wt% or greater while on a moving rover. A second objective is to determine approximate burial depth of enhanced hydrogen-bearing materials up to 1 meter below otherwise dry regolith. The instrument will be carried aboard a landed mobility system at the lunar poles. The instrument operates by measuring the changes in the leakage flux of low energy neutrons out of the regolith. These neutrons are produced by galactic cosmic rays, which are so energetic that they shatter the nuclei in surface materials. The neutrons interact with other nuclei and lose energy, becoming thermalized in the process. Hydrogen is most efficient at thermalizing neutrons owing to protons' similar mass - statistically, neutrons lose half their energy per collision with protons. With hydrogen in the soil, leakage fluxes of neutrons in the 0.5 eV to 500 keV energy range are reduced. A concentration of only1-2 wt% water-equivalent hydrogen results in a decrease in epithermal leakage flux of a factor of two. The leakage flux of thermal neutrons, from 0 to 0.5 eV in energy, can either increase or decrease depending on the hydrogen abundance and stratigraphy. As with the highly successful Lunar Prospector Neutron Spectrometer, the RP NSS detects both thermal and epithermal neutrons by using two helium-3 gas proportional counters, one covered by cadmium and the other uncovered. The former detects only epithermal neutrons with energies above approximately 0.5 eV, the latter detects both thermal (less than 0.5 eV) and epithermal energies (greater than 0.5 eV). When a neutron enters the detector tube and interacts with a helium-3 nucleus, the resulting reaction produces an energetic proton and triton that ionize the gas. The resulting electrons are accelerated toward a high-voltage anode and cascade, amplifying the net charge, which is collected at the anode. The number of electrons produced is proportional to the energy that the triton and proton deposit in the gas. A charge sensitive pre-amplifier converts the total charge to a step voltage output. A shaper amplifier then shapes this step into a uni-polar waveform with peaking time appropriate for the detection depending on the event rate. The integrated shaped waveform, representing the deposited triton/proton energy, is then measured. A histogram, or pulse height analysis, is performed to record the main capture peak and wall effect pulses. A threshold for detection is also required to limit the low amplitude counting rate such as noise floor. The system electronics consists of 2 modules - the Sensor Module (SM) front-end and the Data Processing Module (DPM) back-end circuits. SM is designed as a light-weight and low power front-end housing the two helium-3 proportional counter detectors, preamp and HVPS. It is mounted external to the rover body to detect the thermalized neutron flux with a minimum of host background. The DPM is located inside the rover; it digitizes the SM signals, performs pulse height analysis and accumulates the count rate for each spectral channel. The DPM controls high voltage and thresholding, and sends the science data to the host craft via an RS422 serial asynchronous protocol. The payload host provides all thermal management and control for the SM and DPM.

hydrogen-bearing volatile deposits↗
NASA NTRS2017

What is the Uncertainty in MODIS Aerosol Optical Depth in the Vicinity of Clouds?

MODIS dark-target (DT) algorithm retrieves aerosol optical depth (AOD) using a Look Up Table (LUT) approach. Global comparison of AOD (Collection 6 ) with ground-based sun photometer gives an Estimated Error (EE) of +/-(0.04 + 10%) over ocean. However, EE does not represent per-retrieval uncertainty. For retrievals that are biased high compared to AERONET, here we aim to closely examine the contribution of biases due to presence of clouds and per-pixel retrieval uncertainty. We have characterized AOD uncertainty at 550 nm, due to standard deviation of reflectance in 10 km retrieval region, uncertainty related to gas (H2O, O3) absorption, surface albedo, and aerosol models. The uncertainty in retrieved AOD seems to lie within the estimated over ocean error envelope of +/-(0.03+10%). Regions between broken clouds tend to have higher uncertainty. Compared to C6 AOD, a retrieval omitting observations in the vicinity of clouds (< or = 1 km) is biased by about +/- 0.05. For homogeneous aerosol distribution, clear sky retrievals show near zero bias. Close look at per-pixel reflectance histograms suggests retrieval possibility using median reflectance values.

Patadia, Falguni↗
NASA NTRS2018

Smallholder Crop Area Mapped with Wall-To-Wall WorldView Sub-Meter Panchromatic Image Texture: A Test Case for Tigray, Ethiopia

Global food production in the developing world occurs within sub-hectare fields that are difficult to identify with moderate resolution satellite imagery. Knowledge about the distribution of these fields is critical in food security programs. We developed a semi-automated image segmentation approach using wall-to-wall sub-meter imagery with high-performance computing to map crop area (CA) throughout Tigray, Ethiopia that encompasses over 41,000 km (exp 2). Multiple processing streams were tested to minimize mapping error while applying five unique smoothing kernels to capture differences in land surface texture associated to CA. Typically, very-small fields (mean < 2 ha) have a smooth image roughness compared to natural scrub/shrub woody vegetation at the ~1m scale and these features can be segmented in panchromatic imagery with multi-level histogram thresholding. Multi-temporal very-high resolution (VHR) panchromatic imagery with multi-spectral VHR are sufficient in extracting critical CA information needed in food security programs. A 2011 to 2015 CA map was produced, using over 3000 WorldView-1 panchromatic images wall-to-wall in 1/2 deg mosaics for Tigray, Ethiopia. CA was evaluated with nearly 3000 WorldView-2 2m multispectral 250 X 250 m image subsets by seven expert interpretations, and with in-situ global positioning system photography. CA estimates ranged from 32 to 41% in sub regions of Tigray with median maximum per bin commission and omission errors of 11% and 1% respectively, with most of the error occurring in bins <15%. This empirical, simple, and low direct cost approach via U.S. government license agreement to access commercial VHR data, could be a viable big-data high-performance computing methodology to extract wall-to-wall CA for other regions of the world that have very-small agriculture fields with similar image texture."

Ethiopia↗
NASA NTRS2017

Experimental Visualizations of a Generic Launch Vehicle Flow Field: Time-Resolved Shadowgraph and Infrared Imaging

Time-Resolved shadowgraph and infrared (IR) imaging were performed to investigate off-body and on-body flow features of a generic, 'hammer-head' launch vehicle geometry previously tested by Coe and Nute (1962). The measurements discussed here were one part of a large range of wind tunnel test techniques that included steady-state pressure sensitive paint (PSP), dynamic PSP, unsteady surface pressures, and unsteady force measurements. Image data was captured over a Mach number range of 0.6 less than or equal to M less than or equal to 1.2 at a Reynolds number of 3 million per foot. Both shadowgraph and IR imagery were captured in conjunction with unsteady pressures and forces and correlated with IRIG-B timing. High-speed shadowgraph imagery was used to identify wake structure and reattachment behind the payload fairing of the vehicle. Various data processing strategies were employed and ultimately these results correlated well with the location and magnitude of unsteady surface pressure measurements. Two research grade IR cameras were positioned to image boundary layer transition at the vehicle nose and flow reattachment behind the payload fairing. The poor emissivity of the model surface treatment (fast PSP) proved to be challenging for the infrared measurement. Reference image subtraction and contrast limited adaptive histogram equalization (CLAHE) were used to analyze this dataset. Ultimately turbulent boundary layer transition was observed and located forward of the trip dot line at the model sphere-cone junction. Flow reattachment location was identified behind the payload fairing in both steady and unsteady thermal data. As demonstrated in this effort, recent advances in high-speed and thermal imaging technology have modernized classical techniques providing a new viewpoint for the modern researcher

Shadowgraph↗
NASA NTRS2017

Single Degree-of-Freedom Modeling of SLS Liquid Hydrogen Pre-Valve Flow Guide to Enable Rapid Transient Analysis

A unique single degree-of-freedom approximation technique has been developed to enable rapid application of a temporally-defined multi-spectral semi-narrow-band loading for generation of realistic stress/cycle values compared to a resonant analysis. The technique uses the harmonic analysis at resonance of a high-fidelity finite element model to produce a transfer function, which is then used to calibrate the response of the SDOF model. A standard numerical ordinary differential equation solver is then used to obtain the temporal response, and its histogram is used in a fatigue/fracture model. This technique is related to other SDOF methods used widely in industry, such as Miles' Equation and the Shock Response Spectra, but it is unique in that it produces a realistic time history of the response. The most obvious error in the process, which is the effect of closely-spaced modes, was also assessed using the parallel application of several SDOF models, and the error is shown to be small. The application of this unique and tractable reduced-order methodology has enabled the SLS program to avoid substantial cost and schedule penalties if a redesign or change of material were required. It has also enabled quick analysis of a number of other structures undergoing the same or similar excitation fields, and quick assessment when the excitation and structural configuration has been altered due to design changes in the system.

Brown, Andrew M.↗
NASA NTRS2018

The Kepler Light Curves of AGN: A Detailed Analysis

We present a comprehensive analysis of 21 light curves of Type 1 active galactic nuclei (AGN) from the Kepler spacecraft. First, we describe the necessity and development of a customized pipeline for treating Kepler data of stochastically variable sources like AGN. We then present the light curves, power spectral density functions (PSDs), and flux histograms. The light curves display an astonishing variety of behaviors, many of which would not be detected in ground-based studies, including switching between distinct flux levels. Six objects exhibit PSD flattening at characteristic timescales that roughly correlate with black hole mass. These timescales are consistent with orbital timescales or free-fall accretion timescales. We check for correlations of variability and high-frequency PSD slope with accretion rate, black hole mass, redshift, and luminosity. We find that bolometric luminosity is anticorrelated with both variability and steepness of the PSD slope. We do not find evidence of the linear rms-flux relationships or lognormal flux distributions found in X-ray AGN light curves, indicating that reprocessing is not a significant contributor to optical variability at the 0.1%-10% level.

power spectral density functions(PSDs)↗
NASA NTRS2019

GMAO Seasonal Forecast Ensemble Exploration

GMAO Sub/Seasonal prediction system (S2S) has recently been upgraded. A complete set (1981-2016) of 9-months hindcasts for the previous and current versions (S2S-1.0 and S2S-2.1 respectively) allows for the evaluation of the forecast skill and a study of various characteristics of the ensemble forecasts in particular. We compared the intra-seasonal, interannual and intra-ensemble SST variability of the two systems against the observed. Focusing on the ENSO SST indices, we analyzed the consistency of the forecasts ensembles by studying rank histograms and comparing the ensemble spread with the standard error of the estimate.The S2S-2.1 ensemble appears to be more consistent with observations in Niño1+2 region compared to S2S-1.0, while in the central equatorial Pacific ocean this measure is comparably good for both systems. The S2S-1.0 system tends to be under dispersive, while the new system is under dispersive only at very short lead times, but tends to be over dispersive at long leads and for forecasts verifying in spring in Niño 3.4 region.Overall, the new system has greater skill in predicting ENSO. The evaluation techniques tested here will be applied for testing of the next generation sub/seasonal forecast system under development.

Borovikov, Anna↗
NASA NTRS2019

A Qualitative Investigation of Selected Infrared Flow Visualization Image Processing Techniques

The deployment and integration of high-sensitivity infrared cameras in a transonic wind tunnel testenvironment has resulted in a unique capability to image aerodynamic phenomena in real-time. Multi-camera infrared flow visualization data systems are now routinely utilized at the NASA Ames Unitary Plan Wind Tunnel. The small flow-induced temperature gradients on the surface of the wind tunnel test article coupled with the high bit-depth of the infrared camera sensor makes the processing of the image data critically important. An image processing routine must enhance features of interest with minimal artifacts. Additionally, the production wind tunnel test environment demands that these processed images are made available in a real-time, automatic fashion. Therefore, any image processing routine must be computationally economical and enhance the image data with minimal input from a human operator. The following seeks to qualitatively explore selected image processing techniques by assessing their effectiveness to resolve flow features on a wind tunnel test article. A multi-scale contrast enhancement technique is discussed as well as a new implementation of a multi-scale, non-interpolated adaptive histogram equalization. Finally, a novel method is introduced that demonstrates the ability to resolve flow features imaged on bare-steel test articles possessing low emissivity.This method makes use of dynamic mode decomposition and discrete-time filtering to separate the background reflections that dominate low emissivity surfaces from the aerodynamic driven surface temperature gradients.This process will be shown to resolve the onset of boundary layer transition on a bare metal wing as well as identify and resolve hidden features in the image data. While the implementation of this technique is very preliminary it demonstrates the potential to extend the application of infrared flow-visualization within the wind tunnel test environment.

Image Processeing↗
NASA NTRS2019

Architecture and Dynamics of Kepler's Multi-Transiting Planet Systems.III. Comprehensive Investigation Using All Four Years of Kepler Mission Data

We study the orbital architectures of planetary systems orbiting within ∼1 AU of their stars by analyzing the ensemble of Kepler systems having two or more planet candidates. We use data from the entire Kepler mission, and in many cases we apply improved analysis techniques (e.g., replacing histograms by top-hat Kernel Density Estimators that avoid the loss of information resulting from choosing a particular phase for the bin boundaries) to extend and enhance the studies of Lissauer et al. (2011, ApJS 197, 8) and Fabrycky et al. (2014, ApJ 790, 146).These data show ~ 1700 transiting planet candidates in > 600 multiple-planet systems, far more than were available for our previous two studies. The increased numbers and better information about planetary radii and the properties of stellar hosts made possible by Gaia DR2 allow more statistically-robust analyses of the entire ensemble of Kepler multis as well as independent analyses of subsets of the population. We are thus able to contrast the dynamical configurations of small and large planets, short-period and longer-period planets, and planets orbiting various types of host stars. We reinforce our previous findings that most pairs of planets within the same system are neither in nor near low-order mean motion resonances and that there is a substantial excess of planets having period ratios slightly larger than those of first-order mean-motion resonances. However, neglecting three systems whose planets are locked in 3- body resonances and summing over all first-order mean motion resonances, the deficit of planet pairs with period ratios just narrow of resonance is as large as the excess of planets wide of resonance (within statistical uncertainties), suggesting that overall there is no overall excess of planet pairs in the vicinity of resonance. Other aspects of our study, including estimates of the typical relative inclinations of planetary orbits and their variations as functions of orbital period, planet sizes and stellar properties, are in progress, with results expected to be available for presentation by the time of the conference.

Lissauer, Jack↗
NASA NTRS2019

A New Look at Aubrites: Investigating 3D Modal Mineralogy with X-Ray Computed Tomography

The aubrites (approximately 30 known meteorites) are a unique group of differentiated meteorites that formed on asteroids with oxygen fugacities (ƒO2) from approximately 2 to approximately 6 log units below the iron-wustite buffer. At these highly reduced conditions, elements deviate from the geochemical behavior exhibited at terrestrial ƒO2, forming FeO-poor silicates and exotic sulfides. While previous studies have described the petrology and 2D modal abundances of aubrites, this work investigates the 3D modal mineralogies of silicate, metal, and sulfide phases in aubrite samples, which are then compared to the available 2D data. In addition to 3D modal mineralogies, we have examined the geochemistry of fourteen aubrites, including mineral major-element compositions, bulk-rock compositions, and oxygen isotopic compositions to understand their formation and evolution at extreme ƒO2 conditions. We utilize X-ray computed tomography (XCT) to non-destructively analyze the distribution and abundances of mineral phases in aubrites and locate composite clasts of sulfide grains for future analytical study. In order to better constrain elemental behavior under reduced conditions, we specifically target minerals phases that comprise moderately volatile elements (i.e. oldhamite [CaS], caswellsilverite [NaCrS2] and djerfisherite [K6Na(Fe,Cu,Ni)25S26Cl]) as it has been shown that their geochemical behavior changes as a function of ƒO2. Currently, we have produced 3D scans of the Norton County aubrite. The results of the XCT data have allowed for the determination of the abundances of silicate groundmass (i.e., enstatite, forsterite, albite, and diopside), light (based on electron density) sulfides (i.e. alabandite [MnS] and daubréelite [FeCr2S4]), heavy (based on electron density) sulfides (i.e., troilite [FeS]), and Fe,Ni metal by segmenting a density histogram in Volume Graphics Studio software. XCT scans of additional aubrites are underway. By combining the 3D representation of the exotic phases found in aubrites with existing 2D characterizations, we are able to better determine modal abundances. By integrating 3D and 2D modal abundances and geochemistry, we can ultimately better constrain aubrite petrogenesis and elemental partitioning under reduced conditions. Furthermore, application of this new 3D approach offers the opportunity to identify and select clasts for future study prior to cutting the sample, which will minimize sample loss of this precious material.

Wilbur, Zoe↗
NASA NTRS2019

A Qualitative Review of Selected Infrared Flow Visualization Processing Techniques: Contrast Enhancement and Frequency Domain Analysis

The deployment and integration of high-sensitivity infrared cameras in a transonic wind tunnel test environment has resulted in a unique capability to image aerodynamic phenomena in real-time. Multi-camera infrared flow visualization data systems are now routinely utilized at the NASA Ames Unitary Plan Wind Tunnel. The small flow-induced temperature gradients on the surface of the wind tunnel test article coupled with the high bit-depth of the infrared camera sensor makes the processing of the image data critically important. An image processing routine must enhance features of interest with minimal artifacts. Additionally, the production wind tunnel test environment demands that these processed images are made available in a real-time, automatic fashion. Therefore, any image processing routine must be computationally economical and enhance the image data with minimal input from a human operator. The following seeks to qualitatively explore selected image processing techniques by assessing their effectiveness to resolve flow features on a wind tunnel test article. A multi-scale contrast enhancement technique is introduced as well as a new implementation of a multi-scale, non-interpolated adaptive histogram equalization. Finally, a novel method is introduced that demonstrates the ability to resolve flow features imaged on bare-steel test articles possessing low emissivity. This method merges frequency domain analysis with contrast enhancement and has the potential to extend the application of infrared flow-visualization within the wind tunnel test environment.

Infrared Imaging↗
NASA NTRS2017

Interactive Visualization of High-Dimensional Petascale Ocean Data

We describe an application for interactive visualization of 5 petabytes of time-varying multivariate data from a high-resolution global ocean circulation model. The input data are 10311 hourly (ocean time) time steps of various 2D and 3D fields from a 22-billion point 1/48- degree “lat-lon cap” configuration of the MIT General Circulation Model (MITgcm). We map the global horizontal model domain onto our 128-screen (8x16) tiled display wall to produce a canonical tiling with approximately one MITgcm grid point per display pixel, and using this tiling we encode the entire time series for multiple native and computed scalar quantities at a collection of ocean depths. We reduce disk bandwidth requirements by converting the model’s floating point data to 16-bit fixed point values, and compressing those values with a lossless video encoder, which together allow synchronized playback at 24 time steps per second across all 128 displays. The application allows dynamic assignment of any two encoded tiles to any display, and has multiple interfaces for quickly specifying various orderly arrangements of tiles. All subsequent rendering is done on the fly, with run time control of colormaps, transfer functions, histogram equalization, and labeling. The two data streams on each screen can be rendered independently and combined in various ways, including blending, differencing, horizontal/ vertical wipes, and checkerboarding. The two data streams on any screen can optionally be displayed as a scatterplot in their joint attribute space. All scatterplots and map-view plots from the same x/y location and depth are linked so they all show the current brushable selection. Ocean scientists have used the system, and have found previously unidentified features in the data.

Interactive↗
NASA NTRS2019

Simulation of AMVs for OSSEs

An algorithm to simulate locations of atmospheric motion vectors (AMVs) for use in observing system simulation experiments (OSSEs) is described and demonstrated. It is intended to obviate likely deficiencies in nature run data [standard meteorological parameters and aerosol tracers] if used to produce images for feature tracking. The algorithm employs probabilistic functions that are tuned based on distributions of real observations and histograms of nature run fields. For distinct observation types, the algorithm produces geographical and vertical distributions, time - mean counts, and typical spacings of simulated locations that are, at least, qualitatively similar to those of real observations and are associated with nature run cloud and water vapor fields. It thus appears suitable for generating realistic atmospheric motion vectors for use in observing system simulation experiments.

Errico, R. M.↗