Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “histogram”

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 199 records · Page 11

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↗

Conductance fluctuation of edge-disordered graphene nanoribbons: Crossover from diffusive transport to Anderson localization

Conductance fluctuation of edge-disordered graphene nanoribbons (ED-GNRs) is examined using the non-equilibrium Green's function technique combined with the extended Hückel approximation. The mean free path λ and the localization length ξ of the ED-GNRs are determined to classify the quantum transport regimes. In the diffusive regime where the length L{sub c} of the ED-GNRs is much longer than λ and much shorter than ξ, the conductance histogram is given by a Gaussian distribution function with universal conductance fluctuation. In the localization regime where L{sub c}≫ξ, the histogram is no longer the universal Gaussian distribution but a lognormal distribution that characterizes Anderson localization.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Two-state theory of binned photon statistics for a large class of waiting time distributions and its application to quantum dot blinking

A theoretical method is proposed for the calculation of the photon counting probability distribution during a bin time. Two-state fluorescence and steady excitation are assumed. A key feature is a kinetic scheme that allows for an extensive class of stochastic waiting time distribution functions, including power laws, expanded as a sum of weighted decaying exponentials. The solution is analytic in certain conditions, and an exact and simple expression is found for the integral contribution of “bright” and “dark” states. As an application for power law kinetics, theoretical results are compared with experimental intensity histograms from a number of blinking CdSe/ZnS quantum dots. The histograms are consistent with distributions of intensity states around a “bright” and a “dark” maximum. A gap of states is also revealed in the more-or-less flat inter-peak region. The slope and to some extent the flatness of the inter-peak feature are found to be sensitive to the power-law exponents. Possible models consistent with these findings are discussed, such as the combination of multiple charging and fluctuating non-radiative channels or the multiple recombination center model. A fitting of the latter to experiment provides constraints on the interaction parameter between the recombination centers. Further extensions and applications of the photon counting theory are also discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamics of classical particles in oval or elliptic billiards with a dispersing mechanism

Some dynamical properties for an oval billiard with a scatterer in its interior are studied. The dynamics consists of a classical particle colliding between an inner circle and an external boundary given by an oval, elliptical, or circle shapes, exploring for the first time some natural generalizations. The billiard is indeed a generalization of the annular billiard, which is of strong interest for understanding marginally unstable periodic orbits and their role in the boundary between regular and chaotic regions in both classical and quantum (including experimental) systems. For the oval billiard, which has a mixed phase space, the presence of an obstacle is an interesting addition. We demonstrate, with details, how to obtain the equations of the mapping, and the changes in the phase space are discussed. We study the linear stability of some fixed points and show both analytically and numerically the occurrence of direct and inverse parabolic bifurcations. Lyapunov exponents and generalized bifurcation diagrams are obtained. Moreover, histograms of the number of successive iterations for orbits that stay in a cusp are studied. These histograms are shown to be scaling invariant when changing the radius of the scatterer, and they have a power law slope around −3. The results here can be generalized to other kinds of external boundaries.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Large-deviation joint statistics of the finite-time Lyapunov spectrum in isotropic turbulence

One of the hallmarks of turbulent flows is the chaotic behavior of fluid particle paths with exponentially growing separation among them while their distance does not exceed the viscous range. The maximal (positive) Lyapunov exponent represents the average strength of the exponential growth rate, while fluctuations in the rate of growth are characterized by the finite-time Lyapunov exponents (FTLEs). In the last decade or so, the notion of Lagrangian coherent structures (which are often computed using FTLEs) has gained attention as a tool for visualizing coherent trajectory patterns in a flow and distinguishing regions of the flow with different mixing properties. A quantitative statistical characterization of FTLEs can be accomplished using the statistical theory of large deviations, based on the so-called Cramér function. To obtain the Cramér function from data, we use both the method based on measuring moments and measuring histograms and introduce a finite-size correction to the histogram-based method. We generalize the existing univariate formalism to the joint distributions of the two FTLEs needed to fully specify the Lyapunov spectrum in 3D flows. The joint Cramér function of turbulence is measured from two direct numerical simulation datasets of isotropic turbulence. Results are compared with joint statistics of FTLEs computed using only the symmetric part of the velocity gradient tensor, as well as with joint statistics of instantaneous strain-rate eigenvalues. When using only the strain contribution of the velocity gradient, the maximal FTLE nearly doubles in magnitude, highlighting the role of rotation in de-correlating the fluid deformations along particle paths. We also extend the large-deviation theory to study the statistics of the ratio of FTLEs. The most likely ratio of the FTLEs λ{sub 1} : λ{sub 2} : λ{sub 3} is shown to be about 4:1:−5, compared to about 8:3:−11 when using only the strain-rate tensor for calculating fluid volume deformations. The results serve to characterize the fundamental statistical and geometric structure of turbulence at small scales including cumulative, time integrated effects. These are important for deformable particles such as droplets and polymers advected by turbulence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Path integral Monte Carlo simulations of H{sub 2} adsorbed to lithium-doped benzene: A model for hydrogen storage materials

Finite temperature quantum and anharmonic effects are studied in H{sub 2}–Li{sup +}-benzene, a model hydrogen storage material, using path integral Monte Carlo (PIMC) simulations on an interpolated potential energy surface refined over the eight intermolecular degrees of freedom based upon M05-2X/6-311+G(2df,p) density functional theory calculations. Rigid-body PIMC simulations are performed at temperatures ranging from 77 K to 150 K, producing both quantum and classical probability density histograms describing the adsorbed H{sub 2}. Quantum effects broaden the histograms with respect to their classical analogues and increase the expectation values of the radial and angular polar coordinates describing the location of the center-of-mass of the H{sub 2} molecule. The rigid-body PIMC simulations also provide estimates of the change in internal energy, ΔU{sub ads}, and enthalpy, ΔH{sub ads}, for H{sub 2} adsorption onto Li{sup +}-benzene, as a function of temperature. These estimates indicate that quantum effects are important even at room temperature and classical results should be interpreted with caution. Our results also show that anharmonicity is more important in the calculation of U and H than coupling—coupling between the intermolecular degrees of freedom becomes less important as temperature increases whereas anharmonicity becomes more important. The most anharmonic motions in H{sub 2}–Li{sup +}-benzene are the “helicopter” and “ferris wheel” H{sub 2} rotations. Treating these motions as one-dimensional free and hindered rotors, respectively, provides simple corrections to standard harmonic oscillator, rigid rotor thermochemical expressions for internal energy and enthalpy that encapsulate the majority of the anharmonicity. At 150 K, our best rigid-body PIMC estimates for ΔU{sub ads} and ΔH{sub ads} are −13.3 ± 0.1 and −14.5 ± 0.1 kJ mol{sup −1}, respectively.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Improved neutron-gamma discrimination for a {sup 6}Li-glass neutron detector using digital signal analysis methods

A {sup 6}Li-glass scintillator (GS20) based neutron Anger camera was developed for time-of-flight single-crystal diffraction instruments at Spallation Neutron Source. Traditional Pulse-Height Analysis (PHA) for Neutron-Gamma Discrimination (NGD) resulted in the neutron-gamma efficiency ratio (defined as NGD ratio) on the order of 10{sup 4}. The NGD ratios of Anger cameras need to be improved for broader applications including neutron reflectometers. For this purpose, six digital signal analysis methods of individual waveforms acquired from photomultiplier tubes were proposed using (i) charge integration, (ii) pulse-amplitude histograms, (iii) power spectrum analysis combined with the maximum pulse-amplitude, (iv) two event parameters (a{sub 1}, b{sub 0}) obtained from a Wiener filter, (v) an effective amplitude (m) obtained from an adaptive least-mean-square filter, and (vi) a cross-correlation coefficient between individual and reference waveforms. The NGD ratios are about 70 times those from the traditional PHA method. Our results indicate the NGD capabilities of neutron Anger cameras based on GS20 scintillators can be significantly improved with digital signal analysis methods.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Resolving amorphous solid-liquid interfaces by atomic force microscopy

Recent advancements in liquid atomic force microscopy make it an ideal technique for probing the structure of solid-liquid interfaces. Here, we present a structural study of a two-dimensional amorphous silica bilayer immersed in an aqueous solution utilizing liquid atomic force microscopy with sub-nanometer resolution. Structures show good agreement with atomically resolved ultra-high vacuum scanning tunneling microscopy images obtained on the same sample system, owing to the structural stability of the silica bilayer and the imaging clarity from the two-dimensional sample system. Pair distance histograms of ring center positions are utilized to develop quantitative metrics for structural comparison, and the physical origin of pair distance histogram peaks is addressed by direct assessment of real space structures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Free energy landscape and molecular pathways of gas hydrate nucleation

Despite the significance of gas hydrates in diverse areas, a quantitative knowledge of hydrate formation at a molecular level is missing. The impediment to acquiring this understanding is primarily attributed to the stochastic nature and ultra-fine scales of nucleation events, posing a great challenge for both experiment and simulation to explore hydrate nucleation. Here we employ advanced molecular simulation methods, including forward flux sampling (FFS), p{sub B} histogram analysis, and backward flux sampling, to overcome the limit of direct molecular simulation for exploring both the free energy landscape and molecular pathways of hydrate nucleation. First we test the half-cage order parameter (H-COP) which we developed for driving FFS, through conducting the p{sub B} histogram analysis. Our results indeed show that H-COP describes well the reaction coordinates of hydrate nucleation. Through the verified order parameter, we then directly compute the free energy landscape for hydrate nucleation by combining both forward and backward flux sampling. The calculated stationary distribution density, which is obtained independently of nucleation theory, is found to fit well against the classical nucleation theory (CNT). Subsequent analysis of the obtained large ensemble of hydrate nucleation trajectories show that although on average, hydrate formation is facilitated by a two-step like mechanism involving a gradual transition from an amorphous to a crystalline structure, there also exist nucleation pathways where hydrate crystallizes directly, without going through the amorphous stage. The CNT-like free energy profile and the structural diversity suggest the existence of multiple active transition pathways for hydrate nucleation, and possibly also imply the near degeneracy in their free energy profiles among different pathways. Our results thus bring a new perspective to the long standing question of how hydrates crystallize.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Time filtering of event based neutron scattering data: A pathway to study the dynamic structural responses of materials

Time-resolved diffraction has become a vital tool for probing dynamic responses to an applied stimulus. Such experiments traditionally use hardware solutions to histogram measured data into their respective bin. We will show that a major advantage of event-based data acquisition, which time-stamps measured diffraction data with 100 ns accuracy, is much preferred over hardware histogramming of the data by enabling postprocessing for advanced custom binning using a software solution. This approach is made even more powerful by coupling measured diffraction data with metadata about the applied stimuli and material response. In this work, we present a time-filter approach that leverages the power of event-based diffraction collection to reduce stroboscopic data measured over many hours into equally weighted segments that represent subsets of the response to a single cycle of the applied stimulus. We demonstrate this approach by observing ferroelectric/ferroelastic domain wall motion during electric field cycling of BaTiO 3 . Furthermore, the developed approach can readily be expanded to investigate other dynamic phenomena using complex sample environments.

47 OTHER INSTRUMENTATION↗

NON-THERMAL RESPONSE OF THE CORONA TO THE MAGNETIC FLUX DISPERSAL IN THE PHOTOSPHERE OF A DECAYING ACTIVE REGION

We analyzed Solar Dynamics Observatory line-of-sight magnetograms for a decaying NOAA active region (AR) 11451 along with co-temporal Extreme-Ultraviolet Imaging Spectrometer (EIS) data from the Hinode spacecraft. The photosphere was studied via time variations of the turbulent magnetic diffusivity coefficient, {eta}(t), and the magnetic power spectrum index, {alpha}, through analysis of magnetogram data from the Helioseismic and Magnetic Imager (HMI). These measure the intensity of the random motions of magnetic elements and the state of turbulence of the magnetic field, respectively. The time changes of the non-thermal energy release in the corona was explored via histogram analysis of the non-thermal velocity, v {sub nt}, in order to highlight the largest values at each time, which may indicate an increase in energy release in the corona. We used the 10% upper range of the histogram of v {sub nt} (which we called V {sup upp} {sub nt}) of the coronal spectral line of Fe XII 195 A. A 2 day time interval was analyzed from HMI data, along with the EIS data for the same field of view. Our main findings are the following. (1) The magnetic turbulent diffusion coefficient, {eta}(t), precedes the upper range of the v {sub nt} with the time lag of approximately 2 hr and the cross-correlation coefficient of 0.76. (2) The power-law index, {alpha}, of the magnetic power spectrum precedes V {sup upp} {sub nt} with a time lag of approximately 3 hr and the cross-correlation coefficient of 0.5. The data show that the magnetic flux dispersal in the photosphere is relevant to non-thermal energy release dynamics in the above corona. The results are consistent with the nanoflare mechanism of the coronal heating, due to the time lags being consistent with the process of heating and cooling the loops heated by nanoflares.

79 ASTRONOMY AND ASTROPHYSICS↗

Fluctuation statistics in the scrape-off layer of Alcator C-Mod

We study long time series of the ion saturation current and floating potential, sampled by Langmuir probes dwelled in the outboard mid-plane scrape-off layer and embedded in the lower divertor baffle of Alcator C-Mod. A series of ohmically heated L-mode plasma discharges is investigated with line-averaged plasma density ranging from ${{\bar{n}}_{\text{e}}}/{{n}_{\text{G}}}=0.15$ to 0.42, where $n_ G$ is the Greenwald density. All ion saturation current time series that are sampled in the far scrape-off layer are characterized by large-amplitude burst events. Coefficients of skewness and excess kurtosis of the time series obey a quadratic relationship and their histograms coincide partially upon proper normalization. Histograms of the ion saturation current time series are found to agree well with a prediction of a stochastic model for the particle density fluctuations in scrape-off layer plasmas. The distribution of the waiting times between successive large-amplitude burst events and of the burst amplitudes are approximately described by exponential distributions. The average waiting time and burst amplitude are found to vary weakly with the line-averaged plasma density. Conditional averaging reveals that the radial blob velocity, estimated from floating potential measurements, increases with the normalized burst amplitude in the outboard mid-plane scrape-off layer. For low density discharges, the conditionally averaged waveform of the floating potential associated with large amplitude bursts at the divertor probes has a dipolar shape. In detached divertor conditions the average waveform is random, indicating electrical disconnection of blobs from the sheaths at the divertor targets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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↗

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↗

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)↗

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↗

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↗

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↗