Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “data inference”

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 847 records · Page 47

Effective field theory analysis of the Coulomb breakup of the one-neutron halo nucleus $$^{19}$$C

Abstract We analyse the Coulomb breakup of $$^{19}$$ 19 C measured at 67 A MeV at RIKEN. We use the Coulomb-Corrected Eikonal (CCE) approximation to model the reaction and describe the one-neutron halo nucleus $$^{19}$$ 19 C within Halo Effective Field Theory (Halo EFT). At leading order we obtain a fair reproduction of the measured cross section as a function of energy and angle. The description is insensitive to the choice of optical potential, as long as it accurately represents the size of $${}^{18}$$ 18 C. It is also insensitive to the interior of the $${}^{19}$$ 19 C wave function. Comparison between theory and experiment thus enables us to infer asymptotic properties of the ground state of $$^{19}$$ 19 C: these data put constraints on the one-neutron separation energy of this nucleus and, for a given binding energy, can be used to extract an asymptotic normalisation coefficient (ANC). These results are confirmed by CCE calculations employing next-to-leading order Halo EFT descriptions of $${}^{19}$$ 19 C: at this order the results for the Coulomb breakup cross section are completely insensitive to the choice of the regulator. Accordingly, this reaction can be used to constrain the one-neutron separation energy and ANC of $${}^{19}$$ 19 C.

Physics↗

Bayesian homodyne and heterodyne tomography

Continuous-variable (CV) photonic states are of increasing interest in quantum information science, bolstered by features such as deterministic resource state generation and error correction via bosonic codes. Data-efficient characterization methods will prove critical in the fine-tuning and maturation of such CV quantum technology. Although Bayesian inference offers appealing properties—including uncertainty quantification and optimality in mean-squared error—Bayesian methods have yet to be demonstrated for the tomography of arbitrary CV states. Here we introduce a complete Bayesian quantum state tomography workflow capable of inferring generic CV states measured by homodyne or heterodyne detection, with no assumption of Gaussianity. As examples, we demonstrate our approach on experimental coherent, thermal, and cat state data, obtaining excellent agreement between our Bayesian estimates and theoretical predictions. Our approach lays the groundwork for Bayesian estimation of highly complex CV quantum states in emerging quantum photonic platforms, such as quantum communications networks and sensors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

xRage simulations of magnetic velocity gauges used in gas gun shock initiation experiments

Shock initiation gas gun experiments of PBX 9012 have been performed by Burns and Chiquete [2020] using magnetic velocity gauges as the principal diagnostic. In private communications, the experimenters noted that compatibility issues required a different glue to bond the gauges to PBX 9012 then previously used for experiments with PBX 9501 and PBX 9502. The silicon glue used is more viscous and can result in glue layers up to 50 microns, which is about the gauge thickness. This raises the question on what affect the thicker glue layers have on the gauge response. In particular, the width of the lead shock can affect the accuracy of inferring shock locus points used for the shock pressure of Pop plot data points and for Hugoniot data to calibrate the reactants EOS.

36 MATERIALS SCIENCE↗

SAIL Aerosol Vertical Profiles (SAIL-AVP) Field Campaign Report

Aerosols are central to understanding the water cycle within mountainous regions, but a complete understanding of this role cannot be provided without vertically resolved observations, particularly for aerosol-cloud interactions, since we simultaneously need to know aerosol information and meteorology near cloud bases. The goal of SAIL Aerosol Vertical Profiles (SAIL-AVP) was to perform tethered balloon system (TBS) measurements during different seasons to elucidate process-level understanding of the aerosols and associated meteorological conditions within complex mountainous terrain. Measurements were conducted during the deployment of the U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) user facility’s second Mobile Facility (AMF2) as part of the Surface Atmosphere Integrated Field Laboratory (SAIL) campaign in the East River Watershed of the Upper Colorado River Basin (UCRB) in southwestern Colorado. In mid-latitude continental interior mountain ranges, and especially at SAIL, aerosol vertical profiles and their relationships with meteorology are, almost without exception, inferred from a combination of ground-based and remote-sensing data. These flights provided many more dimensions to aerosol-cloud interactions that are unavailable to those inferential approaches: the TBS collected in situ information in the vertical, optical, chemical, and biological dimensions for aerosol-cloud interactions. Furthermore, the collocation of these flights with longstanding collaborative resources in the region, including the ongoing surface and subsurface hydrologic observations from the DOE’s Watershed Function Science Focus Area (SFA), produced a unique set of atmospheric observations that are complemented by existing land-surface and subsurface (e.g., groundwater) observations.

54 ENVIRONMENTAL SCIENCES↗

Specimen sizing and remaining life sample calculations for high temperature reactor material surveillance

Advanced high-temperature nuclear reactors rely on structural components that will operate for decades under combined mechanical, thermal, and environmental loading. Materials surveillance programs are a promising strategy for managing the resulting uncertainty in long-term structural integrity by monitoring degradation in service using passively actuated mechanical test articles. Previous reports have developed a simplified, spreadsheetimplementable framework for sizing these test articles and for inferring accumulated creep damage and remaining life from ex-situ test data. This report advances that work toward practical deployment by providing a sample-problem book: a collection of worked, end-toend examples in which the ASME Section III, Division 5 design analysis of a representative high-temperature reactor component is carried through to a sized surveillance article, verified with detailed finite-element analysis, and concluded with a remaining-life assessment based on an assumed ex-situ creep-rate measurement on the retrieved specimen. The report also summarizes ongoing ANL engagement with ASTM Committee E10 on Nuclear Technology and Applications toward drafting a standard covering surveillance procedures for advanced reactors.

Barua, Bipul (ORCID:0000000247184113)↗

CyanoCyc cyanobacterial web portal

CyanoCyc is a web portal that integrates an exceptionally rich database collection of information about cyanobacterial genomes with an extensive suite of bioinformatics tools. It was developed to address the needs of the cyanobacterial research and biotechnology communities. The 277 annotated cyanobacterial genomes currently in CyanoCyc are supplemented with computational inferences including predicted metabolic pathways, operons, protein complexes, and orthologs; and with data imported from external databases, such as protein features and Gene Ontology (GO) terms imported from UniProt. Five of the genome databases have undergone manual curation with input from more than a dozen cyanobacteria experts to correct errors and integrate information from more than 1,765 published articles. CyanoCyc has bioinformatics tools that encompass genome, metabolic pathway and regulatory informatics; omics data analysis; and comparative analyses, including visualizations of multiple genomes aligned at orthologous genes, and comparisons of metabolic networks for multiple organisms. CyanoCyc is a high-quality, reliable knowledgebase that accelerates scientists’ work by enabling users to quickly find accurate information using its powerful set of search tools, to understand gene function through expert mini-reviews with citations, to acquire information quickly using its interactive visualization tools, and to inform better decision-making for fundamental and applied research.

59 BASIC BIOLOGICAL SCIENCES↗

Dimensional Reduction for Sampled Priors and Application to Photometric Redshift Distributions

A typical Bayesian inference on the values of some parameters of interest q from some data D involves running a Markov Chain (MC) to sample from the posterior $p$($q$,$n$|$D$) $\propto$ $\mathcal{L}$($D$|$q$,$n$)$p$(q)$p$($n$), where n are some nuisance parameters with a separable prior. In some cases, the nuisance parameters are high-dimensional, and their prior p(n) is itself defined only by a set of samples that have been drawn from some other MC. The MC for the posterior will typically require evaluation of p(n) at arbitrary values of n, i.e., one needs to provide a density estimator over the full n space from the provided samples. But the high dimensionality of n hinders both the density estimation and the efficiency of the MC for the posterior. We describe a solution to this problem: a linear compression of the n space into a much lower-dimensional space u, which projects away directions in n space that cannot appreciably alter $\mathcal{L}$. The algorithm for doing so is a slight modification to principal components analysis, and is less restrictive on p(n) than other proposed solutions to this issue. We demonstrate this “mode projection” technique using the analysis of 2-point correlation functions of weak lensing fields and galaxy density in the Dark Energy Survey, where n is a binned representation of the redshift distribution n(z) of the galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

Constraining the Dense Matter Equation of State with New NICER Mass–Radius Measurements and New Chiral Effective Field Theory Inputs

Pulse profile modeling of X-ray data from the Neutron Star Interior Composition Explorer is now enabling precision inference of neutron star mass and radius. Combined with nuclear physics constraints from chiral effective field theory (χEFT), and masses and tidal deformabilities inferred from gravitational-wave detections of binary neutron star mergers, this has led to a steady improvement in our understanding of the dense matter equation of state (EOS). Here, we consider the impact of several new results: the radius measurement for the 1.42 M ⊙ pulsar PSR J0437−4715 presented by Choudhury et al., updates to the masses and radii of PSR J0740+6620 and PSR J0030+0451, and new χEFT results for neutron star matter up to 1.5 times nuclear saturation density. Using two different high-density EOS extensions—a piecewise-polytropic (PP) model and a model based on the speed of sound in a neutron star (CS)—we find the radius of a 1.4 M ⊙ (2.0 M ⊙ ) neutron star to be constrained to the 95% credible ranges ${12.28}_{-0.76}^{+0.50}$ km (${12.33}_{-1.34}^{+0.70}$ km) for the PP model and ${12.01}_{-0.75}^{+0.56}$ km (${11.55}_{-1.09}^{+0.94}$ km) for the CS model. The maximum neutron star mass is predicted to be ${2.15}_{-0.16}^{+0.14}$M ⊙ and ${2.08}_{-0.16}^{+0.28}$M ⊙ for the PP and CS models, respectively. We explore the sensitivity of our results to different orders and different densities up to which χEFT is used, and show how the astrophysical observations provide constraints for the pressure at intermediate densities. Moreover, we investigate the difference R 2.0 − R 1.4 of the radius of 2 M ⊙ and 1.4 M ⊙ neutron stars within our EOS inference.

gravitational wave sources↗

Characterization of DESI fiber assignment incompleteness effect on 2-point clustering and mitigation methods for DR1 analysis

We present an in-depth analysis of the fiber assignment incompleteness in the Dark EnergySpectroscopic Instrument (DESI) Data Release 1 (DR1). This incompleteness is caused by therestricted mobility of the robotic fiber positioner in the DESI focal plane, which limits thenumber of galaxies that can be observed at the same time, especially at small angular separations.As a result, the observed clustering amplitude is suppressed in a scale-dependent manner, which,if not addressed, can severely impact the inference of cosmological parameters. We discuss themethods adopted for simulating fiber assignment on mocks and data. In particular, we introducethe fast fiber assignment (FFA) emulator, which was employed to obtain the power spectrumcovariance adopted for the DR1 full-shape analysis. We present the mitigation techniques,organised in two classes: measurement stage and model stage. We then use high fidelity mocks as areference to quantify both the accuracy of the FFA emulator and the effectiveness of the differentmeasurement-stage mitigation techniques. This complements the studies conducted in a parallelpaper for the model-stage techniques, namely the θ-cut approach. We find that pairwiseinverse probability (PIP) weights with angular upweighting recover the “true” clustering in allthe cases considered, in both Fourier and configuration space. Notably, we present the first everpower spectrum measurement with PIP weights from real data.

79 ASTRONOMY AND ASTROPHYSICS↗

Trust Your Gut: Comparing Human and Machine Inference from Noisy Visualizations

People commonly utilize visualizations not only to examine a given dataset, but also to draw generalizable conclusions about the underlying models or phenomena. Prior research has compared human visual inference to that of an optimal Bayesian agent, with deviations from rational analysis viewed as problematic. However, human reliance on non-normative heuristics may prove advantageous in certain circumstances. We investigate scenarios where human intuition might surpass idealized statistical rationality. In two experiments, we examine individuals’ accuracy in characterizing the parameters of known data-generating models from bivariate visualizations. Our findings indicate that, although participants generally exhibited lower accuracy compared to statistical models, they frequently outperformed Bayesian agents, particularly when faced with extreme samples. Participants appeared to rely on their internal models to filter out noisy visualizations, thus improving their resilience against spurious data. However, participants displayed overconfidence and struggled with uncertainty estimation. They also exhibited higher variance than statistical machines. Our findings suggest that analyst gut reactions to visualizations may provide an advantage, even when departing from rationality. These results carry implications for designing visual analytics tools, offering new perspectives on how to integrate statistical models and analyst intuition for improved inference and decision-making. The data and materials for this paper are available at https://osf.io/qmfv6

human-machine collaboration↗

Sensitivity Coefficients Calculated for the Prompt Neutron Decay Constant at or Near Delayed Critical [Abstract]

Experimenters at Los Alamos National Laboratory (LANL) measure the prompt neutron decay constant for many experiments at the National Criticality Experiments Research Center (NCERC) to infer reactivity and the effective neutron multiplication factor. These quantities are very important for nuclear criticality safety and validating nuclear data. Uncertainty in measures of criticality of an experimental configuration can be determined prior to physically performing the experiment by applying first order perturbation theory to Monte Carlo codes, such as MCNP®. The first order perturbation theory produces first derivatives of some nuclear parameter to nuclear data (e.g., cross section data). This first derivative is commonly referred to as a sensitivity coefficient. Currently, the MCNP® software has the capability of computing effective neutron multiplication factor sensitivity coefficients to cross section data. This work builds off of this MCNP® capability and the first order perturbation theory to provide a method of calculating sensitivity coefficients for the prompt neutron decay constant at or near delayed critical to cross section data. The prompt neutron decay constant sensitivity coefficient calculated in this work does not depend on any modification of the MCNP® source code. Prompt neutron decay constant sensitivity coefficient calculations can be used to infer reactivity and effective neutron multiplication factor sensitivity coefficient values as well. By investigating the trends of prompt neutron decay constant sensitivity coefficients for nuclide-reaction pairs across energy spectra, experiments can be designed to maximize or minimize the uncertainty in the prompt neutron decay constant in a particular energy region, which can lead to further optimization studies. Prompt neutron decay constant sensitivity coefficients will be calculated for the Jezebel benchmark. Subsequently, these sensitivity coefficients will be used in a data assimilation process to determine if there is or are optimal experiments that can be performed to provide insight into adjustments of uncertain/inaccurate cross section data. Specifically, the effect of the prompt neutron decay constant sensitivity coefficients on the nuclear data-induced uncertainty in the effective neutron multiplication factor will be examined

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

GPU-Accelerated Machine Learning Inference as a Service for Computing in Neutrino Experiments

Machine learning algorithms are becoming increasingly prevalent and performant in the reconstruction of events in accelerator-based neutrino experiments. These sophisticated algorithms can be computationally expensive. At the same time, the data volumes of such experiments are rapidly increasing. The demand to process billions of neutrino events with many machine learning algorithm inferences creates a computing challenge. We explore a computing model in which heterogeneous computing with GPU coprocessors is made available as a web service. The coprocessors can be efficiently and elastically deployed to provide the right amount of computing for a given processing task. With our approach, Services for Optimized Network Inference on Coprocessors (SONIC), we integrate GPU acceleration specifically for the ProtoDUNE-SP reconstruction chain without disrupting the native computing workflow. With our integrated framework, we accelerate the most time-consuming task, track and particle shower hit identification, by a factor of 17. This results in a factor of 2.7 reduction in the total processing time when compared with CPU-only production. For this particular task, only 1 GPU is required for every 68 CPU threads, providing a cost-effective solution.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Selecting Appropriate Model Complexity: An Example of Tracer Inversion for Thermal Prediction in Enhanced Geothermal Systems

Abstract A major challenge in the inversion of subsurface parameters is the ill‐posedness issue caused by the inherent subsurface complexities and the generally spatially sparse data. Appropriate simplifications of inversion models are thus necessary to make the inversion process tractable and meanwhile preserve the predictive ability of the inversion results. In this study, we investigate the effect of model complexity on fracture aperture inversion and thermal performance prediction in a field‐scale EGS model. Principal component analysis was used to map the aperture field to a low‐dimensional latent space. The complexity of the inversion model was quantitatively represented by the percentage of total variance in the original aperture fields preserved by the latent space. Tracer, pressure and flow rate data were used to invert for fracture aperture through an ensemble‐based inversion method, and the inferred aperture field was used to predict thermal performance. With an over‐simplified aperture model, ensemble collapse occurred. The inverted aperture models failed to resolve necessary flow and transport features, leading to a biased thermal performance prediction. A complex aperture model involved excessive features and was prone to overinterpreting the inversion data. Both the tracer/pressure/flow rate data reproduction and thermal prediction showed significant uncertainties, making it difficult to properly estimate long‐term thermal performance. Fortunately, our results indicate that there exists an appropriate model complexity which can simultaneously match inversion data and predict thermal performance with an acceptable uncertainty. The quality of the fit of tracer data appears to be a useful indicator of such an appropriate model complexity.

15 GEOTHERMAL ENERGY↗

Validation of the DESI DR2 Ly$α$ forest full-shape analysis

We present the validation of the Dark Energy Spectroscopic Instrument (DESI) Data Release 2 (DR2) Lyman-$α$ (Ly$α$) forest full-shape analysis. This analysis combines three-dimensional Ly$α$ forest auto-correlations and cross-correlations with quasars to extract information from both the baryon acoustic oscillation (BAO) feature and the broadband clustering signal, with primary emphasis on the Alcock-Paczynski (AP) measurement. Compared to the DESI DR1 analysis, the DR2 validation uses substantially larger and more realistic mock datasets, including CoLoRe 2LPT and AbacusSummit Ly$α$ forest simulations. The modeling framework is also improved through analytic marginalization over small scales ($<10$$h^{-1}$Mpc) and the impact of ultraviolet background fluctuations. The validation program was completed prior to unblinding and defines quantitative requirements for the cosmological parameters of interest, which are evaluated using hundreds of mock realizations. We further test the analysis through independent fits to the auto- and cross-correlations, multiple catalog splits, and a broad suite of analysis and modeling variations applied to both mocks and blinded observational data. We find that the BAO and AP parameters satisfy all validation requirements and remain stable across all tests. In contrast, mock studies reveal a significant bias in the inferred growth-rate parameter $fσ_8$, leading us to exclude this measurement from the final analysis. The consistency across mocks, data splits, and robustness tests demonstrates that the DR2 Ly$α$ full-shape analysis provides a reliable and substantially improved broadband AP measurement over previous Ly$α$ forest studies.

Herbold, M. [Chicago U., KICP; Ohio State U.] (ORC↗

Uncertainty quantification in scientific machine learning: Methods, metrics, and comparisons

Neural networks (NNs) are currently changing the computational paradigm on how to combine data with mathematical laws in physics and engineering in a profound way, tackling challenging inverse and ill-posed problems not solvable with traditional methods. However, quantifying errors and uncertainties in NN-based inference is more complicated than in traditional methods. This is because in addition to aleatoric uncertainty associated with noisy data, there is also uncertainty due to limited data, but also due to NN hyperparameters, overparametrization, optimization and sampling errors as well as model misspecification. Although there are some recent works on uncertainty quantification (UQ) in NNs, there is no systematic investigation of suitable methods towards quantifying the total uncertainty effectively and efficiently even for function approximation, and there is even less work on solving partial differential equations and learning operator mappings between infinite-dimensional function spaces using NNs. In this work, we present a comprehensive framework that includes uncertainty modeling, new and existing solution methods, as well as evaluation metrics and post-hoc improvement approaches. Further, to demonstrate the applicability and reliability of our framework, we present an extensive comparative study in which various methods are tested on prototype problems, including problems with mixed input-output data, and stochastic problems in high dimensions. In the Appendix, we include a comprehensive description of all the UQ methods employed. Further, to help facilitate the deployment of UQ in Scientific Machine Learning research and practice, we present and develop in [1] an open-source Python library (github.com/Crunch-UQ4MI/neuraluq), termed NeuralUQ, that is accompanied by an educational tutorial and additional computational experiments.

11 physics-informed neural networks↗

Ferromagnetic cluster-glass phase in Ca(Co 1-x Ir x ) 2-y As 2 crystals

In this study, single crystals of Ca(Co 1-x Ir x ) 2-y As 2 with 0≤x≤0.35 and 0.10≤y≤0.14 have been grown using the self-flux technique and characterized by single-crystal x-ray diffraction (XRD), energy-dispersive x-ray spectroscopy, magnetization M, and magnetic susceptibility χ measurements versus temperature T, magnetic field H, and time t, and heat-capacity C p (H,T) measurements. The XRD refinements reveal that all the Ir-substituted crystals crystallize in a collapsed-tetragonal structure as does the parent CaCo 2-y As 2 compound. A small 3.3% Ir substitution for Co in CaCo 1.86 As 2 drastically lowers the A-type antiferromagnetic (AFM) transition temperature T N from 52 to 23 K with a significant enhancement of the Sommerfeld electronic heat-capacity coefficient. The A-type AFM structure consists of ab-plane layers of spins ferromagnetically aligned along the c axis with AFM alignment of the spins in adjacent layers along this axis. The positive Weiss temperatures obtained from Curie-Weiss fits to the χ(T>T N ) data indicate that the dominant magnetic interactions are ferromagnetic (FM) for all x. A magnetic phase boundary is inferred to be present between x=0.14 and x=0.17 from a discontinuity in the x dependencies of the effective moment and Weiss temperature in the Curie-Weiss fits. FM fluctuations that strongly increase with increasing x are also revealed from the χ(T) data. Furthermore, the magnetic ground state for x≥0.17 is a spin glass as indicated by hysteresis in χ(T) between field-cooled and zero-field-cooled measurements and from the relaxation of M in a small field that exhibits a stretched-exponential time dependence. The spin glass has a small FM component to the ordering and is hence inferred to be comprised of small FM clusters. The competing AFM and FM interactions along with crystallographic disorder associated with Ir substitution are inferred to be responsible for the development of a FM cluster-glass phase. A logarithmic T dependence of Cp at low T for x=0.14 is consistent with the presence of significant FM quantum fluctuations. This composition is near the T=0 boundary at x≈0.16 between the A-type AFM phase containing ferromagnetically-aligned layers of spins and the FM cluster-glass phase.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Morphostatistical characterization of the spatial galaxy distribution through Gibbs point processes

ABSTRACT This paper proposes a morphostatistical characterization of the galaxy distribution through spatial statistical modelling based on inhomogeneous Gibbs point processes. The galaxy distribution is supposed to exhibit two components. The first one is related to the major geometrical features exhibited by the observed galaxy field, here, its corresponding filamentary pattern. The second one is related to the interactions exhibited by the galaxies. Gibbs point processes are statistical models able to integrate these two aspects in a probability density, controlled by some parameters. Several such models are fitted to real observational data via the ABC shadow algorithm. This algorithm provides simultaneous parameter estimation and posterior-based inference, hence allowing the derivation of the statistical significance of the obtained results.

79 ASTRONOMY AND ASTROPHYSICS↗