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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

A Comparison of the Multiscale Retinex With Other Image Enhancement Techniques

The multiscale retinex with color restoration (MSRCR) has shown itself to be a very versatile automatic image enhancement algorithm that simultaneously provides dynamic range compression, color constancy, and color rendition. A number of algorithms exist that provide one or more of these features, but not all. In this paper we compare the performance of the MSRCR with techniques that are widely used for image enhancement. Specifically, we compare the MSRCR with color adjustment methods such as gamma correction and gain/offset application, histogram modification techniques such as histogram equalization and manual histogram adjustment, and other more powerful techniques such as homomorphic filtering and 'burning and dodging'. The comparison is carried out by testing the suite of image enhancement methods on a set of diverse images. We find that though some of these techniques work well for some of these images, only the MSRCR performs universally well on the test set.

Rahman, Zia-Ur↗
NASA NTRS2005

Analysis of Multi-Scale Radiometric Data Collected during the Cold Land Processes Experiment-1 (CLPX-1)

Brightness temperature maps at 18.7 and 37 GHz collected at the Fraser and North Park Meso-Scale Areas during the Cold Land Processes Experiment by the NOAA Polarimetric Scanning Radiometer (PSWA) airborne sensor are analyzed. The Fraser site is mostly covered by forest with a typical snowpack depth of 1 m while North Park has no forest cover and is characterized by patches of shallow snow. We examine histograms of the brightness temperatures at 500 m resolution for both the Fraser and North Park areas. The histograms can be modelled by a log-normal distribution in the case of the Fraser MSA and by a bi-modal distribution in the case of the North Park MSA. The histograms of the brightness temperatures at coarser resolutions are also plotted to study the effects of sensor resolution on the shape of the distribution, on the values of the average brightness temperatures and standard deviations. Finally, the values of brightness temperatures obtained by re-sampling (aggregating) the data at 25 km resolution are compared with the values of the brightness temperatures collected by the Advanced Microwave Scanning Radiometer (AMSR-E) and Special Sensor Microwave/Imager (SSMII) satellite radiometers. The results show that in both areas for sensor footprint larger than 5000 m, the brightness temperatures show a flat distribution and the memory of the initial distribution is lost. The values of the brightness temperatures measured by the satellite radiometers are in good agreement with the values obtained averaging the airborne data, even if some discrepancies occur.

Tedesco, M.↗
NASA NTRS2005

Automated Counting of Particles To Quantify Cleanliness

A machine vision system, similar to systems used in microbiological laboratories to count cultured microbes, has been proposed for quantifying the cleanliness of nominally precisely cleaned hardware by counting residual contaminant particles. The system would include a microscope equipped with an electronic camera and circuitry to digitize the camera output, a personal computer programmed with machine-vision and interface software, and digital storage media. A filter pad, through which had been aspirated solvent from rinsing the hardware in question, would be placed on the microscope stage. A high-resolution image of the filter pad would be recorded. The computer would analyze the image and present a histogram of sizes of particles on the filter. On the basis of the histogram and a measure of the desired level of cleanliness, the hardware would be accepted or rejected. If the hardware were accepted, the image would be saved, along with other information, as a quality record. If the hardware were rejected, the histogram and ancillary information would be recorded for analysis of trends. The software would perceive particles that are too large or too numerous to meet a specified particle-distribution profile. Anomalous particles or fibrous material would be flagged for inspection.

Rhode, James↗
NASA NTRS2009

Anvil Clouds of Tropical Mesoscale Convective Systems in Monsoon Regions

The anvil clouds of tropical mesoscale convective systems (MCSs) in West Africa, the Maritime Continent and the Bay of Bengal have been examined with TRMM and CloudSat satellite data and ARM ground-based radar observations. The anvils spreading out from the precipitating cores of MCSs are subdivided into thick, medium and thin portions. The thick portions of anvils show distinct differences from one climatological regime to another. In their upper portions, the thick anvils of West Africa MCSs have a broad, flat histogram of reflectivity, and a maximum of reflectivity in their lower portions. The reflectivity histogram of the Bay of Bengal thick anvils has a sharply peaked distribution of reflectivity at all altitudes with modal values that increase monotonically downward. The reflectivity histogram of the Maritime Continent thick anvils is intermediate between that of the West Africa and Bay of Bengal anvils, consistent with the fact this region comprises a mix of land and ocean influences. It is suggested that the difference between the statistics of the continental and oceanic anvils is related to some combination of two factors: (1) the West African anvils tend to be closely tied to the convective regions of MCSs while the oceanic anvils are more likely to be extending outward from large stratiform precipitation areas of MCSs, and (2) the West African MCSs result from greater buoyancy, so that the convective cells are more likely to produce graupel particles and detrain them into anvils

Cetrone, J.↗
NASA NTRS2015

Damage Proxy Map from Interferometric Synthetic Aperture Radar Coherence

A method, apparatus, and article of manufacture provide the ability to generate a damage proxy map. A master coherence map and a slave coherence map, for an area prior and subsequent to (including) a damage event are obtained. The slave coherence map is registered to the master coherence map. Pixel values of the slave coherence map are modified using histogram matching to provide a first histogram of the master coherence map that exactly matches a second histogram of the slave coherence map. A coherence difference between the slave coherence map and the master coherence map is computed to produce a damage proxy map. The damage proxy map is displayed with the coherence difference displayed in a visually distinguishable manner.

Yun, Sang-Ho↗
NASA NTRS2018

Magnetospheric Multiscale Observations of Turbulent Magnetic and Electron Velocity Fluctuations in Earth's Magnetosheath Downstream of a Quasi-Parallel Bow Shock

We present statistical single-spacecraft observations of magnetic and electron velocity fluctuations in Earth's magnetosheath, likely in the vicinity of the magnetopause, downstream of a bow shock immersed in quasi-parallel interplanetary magnetic field conditions, a situation conducive to plasma turbulence in the downstream flow. These fluctuations exhibit scale-dependent behavior, wherein histograms of their Partial Variance of Increments (PVIB or PVIV(sub e)) demonstrate highly non-Gaussian forms at small scales and are reasonably well-described by kappa distributions, albeit with fitted values of the kappa parameter only slightly larger than 1.5, exemplifying their power law nature at large values of PVI. At larger scales, the PVI histograms lose their non-Gaussian nature and are well described by both Gaussian and kappa distributions with large values of the kappa parameter. The PVI histograms furthermore exhibit kurtosis that increases with decreasing scale, a characteristic that is much more prominent in the magnetic fluctuations than in the electron velocity fluctuations. This feature that is not yet explained. In both cases, the results are characteristic of turbulent intermittency.

Pollock, C. J.↗
NASA NTRS2021

Probabilistic Processing of Asynchronous Geiger-Mode Avalanche Photodiode Arrays for Background Rate Measurements

A system-level performance evaluation of Geiger-mode avalanche photodiode (GmAPD) arrays requires accurate measurement and prediction of the background rate of the device due to dark counts and other spurious detection events. Since a GmAPD detector reports only a binary value and timestamp associated with an avalanche event, dark count rates are typically measured by averaging thousands of frames to support a statistically significant measurement. For both synchronous and asynchronous detector, the Poisson distributed background rates are referenced to the time each pixel is armed. Unlike for synchronous GmAPD imagers where all the pixels are armed to an array-wide arm signal, an asynchronous pixel operates independently from its neighboring pixels; requiring the background rates to be calculated using an interarrival histogram. For both types of imagers, the background rate is typically evaluated by fitting an exponential distribution to a fixed window within a measured histogram of time intervals between detection events However, if the statistics of the background rate are insufficient – whether that is due to low population sizes, saturation, or a large dynamic range of population size across the array, the pixel, or array-wide, performance metrics may report results with varying accuracy. This paper reports on an implementation of an algorithm that evaluates GmAPD background rates based on statistical metrics rather than fixed windows. The algorithm functions by determining the appropriate integration window within the interarrival time histogram based on a per-pixel count rate set by a predetermined tolerable measurement error. The implementation of the algorithm allows us to characterize GmAPD arrays with orders of magnitude spread in background rates across the detector using common statistical parameters.

Mumolo, J.↗
DOE OSTI2026

Reconstruction of Six-Dimensional Phase Space

A phase space is a mathematical representation of all possible physical states of a system. Particle beams at Fermilab exist within a six-dimensional (6D) phase space defined by three positional components, (x, y, z) and three momentum components, (px, py, pz). To reconstruct this space implies taking measurement data from detectors and mapping out particle behavior using computational methods. The beam detectors, however, are only able to detect spatial distribution among the events of the beam, therefore being limited to positional data. Also, due to the vast number of events in a particle beam, it is extremely difficult to analyze and differentiate every single one’s behavior. However, with Machine Learning (ML), which can distinguish between patterns and map out particle behavior more efficiently. We first used the particle beam software, G4beamline, to simulate a 10,000-event muon beam, adjusting parameters such as initial momentum magnitude (p¬0) and virtual detector position. Using ten virtual detectors, we analyzed p0 values such that minimum 9,990 events were analyzed by every detector. We then input the data from these beam simulations to a C++ program, that randomly selects 100 events, and creates a 2D histogram based on spatial distribution, detector position, and event intensity. This process is repeated 100 times to create 100 histograms per p0 value. These images were then input to a modified ResNet18 Convolutional Neural Network (CNN) for training, and to predict p0 from some unseen set of histograms. The model was accurate when trained on momentum increments of 5 MeV/c and provided with denser training samples around highly variable test values. These results displayed machine learning being able to accurately predict p0 from being trained on different particle behaviors.

Shirlee, Jermain [Fermilab]↗
DOE OSTI2023

A direct numerical simulation study for confined non-isothermal jet impingement at moderate nozzle-to-plate distances: Capturing jet-to-ambient density effects

A direct numerical simulation (DNS) campaign is deployed for a series of confined downward oriented, non-isothermal turbulent impinging jet configurations. A baseline Reynolds number of 9960 is obtained through a precursor DNS pipe flow simulation (Re τ = 505). Three jet temperature configurations (confinement height to nozzle diameter of three) enter a cylindrical domain that share ambient and impingement plate temperatures (298.15 K). The range of jet temperatures are crafted such that the ratio of inlet to ambient density varies from unity to 0.52, showcasing the effect of density disparity on flow characteristics such as core collapse, radial mixing of momentum and energy, near-wall stagnation behavior, wall-jet profiles, and large-scale vortical structures. Surface quantities provided include mean radial heat flux and wall-shear stress profiles, and heat flux histograms at select radial stations. Results showcase increased radial normal stresses for higher temperature jets that support increased mixing, resulting in large-scale recirculation structures that are smaller, while retaining similar normalized radial wall profiles for shear stress, heat flux and pressure. Furthermore, radial plots for wall shear stress and Nusselt number showcase strong radial decay as compared to previous configurations that share similar jet and ambient temperatures. For the 373.15 K case, a Gaussian-like histogram for heat fluxes at the impingement plate transitions to a log-normal profile as radial distances increase. In contrast, the 573.15 K configuration displays a bi-modal heat flux characteristic at the impingement plate, and in similar manner to the moderate temperature counterpart, transitions to a log-normal profile at larger radial distances.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

Validation of the two-region Rossi-alpha model for reflected assemblies

Nuclear criticality safety, nonproliferation and safeguards, emergency response, and stockpile stewardship utilize estimates of the k eff multiplication factor. The value of k eff cannot be directly measured, but it can be inferred from the prompt neutron period. One modality of measuring the prompt neutron period is the Rossi-alpha method, which is predicated on fitting a histogram of detection time differences due to the nonrandom temporal distribution of same-fission-chain neutrons. Recent works have developed the motivation and theory to expand traditional one-region point kinetic models resulting in one-exponential histogram fits to two-region models that result in two-exponential fits. This paper validates the new two-region model using organic scintillator measurements of copper-reflected weapons-grade plutonium (0:83 ≤ k eff ≤ 0:94) and high-density- polyethylene-reflected highly enriched uranium (0:73 ≤ k eff ≤ 0:95). Furthermore, the results show that more thermal systems have shorter prompt neutron periods in the core region due to increased induced-fission probabilities for moderated neutrons. Here, a new parameter introduced by the two- region model is also shown to be correlated to the amount of reflection, and may be used to infer assembly properties such as type and amount of reflector.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

On the instrument-dependent appearance of ion dissociation events in atom probe tomography mass spectra

The successful application of atom probe tomography (APT) relies on the accurate interpretation of the mass spectrum (i.e. m/z histogram) from a sample. Some materials yield mass spectra that are amenable to a straightforward peak assignment/ranging, however, there are many materials that produce mass spectra with features that defy simple interpretation. One such example is Ga 2 O 3 which yields mass spectra containing several broad and difficult to interpret features. Herein, we study the GaO 2+ → O 1+ + Ga 1+ dissociation and we explain how this dissociation process gives rise to broad and previously unassigned features in the mass spectrum. Trajectory simulations are performed for the dissociation reaction utilizing realistic electrostatic models and compared to experiments using commercially available straight flight and reflectron based local electrode (LE) APT instruments. It is shown that the appearance of these features is strongly dependent on the specific design of the time-of-flight (ToF) mass analyzer. Additionally, we explore how various experimental parameters can affect the appearance of the dissociation process in the one-dimensional (1D) mass spectrum and in the two-dimensional (2D) correlation histogram. While the focus of this work is on a particular dissociation process related to Ga 2 O 3 , the understanding gained in the course of these simulations and experiments should be applicable to the interpretation of dissociation processes in other materials.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Fast correlation function calculator: A high-performance pair-counting toolkit

A novel high-performance exact pair-counting toolkit called fast correlation function calculator (FCFC) is presented. With the rapid growth of modern cosmological datasets, the evaluation of correlation functions with observational and simulation catalogues has become a challenge. High-efficiency pair-counting codes are thus in great demand. We introduce different data structures and algorithms that can be used for pair-counting problems, and perform comprehensive benchmarks to identify the most efficient algorithms for real-world cosmological applications. We then describe the three levels of parallelisms used by FCFC, SIMD, OpenMP, and MPI, and run extensive tests to investigate the scalabilities. Finally, we compare the efficiency of FCFC with alternative pair-counting codes. The data structures and histogram update algorithms implemented in FCFC are shown to outperform alternative methods. FCFC does not benefit greatly from SIMD because the bottleneck of our histogram update algorithm is mainly cache latency. Nevertheless, the efficiency of FCFC scales well with the numbers of OpenMP threads and MPI processes, even though speedups may be degraded with over a few thousand threads in total. FCFC is found to be faster than most (if not all) other public pair-counting codes for modern cosmological pair-counting applications.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

An implementation of neural simulation-based inference for parameter estimation in ATLAS

Neural simulation-based inference (NSBI) is a powerful class of machine-learning-based methods for statistical inference that naturally handles high-dimensional parameter estimation without the need to bin data into low-dimensional summary histograms. Such methods are promising for a range of measurements, including at the Large Hadron Collider, where no single observable may be optimal to scan over the entire theoretical phase space under consideration, or where binning data into histograms could result in a loss of sensitivity. This work develops a NSBI framework for statistical inference, using neural networks to estimate probability density ratios, which enables the application to a full-scale analysis. It incorporates a large number of systematic uncertainties, quantifies the uncertainty due to the finite number of events in training samples, develops a method to construct confidence intervals, and demonstrates a series of intermediate diagnostic checks that can be performed to validate the robustness of the method. As an example, the power and feasibility of the method are assessed on simulated data for a simplified version of an off-shell Higgs boson couplings measurement in the four-lepton final states. This approach represents an extension to the standard statistical methodology used by the experiments at the Large Hadron Collider, and can benefit many physics analyses.

frequentist statistics↗
DOE OSTI2025

Discriminative versus generative approaches to simulation-based inference

Most of the fundamental, emergent, and phenomenological parameters of particle and nuclear physics are determined through parametric template fits. Simulations are used to populate histograms which are then matched to data. This approach is inherently lossy, since histograms are binned and low-dimensional. Deep learning has enabled unbinned and high-dimensional parameter estimation through neural likelihood(-ratio) estimation. We compare two approaches for neural simulation-based inference (NSBI): one based on discriminative learning (classification) and one based on generative modeling. These two approaches are directly evaluated on the same datasets, with a similar level of hyperparameter optimization in both cases. In addition to a Gaussian dataset, we study NSBI using a Higgs boson dataset from the FAIR Universe Challenge. We find that both the direct likelihood and likelihood ratio estimation are able to effectively extract parameters with reasonable uncertainties. For the numerical examples and within the set of hyperparameters studied, we found that the likelihood ratio method is more accurate and/or precise. Both methods have a significant spread from the network training and would require ensembling or other mitigation strategies in practice.

high energy physics↗
DOE OSTI2023

High-energy synchrotron flares powered by strongly radiative relativistic magnetic reconnection: 2D and 3D PIC simulations

ABSTRACT The time evolution of high-energy synchrotron radiation generated in a relativistic pair plasma energized by reconnection of strong magnetic fields is investigated with 2D and 3D particle-in-cell (PIC) simulations. The simulations in this 2D/3D comparison study are conducted with the radiative PIC code OSIRIS, which self-consistently accounts for the synchrotron radiation reaction on the emitting particles, and enables us to explore the effects of synchrotron cooling. Magnetic reconnection causes compression of the plasma and magnetic field deep inside magnetic islands (plasmoids), leading to an enhancement of the flaring emission, which may help explain some astrophysical gamma-ray flare observations. Although radiative cooling weakens the emission from plasmoid cores, it facilitates additional compression there, further amplifying the magnetic field B and plasma density n, and thus partially mitigating this effect. Novel simulation diagnostics utilizing 2D histograms in the n-B space are developed and used to visualize and quantify the effects of compression. The n-B histograms are observed to be bounded by relatively sharp power-law boundaries marking clear limits on compression. Theoretical explanations for some of these compression limits are developed, rooted in radiative resistivity or 3D kinking instabilities. Systematic parameter-space studies with respect to guide magnetic field, system size, and upstream magnetization are conducted and suggest that stronger compression, brighter high-energy radiation, and perhaps significant quantum electrodynamic effects such as pair production, may occur in environments with larger reconnection-region sizes and higher magnetization, particularly when magnetic field strengths approach the critical (Schwinger) field, as found in magnetar magnetospheres.

Schoeffler, K. M. (ORCID:0000000276447381)↗
DOE OSTI2021

Geometrical-Based Generative Adversarial Network to Enhance Digital Rock Image Quality

X-ray microcomputed tomography (micro-CT) is a common tool for the study of porous media structures and properties. High-quality micro-CT data are required to accurately capture pore structures. Acquiring high-quality micro-CT data, however, is not always possible, owing to application limitations and experimental constraints. Therefore, we propose a geometrical-based generative adversarial network (GAN) to rapidly restore noisy micro-CT images to their clean counterparts. The training data and related ground-truth (GT) data are scanned for 7 min and 9.5 h, respectively. To evaluate the performance of the geometrical-based GAN, a 6003 voxel image that has never been used for training is reconstructed and compared with the corresponding GT image. Histogram matching and linear normalization are implemented to adjust the histogram of the reconstructed image to that of the GT image. A watershed-based segmentation method is then applied to delineate pore and solid phases. Lastly, we measure the Minkowski functionals and petrophysical properties, including absolute permeability, pore size distribution, drainage capillary pressure-saturation curve, and imbibition relative permeability, to estimate the physical accuracy of the denoised image. The results show that the proposed geometrical-based GAN can accurately restore noisy micro-CT data. By reducing the scanning time from 9.5 h to 7 min, the expenditure of collecting micro-CT can be decreased significantly. This is particularly important for applications where time-lapse images of a dynamic process are required, high-throughput imaging is necessary for real-time data analysis, or where the quantification of large sample volumes is required.

58 GEOSCIENCES↗
DOE OSTI2023

MILK : a Python scripting interface to MAUD for automation of Rietveld analysis

Modern diffraction experiments ( e.g. in situ parametric studies) present scientists with many diffraction patterns to analyze. Interactive analyses via graphical user interfaces tend to slow down obtaining quantitative results such as lattice parameters and phase fractions. Furthermore, Rietveld refinement strategies ( i.e. the parameter turn-on-off sequences) tend to be instrument specific or even specific to a given dataset, such that selection of strategies can become a bottleneck for efficient data analysis. Managing multi-histogram datasets such as from multi-bank neutron diffractometers or caked 2D synchrotron data presents additional challenges due to the large number of histogram-specific parameters. To overcome these challenges in the Rietveld software Material Analysis Using Diffraction ( MAUD ), the MAUD Interface Language Kit ( MILK ) is developed along with an updated text batch interface for MAUD . The open-source software MILK is computer-platform independent and is packaged as a Python library that interfaces with MAUD . Using MILK , model selection ( e.g. various texture or peak-broadening models), Rietveld parameter manipulation and distributed parallel batch computing can be performed through a high-level Python interface. A high-level interface enables analysis workflows to be easily programmed, shared and applied to large datasets, and external tools to be integrated with MAUD . Through modification to the MAUD batch interface, plot and data exports have been improved. The resulting hierarchical folders from Rietveld refinements with MILK are compatible with Cinema: Debye–Scherrer , a tool for visualizing and inspecting the results of multi-parameter analyses of large quantities of diffraction data. In this manuscript, the combined Python scripting and visualization capability of MILK is demonstrated with a quantitative texture and phase analysis of data collected at the HIPPO neutron diffractometer.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Bimodal Visualization of Industrial X-Ray and Neutron Computed Tomography Data

Advanced manufacturing creates increasingly complex objects with material compositions that are often difficult to characterize by a single modality. Our collaborating domain scientists are going beyond traditional methods by employing both X-ray and neutron computed tomography to obtain complementary representations expected to better resolve material boundaries. However, the use of two modalities creates its own challenges for visualization, requiring either complex adjustments of bimodal transfer functions or the need for multiple views. Together with experts in nondestructive evaluation, we designed a novel interactive bimodal visualization approach to create a combined view of the co-registered X-ray and neutron acquisitions of industrial objects. Using an automatic topological segmentation of the bivariate histogram of X-ray and neutron values as a starting point, the system provides a simple yet effective interface to easily create, explore, and adjust a bimodal visualization. Here, we propose a widget with simple brushing interactions that enables the user to quickly correct the segmented histogram results. Our semiautomated system enables domain experts to intuitively explore large bimodal datasets without the need for either advanced segmentation algorithms or knowledge of visualization techniques. We demonstrate our approach using synthetic examples, industrial phantom objects created to stress bimodal scanning techniques, and real-world objects, and we discuss expert feedback.

image segmentation↗