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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Standardized Data to Improve Understanding and Modeling of Soil Nitrogen at Continental Scale

Abstract Nitrogen (N) is a key limiting nutrient in terrestrial ecosystems, but there remain critical gaps in our ability to predict and model controls on soil N cycling. This may be in part due to lack of standardized sampling across broad spatial–temporal scales. Here, we introduce a continentally distributed, publicly available data set collected by the National Ecological Observatory Network (NEON) that can help fill these gaps. First, we detail the sampling design and methods used to collect and analyze soil inorganic N pool and net flux rate data from 47 terrestrial sites. We address methodological challenges in generating a standardized data set, even for a network using uniform protocols. Then, we evaluate sources of variation within the sampling design and compare measured net N mineralization to simulated fluxes from the Community Earth System Model 2 (CESM2). We observed wide spatiotemporal variation in inorganic N pool sizes and net transformation rates. Site explained the most variation in NEON’s stratified sampling design, followed by plots within sites. Organic horizons had larger pools and net N transformation rates than mineral horizons on a sample weight basis. The majority of sites showed some degree of seasonality in N dynamics, but overall these temporal patterns were not matched by CESM2, leading to poor correspondence between observed and modeled data. Looking forward, these data can reveal new insights into controls on soil N cycling, especially in the context of other environmental data sets provided by NEON, and should be leveraged to improve predictive modeling of the soil N cycle.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Psilocybin: crystal structure solutions enable phase analysis of prior art and recently patented examples

Psilocybin {systematic name: 3-[2-(dimethylamino)ethyl]-1 H -indol-4-yl dihydrogen phosphate} is a zwitterionic tryptamine natural product found in numerous species of fungi known for their psychoactive properties. Following its structural elucidation and chemical synthesis in 1959, purified synthetic psilocybin has been evaluated in clinical trials and has shown promise in the treatment of various mental health disorders. In a recent process-scale crystallization investigation, three crystalline forms of psilocybin were repeatedly observed: Hydrate A, Polymorph A, and Polymorph B. The crystal structure for Hydrate A was solved previously by single-crystal X-ray diffraction. This article presents new crystal structure solutions for the two anhydrates, Polymorphs A and B, based on Rietveld refinement using laboratory and synchrotron X-ray diffraction data, and density functional theory (DFT) calculations. Utilizing the three solved structures, an investigation was conducted via Rietveld method (RM) based quantitative phase analysis (QPA) to estimate the contribution of the three different forms in powder X-ray diffraction (PXRD) patterns provided by different sources of bulk psilocybin produced between 1963 and 2021. Over the last 57 years, each of these samples quantitatively reflect one or more of the hydrate and anhydrate polymorphs. In addition to quantitatively evaluating the composition of each sample, this article evaluates correlations between the crystal forms present, corresponding process methods, sample age, and storage conditions. Furthermore, revision is recommended on characterizations in recently granted patents that include descriptions of crystalline psilocybin inappropriately reported as a single-phase `isostructural variant.' Rietveld refinement demonstrated that the claimed material was composed of approximately 81% Polymorph A and 19% Polymorph B, both of which have been identified in historical samples. In this article, we show conclusively that all published data can be explained in terms of three well-defined forms of psilocybin and that no additional forms are needed to explain the diffraction patterns.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2020

Multi Model Monte Carlo with Python (MXMCPy)

Multi Model Monte Carlo with Python (\mxmc {}) is a software package developed as a general capability for computing the statistics of outputs from an expensive, high-fidelity model by leveraging faster, low-fidelity models for speedup. Motivated by uncertainty propagation problems where classical Monte Carlo (MC) simulation is computationally intractable, various multi-model MC approaches have recently emerged that yield unbiased estimators with significantly reduced variance relative to MC for the same cost. These existing methods include multi-level Monte Carlo (MLMC), multi-fidelity Monte Carlo (MFMC), and approximate control variates (ACV). Given a fixed computational budget and a collection of models with varying cost/accuracy, each method seeks a sample allocation strategy across the models that results in an estimator with optimal variance reduction. \mxmc {} is a versatile tool that enables convenient access to many existing multi-model MC approaches within one modular and extensible package. With \mxmc {}, users can easily compare existing methods to determine the best choice for their particular problem, while developers have a basis for implementing and sharing new variance reduction approaches. This report introduces the \mxmc {} software, providing a summary of the problem-solving workflow for users as well as a brief overview of the code layout for developers.

Geoffrey F Bomarito↗
DOE OSTI2025

Site-specific plan-view (S)TEM sample preparation from thin films using a dual-beam FIB-SEM

To fully evaluate the atomic structure, and associated properties of materials using transmission electron microscopy, examination of samples from three non-collinear orientations is needed. This is particularly challenging for thin films and nanoscale devices built on substrates due to limitations with plan-view sample preparation. In this work, a new method for preparation of high-quality, site-specific, plan-view TEM samples from thin-films grown on substrates, is presented and discussed. Here, it is based on using a dual-beam focused ion beam scanning electron microscope (FIB-SEM) system. To demonstrate the method, the samples were prepared from thin films of perovskite oxide BaSnO 3 grown on a SrTiO 3 substrate and metal oxide IrO 2 on a TiO 2 substrate, ranging from 20–80 nm in thicknesses using molecular beam epitaxy. While the method is optimized for the thin films, it can be extended to other site-specific plan-view samples and devices build on wafers. Aberration-corrected STEM was used to evaluate the quality of the samples and their applicability for atomic-resolution imaging and analysis.

BaSnO3↗
DOE OSTI2024

Complete resolution across the neodymium/samarium isotopic envelope with a liquid sampling‐atmospheric pressure glow discharge — Orbitrap mass spectrometer

Rationale Nd and Sm isotope ratios play an important role in geological dating and as nuclear forensic signatures; however, the overlap of the respective 144, 148, 150 Nd/Sm isobars requires prior separations to be performed before analysis on typical MS platforms. The work presented here overcomes these isobaric interferences using ultrahigh‐mass resolution to alleviate interference without prior chemical separations. Methods A liquid sampling‐atmospheric pressure glow discharge ion source was coupled to a standard, QExactive Focus Orbitrap mass spectrometer, providing a mass resolution of ~80 k. A Spectroswiss FTMS booster X2 data acquisition package was used to collect extended transients, providing much higher mass resolution; ~230 k and ~600 k are employed here for Nd and Sm isotopes. Results While the standard Orbitrap resolution is far greater than typical “atomic” MS platforms, it was insufficient to alleviate all isobars. The use of a resolution of ~230 k resulted in baseline separation across the entire isotopic envelope for both Nd and Sm. Isotope ratios obtained from Nd:Sm mixtures using high‐resolution were equivalent to those found for individual‐element solutions, while isotope ratios obtained at a resolution of ~80 k (standard for the OEM data system) showed large deviations. Conclusions Use of ultrahigh‐resolution is an attractive alternative to extensive chemical separations to alleviate severe isobaric interferences. Sufficient mass resolution greatly reduces/eliminates the need for sample manipulations (separations) before analysis while reducing costs and total analysis times.

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

Matrix-Assisted Ionization of Molecular Uranium Species

Matrix-assisted ionization (MAI) demonstrates high sensitivity for a variety of organic compounds; however, few studies have reported the application of MAI for the detection and characterization of inorganic analytes. Trace-level uranium analysis is important in the realms of nuclear forensics, nuclear safeguards, and environmental monitoring. Traditional mass spectrometry methods employed in these fields require combinations of extensive laboratory chemistry sample preparation and destructive ionization methods. There has been recent interest in exploring ambient mass spectrometry methods that enable timely sample analysis and higher sensitivity than what is attainable by field-portable radiation detectors. Rapid characterization of uranium at nanogram levels is demonstrated in this study using MAI techniques. Mass spectra were collected on an atmospheric pressure mass spectrometer for solutions of uranyl nitrate, uranyl chloride, uranyl acetate, and uranyl oxalate utilizing 3-nibrobenzonitrile as the ionization matrix. The uranyl complexes investigated were detectable, and the chemical speciation was preserved. Sample analysis was accomplished in a matter of seconds, and limits of detection of 5 ng of uranyl nitrate, 10 ng of uranyl oxalate, 100 ng of uranyl chloride, and 200 ng of uranyl acetate were achieved. The observed gas-phase speciation was similar to negative-ion electrospray ionization of uranyl compounds with notable differences. Six matrix-derived ions were detected in all negative-ion mass spectra, and some of these ions formed adducts with the uranyl analyte. Subsequent analysis of the matrix suggests that these molecules are not matrix contaminants and are instead created during the ionization process.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2022

U–Pb Geochronology and Stable Isotope Geochemistry of Terrestrial Carbonates, Lower Cretaceous Cedar Mountain Formation, Utah: Implications for Synchronicity of Terrestrial and Marine Carbon Isotope Excursions

The terrestrial Lower Cretaceous Cedar Mountain Formation, Utah, is a critical archive of paleoclimate, tectonics, and vertebrate ecology and evolution. Early Cretaceous carbon cycle perturbations associated with ocean anoxia have been interpreted from this succession, as expressed in stable carbon isotopes. However, refining the timing of the observed stable isotope excursions remains a key challenge in understanding how marine anoxia affects the Earth system, and is ultimately recorded in the terrestrial realm. The geochronology and geochemistry of a terrestrial carbonate near the base of this succession, which potentially records the Ap7 global carbon isotope excursion, is studied here. Petrographic and geochemical analyses are used to test plausible mechanisms for U incorporation into the calcite lattice in this sample. Using these methods, the hypothesis that the incorporation of U was at or close to the timing of carbonate precipitation is evaluated. U–Pb geochronology of calcite indicates a plausible Early Cretaceous age. However, comparison of the new U–Pb ages of calcite with detrital zircon maximum depositional ages immediately beneath the studied sample indicates a disparity in the apparent sedimentation rates if both types of geochronologic information are interpreted as reflecting the timing of sediment deposition. The totality of data supports an early, and high-temperature, diagenetic timing of U incorporation, with potential for minor leaching of U in subsequent fluid–rock interaction. The most likely mechanism for U transport and immobilization in these samples is hydrothermal fluid–rock interaction. Therefore, the radiometric ages, and corresponding stable isotope composition of U-bearing carbonate domains in this sample, indicate early subsurface fluid–rock interactions and not a record of atmosphere–soil geochemical reactions.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

Application of a Bayesian Framework for Plasticity Model Selection

Interpretable Machine Learning (IML) has performed well when tasked with deriving constitutive material models. However, IML has been shown to prefer models that overfit noise in data, which tends to lead to bloat and a decrease in interpretability. Due to these issues, the ability of IML to reliably derive models that fit the data and are both interpretable and generalizable is limited. A method developed recently has shown promise to improve upon traditional IML by using a Bayesian fitness definition for the evolution of free-form models with non-deterministic parameters. This framework was developed for genetic-programming-based symbolic regression(GPSR) and involves model parameter estimation using Sequential Monte Carlo sampling (SMC).The method has demonstrated a reduction in bloat when dealing with noisy data in comparison to conventional GPSR. The results of this framework applied to stress-strain data for copper show models that more effectively predict the experimental data better than was previously shown with GPSR.

plasticity↗
DOE OSTI2024

Variational inference of effective range parameters for 3 He− 4 He scattering

We use two different methods, Monte Carlo sampling and variational inference (VI), to perform a Bayesian calibration of the effective-range parameters in 3 He– 4 He elastic scattering. The parameters are calibrated to data from a recent set of 3 He– 4 He elastic scattering differential cross section measurements. Analysis of these data for E lab ≤ 4.3 MeV yields a unimodal posterior for which both methods obtain the same structure. However, the effective-range expansion amplitude does not account for the 7/2 − state of 7 Be so, even after calibration, the description of data at the upper end of this energy range is poor. The data up to E lab = 2.6 MeV can be well described, but calibration to this lower-energy subset of the data yields a bimodal posterior. After adapting VI to treat such a multi-modal posterior we find good agreement between the VI results and those obtained with parallel-tempered Monte Carlo sampling.

effective field theory↗
NASA NTRS1971

Lunar rilles: A catalog and method of classification

A representative sample of lunar rilles observed on Lunar Orbiter 4 and 5 photographs is presented. A method of classification of rilles is described which will provide a means of quantitatively grouping rilles of similar planemetric shape. The method is based on a finite Fourier approximation of the edge of the rille. It provides an adequate approximation to all planemetric characteristics of the rille except rille recurvature.

Oberbeck, V. R.↗
NASA NTRS1981

A remark on the theory of measuring thermal diffusivity by the modified Angstrom's method

A theory of the measurement of the thermal diffusivity of a sample by the modified Angstrom method is developed for the case in which radiative heat loss from the end surface of the sample is not negligible, and applied to measurements performed on lunar samples. Formulas allowing sample thermal diffusivity to be determined from the amplitude decay and phase lag of a temperature wave traveling through the sample are derived for a flat disk sample for which only heat loss from the end surface is important, and a sample of finite diameter and length for which heat loss through the end and side surfaces must be considered. It is noted that in the case of a flat disk, measurements at a single angular frequency of the temperature wave are sufficient, while the sample of finite diameter and length requires measurements at two discrete angular frequencies. Comparison of the values of the thermal diffusivities of two lunar samples of dimensions approximately 1 x 1 x 2 cm derived by the present methods and by the Angstrom theory for a finite bar reveals them to differ by not more than 5%, and indicates that more refined data are required as the measurement theory becomes more complicated.

Horai, K.-I.↗
NASA NTRS2005

Improved Compression of Wavelet-Transformed Images

A recently developed data-compression method is an adaptive technique for coding quantized wavelet-transformed data, nominally as part of a complete image-data compressor. Unlike some other approaches, this method admits a simple implementation and does not rely on the use of large code tables. A common data compression approach, particularly for images, is to perform a wavelet transform on the input data, and then losslessly compress a quantized version of the wavelet-transformed data. Under this compression approach, it is common for the quantized data to include long sequences, or runs, of zeros. The new coding method uses prefixfree codes for the nonnegative integers as part of an adaptive algorithm for compressing the quantized wavelet-transformed data by run-length coding. In the form of run-length coding used here, the data sequence to be encoded is parsed into strings consisting of some number (possibly 0) of zeros, followed by a nonzero value. The nonzero value and the length of the run of zeros are encoded. For a data stream that contains a sufficiently high frequency of zeros, this method is known to be more effective than using a single variable length code to encode each symbol. The specific prefix-free codes used are from two classes of variable-length codes: a class known as Golomb codes, and a class known as exponential-Golomb codes. The codes within each class are indexed by a single integer parameter. The present method uses exponential-Golomb codes for the lengths of the runs of zeros, and Golomb codes for the nonzero values. The code parameters within each code class are determined adaptively on the fly as compression proceeds, on the basis of statistics from previously encoded values. In particular, a simple adaptive method has been devised to select the parameter identifying the particular exponential-Golomb code to use. The method tracks the average number of bits used to encode recent runlengths, and takes the difference between this average length and the code parameter. When this difference falls outside a fixed range, the code parameter is updated (increased or decreased). The Golomb code parameter is selected based on the average magnitude of recently encoded nonzero samples. The coding method requires no floating- point operations, and more readily adapts to local statistics than other methods. The method can also accommodate arbitrarily large input values and arbitrarily long runs of zeros. In practice, this means that changes in the dynamic range or size of the input data set would not require a change to the compressor. The algorithm has been tested in computational experiments on test images. A comparison with a previously developed algorithm that uses large code tables (generated via Huffman coding on training data) suggests that the data-compression effectiveness of the present algorithm is comparable to the best performance achievable by the previously developed algorithm.

Kiely, Aaron↗
DOE OSTI2022

Reconstruction of Fission Product Distribution from Tomographic Scans in TRISO Fuel Graphitic Matrix and Nuclear Grade Graphites

An image reconstruction method was developed to rectify shortcomings of earlier methods that became apparent as destructively sampled data became available. This reconstruction method was applied to the tomographic gamma scans of nuclear graphite and graphitic matrix samples from AGR-3/4. There is generally agreement between profiles measured via destructive sampling of these rings and profiles from tomographic reconstruction.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Reconstruction of Fission Product Distribution from Tomographic Scans in TRISO Fuel Graphitic Matrix and Nuclear Grade Graphites

An image reconstruction method was developed to rectify shortcomings of earlier methods that became apparent as destructively sampled data became available. This reconstruction method was applied to the tomographic gamma scans of nuclear graphite and graphitic matrix samples from AGR-3/4. There is generally agreement between profiles measured via destructive sampling of these rings and profiles from tomographic reconstruction.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Bootstrap-determined p values in lattice QCD

We present a general method to determine the probability that stochastic Monte Carlo data, in particular those generated in a lattice QCD calculation, would have been obtained were that data drawn from the distribution predicted by a given theoretical hypothesis. Such a probability, or p -value, is often used as an important heuristic measure of the validity of that hypothesis. The proposed method offers the benefit that it remains usable in cases where the standard Hotelling T 2 methods based on the conventional χ 2 statistic do not apply, such as for uncorrelated fits. Specifically, we analyze q 2 , defined as the correlated χ 2 statistic obtained using an arbitrary covariance matrix estimator, and show how to use the bootstrap as a data-driven method to determine the expected distribution of q 2 for a given hypothesis with minimal assumptions. This distribution can then be used to determine the p -value for a fit to the data. We also describe a bootstrap approach for quantifying the impact upon this p -value of estimating population parameters from a single ensemble of N samples. The overall method is accurate up to a 1 / N bias which we do not attempt to quantify. Published by the American Physical Society 2025

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Application of end-of-shift respirable crystalline silica monitoring to construction

A pilot project was conducted to determine the effect of common construction dusts as interferences in a new portable end-of-shift (EoS), direct-on-filter (DoF) sampling and analysis method for respirable crystalline silica (RCS), in this case, quartz. Construction dusts were prepared from plaster, drywall, cement and brick by grinding, aerosolizing, and collecting respirable dust with high flow rate cyclones. Filters were loaded with different levels of commercial α-quartz powder Min-u-Sil 5, and different levels of interfering dusts, singly and in combination. Samples were analyzed by Fourier Transform Infrared Spectroscopy (FTIR). Good correlations were found between nominal quartz loading (0 µg, 25 µg, 50 µg, and 100 µg) adjusted for quartz in the interfering dust and FTIR absorbance alone and in the presence of all interfering dusts. The slopes of the correlations were similar whether the loading was quartz without interference, or with plaster, drywall, and cement dusts, regardless of quantity. The results show that (a) plaster and drywall dusts do not interfere substantially; (b) cement does not interfere, but a change in the intercept suggests an effect on the background absorbance of the filter; and (c) in addition to having a substantial quartz content, brick dust contains an additional material, probably a silicate mineral, which interferes with the quartz peak. Thus, the presence of cement leads to lower quartz values and brick leads to higher values, but overall, 83% of the quartz contents predicted from the calibration data agreed within 50% of the adjusted nominal loadings within the range 20–110 µg. This result is encouraging given the high levels of interfering dusts. Nine samples loaded with smaller amounts of all four dusts together gave results within 25% of the adjusted nominal loadings. A single mixture addition of the dusts to the filter gave tighter variance in results than sequential additions. Unexpectedly, the two Certified Reference Materials (CRMs) 1878a and 1878b, gave different results when used to calibrate XRD analysis of Min-u-Sil 5.

42 ENGINEERING↗
DOE OSTI2025

Comparison of measurement techniques and sorption of radium-226 in low and high salinity aqueous samples

Human activities have the potential to redistribute radium (Ra) in the marine environment in a manner that may necessitate monitoring or management of subsequent human or environmental exposures. There is therefore a need to identify accurate and accessible techniques for Ra measurement in high salinity samples and to describe the distribution of Ra in estuarine and marine environments, but most efforts in these areas have focused on low salinity matrices. In addition, rapid and reliable measurements are crucial for time-sensitive samples such as short-lived isotopes or emergency situations. The objective of this study is to describe the limits of detection, cost, and relative ease for measurement of Ra in both low and high salinity aqueous samples via three analytical methods: liquid scintillation counting (LSC), high purity germanium (HPGe) gamma spectrometry, and inductively coupled plasma mass spectrometry (ICP-MS). To contextualize these measurements for real-world scenarios, the partitioning of 226 Ra to substrates relevant to the marine environment was also characterized. Although HPGe detection with solid phase extraction had the lowest limit of detection for low salinity samples (0.27 Bq L −1 ), poor 226 Ra recovery for high salinity samples and high materials costs make this method prohibitive for many users. Limits of detection for high salinity samples were lower for LSC (1.28 Bq L −1 ) than for ICP-MS without dilution (11.4 Bq L −1 ), but significant and unexpected degradation of the high salinity LSC standards was observed after six months. Furthermore, our preferred measurement method for high salinity Ra samples is ICP-MS with sample dilution as necessary to reduce matrix effects.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2022

Testing the role of AGN on the star formation and metal enrichment of ‘twin galaxies’

Here, we explore the effect of AGN activity on the star formation history of galaxies by analysing the stellar population properties of 10 pairs of nearby twin galaxies – selected as being visually similar except for the presence of an AGN. The selection of such twin samples represents a method to study AGN feedback, as recently proposed by del Moral Castro et al. We use integral field unit (IFU) data from CALIFA, stacked within three fixed apertures. AGN galaxies in a twin pair suggest more evolved stellar populations than their non-AGN counterpart 90 per cent of the time, regardless of aperture size. A comparison with a large sample from SDSS confirms that most twins are representative of the general population, but in each twin the differences between twin members is significant. A set of targeted line strengths reveal the AGN member of a twin pair is older and more metal-rich than the non-AGN galaxy, suggesting AGN galaxies in our sample may either have an earlier formation time or follow a different star formation and chemical enrichment history. These results are discussed within two simple, contrasting hypotheses for the role played by AGN in galaxy evolution, which can be tested in the future at a greater detail with the use of larger data sets.

79 ASTRONOMY AND ASTROPHYSICS↗