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

CIRCLEZ : Reliable photometric redshifts for active galactic nuclei computed solely using photometry from Legacy Survey Imaging for DESI

Photometric redshifts for galaxies hosting an accreting supermassive black hole in their center, known as active galactic nuclei (AGNs), are notoriously challenging. At present, they are most optimally computed via spectral energy distribution (SED) fittings, assuming that deep photometry for many wavelengths is available. However, for AGNs detected from all-sky surveys, the photometry is limited and provided by a range of instruments and studies. This makes the task of homogenizing the data challenging, presenting a dramatic drawback for the millions of AGNs that wide surveys such as SRG/eROSITA are poised to detect. This work aims to compute reliable photometric redshifts for X-ray-detected AGNs using only one dataset that covers a large area: the tenth data release of the Imaging Legacy Survey (LS10) for DESI. LS10 provides deep grizW1-W4 forced photometry within various apertures over the footprint of the eROSITA-DE survey, which avoids issues related to the cross-calibration of surveys. We present the results from CIRCLEZ, a machine-learning algorithm based on a fully connected neural network. CIRCLEZ is built on a training sample of 14 000 X-ray-detected AGNs and utilizes multi-aperture photometry, mapping the light distribution of the sources. The accuracy (σNMAD) and the fraction of outliers (η) reached in a test sample of 2913 AGNs are equal to 0.067 and 11.6%, respectively. The results are comparable to (or even better than) what was previously obtained for the same field, but with much less effort in this instance. We further tested the stability of the results by computing the photometric redshifts for the sources detected in CSC2 and Chandra-COSMOS Legacy, reaching a comparable accuracy as in eFEDS when limiting the magnitude of the counterparts to the depth of LS10. The method can be applied to fainter samples of AGNs using deeper optical data from future surveys (for example, LSST, Euclid), granting LS10-like information on the light distribution beyond the morphological type. Along with this paper, we have released an updated version of the photometric redshifts (including errors and probability distribution functions) for eROSITA/eFEDS.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Sampling two-dimensional isometric tensor network states

Sampling a quantum system’s underlying probability distributions is an important computational task, e.g., for quantum advantage experiments and quantum Monte Carlo algorithms. Tensor networks are an invaluable tool for efficiently representing states of large quantum systems with limited entanglement. Algorithms for sampling one-dimensional (1D) tensor networks are well-established and utilized in several 1D tensor network methods. In this paper we introduce two novel sampling algorithms for two-dimensional (2D) isometric tensor network states (isoTNS) that generalize existing 1D tensor network sampling algorithms. Our first proposed algorithm performs independent sampling and yields a single configuration together with its associated probability. The second algorithm employs a greedy search strategy to identify high-probability configurations and their corresponding probabilities. Numerical results demonstrate the effectiveness of these algorithms across quantum states with varying entanglement and system size.

Dumitrescu, Eugene [ORNL] (ORCID:0000000158519567)↗
DOE OSTI2022

Combining Spatial and Temporal Resolution in Cryo-TEM of Device Materials: Preprint

Cryogenic transmission electron microscopy (cryo-TEM) has been instrumental to reducing the effects of electron-beam sample interactions and enabling higher spatial resolution studies of beam sensitive materials. However, the process of thinning and freezing samples for cryo-TEM characterization is often complex and frequently involves removing materials from their environment or stimuli of interest for extended times. This can allow for structural relaxation, diffusion, and other dynamic processes to occur prior to characterization; for energy materials like battery electrodes, this makes it difficult to precisely correlate the imaged structure with native structures that evolve in-situ during device processing, operation, or aging. Here, we present methods to prepare cryo-TEM samples from energy materials during these critical processes in the relevant device lifetimes. This offers insights into improving the temporal resolution of cryo-TEM, while still exploiting its ability to achieve high spatial resolutions in characterizing beam sensitive materials.

BASIC BIOLOGICAL SCIENCES,MATERIALS SCIENCE↗
DOE OSTI2025

LDRD Abbreviated report: High-Order General-Discrete-Ordinates Method Enabling Efficient Deterministic Transport in Hydrodynamic Simulations

Deterministic transport simulations for national-security and energy applications often operate in high-dimensional phase-space, where accuracy and cost both become major challenges. A common numerical artifact in such problems is the “ray-effect,” which appears as unphysical streaks. Beyond misinterpretation, these artifacts can contaminate tightly coupled physics, such as fluid dynamics, radiation-hydrodynamics, and laser-plasma interactions, eroding the predictive capability of entire multiphysics workflows. Our objective was to make high-dimension studies practical on modern hardware while mitigating the ray-effect without relying on prohibitively expensive sampling approaches such as Monte Carlo methods. We developed the Generic Discretization Library (GenDiL), a Graphics Processing Unit (GPU)-first framework that uses high-order Discontinuous Galerkin (DG) methods and matrix-free algorithms to reduce memory usage and improve computational efficiency, critical for phase-space simulations. GenDiL supports phase-space adaptivity in both mesh size and polynomial order (hp-adaptivity) to place resolution only where it is needed. A central capability is Local Dimensional Refinement (LDR), which couples lower-dimension continuum models to higher-dimension kinetic models through stable and conservative interfaces, so that high-fidelity physics is applied only in regions where it is essential. Building on the GenDiL framework, we developed the General SN (GSN) family of algorithms as a true generalization of the polar SN approach (discrete ordinates, often denoted SN). Rather than tying discrete ordinates to a specific polar change of coordinates, GSN formulates transport on an arbitrary change of coordinates chosen to reduce ray-effect. We studied two complementary variants: an analytic variant, where the coordinate map is prescribed in advance by a closed-form function; and a data-driven variant, where a quantity of interest, such as the net flux, guides the coordinate system. GenDiL provides the library infrastructure for efficient GPU execution, but the GSN concept is algorithmic and independent of any one library. Across representative high-dimension tests, including non-symmetric solutions, both variants delivered strong ray-effect mitigation at practical cost, moving four- to six-dimensional analysis toward repeatable, routine studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Computing Potential of the Mean Force Profiles for Ion Permeation Through Channelrhodopsin Chimera, C1C2

Umbrella sampling, coupled with a weighted histogram analysis method (US-WHAM), can be used to construct potentials of mean force (PMFs) for studying the complex ion permeation pathways of membrane transport proteins. Despite the widespread use of US-WHAM, obtaining a physically meaningful PMF can be challenging. Here, we provide a protocol to resolve that issue. Then, we apply that protocol to compute a meaningful PMF for sodium ion permeation through channelrhodopsin chimera, C1C2, for illustration.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Elastic Modulus Measurement at High Temperatures for Miniature Ceramic Samples Using Laser Micro-Machining and Thermal Mechanical Analyzer

In this paper, we demonstrate a method of measuring the flexural elastic modulus of ceramics at an intermediate (~millimeter) scale at high temperatures. We used a picosecond laser to precisely cut microbeams from the location of interest in a bulk ceramic. They had a cross-section of approximately 100 μm × 300 μm and a length of ~1 cm. They were then tested in a thermal mechanical analyzer at room temperature, 500 °C, 800 °C, and 1100 °C using the four-point flexural testing method. We compared the elastic moduli of high-purity Al2O3 and AlN measured by our method with the reported values in the literature and found that the difference was less than 5% for both materials. This paper provides a new and accurate method of characterizing the high-temperature elastic modulus of miniature samples extracted from representative/selected areas of bulk materials.

Chemistry↗
DOE OSTI2024

Inactivation of highly transmissible livestock and avian viruses including influenza A and Newcastle disease virus for molecular diagnostics

There is a critical need for an inactivation method that completely inactivates pathogens at the time of sample collection while maintaining the nucleic acid quality required for diagnostic PCR testing. This inactivation method is required to alleviate concerns about transmission potential, minimize shipping complications and cost, and enable testing in lower containment laboratories, thereby enhancing disease diagnostics through improved turn-around time. This study evaluated a panel of 10 surrogate viruses that represent highly pathogenic animal diseases. These results showed that a commercial PrimeStore® molecular transport media (PSMTM) completely inactivated all viruses tested by >99.99%, as determined by infectivity and serial passage assays. However, the detection of viral nucleic acid by qRT-PCR was comparable in PSMTM and control-treated conditions. These results were consistent when viruses were evaluated in the presence of biological material such as sera and cloacal swabs to mimic diagnostic sample conditions for non-avian and avian viruses, respectively. The results of this study may be utilized by diagnostic testing laboratories for highly pathogenic agents affecting animal and human populations. These results may be used to revise guidance for select agent diagnostic testing and the shipment of infectious substances.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Learning the temporal evolution of multivariate densities via normalizing flows

In this work, we propose a method to learn multivariate probability distributions using sample path data from stochastic differential equations. Specifically, we consider temporally evolving probability distributions (e.g., those produced by integrating local or nonlocal Fokker–Planck equations). Here, we analyze this evolution through machine learning assisted construction of a time-dependent mapping that takes a reference distribution (say, a Gaussian) to each and every instance of our evolving distribution. If the reference distribution is the initial condition of a Fokker–Planck equation, what we learn is the time-T map of the corresponding solution. Specifically, the learned map is a multivariate normalizing flow that deforms the support of the reference density to the support of each and every density snapshot in time. We demonstrate that this approach can approximate probability density function evolutions in time from observed sampled data for systems driven by both Brownian and Lévy noise. We present examples with two- and three-dimensional, uni- and multimodal distributions to validate the method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Prediction of tissue optical properties using the Monte Carlo modeling of photon transport in turbid media and integrating spheres

Monte Carlo methods are an established technique for simulating light transport in biological tissue. Integrating spheres make experimental measurements of the reflectance and transmittance of a sample straightforward and inexpensive. This work presents an extension to existing Monte Carlo photon transport methods to simulate integrating sphere experiments. Crosstalk between spheres in dual-sphere experiments is accounted for in the method. Analytical models, previous works on Monte Carlo photon transport, and experimental measurements of a synthetic tissue phantom validate this method. We present two approaches for using this method to back-calculate the optical properties of samples. Experimental and simulation uncertainties are propagated through both methods. Both back-calculation methods find the optical properties of a sample accurately and precisely. Our model is implemented in standard Python 3 and CUDA C++ [J. Nickolls, I. Buck, M. Garland, and K. Skadron, ACM Queue 6 , 40 ( 2008 ) ] and is publicly available in Code 1.

Cook, Patrick D. (ORCID:0000000279345428)↗
DOE OSTI2026

Score-based deterministic density sampling

We propose a deterministic sampling framework using Score-Based Transport Modeling for sampling an unnormalized target density π given only its score ∇ log π. Our method approximates the Wasserstein gradient flow on KL($f_t$∥π) by learning the time-varying score ∇ log $f_t$ on the fly using score matching. While having the same marginal distribution as Langevin dynamics, our method produces smooth deterministic trajectories, resulting in monotone noise-free convergence. We prove that our method dissipates relative entropy at the same rate as the exact gradient flow, provided sufficient training. Numerical experiments validate our theoretical findings: our method converges at the optimal rate, has smooth trajectories, and is often more sample efficient than its stochastic counterpart. Experiments on high-dimensional image data show that our method produces high-quality generations in as few as 15 steps and exhibits natural exploratory behavior. The memory and runtime scale linearly in the sample size.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

The first high-redshift cavity power measurements of cool-core galaxy clusters with the International LOFAR Telescope

Radio-mode feedback associated with the active galactic nuclei (AGNs) at the cores of galaxy clusters injects a large amount of energy into the intracluster medium (ICM), offsetting radiative losses through X-ray emission. This mechanism prevents the ICM from rapidly cooling down and fueling extreme starburst activity as it accretes onto the central galaxies, and it is therefore a key ingredient in the evolution of galaxy clusters. However, the influence and mode of feedback at high redshifts (z ∼ 1) remains largely unknown. Low-frequency sub-arcsecond-resolution radio observations taken with the International LOFAR Telescope have demonstrated their ability to assist X-ray observations with constraining the energy output from the AGNs (or “cavity power”) in galaxy clusters, thereby enabling research at higher redshifts than before. In this pilot project, we tested this hybrid method on a high-redshift (0.6 < z < 1.3) sample of 13 galaxy clusters for the first time with the aim of verifying the performance of this method at these redshifts and providing the first estimates of the cavity power associated with the central AGN for a sample of distant clusters. We were able to detect clear radio lobes in three out of 13 galaxy clusters at redshifts of 0.7 < z < 0.9, and we used these detections in combination with ICM pressures surrounding the radio lobes obtained from standard profiles to calculate the corresponding cavity powers of the AGNs. Combining our results with the literature, the current data appear to suggest that the average cavity power peaked at a redshift ofz ∼ 0.4 and slowly decreases toward higher redshifts. However, we require more and tighter constraints on the cavity volume and a better understanding of our observational systematics to confirm any deviation of the cavity power trend from a constant level.

Astronomy & Astrophysics↗
DOE OSTI2025

Resolving Mixtures of Soot Characterized by SP-AMS Spectra Using a Latent Dirichlet Allocation Model

Soot produced by detonation or combustion events exhibits different chemical properties depending on the fuel, device construction, and environmental conditions in which the event occurs. These properties can be useful for defining relevant signatures for probabilistically identifying the different types of events that occurred, based on the soot that is produced from these events. However, it is rare to observe samples of soot from a detonation or combustion that are not contaminated by outside particles. In this paper, we present a method for resolving mixtures of soot to determine the contributions of sources that may be present in samples of recovered soot. We use Latent Dirichlet Allocation to describe the generative process for a sample of recovered soot, and use Variational Bayesian Inference to learn about the parameters associated with the generative model. We demonstrate the utility of this method by considering real samples of mixtures of soot under various frameworks to show that the model is able to identify the different components present in a sample of soot as well as their mixing proportions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Towards resolving protein structures at the atomic scale using atom probe tomography

In the field of Structural Biology, Atom Probe Tomography (APT) is in the nascent stages of development wherein coarse-grained visuals of proteins have been captured. The characterization of organic samples or biomolecules through the technique is currently limited to the detection of a few dominant signatures. The problem of indecipherable characterization can inherently be traced back to multiple forms of technique-specific responses to organic samples and consequent triggers leading to organic-sample and sample-medium interactions. While it is possible for captured manifestations of protein reconstructions to seemingly appear intact from a basic visual purview, the parameter-protein associative responses throughout the structure as a direct consequence of the inherent workings of the technique (until harmonized with organic sample complexity and behavior) and field evaporation-based factors make non-aberrative atomic associations infeasible. The work focused on identifying and theorizing the (above stated and other) fundamental mechanisms that stronghold the study of intricate atomic to higher order associations in proteins through APT. Attempts at structure elucidation of the cryogenic sample under study, through indirect associations (and methods) based on other imaging techniques, further revealed the distinct and highly distortive nature at the atomic scale deterring structural tunability and thus characterization of APT based cryogenic samples under analysis. As a direct counter to the atomic scale characterization problem, by taking the experiment-specific uncertainties, and probable APT-centric organic sample-based variabilities into account, a basic result is extracted and presented. Through mass-spectrometric and computational analysis, specific individual amino acids (Sulfur-containing protein-bound amino acids) in proteins and aspects of protein structure (probable backbone fragments, partial sequence - partial backbone portions) have been identified and characterized. Under analysis considerations, a few of the simplest known and easily inferable segments that favor structural deteriorations in the reconstructions are stated. Additionally, to overcome technique-specific deterrents to the characterization of biomolecules in cryogenic sample medium, the development of a protein-labeling strategy tailored to APT is suggested.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

An Integrated Framework for Memory-Centric Analysis: From Trace Collection to Co-Design

The memory wall phenomenon—where advances in processor performance significantly outpace those in memory subsystems—poses a fundamental challenge for contemporary computing systems. In memory-bound applications, memory subsystem behavior dominates performance, yet existing analysis approaches present significant limitations: detailed microarchitectural simulators require days to weeks to simulate modest workloads; hardware performance counters provide only aggregate statistics that obscure temporal and spatial access patterns; and scaled simulation approaches face challenges in capturing certain behaviors that emerge at larger scales. These limitations reflect a processor-centric design philosophy increasingly misaligned with memory-bound workloads where detailed understanding of memory access patterns, cache hierarchy interactions, and contention is critical for effective optimization. This paper presents an integrated framework for memory-centric analysis that enables effective hardware-software co-design. We describe practical trace collection techniques, including hardware-assisted processor tracing with minimal overhead and portable software-based instrumentation with statistical sampling. We present multi-perspective analysis methods that examine memory behavior from temporal, sequential, spatial, and relational viewpoints, revealing distinct optimization opportunities invisible in aggregate metrics. We detail an architectural modeling framework that uses sampled traces with temporal interpolation and confidence-based filtering to evaluate cache and memory configurations. Evaluation on representative benchmarks demonstrates that this framework achieves practical accuracy (L2 cache errors of 2.64\%, confidence-filtered L3 errors of 9.92\%, bandwidth errors of 7.33\%) while providing substantial speedup (26.8×) over cycle-accurate simulation, enabling rapid design space exploration. We demonstrate how this integrated framework enables systematic identification of both hardware optimizations (memory controller tuning, bank partitioning, NUMA configuration) and software optimizations (data layout restructuring, prefetching strategies, memory-aware scheduling). Through this comprehensive treatment of the memory-centric analysis pipeline—from trace collection through architectural modeling to co-design application—we provide researchers and practitioners with practical techniques for addressing memory bottlenecks in contemporary computing systems.

Gajaria, Dhruv Mayur↗
DOE OSTI2020

STOCHASTIC OPTIMIZATION FOR LONG TERM CAPITAL STRUCTURES, SYSTEMS, AND COMPONENTS REFURBISHMENT AND REPLACEMENT

As commercial nuclear power plants (NPPs) pursue extended plant operations in the form of Second License Renewals (SLRs), opportunities exist for these plants to provide capital investments to ensure long-term, safe, and economic performance. Several utilities have already announced their intention to pursue extended operations for one or more of their NPPs via SLR2. The goal of this research is to develop a riskinformed approach to evaluate and prioritize plant capital investments made in preparation for, and during the period of, extended plant operations to support decisions in NPP operations. In order to prioritize project selection via a riskinformed approach we developed a single decision-making tool that integrates safety/reliability, cost, and stochastic optimization models to provide users with data analysis capabilities to more cost effectively manage plant assets. Both stochastic analysis methods—such as Monte Carlo-based sampling strategies—and multi-stage stochastic optimization strategies are employed to provide priority lists to decisionmakers in support of risk-informed decisions. We applied the proposed method to a trial application of projected replacement/refurbishment expenditures for plant capital assets (i.e., Structures, Systems, and Components [SSCs]). The objective is to optimize the SSC replacement/refurbishment schedule in terms of economic constraints, data uncertainties, and SSC reliability data, as well to generate a priority list for maximizing returns on investment.

42 ENGINEERING↗
DOE OSTI2021

Determination of Seven Organic Impurities in FD&C Yellow No. 6 by Ultra-High-Performance Liquid Chromatography

The U.S. Food and Drug Administration (FDA) batch certifies FD&C Yellow No. 6 (Y6) to ensure that the color additive meets requirements published in the Code of Federal Regulations (CFR), including specifications for seven organic manufacturing impurities consisting of two intermediates, a reaction by-product, an impurity originating from an intermediate, and three subsidiary colors. An ultra-high-performance liquid chromatography (UHPLC) method was developed and validated for determining seven organic impurities in Y6. For comparison, the currently used high-performance liquid chromatography (HPLC) method was also validated. The new UHPLC method uses a 1.7 µm particle size biphenyl column with aqueous ammonium formate and methanol as eluants. Analytes are identified by comparing their retention times and UV-visible spectra to those of reference standards. Calibration is performed in the presence of the dye matrix and analyte levels are determined from their peak areas. UHPLC and HPLC validation studies obtained linear calibration curves and excellent values for limits of detection, limits of quantitation, recovery, repeatability, and intermediate precision for all analytes. Survey analyses of 30 samples using the UHPLC and HPLC methods yielded results that were consistent within experimental error; however, the UHPLC method provided more accurate results for two analytes that coelute using the HPLC method. The UHPLC method satisfies the accuracy and precision requisites for routine certification of Y6. The UHPLC method is faster and generates less waste than the HPLC method while delivering better sensitivity, separation, and accuracy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗