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 253 records · Page 14

Prediction and uncertainty quantification of shale well performance using multifidelity Monte Carlo

Uncertainty quantification is an integral component of reservoir management, especially considering the inherent uncertainty in subsurface systems. While a standard practice to estimate the uncertainty, Monte Carlo (MC) simulation is computationally intense when the sampling population comprises high-fidelity simulations. Alternatively, the Multi-fidelity Monte Carlo (MFMC) simulation overcomes this computational intensity by integrating low- and high-fidelity simulations. Our goal is to minimize the number of expensive high-fidelity simulations while maintaining accuracy and using numerous fast and cheap low-fidelity simulations to efficiently sample to input parameter space of interest. We selected gas production from unconventional wells to demonstrate the potential speedups and accuracy of the MFMC approach. The model fidelity usually determines the trade-off between accuracy and efficiency. While the high-fidelity model is more accurate, the low-fidelity model is more efficient. Our high-fidelity simulation includes reservoir simulations of a hydraulically fractured well. On the other hand, our low-fidelity model comprises the parallel-plate flow model. We used differential programming to efficiently solve the 1D flow model, where automatic differentiation is used to efficiently compute the gradients. We matched the production profile of high-fidelity simulations with our low-fidelity simulations. Then, we used a support vector regression to map the high- and low-fidelity input parameters. The mapping function is essential to tune the low-dimensional parameter space of the low-fidelity model to the high-dimensional parameter space of the high-fidelity model. We found that we can use a combination of 9 high fidelity and 10,000 low fidelity simulations to efficiently and accurately simulate pressure management. This method is at least two orders of magnitude faster than only using high-fidelity simulations. Finally, from a broader perspective, MFMC could efficiently estimate the uncertainty of various systems and models, integrating low- and high-fidelity models.

04 OIL SHALES AND TAR SANDS↗

Control of magnon–photon coupling by spin torque

We demonstrate the influence of damping and field-like torques in the magnon–photon coupling process by classically integrating the generalized Landau–Lifshitz–Gilbert equation with RLC equation in which a phase correlation between dynamic magnetization and microwave current through combined Ampère and Faraday effects are considered. We show that the gap between two hybridized modes can be controlled in samples with damping parameter in the order of 10 –3 by changing the direction of the dc current density J if a certain threshold is reached. Furthermore, our results suggest that an experimental realization of the proposed magnon–photon coupling control mechanism is feasible in yttrium iron garnet/Pt hybrid structures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Exploring the potential of using L-Band InSAR for the mapping of flooded vegetation in tropical wetlands

Wetlands play a critical role in global water and carbon cycles, yet monitoring their water extent remains difficult, particularly beneath dense vegetation. SAR-based techniques such as backscatter thresholding are limited by complex scattering mechanisms, while fully polarimetric SAR (PolSAR) data capable of detecting doublebounce scattering remain scarce. To address these challenges, this study evaluates the potential of Interferometric SAR (InSAR) for mapping water surfaces beneath vegetation, termed flooded vegetation, using the Atrato floodplain in Colombia as a case study. We develop an automated workflow combining InSAR fringe detection with local phase homogeneity analysis and random sampling of processing parameters to generate probabilistic flooded vegetation maps. Applied to ALOS PALSAR-1 L-band image pairs from 2007–2011, the workflow captures seasonal fluctuations in flooded extent ranging from 500 to 1,500 km2. Compared to other L-band SAR inundation products, the InSAR-based maps identify broader flooded areas, with ~70% agreement in pairwise comparisons. Around 84% of detections align with existing wetland inventories and seasonal changes correspond with regional hydrological indicators, including terrestrial water storage anomalies and water gauge measurements. PolSAR analysis shows that InSAR complements backscatter-based methods by detecting inundation in areas with weak double-bounce signals. These findings suggest that combining InSAR with backscatter-based methods can improve detection of flooded vegetation, which is especially relevant for the upcoming NISAR mission that will offer frequent global L-band observations.

Coastal inundation↗

Constraining Bedrock Groundwater Residence Times in a Mountain System With Environmental Tracer Observations and Bayesian Uncertainty Quantification

Groundwater residence time distributions provide fundamental insights on the hydrological processes within watersheds. Yet, observations that can constrain groundwater residence times over broad timescales remain scarce in mountain catchment studies. We use environmental tracers (CFC-12, SF 6 , 3 H, and 4 He) to investigate groundwater residence times along a hillslope in the East River Watershed, Colorado, USA. We develop a Bayesian inference framework that applies a Markov-chain Monte Carlo (MCMC) approach to estimate noble gas recharge temperature, elevation, and excess-air parameters and the resulting environmental tracer concentrations. MCMC is then used to propagate the environmental tracer uncertainties to estimates of groundwater mean residence times inferred with lumped parameter models. All samples contain 3 H, CFC-12, and SF 6 in addition to terrigenic 4 He, suggesting a mixture of water characterized by modern and premodern residence times. 4He exponential mean residence times range from hundreds of years at the upslope well to thousands of years at the toe-slope well assuming average crustal production rates. We find that binary mixing residence time distributions with separate young and old mixing fractions are needed to predict the 4 He, CFC-12, SF 6 , and 3 H observations, supporting the importance of flow path mixing in this bedrock system. Our findings that the fractured bedrock hosts groundwater with a mixture of residence times ranging from decades to millennia suggest variable recharge dynamics and flow path mixing along the hillslope and highlight the importance of characterizing groundwater systems with observations that are sensitive to transport over a broad range of residence times.

54 ENVIRONMENTAL SCIENCES↗

High Mg activation in implanted GaN by high temperature and ultrahigh pressure annealing

Not provided.We demonstrate high p-type conductivity and hole concentrations >10 18 cm -3 in Mg-implanted GaN. The implantation was performed at room temperature and by post-implantation annealing at 1 GPa of N 2 and in a temperature range of 1200–1400 °C. The high pressure thermodynamically stabilized the GaN surface without the need of a capping layer. We introduce a “diffusion budget,” related to the diffusion length, as a convenient engineering parameter for comparing samples annealed at different temperatures and for different times. Although damage recovery, as measured by XRD, was achieved at relatively low diffusion budgets, these samples did not show p-type conductivity. Further analyses showed heavy compensation by the implantation-induced defects. Higher diffusion budgets resulted in a low Mg ionization energy (~115 meV) and almost complete Mg activation. For even higher diffusion budgets, we observed significant loss of Mg to the surface and a commensurate reduction in the hole conductivity. High compensation at low diffusion budgets and loss of Mg at high diffusion budgets present a unique challenge for shallow implants. A direct control of the formation of compensating defects arising from the implantation damage may be necessary to achieve both hole conductivity and low Mg diffusion.

36 MATERIALS SCIENCE↗

Fast physics-based launcher optimization for electron cyclotron current drive

With the increased urgency to design fusion pilot plants, fast optimization of electron cyclotron current drive (ECCD) launchers is paramount. Traditionally, this is done by coarsely sampling the 4D parameter space of possible launch conditions consisting of (1) the launch location (constrained to lie along the reactor vessel), (2) the launch frequency, (3) the toroidal launch angle, and (4) the poloidal launch angle. For each initial condition, a ray-tracing simulation is performed to evaluate the ECCD efficiency. Unfortunately, this approach often requires a large number of simulations (sometimes millions in extreme cases) to build up a dataset that adequately covers the plasma volume, which must then be repeated every time the design point changes. Here we adopt a different approach. Rather than launching rays from the plasma periphery and hoping for the best, we instead directly reconstruct the optimal ray for driving current at a given flux surface using a reduced physics model coupled with a commercial ray-tracing code. Repeating this throughout the plasma volume requires only hundreds of simulations, constituting a significant speedup. The new method is validated on two separate example tokamak profiles, and is shown to reliably drive localized current at the specified flux surface with the same optimal efficiency as obtained from the traditional approach.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Towards automating structural discovery in scanning transmission electron microscopy *

Abstract Scanning transmission electron microscopy is now the primary tool for exploring functional materials on the atomic level. Often, features of interest are highly localized in specific regions in the material, such as ferroelectric domain walls, extended defects, or second phase inclusions. Selecting regions to image for structural and chemical discovery via atomically resolved imaging has traditionally proceeded via human operators making semi-informed judgements on sampling locations and parameters. Recent efforts at automation for structural and physical discovery have pointed towards the use of ‘active learning’ methods that utilize Bayesian optimization with surrogate models to quickly find relevant regions of interest. Yet despite the potential importance of this direction, there is a general lack of certainty in selecting relevant control algorithms and how to balance a priori knowledge of the material system with knowledge derived during experimentation. Here we address this gap by developing the automated experiment workflows with several combinations to both illustrate the effects of these choices and demonstrate the tradeoffs associated with each in terms of accuracy, robustness, and susceptibility to hyperparameters for structural discovery. We discuss possible methods to build descriptors using the raw image data and deep learning based semantic segmentation, as well as the implementation of variational autoencoder based representation. Furthermore, each workflow is applied to a range of feature sizes including NiO pillars within a La:SrMnO 3 matrix, ferroelectric domains in BiFeO 3 , and topological defects in graphene. The code developed in this manuscript is open sourced and will be released at github.com/nccreang/AE_Workflows .

47 OTHER INSTRUMENTATION↗

A panoramic view of the Local Group dwarf galaxy NGC 6822

ABSTRACT We present a panoramic survey of the isolated Local Group dwarf irregular galaxy NGC 6822. Our photometry reaches ∼2–3 mag deeper than most previous studies and spans the widest area around the dwarf compared to any prior work. We observe no stellar overdensities in the outskirts of NGC 6822 to V ∼ 30 mag arcsec−2 and a projected radius of 16.5 kpc. This indicates that NGC 6822 has not experienced any recent interaction with a companion galaxy, despite previous suggestions to the contrary. Similarly, we find no evidence for any dwarf satellites of NGC 6822 to a limiting luminosity MV ≈ −5. NGC 6822 contains a disc of H i gas and young stars, oriented at ∼60○ to an extended spheroid composed of old stellar populations. We observe no correlation between the distribution of young stars and spheroid members. Our imaging allows us to trace the spheroid to nearly 11 kpc along its major axis, commensurate with the extent of the NGC 6822 globular cluster system. We find that the spheroid becomes increasingly flattened at larger radii, and its position angle twists by up to 40○. We use Gaia EDR3 astrometry to measure a proper motion for NGC 6822, and then sample its orbital parameter space. While this galaxy has spent the majority of its life in isolation, we find that it likely passed within the virial radius of the Milky Way ∼3–4 Gyr ago. This may explain the apparent flattening and twisting observed in the outskirts of its spheroid.

79 ASTRONOMY AND ASTROPHYSICS↗

Local primordial non-Gaussianity from the large-scale clustering of photometric DESI luminous red galaxies

ABSTRACT We use angular clustering of luminous red galaxies from the Dark Energy Spectroscopic Instrument (DESI) imaging surveys to constrain the local primordial non-Gaussianity parameter fNL. Our sample comprises over 12 million targets, covering 14 000 deg2 of the sky, with redshifts in the range 0.2 < z < 1.35. We identify Galactic extinction, survey depth, and astronomical seeing as the primary sources of systematic error, and employ linear regression and artificial neural networks to alleviate non-cosmological excess clustering on large scales. Our methods are tested against simulations with and without fNL and systematics, showing superior performance of the neural network treatment. The neural network with a set of nine imaging property maps passes our systematic null test criteria, and is chosen as the fiducial treatment. Assuming the universality relation, we find $f_{\rm NL} = 34^{+24(+50)}_{-44(-73)}$ at 68 per cent (95 per cent) confidence. We apply a series of robustness tests (e.g. cuts on imaging, declination, or scales used) that show consistency in the obtained constraints. We study how the regression method biases the measured angular power spectrum and degrades the fNL constraining power. The use of the nine maps more than doubles the uncertainty compared to using only the three primary maps in the regression. Our results thus motivate the development of more efficient methods that avoid overcorrection, protect large-scale clustering information, and preserve constraining power. Additionally, our results encourage further studies of fNL with DESI spectroscopic samples, where the inclusion of 3D clustering modes should help separate imaging systematics and lessen the degradation in the fNL uncertainty.

79 ASTRONOMY AND ASTROPHYSICS↗

Visual Tool for Assessing Stability of DER Configurations on Three-Phase Radial Networks

Here, we present a method and tool for evaluating the placement of Distributed Energy Resources (DER) on distribution circuits in order to control voltages and power flows. Our previous work described Phasor-Based Control (PBC), a novel control framework where DERs inject real and reactive power to track voltage magnitude and phase angle targets. Here, we employ linearized power flow equations and integral controllers to develop a linear state space model for PBC acting on a three-phase unbalanced network. We use this model to evaluate whether a given inverter-based DER configuration admits a stable set of controller gains, which cannot be done by analyzing controllability nor by using the Lyapunov equation. Instead, we sample over a parameter space to identify a stable set of controller gains. Our stability analysis requires only a line impedance model and does not entail simulating the system or solving an optimization problem. We incorporate this assessment into a publicly available visualization tool and demonstrate three processes for evaluating many control configurations on the IEEE 123-node test feeder (123NF).

24 POWER TRANSMISSION AND DISTRIBUTION↗

Out of Distribution Detection with Neural Network Anchoring

This is code to reproduce and build on OOD detection from the paper "Out of Distribution Detection with Neural Network Anchoring". Our goal here is to exploit heteroscedastic temperature scaling as a calibration strategy for out of distribution (OOD) detection. Heteroscedasticity here refers to the fact that the optimal temperature parameter for each sample can be different, as opposed to conventional approaches that use the same value for the entire distribution. To enable this, we propose a new training strategy called anchoring that can estimate appropriate temperature values for each sample, leading to state-of-the-art OOD detection performance across several benchmarks. Using NTK theory, we show that this temperature function estimate is closely linked to the epistemic uncertainty of the classifier, which explains its behavior. In contrast to some of the best-performing OOD detection approaches, our method does not require exposure to additional outlier datasets, custom calibration objectives, or model ensembling. Through empirical studies with different OOD detection settings - far OOD, near OOD, and semantically coherent OOD - we establish a highly effective OOD detection approach.

Thiagarajan, Jayaraman↗

GP Cosmology Surrogate v1.0

GP Cosmology Surrogate is a Python library for building and training a generalized multi-output Gaussian process (GP) framework of @takhtaganov2021cosmic. In this approach, the surrogate is constructed sequentially, guided by a Bayesian optimization acquisition function that targets reduction of emulation error in the regions most consistent with the observational data. This adaptive design concentrates computational resources where they have the greatest impact on inference accuracy. The library supports efficient training for separable GP kernels, which allows the use of Kronecker algebra to handle high-dimensional input spaces and large numbers of correlated outputs. This makes it well suited for applications such as modeling cosmological power spectra, large-scale physical simulations, and multi-output hyperparameter tuning. By combining scalable multi-output GP modeling with data-driven adaptive sampling, GPsurrogate enables parameter inference and optimization with substantially fewer simulations than conventional space-filling designs.

Lukic, Zarija [Lawrence Berkeley National Laborato↗

Measuring QCD Splittings with Invertible Networks

QCD splittings are among the most fundamental theory concepts at the LHC. We show how they can be studied systematically with the help of invertible neural networks. These networks work with sub-jet information to extract fundamental parameters from jet samples. Our approach expands the LEP measurements of QCD Casimirs to a systematic test of QCD properties based on low-level jet observables. Starting with an toy example we study the effect of the full shower, hadronization, and detector effects in detail.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

In-situ spectroscopic ellipsometry for real-time characterization of the effects of high-flux helium plasmas on tungsten surfaces

Tungsten samples were exposed to He plasmas generated by an RF source (Γ i = 3.5 x 10 16 He cm -2 s -1 , ion energy = 92 eV.) The range of exposure conditions selected here is conducive to the growth of nearsurface He bubbles, and at higher fluence, the formation of W nanotendrils ranging between 50 — 100 nm in diameter. The evolution of these surface features was probed using a fixed-angle ellipsometer (280 — 1000 nm wavelength range) with direct line-of-sight to the sample. Over the parameter space explored here, changes in the two angles (p, S) that define the polarization of the reflected light followed a distinct trajectory with increasing plasma fluence. Ex-situ ellipsometry of 22 additional tungsten specimens tested at a wide range of plasma fluences and temperatures mapped onto these in-situ results well. We used helium ion microscopy and focused ion beam profiling to provide a direct calibration of the ellipsometry measurements. Our results indicate that for a reproducible process such as the growth helium-induced surface morphologies, ellipsometry is a practical in-situ diagnostic to study how fusion plasmas modify materials. To study more general effects of plasmas on surfaces, including co-deposition and sputtering, different approaches to modelling the optical properties of the exposed surfaces are also considered.

36 MATERIALS SCIENCE↗

How Low Can You Go? Using Synthetic 3D Imagery to Drastically Reduce Real-World Training Data for Object Detection

Deep convolutional neural networks (DCNNs) currently provide state-of-the-art performance on image classification and object detection tasks, and there are many global security mission areas where such models could be extremely useful. Crucially, the success of these models is driven in large part by the widespread availability of high-quality open source data sets such as Image Net, Common Objects in Context (COCO), and KITTI, which contain millions of images with thousands of unique labels. However, global security relevant objects-of-interest can be difficult to obtain: relevant events are low frequency and high consequence; the content of relevant images is sensitive; and adversaries and proliferators seek to obscure their activities. For these cases where exemplar data is hard to come-by, even fine-tuning an existing model with available data can be effectively impossible. Recent work demonstrated that models can be trained using a combination of real-world and synthetic images generated from 3D representations; that such models can exceed the performance of models trained using real-world data alone; and that the generated images need not be perfectly realistic (Tremblay, et al., 2018). However, this approach still required hundreds to thousands of real-world images for training and fine tuning, which for sparse, global security-relevant datasets can be an unrealistic hurdle. In this research, we validate the performance and behavior of DCNN models as we drive the number of real-world images used for training object detection tasks down to a minimal set. We perform multiple experiments to identify the best approach to train DCNNs from an extremely small set of real-world images. In doing so, we: Develop state-of-the-art, parameterized 3D models based on real-world images and sample from their parameters to increase the variance in synthetic image training data; Use machine learning explainability techniques to highlight and correct through targeted training the biases that result from training using completely synthetic images; and Validate our results by comparing the performance of the models trained on synthetic data to one another, and to a control model created by fine-tuning an existing ImageNet-trained model with a limited number (hundreds) of real-world images.

97 MATHEMATICS AND COMPUTING↗

Complete Round-Robin Hydrogen Gas Analysis Capability Comparison (Milestone 2.6)

The Department of Energy (DOE) is currently evaluating strategies for the extended dry storage of aluminum-clad spent nuclear fuel (ASNF). Part of this assessment concerns the extent of radiolytic molecular hydrogen (H2) generation from the aluminum cladding’s oxyhydroxide corrosion layers. Understanding this radiation-induced process and the factors effecting it (e.g., system conditions such as temperature and gaseous environment) are essential for the development of predictive computer models to support the Technical Considerations and Challenges for Extended (> 50 yrs) Dry Storage of ASNF program. To achieve this goal and ensure that the experimental data gathered by Task 2 (Oxyhydroxide Layer Radiolytic Gas Generation Resolution) research groups (Idaho National Laboratory and Savannah River National Laboratory) are consistent, a round-robin H2 analysis capability comparison was initiated. Here we present the results from said round-robin and conclude that despite differences in sample preparation, irradiation parameters, and analytical procedures, the measured data are sufficiently consistent between the two laboratories (= 15%).

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Improving the Parameterization of Cloud and Rain Microphysics in E3SM using Novel Observationally-Constrained Bayesian Approach (Final Technical Report)

In this project, we sought to develop new cloud and rain microphysics frameworks within the Energy Exascale Earth System Model (E3SM). This work encompassed two primary avenues of research: 1) Further development of a Bayesian-based scheme called BOSS (Bayesian Observationally-constrained Statistical-physical Scheme) to represent cloud and rain microphysics, testing it in realistic high-resolution cloud models, and implementing it in E3SM; 2) Development of a methodology utilizing machine learning to enable computationally tractable use of tractable use of Markov chain Monte Carlo sampling for Bayesian parameter estimation in Earth system and cloud models. In this project, we adapted the BOSS microphysics scheme, originally formulated for rain-only, to include all liquid-phase microphysical processes for cloud and rain, in particular the processes that mediate between these two categories, for example the conversion from cloud to rain through collision and coalescence of drops. We constrained the scheme via comparison and testing against a detailed model that explicitly represents the evolution of cloud and rain particles, called a bin microphysics scheme.

54 ENVIRONMENTAL SCIENCES↗

LYNM Containment Hydrogeology Data Update for FY25

This report presents new two-phase flow parameters measured in samples of vitric tuff units located stratigraphically above the PE1 A chamber, which intersect or are near the top of the B and DL chambers. These two-phase flow data show the permeability, porosity, and mean pore size of these units to be higher than the stratigraphic units below them. These data indicate the rocks support fast gas transport (high porosity and permeability) across these units, depending on their water content.

58 GEOSCIENCES↗