Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “High resolution data analysis”

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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Grand Challenges of Hydrologic Modeling for Food-Energy-Water Nexus Security in High Mountain Asia

Climate-influenced changes in hydrology affect water-food-energy security that may impact up to two billion people downstream of the High Mountain Asia (HMA) region. Changes in water supply affect energy, industry, transportation, and ecosystems (agriculture, fisheries) and as a result, also affect the region's social, environmental, and economic fabrics. Sustaining the highly interconnected food-energy-water nexus (FEWN) will be a fundamental and increasing challenge under a changing climate regime. High variability in topography and distribution of glaciated and snow-covered areas in the HMA region, and scarcity of high resolution (in-situ) data make it difficult to model and project climate change impacts on individual watersheds. We lack basic understanding of the spatial and temporal variations in climate, surface impurities in snow and ice such as black carbon and dust that alter surface albedo, and glacier mass balance and dynamics. These knowledge gaps create challenges in predicting where and when the impact of changes in river flow will be the most significant economically and ecologically. In response to these challenges, the United States National Aeronautics and Space Administration (NASA) established the High Mountain Asia Team (HiMAT) in 2016 to conduct research to address knowledge gaps. This paper summarizes some of the advances HiMAT made over the past 5 years, highlights the scientific challenges in improving our understanding of the hydrology of the HMA region, and introduces an integrated assessment framework to assess the impacts of climate changes on the FEWN for the HMA region. The framework, developed under a NASA HMA project, links climate models, hydrology, hydropower, fish biology, and economic analysis. The framework could be applied to develop scientific understanding of spatio-temporal variability in water availability and the resultant downstream impacts on the FEWN to support water resource management under a changing climate regime.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Air quality impact of the Northern California Camp Fire of November 2018

The Northern California Camp Fire that took place in November 2018 was one of the most damaging environmental events in California history. Here, we analyze ground-based station observations of airborne particulate matter that has a diameter <2.5 µm (PM 2.5 ) across Northern California and conduct numerical simulations of the Camp Fire using the Weather Research and Forecasting model online coupled with chemistry (WRF-Chem). Simulations are evaluated against ground-based observations of PM 2.5 , black carbon, and meteorology, as well as satellite measurements, such as Tropospheric Monitoring Instrument (TROPOMI) aerosol layer height and aerosol index. The Camp Fire led to an increase in Bay Area PM 2.5 to over 50 µg m –3 for nearly 2 weeks, with localized peaks exceeding 300 µg m –3 . Using the Visible Infrared Imaging Radiometer Suite (VIIRS) high-resolution fire detection products, the simulations reproduce the magnitude and evolution of surface PM 2.5 concentrations, especially downwind of the wildfire. The overall spatial patterns of simulated aerosol plumes and their heights are comparable with the latest satellite products from TROPOMI. WRF-Chem sensitivity simulations are carried out to analyze uncertainties that arise from fire emissions, meteorological conditions, feedback of aerosol radiative effects on meteorology, and various physical parameterizations, including the planetary boundary layer model and the plume rise model. Downwind PM 2.5 concentrations are sensitive to both flaming and smoldering emissions over the fire, so the uncertainty in the satellite-derived fire emission products can directly affect the air pollution simulations downwind. Our analysis also shows the importance of land surface and boundary layer parameterization in the fire simulation, which can result in large variations in magnitude and trend of surface PM 2.5 . Inclusion of aerosol radiative feedback moderately improves PM 2.5 simulations, especially over the most polluted days. Results of this study can assist in the development of data assimilation systems as well as air quality forecasting of health exposures and economic impact studies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Quantification of the REE 3+ aqua ions and chloride species in aqueous fluids by in situ Raman spectroscopy using perturbations of the water band

Acidic NaCl-rich aqueous fluids play a crucial role in forming hydrothermal rare earth elements (REE) mineral deposits. Aqueous REE mobility is mostly controlled by the stabilities of REE 3+ and REE chloride species. Our current knowledge of REE speciation is based on solubility data, thermodynamic models and in situ spectroscopic measurements, sometimes coupled with molecular simulations. Here, in this study, we investigate Nd and Yb speciation in pH2 Cl-bearing solutions at 25 °C and 0.1 MPa with variable Cl/REE ratios using Raman Spectroscopy in solutions with 0.1 to 0.6 mol/kg NdCl 3 or YbCl 3 and 0.2 to 3.2 mol/kg NaCl. Due to the challenges in resolving the REE-Cl band, we developed a new method using the water vibrational mode and multivariate curve resolution (MCR) analysis. The Raman spectra for the vibrational band of water (2700 to 3900 cm –1 ) were collected at 25 °C and fitted by three Gaussian sub-peaks, then quantified using MCR analysis to de-convolute the water band into bulk H 2 O and the perturbations caused by of Cl – , REE 3+ , and REE chloride species. REE speciation based on the perturbations of the water band indicates that REE 3+ aqua ions dominate acidic solutions at 25 °C, but up to ~20 mol% YbCl 2+ forms at high YbCl 3 concentrations. The new method is promising for quantifying in situ speciation of the REE 3+ aqua ions and REE chloride species in aqueous fluids while providing information on the hydration of ions. This method improves our molecular level understanding of REE aqueous species stability and their role in REE mobilization during fluid-rock interaction.

58 GEOSCIENCES↗
DOE OSTI2023

A novel integrated time-resolved array avalanche photodiode detection system for nuclear resonant scattering measurements

Here, the nuclear resonant scattering (NRS) experiment requires photon-counting detectors with high time resolution, short dead time, large dynamic range, low noise, and large detection area. An 8-channel avalanche photodiode (APD) array detector system with high integrity, flexibility, and reliability has been developed to adapt to the demands of NRS experiments. The detector system mainly consists of four key parts: (i) an array-APD sensor, (ii) 8-channel integrated fast preamplifiers, (iii) the time-to-digital converter readout electronics, and (iv) a data acquisition system and EPICS support software. Remarkably, the system exhibits a time resolution of better than 500 ps and has a sufficiently low noise level, allowing for the lowest detection energy threshold of 4 keV. The performance of the new array-APD system as well as its real application in nuclear forward scattering (NFS) and nuclear resonant inelastic x-ray scattering (NRIXS) experiments was tested in two synchrotron facilities. With the new system, the NFS signal very close to the prompt electronic scattering signal can be extracted. Thanks to the customized EPICS-areaDetector-based control software, NRIXS spectra can be readily measured with time and energy information of the NRIXS signal stored in the raw data, which is promising for developing NRIXS data analysis in the time domain. The array-APD detector can be deployed for nuclear resonant scattering experiments at various synchrotron radiation facilities.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Asymptotically Compatible Reproducing Kernel Collocation and Meshfree Integration for Nonlocal Diffusion

Reproducing kernel (RK) approximations are meshfree methods that construct shape functions from sets of scattered data. We present an asymptotically compatible (AC) RK collocation method for nonlocal diffusion models with Dirichlet boundary condition. The numerical scheme is shown to be convergent to both nonlocal diffusion and its corresponding local limit as nonlocal interaction vanishes. The analysis is carried out on a special family of rectilinear Cartesian grids for a linear RK method with designed kernel support. The key idea for the stability of the RK collocation scheme is to compare the collocation scheme with the standard Galerkin scheme, which is stable. In addition, assembling the stiffness matrix of the nonlocal problem requires costly computational resources because high-order Gaussian quadrature is necessary to evaluate the integral. We thus provide a remedy to the problem by introducing a quasi-discrete nonlocal diffusion operator for which no numerical quadrature is further needed after applying the RK collocation scheme. The quasi-discrete nonlocal diffusion operator combined with RK collocation is shown to be convergent to the correct local diffusion problem by taking the limits of nonlocal interaction and spatial resolution simultaneously. The theoretical results are then validated with numerical experiments. We additionally illustrate a connection between the proposed technique and an existing optimization based approach based on generalized moving least squares.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Characterization of HRRR-simulated rotor layer wind speeds and clouds along the coast of California

Stratocumulus clouds, with their low cloud base and top, affect the atmospheric boundary layer wind and turbulence profile, thereby modulating wind energy resources. GOES satellite data reveal an abundance of stratocumulus clouds in the late spring and summer months off the coast of northern and central California, where there are active plans to deploy floating offshore wind farms at two lease areas (near Morro Bay and Humboldt). Since the fall of 2020, two buoys equipped with multiple instruments, including Doppler lidar, have been deployed for about 1 year in these wind farm lease areas to assess the rotor layer wind conditions in these locations. The objective of this study is to evaluate how well the High-Resolution Rapid Refresh (HRRR) model represents stratocumulus cloud characteristics and turbine-relevant rotor layer winds (surface to 300 m) by comparing HRRR simulations with buoy and satellite observations. We first find that the HRRR model reproduces the seasonal cycle of cloud top height reasonably well in these regions. However, during the warm season – especially at Morro Bay – the HRRR-simulated stratocumulus clouds tend to have lower tops by about 150 m and exhibit weaker diurnal cycles than satellite observations. Our analysis also shows that rotor layer wind speeds and vertical shear are stronger at Humboldt than at Morro Bay, and both are generally stronger under clear-sky conditions. Finally, the HRRR model bias in rotor layer wind speed is small under cloudy conditions but larger and dependent on observed wind speed under clear skies. Specifically, HRRR underestimates wind speeds at Morro Bay and overestimates them at Humboldt under clear-sky conditions.

17 WIND ENERGY↗
DOE OSTI2021

Omics Lethal Human Viruses Project Profiling of the Host Response to MERS-CoV Infection, Processed Experimental Dataset Catalog

Middle East Respiratory Syndrome coronavirus (MERS-CoV) is classified as a Category C priority pathogen (Coronaviridae) by the National Institute of Allergy and Infectious Diseases (NIAID), and is known to cause severe respiratory disease with high mortality rates in humans. Lethal host-pathogen invasion mechanisms and the cellular intricacies behind these fatal infections still remain unclear. The NIAID Modeling Host Responses to Understand Severe Human Virus Infections Research Program project (2013-2018) aimed to develop an improved comprehensive understanding of the host response to a suite of viruses causing lethal infections leveraging a systems biology approach. Herein, PNNL sub-projects provide a never before released comprehensive infectious disease collection of primary and secondary transformation multi-Omics data profiling a series of priority pathogen primary experimental studies for enhanced open-access to viral Omics datasets and project lifecycle metadata. Secondary host-pathogen viral dataset downloads contain one or more statistically processed (normalization data transformation) quantitative dataset collections resulting in qualitative expression analyses of primary host-pathogen experimental study designs. Leveraging unique high-resolution Omics capabilities for proteomics (P), metabolomics (M), lipidomics (L), and transcriptomics (T) dataset downloads each have a direct relationship to a primary sample submission corresponding to a specific MERS-CoV [NCBITAXON:1335626] experimental infection study. Host sample types include human lung adenocarcinoma cells ["Calu-3", BTO:0002750], human bronchial epithelial cells ["Calu-3 clone 2B4"; BTO:0002022], primary human fibroblasts ["FB"; BTO:0000452], primary human airway epithelial cells ["HAE"; BTO:0005571], human microvascular endothelial cells ["HMVE"; BTO:0003123], and whole mouse lung [BTO:0000763] tissue collections.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Separate Effects Validation for Subcooled Boiling in CTF

In 2010, the U.S. Department of Energy created its first Energy Innovation Hub, which is focused on developing high-fidelity and high-resolution modeling and simulation (M&S) tools for modeling of light water reactors (LWRs). This hub, the Consortium for Advanced Simulation of LWRs (CASL), has developed an LWR simulation tool called the Virtual Environment for Reactor Applications (VERA). The multi-physics capability of VERA is achieved through the coupling of single-physics codes, including CTF (the CASL version of Coolant Boiling in Rod Arrays— Three Field (COBRA-TF)), Michigan Parallel Characteristics Transport (MPACT), BISON, and Materials Performance and Optimization (MPO) Advanced Model for Boron Analysis (MAMBA). As part of its M&S efforts, CASL has identified various challenge problems, including Crud Induced Power Shift (CIPS), Crud-Induced Localized Corrosion (CILC), Pellet-Cladding Interaction (PCI), and Departure from Nucleate Boiling (DNB). This work addresses CASL milestone L2:VVI.P19.03, which focuses on uncertainty quantification of crud, which is relevant to both CIPS and CILC. This is achieved through an analysis and separate effects validation of the thermal hydraulic phenomenon known as subcooled boiling. As part of this work, various sources of experimental data are examined and compared to different options for empirical modeling of subcooled boiling. Through this analysis, a complete understanding of the underlying models and their implementation details are understood. A subset of these data are incorporated into a separate effects validation study of CTF. The Westinghouse Advanced Loop Tester (WALT) and Rohsenow experiments are modeled, and it is shown that the newly-implemented Gorenflo correlation is more accurate than the existing Chen and Thom correlations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Single-nuclei transcriptome analysis of the shoot apex vascular system differentiation in Populus

ABSTRACT Differentiation of stem cells in the plant apex gives rise to aerial tissues and organs. Presently, we lack a lineage map of the shoot apex cells in woody perennials – a crucial gap considering their role in determining primary and secondary growth. Here, we used single-nuclei RNA-sequencing to determine cell type-specific transcriptomes of the Populus vegetative shoot apex. We identified highly heterogeneous cell populations clustered into seven broad groups represented by 18 transcriptionally distinct cell clusters. Next, we established the developmental trajectories of the epidermis, leaf mesophyll and vascular tissue. Motivated by the high similarities between Populus and Arabidopsis cell population in the vegetative apex, we applied a pipeline for interspecific single-cell gene expression data integration. We contrasted the developmental trajectories of primary phloem and xylem formation in both species, establishing the first comparison of vascular development between a model annual herbaceous and a woody perennial plant species. Our results offer a valuable resource for investigating the principles underlying cell division and differentiation conserved between herbaceous and perennial species while also allowing us to examine species-specific differences at single-cell resolution.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Tracking seasonal variability in plant traits from spaceborne PRISMA and NEON AOP across forest types and ecoregions

Plant traits serve as critical indicators of how plants adapt to environmental changes and influence ecosystem functions. While airborne hyperspectral remote sensing effectively maps plant traits through detailed reflectance properties, it is limited by cost and scale, making large-scale and temporal studies challenging. The recently launched spaceborne hyperspectral imager, PRecursore IperSpettrale della Missione Applicativa (PRISMA), offers frequent, large scale and high-fidelity observations on a spatial resolution of 30 m and a revisit time of around 29 days, making it suitable for large-scale seasonal trait mapping. However, their potential remains largely unexplored. This study developed a multi-stage framework by leveraging the PRISMA spaceborne hyperspectral data and National Ecological Observatory Network (NEON) Airborne Observation Platform (AOP) hyperspectral data to investigate the seasonal dynamics of four key plant traits — chlorophyll content, carotenoid content, equivalent water thickness, and nitrogen content — across eleven NEON sites representing diverse forest types and ecoregions in the contiguous U.S. Our results demonstrated that PRISMA hyperspectral data can reliably track seasonal variability in plant traits, achieving overall R 2 values ranging from 0.78 to 0.88 and normalized root mean square error (NRMSE) values ranging from 5.4% to 8.4% for the four traits. Seasonal patterns revealed bell-shaped trajectories for chlorophyll and carotenoids, while equivalent water thickness decreased steadily across most sites, driven by structural changes during leaf maturation and senescence. Nitrogen content exhibited less pronounced seasonal variation but followed expected nutrient resorption patterns. Analysis of environmental drivers showed that seasonal variability is primarily controlled by solar radiation and day length in northern sites, vapor pressure in semi-arid regions, and temperature in mid-southeastern sites. Spatial variability, meanwhile, was primarily driven by soil properties, particularly during the peak growing season. However, the influence of soil variables slightly declines toward the end of the season at several sites, as climatic factors become more prominent. This study highlights the capability of PRISMA, and potentially other similar spaceborne hyperspectral data for large-scale, time-series plant trait mapping and provides valuable insights into the interactions between plant traits and environmental factors. In conclusion, these findings contribute to advancing our understanding of plant functional ecology and improving predictions of ecosystem responses to environmental changes.

Environmental drivers↗
DOE OSTI2024

Projected income data under different shared socioeconomic pathways for Washington state

Abstract High-resolution income projections under different Shared Socioeconomic Pathways (SSPs) are essential for the climate change research communities to devise climate change adaptation and mitigation strategies. To generate income projections for Washington state, we obtain state-level GDP per capita projections and convert them into projected annual household income. The resulting state-level income projections are subsequently downscaled to the census block-level based on the Longitudinal Origin-Destination Employment Statistics (LODES) dataset. For accuracy assessment, we downscale historical income data from state- level to block- and block group-level and compare the downscaled results against the actual income data from LODES. County-level accuracy assessment is also conducted based on American Community Survey. The results demonstrate a good agreement (Average R 2 of 0.67, 0.8, and 0.99 for block-, block group-, and county-level, respectively) between the downscaled income data and the reference data, thereby validating the methodology employed. Our approach is applicable to other states for income projections, which can be utilized by a broader audience, including those involved in demographic analysis, economic research, and urban planning.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2020

Neutral Pion Electroproduction and development of a Neutral Particle Spectrometer

Protons and neutrons, i.e. nucleons, are the basic building blocks of the matter in the visible universe. The strong force binds the nucleons to form nuclei. The electromagnetic force forms the atoms by binding the electrons with the nuclei. The electromagnetic interaction is well understood by Quantum Electrodynamics (QED), which shows the most precise predictability amongst all the theories in physics. In QED, charges interact with each other by exchanging photons. The Quantum Chromodynamics (QCD) describes the strong interaction. Its degrees of freedom are quarks and gluons, the fundamental constituents of the nucleons. The quarks interact with each other by exchanging gluons. However, unlike QED, the gluons interact amongst themselves. This feature of self bindings of the gluons confines the quarks and gluons in the nucleons/hadrons, never to be seen as free. In order to study some of the features of QCD, such as confinement or the structure of the nucleon, one usually needs to rely on experiments. Electromagnetic probes, governed by the well-understood QED, are excellent tools to probe the nucleon. In general, different scales, e.g. electron beam energies, probe different regions of the nucleon. At low energy, of the order of a few GeV, the electron probes the nucleon in the valence quark region. As its energy increases, the electron probes the sea quark and gluon regions. The study of the nucleon structure in all these regions is needed to fully understand QCD. Form factors and parton distribution functions measured from elastic scattering and deep inelastic scattering of leptons off nucleons have provided a partial view of the internal structure of the nucleon. In the mid-1990s, Generalized Parton Distributions (GPDs) were developed. These new objects are a generalization of the form factors and parton distribution functions, but contain richer information on the nucleon internal structure. GPDs are accessible experimentally by deep exclusive reactions. Deeply virtual Compton scattering (DVCS) and deeply virtual meson production (DVMP) are some examples. The first dedicated DVCS/DVMP experiment took place in 2004 in Hall A at The Thomas Jefferson National Accelerator Facility, i.e. Jefferson Lab, in Virginia, U.S.A. A new DVCS/DVMP experiment, after the beam energy upgrade of Jefferson Lab, was carried out in Hall A in a wider kinematic range. Its data were taken from 2014 to 2016. In Hall C at Jefferson Lab, the next DVCS/DVMP experiment will take place. The Hall C experiment will further exploit the kinematic range with higher precision. A Neutral Particle Spectrometer (NPS) is in development to measure DVCS/DVMP events under high background conditions. Jefferson Lab will provide the highest precision data in the valence quark region for various exclusive reactions. The Electron-Ion Collider (EIC) is a future experimental facility currently planned to start operations around 2030 in the U.S.A. Its high energy and high luminosity will probe the sea quark and gluon regions providing answers to the outstanding questions of QCD, in particular in the region where matter is dominated by gluons. First of all, this document describes the data analysis and results of the Hall A neutral pion electroproduction off the proton, from the data taken in 2014-2016. Later, some of the developments towards the construction of the electromagnetic calorimeter of the NPS for the upcoming DVCS/DVMP experiment in Hall C are presented. Finally, one of the candidate materials for the EIC calorimeter, a glass scintillator, will be briefly introduced. I have participated to all these projects, in collaboration with many colleagues. I present in this thesis my contributions to each of these projects. My contributions to the neutral pion data analysis were focused on background subtractions on the calorimeter, acceptance calculations, and the estimation of the systematic uncertainty associated to the event selection cuts. Some necessary information on calibrations of the detectors and data analysis methods are also described. In the NPS project, I performed background dose calculations and energy and position resolution studies of the calorimeter, all using Monte Carlo simulations, with realistic geometries of the experimental apparatus. Characterization of the crystals of the calorimeter was also done. Additionally, I measured the radiation hardness of some glass scintillator in its early stage of development. In order to have a future reference when the glass calorimeter prototype will be tested, I simulated the energy resolution of the prototype.

Ko, Ho-San↗
DOE OSTI2024

Optimal Spectral Resolution for Infrared Studies of Solids and Liquids

Due to a legacy originating in the limited capability of early computers, the spectroscopic resolution used in Fourier transform infrared spectroscopy and other systems has largely been implemented using only powers of two for more than 50 years. In this study, we investigate debunking the spectroscopic lore of, e.g., using only 2, 4, 8, or 16 cm −1 resolution and determine the optimal resolution in terms of both (i) a desired signal-to-noise ratio and (ii) efficient use of acquisition time. The study is facilitated by the availability of solids and liquids reference spectral data recorded at 2.0 cm −1 resolution and is based on an examination in the 4000–400 cm −1 range of 61 liquids and 70 solids spectra, with a total analysis of 4237 peaks, each of which was also examined for being singlet/multiplet in nature. Of the 1765 liquid bands examined, only 27 had widths <5 cm −1 . Of the 2472 solid bands examined, only 39 peaks have widths <5 cm −1 . For both the liquid and solid bands, a skewed distribution of peak widths was observed: For liquids, the mean peak width was 24.7 cm −1 but the median peak width was 13.7 cm −1 , and, similarly, for solids, the mean peak width was 22.2 cm −1 but the median peak width was 11.2 cm −1 . While recognizing other studies may differ in scope and limiting the analysis to only room temperature data, we have found that a resolution to resolve 95% of all bands is 5.7 cm −1 for liquids and 5.3 cm −1 for solids; such a resolution would capture the native linewidth (not accounting for instrumental broadening) for 95% of all the solids and liquid bands, respectively. After decades of measuring liquids and solids at 4, 8, or 16 cm −1 resolution, we suggest that, when accounting only for intrinsic linewidths, an optimized resolution of 6.0 cm −1 will capture 91% of all condensed-phase bands, i.e., broadening of only 9% of the narrowest of bands, but yielding a large gain in signal-to-noise with minimal loss of specificity.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Machine learning changes the rules for flux limiters

Learning to integrate non-linear equations from highly resolved direct numerical simulations has seen recent interest for reducing the computational load for fluid simulations. Here, we focus on determining a flux-limiter for shock capturing methods. Focusing on flux limiters provides a specific plug-and-play component for existing numerical methods. Since their introduction, an array of flux limiters has been designed. Using the coarse-grained Burgers' equation, we show that flux-limiters may be rank-ordered in terms of their log-error relative to high-resolution data. We then develop a theory to find an optimal flux-limiter and present flux-limiters that outperform others tested for integrating Burgers' equation on lattices with [Formula: see text], and 2x, 3x, 4x, and 8x coarse-grainings. We train a continuous piecewise linear limiter by minimizing the mean-squared misfit to six-grid point segments of high-resolution data, averaged over all segments. While flux limiters are generally designed to have an output of φ(r) = 1 at a flux ratio of r = 1, our limiters are not bound by this rule and yet produce a smaller error than standard limiters. Here we find that our machine learned limiters have distinctive features that may provide new rules-of-thumb for the development of improved limiters. Additionally, we use our theory to learn flux-limiters that outperform standard limiters across a range of values (as opposed to at a specific fixed value) of coarse-graining, number of discretized bins, and diffusion parameter. This demonstrates the ability to produce flux limiters that should be more broadly useful than standard limiters for general applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Initial Characterization and Optimization of the Liquid Sampling-Atmospheric Pressure Glow Discharge Ionization Source Coupled to an Orbitrap Mass Spectrometer for the Determination of Plutonium

Plutonium measurements are essential to the nuclear forensics and safeguards community. The liquid sampling-atmospheric pressure glow discharge (LS-APGD) microplasma ionization source coupled with an Orbitrap mass spectrometer is a proven platform for uranium isotope ratio determinations. Here, this work expands the LS-APGD-Orbitrap platform capabilities by reporting the first-ever analysis of plutonium with the LS-APGD and the first-ever measurement of elemental plutonium with an Orbitrap mass spectrometer. This coupling has the potential to dramatically reduce the complex sample manipulations required for traditional analysis techniques employed for actinide isotope ratio determinations. As a first step toward the goal of simultaneous uranium and plutonium isotope ratio determinations, the initial characterization and optimization of the platform for the detection of plutonium are reported. Collision-induced dissociation modality settings were optimized to reduce water-related and other molecular clusters containing plutonium, maximizing 242 Pu 16 O 2 + responses. A design of experiments study was conducted to optimize the discharge conditions of the dual-electrode LS-APGD toward the responsivity of 242 Pu 16 O 2 +. The measurement sensitivity was determined from a Pu response curve, yielding a limit of detection of 10 fg (absolute) of total analyte when data was collected and processed with a Spectroswiss FTMS Booster X2 data acquisition system. Additionally, plutonium and uranium were measured in a simultaneous acquisition, and each analyte remained unaffected by the other. It is believed that the LS-APGD-Orbitrap platform could be a valuable addition to the nuclear forensics’ toolbox and, indeed, other scientific disciplines and regulatory communities in which rapid, high-resolution plutonium determinations are paramount.

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

K-space algorithmic reconstruction (KAREN): a robust statistical methodology to separate Bragg and diffuse scattering

Diffuse scattering occurring in the Bragg diffraction pattern of a long-range-ordered structure represents local deviation from the governing regular lattice. However, interpreting the real-space structure from the diffraction pattern presents a significant challenge because of the dramatic difference in intensity between the Bragg and diffuse components of the total scattering function. In contrast to the sharp Bragg diffraction, the diffuse signal has generally been considered to be a weak expansive or continuous background signal. In this paper, using 1D and 2D models, it is demonstrated that diffuse scattering in fact consists of a complex array of high-frequency features that must not be averaged into a low-frequency background signal. To evaluate the actual diffuse scattering effectively, an algorithm has been developed that uses robust statistics and traditional signal processing techniques to identify Bragg peaks as signal outliers which can be removed from the overall scattering data and then replaced by statistically valid fill values. This method, described as a 'K-space algorithmic reconstruction' (KAREN), can identify Bragg reflections independent of prior knowledge of a system's unit cell. KAREN does not alter any data other than that in the immediate vicinity of the Bragg reflections, and reconstructs the diffuse component surrounding the Bragg peaks without introducing discontinuities which induce Fourier ripples or artifacts from underfilling `punched' voids. The KAREN algorithm for reconstructing diffuse scattering provides demonstrably better resolution than can be obtained from previously described punch-and-fill methods. The superior structural resolution obtained using the KAREN method is demonstrated by evaluating the complex ordered diffuse scattering observed from the neutron diffraction of a single plastic crystal of CBr 4 using pair distribution function analysis.

36 MATERIALS SCIENCE↗
DOE OSTI2020

TOI-1235 b: A Keystone Super-Earth for Testing Radius Valley Emergence Models around Early M Dwarfs

Small planets on close-in orbits tend to exhibit envelope mass fractions of either effectively zero or up to a few percent depending on their size and orbital period. Models of thermally driven atmospheric mass loss and of terrestrial planet formation in a gas-poor environment make distinct predictions regarding the location of this rocky/nonrocky transition in period–radius space. Here we present the confirmation of TOI-1235 b (P = 3.44 days, r{sub p}=1.738{sub −0.076}{sup +0.087} R{sub ⊕}), a planet whose size and period are intermediate between the competing model predictions, thus making the system an important test case for emergence models of the rocky/nonrocky transition around early M dwarfs (R {sub s} = 0.630 ± 0.015 R{sub ⊙}, M {sub s} = 0.640 ± 0.016 M{sub ⊙}). We confirm the TESS planet discovery using reconnaissance spectroscopy, ground-based photometry, high-resolution imaging, and a set of 38 precise radial velocities (RVs) from HARPS-N and HIRES. We measure a planet mass of 6.91{sub −0.85}{sup +0.75} M{sub ⊕}, which implies an iron core mass fraction of 20{sub −12}{sup +15}% in the absence of a gaseous envelope. The bulk composition of TOI-1235 b is therefore consistent with being Earth-like, and we constrain an H/He envelope mass fraction to be <0.5% at 90% confidence. Our results are consistent with model predictions from thermally driven atmospheric mass loss but not with gas-poor formation, suggesting that the former class of processes remains efficient at sculpting close-in planets around early M dwarfs. Our RV analysis also reveals a strong periodicity close to the first harmonic of the photometrically determined stellar rotation period that we treat as stellar activity, despite other lines of evidence favoring a planetary origin (P=21.8{sub −0.8}{sup +0.9} days, m{sub p}sin i=13.0{sub −5.3}{sup +3.8} M{sub ⊕}) that cannot be firmly ruled out by our data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

CERF: IM3 Projected Western US Power Plant Locations

Overview The Capacity Expansion Regional Feasibility (CERF) model is an open-source geospatial python package that provides new power plant locations at a 1km resolution. The model ingests U.S. state or regional-scale electricity system capacity expansion plans, such as those produced by the Global Change Analysis Model (GCAM-USA), and identifies feasible, site-specific locations for individual new power plants (renewable and non-renewable). CERF combines high-resolution geospatial suitability analyses with an economic algorithm that selects individual plant siting locations based on grid interconnection costs and the locational marginal value of new generation. The model incorporates a wide range of dynamic constraints and opportunities, such as protected lands, population density, existing infrastructure, and water availability. This dataset provides CERF power plant siting results for IM3 Phase 2 simulations across eight different scenarios for the Western US through 2055. The scenarios include combinations of two Shared Socioeconomic Pathways (SSP3 and SSP5) with four high-resolution climate projections specific to the United States (see, https://tgw-data.msdlive.org/). These climate projections include "hotter" and "cooler" variants for two Representative Concentration Pathways (RCP4.5 and RCP8.5). The resulting eight simulations are: rcp45cooler_ssp3 rcp45cooler_ssp5 rcp45hotter_ssp3 rcp45hotter_ssp5 rcp85cooler_ssp3 rcp85cooler_ssp5 rcp85hotter_ssp3 rcp85hotter_ssp5 CERF siting results in this dataset correspond to capacity expansion plans in the GCAM-USA IM3 Phase 2 simulation data and are available for each of the above scenarios. Data Details Temporal Range: 2015-2055 in 5-year timesteps. Note that 2015 is the experiment base year and 2020 and beyond represent model simulation years. Spatial Range: Plant locations are provided for the eleven states in the Western US including Arizona, California, Colorado, Idaho, Montana, New Mexico, Nevada, Oregon, Utah, Washington, and Wyoming. Spatial Resolution: 1 km-squared, provided in x and y coordinates Geospatial Projection: Albers Equal Area Conic (ESRI:102003) File Type: csv The dataset contains subdirectories for each of the eight scenarios described in the overview. Each scenario folder contains two subfolders with the following information: 1. Power Plant Data This directory contains a single .csv file of power plant locations for both pre-existing (non-CERF sited plants in operation in 2015) and new (CERF-sited) power plants across the temporal range along with additional CERF model output parameters for CERF-sited plants. Plant with a siting year earlier than 2020 correspond to facilities that are operational leading into the first timestep CERF simulation. For a more detailed description of CERF model output parameters, see the CERF model documentation. Note that the cerf_plant_id parameter is unique within each scenario file but not across scenario files. Parameter Descriptions scenario - Name of scenario cerf_plant_id - Unique siting identifier cerf_sited - If True, indicates that plant was sited by CERF model. If False, indicates pre-existing facility region_name - Name of region (state) tech_id - Technology ID tech_name - Full generation technology name inclusive of cooling type (if applicable) and additional characteristics tech_simple - Simplified generation technology type unit_size_mw - Power plant unit size (MW) xcoord - X coordinate in the default CRS (meters) ycoord - Y coordinate in the default CRS (meters) index - Index position in the flattend 2D array buffer_in_km - Exclusion buffer around site (km) sited_year - Year of siting retirement_year - Year of retirement lmp_zone - Locational marginal price (LMP) zone ID locational_marginal_price_usd_per_mwh - Locational marginal price ($/MWh) generation_mwh_per_year - Generation output (MWh/yr) operating_cost_usd_per_year - Cost of plant operations ($/yr) net_operational_value - Net operational value based on LMP and and operating costs ($/yr) interconnection_cost - Cost of interconnection for transmission & gas pipeline (if applicable) net_locational_cost -- Difference of interconnection cost and operating value ($/yr) capacity_factor_fraction - Capacity factor (fraction) carbon_capture_rate_fraction - Carbon capture rate (fraction) fuel_co2_content_tons_per_btu - Fuel CO2 content (tons/Btu) fuel_price_usd_per_mmbtu - Fuel price ($/MMBtu) fuel_price_esc_rate_fraction - Fuel price escalation rate (fraction) heat_rate_btu_per_kWh - Heat rate (Btu/kWh) lifetime_yrs - Technology lifetime for annuity (years) operational_life_yrs - Operational lifetime for retirement (years) variable_om_usd_per_mwh - Variable operation and maintenance costs of yearly capacity use ($/MWh) variable_om_esc_rate_fraction - Variable operation and maintenance costs escalation rate (fraction) carbon_tax_usd_per_ton - Carbon tax ($/ton) carbon_tax_esc_rate_fraction - Carbon tax escalation rate (fraction) 2. Storage Data This directory contains information on new and pre-existing energy storage facilities operational in each timestep along with various storage operational parameters. The 2015 timestep provides pre-existing energy storage data and corresponds with facilities that are operational leading into the first model simulation timestep. Note that coordinates in the storage files correspond to the interconnection point on the grid (substation location), not individual energy storage locations. Energy storage is added in a cumulative process at each given interconnection point. That is, each individual file provides the total operational storage capacity interconnected to the specified substation for the given timestep, inclusive of previously installed storage at that location and new storage installed in that timestep at that location. Parameters scenario - Name of scenario timestep - Simulation timestep name - Unique storage identifier s_typ - Type of energy storage technology (battery or pumped storage hydro) s_node - Node ID of interconnecting substation xcoord - X coordinate in the default CRS (meters) ycoord - Y coordinate in the default CRS (meters) charge_rate - Maximum charge rate (power capacity) of storage system (MW) discharge_rate - Maximum discharge rate (power capacity) of storage system (MW) duration - Duration of storage system (hours) max_SoC - Allowed maximum state of charge (energy capacity) of storage system (MWh) min_SoC -Allowed minimum state of charge (energy capacity) of storage system (MWh) charge_eff - Efficiency of charge (fraction between 0 and 1) discharge_eff - Efficiency of discharge (fraction between 0 and 1) Acknowledgment IM3 is a multi-institutional effort led by Pacific Northwest National Laboratory and supported by the U.S. Department of Energy's Office of Science as part of research in MultiSector Dynamics, Earth and Environmental Systems Modeling Program.

CERF↗