Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sampling 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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Modeling PAH Mixture Interactions in a Human In Vitro Organotypic Respiratory Model

One of the most significant challenges in human health risk assessment is to evaluate hazards from exposure to environmental chemical mixtures. Polycyclic aromatic hydrocarbons (PAHs) are a class of ubiquitous contaminants typically found as mixtures in gaseous and particulate phases in ambient air pollution associated with petrochemicals from Superfund sites and the burning of fossil fuels. However, little is understood about how PAHs in mixtures contribute to toxicity in lung cells. To investigate mixture interactions and component additivity from environmentally relevant PAHs, two synthetic mixtures were created from PAHs identified in passive air samplers at a legacy creosote site impacted by wildfires. The primary human bronchial epithelial cells differentiated at the air–liquid interface were treated with PAH mixtures at environmentally relevant proportions and evaluated for the differential expression of transcriptional biomarkers related to xenobiotic metabolism, oxidative stress response, barrier integrity, and DNA damage response. Component additivity was evaluated across all endpoints using two independent action (IA) models with and without the scaling of components by toxic equivalence factors. Both IA models exhibited trends that were unlike the observed mixture response and generally underestimated the toxicity across dose suggesting the potential for non-additive interactions of components. Overall, this study provides an example of the usefulness of mixture toxicity assessment with the currently available methods while demonstrating the need for more complex yet interpretable mixture response evaluation methods for environmental samples.

3D in vitro models↗
DOE OSTI2021

Using Multigrain Crystallography to Explore the Microstructural Evolution of the α-Olivine to γ-Ringwoodite Transformation and ε-Mg2SiO4 at High Pressure and Temperature

The introduction of multigrain crystallography (MGC) applied in a laser-heated diamond anvil cell (LH-DAC) using synchrotron X-rays has provided a new path to investigate the microstructural evolution of materials at extreme conditions, allowing for simultaneous investigations of phase identification, strain state determination, and orientation relations across phase transitions in a single experiment. Here, we applied this method to a sample of San Carlos olivine beginning at ambient conditions and through the α-olivine → γ-ringwoodite phase transition. At ambient temperatures, by measuring the evolution of individual Bragg reflections, olivine shows profuse angular streaking consistent with the onset of yielding at a measured stress of ~1.5 GPa, considerably lower than previously reported, which may have implications for mantle evolution. Furthermore, γ-ringwoodite phase was found to nucleate as micron to sub-micron grains imbedded with small amounts of a secondary phase at 15 GPa and 1000 °C. Using MGC, we were able to extract and refine individual crystallites of the secondary unknown phase where it was found to have a structure consistent with the ε-phase previously described in chondritic meteorites.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Light Element Abundances and Multiple Populations in M53

We present results from a study of 94 red giant stars in the globular cluster M53. We use low-resolution spectra to measure the strength of CN and CH features at ∼3800 and 4300 Å, respectively. The strengths of these features are used to classify stars into a CN-enhanced and CN-normal population and to measure C and N abundances in all 94 stars. We find the red giant branch stars to be evenly split between the two populations identified, and observe the presence of CN-enhanced stars on the asymptotic giant branch. In addition, we identify five CH star candidates based on the strength of their CN and CH band features, and the presence of a P-branch in their CH band. We compare our identification of multiple populations to those based on the Na–O anticorrelation and pseudo-color indices in Hubble Space Telescope UV photometry, and find general agreement between all three methods. Our large sample size also allows us to study the radial distribution of each population, and we find that the CN-enhanced population is more centrally concentrated. We use our C and N measurements to compare the evolutionary changes in these elements as a function of magnitude between the two populations, and show that both populations experience similar evolutionary changes to the surface abundances of C and N. Finally, we calculate C+N+O abundances for each population and compare them to similar measurements made in M10; we find that in both clusters, CN-enhanced stars have a slightly enhanced C+N+O (Δ(C+N+O) ∼ 0.2 dex).

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Early Ultraviolet Observations of Type IIn Supernovae Constrain the Asphericity of Their Circumstellar Material

Here we present a survey of the early evolution of 12 Type IIn supernovae (SNe IIn) at ultraviolet and visible light wavelengths. We use this survey to constrain the geometry of the circumstellar material (CSM) surrounding SN IIn explosions, which may shed light on their progenitor diversity. In order to distinguish between aspherical and spherical CSM, we estimate the blackbody radius temporal evolution of the SNe IIn of our sample, following the method introduced by Soumagnac et al. We find that higher-luminosity objects tend to show evidence for aspherical CSM. Depending on whether this correlation is due to physical reasons or to some selection bias, we derive a lower limit between 35% and 66% for the fraction of SNe IIn showing evidence for aspherical CSM. This result suggests that asphericity of the CSM surrounding SNe IIn is common - consistent with data from resolved images of stars undergoing considerable mass loss. It should be taken into account for more realistic modeling of these events.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Real-time Ridesharing for Transportation Hubs with Demand and Supply Uncertainty

Transportation hubs in major cities generate a significant amount of trips by taxis and for-hail vehicles (FHV), with many of the trips sharing similar destinations. This suggests promising opportunities to leverage the collective travel needs with dedicated ridesharing solutions to reduce the externalities of excessive traffic at transportation hubs. In this study, we develop a novel dynamic ridesharing approach to serve trips from the transportation hub by considering (1) demand (new passengers) and supply (newly available vehicles) in the near future and (2) the uncertainty of future predictions. Our approach consists of two stages. In the first stage, we develop a data structure called hub mobility tree to generate potential combinations of shareable trips as candidate schedules efficiently. Then the generated schedules are used in the second stage to formulate the stochastic hub-based ridesharing problem (SHRP), which is a stochastic integer programming problem with the objective to maximize the total expected ridesharing profit over time. Due to the prohibitive number of shareable trips, we then approximately solve SHRP by the sample average approximate method (SAA), and a dual Lagrangian technique is implemented to further improve the scalability of the solution approach. We demonstrate the performances of the proposed method by simulating the ridesharing service at JFK airport using NYC taxi and FHV data. The results indicate that the proposed method outperforms the myopic ridesharing (maximize profit for a single time step) and the rolling horizon method with point estimation of future demand and supply.

dynamic ridesharing↗
DOE OSTI2020

The Effects of Contaminant-Aging on Decontamination Efficacy for Rapid Remediation of Concrete Surfaces - 20412

Understanding the relationship between decontamination efficacy and contaminant aging is imperative for developing effective remediation strategies following a large-scale contamination event. Non-destructive decontamination of concrete following a contamination incident becomes more difficult with increased time between contamination and decontamination because contaminants may chemically bind to concrete and/or penetrate into the concrete subsurface. In this work, we evaluated the decontamination efficacy of two decontamination methods on concrete samples contaminated with soluble Cesium-137 (Cs-137) and silica particles with 0.5 μm and 2 μm diameters. Concrete samples were aged between 1 and 59 days, with half of the coupons receiving 1 mL of artificial rainwater about once every three days. After aging, coupons were either decontaminated or analyzed to determine contaminant penetration depths. Coupons were decontaminated with 0.1 M Potassium Chloride (KCI) solution applied either by flowing solution across the contaminated coupon face or non-destructive power washing. Contaminant depth profiles were created by removing the top surface of the coupon at least 15 times, measuring the activity in each removed layer, normalizing the total activity removed to the measured decontamination efficacy following the last layer removal, and determining the penetration thickness using the coupon dimensions, bulk density, and mass of material removed for each layer. Pressurized washing effectively removed particles, but not soluble cesium, despite the majority of Cs-137 being located within the first millimeter of the subsurface. Soluble cesium removal by flowing 0.1 M KCI solution across the contaminated surface dropped to less than 12% within the first ten days. This drop in Cs-137 removal efficacy occurred within the first five days for coupons experiencing artificial rainfall events. Flowing 0.1 M KCI solution across the contaminated face was ineffective at removing particles, likely because particles had localized in surface depressions. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Rapid prediction of fuel research octane number and octane sensitivity using the AFIDA constant-volume combustion chamber

Current research octane number (RON) and motor octane number (MON) gasoline performance characterization techniques use dated, complex engine testing methodology and limit researchers’ ability to easily characterize small volumes of experimental fuels. A novel methodology is presented that correlates measured ignition delay (ID) time to RON in an Advanced Fuel Ignition Delay Analyzer (AFIDA) constant-volume combustion chamber device at a single pressure/temperature condition, with an r2 of 0.99 and standard error (SE) of 1.0. The correlation of the slope of the ID time between two additional temperature points to octane sensitivity (S) produces an r2 of 0.97 and SE of 0.69; however, fuels with S>12 are indistinguishable. These results are based on methodology calibration using 31 primary and toluene reference fuels containing 0%-40% ethanol with RON values ranging from 85 to 113. Validation of these methods using a 102-sample fuel matrix spanning an array of base fuels and additive chemistry designed to test the robust applicability of the method, along with pump gasoline and high-octane surrogate blend samples, demonstrates an r2 of 0.94 and SE of 1.3 for the RON correlation over all samples, whereas the equivalent S correlation produces an r2 of 0.78 and SE of 1.2 by excluding two additives, 3-pentanone and diisobutylene, which displayed poor S correlation results. This novel AFIDA analysis method can be performed in 1 h and with 40 mL of fuel, offering significant improvements in time and volume requirements over traditional techniques.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Primal-Dual Differentiable Programming for Distribution System Critical Load Restoration: Preprint

Swift and reliable critical load restoration (CLR) can help make a distribution system resilient towards extreme events. To optimally achieve that, alongside practical concerns such as limiting online computational burden, some studies leverage model-free reinforcement learning (RL) to train control policies. Despite the advantages provided by RL algorithms, these approaches suffer from two issues: 1) the lack of a proper mechanism for constraint enforcement, and 2) poor sample efficiency. Therefore, in this paper, a primal-dual differentiable programming (PDDP) method is developed for guiding the training leading to a constraint-satisfying policy. Additionally, the model-based nature of the proposed method aims at improving sample efficiency. The experiment on a CLR problem demonstrates that PDDP can effectively train a control policy that both achieves desirable performance and satisfies required constraints.

differentiable programming↗
DOE OSTI2025

Analysis of semivolatile organics in liquid radioactive residue using sorptive stir bar and solvent back-extraction

Current methods for semivolatiles analysis in radioactive samples can produce large volumes of radioactive solvent residue. A method utilizing stir-bar sorptive extraction has been explored in this work for its applicability to radioactive waste samples. This low solvent analytical method may accelerate remediation, minimize hazardous solvent waste, and reduce exposure risk to workers. Organic compounds (polyaromatic hydrocarbons, chlorinated aromatics, and phenolic compounds) were chosen as surrogates for common Liquid Waste System (LWS) contaminants at Savannah River Site (Aiken, SC). Stir-bar extraction parameters (extraction time, matrix modification, and effective pH range) and solvent back extraction parameters (solvent type, volume, and extraction time) were optimized experimentally for the chosen compounds. Affinity of the stir-bar extraction polymer to radionuclides Cs-137 and Am-241 was observed to determine radionuclide concentration effects. The stir-bar method achieved mean recovery of 100 ± 0.7% (1σ), relative to 114 ± 7% using solvent extraction, while reducing weekly method hands-on time by 93.4% and solvent volume consumption by 99.3%. Sensitivity was improved by 378% in simulated tank waste and 278% in real-world LWS matrix, relative to solvent extraction. This work has produced a safe and optimized method for the low solvent analysis of organics in legacy radioactive tank waste by stir-bar sorptive extraction.

GC-MS↗
DOE OSTI2020

Modern Monte Carlo methods for efficient uncertainty quantification and propagation: A survey

Uncertainty quantification (UQ) includes the characterization, integration, and propagation of uncertainties that result from stochastic variations and a lack of knowledge or data in the natural world. Monte Carlo (MC) method is a sampling-based approach that has widely used for quantification and propagation of uncertainties. However, the standard MC method is often time-consuming if the simulation-based model is computationally intensive. This article gives an overview of modern MC methods to address the existing challenges of the standard MC in the context of UQ. Specifically, multilevel Monte Carlo (MLMC) extending the concept of control variates achieves a significant reduction of the computational cost by performing most evaluations with low accuracy and corresponding low cost, and relatively few evaluations at high accuracy and corresponding high cost. Multifidelity Monte Carlo (MFMC) accelerates the convergence of standard Monte Carlo by generalizing the control variates with different models having varying fidelities and varying computational costs. Multimodel Monte Carlo method (MMMC), having a different setting of MLMC and MFMC, aims to address the issue of UQ and propagation when data for characterizing probability distributions are limited. Multimodel inference combined with importance sampling is proposed for quantifying and efficiently propagating the uncertainties resulting from small data sets. All of these three modern MC methods achieve a significant improvement of computational efficiency for probabilistic UQ, particularly uncertainty propagation. Finally, an algorithm summary and the corresponding code implementation are provided for each of the modern MC methods. The extension and application of these methods are discussed in detail.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Randomized Algorithms for Symmetric Nonnegative Matrix Factorization

Symmetric Nonnegative Matrix Factorization (SymNMF) is a technique in data analysis and machine learning that approximates a matrix with a product of a nonnegative, low-rank matrix and it transpose. To design faster and more scalable algorithms for SymNMF we develop two randomized algorithms for its computation. The first method uses randomized matrix sketching to compute an initial low-rank approximation to the input matrix and proceeds to uses this as a low-rank input to rapidly compute a SymNMF. The second methods uses randomized leverage score sampling to approximately solve constrained least squares problems. Many successful methods for SymNMF rely on (approximately) solving sequences of constrained least squares problems. Here, we prove theoretically that leverage score sampling can approximately solve constrained least squares problems to e-accuracy. Finally we demonstrate both methods work in practice by applying them to graph clustering tasks on large real world data sets. These experiments show that our methods approximately maintain solution quality and achieve significant speed ups for both large dense and large sparse problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Binning is Sinning: Redemption for Hubble Diagram Using Photometrically Classified Type Ia Supernovae

Bayesian Estimation Applied to Multiple Species (BEAMS) is implemented in the BEAMS with Bias Corrections (BBC) framework to produce a redshift-binned Hubble diagram (HD) for Type Ia supernovae (SNe Ia). BBC corrects for selection effects and non–SN Ia contamination, and systematic uncertainties are described by a covariance matrix with dimension matching the number of BBC redshift bins. For spectroscopically confirmed SN Ia samples, a recent "Binning is Sinning" article showed that an unbinned HD and covariance matrix reduces the systematic uncertainty by a factor of ~1.5 compared to the binned approach. Here we extend their analysis to obtain an unbinned HD for a photometrically identified sample processed with BBC. To test this new method, we simulate and analyze 50 samples corresponding to the Dark Energy Survey (DES) with a low-redshift anchor; the simulation includes SNe Ia, and contaminants from core-collapse SNe and peculiar SNe Ia. The analysis includes systematic uncertainties for calibration and measures the dark energy equation of state parameter (w). Compared to a redshift-binned HD, the unbinned HD with nearly 2000 events results in a smaller systematic uncertainty, in qualitative agreement with BHS21, and averaging results among the 50 samples we find no evidence for a w-bias. To reduce computation time for fitting an unbinned HD with large samples, we propose an HD-rebinning method that defines the HD in bins of redshift, color, and stretch; the rebinned HD results in similar uncertainty as the unbinned case, and shows no evidence for a w-bias.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Development of the uncertainty quantification toolkit's python interface and surrogate construction tutorial

The uncertainty quantification toolkit (UQTk) is a collection of c++ libraries that assess the confidence of numerical models. Surrogate approximations, often polynomial chaos expansions (PCEs), lessen the computational cost of these assessments. I developed a Python interface in UQTk for regression and Bayesian compressive sensing to add to the existing Galerkin projection method. These methods receive an object containing the polynomial basis information and NumPy arrays of sample points, call c++ methods, and return the PCE coefficients in a NumPy array. To demonstrate these methods, I wrote a tutorial in which I use them to construct surrogates for Genz functions and calculate the resulting error.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Variance Reduction within Implicit Monte Carlo Thermal Radiation Transport using the Local Importance Function Transform [Slides]

System of equations derived by Fleck and Cummings to model time-dependent thermal radiative transfer (TRT) problems: Solution of radiation specific intensity (I), non-linearly coupled to material internal energy (U); Absorption and re-emission approximated by effective scattering events. Useful for high-energy density physics simulations: i.e., astrophysics, inertial confinement fusion (ICF). Whereas deterministic methods are fully discretized and free of statistical noise, Monte Carlo methods allow for dynamic sampling of the phase space at the price of statistical noise; Monte Carlo methods are also characterized by a slow $\frac{1}{√Ν}$ convergence rate.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Temporal Study 2021-2022: Sample-Based Surface Water Chemistry and Organic Matter Characterization across Watersheds in the Yakima River Basin, Washington, USA (v3)

This dataset supports a broader study examining the drivers of temporal variability in sediment respiration rates in the Yakima River Basin. The dataset provides geochemistry and organic matter characterization data generated from samples collected at weekly or bi-weekly intervals at six sites across the Yakima River Basin in Washington, USA. Related sensor data will be published separately and can be used to link sediment respiration rates to biogeochemical processing rates. The data package was originally published in November 2022. It was updated in November 2023 (v2; modified files) and April 2025 (new and modified files). See the change history section in readme for more details. For details on how to navigate data packages generated by this project, see https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About. This dataset is comprised of one main data folder containing (1) file-level metadata; (2) data dictionary; (3) field metadata; (4) dissolved inorganic carbon (DIC), dissolved organic carbon (DOC; reported as non-purgeable organic carbon; NPOC), total nitrogen (TN), total suspended solids (TSS), and ions; (5) averaged values from water chemistry data; (6) surface water sampling protocol; (7) sensor protocol; (8) readme; (9) methods codes; (10) international generic sample number (IGSN) mapping file; and (11) folder of high resolution characterization of organic via 12 Tesla Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS) through the Environmental Molecular Sciences Laboratory (EMSL; https://www.pnnl.gov/environmental-molecular-sciences-laboratory). This folder contains two subfolders, one containing the .xml data files and the other containing instructions for using Formularity (https://omics.pnl.gov/software/formularity) and an R script to process the data based on the user's specific needs. All files are .csv, .pdf, .R, .ref, or .xml. We acknowledge the Yakama Nation as owners and caretakers of the lands where we collected our data. We thank the Confederated Tribes and Bands of the Yakama Nation Tribal Council and Yakama Nation Fisheries for working with us to facilitate sample collection and optimization of data usage according to their values and worldview.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Spontaneous Imbibition Tests and Parameter Estimation in Volcanic Tuff

We present a dynamic laboratory spontaneous imbibition test and interpretation method, demonstrated on volcanic tuff samples from the Nevada National Security Site. The method includes numerical inverse modeling to quantify uncertainty of estimated two-phase fluid flow properties. As opposed to other approaches requiring multiple different laboratory instruments, the dynamic imbibition method simultaneously estimates capillary pressure and relative permeability from one test apparatus.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Improving cellulose attribution by selectively removing yeast glucans from grain fermentation intermediates

Abstract This article presents an industry-relevant method for quantifying cellulose in mixed substrate samples. We built upon the cellulosic glucan measurement proposed by Sluiter et al. (2021) to investigate significant cellulose loss under cold caustic conditions, which has hindered the establishment of an industrially relevant method. To overcome this issue, we used dimethyl sulfoxide (DMSO) as an alternative solvent, which avoids mercerization reactions (Budtova and Navard 2015) but may leave some resistant starch in the sample. Treatment via an enzyme mixture removed starch and conformed to US Environmental Protection Agency (EPA) guidance on cellulose measurement by establishing quantitative de-starching via nuclear magnetic resonance (EPA 2022). The selective removal of yeast beta-glucan was accomplished using a commercially available Zymolyase. Our results demonstrated excellent reproducibility, with coefficients of variance of 7.14% or less, when measuring cellulose in low cellulose content samples. The method was tested on relevant lab and plant samples, showing an average 0.9% ethanol derived from the conversion of cellulose when cellulase was added to the process and zero response when no cellulase was added. These findings indicate that a cellulose method using DMSO and pullulanase provides a complete and accurate view of cellulose content, composition, and conversion in industrial fermentation processes.

Sluiter, Justin [National Renewable Energy Laborat↗
DOE OSTI2023

Importance of nonuniform Brillouin zone sampling for ab initio Bethe-Salpeter equation calculations of exciton binding energies in crystalline solids

Excitons are prevalent in semiconductors and insulators, and their binding energies are critical for optoelectronic applications. The state-of-the-art method for first-principles calculations of excitons in extended systems is the ab initio GW-Bethe-Salpeter equation (BSE) approach, which can require a fine sampling of reciprocal space to accurately resolve solid-state exciton properties. Here, in this study, we show, for a range of semiconductors and insulators, that the commonly employed approach of uniformly sampling the Brillouin zone can lead to underconverged exciton binding energies, as impractical grid sizes are required to achieve adequate convergence. We further show that nonuniform sampling of the Brillouin zone, focused on the region of reciprocal space where the exciton wave function resides, enables efficient rapid numerical convergence of exciton binding energies at a given level of theory. We propose a well-defined convergence procedure, which can be carried out at relatively low computational cost and which in some cases leads to a correction of previous best theoretical estimates by almost a factor of 2, qualitatively changing the predicted exciton physics. These results call for the adoption of nonuniform sampling methods for ab initio GW-BSE calculations and for revisiting previously computed values for exciton binding energies of many systems.

36 MATERIALS SCIENCE↗