Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “clustering statistics”

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 127 records · Page 7

Producing two-dimensional dust clouds and clusters using a movable electrode for complex plasma and fundamental physics experiments

We report a Bidirectional Electrode Control Arm Assembly (BECAA) for precisely manipulating dust clouds levitated above the powered electrode in RF plasmas. The reported techniques allow the creation of perfectly 2D dust layers by eliminating off-plane particles by moving the electrode from outside the plasma chamber without altering the plasma conditions. Here, the tilting and moving of electrodes using BECAA also allows the precise and repeatable elimination of dust particles one by one to achieve any desired number of grains N without trial and error. Simultaneously acquired top and side view images of dust clusters show that they are perfectly planar or 2D. A demonstration of clusters with N = 1–28 without changing the plasma conditions is presented to show the utility of BECAA for complex plasma and statistical physics experimental design. Demonstration videos and 3D printable part files are available for easy reproduction and adaptation of this new method to repeatably produce 2D clusters in existing RF plasma chambers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Machine-learning-assisted analysis of transition metal dichalcogenide thin-film growth

In situ reflective high-energy electron diffraction (RHEED) is widely used to monitor the surface crystalline state during thin-film growth by molecular beam epitaxy (MBE) and pulsed laser deposition. With the recent development of machine learning (ML), ML-assisted analysis of RHEED videos aids in interpreting the complete RHEED data of oxide thin films. The quantitative analysis of RHEED data allows us to characterize and categorize the growth modes step by step, and extract hidden knowledge of the epitaxial film growth process. In this study, we employed the ML-assisted RHEED analysis method to investigate the growth of 2D thin films of transition metal dichalcogenides (ReSe2) on graphene substrates by MBE. Principal component analysis (PCA) and K-means clustering were used to separate statistically important patterns and visualize the trend of pattern evolution without any notable loss of information. Using the modified PCA, we could monitor the diffraction intensity of solely the ReSe2 layers by filtering out the substrate contribution. These findings demonstrate that ML analysis can be successfully employed to examine and understand the film-growth dynamics of 2D materials. Further, the ML-based method can pave the way for the development of advanced real-time monitoring and autonomous material synthesis techniques.

36 MATERIALS SCIENCE↗

Many-body Nuclear Dynamics

The principal goal of this work was to better understand the fusion of neutron-rich nuclei a topic relevant to the fields of both nuclear physics and nuclear astrophysics. The work provides insight into the structure and reactions of neutron-rich nuclei namely the extent of their neutron density distribution and its polarizability. By comparing the fusion excitation functions for a chain of isotopes with a common target nucleus changes in the attractive nuclear potential are assessed. This change in the attractive potential is related to changes in the neutron density distribution with increasing number of neutrons or changes in fusion dynamics with increasing neutron number. The impact of the pairing of valence neutrons and protons on the fusion cross-section is also examined. Measurement of an isotopic chain is a powerful tool to address this topic. The experimental program made use of several different accelerator facilities. The core of the experimental work involved experiments at the ReA3 accelerator and the Facility for Rare Isotope Beams (FRIB) situated at Michigan State University, at GANIL, the French national nuclear physics laboratory, and the University of Notre Dame. At ReA3 the degree to which α-clusters associated with fusion of 28,30,32 Si + 28 Si result from the collision dynamics or reflect an initial α-cluster structure was explored. Alpha clusters observed in fusion reactions exceed the predictions of the standard statistical model. This experiment provides a measure of how clusterization is impacted by the increasing neutron-richess of the system. In the experiment at GANIL we investigated fusion in 19 O + 12 C and 20 O + 12 C. This experiment utilized the recently developed active-target detector MuSIC@Indiana. Using this proven, highly efficient, active-target detector we measured the fusion excitation function for these neutron-rich systems. Grounded by these experimental measurements, theoretical models were used to examine the role of unpaired valence neutrons on the fusion cross-section at energies just above the fusion barrier.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Data augmentation for disruption prediction via robust surrogate models

The goal of this work is to generate large statistically representative data sets to train machine learning models for disruption prediction provided by data from few existing discharges. Such a comprehensive training database is important to achieve satisfying and reliable prediction results in artificial neural network classifiers. Here, we aim for a robust augmentation of the training database for multivariate time series data using Student t process regression. We apply Student t process regression in a state space formulation via Bayesian filtering to tackle challenges imposed by outliers and noise in the training data set and to reduce the computational complexity. Thus, the method can also be used if the time resolution is high. We use an uncorrelated model for each dimension and impose correlations afterwards via colouring transformations. We demonstrate the efficacy of our approach on plasma diagnostics data of three different disruption classes from the DIII-D tokamak. To evaluate if the distribution of the generated data is similar to the training data, we additionally perform statistical analyses using methods from time series analysis, descriptive statistics and classic machine learning clustering algorithms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Robust System Separation Strategy Considering Online Wide-Area Coherency Identification and Uncertainties of Renewable Energy Sources

With the fast growth of renewable energy sources (RES), more and more uncertainties are involved and influencing the stable operation of power systems. Controlled islanding is the last measure to prevent power system blackouts, thus this paper aims to propose a novel model of system separation based on Online Coherency Identification and Adjustable Robust Optimization Programming (OCI-AROP) for minimizing load shedding considering the uncertainties of RES. First, Fuzzy C-Means (FCM) clustering method with F-statistics is utilized to identify the coherent generator groups with the frequency data measured by Phasor Measurement Units (PMUs). Then, the OCI-AROP model considering coherent group constraints, connectivity constraints and robustness constraints about RES are presented. Finally, the case studies on IEEE-39 bus system and WECC-179 bus system are employed to demonstrate the effectiveness of the proposed OCI-AROP model, and comparisons among the OCI-AROP model and the other models are also given to show its superiority.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Data augmentation for disruption prediction via robust surrogate models

The goal of this work is to generate large statistically representative datasets to train machine learning models for disruption prediction provided by data from few existing discharges. Such a comprehensive training database is important to achieve satisfying and reliable prediction results in artificial neural network classifiers. Here, we aim for a robust augmentation of the training database for multivariate time series data using Student-t process regression. We apply Student-t process regression in a state space formulation via Bayesian filtering to tackle challenges imposed by outliers and noise in the training data set and to reduce the computational complexity. Thus, the method can also be used if the time resolution is high. We use an uncorrelated model for each dimension and impose correlations afterwards via coloring transformations. We demonstrate the efficacy of our approach on plasma diagnostics data of three different disruption classes from the DIII-D tokamak. To evaluate if the distribution of the generated data is similar to the training data, we additionally perform statistical analyses using methods from time series analysis, descriptive statistics, and classic machine learning clustering algorithms.

97 MATHEMATICS AND COMPUTING↗

Modeling of H2 Dispersion at ARIES

Hydrogen is a versatile and clean energy carrier that can be produced from various renewable sources such as wind, solar, and hydropower and help decarbonize electricity grids, industry, and transportation. Using the Hydrogen Research Facility under Advanced Research on Integrated Energy Systems (ARIES) at the National Renewable Energy Laboratory's (NREL) Flatirons campus as a test bench, the study examines the feasibility, useability, and value of using computational fluid dynamics (CFD) techniques to model hydrogen dispersion. The ARIES facility was chosen because controlled hydrogen releases can be performed at a rate of 27 kg-H2/hr. Site-specific atmospheric and weather condition data such as wind speed and temperature were used as inputs to the model. The results show statistical distributions and ranges of hydrogen concentrations at locations throughout the domain. Wind conditions are found to significantly impact the release behavior, including the hydrogen cloud's direction and concentrations. At low wind speeds (below 1 mph), hydrogen forms a cloud and at higher wind speeds (> 2-4 mph) hydrogen plume stretches in the direction of wind momentum. From >100 simulations for ARIES site-specific conditions, statistical quantities combined with a clustering algorithm were used to propose sensor location at various elevations from ground.

dispersion↗

Segmental Dynamics and Vitrification in Associating Copolymer Melts: Role of Cluster Formation, Microdomains, and Cross-Linking

A microscopic statistical mechanical theory of the structure, self-assembly, and activated segmental relaxation is employed to study associating copolymer melts with high attractive sticker fractions, local clustering, and disordered microphase ordering. The stickers are dynamically pinned in a manner that does not affect equilibrium structure which mimics the much slower physical bond breaking process or postassembly cross-linking of sticky monomers. Local sticker clustering and microdomain spatial correlations significantly modify the activated relaxation of nonstickers and glass transition temperature, T g . A re-entrant glass-melting feature is predicted as sticker attraction strength is initially increased corresponding to a speed up of segmental relaxation, and hence reduction of T g relative to the cross-linked homopolymer network. A mechanistic analysis reveals three competing effects: a purely kinetic slowing down of nonstickers down due to cross-linking, disordering of the nonsticker local cage and weakening of effective forces they experience due to sticker physical clustering, and a longer range impact of microdomain scale correlations that results in nonmonotonic dynamical effects. At high enough attraction strength, a qualitative change emerges corresponding to a sticker fraction dependent elevation of T g , which eventually surpasses that of the cross-linked homopolymer network. Dynamically, the new physics arises from a complex evolution of the amplitude of the collective elastic field that dresses the large amplitude mobile segment hopping within a coupled local-nonlocal description of the alpha relaxation. Furthermore, the results are qualitatively consistent with recent experiments on associating PDMS and PPG telechelics of fixed sticker fraction but with chemically different end groups of variable attraction strengths. Finally, possible tests using simulation and the influence of material or model specific interaction potentials and other real world complications are discussed.

Copolymers↗

Subseasonal Clustering of Atmospheric Rivers Over the Western United States

Abstract The serial occurrence of atmospheric rivers (ARs) along the US West Coast can lead to prolonged and exacerbated hydrologic impacts, threatening flood‐control and water‐supply infrastructure due to soil saturation and diminished recovery time between storms. Here a statistical approach for quantifying subseasonal temporal clustering among extreme events is applied to a 41‐year (1979–2019) wintertime AR catalog across the western United States (US). Observed AR occurrence, compared against a randomly distributed AR timeseries with the same average event density, reveals temporal clustering at a greater‐than‐random rate across the western US with a distinct geographical pattern. Compared to the Pacific Northwest, significant AR clusters over the northern Coastal Range of California and Sierra Nevada are more frequent and occur over longer time periods. Clusters along the California Coastal Range typically persist for 2 weeks, are composed of 4–5 ARs per cluster, and account for over 85% of total AR occurrence. Across the northwest Coast‐Cascade Ranges, clusters account for ∼50% of total AR occurrence, typically last 8–10 days, and contain 3–4 individual AR events. Based on precipitation data from a high‐resolution dynamical downscaling of reanalysis, the fractions of total and extreme hourly precipitation attributable to AR clusters are largest along the northern California coast and in the Sierra Nevada. Interannual variability among clusters highlights their importance for determining whether a particular water year is anomalously wet or dry. The mechanisms behind this unusual clustering are unclear and require further research.

Meteorology & Atmospheric Sciences↗

Measuring Nuclear Clusters in the Short-Baseline Near Detector

The Short-Baseline Near Detector (SBND) is the first Liquid Argon Time Projection Chamber (LArTPC) with high enough resolution and large enough neutrino flux to measure the cross-section for production of heavier-than-proton fragments in neutrino interactions. SBND is an 112-ton LArTPC, and lies 110 m downstream from the Booster Neutrino Beam target, where it is collecting more than two million neutrino interactions per year. As a result of both the intranuclear cascade and nuclear de-excitation, neutrino interactions with argon can produce nuclear clusters such as deuterons, tritons, helions and alpha particles, which have so far only been seen as vertex activity or mis-reconstructed as protons. With LArTPC technology's excellent reconstruction capabilities and SBND's unprecedented neutrino statistics, a measurement of these nuclear clusters could distinguish between nuclear models and improve neutrino energy reconstruction. This poster presents the status of simulation and reconstruction of nuclear clusters in SBND.

Beever, Anna L. [Unlisted, US] (ORCID:000900069339↗

Beyond Limber: efficient computation of angular power spectra for galaxy clustering and weak lensing

Angular two-point statistics of large-scale structure observables are important cosmological probes. To reach the high accuracy required by the statistical precision of future surveys, some of these statistics may need to be computed without the commonly employed Limber approximation; the exact computation however requires integration over Bessel functions, and a brute-force evaluation is slow to converge. Here, we present a new method based on our generalized FFTLog algorithm for the efficient computation of angular power spectra beyond the Limber approximation. The new method significantly simplifies the calculation and improves the numerical speed and stability. It is easily extended to handle integrals involving derivatives of Bessel functions, making it equally applicable to numerically more challenging cases such as contributions from redshift-space distortions and Doppler effects. We implement our method for galaxy clustering and galaxy-galaxy lensing power spectra. We find that using the Limber approximation for galaxy clustering in future analyses like LSST Year 1 and DES Year 6 may cause significant biases in cosmological parameters, indicating that going beyond the Limber approximation is necessary for these analyses.

weak gravitational lensing↗

Effects of cluster expansion on the locations of phase transition boundary as a first step to quantify uncertainty in first principles statistical mechanics framework

Predicting phase diagrams from first principle calculations eliminates the need of tedious experimental trials and errors. Fully automating first principle phase diagram calculations without any sort of human intervention has been a long daunting task and troubling scientific communities for decades. This grand problem remains not fully resolved, largely due to the vastly high-dimensional parameter space associated with density functional theory, cluster expansion, lattice Monte Carlo, and the substantial uncertainty propagating through a set of complex simulations. As a first step to tackle this grand problem, we reported a first demonstration of how sensitive phase boundary locations can be to various cluster expansion fittings and input DFT training data. To the best knowledge of the authors, this study reported the first ever attempt to quantify uncertainty in first principles statistical mechanics framework. In addition, a semi-automated phase transition detection algorithm has been devised in this paper to deal with the associated statistical errors and uncertainties from Monte Carlo method and its predecessor DFT calculations and cluster expansions. This algorithm has been applied in a classical cluster expansion Mg-Cd binary alloy system to detect phase transitions at various locations of phase diagram using different cluster expansions to demonstrate its predictive power and quantify uncertainties. The results suggested that using Chebyshev basis function shifted the transition locations toward the dilute solid solution phase and expanded the phase stability of concentrated ordered phases to wider composition range. The addition of perturbed defect configurations into training data set lowered the order-disorder transition temperature to be away from true transition temperature, suggesting the transition nature is indeed configurational disorder dominated rather than defect assisted. Finally, we found that the weighting has negligible effect on the transition locations, except for the case of using Chebyshev basis function to fit all non-weighted configurations that can be susceptible to Monte Carlo sampling hysteresis.

36 MATERIALS SCIENCE↗

Does a Free Electron Laser Exhibit Non-Standard Statistics? [Slides and Report]

The originally proposed work under this award sought to develop a theoretical basis for previous experimental observations in free-electron laser (FEL) physics involving certain non-classical effects in the light produced by these devices. A sound and valid theoretical interpretation of these experiments would represent a transformational understanding of free-electron lasers, and of electron-photon interactions in general. The specific objective of the completed work, as described in the attached progress report by Jeongwan Park, was to perform a detailed analysis of an experiment by Chen and Madey (Physical Review Letters; volume 86 number 26, 2001) concerning the purported sub-Poissonian photon statistics in an FEL. Since that published work contradicted FEL theory as presently understood, the purpose of the work under the present subcontract was to obtain a better understanding of the quantum nature of the FEL, to investigate the possibility of theoretical and experimental evidence for the non-standard photon statistics of the FEL, and to critically examine the experimental evidence through a re-evaluation of the data analysis of Chen and Madey. The results of this study showed that there were numerous experimental conditions under which the observed photon statistics could be explained by combining both the photon clustering property and the dead-time effect, and consequently, that the observed photon statistics in the Chen-Madey experiment did not require sub-Poissonian photon statistics for their explanation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Investigating the influence of particle distribution on force and torque statistics using hierarchical machine learning

An accurate representation of hydrodynamic force and torque experienced by every particle in a distribution can be obtained from particle resolved (PR) simulations. These unique quantities are influenced by the deterministic position of surrounding particles. However, systems simulated with this methodology are typically limited to particles due to the involved computational cost. This resource requirement is a major bottleneck in analyzing the effect of variations in particle distribution. Here, this article attempts to address this bottleneck by availing relatively inexpensive deep learning models. The surrogate models that we employ in this article use a physics‐based hierarchical framework and symmetry‐preserving neural networks to achieve robustness with limited training data. This article first performs additional generalizability tests on PR data of distinct distributions that are not involved in the training process. The models are then deployed on several different particle distributions. Impact of clustering and structure on the observed statistics are investigated.

42 ENGINEERING↗

SUNBIRD : a simulation-based model for full-shape density-split clustering

Combining galaxy clustering information from regions of different environmental densities can help break cosmological parameter degeneracies and access non-Gaussian information from the density field that is not readily captured by the standard two-point correlation function (2PCF) analyses. However, modelling these density-dependent statistics down to the non-linear regime has so far remained challenging. We present a simulation-based model that is able to capture the cosmological dependence of the full shape of the density-split clustering (DSC) statistics down to intra-halo scales. Our models are based on neural-network emulators that are trained on high-fidelity mock galaxy catalogues within an extended-ΛCDM framework, incorporating the effects of redshift-space, Alcock–Paczynski distortions, and models of the halo–galaxy connection. Our models reach sub-percent level accuracy down to $1 \, h^{-1}\text{Mpc}$ and are robust against different choices of galaxy–halo connection modelling. When combined with the galaxy 2PCF, DSC can tighten the constraints on ω cdm , σ 8 , and n s by factors of 2.9, 1.9, and 2.1, respectively, compared to a 2PCF-only analysis. DSC additionally puts strong constraints on environment-based assembly bias parameters.

79 ASTRONOMY AND ASTROPHYSICS↗

Large‐Scale Statistically Meaningful Patterns (LSMPs) Associated With Precipitation Extremes Over Northern California

Abstract We analyze large‐scale statistically meaningful patterns (LSMPs) that precede extreme precipitation (PEx) events over Northern California (NorCal). We find LSMPs by applying k‐means clustering to the two leading principal components of daily 500 hPa geopotential height anomalies two days before the onset, from October to March during 1948–2015. Statistical significance testing based on Monte Carlo simulations suggests a minimum of four statistically distinguished LSMP clusters. The four LSMP clusters are characterized as Northwest continental negative height anomaly, Eastward positive “Pacific‐North American Pattern (PNA),” Westward negative “PNA,” and Prominent Alaskan ridge. These four clusters, shown in multiple variables, evolve very differently and have differing links to the Arctic and tropical Pacific regions. Using binary forecast skill measures and a new copula‐based framework for predicting PEx events, we find LSMP indices that are useful predictors of NorCal PEx events, with moisture‐based variables being the best predictors of PEx events at least 6 days before the onset, and the lower atmospheric variables being better than their upper atmospheric counterparts any day in advance tested. To ensure statistical rigor, the LSMPs analyzed here (with the modified acronym) include local tests of both significance and consistency, which are not always featured in the literature on large‐scale meteorological patterns.

54 ENVIRONMENTAL SCIENCES↗

HFTS-2 Completions Design and State-of-the-Art Diagnostics Results

Various completion schemes were tested as part of the Hydraulic Fracture Test Site-2 (HFTS-2) project. The wells were instrumented with permanent fiber optic (FO) cable, and unique near-wellbore and far-field region surveys were acquired to evaluate spatial and temporal hydraulic fracture characteristics. The state-of-the-art diagnostics included substantial monitoring of both the well stimulation and longer-term production performance. The objective was to evaluate and improve stimulation distribution effectiveness (SDE) for different completion designs based on the diagnostics results. A hydraulic fracturing (HF) profile (i.e., proppant/fluid allocation among clusters) was obtained using near-wellbore Distributed Acoustic Sensing (DAS) data acquired during stimulation. Statistical analysis provided a quantitative assessment of SDE for the different completion designs. A novel workflow was developed to assess time-lapse fracture property changes at the cluster and stage levels based on Distributed Strain Sensing (DSS) data acquired during the production phase. These time-lapse fracture property changes were compared with the HF profile. Based on the statistical analysis, Rate Per Cluster (RPC) during proppant slurry placement appears to be a primary completion design variable, rather than the typical Rate Per Perf (RPP). If RPC was not maintained at a certain level, the likelihood of uneven distribution increased dramatically. The stage configurations that created low RPC while maintaining high RPP did not show improved SDE. The strain change results based on DSS data during the production phase indicated near-wellbore fracture property changes (e.g., aperture changes during well shut-in and flow periods) and aligned well with near-wellbore DAS signals acquired during stimulation. More importantly, fracture property changes based on DSS data acquired during two tests over a 7-month period show high consistency and good alignment with the near-wellbore HF profile. This workflow quantitatively evaluates SDE for different completion designs using diagnostic results and could help optimize completion strategy for future development projects. Here, the analysis and interpretation of results demonstrate the technical feasibility of using permanent fiber optic cable for stimulation and production monitoring and indic

58 GEOSCIENCES↗

Multiscale spatial analysis of fracture arrangement and pattern reconstruction using Ripley's K-function

This work presents novel multiscale spatial data analytics using Ripley's K-function, as a measure of spatial interaction, to study one-dimensional arrangement of fractures. Fracture spatial arrangements are classified into clustered, anticlustered, or indistinguishable from random by testing statistical significance of the calculated Ripley's K-function. Characterizations of fracture arrangements are performed as a function of length scale and position. Analysis of the K-function along the study interval identifies where fracture clustering and anticlustering occur. Here, a simulation technique is also introduced here to statistically reconstruct spatial arrangements and to generate fracture realizations that are spatially similar to the fractures observed in the field. With this simulation technique, one can also fill spatial gaps in fracture measurements where data are absent, unreliable, or unused. Synthetic as well as field-measured 1D fracture datasets are used for testing and demonstration. Methods introduced in this work can be readily applied to fracture datasets observed in outcrops, borehole image logs, and cores.

58 GEOSCIENCES↗