Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mathematical physics methods”

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 1,603 records · Page 89
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Disentangling the Black Hole Mass Spectrum with Photometric Microlensing Surveys

Abstract From the formation mechanisms of stars and compact objects to nuclear physics, modern astronomy frequently leverages surveys to understand populations of objects to answer fundamental questions. The population of dark and isolated compact objects in the Galaxy contains critical information related to many of these topics, but is only practically accessible via gravitational microlensing. However, photometric microlensing observables are degenerate for different types of lenses, and one can seldom classify an event as involving either a compact object or stellar lens on its own. To address this difficulty, we apply a Bayesian framework that treats lens type probabilistically and jointly with a lens population model. This method allows lens population characteristics to be inferred despite intrinsic uncertainty in the lens class of any single event. We investigate this method’s effectiveness on a simulated ground-based photometric survey in the context of characterizing a hypothetical population of primordial black holes (PBHs) with an average mass of 30 M ⊙ . On simulated data, our method outperforms current black hole (BH) lens identification pipelines and characterizes different subpopulations of lenses while jointly constraining the PBH contribution to dark matter to ≈25%. Key to robust inference, our method can marginalize over population model uncertainty. We find the lower mass cutoff for stellar origin BHs, a key observable in understanding the BH mass gap, particularly difficult to infer in our simulations. This work lays the foundation for cutting-edge PBH abundance constraints to be extracted from current photometric microlensing surveys.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Metrics and Methods for Radiation Detection Algorithm Characterization for Nuclear/Radiological Source Search

This report presents a series of recommendations for data to train and evaluate radiation detection algorithms and performance metrics to evaluate these algorithms. These recommendations were formed through a community consensus approach through the Detection Radiation Algorithms Group (DRAG), a multi-institution collaboration spanning eight Department of Energy laboratories and John Hopkins Applied Physics Laboratory. This report includes recommendations on background data variability, and metrics to quantify variability, sources and shielding configurations to include in data collection campaigns and detector response variability. In addition, this report describes several anomaly detection and identification algorithms and recommends metrics to report their performance. Finally, this report ends with a discussion on machine learning algorithms.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Comparison of Numerical Methods for Modeling the Wave Field Effects Generated by Individual Wave Energy Converters and Multiple Converter Wave Farms

This numerical study compares the wave field generated by the spectral wave action balance code, SNL-SWAN, to the linear-wave boundary-element method (BEM) code, WAMIT. The objective of this study is to assess the performance of SNL-SWAN for modeling wave field effects produced by individual wave energy converters (WECs) and wave farms comprising multiple WECs by comparing results from SNL-SWAN with those produced by the BEM code WAMIT. BEM codes better model the physics of wave-body interactions and thus simulate a more accurate near-field wave field than spectral codes. In SNL-SWAN, the wave field’s energy extraction is modeled parametrically based on the WEC’s power curve. The comparison between SNL-SWAN and WAMIT is made over a range of incident wave conditions, including short-, medium-, and long-wavelength waves with various amounts of directional spreading, and for three WEC archetypes: a point absorber (PA), a pitching flap (PF) terminator, and a hinged raft (HR) attenuator. Individual WECs and wave farms of five WECs in various configuration were studied with qualitative comparisons made of wave height and spectra at specific locations, and quantitative comparisons of the wave fields over circular arcs around the WECs as a function of radial distance. Results from this numerical study demonstrate that in the near-field, the difference between SNL-SWAN and WAMIT is relatively large (between 20% and 50%), but in the far-field from the array the differences are minimal (between 1% and 5%). The resultant wave field generated by the two different numerical approaches is highly dependent on parameters such as: directional wave spreading, wave reflection or scattering, and the WEC’s power curve.

16 TIDAL AND WAVE POWER↗
DOE OSTI2026

Quantum Filtering and Analysis of Multiplicities in Eigenvalue Spectra

Fine-grained spectral properties of quantum Hamiltonians, including both eigenvalues and their multiplicities, provide useful information for characterizing many-body quantum systems as well as for understanding phenomena such as topological order. Extracting such information with small additive error is #BQP-complete in the worst case. In this work, we introduce QFAMES (quantum filtering and analysis of multiplicities in eigenvalue spectra), a quantum algorithm that efficiently identifies clusters of closely spaced dominant eigenvalues and determines their multiplicities under physically motivated assumptions, which allows us to bypass worst-case complexity barriers. QFAMES also enables the estimation of observable expectation values within targeted energy clusters, providing a powerful tool for studying quantum phase transitions and other physical properties. We validate the effectiveness of QFAMES through numerical demonstrations, including its applications to characterizing quantum phases in the transverse-field Ising model and estimating the ground-state degeneracy of a topologically ordered phase in the two-dimensional toric code model. We also generalize QFAMES to the setting of mixed initial states. Our approach offers rigorous theoretical guarantees and significant advantages over existing subspace-based quantum spectral analysis methods, particularly in terms of the sample complexity and the ability to resolve degeneracies.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

G2Aero: A Python package for separable shape tensors

G2Aero is a Python package for the design and deformation of discrete planar curves and tubular surfaces using a geometric data-driven approach. G2Aero utilizes a topology of product manifolds: the Grassmannian, $\mathcal{G}$($\mathcal{n}$, 2) - the set of 2-dimensional subspaces in $\mathbb{R}$ $\mathcal{n}$ - and the symmetric positive-definite (SPD) manifold, $\mathcal{S}^{2}_{++}$ - the set of 2x2 SPD matrices. The package provides a novel framework for representing separable deformations to shapes, which consist of stretching, scaling, rotating, and translating - also known as affine deformations - and a set of complementary deformations - which we refer to as undulation-type deformations. We focus on airfoil and blade design applications to emphasize the utility of the methods in an environment where the separation of affine and undulation-type deformations is critical. Notable functionalities of the framework for blade design include: 1) generating novel 2D (airfoil) shapes informed by a database of physically relevant airfoils, 2) building 3D blades by interpolating sequences of 2D airfoil cross-sections, and 3) generating blades with consistent perturbations along the blade span. We discuss the framework and provide examples in the context of wind energy applications, specifically wind turbine blade design. Figure 1 shows the wire frame obtained by interpolating airfoils defining the IEA 15-MW wind turbine blade and applying affine transformations corresponding to twist, chordal scaling, and bending. This, and all other figures in the paper, can be reproduced following examples and referencing supporting documentation provided in the G2Aero package.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Developing and applying quantifiable metrics for diagnostic and experiment design on Z

This project applies methods in Bayesian inference and modern statistical methods to quantify the value of new experimental data, in the form of new or modified diagnostic configurations and/or experiment designs. We demonstrate experiment design methods that can be used to identify the highest priority diagnostic improvements or experimental data to obtain in order to reduce uncertainties on critical inferred experimental quantities and select the best course of action to distinguish between competing physical models. Bayesian statistics and information theory provide the foundation for developing the necessary metrics, using two high impact experimental platforms on Z as exemplars to develop and illustrate the technique. We emphasize that the general methodology is extensible to new diagnostics (provided synthetic models are available), as well as additional platforms. We also discuss initial scoping of additional applications that began development in the last year of this LDRD.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Rapid data-driven model reduction of nonlinear dynamical systems including chemical reaction networks using ℓ 1 -regularization

We develop a new data-driven paradigm for efficient model reduction of a broad class of nonlinear dynamical systems. Our model reduction method directly enables the interpretation of key components of the dynamical system, unlike traditional projection-based model reduction methods that focus on reducing computational complexity more than interpretability. Our method is not application specific and is simple to implement on nonlinear dynamical systems arising from a variety of different fields. It requires minimal parameterization using a single parameter to trade-off between model complexity and estimation error. We use a data-driven paradigm to formulate model reduction as an efficient convex optimization problem that scales polynomially in the original size of the complex system, enabling systems with as many as thousands of components to be reduced in a matter of minutes.

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

Neural Density Estimation and Uncertainty Quantification for ChemCam Spectra [Slides]

The ChemCam instrument of Curiosity uses laser-induced breakdown spectroscopy (LIBS). It fires a laser at target and vaporizes rock surfaces, creating a plasma. Three spectrographs divide the plasma light into wavelengths for chemical analysis: ultraviolet, violet, and visible near-infrared. Regression methods (SVR, PCR, CNN) have been employed for calibration (prediction of the elemental composition of samples); however, labeled ChemCam samples are limited. Here, we focus on unsupervised learning and employ generative models from ChemCam analysis. Further, we use labels (supervised) in combination to the generative model to compute uncertainties related to predictions. We report generative modeling can be successfully applied to model real-world data. Normalizing flow models can be efficiently constructed on latent spaces for fast downstream inference. Unsupervised and supervised learning can be combined to form an uncertainty quantification framework.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

SAM User’s Guide

The System Analysis Module (SAM) is a modern system analysis tool being developed at Argonne National Laboratory for advanced non-LWR safety analysis. It aims to provide fast-running, whole-plant transient analyses capability with improved-fidelity for Sodium-cooled Fast Reactors (SFR), Lead-cooled Fast Reactors (LFR), and Molten Salt Reactors (MSR) or Fluoride-cooled High-temperature Reactors (FHR). SAM takes advantage of advances in physical modeling, numerical methods, and software engineering to enhance its user experience and usability. It utilizes an object-oriented application framework (MOOSE), and its underlying meshing and finite-element library (libMesh) and linear and non-linear solvers (PETSc), to leverage the modern advanced software environments and numerical methods. This document provides a user’s guide, which will help users understand the input description and core capabilities of the SAM code. A brief overview of the code is presented, as well as how to obtain and run it. The input syntax for various parts of the code is provided. Additionally, a number of example problems, starting with simple unit component problems to problems with increasing complexity, are provided. Because the code is still under active development, this SAM User’s Guide will evolve with periodic updates.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Measurement of the branching ratio of 16 N , 15 C , 12 B , and 13 B isotopes through the nuclear muon capture reaction in the Super-Kamiokande detector

The Super-Kamiokande detector has measured solar neutrinos for more than 25 years. The sensitivity for solar neutrino measurement is limited by the uncertainties of energy scale and background modeling. Decays of unstable isotopes with relatively long half-lives through nuclear muon capture, such as 16 N, 15 C, 12 B, and 13 B, are detected as background events for solar neutrino observations. Here, in this study, we developed a method to form a pair of stopping muon and decay candidate events and evaluated the production rates of such unstable isotopes. We then measured their branching ratios considering both their production rates and the estimated number of nuclear muon capture processes as Br⁡( 16 N) = (9.0 ± 0.1)%, Br⁡( 15 C) = (0.6 ± 0.1)%, Br⁡( 12 B) = (0.98 ± 0.18)%, Br⁡( 13 B) = (0.14 ± 0.12)%, respectively. The result for 16 N has world-leading precision at present and the results for 15 C, 12 B, and 13 B are the first branching ratio measurements for those isotopes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2018

Remote Sensing, Uncertainty Quantification, and a Theory of Data Systems; Workshop Report

The purpose of the workshop was to invite statisticians, applied mathematicians, computer scientists, data system architects, experts in remote sensing technology, and Climate and Earth System scientists to review, discuss, and plan research on issues related to large-scale, efficient analysis of distributed data using spatial statistical methods. Our motivation in organizing this event was to catalyze interchange among experts on the fast-emerging problem of analysis of distributed data. As part of SAMSI's 2017-2018 Program on Mathematical and Statistical Methods for Climate and the Earth System, a Working Group on Remote Sensing was established to address statistical and mathematical research problems in the analysis of remote sensing data. The Working Group has five subgroups: 1) Spatial Retrieval Methodology (the so-called \Spatial-X" subgroup); 2) Spatial Analysis for Hyperspectral Data (the so-called \Spatial-Y" subgroup); 3) Emulators for Complex Forward Models; 4) Optimization for Remote Sensing Retrievals; and 5) Theory of Data Systems (ToDS). The ToDS subgroup spent the first half of this academic year formulating a framework in which to consider the joint problem of a) optimizing statistical methods for environments where data are distributed and too large to move to a central location, and b) the design of data system infrastructures within which to implement those statistical methods. To x ideas, the Workshop focused on spatial statistical methods. To date there are many new spatial statistical methods designed with massive data sets in mind, in the literature. However, very few have been implemented for remote sensing data, and none have been implemented in operational settings like those used by NASA and NOAA. A major impediment to their use in these cases is that the data are not only massive, but are stored in different physical locations. These data must be brought together in some way in order to estimate spatial covariance functions, but moving data to a central location for analysis is tedious at best and impossible at worst. Some remote data reduction is almost certainly necessary, but how much? What are the consequences for inference? The fundamental issue underlying these questions is how to navigate the trade-space between costs and uncertainty in the estimates or inferences that are ultimately produced.

Braverman, Amy↗
DOE OSTI2024

A quasi-meshfree method for constructing boundary-aware reproducing bases on geometrically complex domains using manifold geodesics

In many applications, physical domains are geometrically complex making it challenging to perform coarse-scale approximation. A defeaturing process is often used to simplify the domain in preparation for approximation and analysis at the coarse scale. In this paper, a methodology is presented for constructing a coarse-scale reproducing basis on geometrically complex domains given an initial fine-scale mesh of the fully featured domain. The initial fine-scale mesh can be of poor quality and extremely refined. The construction of the basis functions begins with a coarse-scale covering of the domain and generation of weighting functions with local support. Manifold geodesics are used to define distances within the local support for general applicability to non-convex domains. Conventional moving least squares is used to construct the coarse-scale reproducing basis. Applications in quasi-interpolation and linear elasticity are presented.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Simulation of gas mixture dynamics in a pipeline network using explicit staggered-grid discretization

Here we develop an explicit staggered finite difference discretization scheme for simulating the transport of highly heterogeneous gas mixtures through pipeline networks. This study is motivated by the proposed blending of hydrogen into natural gas pipelines to reduce end use carbon emissions while using existing pipeline systems throughout their planned lifetimes. Our computational method accommodates an arbitrary number of constituent gases with very different physical properties that may be injected into a network with significant spatiotemporal variation. In this setting, the gas flow physics are highly location- and time- dependent, so that local composition and nodal mixing must be accounted for. The resulting conservation laws are formulated in terms of pressure, partial densities and flows, and volumetric and mass fractions of the constituents. We include non-ideal equations of state that employ linear approximations of gas compressibility factors, so that the pressure dynamics propagate locally according to a variable wave speed that depends on mixture composition and density. We derive compatibility relationships for network edge boundary values that are more complex than for a homogeneous gas. The simulation method is evaluated on initial boundary value problems for a single pipe and a small network, is cross-validated with a lumped element simulation, and used to demonstrate a local monitoring and control policy for maintaining allowable concentration levels.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Summary Report Of The FY25 Reactor Physics Verification And Validation Exercises In The Advanced Reactor Technologies - Gas-cooled Reactor Program

Valdiation and verification of numerical tools is critical for ensuring reasonable predictions for design scoping, licensing, and safety analsyis. In this report, two reactor physics verification and validation exercises are presented. The first of these exercises focuses on burnup analysis with data from the Advanced Gas Reactor (AGR) program. Simulations are performed with Monte Carlo N-Particle (MCNP) and are compared with the experimental measurements for the AGR 1 and 2 experiments that utilize both UCO and UO2 fuel. The second exercises utilizes data from the HTR-Proteus experiments to perform reactor physics validation. Specifications of the experimental facility are provdied, along with a demonstration of initial modeling efforts in Serpent for one of the determistic packing experiments. Both cases are part of the Generation-IV international forum (GIF) Very High-Temperature Reactor (VHTR) Computational Methods, Validation, and Benchmarking (CMVB) program, an international collaborative organization dedicated to the verification and validation of High-Temperature Gas-Cooled Reactor (HTGR) analysis. Participation in the CMVB allows the US Department of Energy (DOE) to leverage these existing validation activities to provide extra value through benchmarking activities with other CMVB members.

and Benchmarking (CMVB) program↗
DOE OSTI2020

Investigation of Mass Transfer and Sorption in CO 2 /Brine/Rock Systems via In Situ FT-IR

CO 2 geological storage in deep saline formations is considered a promising method to mitigate anthropogenic CO 2 emissions and, thereby, minimize changes to the Earth’s atmosphere. A fundamental understanding of CO 2 mass transfer and sorption phenomena in brine-saturated reservoir formations is necessary to understand the long-term fate of injected CO 2 as it is subjected to different (physical, dissolution, and mineral) trapping mechanisms. In this work, we investigate CO 2 sorption in brine-saturated and dry Mt. Simon sandstone samples via in-situ Fourier Transform infrared spectroscopy (FT-IR) at elevated pressures, ranging from 0.3 MPa to 8.3 MPa, at a temperature of 50 ⁰C. The FT-IR spectra of bulk-phase CO 2 were simultaneously recorded under the same conditions. For bulk-phase CO 2 , we observed, in agreement with past studies, a doublet peak at 2361 cm -1 and 2336 cm -1 and another peak (ν 2 bending mode) at 667 cm -1 . With increasing pressure, the position of the peak at 667 cm -1 remains invariant, however, when crossing into the supercritical region the doublet peak degenerates onto a single peak at 2336 cm -1 with a barely visible shoulder at 2361 cm -1 . The bulk CO 2 data provide a perfect fit for Beer’s law for the whole range of pressure conditions. For the dry sample, the IR spectrum is experimentally indistinguishable from the bulk CO 2 spectrum, signifying that if physical adsorption occurs to any significant extent, the adsorbed CO 2 molecules are not substantially more rotationally constrained than the dense bulk CO 2 molecules. For the brine-saturated sample, we observe a strong band centered at 2342 cm -1 and a small companion peak at 2360-2361 cm -1 that degenerates into a barely visible shoulder peak at the higher pressures. The 2342 cm -1 band has been previously observed by other investigators for CO 2 dissolved in bulk water/brine as well during its adsorption on a variety of other wet natural porous media. We observe no peaks corresponding to bicarbonate or carbonate bulk species, which correlates well with the prior literature on similar low-pH aqueous solutions. The integrated peak area for the CO 2 sorbed in the brine-saturated sample correlates linearly with its solubility in the same bulk brine, as measured separately via a PVT-cell approach. Here, this validates the accuracy of both techniques, and the potential of the FT-IR method to be used in the study of mass transfer and adsorption in such systems. To that effect, a simple mathematical model is presented to analyze the FT-IR data to determine the CO 2 effective diffusivity in the brine-saturated sandstone sample.

58 GEOSCIENCES↗
DOE OSTI2024

Diffusion quantum Monte Carlo approach to the polaritonic ground state

Making and using polaritonic states (i.e., hybrid electron-photon states) for chemical applications has recently become one of the most prominent and active fields that connects the communities of chemistry and quantum optics. Modeling of such polaritonic phenomena using ab initio approaches calls for new methodologies, leading to the reinvention of many commonly used electronic structure methods, such as Hartree-Fock, density functional, and coupled cluster theories. Here in this work, we explore the formally exact diffusion quantum Monte Carlo approach to obtain numerical solutions to the polaritonic ground state during the dissociation of the H 2 molecular system. Here, we examine various electron-nuclear-photon properties throughout the dissociation, such as changes to the minimum of the cavity Born-Oppenheimer surface, the localization of the electronic wave function, and the average mode occupation. Finally, we directly compare our results to that obtained with state-of-the-art, yet approximate, polaritonic coupled cluster approaches.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Machine Learning-driven Molecular Design for Therapeutic Discovery

The ongoing novel coronavirus pandemic (COVID-19) has highlighted the need for new therapeutics to counter the threat of emerging viral pathogens. The main proteases are a promising target for developing antiviral inhibitors. In this work, we utilized a novel combination of artificial intelligence-driven iterative design of covalent inhibitor candidates, physics-based computational modeling of protein-inhibitor interactions, and “All in One” Native MS biophysical assay screening and characterization of therapeutic candidates. With our existing expertise in hit generation using a particular scaffold as a starting point, we first generated tens of thousands of compounds that preserve the key scaffold. In order to optimize the candidates, we calculated about 136 descriptors consisting of 2D and 3D features for molecules targeting the SARS-CoV-2 Main protease (Mpro). These compounds were initially filtered according to properties and further sorted by predicted binding affinity using our automated docking modeling and machine learning methods. We tested a handful of candidates and identified two as inhibitors of Mpro with micromolar affinities.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

New charged-particle transport computational capability: the SIT code An L-4 milestone

We have developed a new high-fidelity code for direct transport of charged particles using the simple integral transport method. The code can be coupled to the outputs of any hydrodynamical simulation code in 1-D, 2-D or 3-D. In this report we summarize the formalism involved in treating complex transport problems. We present physical examples wherein we have used the code to calculate the transport of alpha particles. Future work is planned to study the sensitivity of hydrodynamical mix to charged-particle radiochemistry and reaction-in-flight neutrons for the complex inertial confinement fusion problems encountered at NIF and at the Z-machine.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗