Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Random Effects”

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 739 records · Page 41

W$^+$W$^-$ boson pair production in proton-proton collisions at $\sqrt{s} =$ 13 TeV

A measurement of the W+W- boson pair production cross section in proton-proton collisions at s=13 TeV is presented. The data used in this study are collected with the CMS detector at the CERN LHC and correspond to an integrated luminosity of 35.9 fb-1. The W+W- candidate events are selected by requiring two oppositely charged leptons (electrons or muons). Two methods for reducing background contributions are employed. In the first one, a sequence of requirements on kinematic quantities is applied allowing a measurement of the total production cross section, 117.6±6.8 pb, which agrees well with the theoretical prediction. Fiducial cross sections are also reported for events with zero or one jet, and the change in the zero-jet fiducial cross section with the jet transverse momentum threshold is measured. Normalized differential cross sections are reported within the fiducial region. A second method for suppressing background contributions employs two random forest classifiers. The analysis based on this method includes a measurement of the total production cross section and also a measurement of the normalized jet multiplicity distribution in W+W- events. Finally, a dilepton invariant mass distribution is used to probe for physics beyond the standard model in the context of an effective field theory, and constraints on the presence of dimension-6 operators are derived.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Estimating Subhourly Inverter Clipping Loss From Satellite-Derived Irradiance Data

Photovoltaic system production simulations are conventionally run using hourly weather datasets. Hourly simulations are sufficiently accurate to predict the majority of long-term system behavior but cannot resolve high-frequency effects like inverter clipping caused by short-duration irradiance variability. Direct modeling of this subhourly clipping error is only possible for the few locations with high-resolution irradiance datasets. This paper describes a method of predicting the magnitude of this error using a machine learning regressor ensemble model, comprised of a random forest and an XGBoost model, and 30-minute satellite irradiance data. The method predicts a correction for each 30-minute interval with the potential to roll up into 60-minute corrections to match an hourly energy model. The model is trained and validated at locations where the error can be directly simulated from 1-minute ground data. The validation shows low bias at most ground station locations. The model is also applied to gridded satellite irradiance to produce a heatmap of the estimated clipping error across the United States. Finally, the relative importance of each predictor satellite variable is retrieved from the model and discussed.

41 EE - Solar Energy Technologies Office (EE-4S)↗

Using Behavioral Science to Target LMI and High-Value Solar Installations (Final Technical Report)

This project is a data-driven analysis of approaches to facilitate solar adoption in LMI communities and areas with high potential benefits to the electric grid. The cornerstone of the project was two waves of randomized field experiments. The first wave involved novel approaches to reaching LMI communities based on insights from literature. The second wave involved areas of high potential of solar to the electricity grid, as determined in collaboration with the local utilities in Connecticut. Both waves were followed by surveys to develop insight into the effectiveness of different approaches. A third key element of the project is the development of new agent-based models of solar energy diffusion that explicitly incorporate the real spatial and temporal constraints that influence solar diffusion, based on real-world data on solar adoption. The project also involved outreach and dissemination of key findings to relevant stakeholders, including a guidebook on strategies to accelerate solar deployment to low and moderate income communities, multiple academic publications, and numerous presentations of the research findings.

14 SOLAR ENERGY↗

Evolution of an Asteroid Family under YORP, Yarkovsky, and Collisions

Any population of asteroids, like asteroid families, will disperse in semimajor axis due to the Yarkovsky effect. The amount of drift is modulated by the asteroid spin state evolution, which determines the balance between the diurnal and seasonal Yarkovsky forces. The asteroid’s spin state is, in turn, controlled in part by the Yarkovsky–O’Keefe–Radzievskii–Paddack (YORP) effect. The otherwise smooth evolution of an asteroid can be abruptly altered by collisions, which can cause impulsive changes in the spin state and can move the asteroid onto a different YORP track. In addition, collisions may also alter the YORP parameters by changing the superficial features and overall shape of the asteroid. Thus, the coupling between YORP and Yarkovsky is also strongly affected by the impact history of each body. To investigate this coupling, we developed a statistical code modeling the time evolution of semimajor axis under YORP–Yarkovsky coupling. It includes the contributions of NYORP (normal YORP), TYORP (tangential YORP), and collisions whose effects are deterministically calculated and not added in a statistical way. We find that both collisions and TYORP increase the dispersion of a family in semimajor axis by making the spin axis evolution less smooth and regular. We show that the evolution of a family’s structure with time is complex and collisions randomize the YORP evolution. In our test families, we do not observe the formation of a “YORP-eye” in the semimajor axis versus diameter distribution, even after a long period of time. If present, the “YORP-eye” might be a relic of an initial ejection velocity pattern of the collisional fragments.

79 ASTRONOMY AND ASTROPHYSICS↗

On the Marriage of Asynchronous Many Task Runtimes and Big Data: A Glance

The rise of the accelerator-based architectures and reconfigurable computing have showcased the weakness of software stack toolchains that still maintain a static view of the hardware instead of relying on a symbiotic relationship between static (e.g., compilers) and dynamic tools (e.g., runtimes). In the past decades, this need has given rise to adaptive runtimes with increasingly finer computational tasks. These finer tasks help to take advantage of the hardware by switching out when a long latency operation is encountered (because of the deeper memory hierarchies and new memory technologies that might target streaming instead of random access), thus trading off idle time for unrelated work. Examples of these finer task runtimes are Asynchronous Many Task (AMT) runtimes, in which highly efficient computational graphs run on a variety of hardware. Due to its inherent latency tolerant characteristics, Latency-sensitive applications, such as Graph Analytics and Big Data can effectively use these runtimes. This paper aims to present an example of how the careful design of an AMT can exploit the hardware substrate when faced with high latency applications such as the ones given in the Big Data domain. Moreover, with its introspection and adaptive capabilities, we aim to show the power of these runtimes when facing the changing requirements of the application workloads. We use the Performance Open Community Runtime (P-OCR) as our vehicle to demonstrate the concepts presented here.

adaptive runtime, big data analysis↗

Synthetic High Impedance Fault Data through Deep Convolutional Generated Adversarial Network

High impedance faults (HIFs) have always been significant challenge in the power grids. Researchers have developed some advanced protective methods to detect the HIFs. To test and validate these methods, large amounts of HIF data are required. This paper presents a synthetic HIF data generating method using the deep convolutional generated adversarial network (DCGAN). The DCGAN includes a generator module to create synthetic HIF waveform from random noises; and a discriminator module to identify the flaws of those synthetic data, which ultimately help improve the quality of the synthetic data created by the generator. To test the fidelity of the generated synthetic HIF data, two different HIF-detection methods have been applied. Extensive simulation results have validated the effectiveness of using the DCGAN to create synthetic HIF data.

Yang, Kun↗

Comprehensive Analysis of the Neutrino Process in Core-collapsing Supernovae

We investigate the neutrino flavor change effects due to neutrino self-interaction and shock wave propagation, as well as the matter effects on the neutrino process in core-collapsing supernovae (CCSNe). For the hydrodynamics, we use two models: a simple thermal bomb model and a specified hydrodynamics model for SN1987A. For the presupernova model, we take an updated model, adjusted to explain SN1987A, which employs recent developments in the (n, γ) reaction rates for nuclei near the stability line (A ~ 100). As for the neutrino luminosity, we adopt two different models: equivalent neutrino luminosity and nonequivalent luminosity models. The latter is taken from a synthetic analysis of CCSN simulation data, which quantitatively presented the results obtained by various neutrino transport models. Relevant neutrino-induced reaction rates are calculated using a shell model for light nuclei and a quasiparticle random phase approximation model for heavy nuclei. For each model, we present abundances of the light nuclei ( 7 Li, 7 Be, 11 B, and 11 C) and the heavy nuclei ( 92 Nb, 98 Tc, 138 La, and 180 Ta) produced by the neutrino process. The light nuclei abundances turn out to be sensitive to the Mikheyev–Smirnov–Wolfenstein (MSW) region around O-Ne-Mg layer while the heavy nuclei are mainly produced prior to the MSW region. Through detailed analyses, we find that neutrino self-interaction becomes a key ingredient, in addition to the MSW effect, for understanding the neutrino process and the relevant nuclear abundances. The normal mass hierarchy is shown to be more compatible with the meteorite data. The main nuclear reactions for each nucleus are also investigated in detail.

79 ASTRONOMY AND ASTROPHYSICS↗

Short-Depth QAOA circuits and Quantum Annealing on Higher-Order Ising Models (Rev.2)

The Quantum Alternating Operator Ansatz (QAOA) and Quantum Annealing (QA) are quantum algorithms that are both based on the adiabatic theorem and both have the goal of sampling the optimal solution(s) of combinatorial optimization problems. Quantum annealing has been physically instantiated on D-Wave devices using superconducting flux qubits, and QAOA can be programmed on digital gate-model quantum computers such as the programmable superconducting transmon qubits devices of the IBMQ series, for instance ibm washington. QAOA and QA address the same types of problems, but it is unclear how they will scale to large problem sizes and to larger and higher-fidelity quantum computers. In this article, we present a direct comparison between QAOA, one and two rounds, run on all 127 qubits of ibm washington and QA run on D-Wave Advantage system4.1 and Advantage system6.1. The problems which allow for this comparison are random Ising model problems whose connectivity matches the heavy hexagonal lattice topology of ibm washington and the Pegasus graph connectivity of the two D-Wave devices. We create two classes of problem instances for this comparison: one with higher order terms (ZZZ variable interactions), linear terms, and quadratic terms, and a separate problem type with only linear and quadratic terms. Our QAOA circuits are novel and extremely short depth, with a CNOT depth of 6 per round, which allows whole chip usage of ibm washington’s heavy hexagonal lattice and can be applied to future heavy-hex chips. We also test the effectiveness of the error suppression technique digital dynamical decoupling on the QAOA circuits. The QAOA circuits compiled to ibm washington are composed of several thousand circuit instructions, approximately 3, 000 depending on the details of the circuit, making these some the largest quantum circuits ever executed on a digital quantum processor. QAOA and QA are compared against the classical heuristic algorithm of simulated annealing and all problem instances are exactly solved using CPLEX in order to evaluate which samplers, if any, correctly found the ground state solution(s) of the problem instances. We find that (i) QA outperforms QAOA on all problem instances, (ii) QAOA samples the problems better than random sampling, and (iii) QAOA angle computation exhibits clear parameter concentration across the ensemble of Ising models.

127 Qubits↗

Unveiling the drivers contributing to global wheat yield shocks through quantile regression

Sudden reductions in crop yield (i.e., yield shocks) severely disrupt the food supply, intensify food insecurity, depress farmers' welfare, and worsen a country's economic conditions. Here, we study the spatiotemporal patterns of wheat yield shocks, quantified by the lower quantiles of yield fluctuations, in 86 countries over 30 years. Furthermore, we assess the relationships between shocks and their key ecological and socioeconomic drivers using quantile regression based on statistical (linear quantile mixed model) and machine learning (quantile random forest) models. Using a panel dataset that captures spatiotemporal patterns of yield shocks and possible drivers in 86 countries, we find that the severity of yield shocks has been increasing globally since 1997. Moreover, our cross-validation exercise shows that quantile random forest outperforms the linear quantile regression model. Despite this performance difference, both models consistently reveal that the severity of shocks is associated with higher weather stress, nitrogen fertilizer application rate, and gross domestic product (GDP) per capita (a typical indicator for economic and technological advancement in a country). While the unexpected negative association between more severe wheat yield shocks and higher fertilizer application rate and GDP per capita does not imply a direct causal effect, they indicate that the advancement in wheat production has been primarily on achieving higher yields and less on lowering the possibility and magnitude of sharp yield reductions. Hence, in the context of growing extreme weather stress, there is a critical need to enhance the technology and management practices that mitigate yield shocks to improve the resilience of the world food systems.

60 APPLIED LIFE SCIENCES↗

Exact bounds for dynamical critical exponents of transverse-field Ising chains with a correlated disorder

Highlights: • Analytic results in disordered quantum systems are obtained. • Correlated disorder changes the universality class of quantum phase transitions. • The lower and upper bounds of the dynamical critical exponent are obtained. This study investigates the dynamical critical exponent of disordered Ising chains under transverse fields to examine the effect of a correlated disorder on quantum phase transitions. The correlated disorder, where the on-site transverse field depends on the nearest-neighbor coupling strengths connecting the site, gives a qualitatively different result from the uncorrelated disorder. In the uncorrelated disorder cases where the transverse field is either homogeneous over sites or random independently of the nearest-neighbor coupling strengths, the dynamical critical exponent is infinite. In contrast, in the presence of the correlated disorder, we analytically show that the dynamical critical exponent is finite. We also show that the dynamical critical exponent depends on the tuning process of the transverse field strengths.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Land surface dynamics and meteorological forcings modulate land surface temperature characteristics

This study examines the effect of land cover, vegetation health, climatic forcings, elevation heat loads, and terrain characteristics (LVCET) on land surface temperature (LST) distribution over West Africa (WA). We employ fourteen machine-learning models, which preserve nonlinear relationships, to downscale LST and other predictands while preserving the geographical variability of WA. Our results showed that the random forest model performs best in downscaling predictands. This is important for the sub-region since it has limited access to mainframes to power multiplex machine-learning algorithms. In contrast to the northern regions, the southern regions consistently exhibit healthy vegetation. Also, areas with unhealthy vegetation coincide with hot LST clusters. The positive Normalized Difference Vegetation Index (NDVI) trends in the Sahel underscore rainfall recovery and subsequent Sahelian greening. The southwesterly winds cause the upwelling of cold waters, lowering LST in southern WA and highlighting the cooling influence of water bodies on LST. Identifying regions with elevated LST is paramount for prioritizing greening initiatives, and our study underscores the importance of considering LVCET factors in urban planning. Topographic slope-facing angles, heat loads, and diurnal anisotropic heat all contribute to variations in LST, emphasizing the need for a holistic approach when designing resilient and sustainable landscapes.

54 ENVIRONMENTAL SCIENCES↗

Scaling Field-Theoretic Simulation for Multicomponent Mixtures with Neural Operators

Multicomponent polymer mixtures are ubiquitous in biological self-organization but are notoriously difficult to study computationally. Plagued by both slow single molecule relaxation times and slow equilibration within dense mixtures, molecular dynamics simulations are typically infeasible at the spatial scales required to study the stability of mesophase structure. Polymer field theories offer an attractive alternative, but analytical calculations are only tractable for mean-field theories and nearby perturbations, constraints that become especially problematic for fluctuation-induced effects such as coacervation. Here, we show that a recently developed technique for obtaining numerical solutions to partial differential equations based on operator learning, neural operators, lends itself to a highly scalable training strategy by parallelizing per-species operator maps. We illustrate the efficacy of our approach on six-component mixtures with randomly selected compositions and that it significantly outperforms the state-of-the-art pseudospectral integrators for field-theoretic simulations, especially as polymer lengths become long.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Enhanced Estimation of Quantum Properties with Common Randomized Measurements

We present a technique for enhancing the estimation of quantum state properties by incorporating approximate prior knowledge about the quantum state. This consists in performing randomized measurements on a quantum processor and comparing the results with those obtained from a classical computer that stores an approximation of the quantum state. We provide unbiased estimators for expectation values of multicopy observables and present performance guarantees in terms of variance bounds that depend on the prior knowledge accuracy. We demonstrate the effectiveness of our approach through experimental and numerical examples detecting mixed-state entanglement, and estimating polynomial approximations of the von Neumann entropy and state fidelities. Published by the American Physical Society 2024

Physics↗

Effective Local Geometry Descriptor for 29Si NMR Q4 Anisotropy

The nuclear shielding anisotropy, ζ, is a useful nuclear magnetic resonance (NMR) shielding tensor parameter in describing the extent of electron cloud distortion about an atom. Despite the advantages afforded by NMR in structural characterization, the relationship between ζ and local structure of an atom in high-symmetry environments, such as Si-Q 4 sites, is poorly understood. Here, we use a data-driven approach combining random forest feature ranking and the Sure Independence Screening and Sparsifying Operator (SISSO) approach to derive a simple and accurate geometric descriptor for ζ with a root-mean-squared prediction error of 6.77 ppm and an R 2 of 0.761. We then apply this descriptor to describe the local geometric distortion of zeolites Sigma-2 and silica-ZSM-5 whose chemical shift anisotropy tensor has been reported. We envision that this geometric descriptor will allow for structural description and refinement in previously difficult-to-describe materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Bootstrap estimation of the effect of instrument response function uncertainty on the reconstruction of fusion neutron sources

Neutron imagers are important diagnostics for the inertial confinement fusion implosions at the National Ignition Facility. They provide two- and three-dimensional reconstructions of the neutron source shape that are key indicators of the overall performance. To interpret the shape results properly, it is critical to estimate the uncertainty in those reconstructions. There are two main sources of uncertainties: limited neutron statistics, leading to random errors in the reconstructed images, and incomplete knowledge of the instrument response function (the pinhole-dependent point spread function). While the statistical errors dominate the uncertainty for lower yield deuterium-tritium (DT) shots, errors due to the instrument response function dominate the uncertainty for DT yields on the order of 10 16 neutrons or higher. In this work, a bootstrapping method estimates the uncertainty in a reconstructed image due to the incomplete knowledge of the instrument response function. The main reconstruction is created from the fixed collection of pinhole images that are best aligned with the neutron source. Additional reconstructions are then built using subsets of that collection of images. Variations in the shapes of these additional reconstructions originate solely from uncertainties in the instrument response function, allowing us to use them to provide an additional systematic uncertainty estimate.

47 OTHER INSTRUMENTATION↗

Molecular dynamics studies of helium bubble effects on grain boundary fracture vulnerabilities in an Fe 70 Ni 11 Cr 19 –1%H austenitic stainless steel

Comprehensive molecular dynamics tensile test simulations have been performed to study the delamination processes of seven different grain boundaries / cleavage planes (Σ1{111}, Σ3{111}, Σ5{100}, Σ7{111}, Σ9{411}, Σ11{311}, and R{100}/{411}) containing a helium bubble. Combinations of a variety of conditions are explored including different strain rates, system dimensions, bubble density, bubble radius, bubble pressure, and temperature. We found that in general, grain boundaries absorb less energies with decreasing strain rate but increasing bubble areal density, bubble pressure, bubble radius, and temperature. The propensity of grain boundary delamination is sensitive to grain boundary type: The random grain boundary R{100}/{411} is one of the most brittle boundaries whereas the Σ1{111} cleavage plane and the Σ3{111} twin boundary are two of the toughest boundaries. The sorted list of grain boundary fracture vulnerability obtained from our dynamic tensile test simulations differs from the one obtained from our decohesion energy calculations, confirming the important role of plastic deformation during fracture. Finally, detailed mechanistic analyses are performed to interpret the simulated results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

gLaSDI: Parametric physics-informed greedy latent space dynamics identification

A parametric adaptive physics-informed greedy Latent Space Dynamics Identification (gLaSDI) method is proposed for accurate, efficient, and robust data-driven reduced-order modeling of high-dimensional nonlinear dynamical systems. In the proposed gLaSDI framework, an autoencoder discovers intrinsic nonlinear latent representations of high-dimensional data, while dynamics identification (DI) models capture local latent-space dynamics. Here, an interactive training algorithm is adopted for the autoencoder and local DI models, which enables identification of simple latent-space dynamics and enhances accuracy and efficiency of data-driven reduced-order modeling. To maximize and accelerate the exploration of the parameter space for the optimal model performance, an adaptive greedy sampling algorithm integrated with a physics-informed residual-based error indicator and random-subset evaluation is introduced to search for the optimal training samples on the fly. Further, to exploit local latent-space dynamics captured by the local DI models for an improved modeling accuracy with a minimum number of local DI models in the parameter space, a -nearest neighbor convex interpolation scheme is employed. The effectiveness of the proposed framework is demonstrated by modeling various nonlinear dynamical problems, including Burgers equations, nonlinear heat conduction, and radial advection. The proposed adaptive greedy sampling outperforms the conventional predefined uniform sampling in terms of accuracy. Compared with the high-fidelity models, gLaSDI achieves 17 to 2,658× speed-up with 1 to 5% relative errors.

97 MATHEMATICS AND COMPUTING↗