Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sampled value”

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,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Superconducting gap of H 3 S measured by tunnelling spectroscopy

Several hydrogen-rich superconductors have been found to show unprecedentedly high critical temperatures, stimulating investigations into the nature of the superconductivity in these materials. Although their macroscopic superconducting properties are established, microscopic insights into the pairing mechanism remains unclear. Here we characterize the superconducting gap structure in the high-temperature superconductor H 3 S and its deuterium counterpart D 3 S by performing tunnelling spectroscopy measurements. The tunnelling spectra reveal that H 3 S and D 3 S both have a fully gapped structure, which could be well described by a single s-wave Dynes model, with gap values 2Δ of approximately 60 meV and 44 meV, respectively. Furthermore, we observed gap features of another likely H-depleted H x S superconducting phase in a poorly synthesized hydrogen sulfide sample. Our work offers direct experimental evidence for superconductivity in the hydrogen-rich superconductor H 3 S from a microscopic perspective. It validates the phonon-mediated mechanism of superconducting pairing and provides a foundation for further understanding the origins of high-temperature superconductivity in hydrogen-rich compounds.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Nanoscale control of LaAlO 3 /SrTiO 3 conductance beyond the polar catastrophe

Conductive atomic force microscope (c-AFM) lithography can be utilized to create a wide range of LaAlO 3 /SrTiO 3 (LAO/STO)-based nanoelectronic devices in a reconfigurable manner. Experiments were generally performed with intrinsically insulating LAO/STO heterostructures, with LAO thickness less than the critical value at which a polar catastrophe takes place [<4 unit cell (u.c.)]. Here, we use inductively coupled plasma reactive ion etching (ICPRIE) to fabricate c-AFM “canvases” on intrinsically conducting LAO/STO samples with ≥4 u.c. LAO. We observe that its interfacial two-dimensional electron gas (2DEG) can be pinched off and then switched back on by c-AFM lithography. Nanowires created with initially conductive LAO/STO interfaces have an order-of-magnitude longer lifetime in ambient conditions, when compared to an identically created 3.4 u.c. LAO/STO nanowire. We also demonstrate key nanoscale properties such as ballistic transport in a quasi-one-dimensional electron waveguide at a 5 u.c. LAO/STO interface. This approach frees c-AFM-written nanodevice designs from time constraints in air associated with <4 u.c. LAO/STO heterostructures.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Electronic transport, thermal transport, thermal expansion, and magnetization in the strongly correlated metal LaNi⁢O 3

Perovskite structured LaNiO 3 is a strongly correlated metal with intriguing thermal and magnetic properties. The volume dependence of calculated and measured physical properties can add additional critical information to develop a more in-depth understanding of this strongly correlated phenomenon. Taking advantage of recent single crystal LaNiO 3 growth using the floating-zone method, we have measured the thermal expansion, the magnetostriction, and the pressure dependence of the magnetic susceptibility, which then allows derivation of the Grüneisen parameters γ e = $\frac{dlnN(E_F)}{d lnV}$, γ χ = $\frac{dlnχ}{d lnV}$, as well as of electric and thermal transport properties. We simulate the volume dependence of structural and magnetic properties using Density Functional Theory calculations at the Generalized Gradient Approximation level. A large discrepancy between experimental values and calculated ones suggests that strong correlations are likely to be dynamic in nature. This study also provides a side-by-side comparison of measurements in single crystal and polycrystalline samples of LaNiO 3 to elucidate intrinsic materials properties. A broad hump at high temperatures in the temperature dependence of magnetization found in the single crystal sample of LaNiO 3 has been rationalized by a model that includes the influence of electron correlations on the Landau diamagnetism.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Improved $^{95}\mathrm{Mo}$ neutron resonance parameters and astrophysical reaction rates

We report improved 95 Mo neutron resonance parameters and reaction rates are important for nuclear astrophysics, testing nuclear models, and nuclear criticality safety. However, despite many previous neutron-capture and total cross-section measurements on this nuclide, there still is much room for improvement as well as several discrepancies. For example, there are very few firm resonance spin and parity assignments; average resonance parameters are available only for each parity, the currently recommended astrophysical reaction rate results in disagreements between stellar models and meteoric isotopic anomalies, and there are substantial disagreements in the neutron-capture cross section at low energies important for nuclear criticality safety. To obtain an improved set of neutron resonance parameters and astrophysical reaction rates for 95 Mo. High-resolution neutron-capture and transmission data were measured at the Oak Ridge Electron Linear Accelerator (ORELA) using highly isotopically enriched 95 Mo samples. The neutron-capture apparatus, data reduction, and analysis were improved so that information contained in the γ-ray cascade following neutron capture were used to assign resonance J π values. Following this, simultaneous analysis of the new neutron-capture and transmission data was used to obtain resonance energies, gamma widths, and neutron widths and their uncertainties to a maximum energy of 10 keV. Accurate neutron-capture cross sections also were obtained for the unresolved resonance region to a maximum energy of 500 keV and, together with the new resonance parameters, used to calculate the astrophysical reaction rates in the temperature range from 5 to 30 keV. A vastly improved set of 95Mo neutron resonance parameters and an astrophysical reaction rate accurate to about 3% were obtained. Firm J π assignments were determined for 261 of the 314 observed resonances. This is a very large improvement over the previously published 32 firm J π assignments for 108 resonances. Also, the number of resonances having both firm J π assignments and Γ γ values was increased by almost a factor of 24—from 11 to 261. Neutron- and total-radiation-width distributions and average resonance spacings, average total radiation widths, and neutron strength functions were obtained for the six different s- and p-wave possibilities. Parameters for the lowest s-wave resonance, which is most important for criticality benchmarks, were obtained with high accuracy. Simple modification of the neutron-capture apparatus and expansion and improvement of data analysis techniques led to a large increase in firm J π assignments for 95 Mo neutron resonances. The resulting astrophysical reaction rate is 20%–30% larger than the currently recommended rate at s-process temperatures, which should lead to better agreement between stellar models and meteoric isotopic anomalies. The neutron-capture cross section at low energies is substantially larger than recommended in the latest evaluation, which is problematical for criticality benchmarks. The average resonance spacing as a function of spin and parity is significantly different from current models. The total-radiation-width distributions are significantly broader than predicted by theory and show significant departures from the expected Gaussian shapes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Inverse Design of Two-Dimensional Airfoils Using Conditional Generative Models and Surrogate Log-Likelihoods

Abstract This paper shows how to use conditional generative models in two-dimensional (2D) airfoil optimization to probabilistically predict good initialization points within the vicinity of the optima given the input boundary conditions, thus warm starting and accelerating further optimization. We accommodate the possibility of multiple optimal designs corresponding to the same input boundary condition and take this inversion ambiguity into account when designing our prediction framework. To this end, we first employ the conditional formulation of our previous work BézierGAN–Conditional BézierGAN (CBGAN)—as a baseline, then introduce its sibling conditional entropic BézierGAN (CEBGAN), which is based on optimal transport regularized with entropy. Compared with CBGAN, CEBGAN overcomes mode collapse plaguing conventional GANs, improves the average lift-drag (Cl/Cd) efficiency of airfoil predictions from 80.8% of the optimal value to 95.8%, and meanwhile accelerates the training process by 30.7%. Furthermore, we investigate the unique ability of CEBGAN to produce a log-likelihood lower bound that may help select generated samples of higher performance (e.g., aerodynamic performance). In addition, we provide insights into the performance differences between these two models with low-dimensional toy problems and visualizations. These results and the probabilistic formulation of this inverse problem justify the extension of our GAN-based inverse design paradigm to other inverse design problems or broader inverse problems.

Engineering↗
DOE OSTI2021

Comprehensive Material Characterization and Simultaneous Model Calibration for Improved Computational Simulation Credibility

Computational simulation is increasingly relied upon for high-consequence engineering decisions, and a foundational element to solid mechanics simulations is a credible material model. Our ultimate vision is to interlace material characterization and model calibration in a real-time feedback loop, where the current model calibration results will drive the experiment to load regimes that add the most useful information to reduce parameter uncertainty. The current work investigated one key step to this Interlaced Characterization and Calibration (ICC) paradigm, using a finite load-path tree to incorporate history/path dependency of nonlinear material models into a network of surrogate models that replace computationally-expensive finite-element analyses. Our reference simulation was an elastoplastic material point subject to biaxial deformation with a Hill anisotropic yield criterion. Training data was generated using either a space-filling or adaptive sampling method, and surrogates were built using either Gaussian process or polynomial chaos expansion methods. Surrogate error was evaluated to be on the order of 10 ⁻5 and 10 ⁻3 percent for the space-filling and adaptive sampling training data, respectively. Direct Bayesian inference was performed with the surrogate network and with the reference material point simulator, and results agreed to within 3 significant figures for the mean parameter values, with a reduction in computational cost over 5 orders of magnitude. These results bought down risk regarding the surrogate network and facilitated a successful FY22-24 full LDRD proposal to research and develop the complete ICC paradigm.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Molecular level characterization of wildfire aerosol induced changes to plant health and value

The potential impacts of wildfire smoke on plant health and value are of increasing concern, due both to increasing awareness, and more frequent wildfire events. Wildfires produce substantial amounts of atmospheric pollution in smoke and particulate matter which can travel thousands of miles. These smoke events can blanket entire agricultural regions and cause impacts to plant development. One notable example of this is ‘smoke taint’ in wine, where vines, exposed to smoke, absorb many small volatile phenolic (VP) compounds. These VPs can be metabolized (glycosylated) and transported throughout the plant. Later, during downstream processing (fermentation of the grape musts), these glycosylated VPs can be hydrolyzed to re-release and volatilize the aroma-active phenolic compounds. The final product, thus, smells of smoke. Despite the obvious commercial and general research interests here, the nature of wildfire impacts on plant health are incompletely understood. Using established metabolomics and organic matter (aerosol) characterization methods, the overall aim of this proposal is thus to study wildfire smoke impacts on plant health and value, in the context of a model system (V. vinifera). This project closed early, and this report reflects some of the limited work conducted prior to conclusions. Specifically, this report discusses some sample preparation and initial MALDI-mass spectrometry analysis.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A Bayesian model for multivariate discrete data using spatial and expert information with application to inferring building attributes

When modeling sparsely observed multivariate data, strong prior information elicited from experts can be used to bolster predictive accuracy and counteract sampling bias. Similarly, modeling autocorrelation in space can help make use of co-occurrence patterns present in many types of spatial data. To make use of both expert prior information and spatial structure, we propose a novel graphical model for a spatial Bayesian network developed specifically to address challenges in inferring the attributes of buildings from geographically sparse observational data. This model is implemented as the sum of a spatial multivariate Gaussian random field and a tabular conditional probability function in real-valued space prior to projection onto the probability simplex. This modeling form is especially suitable for the usage of prior information in the form of sets of atomic rules obtained from experts. To perform inference with missing data, we implement a Markov chain Monte Carlo scheme composed of alternating steps of Gibbs sampling of missing entries and Hamiltonian Monte Carlo for model parameters. A case study in building attribution is presented to highlight the advantages and limitations of this approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Linking Extragalactic Transients and Their Host Galaxy Properties: Transient Sample, Multiwavelength Host Identification, and Database Construction

Understanding the preferences of transient types for host galaxies with certain characteristics is key to studies of transient physics and galaxy evolution, as well as to transient identification and classification in the LSST era. Here we describe a value-added database of extragalactic transients—supernovae, tidal disruption events, gamma-ray bursts, and other rare events—and their host galaxy properties. Based on reported coordinates, redshifts, and host galaxies (if known) of events, we cross-identify their host galaxies or most likely host candidates in various value-added or survey catalogs, and compile the existing photometric, spectroscopic, and derived physical properties of the host galaxies in these catalogs. This new database covers photometric measurements from the far-ultraviolet to mid-infrared. Spectroscopic measurements and derived physical properties are also available for a smaller subset of hosts. For our 36,333 unique events, we have cross-identified 13,753 host galaxies using host names, plus 4480 using host coordinates. Besides those with known hosts, there are 18,100 transients with newly identified host candidates. This large database will allow explorations of the connections of transients to their hosts, including a path toward transient alert filtering and probabilistic classification based on host properties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Galaxy Clusters Discovered via the Thermal Sunyaev-Zel’dovich Effect in the 500-square-degree SPTpol Survey

We present a catalog of 689 galaxy cluster candidates detected at significance ξ > 4 via their thermal Sunyaev-Zel’dovich (SZ) effect signature in 95 and 150 GHz data from the 500-square-degree SPTpol survey. We use optical and infrared data from the Dark Energy Camera and the Wide-field Infrared Survey Explorer (WISE) and Spitzer satellites, to confirm 544 of these candidates as clusters with ~ 94% purity. The sample has an approximately redshift-independent mass threshold at redshift z > 0.25. The confirmed sample spans 1.5 × 10 14 < M 500 c < 9 × 10 14 M ⊙ /h 70 and 0.03 < z ≲ 1.6 in mass and redshift, respectively, with a median mass of 2.5×10 14 M ⊙ /h 70 and median redshift z = 0.7; 21% of the confirmed clusters are at z > 1. We use external radio data from the Sydney University Molonglo Sky Survey (SUMSS) to estimate contamination to the SZ signal from synchrotron sources. The contamination reduces the recovered ξ by a median value of 0.032, or ~ 0.8% of the ξ = 4 threshold value, and ~ 7% of candidates have a predicted contamination greater than Δξ = 1. With the exception of a small number of systems (< 1%), an analysis of clusters detected in single-frequency 95 and 150 GHz data shows no significant contamination of the SZ signal by emission from dusty or synchrotron sources. This cluster sample, representing the deepest SZ-selected cluster sample to date, will be a key component in upcoming astrophysical and cosmological analyses of clusters. In addition to the cluster catalog, we also release the millimeter-wave maps and associated data products used to produce this sample. These maps have depths of 5.3 (11.7) µK CMB -arcmin at 150 (95) GHz and an effective angular resolution of 1. '2 (1. '7). The SPTpol products are available at https:// pole.uchicago.edu/public/data/sptpol_500d_clusters/index.html, and the NASA LAMBDA website. An interactive sky server with the SPTpol maps and Dark Energy Survey data release 2 images is also available at NCSA https://skyviewer.ncsa.illinois.edu.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Forecasts on the Dark Matter Density Profiles of Dwarf Spheroidal Galaxies with Current and Future Kinematic Observations

We forecast parameter uncertainties on the mass profile of a typical Milky Way dwarf spheroidal galaxy (dSph) using the spherical Jeans equation and Fisher matrix formalism. For a Draco-like system we show that radial velocity measurements for 1000 individual stars can constrain the mass contained within the effective radius of a dSph to within 5%. This is consistent with constraints extracted from current observational data. We compare two systems, a cusp and core, and demonstrate that a minimum sample of 100,000 (10,000) stars with both radial and proper motions measurements is required to disentangle their inner slopes at the 2σ (1σ) level. If using the log-slope measured at the half-light radius as a proxy for differentiating between a core or cusp slope, only 1000 line-of-sight and proper motions measurements are required; however, we show this choice of radius does not always unambiguously differentiate between core and cusped profiles. Once observational errors are below half the value of the intrinsic dispersion, improving the observational precision yields little change in the density profile uncertainties. The choice of priors in our profile shape analysis plays a crucial role when the number of stars in a system is less than 100 but does not affect the resulting uncertainties for larger kinematic samples. Our predicted 2D confidence regions agree well with those from a full likelihood analysis run on a mock kinematic data set taken from the Gaia Challenge, validating our Fisher predictions. Our methodology is flexible, allowing us to predict density profile uncertainties for a wide range of current and future kinematic data sets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Assessing high fidelity multi-component models to facilitate safeguards at Gas Centrifuge Enrichment Plants

We report that the International Atomic Energy Agency (IAEA) inspectors routinely carry out environmental sampling (ES) as a verification method. Collection of environmental swipe samples at various locations in Gas Centrifuge Enrichment Facilities (GCEPs) is an important process in detecting misuse of a declared facility and possibly the existence of undeclared nuclear material. These samples are measured for isotopic composition in uranium containing particles by Thermal Ionization Mass Spectrometry (TIMS) or Inductively Coupled Plasma Mass Spectrometry (ICP-MS). Even though ES is highly effective in detecting the absolute value of enrichments and their deviations from the declared values, it cannot explain the cause of those changes. Several potential explanations can serve as possibilities for particles detected above or other than the declared enrichment. These include normal and non-malicious events such as the design of the enrichment cascades, unintentional failures of the machines, or overshoot during startup of a cascade. It can also be the result of deliberate misuse by the facility operators. The primary objective of this work is to understand how these factors affect the enrichments produced by a cascade and quantify anticipated multi-isotopic concentrations for each case. The following methodology is employed to determine signatures at a particular facility. 1) Utilize a new two-dimensional multi-component diffusion code to obtain centrifuge performance data and use that information to design and perform cascade analysis. Compare and contrast the results with previous 1-D radially averaged solutions from the Pancake code. 2) Design a GCEP cascade with the production goal of 19.75% 235 U for each set of machine data above. Investigate two cascade scenarios that include enrichment of natural uranium (NU) feed to 19.75% 235U in a single cascade compared to a two-step process of NU to 5% and then 5% to 19.75%. 3) Simulate the intentional vs. unintentional off-normal scenarios in the cascades to assess the differences in isotopic concentrations. A non-ideal squared-off cascade model developed at the University of Virginia is used to calculate flow rates and isotopic concentrations of the process gas. The analysis is performed using the Rome machine model operated at 600 m/s rotor speed. The upper and lower bounds of normal and abnormal enrichments in a typical facility are used in conjunction with ES results to understand the root causes of such observations.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2024

SPRUCE Xylem Native Embolism and Leaf Traits of Picea mariana and Larix laricina, 2019

This data set contains measurements of native embolism in branch xylem and associated branch and leaf traits from Picea mariana (Black Spruce) and Larix laricina (Tamarack) from September-October 2019 at the SPRUCE experiment (Hanson et al. 2017). Data are presented in one comma-separated (*.csv) file. Native embolism, a measurement of in-situ embolism in the xylem tissue that blocks water movement, measurements were conducted at the end of the growing season on cut branches using the hydraulic pipette method (see Peters et al. 2023 for full method). In short, branch segments were connected to hydraulic apprentice and flow rates of perfusion liquid were measured using graduated pipettes and stopwatch. After initial conductance measurements, branch segments were flushed using vacuum infiltration and hydraulic conductance was remeasured to calculate the percent loss in conductance due to embolism (PLC). Three branch segments (distal, middle, and proximal) were measured from each branch representing different diameters size classes. This dataset also contains leaf area associated with each measured branch segment, calculate leaf mass per area (LMA), sapwood specific conductivity (Ks), leaf area specific conductivity (Kl) and the sapwood area to leaf area ratio (Huber value). Measurements were made on mature trees from all ten treatment enclosures. One branch from each of five trees per species were used from each enclosure where available. Some plots do not contain five individual Larix laricina, in which case all available trees were sampled. All measurements were conducted in late September/early October 2019, at the end of growing season but before Larix laricina needle senescence. There are 10 experimental plots at SPRUCE: five temperature treatments (+0, +2.25, +4.5, +6.75, +9°C) at ambient CO2, and the same five temperature treatments at elevated CO2 (+500 ppm). These data were collected after the treatments had been running in full for three years, meaning much of the material measured was grown under treatment conditions.

Huber value↗
DOE OSTI2021

Changes in Hydrogen Concentration and Defect State Density at the Poly-Si/SiOx/c-Si Interface Due to Firing

We determined the density of defect states of poly-Si/SiOx/c-Si junctions featuring a wet chemical interfacial oxide from lifetime measurements using the MarcoPOLO model to calculate recombination and contact resistance in poly-Si/SiOx/c-Si-junctions. In samples that did not receive any hydrogen treatment, the Dit,cSi is about 2 × 1012 cm-2 eV1 before firing and rises to 3–7 × 1012 cm2 eV1 during firing at measured peak temperatures between 620 °C and 863 °C. To address the question of why AlOx/SiNy stacks in contrast to pure SiNy layers for hydrogenation during firing provides better passivation quality, we have measured the hydrogen concentrations at the poly-Si/SiOx/c-Si interface as a function of AlOx layer thickness and compared these to J0 and calculated Dit,c-Si values. We observe an increase of the hydrogen concentration at the SiOx/c-Si interface upon firing as a function of the firing temperature that exceeds the defect concentrations at the interface several times. However, the AlOx layer thickness appears to cause an increase in hydrogen concentration at the SiOx/c-Si interface in these samples rather than exhibiting a hydrogen blocking property.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2025

Ultra-high temperature testing and performance of L-PBF C103

Additive Manufacturing (AM) of refractory alloys is gaining traction as a materials processing route for components subject to extreme temperature environments. Due to the low oxidation resistance of refractory alloys, novel methods for evaluating their elevated temperature performance must be developed. In this work, a Gleeble® 3800 thermomechanical load frame was modified to evaluate the mechanical properties of laser powder bed fusion (L-PBF) consolidated niobium alloy C103 ranging from room temperature (RT) to 1400 °C. The fixturing and sample geometry were designed to accommodate Joule heating and prevent damage to the test chamber. Oxidation of the samples was minimized via testing in vacuum level of 1E-5 Torr. Ultimate tensile strength (UTS), yield strength (YS), elongation, and strain-hardening behavior were determined as a function of temperature. L-PBF C103 presented an average UTS of ∼650 MPa and over 25 % elongation at RT. Above RT, the UTS and YS dropped then leveled off from 500 °C to 1000 °C with values ranging from ∼400 MPa to ∼460 MPa, which is consistent with dynamic strain aging observed in this class of alloys. The strength rapidly declined after 1200 °C to ∼150 MPa at 1400 °C. Fractography indicated ductile fractures for the C103 at all test temperatures, and Electron Backscatter Diffraction (EBSD) analysis revealed a textured microstructure and the presence of dynamic recrystallization within the necked region of the sample tested at 1400 °C.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

IUPAC-IUGS recommendation on the half-lives of 147 Sm and 146 Sm

In this paper, the IUPAC-IUGS joint Task Group “Isotopes in Geosciences” recommends a value of (106.25 ± 0.38) Ga for the half-life of 147 Sm, and a corresponding decay constant λ 147 = (6.524 ± 0.024) × 10 –12 a –1 , both with a coverage factor k = 2. For the extinct radionuclide 146 Sm two very different half-lives are used in the scientific community ( c . 68 and 103 Ma), to such a degree that no consensus value can be endorsed at present by the Task Group. Pending dedicated re-investigations it is recommended that papers using the 146 Sm decay to quantify the cosmo/geological evolution of (extra)terrestrial samples perform a twin set of calculations using both proposed half-lives.

146Sm↗
DOE OSTI2021

Inorganic characterization of switchgrass biomass using laser-induced breakdown spectroscopy

The inorganic characterization of 74 samples of switchgrass using laser-induced breakdown spectroscopy (LIBS) was undertaken. Determination of ash and inorganic elements content in biomass materials is vital for feedstock screening for bioconversion processes. Hierarchical models using principal component analysis (PCA) and partial least square analysis (PLS) were used to determine the presence of specific elemental micronutrients that are important in determining plant health for robust biomass production. LIBS uses a 532 nm laser with 45 mJ of laser power to excite the samples of switchgrass plant material and the emission of all the elements present in the plant samples were recorded in single spectra with a wide wavelength range of 200–800 nm. The results were compared to the laboratory standard technique, e.g., ICP-OES technique, to determine the true values for major micronutrients such as, silicon (Si), potassium (K), calcium (Ca), magnesium (Mg), phosphorus (P), and sulfur (S). Overall, our objectives were: 1) To determine the spectral features of switchgrass containing different amounts of these elements and 2) To examine the viability of this technique for determining the quality of the feedstock in terms of its inorganic composition. Cross-validation results showed that the broad-based model developed is promising for inorganics prediction in switchgrass. The LIBS validation prediction for the micronutrient elements mentioned here have been obtained. The regression coefficients for Si, were obtained to be 0.995, 0.994 for calibration and validation respectively, in case of Ca the regression coefficients were, 0.994 and 0.992 for calibration and validation. Similarly, in the case of Mg and K these were calculated to be 0.992 and 0.985, and 0.994 and 0.993 respectively. The regression coefficients are not as good as those for the elements mentioned, in case of the two elements S and P. They are 0.957, and 0.878, and 0.952 and 0.894 respectively for calibration, validation for the two elements. This demonstrates that LIBS-based techniques are inherently well suited for diverse environmental applications. Furthermore, LIBS along with PLS model can show capability in determining the viability of switchgrass as a biomass in the production of biofuels and survivability of switchgrass in processes associated with climate change. LIBS can help determining which switchgrass would be appropriate for a specific conversion process that favors low ash content overall or low value of specific inorganics.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Emulator-based Bayesian calibration of a subglacial drainage model

Subglacial drainage models, often motivated by the relationship between hydrology and ice flow, sensitively depend on numerous unconstrained parameters. We explore using borehole water-pressure time series to calibrate the uncertain parameters of a popular subglacial drainage model, taking a Bayesian perspective to quantify the uncertainty in parameter estimates and in the calibrated model predictions. To reduce the computation time associated with Markov Chain Monte Carlo sampling, we construct a fast Gaussian process emulator to stand in for the subglacial drainage model. We first carry out a calibration experiment using synthetic observations consisting of model simulations with hidden parameter values as a demonstration of the method. Using real borehole water pressures measured in western Greenland, we find meaningful constraints on four of the eight model parameters and a factor-of-three reduction in uncertainty of the calibrated model predictions. These experiments illustrate Gaussian process-based Bayesian inference as a useful tool for calibration and uncertainty quantification of complex glaciological models using field data. However, significant differences between the calibrated model and the borehole data suggest that structural limitations of the model, rather than poorly constrained parameters or computational cost, remain the most important constraint on subglacial drainage modelling.

58 GEOSCIENCES↗