Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “statistics of extremes”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Fluctuations of subsystem entropies at late times

We study the fluctuations of subsystem entropies in closed quantum many-body systems after thermalization. Using a combination of analytics and numerics for both random quantum circuits and Hamiltonian dynamics, we find that the statistics of such entropy fluctuations is drastically different than in the classical setting. For instance, shortly after a system thermalizes, the probability of entropy fluctuations for a subregion is suppressed in the dimension of the Hilbert space of the complementary subregion. This suppression becomes increasingly stringent as a function of time, ultimately depending on the exponential of the Hilbert space dimension, until extremely late times when the amount of suppression saturates. We also use our results to estimate the total number of rare fluctuations at large timescales. We find that the “Boltzmann brain” paradox is largely ameliorated in quantum many-body systems, in contrast with the classical setting.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Microphysical Sensitivity of Superparameterized Precipitation Extremes in the Contiguous United States Due to Feedbacks on Large-Scale Circulation

Superparameterized (SP) global climate models have been shown to better simulate various features of precipitation relative to conventional models, including its diurnal cycle as well as its extremes. While various studies have focused on the effect of differing microphysics parameterizations on precipitation within limited-area cloud-resolving models, we examine here the effect on contiguous U.S. (CONUS) extremes in a global SP model. We vary the number of predicted moments for hydrometeor distributions, the character of the rimed ice species, and the representation of raindrop self-collection and breakup. Using a likelihood ratio test and accounting for the effects of multiple hypothesis testing, we find that there are some regional differences, particularly during spring and summer in the Southwest and the Midwest, in both the current climate and a warmer climate with uniformly increased sea surface temperatures. These differences are most statistically significant and widespread when the number of moments is changed. To determine whether these results are due to (fast) local effects of the different microphysics or the (slower) ensuing feedback on the large-scale atmospheric circulation, we run a series of short, 5-day simulations initialized from reanalysis data. We find that the differences largely disappear in these runs and therefore infer that the different parameterizations impact precipitation extremes indirectly via the large-scale circulation. Finally, we compare the present-day results with hourly rain gauge data and find that SP underestimates extremes relative to observations regardless of which microphysics scheme is used given a fixed model configuration and resolution.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Using statistical modeling to predict and understand fusion experiments

Over 300 cryogenic layered direct-drive inertial confinement fusion implosions have been successfully executed on the OMEGA Laser System in the last decade. However, extracting sufficient understanding from these experiments to develop new designs or to identify or mitigate degradation sources remains challenging. Recently, a statistical modeling approach was developed to successfully design and predict improved implosion experiments on OMEGA. Here, we show that one-dimensional simulations can be used to predict the outcomes of systematically perturbed three-dimensional simulations and that this statistical modeling approach can be used to identify or rule out physical mechanisms for some of the degradation sources observed on the OMEGA Laser System for direct-drive cryogenic inertial confinement fusion. In this instance, we investigate the fusion yield dependencies on the ion temperature asymmetries and laser beam size observed in experiments by comparing with trends in three-dimensional synthetic simulation databases. Using the statistical model on these systematically perturbed simulations, we find that the statistically inferred dependency on the measured ion temperature asymmetries is well explained by imposed ℓ =1 modes. However, we find that the dependency on the laser beam size is only dominated by the illumination non-uniformity for some extreme cases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Open-Loop Control of Adjustable Tuned Mass Dampers for Floating Wind Turbine Platforms: Preprint

Floating offshore wind turbines have several advantages over their land-based counterparts, including access to stronger, more consistent winds and proximity to population centers. However, increased costs from the platform have been a challenge inhibiting their widespread adoption. New damping technologies with tunable target frequencies and damping ratios promise a greater degree of control over platform motions, allowing hulls to be designed smaller and reducing the overall cost of the turbines. In this work, a control strategy is proposed for these adjustable tuned mass dampers (TMDs). A frequency-domain model of the system is developed, from which response amplitude operators (RAOs) for the platform's rigid-body modes can be obtained for different natural frequencies and damping ratios of the dampers. Using these RAOs, and JONSWAP spectrums of operational and extreme sea states, the performance of the various damper settings are compared by evaluating the standard deviation of a rigid-body motion of interest (platform pitch or heave). By finding the optimal damper setting for a range of representative design load cases (DLCs), a lookup function is made to return the damper target frequency given the peak period of the sea state. To determine the current sea state, we propose a simple estimator that uses wave measurements (e.g., from a measurement buoy) to determine the peak period and significant wave height. The significant wave height is determined using the statistics of the past 100 seconds of wave elevation information and the peak period is computed using a frequency locked loop. Initial results show that we can use these estimated values to control the TMD natural frequency with an open loop controller. The controlled system is tested in a range of environmental conditions determined by the International Electrotechnical Commission design load cases (DLCs), which include normal and extreme wind and wave models. In these tests, we compare the performance of 4 cases: (1) no TMD, (2) a constant TMD based on the worst case DLC, (3) a controlled TMD based on known wind and wave environments, and (4) a real-time controlled TMD using estimated wind and wave environments. The effect of hull-based TMD control is also compared to changes in traditional wind turbine control via blade pitch.

control↗
DOE OSTI2025

Temperature measurement of Quark-Gluon plasma at different stages

In a Quark-Gluon Plasma (QGP), the fundamental building blocks of matter, quarks and gluons, are under extreme conditions of temperature and density. A QGP could exist in the early stages of the Universe, and in various objects and events in the cosmos. The thermodynamic and hydrodynamic properties of the QGP are described by Quantum Chromodynamics (QCD) and can be studied in heavy-ion collisions. Despite being a key thermodynamic parameter, the QGP temperature is still poorly known. Thermal lepton pairs (e + e − and μ + μ − ) are ideal penetrating probes of the true temperature of the emitting source, since their invariant-mass spectra suffer neither from strong final-state interactions nor from blue-shift effects due to rapid expansion. Here we measure the QGP temperature using thermal e+e− production at the Relativistic Heavy Ion Collider (RHIC). The average temperature from the low-mass region (in-medium ρ 0 vector-meson dominant) is (2.01 ± 0.23) × 10 12 K, consistent with the chemical freeze-out temperature from statistical models and the phase transition temperature from Lattice QCD. The average temperature from the intermediate mass region (above the ρ0 mass, QGP dominant) is significantly higher at (3.25 ± 0.60) × 10 12 K. This work provides essential experimental thermodynamic measurements to map out the QCD phase diagram and understand the properties of matter under extreme conditions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Spatiotemporal Learning in Power Modules: Wavelet-Enhanced Forecasting of Thermomechanical Degradation

Detecting internal defects in power electronics packages is critical for their performance and reliability, especially under extreme operating conditions, as these defects can lead to catastrophic failure if not properly addressed. Confocal scanning acoustic microscopy (C-SAM) plays a key role in the nondestructive evaluation of bond layer degradation within a power electronics package by detecting defects such as delamination, voids, and cracks. However, accurately quantifying and predicting these defects from C-SAM images remains a significant challenge due to the low noise-to-signal ratio, which typically arises from both imaging process and bond patterns itself. In this paper, we explore machine learning strategies for processing C-SAM images and providing predictive models of defect growth. We use C-SAM images of sintered copper and sintered silver samples, which are obtained under accelerated thermal experiments, as the representative dataset for our study. We investigate the effect of Fourier transforms and wavelet transforms on these datasets to remove high-frequency noise and address noise across multiple scales with histogram equalization to enhance the contrast and improve the visibility of defects. As a result, defect boundaries can be clearly distinguished, enabling more accurate tracking of their growth over time. We then employ different time-series forecasting algorithms on the denoised images to formulate an image-based lifetime prediction model. Statistical models and deep-learning techniques are trained on images obtained in the early stages of thermal shock, and defect growth in the later stages is predicted. Our work serves as a preliminary attempt to improve the accuracy of lifetime prediction models of power electronics packages, which is critical under extreme operating environments.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Exploring topography-based methods for downscaling subgrid precipitation for use in Earth System Models

Topography exerts major control on land surface processes. To improve representation of topographic impacts on land surface processes, a new topography-based subgrid structure has been introduced to the Energy Exascale Earth System Model to represent the subgrid heterogeneity of surface elevation. To take advantage of the new subgrid structure for improving land surface modeling, this study explores four variations of topography-based methods for downscaling grid precipitation to the corresponding subgrids. In the first three methods, the deviation of the subgrid precipitation from the grid’s is equal to the grid precipitation multiplied by the ratio of the elevation difference between the subgrid and grid mean to a specified elevation equals to the grid elevation, the difference between the maximum and minimum subgrid elevation, and the maximum subgrid elevation, respectively. The second method limits the ratio to 0.5 to avoid extreme values on mountains and the third method accounts for the slope effect. The fourth method is similar to the third method except that the Froude Number is used to limit the numerator in a blocking regime of the ambient flow. The downscaled precipitation is evaluated using the PRISM precipitation data over the U.S. using statistical metrics. Results show that by accounting for topographic slope besides elevation, the third and fourth methods show clear advantages over the first and second methods. Furthermore, introducing the Froude Number in the fourth method improves downscaling skill and shows consistent advantages over the third method in areas with larger subgrid heterogeneity across different grids sizes.

Tesfa, Teklu K.↗
DOE OSTI2023

An upper bound for extreme temperatures over midlatitude land

Heatwaves damage societies worldwide and are intensifying with global warming. Several mechanistic drivers of heatwaves, such as atmospheric blocking and soil moisture-atmosphere feedback, are well-known for their ability to raise surface air temperature. However, what limits the maximum surface air temperature in heatwaves remains unclear; this became evident during recent Northern Hemisphere heatwaves which achieved temperatures far beyond the upper tail of the observed statistical distribution. Here, we present evidence for the hypothesis that convective instability limits annual maximum surface air temperatures (TXx) over midlatitude land. We provide a theory for the corresponding upper bound of midlatitude temperatures, which accurately describes the observed relationship between temperatures at the surface and in the midtroposphere. We show that known heatwave drivers shift the position of the atmospheric state in the phase space described by the theory, changing its proximity to the upper bound. This theory suggests that the upper bound for midlatitude TXx should increase 1.9 times as fast as 500-hPa temperatures at the time and location of TXx occurrences. Using empirical 500-hPa warming, we project that the upper bound of TXx over Northern Hemisphere midlatitude land (40°N to 65°N) will increase about twice as fast as global mean surface air temperature, and TXx will increase faster than this bound over regions that dry on the hottest days.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Self-Organizing Map-Based Resilience Quantification and Resilient Control of Distribution Systems Under Extreme Events

Due to climate change, extreme weather events are occurring more frequently and with increasing impact. This trend poses a significant challenge for distribution system operators (DSO) to ensure that there is uninterrupted power supply to critical loads in their networks. To embed resilience into DSO's decision-making, resilience needs to be first quantified and then integrated into the system-level optimization. Therefore, this paper first develops a novel self-organizing map (SOM) based method (called SomRes) to quantify the time-varying resilience index of a system that can leverage the powerful classification property of SOMs and removes some of the disadvantages of subjective weight assignment methods. Using SomRes, a resilient resource allocation and operational dispatch algorithm is further developed to enhance system resilience against extreme events by considering the SomRes resilience index directly as the feedback. Here, the proposed resilience quantification approach is benchmarked with a state-of-the-art approach and the efficacy of the proposed resilient dispatch algorithm is demonstrated through several deterministic and statistical case studies on the IEEE 123-bus distribution system. Simulation studies show that the proposed SomRes quantification method is an appropriate indicator of system resilience, and the resilient resource allocation and dispatch strategy can significantly reduce critical load shedding under varying event propagation scenarios.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Optimizing enzymes for plastic upcycling using machine learning design and high throughput experiments

Plastic use is ubiquitous in the modern world, and polyethylene terephthalate (PET) is one of the most abundantly produced plastics (and the most highly produced polyester), with ~65 million metric tons manufactured annually. To the consumer, PET is likely most recognizable as the plastic used to make beverage bottles. Like many plastics, traditional mechanical or chemical means of PET deconstruction and upcycling are costly and inefficient. Because of these challenges, recycled plastic is generally of lower quality and is more expensive to produce than virgin plastic derived from petroleum. Ultimately, this results in most plastic ending up as waste. We view plastic waste as an underutilized resource which, with the development of more efficient and high-quality recycling processes, could (1) generate significant economic value while (2) decreasing petroleum usage and greenhouse gas emissions, as well as (3) minimizing its negative environmental and health impacts. Biocatalytic recycling, or biomanufacturing the basic building blocks of new plastic from plastic waste, is a promising approach to plastic reuse that complements existing recycling technologies. Recently, biological enzymes capable of breaking down PET have garnered significant attention as an attractive means of dealing with the plastic problem. These enzymes are currently undergoing pilot studies for implementation in industrial-scale enzyme-based recycling. However, there are significant limitations to current enzymes, including the need to perform costly pre-processing of the plastic waste before the enzymes are able to work. Further optimization of these enzymes is necessary to make these technologies competitive, and ultimately incentivise industry-wide adoption of this biology-based green recycling technology. n this work we demonstrate a means to design and generate performant biological enzymes, capable of efficiently deconstructing plastic waste. Specifically, we applied recent advances in artificial intelligence, machine learning, and statistical analysis to design new versions and discover natural enzymes capable of breaking down PET. We focused on optimizing key properties that are important for industrial-scale enzymatic recycling such as pH and thermotolerance. Normal testing of enzymatic plastic-deconstruction is extremely labor intensive and so through this work we also developed a robotic-assisted experimental pipeline capable of characterizing thousands of candidate enzymes. The results of this iterative, AI-guided, multi-discipline approach have led to increases in enzymatic breakdown of over 150X over starting enzymes. This work supports the rapidly developing and transformative field of biocatalytic solutions to environmental problems beyond the discovery and predictive understanding of enzymes for polymer recycling, and has wide implications for tackling numerous energy problems such as carbon capture and fixation (e.g., engineering carbon monoxide dehydrogenase and the rubisco-pathway), biomining (e.g., design of lanthanide-binding proteins) and biomanufacturing (e.g., lignin-deconstruction enzymes).

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Simulating Impacts of Extreme Events on Grids with High Penetrations of Wind Power Resources

As extreme weather events become more frequent and intense, the demand for connecting grid operation and infrastructure planning with extreme event models will increase as well. We present a methodology for creating damage contingencies and scenarios for electric transmission grids during a hurricane strike. Included is an example case study: Hurricane Dolly damaging a synthetic 2000 bus test system during its landing in Southern Texas. Using WIND Toolkit meteorological data in conjunction with fragility curves for various electric grid elements, we generate stochastic damage scenarios that can be used for short- and long-term planning problems including emergency asset management. We perform statistical analysis of damages and quantify topological effects on example synthetic grid. Also, we investigate loss-of-load events during two-day economic dispatch experiment. Finally, we point out various shortcomings of our method and suggest how it can be improved.

data-driven forecasting↗
DOE OSTI2021

Considering uncertainties expands the lower tail of maize yield projections

Crop yields are sensitive to extreme weather events. Improving the understanding of the mechanisms and the drivers of the projection uncertainties can help to improve decisions. Previous studies have provided important insights, but often sample only a small subset of potentially important uncertainties. Here we expand on a previous statistical modeling approach by refining the analyses of two uncertainty sources. Specifically, we assess the effects of uncertainties surrounding crop-yield model parameters and climate forcings on projected crop yield. We focus on maize yield projections in the eastern U.S.in this century. We quantify how considering more uncertainties expands the lower tail of yield projections. We characterized the relative importance of each uncertainty source and show that the uncertainty surrounding yield model parameters is the main driver of yield projection uncertainty.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Data Mining and Machine Learning for Power System Monitoring, Understanding, and Impact Evaluation

This chapter presents results from the Big Data analysis framework to improve power system situational awareness and system reliability. For this purpose, a dataset with real-world phasor measurement unit data and historical transmission system outage data has been created and used to carry out the analysis. Several statistical analysis and machine learning methods have been developed and implemented for event and anomaly detection and modeling. Detection and analysis results for actual examples of power system events are presented. Finally, data-driven characterization and risk assessment methods for weather-related extremes in power systems are developed and demonstrated on the Bonneville Power Administration system. These applications demonstrate the capability of Machine Learning (ML) methods to monitor system abnormalities, to predict system events, and to characterize the impact of extreme events on power grid

data mining, power grid, machine learning, anomaly↗
DOE OSTI2022

Chemomechanical weakening of muscovite quantified with in situ liquid nanoindentation

Abstract Chemomechanical weakening of layered phyllosilicate muscovite mica was studied as a function of chemical environment via in situ liquid-nanoindentation under four conditions (dry, deionized water, and two NaCl solutions of different pH). While traditional indentation analyses of layered materials with extreme mechanical anisotropy have been limited due to displacement bursts (pop-ins), here the bursts were used as proxies for delamination, fracture, and spalling events. Since displacement bursts during an indent represent a stochastic process, 120 indents were conducted for each condition to generate statistically significant amounts of data. In total, over 9000 bursts were assessed using a load–displacement threshold criterion, classifying this as a high-throughput nanoscale fracture technique. For each burst, initiation load, initiation displacement, plastic zone volume at initiation, and energy dissipation were analyzed. A power-law relationship between the burst load and burst energy was noted which separated the bursts into two continuous distributions: (1) bursts due only to the mechanics of the indent and (2) bursts due to both the mechanics of the indent and the environment. By using a cumulative probability distribution, it was found that the NaCl solutions decreased the minimum plastic zone volume necessary to initiate a displacement burst by an order of magnitude relative to the dry condition. Finally, the underlying mechanisms explaining the trends in initiation volume as a function of environment were discussed, with a focus on the chemomechanical degradation processes via chemical attack and cation exchange.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Enhancing scalability of a matrix-free eigensolver for studying many-body localization

We propose several techniques to enhance the parallel scalability of a matrix-free eigensolver designed for studying many-body localization (MBL) of quantum spin chain models with nearest-neighbor interactions and on-site disorder. This type of problem is computationally challenging because the dimension of the associated Hamiltonian matrix grows exponentially with respect to the number of spins L, and we need to average over different realizations of the random disorder to obtain relevant statistical behavior. For each disorder realization, we need to compute eigenvalues from different regions of the spectrum and their corresponding eigenvectors. In previous work, the interior eigenstates for a single eigenvalue problem are computed via the shift-and-invert Lanczos algorithm. Due to the extremely high memory footprint of the LU factorizations, this technique is not well suited for large L’s. For example, we need thousands of compute nodes on modern high performance computing infrastructures to go beyond L = 24. The matrix-free approach does not suffer from this memory bottleneck, however, its scalability is limited by a computation and communication load imbalance. To reduce this imbalance and to significantly enhance the scalability of the matrix-free eigensolver, we reorder the matrix and leverage the consistent space runtime, CSPACER. We also show its efficiency in managing irregular communication patterns at scale compared to optimized MPI non-blocking two-sided and one-sided RMA implementation variants. This effort enables us to study MBL for spin chains with a larger number of spins. The efficiency and effectiveness of the proposed algorithm is demonstrated by computing eigenstates on a massively parallel many-core high performance computer.

METIS↗
DOE OSTI2020

The ongoing need for high-resolution regional climate models: Process understanding and stakeholder information

Regional climate modeling addresses our need to understand and simulate climatic processes and phenomena unresolved in global models. High resolution models are generally more skillful in simulating extremes, such as heavy precipitation, strong winds, and severe storms. In addition, research has shown that fine-scale features such as mountains, coastlines, lakes, irrigation, land use, and urban heat islands can substantially influence a region’s climate and its response to changing forcings. Regional climate simulations explicitly simulating convection are now being performed, providing an opportunity to illuminate new physical behavior that previously was represented by parameterizations with large uncertainties. Regional and global models are both advancing toward higher resolution, as computational capacity increases. However, the resolution and ensemble size necessary to produce a sufficient statistical sample of these processes in global models has proven too costly for contemporary supercomputing systems. Regional climate models are thus indispensable tools that complement global models for understanding regional climate variability and change, and are critical for supporting societal responses to changing climate.

Gutowski, William↗
DOE OSTI2023

Modeling aqueous association constants and mineral solubilities at subcritical and supercritical temperatures

The need for sustainable power generation has increased interest in the use of hydrothermal fluids for industrial applications. New high-enthalpy geothermal systems and biowaste-to-fuel processes are two relevant examples that employ supercritical fluids which require an in-depth understanding of complex chemical reactions occurring near the supercritical temperature of water (374 °C). As these processes operate in thermodynamic regimes that are not currently covered by a standard molar Gibbs energy of formation model, only empirical fits for single reaction systems are available which limit the use of multi-component phase equilibria calculations that are standard practice for less extreme environments. Here, in this paper, we advance a standard molar Gibbs energy of formation model able to operate in these otherwise inaccessible thermodynamic states to include species needed for key mineral solubility systems and ion association reactions. This work extends a model based on molecular statistical thermodynamics (MST) into four new systems (Na 3 PO 4 -H 2 O, LiOH-H 2 O, KOH-H 2 O, and BaSO 4 -H 2 O) by extending the model to cover 10 new species. For each of these systems, model predictions were consistently within the experimental uncertainties for the new systems covered. A breakdown of MST contributions to the model revealed that electrostatic and hard sphere contributions were key to reproducing density dependencies of standard molar Gibbs energy of formation values around the critical point of water.

42 ENGINEERING↗
DOE OSTI2020

Genome-Wide Analysis of RNA Decay in the Cyanobacterium Synechococcus sp. Strain PCC 7002

ABSTRACT RNA degradation is an important process that influences the ultimate concentration of individual proteins inside cells. While the main enzymes that facilitate this process have been identified, global maps of RNA turnover are available for only a few species. Even in these cases, there are few sequence elements that are known to enhance or destabilize a native transcript; even fewer confer the same effect when added to a heterologous transcript. To address this knowledge gap, we assayed genome-wide RNA degradation in the cyanobacterium Synechococcus sp. strain PCC 7002 by collecting total RNA samples after stopping nascent transcription with rifampin. We quantified the abundance of each position in the transcriptome as a function of time using RNA-sequencing data and later analyzed the global mRNA decay map using machine learning principles. Half-lives, calculated on a per-ORF (open reading frame) basis, were extremely short, with a median half-life of only 0.97 min. Despite extremely rapid turnover of most mRNA, transcripts encoding proteins involved in photosynthesis were both highly expressed and highly stable. Upon inspection of these stable transcripts, we identified an enriched motif in the 3′ untranslated region (UTR) that had similarity to Rho-independent terminators. We built statistical models for half-life prediction and used them to systematically identify sequence motifs in both 5′ and 3′ UTRs that correlate with stabilized transcripts. We found that transcripts linked to a terminator containing a poly(U) tract had a longer half-life than both those without a poly(U) tract and those without a terminator. IMPORTANCE RNA degradation is an important process that affects the final concentration of individual mRNAs, affecting protein expression and cellular physiology. Studies of how RNA is degraded increase our knowledge of this fundamental process as well as enable the creation of genetic tools to manipulate RNA stability. By studying global transcript turnover, we searched for sequence elements that correlated with transcript (in)stability and used these sequences to guide tool design. This study probes global RNA turnover in a cyanobacterium, Synechococcus sp. strain PCC 7002, that both has a unique array of RNases that facilitate RNA degradation and is an industrially relevant strain that could be used to convert CO 2 and sunlight into useful products.

59 BASIC BIOLOGICAL SCIENCES↗