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 271 records · Page 15

Universal construction of black hole microstates

We refine and extend a recent construction of sets of black hole microstates with semiclassical interiors that span a Hilbert space of dimension e S , where S is the black hole entropy. We elaborate on the definition and properties of microstates in statistical and black hole mechanics. The gravitational description of microstates employs matter shells in the interior of the black hole, and we argue that in the limit where the shells are very heavy, the construction acquires universal validity. To this end, we show it for very wide classes of black holes; we first extend the construction to rotating and charged black holes, including extremal and near-extremal solutions, with or without supersymmetry, and we sketch how the construction of microstates can be embedded in string theory. We then describe how the approach can include general quantum corrections, near or far from extremality. For supersymmetric black holes, the microstates we construct differ from other recent constructions in that the interior excitations are not confined within the near-extremal throat. Published by the American Physical Society 2024

Astronomy & Astrophysics↗

Atmospheric temperature, water vapour and liquid water path from two microwave radiometers during MOSAiC

The microwave radiometers HATPRO (Humidity and Temperature Profiler) and MiRAC-P (Microwave Radiometer for Arctic Clouds - Passive) continuously measured radiation emitted from the atmosphere throughout the Multidisciplinary drifting Observatory for the Study of the Arctic Climate (MOSAiC) expedition on board the research vessel Polarstern. From the measured brightness temperatures, we have retrieved atmospheric variables using statistical methods in a temporal resolution of 1 s covering October 2019 to October 2020. The integrated water vapour (IWV) is derived individually from both radiometers. In addition, we present the liquid water path (LWP), temperature and absolute humidity profiles from HATPRO. To prove the quality and to estimate uncertainty, the data sets are compared to radiosonde measurements from Polarstern. The comparison shows an extremely good agreement for IWV, with standard deviations of 0.08–0.19 kg m -2 (0.39–1.47 kg m -2 ) in dry (moist) situations. The derived profiles of temperature and humidity denote uncertainties of 0.7–1.8 K and 0.6–0.45 gm -3 in 0–2 km altitude.

54 ENVIRONMENTAL SCIENCES↗

Regression of exchangeable relational arrays

Relational arrays represent measures of association between pairs of actors, often in varied contexts or over time. Trade flows between countries, financial transactions between individuals, contact frequencies between school children in classrooms and dynamic protein-protein interactions are all examples of relational arrays. Elements of a relational array are often modelled as a linear function of observable covariates. Uncertainty estimates for regression coefficient estimators, and ideally the coefficient estimators themselves, must account for dependence between elements of the array, e.g., relations involving the same actor. Existing estimators of standard errors that recognize such relational dependence rely on estimating extremely complex, heterogeneous structure across actors. This paper develops a new class of parsimonious coefficient and standard error estimators for regressions of relational arrays. Here we leverage an exchangeability assumption to derive standard error estimators that pool information across actors, and are substantially more accurate than existing estimators in a variety of settings. This exchangeability assumption is pervasive in network and array models in the statistics literature, but not previously considered when adjusting for dependence in a regression setting with relational data. We demonstrate improvements in inference theoretically, via a simulation study, and by analysis of a dataset involving international trade.

97 MATHEMATICS AND COMPUTING↗

Collaborative: in situ visual analytics technologies for extreme scale combustion simulations

This project aims to drastically enhance the usability of in situ analysis and visualization for extreme-scale scientific simulations. Current exascale computing capabilities promise to offer greater predictive ability of simulations and to further push the frontiers of science and technology. However, to validate the simulation output at extreme scale, examine the modeled phenomena, and discover previously unknowns from the output data, the output must be reduced or transformed in situ as it is being generated during the simulation such that the amount of data to examine and store is kept to a minimum. Such in situ approaches allow us to process and analyze the data and any embedded geometry to an extent that would be prohibitively expensive, if not impossible, to perform as a post hoc task. While in situ processing has been demonstrated to be a feasible and promising approach, its full potential has not yet been leveraged. In this project, we have developed comprehensive enhancements to in situ technology based on probability distributions in data. Our research focuses on jointly developing new ways of interacting with massive statistical samples while creatively utilizing new state-of-the-art computational resources to push the boundaries of in situ exploration. Moreover, we have developed new time-dependent techniques to enable previously unattainable capabilities in areas such as intelligent simulation steering and precise feature identification. We have experimentally studied our design and implementation at NERSC and OLCF, and are able to leverage existing in situ infrastructures whenever possible. While the exemplar in this project is combustion, many other fields for which turbulent transport is important, e.g., fusion, climate, astrophysics among others, encounter similar issues as simulations scale up to the exascale. This project shows its potential to generate high impact on DOE missions since the resulting technology promises to improve scientists’ ability to rapidly and correctly interpret and tune extreme-scale simulations, leading to new scientific understanding and advancements.

97 MATHEMATICS AND COMPUTING↗

Changes in Climate and Its Effect on Timing of Snowmelt and Intensity-Duration-Frequency Curves

Snow is a critical water resource for much of the U.S. and failure to account for changes in climate could deleteriously impact military assets. In this study, we produced historical and future snow trends through modeling at three military sites (in Washington, Colorado, and North Dakota) and the Western U.S. For selected rivers, we performed seasonal trend analysis of discharge extremes. We calculated flood frequency curves and estimated the probability of occurrence of future annual maximum daily rainfall depths. Additionally, we generated intensity-duration-frequency curves (IDF) to find rainfall intensities at several return levels. Generally, our results showed a decreasing trend in historical and future snow duration, rain-on-snow events, and snowmelt runoff. This decreasing trend in snowpack could reduce water resources. A statistically significant increase in maximum streamflow for most rivers at the Washington and North Dakota sites occurred for several months of the year. In Colorado, only a few months indicated such an increase. Future IDF curves for Colorado and North Dakota indicated a slight increase in rainfall intensity whereas the Washington site had about a twofold increase. This increase in rainfall intensity could result in major flood events, demonstrating the importance of accounting for climate changes in infrastructure planning.

58 GEOSCIENCES↗

Dark Energy Survey Year 3 results: optimized $w$CDM simulation-based inference with weak lensing map-level hybrid statistics

We present cosmological constraints from the Dark Energy Survey Year 3 (DES Y3) weak lensing data using hierarchical hybrid statistics within a Bayesian simulation-based inference framework that is based on the Gower Street simulations. To maximize the precision of the inference, we have developed a new, information-theory based, data compression of the weak lensing maps to just seven highly informative summary statistics. The hybrid scheme exploits the high information content of the power spectrum, compressing both the power spectrum and neural-based summaries that are designed to extract further information. Our simulation-based approach enables principled forward modelling of all major sources of systematic uncertainty and survey properties into realistic mock observations, including the survey mask, photometric redshift uncertainties, intrinsic galaxy alignments, multiplicative shear calibration bias, source galaxy clustering, non-Gaussian shape noise, and non-linear structure formation. The summary statistics are then used in a Bayesian simulation-based inference pipeline. The inference is validated through coverage tests and checks for robustness against baryonic feedback. Assuming a $w$CDM cosmology, our analysis yields $S_8 = 0.808 \pm 0.017$, $Ω_{\rm m} = 0.325 \pm 0.024$, and $w < -0.766$ (marginalized posterior 68 per cent credible intervals). This rigorous combination of information theory, physics- and neural network-based extreme data compression, and principled Bayesian analysis improves the figure of merit for $(Ω_{\rm m}, S_8, w)$ by 60 per cent over the previous state-of-the-art, and by almost a factor of 3 over two-point analyses of the same data. They are the most precise joint constraints on $(Ω_{\rm m}, S_8, w)$ from weak gravitational lensing data alone of any survey to date. We intend to apply this analysis to the more recent DES Y6 data.

Williamson, J. [University Coll. London]↗

A simple framework to characterize land aridity based on surface energy partitioning regimes

Land aridity is often characterized by the aridity index (AI), which does not account for land surface water-energy interactions that are crucially important in determining regional climate. Such interactions can be captured by the evaporative fraction (EF, ratio of evapotranspiration to available energy) regimes. As EF is subject to energy and water limitations in humid and dry areas, respectively, EF regimes may be used to characterize land aridity to account for the influence of complex land characteristics and their impact on water availability. Here, we propose a simple framework to characterize land aridity by statistically ranking the coupling strength between EF and surface energy and water terms. The framework is demonstrated using gridded data and compared with AI over the U.S. and China. Results show that regionalization of aridity zones based on EF regimes and a two-tiered classification scheme may provide information such as surface energy and water variability complementary to the background aridity depicted by AI, with implications for extreme events.

54 ENVIRONMENTAL SCIENCES↗

Heat transfer and transport property contrast effects on the compressible Rayleigh-Taylor instability

In extreme environments, the Rayleigh-Taylor instability (RTI) may occur under large variations in density and temperature and with fluid transport properties strongly dependent on temperature. Direct numerical simulations of the 3D fully compressible RTI are conducted, examining the idealized configuration of a hotter, less dense fluid pushing against a colder, denser fluid. Various temperature ratios and transport property configurations are explored to examine how heat conduction, large variations in transport properties, and sudden changes in transport properties can affect the evolution of the mixing layer. Nonuniform fluid expansion and contraction induced by heat transfer can significantly affect local density differences and overall instability growth, causing profile asymmetries about the initial interface for flow and mixing statistics. The departures from classical self-similar development of the instability along with misalignment between regions of mixing and regions of most intense turbulent activity caused by both heat transfer and transport property contrasts are examined. After sudden changes in fluid transport properties, which may occur as a result of rapid heating (e.g., in inertial confinement fusion), the flow quickly responds and begins to relax towards quasi-self-similar late-time evolution. For many dynamical quantities such as vorticity and dissipation, this late-time evolution resembles that of the configuration that already started with the final transport property magnitudes, suggesting that these quantities depend only on the transport properties and not on past flow history, provided that the density field distributions for the flows remain similar. On the other hand, the mixing evolution after the transport property change is unique, implying that both property magnitudes and previous history are impactful on the mixing. In conclusion, these simulations demonstrate how various temperature-related effects are extremely important to consider in compressible RTI flows with large temperature variations.

42 ENGINEERING↗

Time-dependent THMC properties and microstructural evolution of damaged rocks in excavation damage zone

Modeling coupled thermo-hydro-mechanical-chemical (THMC) processes in host rocks near high-level nuclear waste (HLW) repositories at various time scales is an extremely challenging task. The current study integrates experimental, theoretical, and numerical methods in assessing the evolution of excavation damage zone (EDZ) over time and its implication on the long-term migration of hazardous species. Argillite and rock salt and are the focus of this study. The first part of the report presents a novel time-dependent directional microcrack damage theory for generic brittle rocks. It features detailed statistical description of the microcracks within a damaged solid, permitting a direct upscaling of microscale processes such as crack growth kinetics, crack closure/opening, sliding friction to explain the macroscopic creep, nonlinear elasticity, shear dilation, and loading-unloading hysteresis. This provides a basic platform for describing the anisotropic mechanical and transport properties of damaged rocks during excavation and subsequent THMC loadings. The model is validated and numerically implemented to Finite Element (FE) package ABAQUS through the user-defined material (UMAT) interface and have demonstrated great potential in resolving the time-dependent and anisotropic evolution of damage in the EDZ. The second part of the report presents a multi-scale experimental effort in characterizing the thermal, hydraulic, and mechanical properties of Mancos shale and Avery Island salt. For the Mancos shale, triaxial compression tests are performed at different confining pressures and temperatures to probe its thermomechanical properties relevant to HLW repositories. The obtained stress-strain data are interpreted using the proposed directional damage theory. Post-test specimens are subjected to gas permeability tests to reveal the correlation between permeability and the degree of microcracking. At microscale, temperature-controlled nanoindentation tests are performed and found a linear correlation between fracture toughness and elastic modulus from 25°C to 300°C. For Avery Island salt, we have designed and manufactured a novel relative-humidity controlled uniaxial creep device. Long-term creep tests at low stresses (< 5 MPa) are performed at different levels of relative humidity (RH). Besides confirming the much higher creep rates as one would expect through extrapolating the high-stress creep data, the results reveal that the steady-state creep rate of rock salt is strongly dependent on the ambient RH, an aspect that is often neglected in the literature. Both behaviors can be attributed to the pressure-solution creep mechanism which dominates at low stress and high RH levels. The third part of the report explores a set of numerical strategies in modeling the THMC behavior of porous geomaterials. A fully implicit, monolithic FE solution that can flexibly interface with different material models and coupling mechanisms for THMC problems is developed and verified through the ABAQUS user-defined element (UEL) interface. The scheme is used to study the THM response of a hypothetical HLW storage site with reference to an existing in-situ heater test. Strategies for integrating the UEL and the microcrack UMAT are suggested. Another numerical endeavor of this study is to implement a higher-order asymptotic homogenization method to account for the heterogeneous porous structures. The same method is then extended to perform microstructure-informed thermo-mechanical modeling of generalized continua. The above outcomes of this project provide a strong thrust towards enhancing the fundamental understanding and modeling capability of the long-term evolution of host rocks in EDZ, thus helping achieve the design goal of 1-million-year isolation of high-level nuclear wastes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Development of A High-Resolution Dataset for Solar Resource Adequacy Studies

High-resolution, long-term solar dataset is essential for characterizing the variability of solar energy resources and for informing strategies that ensure grid reliability and resilience in grid systems with high levels of solar energy integration. We present the development of a new 4-km, hourly Earth system dataset for the contiguous United States (CONUS), using a statistical downscaling approach that integrates the National Solar Radiation Database (NSRDB) with regional Earth system model projections. The new high-resolution Earth system dataset includes key variables - GHI, DNI, DHI, surface air temperature, and wind speed - under two future scenarios. Preliminary results show a reasonable agreement with NSRDB observations, with nBias less than 1% for GHI across CONUS. The dataset is expected to support in-depth analyses of extreme weather impacts and provide input to resource adequacy for future energy systems with diverse generation sources.

14 SOLAR ENERGY↗

Contaminant Transport Modeling for Technology Evaluation and Long-Term Monitoring in the Tims Branch Testbed, SC - 20343

Studies conducted at U.S. DOE sites have shown the presence of heavy metals, radionuclides, and volatile organic compounds in surface water, groundwater, and soil as a result of nuclear activity in the Cold War era. Since the 1990's, innovative cleanup methods have been implemented in the Tims Branch watershed at Savannah River Site (SRS) to limit the contaminant flux to the stream that have reduced the contaminant concentrations to acceptable regulatory levels in the dissolved phase. A tin-based treatment which effectively eliminated all local anthropogenic mercury inputs to this ecosystem resulted in a known step function addition of inert tin oxide particles which now serve as a potential tracer for sedimentation and particle transport processes in the stream. The long-term effectiveness of this and other remediation techniques and the potential for remobilization of adsorbed contaminant in sediment during extreme hydrologic conditions however remains unclear. It is therefore important to understand not only the fate and transport of dissolved contaminants, but also the movement of sediment and the relevant interactions with dissolved contaminant. To narrow this knowledge gap, a study is being conducted using the Tims Branch watershed as a stream-scale ecosystem test-bed to identify the primary transport processes of major contaminants of concern (such as mercury, nickel and uranium) with an emphasis on interactions with sediment transport. This involves the development of a fully distributed hydrologic watershed model of the Tims Branch watershed to predict streamflow under extreme weather conditions, as well as the development of a comprehensive contaminant transport model that can properly account for coupled contaminant and sediment transport. Review of relevant research reports and peer-reviewed journals revealed that aside from advection-dispersion transport of dissolved contaminants, adsorption and desorption with suspended solids and bed sediment also play an important role in the transport of those contaminants of concern. To develop the fully distributed hydrologic watershed model, the MIKE SHE 2-dimensional (2D) land surface/3D groundwater model that simulates surface/subsurface hydrologic processes (such as overland flow, evapotranspiration, and infiltration) was coupled with a 1D streamflow model that accounts for stream water hydraulics (such as hydraulic structures, cross-sections, and network). To model contaminant transport, the MIKE 11 streamflow component was coupled with the MIKE 11 AD module that simulates solute transport through advection and dispersion, and MIKE ECO Lab module that accounts for both sediment transport and interactions with dissolved contaminant. At this stage, the development and optimization of the fully distributed hydrologic model has been completed, achieving satisfactory statistical results between observed and predicted discharge as indicated by a root mean square error (RMSE) of 0.039 cms and a Nash-Sutcliffe efficiency coefficient (NSE) of 0.764. The ongoing development of the contaminant transport model has also yielded realistic results from preliminary tests. Results from this study are a key to evaluating the effectiveness of tin (II)-based mercury treatment of wetlands at the SRS site, and are also relevant to evaluating the potential of using this type of novel remediation technology in other mercury-contaminated stream systems. Knowledge acquired from this research will also support interpretation of historical data on the trends of contaminant concentration distribution in Tims Branch, particularly considering the effect of extreme hydrological events on the stream flow and pollutant transport. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Mechanism of Stoichiometrically Governed Titanium Oxide Brownian Tree Formation on Stepped Au(111)

Previously observed formation of substoichiometric titanium oxide dendritic structures across terraces of Au(111) is computationally studied and shown to follow the classical fractal formation mechanism of diffusion-limited aggregation (DLA). Globally optimized gas-phase oxide cluster structures are sampled in a variety of landing formations on gold surfaces and shown to favor isomers driving polymerization to Brownian tree fractal structures. Mobility of Ti 3 O 5 monomers is shown to be extremely high, with diffusion barriers of 0.21 eV or less. Through bonding stabilization, polymerization of these monomers is energetically favorable and irreversible on the 111 terrace but geometrically impossible to propagate along the step edge. Simulated scanning tunneling microscopy (STM) images show strong similarity to experiment. By contrast, observation of Ti 3 O 6 aggregating as wires along step edges is explained by the affinity of oxygen to step edges and statistical arguments for aggregation entropy at the step, in addition to low barriers for monomer diffusion and polymerization.

36 MATERIALS SCIENCE↗

Code Release for “Unlocking Extreme Space Weather through Advanced Modeling of Legacy Vela Spacecraft” ER

The software being developed for this project has two mains aims. First, a Bayesian Model Calibration (BMC) procedure is being developed to calibrate a spallation model that simulates protons hitting a spacecraft orbiting earth to real data. The procedure will be developed for general data (there is no data release requested as part of this code release). Second, an inverse physics modeling task is being undertaken to map the number of resulting neutrons observed from this process to the expected number of protons that hit the model. This second task is of statistical interest; to publish on it, the code will need to be open source.

Murph, Alexander↗

Comparison of the spatial statistics of random and defined-sequence photoresist films

The resolution-line edge roughness-sensitivity tradeoff has motivated the exploration of potential improvements using defined sequence polymers and polymer-bound photoacid generators and quenchers. We characterize the internal structures of positive tone photoresist polymer films formed from defined sequence polymers and compare them with random copolymers of the same composition. We model their imaging to connect initially to developable film structures. We use a polymer packing algorithm to simulate films of diverse compositions and locations of photoacid generators and quenchers, using the composition of an ESCAP photoresist. We use a simple extreme ultraviolet exposure-deprotection algorithm to model developable image formation within them. In all cases, the spatial distribution of chemical moieties in the film for defined sequence polymers is nearly indistinguishable from random copolymers. We evaluate several exposure-deprotection scenarios and find that a defined sequence copolymer has a distinctive developable image under certain circumstances. The use of defined sequence polymers within a photoresist layer does not automatically result in improved imaging; however, they do have some characteristics different from random polymers of the same composition. Further study of these characteristics may provide a route to improved control over the nanoscale imaging process.

36 MATERIALS SCIENCE↗

Experimental investigation of longitudinal scraping of H- bunches via photo-detachment

Longitudinal emittance growth is a significant challenge in RF linacs, especially for poorly bunched beams. This stems from particles occupying outer synchrotron oscillation orbits in the LBET, causing unwanted bunch-bunch interactions and degraded beam quality. To address this, we proposed using temporally spaced laser pulses to selectively photo-detach electrons from the longitudinal head and tail regions of H- ion bunches. This approach aims to reduce particle density in extreme orbits, enhancing beam uniformity and limiting emittance growth. Our experiments employed Fermilab's 'LaserNotcher' system at the font end of the linac, delivering 1.6 MW peak power with sub-nanosecond precision. By neutralizing the first and last half-nanosecond of several H- bunches, we measured their propagation injection into the booster. Measurements of pulse width, average height, and temporal spacing over booster cycles were compared between the scraped and unscraped bunches. Statistical analysis evaluated the results’ significance, highlighting the feasibility of laser-based scraping for future linac designs to achieve higher beam energies with improved emittance control.

Landon, Parker [Boston U.]↗

Experimental investigation of longitudinal scraping of $H-$ bunches via photo-detachment

Longitudinal emittance growth is a significant challenge in RF linacs, especially for poorly bunched beams. This stems from particles occupying outer synchrotron oscillation orbits in the LBET, causing unwanted bunch-bunch interactions and degraded beam quality. To address this, we proposed using temporally spaced laser pulses to selectively photo-detach electrons from the longitudinal head and tail regions of H- ion bunches. This approach aims to reduce particle density in extreme orbits, enhancing beam uniformity and limiting emittance growth. Our experiments employed Fermilab's 'LaserNotcher' system at the font end of the linac, delivering 1.6 MW peak power with sub-nanosecond precision. By neutralizing the first and last half-nanosecond of several H- bunches, we measured their propagation injection into the booster. Measurements of pulse width, average height, and temporal spacing over booster cycles were compared between the scraped and unscraped bunches. Statistical analysis evaluated the results’ significance, highlighting the feasibility of laser-based scraping for future linac designs to achieve higher beam energies with improved emittance control.

Landon, P. [Boston U.; Fermilab]↗

The Forward Physics Facility: Sites, Experiments, and Physics Potential

The Forward Physics Facility (FPF) is a proposal to create a cavern with the space and infrastructure to support a suite of far-forward experiments at the Large Hadron Collider during the High Luminosity era. Located along the beam collision axis and shielded from the interaction point by at least 100 m of concrete and rock, the FPF will house experiments that will detect particles outside the acceptance of the existing large LHC experiments and will observe rare and exotic processes in an extremely low-background environment. In this work, we summarize the current status of plans for the FPF, including recent progress in civil engineering in identifying promising sites for the FPF; the FPF experiments currently envisioned to realize the FPF’s physics potential; and the many Standard Model and new physics topics that will be advanced by the FPF, including searches for long-lived particles, probes of dark matter and dark sectors, high-statistics studies of TeV neutrinos of all three flavors, aspects of perturbative and non-perturbative QCD, and high-energy astroparticle physics.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Seeing is Believing: Autonomous Microscopy and the Data Revolution in Materials Science

This presentation explores the transformative potential of autonomous electron microscopy and artificial intelligence (AI) in accelerating materials science discovery, particularly for energy applications and materials operating in extreme environments. We discuss pioneering self-driving laboratories at NREL designed to intelligently probe material synthesis and degradation across multiple scales, aiming to rapidly bridge the gap between atomic-level understanding and the development of high-performance, reliable materials. Utilizing advanced machine learning techniques, such as few-shot learning and multimodal analysis integrating imaging and spectroscopy, we demonstrate methods to extract actionable descriptors for material behavior, quantify complex microstructural evolution, and statistically link synthesis parameters to defect populations. This AI-driven approach promises to accelerate the creation of predictive materials tailored for specific missions, enabling faster development cycles and enhanced material assurance.

36 MATERIALS SCIENCE↗