Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parameter sampling”

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

A Bayesian Approach for Quantifying Data Scarcity when Modeling Human Behavior via Inverse Reinforcement Learning

Computational models that formalize complex human behaviors enable study and understanding of such behaviors. However, collecting behavior data required to estimate the parameters of such models is often tedious and resource intensive. Thus, estimating dataset size as part of data collection planning (also known as Sample Size Determination) is important to reduce the time and effort of behavior data collection while maintaining an accurate estimate of model parameters. In this paper, we present a sample size determination method based on Uncertainty Quantification (UQ) for a specific Inverse Reinforcement Learning (IRL) model of human behavior, in two cases: 1) pre-hoc experiment design—conducted in the planning stage before any data is collected, to guide the estimation of how many samples to collect; and 2) post-hoc dataset analysis—performed after data is collected, to decide if the existing dataset has sufficient samples and whether more data is needed. Here, we validate our approach in experiments with a realistic model of behaviors of people with Multiple Sclerosis (MS) and illustrate how to pick a reasonable sample size target. Our work enables model designers to perform a deeper, principled investigation of effects of dataset size on IRL.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Photometric redshift uncertainties in weak gravitational lensing shear analysis: models and marginalization

ABSTRACT Recovering credible cosmological parameter constraints in a weak lensing shear analysis requires an accurate model that can be used to marginalize over nuisance parameters describing potential sources of systematic uncertainty, such as the uncertainties on the sample redshift distribution n(z). Due to the challenge of running Markov chain Monte Carlo (MCMC) in the high-dimensional parameter spaces in which the n(z) uncertainties may be parametrized, it is common practice to simplify the n(z) parametrization or combine MCMC chains that each have a fixed n(z) resampled from the n(z) uncertainties. In this work, we propose a statistically principled Bayesian resampling approach for marginalizing over the n(z) uncertainty using multiple MCMC chains. We self-consistently compare the new method to existing ones from the literature in the context of a forecasted cosmic shear analysis for the HSC three-year shape catalogue, and find that these methods recover statistically consistent error bars for the cosmological parameter constraints for predicted HSC three-year analysis, implying that using the most computationally efficient of the approaches is appropriate. However, we find that for data sets with the constraining power of the full HSC survey data set (and, by implication, those upcoming surveys with even tighter constraints), the choice of method for marginalizing over n(z) uncertainty among the several methods from the literature may modify the 1σ uncertainties on Ωm–S8 constraints by ∼4 per cent, and a careful model selection is needed to ensure credible parameter intervals.

Zhang, Tianqing (ORCID:000000025596198X)↗
DOE OSTI2026

Quantum-Inspired Bayesian Sampling for Uncertainty Quantification and Machine Learning (Final Technical Report)

With increasing simulation and measurement data, machine learning and artificial intelligence have been widely used in computational decision-making of complex engineering systems. The resulting tools, such as uncertainty quantification solvers, reinforcement learning, and physics-informed machine learning, have achieved great success in critical DOE tasks such as material discovery and design, energy system modeling and control, and numerical weather and climate prediction. A core topic in scientific machine learning and artificial intelligence is Bayesian inference: given an observed data set, people want to estimate the posterior distribution of a (possibly large) number of hidden parameters. Due to the flexibility and weak assumptions, Bayesian sampling has been the mainstream Bayesian inference solvers despite the rapid progress of approximate Bayesian inference. Classical Bayesian sampling methods such as Markov-chain Monte Carlo suffer from a low-acceptance rate due to the random walk nature, therefore state-of-the-art techniques use Hamiltonian Monte Carlo and its variants to efficiently draw posterior samples in a high dimension. The key idea of Hamiltonian Monte Carlo and its variants is to simulate the Hamiltonian dynamics of a classical particle with a fixed mass, and their performance significantly degrades when the posterior distribution is highly spiky or has multiple modes. Leveraging the idea of quantum physics, this project has investigated new theory, algorithms and applications of Bayesian inference (especially Bayesian sampling). The main results include: (1) novel quantum-inspired Bayesian sampling methods that can lead to better accuracy for challenging multi-modal or spiky distributions, (2) more scalable machine learning framework leveraging tensor-compressed Bayesian inference, and (3) Bayesian and sampling approaches for verifying the robustness of continuous and binary neural networks.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A forward modeling approach to analyzing galaxy clustering with S IM BIG

We present cosmological constraints from a simulation-based inference (SBI) analysis of galaxy clustering from the SimBIG forward modeling framework. SimBIG leverages the predictive power of high-fidelity simulations and provides an inference framework that can extract cosmological information on small nonlinear scales. In this work, we apply SimBIG to the Baryon Oscillation Spectroscopic Survey (BOSS) CMASS galaxy sample and analyze the power spectrum, P ℓ (k), to k max = 0.5 h/Mpc. We construct 20,000 simulated galaxy samples using our forward model, which is based on 2,000 high-resolution Quijote N -body simulations and includes detailed survey realism for a more complete treatment of observational systematics. We then conduct SBI by training normalizing flows using the simulated samples and infer the posterior distribution of ΛCDM cosmological parameters: Ω m , Ω b , h, n s , σ 8 . We derive significant constraints on Ω m and σ 8 , which are consistent with previous works. Our constraint on σ 8 is 27% more precise than standard P ℓ analyses because we exploit additional cosmological information on nonlinear scales beyond the limit of current analytic models, k > 0.25 h/Mpc. This improvement is equivalent to the statistical gain expected from a standard P ℓ analysis of galaxy sample ~ 60% larger than CMASS. While we focus on P ℓ in this work for validation and comparison to the literature, SimBIG provides a framework for analyzing galaxy clustering using any summary statistic. We expect further improvements on cosmological constraints from subsequent SimBIG analyses of summary statistics beyond P ℓ .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Microstructural evolution, defect mitigation, and precipitation behavior in AA6061 via laser powder bed fusion with high-temperature substrate heating

Defects, particularly solidification cracking, remain persistent challenges in the laser powder bed fusion (LPBF) processing of AA6061 aluminum alloy. This study systematically investigates defect mitigation, microstructural evolution, and mechanical properties associated with high-temperature substrate preheating at 500 °C. Comprehensive microstructural analyses, including characterization of defects, grain structures, and precipitation behavior, were performed on samples in both as-built and T6 heat-treated states. Elevated preheating substantially reduced solidification cracking across a wide processing window, while demonstrating decreased crack sensitivity to laser parameters. Columnar cracks along the build direction were observed despite substrate preheating. Lack-of-fusion and keyhole porosity were effectively eliminated, though gas-induced microporosity persisted at higher powers. In-depth characterization of two distinct laser power and speed conditions confirmed the formation of micron-sized, non-coherent Mg 2 Si precipitates under heated substrate conditions, alongside α-AlFeCrMnSi intermetallic phases indicative of in-situ thermal effects during fabrication. Subsequent T6 heat treatment revealed the formation of fine, coherent needle-shaped β″ strengthening precipitates. Despite substantial differences in processing parameters, comparable mechanical properties were measured in the as-built samples (∼52 MPa yield strength, ∼130 MPa tensile strength), primarily due to reduced strain hardening effects and consistent precipitation characteristics. Meanwhile, the T6 heat treatment led to significant improvement in properties, enhancing yield strength by over 400%, aligning closely with the performance of conventional wrought AA6061-T6. These findings underscore that high-temperature substrate preheating offers an effective means to suppress cracks, control precipitation, and enhance mechanical performance in LPBF-processed AA6061.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Diffusion kinetics of 3 He in pyroxene and plagioclase and applications to cosmogenic exposure dating and paleothermometry in mafic rocks

In this study, we investigate the diffusivity of cosmogenic 3 He in a variety of plagioclase and pyroxene compositions, and its application to paleothermometry and exposure dating in these minerals, through stepwise degassing experiments. While cosmogenic 3 He has been utilized for exposure dating in pyroxene for decades due to its retentivity, plagioclase, often found along with pyroxene in mafic rocks, is generally less retentive of cosmogenic noble gas. However, the diffusivity of 3 He in either plagioclase or pyroxene has not yet been measured quantitatively. A challenge in measuring diffusion kinetics by step-degassing experiments in poorly retentive minerals is the fact that significant amounts of He can be lost prior to the experiment. To address this issue, we apply a forward “multiple diffusion domain” (MDD) inversion model that includes model predictions of initial gas loss during irradiation and storage of the samples to account for this observation and add constraints to the diffusion parameters. We find that 3 He diffusivity in plagioclase appears to be highly variable. This variability can be explained by the MDD inversion models' inability to constrain the diffusion parameters when significant gas has been lost during irradiation and/or prolonged storage prior to experiment analysis, resulting in an overestimation of 3 He retentivity. Plagioclase samples that were kept frozen after irradiation, to limit this initial gas loss, yielded the most reliable estimate of diffusion kinetics. We find that 3 He in plagioclase is diffusively lost at Earth's surface temperatures on a timescale of 100 years and is therefore unsuitable for surface temperature paleothermometry. On the contrary, we find cosmogenic 3 He in pyroxene to be retentive at Earth's surface temperatures on a 1-million-year timescale.

Bergelin, Marie [Berkeley Geochronology Center (BG↗
DOE OSTI2024

Learning new physics from data: A symmetrized approach

Thousands of person years have been invested in searches for new physics (NP), the majority of them motivated by theoretical considerations. Yet, no evidence of beyond the Standard Model physics has been found. This suggests that model-agnostic searches might be an important key to explore NP, and help discover unexpected phenomena which can inspire future theoretical developments. A possible strategy for such searches is identifying asymmetries between data samples that are expected to be symmetric within the Standard Model. We propose exploiting neural networks (NNs) to quickly fit and statistically test the differences between two samples. Our method is based on an earlier work, originally designed for inferring the deviations of an observed dataset from that of a much larger reference dataset. We present a symmetric formalism, generalizing the original one, avoiding fine-tuning of the NN parameters and any constraints on the relative sizes of the samples. Our formalism could be used to detect small symmetry violations, extending the discovery potential of current and future particle physics experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Lens Modeling of STRIDES Strongly Lensed Quasars Using Neural Posterior Estimation

Strongly lensed quasars can be used to constrain cosmological parameters through time-delay cosmography. Models of the lens masses are a necessary component of this analysis. To enable time-delay cosmography from a sample of $\mathcal{O}(10^3)$ lenses, which will soon become available from surveys like the Rubin Observatory’s Legacy Survey of Space and Time and the Euclid Wide Survey, we require fast and standardizable modeling techniques. To address this need, we apply neural posterior estimation (NPE) for modeling galaxy-scale strongly lensed quasars from the Strong Lensing Insights into the Dark Energy Survey (STRIDES) sample. NPE brings two advantages: speed and the ability to implicitly marginalize over nuisance parameters. We extend this method by employing sequential NPE to increase precision of mass model posteriors. We then fold individual lens models into a hierarchical Bayesian inference to recover the population distribution of lens mass parameters, accounting for out-of-distribution shift. After verifying our method using simulated analogs of the STRIDES lens sample, we apply our method to 14 Hubble Space Telescope single-filter observations. We find the population mean of the power-law elliptical mass distribution slope, γ lens , to be $\mathcal{M}_γ$ lens = 2.13 ± 0.06. Our result represents the first population-level constraint for these systems. This population-level inference from fully automated modeling is an important stepping stone toward cosmological inference with large samples of strongly lensed quasars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Combined full shape analysis of BOSS galaxies and eBOSS quasars using an iterative emulator

ABSTRACT Standard full-shape clustering analyses in Fourier space rely on a fixed power spectrum template, defined at the fiducial cosmology used to convert redshifts into distances, and compress the cosmological information into the Alcock–Paczynski parameters and the linear growth rate of structure. In this paper, we propose an analysis method that operates directly in the cosmology parameter space and varies the power spectrum template accordingly at each tested point. Predictions for the power spectrum multipoles from the TNS model are computed at different cosmologies in the framework of $\Lambda \rm {CDM}$. Applied to the final eBOSS QSO and LRG samples together with the low-z DR12 BOSS galaxy sample, our analysis results in a set of constraints on the cosmological parameters Ωcdm, H0, σ8, Ωb, and ns. To reduce the number of computed models, we construct an iterative process to sample the likelihood surface, where each iteration consists of a Gaussian process regression. This method is validated with mocks from N-body simulations. From the combined analysis of the (e)BOSS data, we obtain the following constraints: σ8 = 0.877 ± 0.049 and $\Omega _{\rm m}=0.304^{+0.016}_{-0.010}$ without any external prior. The eBOSS quasar sample alone shows a 3.1σ discrepancy compared to the Planck prediction.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Deducing neutron star equation of state from telescope spectra with machine-learning-derived likelihoods

The interiors of neutron stars reach densities and temperatures beyond the limits of terrestrial experiments, providing vital laboratories for probing nuclear physics. While the star's interior is not directly observable, its pressure and density determine the star's macroscopic structure which affects the spectra observed in telescopes. The relationship between the observations and the internal state is complex and partially intractable, presenting difficulties for inference. Previous work has focused on the regression from stellar spectra of parameters describing the internal state. We demonstrate a calculation of the full likelihood of the internal state parameters given observations, accomplished by replacing intractable elements with machine learning models trained on samples of simulated stars. Our machine-learning-derived likelihood allows us to perform maximum a posteriori estimation of the parameters of interest, as well as full scans. We demonstrate the technique by inferring stellar mass and radius from an individual stellar spectrum, as well as equation of state parameters from a set of spectra. Our results are more precise than pure regression models, reducing the width of the parameter residuals by 11.8% in the most realistic scenario. The neural networks will be released as a tool for fast simulation of neutron star properties and observed spectra.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Joint iterative reconstruction and 3D rigid alignment for X-ray tomography

X-ray tomography is widely used for three-dimensional structure determination in many areas of science, from the millimeter to the nanometer scale. The resolution and quality of the 3D reconstruction is limited by the availability of alignment parameters that correct for the mechanical shifts of the sample or sample stage for the images that constitute a scan. In this paper we describe an algorithm for marker-free, fully automated and accurately aligned and reconstructed X-ray tomography data. Our approach solves the tomographic reconstruction jointly with projection data alignment based on a rigid-body deformation model. We demonstrate the robustness of our method on both synthetic phantom and experimental data and show that our method is highly efficient in recovering relatively large alignment errors without prior knowledge of a low resolution approximation of the 3D structure or a reasonable estimate of alignment parameters.

36 MATERIALS SCIENCE↗
DOE OSTI2020

NNSS Soils Monitoring: Plutonium Valley (CAU 366) FY2019

Desert Research Institute (DRI) conducted a field assessment of the potential for radionuclide-contaminated soil to be transported from the Plutonium Valley Dispersion Sites Contamination Area (CA) because of both wind and storm water runoff. This activity supported U.S. Department of Energy (DOE) Environmental Management Nevada Program (EM NV) efforts to establish post-closure monitoring plans for the Plutonium Valley Dispersion Sites Corrective Action Unit (CAU) 366. The DRI task was intended to identify the likely mechanism(s) of transport and determine the meteorological conditions that might cause the movement of radionuclide-contaminated soils. The emphasis of the work was on collecting sediment transported by channelized storm runoff and measuring airborne dust concentrations and associated wind conditions. These data will facilitate an appropriate post-closure monitoring program. In 2011, DRI installed meteorological monitoring stations north and south of the Plutonium Valley CA—known as Pu Valley North and Pu Valley South, respectively—as well as a fluvial sediment sampling station within the CA. Since installation, temperature, wind speed, wind direction, relative humidity, precipitation, solar radiation, barometric pressure, soil temperature, volumetric soil moisture content, and airborne particulate concentrations have been collected at both meteorological stations. The maximum, minimum, and average or total (as appropriate) for each of these parameters were recorded for each 10-minute interval. The sediment sampling station included an automatically activated sampling pump with collection bottles for suspended sediment, which was activated when sufficient flow was present in the channel, and passive traps for bedload material that was transported down the channel during runoff events. This report presents data collected from these stations during fiscal year (FY) 2019. During the FY2019 (October 1, 2018, through September 30, 2019) reporting period, the warmest month was July and the coldest month was February. Monthly total solar radiation was highest in June or July (depending on the station) and lowest in December. Monthly mean wind speeds were highest during February. At the Pu Valley North station, winds were commonly northerly, northwesterly, and southerly. At the Pu Valley South station, the wind direction was more variable and winds frequently came from most directions, excluding westerly and northwesterly winds. Winds above 15.0 miles per hour (mph) (24.1 kilometers per hour [km/hr]) were frequently southerly at both stations, and less-frequent strong winds were northerly (at the Pu Valley South station) or northwesterly (at the Pu Valley North station). Monthly average relative humidity ranged from 18 percent in August to 62 percent in January and February. Monthly total precipitation ranged from zero at both stations in September to 2.23 inches (in) (56.6 millimeters [mm]) during February at the Pu Valley South station and 2.09 in (53.1 mm) during March at the Pu Valley North station. From October 1, 2018, through September 30, 2019, the total precipitation was 8.23 in (209 mm) and 8.15 in (207 mm) at the Pu Valley South and Pu Valley North stations, respectively. The largest daily precipitation totals of 1.46 in (37.1 mm) and 1.14 in (29.0 mm) occurred at the Pu Valley South and Pu Valley North stations, respectively, on March 6, 2019. However, the ISCO autosampler was not activated during any storm event in FY2019 because the channel water depth measurements were insufficient for suspended sediment sample collection. Additionally, because no significant flow through the channel was indicated, no bedload samples were collected. Light breezes of 0.0 mph to 5.0 mph (0.0 km/hr to 8.0 km/hr) occurred most frequently (approximately 55 percent to 57 percent of the time). The frequency of occurrence diminished exponentially as the wind speed increased such that winds in excess of 25.0 mph (40.2 km/hr) occurred less than 0.01 percent of the time. Winds in excess of 15.0 mph (24.1 km/hr) were most commonly southerly, which is likely controlled by the topography of the valley because topographic highs define the east and west sides of the valley and converge toward the north.The concentrations of PM2.5 (particulate matter with an aerodynamic diameter ≤2.5 micrometers [μm]) and PM10 (particulate matter with an aerodynamic diameter ≤10 μm) in the air generally increase as wind speed increases. However, the particle profiler at the Pu Valley South station regularly malfunctioned, which resulted in poor data quality and a limited period of record within FY2019. Therefore, discussion of the annual dust observations is mostly limited to the Pu Valley North station. At the Pu Valley North station, significant increases in windblown dust concentrations were observed when wind speeds exceeded 15.0 mph (24.1 km/hr), especially during southerly winds. When all wind speeds were considered, high dust concentrations at the Pu Valley North station were frequently associated with northerly, northwesterly, or southerly winds. The majority of PM10 transport at the Pu Valley North station occurred when winds were blowing northerly and northwesterly because of the high frequency of winds up to 15.0 mph (24.1 km/hr). The ratio of PM10 to PM2.5 is a qualitative indicator of the proximity of dust sources to the observation point and the values of the PM10 to PM2.5 ratio tend to be higher nearer to the source. The PM10 to PM2.5 ratio increased from 1.2 to 2.8 with higher wind speeds at the Pu Valley North station. This increase suggests that some locally sourced PM10 was present in the air when local winds exceeded 15.0 mph (24.1 km/hr). Local dust suspension is interpreted during major dust events when peaks in the PM10 to PM2.5 ratio occur simultaneously with the maximum wind speed and the peak PM10 concentration. Analysis of the five major dust events in FY2019 revealed that northerly winds were dominant in three events and southerly winds were dominant in two events. Local suspension of dust within Plutonium Valley is suspected to have occurred during two of the five major events (April 9, 2019, and September 28 to 29, 2019), and a nonlocal source of dust is presumed for the three other major dust events (October 14, 2018; May 16, 2019; and September 2, 2019). The three suspected regional dust transport events exhibited winds that typically did not exceed the 15.0 mph (24.1 km/hr) wind speed threshold for local dust transport and featured PM10 to PM2.5 ratios that indicate large-scale regional dust transport. Similar event characteristics (e.g., wind speed and PM10 concentration) were observed at air monitoring stations approximately 70 mi (113 km) to the northwest within the Tonopah Test Range (TTR). Therefore, these three regional events may not reflect transport from the Plutonium Valley CA. Routine monitoring of meteorological and hydrologic parameters at Plutonium Valley has been discontinued at the request of EM NV, effective at the conclusion of FY2019. Environmental data acquisition has been completed, and both Pu Valley North and Pu Valley South meteorological stations within Plutonium Valley as well as the ISCO autosampler station within the CA were removed on October 14, 2019.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Simulation Study of High-Precision Characterization of MeV Electron Interactions for Advanced Nano-Imaging of Thick Biological Samples and Microchips

The resolution of a mega-electron-volt scanning transmission electron microscope (MeV-STEM) is primarily governed by the properties of the incident electron beam and angular broadening effects that occur within thick biological samples and microchips. A precise understanding and mitigation of these constraints require detailed knowledge of beam emittance, aberrations in the STEM column optics, and energy-dependent elastic and inelastic critical angles of the materials being examined. This simulation study proposes a standardized experimental framework for comprehensively assessing beam intensity, divergence, and size at the sample exit. This framework aims to characterize electron-sample interactions, reconcile discrepancies among analytical models, and validate Monte Carlo (MC) simulations for enhanced predictive accuracy. Our numerical findings demonstrate that precise measurements of these parameters, especially angular broadening, are not only feasible but also essential for optimizing imaging resolution in thick biological samples and microchips. By utilizing an electron source with minimal emittance and tailored beam characteristics, along with amorphous ice and silicon samples as biological proxies and microchip materials, this research seeks to optimize electron beam energy by focusing on parameters to improve the resolution in MeV-STEM/TEM. This optimization is particularly crucial for in situ imaging of thick biological samples and for examining microchip defects with nanometer resolutions. Our ultimate goal is to develop a comprehensive mapping of the minimum electron energy required to achieve a nanoscale resolution, taking into account variations in sample thickness, composition, and imaging mode.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Is the Local Ion Density Sufficient to Drive NaCl Nucleation from the Melt and Aqueous Solution?

Even though nucleation is ubiquitous in different science and engineering problems, investigating nucleation is extremely difficult due to the complicated ranges of time and length scales involved. In this work, we simulate NaCl nucleation in both molten and aqueous environments using enhanced sampling of all-atom molecular dynamics with deep-learning-based estimation of reaction coordinates. By incorporating various structural order parameters and learning the reaction coordinate as a function thereof, we achieve significantly improved sampling relative to traditional ad hoc descriptions of what drives nucleation, particularly in an aqueous medium. Our results reveal a one-step nucleation mechanism in both environments, with reaction coordinate analysis highlighting the importance of local ion density in distinguishing solid and liquid states. However, although fluctuations in the local ion density are necessary to drive nucleation, they are not sufficient. Our analysis shows that near the transition states, descriptors such as enthalpy and local structure become crucial. Furthermore, our protocol proposed here enables robust nucleation analysis and phase sampling and could offer insights into nucleation mechanisms for generic small molecules in different environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Application of sensitivity analysis in DYMOND/Dakota to fuel cycle transition scenarios

The ability to perform sensitivity analysis has been enabled for the nuclear fuel cycle simulator DYMOND through its coupling with the design and analysis toolkit Dakota. To test and demonstrate these new capabilities, a transition scenario and multi-parameter study were devised. The transition scenario represents a partial transition from the US nuclear fleet to a closed fuel cycle with small modular LWRs and fast reactors fueled by reprocessed used nuclear fuel. Four uncertain parameters in this transition were studied – start date of reprocessing, total reprocessing capacity, the nuclear energy demand growth, and the rate at which the fast reactors are deployed – with respect to their impact on four response metrics. The responses – total natural uranium consumed, maximum annual enrichment capacity required, total disposed mass, and total cost of the nuclear fuel cycle – were chosen based on measures known to be of interest in transition scenarios and to be significantly impacted by the varying parameters. Furthermore, analysis of this study was performed both from the direct sampling and through surrogate models developed in Dakota to calculate the global sensitivity measures Sobol’ indices. This example application of this new capability showed that the most consequential parameter to most metrics was the share of new build capacity that is fast reactors. However, for the cost metric, the scaling factor of the energy demand growth was significant and had synergistic behavior with the fast reactor new build share.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Efficient Distance-based Global Sensitivity Analysis for Terrestrial Ecosystem Modeling

Sensitivity analysis in terrestrial ecosystem modeling is important for understanding controlling processes, guiding model development, and targeting new observations to reduce parameter and prediction uncertainty. Complex and computationally expensive terrestrial ecosystem models (TEM) limit the number of ensemble simulations, requiring sophisticated and efficient methods to analyze sensitivities of multiple model responses to different types of parameter uncertainties. In this study, we propose a distance-based global sensitivity analysis (DGSA) method. DGSA first classifies model response samples into a small set of discrete classes and then calculates the distance between parameter frequency distributions in different classes to measure the parameter sensitivity. The principle is that, if the parameter distribution is the same in each class, then the model response is insensitive to the parameter, while a large difference in the distributions indicates the parameter is influential to the response. Built on this idea, DGSA can be applied to analyze sensitivity of a single and a group of responses to different kinds of parameter uncertainties including continuous, discrete and even stochastic. Besides the main-effect sensitivity from a single parameter, DGSA can also quantify the sensitivity from parameter interactions. Additionally, DGSA is computationally efficient which can use a small number of model evaluations to obtain an accurate and statistically significant result. We applied DGSA to two TEMs, one having eight parameters and three kinds of model responses, and the other having 47 parameters and a long-period response. We demonstrated that DGSA can be used for sensitivity problems with multiple responses and high-dimensional parameters efficiently.

Lu, Dan↗
DOE OSTI2023

Habitat severity characteristics structure soil communities at regional and local spatial scales along the Antarctica Peninsula

Antarctic soils provide an excellent setting to test biogeographical patterns across spatial and environmental scales given their relatively simple communities and the dominance of physical factors that create strong environmental gradients. Additional urgency is given by the fact that their unique terrestrial communities are the subject of conservation efforts in a rapidly changing environment. We investigated relationships of soil community assembly and alpha and beta diversity with climatic and environmental parameters across regional and local scales in Maritime Antarctica. We sampled from a regional gradient of sites that differ in habitat severity, ranging from relatively favourable to harsher physicochemical conditions. At the regional scale, bacterial community characteristics and microarthropod abundance varied along this severity gradient, but most measures of fungal communities did not. Microarthropod and microbial communities differed in which soil and climate parameters were most influential, and the specific parameters that influenced each taxon differed across broad and fine spatial scales. This suggests that conservation efforts will need to focus on a large variety of habitat characteristics to successfully encompass diversity across taxa. Because beta diversity was the result of species turnover, conservation efforts also cannot focus on only the most biodiverse sites to effectively preserve all aspects of biodiversity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Interplay between magnetism and superconductivity in UTe 2

Time-reversal symmetry breaking (TRSB) in UTe 2 was inferred from observations of a spontaneous Kerr response in the superconducting state after cooling in zero magnetic field, while a finite c-axis magnetic field training was further used to determine the nature of the nonunitary composite order parameter of this material. Here, we present an extensive study of the magnetic-field-trained Kerr effect, which unveils a unique critical state of pinned “ferromagnetic vortices.” We show that a remanent Kerr signal appears following the removal of a training magnetic field, which reflects the response to the TRSB order parameter and the external magnetic field through the paramagnetic susceptibility. Furthermore, this unambiguously demonstrates the importance of the magnetic fluctuations and their intimate relation to the composite order parameter. Focusing the beam on the center of the sample, we are able to accurately determine the maximum field that is screened by the critical state and the respective critical current. Measurements in the presence of magnetic field show the tendency of the superconductor to produce shielding currents that oppose the increase in vortex-induced magnetization due to the diverging paramagnetic susceptibility.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗