SEARCH · Engineering Papers

Results for “Importance sampling”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Predictive Design of Hybrid Improper Ferroelectric Double Perovskite Oxides

The computational design of suitable multiferroic double perovskite oxides requires finding materials that exhibit sizable polarization, magnetization, and coupling between them. Oxides with the chemical formula of AA'BB'O 6 with building blocks of ABO 3 single perovskite oxides in centrosymmetric Pnma symmetry are strong candidates that have been reported to satisfy such criteria. The system lowers to noncentrosymmetric, polar P2 1 symmetry if A/A' layered and B/B' rocksalt cation orderings are imposed. A detailed compositional search over a variety of chemical spaces followed by evaluating their polarization may lead to the identification of more of these compounds with ferroelectric ordering. The standard density functional theory practices to estimate polarization within the Berry phase formalism require the systems to be perfectly insulating. The number of compounds that can be evaluated using this method is therefore limited. In this work, we introduce a predictive learning strategy based on importance sampling to build a series of machine learning models using results from first-principles simulations to predict polarization and the corresponding switching barrier. The geometry-driven features related to charge states and cationic radii play key roles in predicting the switching barrier with complementary contributions from the key structural mode-based order parameters. These modes become important to draw reasonable predictions of polarization components from machine learning models. In conclusion, our predictive models identify candidates with high polarizations and low switching barriers from a pool of double perovskite oxides, suitable for future investigation for their potential applications in spintronic devices.

36 MATERIALS SCIENCE↗

NASA NTRS2012

The Violent Early Solar System, as Told by Lunar Sample Geochronology

One of the legacies of the samples collected by the Apollo and Luna missions is the link forged between radiometric ages of rocks and relative ages according to stratigraphic relationships and impact crater size-frequency distributions. Our current understanding of the history of the inner solar system is based on the relative chronology of individual planets, tied to the absolute geochronology of the Moon via these important samples. Sample ages have enabled us to infer that impact-melt breccias from Apollo 14 and 15 record the formation of the Imbrium Basin, those from the highland massifs at Apollo 17 record the age of Serenitatis, those from the KREEP-poor Apollo 16 site record the age of Nectaris, and materials from Luna 24 record the age of Crisium. Ejecta from smaller and younger craters Copernicus and Tycho were sampled at Apollo 12 and 17, respectively, and local craters such as Cone at Apollo 14, and North Ray and South Ray at Apollo 16 were also sampled and ages determined for those events. Much of what we understand about the lunar impact flux is based on these ages. Samples from these nearside locations reveal a preponderance of impact-disturbed or recrystallized ages between 3.75 and 3.95 billion years. Argon and lead loss (and correlated disturbances in the Rb-Sr system) have been attributed to metamorphism of the lunar crust by an enormous number of impacts in a brief pulse of time, called the Lunar Cataclysm or Late Heavy Bombardment. Subsequent high-precision geochronometric analyses of Apollo samples and lunar highlands meteorites show a wider range of ages, but very few older than 4 Ga. The paucity of ancient impact melt rocks has been interpreted to mean that either that most impact basins formed at this time, or that ejecta from the large, near-side, young basins dominates the Apollo samples.

Cohen, Barbara A.↗

NASA NTRS2012

The Violent Early Solar System, as Told by Lunar Sample Geochronology

One of the legacies of the samples collected by the Apollo and Luna missions is the link forged between radiometric ages of rocks and relative ages according to stratigraphic relationships and impact crater size-frequency distributions. Our current understanding of the history of the inner solar system is based on the relative chronology of individual planets, tied to the absolute geochronology of the Moon via these important samples. Samples from these nearside locations reveal a preponderance of impact-disturbed or recrystallized ages between 3.75 and 3.95 billion years. Argon and lead loss (and correlated disturbances in the Rb-Sr system) have been attributed to metamorphism of the lunar crust by an enormous number of impacts in a brief pulse of time, called the Lunar Cataclysm or Late Heavy Bombardment. Subsequent high-precision geochronometric analyses of Apollo samples and lunar highlands meteorites show a wider range of ages, but very few older than 4 Ga. The paucity of ancient impact melt rocks has been interpreted to mean that either that most impact basins formed at this time, or that ejecta from the large, near-side, young basins dominates the Apollo samples. Selenochronology is getting more complicated: new results question meaning of sample ages, crater counts, crater production functions, and the solar system itself. Improved geological mapping of lunar geologic units and boundaries using multiple remote sensing datasets. High-resolution image-based crater counting of discrete geologic units and relating them to location. Improved understanding of the regolith thickness and its global variation (GRAIL). Tying the sampling of impact-melt rocks to the lunar impact flux. Using improved techniques (magnetic fields, diffusion studies, isotopic analysis) on existing samples. New sample return from benchmark craters, particularly SPA, which appears in 2013 Decadal Survey.

Cohen, Barbara↗

NASA NTRS2013

The Violent Early Solar System, as Told by Sample Geochronology

Cohen, Barbara↗

DOE OSTI2025

Validating sequential Monte Carlo for gravitational-wave inference

Nested sampling (NS) is the preferred stochastic sampling algorithm for gravitational-wave inference for compact binary coalescences. It can handle the complex nature of the gravitational-wave likelihood surface and provides an estimate of the Bayesian model evidence. However, there is another class of algorithms that meets the same requirements, but has not been used for gravitational-wave analyses: sequential Monte Carlo (SMC), an extension of importance sampling that maps samples from an initial density to a target density via a series of intermediate densities. In this work, we validate a type of SMC algorithm, called persistent sampling (PS), for gravitational-wave inference. We consider a range of different scenarios including binary black holes and binary neutron stars and real and simulated data and show that PS produces results that are consistent with NS whilst being, on average, 2 times more efficient and 2.74 times faster. This demonstrates that PS is a viable alternative to NS that should be considered for future gravitational-wave analyses.

black hole mergers↗

DOE OSTI2022

Bilayer Hubbard model: Analysis based on the fermionic sign problem

The bilayer Hubbard model describes the antiferromagnet to spin singlet transition and, potentially, aspects of the physics of unconventional superconductors. Despite these important applications, significant aspects of its phase diagram in the interplane hopping t ⊥ versus on-site interaction U parameter space, at half filling, are largely in disagreement. Here we provide an analysis making use of the average sign of weights over the course of the importance sampling in quantum Monte Carlo simulations to resolve several central open questions. Specifically, this metric of the weights clarifies the finite-sized metallic regimes at small U. Furthermore, at strong interactions, it points to the existence of a crossover from a correlated to uncorrelated band insulator not yet explored in a variety of existing, unbiased numerical methods. Our paper demonstrates the versatility of using properties of the weights in quantum Monte Carlo simulations to reveal important physical characteristics of the models under study.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2018

Improved Hybrid Modeling of Spent Fuel Storage Facilities

This work developed a new computational method for improving the ability to calculate the neutron flux in deep-penetration radiation shielding problems that contain areas with strong streaming. The “gold standard” method for radiation transport is Monte Carlo (MC) as it samples the physics exactly and requires few approximations. Historically, however, MC was not useful for shielding problems because of the computational challenge of following particles through dense shields. Instead, deterministic methods, which are superior in term of computational effort for these problems types but are not as accurate, were used. Hybrid methods, which use deterministic solutions to improve MC calculations through a process called variance reduction, can make it tractable from a computational time and resource use perspective to use MC for deep-penetration shielding. Perhaps the most widespread and accessible of these methods are the Consistent Adjoint Driven Importance Sampling (CADIS) and Forward-Weighted CADIS (FW-CADIS) methods. For problems containing strong anisotropies, such as power plants with pipes through walls, spent fuel cask arrays, active interrogation, and locations with small air gaps or plates embedded in water or concrete, hybrid methods are still insufficiently accurate. In this work, a new method for generating variance reduction parameters for strongly anisotropic, deep penetration radiation shielding studies was developed. This method generates an alternate form of the adjoint scalar flux quantity, Φ Ω , which is used by both CADIS and FW-CADIS to generate variance reduction parameters for local and global response functions, respectively. The new method, called CADIS-Ω, was implemented in the Denovo/ADVANTG software. Results indicate that the flux generated by CADIS-Ω incorporates localized angular anisotropies in the flux more effectively than standard methods. CADIS-Ω outperformed CADIS in several test problems. This initial work indicates that CADIS- may be highly useful for shielding problems with strong angular anisotropies. This is a benefit to the public by increasing accuracy for lower computational effort for many problems that have energy, security, and economic importance.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2025

Dark Energy Survey Year 6 Results: Redshift Calibration of the MagLim++ Lens Sample

In this work, we derive and calibrate the redshift distribution of the MagLim++ lens galaxy sample used in the Dark Energy Survey Year 6 (DES Y6) 3x2pt cosmology analysis. The 3x2pt analysis combines galaxy clustering from the lens galaxy sample and weak gravitational lensing. The redshift distributions are inferred using the SOMPZ method - a Self-Organizing Map framework that combines deep-field multi-band photometry, wide-field data, and a synthetic source injection (Balrog) catalog. Key improvements over the DES Year 3 (Y3) calibration include a noise-weighted SOM metric, an expanded Balrog catalogue, and an improved scheme for propagating systematic uncertainties, which allows us to generate O($10^8$) redshift realizations that collectively span the dominant sources of uncertainty. These realizations are then combined with independent clustering-redshift measurements via importance sampling. The resulting calibration achieves typical uncertainties on the mean redshift of 1-2%, corresponding to a 20-30% average reduction relative to DES Y3. We compress the $n(z)$ uncertainties into a small number of orthogonal modes for use in cosmological inference. Marginalizing over these modes leads to only a minor degradation in cosmological constraints. This analysis establishes the MagLim++ sample as a robust lens sample for precision cosmology with DES Y6 and provides a scalable framework for future surveys.

Giannini, G. [Chicago U., Astron. Astrophys. Ctr.;↗

DOE OSTI2026

Dark Energy Survey Year 6 results: Redshift calibration of the MagLim++ lens sample

In this work, we derive and calibrate the redshift distribution of the MagLim++ lens galaxy sample used in the Dark Energy Survey Year 6 (DES Y6) 3 x 2pt cosmology analysis. The 3 x 2pt analysis combines galaxy clustering from the lens galaxy sample and weak gravitational lensing. The redshift distributions are inferred using the SOMPZ method - a Self-Organizing Map framework that combines deep-field multi-band photometry, wide-field data, and a synthetic source injection ( B alrog) catalog. Key improvements over the DES Year 3 (Y3) calibration include a noise-weighted SOM metric, an expanded Balrog catalogue, and an improved scheme for propagating systematic uncertainties, which allows us to generate O(10 8 ) redshift realizations that collectively span the dominant sources of uncertainty. These realizations are then combined with independent clustering-redshift measurements via importance sampling. The resulting calibration achieves typical uncertainties on the mean redshift of 1-2%, corresponding to a 20-30% average reduction relative to DES Y3. We compress the n(z) uncertainties into a small number of orthogonal modes for use in cosmological inference. Marginalizing over these modes leads to only a minor degradation in cosmological constraints. Here, this analysis establishes the MagLim++ sample as a robust lens sample for precision cosmology with DES Y6 and provides a scalable framework for future surveys.

dark energy↗

NASA NTRS1993

Experimental analysis of computer system dependability

This paper reviews an area which has evolved over the past 15 years: experimental analysis of computer system dependability. Methodologies and advances are discussed for three basic approaches used in the area: simulated fault injection, physical fault injection, and measurement-based analysis. The three approaches are suited, respectively, to dependability evaluation in the three phases of a system's life: design phase, prototype phase, and operational phase. Before the discussion of these phases, several statistical techniques used in the area are introduced. For each phase, a classification of research methods or study topics is outlined, followed by discussion of these methods or topics as well as representative studies. The statistical techniques introduced include the estimation of parameters and confidence intervals, probability distribution characterization, and several multivariate analysis methods. Importance sampling, a statistical technique used to accelerate Monte Carlo simulation, is also introduced. The discussion of simulated fault injection covers electrical-level, logic-level, and function-level fault injection methods as well as representative simulation environments such as FOCUS and DEPEND. The discussion of physical fault injection covers hardware, software, and radiation fault injection methods as well as several software and hybrid tools including FIAT, FERARI, HYBRID, and FINE. The discussion of measurement-based analysis covers measurement and data processing techniques, basic error characterization, dependency analysis, Markov reward modeling, software-dependability, and fault diagnosis. The discussion involves several important issues studies in the area, including fault models, fast simulation techniques, workload/failure dependency, correlated failures, and software fault tolerance.

Iyer, Ravishankar, K.↗

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 OSTI2018

Spectroscopic and Theoretical Studies of Fluorescence Effects in 2-Methylamino-5-(2,4-dihydroxyphenyl)-1,3,4-thiadiazole Induced by Molecular Aggregation

The article presents the results of fluorescence analyses of 2-methylamino-5-(2,4-dihydroxyphenyl)-1,3,4-thiadiazole (MDFT) in an aqueous environment. MDFT dissolved in aqueous solutions with a pH value in the range from 1 to 4.5 yielded an interesting effect of two clearly separated fluorescence emissions. In turn, a single fluorescence was observed in MDFT dissolved in water solutions with a pH value from 4.5 to 12. As it was suggested in the previous investigations of other 1,3,4-thiadiazole compounds, these effects may be associated with conformational changes in the structure of the analysed molecule accompanied by aggregation effects. Crystallographic data showed that the effect of the two separated fluorescence emissions occurred in a conformation with the –OH group in the resorcyl ring bound on the side of the sulphur atom from the 1,3,4-thiadiazole ring. The hypothesis of aggregation as the mechanism involved in the change in the spectral properties at low pH is supported by the results of (Time-Dependent) Density Functional Theory calculations. The possibility of rapid analysis of conformational changes with the fluorescence spectroscopy technique may be rather important outcome obtained from the spectroscopic studies presented in this article. Additionally, the presented results seem to be highly important as they can be easily observed in solutions and biologically important samples.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2023

Removal of fission product tellurium and iodine from recently irradiated HEU using copper metal

Gamma detection of select fission products or activation products such as 153 Sm and 239 Np in a fission sample is difficult due to gamma interferences from radiotellurium and radioiodine. To aid in the detection of these difficult-to-detect nuclides, the radiotellurium and radioiodine were removed using Cu metal turnings. This method is highly effective at the removal of the interfering radionuclides removing > 98% of the radiotellurium and radioiodine. A mechanistic study was conducted. Determining the deposition mechanism, using a number of techniques including XANES and electrochemistry providing an understanding of the electrochemical path for the Te deposition. Initial work focused on the use of stable Te, further expanding into activated TeO 2 , and a fresh fission sample. Importantly, this method presents a rapid method for the direct detection of select fission products that could have a significant impact on nuclear related fields including nuclear forensics.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

DOE OSTI2019

Particle filter-based modulation domain infrared targets tracking

Faced with problems of low contrast, poor SNR, and relatively complicated tracking environment, stable infrared target tracking is worth researching for its many potential applications. In this paper, instead of traditional target tracking in the pixel domain, we propose a sampling importance resampling (SIR) particle filter method with indirect velocity measurements to track infrared targets in the modulation domain. The dominant amplitude modulation (AM) features used for tracking is extracted by decomposing the input image using an 18-channel Gabor filter bank followed by the application of the dominant component analysis approach. The dominant AM modulation features provide a significant partial texture characteristic of the target which can be separated from background with better discrimination. To take advantage of observed kinematics, we utilize the augmented state vector with indirect velocity information via combining the measurements of velocity in adjacent frames to the SIR particle filter framework, which weakens weights of particles with bad velocity estimates but still having association with the cluttered background or other moving objects. A dynamic template update strategy is also provided to prevent the tracker from appearance model drift. Experiments indicate that the proposed method is effective for raising the tracking accuracy compared with other tracking methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2018

Radiation dose rate distributions of spent fuel dry casks estimated with MAVRIC based on detailed geometry and continuous-energy models

This study presents a detailed comparison of the dose rate distributions of a dry TN–32 fuel cask with two geometry models and two cross-sectional datasets. The accuracies of radiation dose rate estimation and computational efficiencies of each geometry model with two cross-sectional data-sets are compared. The use of automated variance reduction techniques can significantly improve the computational efficiency of such a realistic, deep penetration problem that involves radiation transport from different volumetric sources, thereby eliciting only a small statistical error. Monaco with Automated Variance Reduction using Importance Calculations (MAVRIC) is a computational sequence within the SCALE 6.2 code package based on consistent adjoint driven importance sampling (CADIS), a type of automated variance reduction technique. Homogenous and full fuel assembly models are built herein, and two nuclear cross-section libraries (V7–200N47G and continuous energy) are applied in this work. Based on the detailed comparisons, we found that neutron dose rate estimation is more dependent on geometry modeling than on cross-section data. For neutron-induced gamma rays, the dose rate distribution depends on both the spatial self-shielding effect and the cross-section library. The primary gamma rays respectively contribute to the total dose rate by ~91% and ~ 99% on side and top surfaces, and the dose rate accuracy is more dependent on the cross-section library than on geometry modeling. In terms of the computation efficiency and efforts spent on geometry modeling, the homogenous fuel assembly model with the MG library can produce an acceptable dose rate distribution, but the detailed fuel assembly model with the continuous-energy library is required for more precise dose rate estimation.

61 RADIATION PROTECTION AND DOSIMETRY↗

DOE OSTI2019

Bayesian modeling of inconsistent plastic response due to material variability

The advent of fabrication techniques such as additive manufacturing has focused attention on the considerable variability of material response due to defects and other microstructural aspects. This variability motivates the development of an enhanced design methodology that incorporates inherent material variability to provide robust predictions of performance. In this work, we develop plasticity models capable of representing the distribution of mechanical responses observed in experiments using traditional plasticity models of the mean response and recently developed uncertainty quantification (UQ) techniques. To account for material response variability through variations in physical parameters, we adapt a recent Bayesian embedded modeling error calibration technique. We use Bayesian model selection to determine the most plausible of a variety of plasticity models and the optimal embedding of parameter variability. To expedite model selection, we develop an adaptive importance-sampling-based numerical integration scheme to compute the Bayesian model evidence. In conclusion, we demonstrate that the new framework provides predictive realizations that are superior to more traditional ones, and how these UQ techniques can be used in model selection and assessing the quality of calibrated physical parameters.

36 MATERIALS SCIENCE↗

DOE OSTI2020

On the quantification and efficient propagation of imprecise probabilities with copula dependence

This paper addresses the problem of quantification and propagation of uncertainties associated with dependence modeling when data for characterizing probability models are limited. Practically, the system inputs are often assumed to be mutually independent or correlated by a multivariate Gaussian distribution. However, this subjective assumption may introduce bias in the response estimate if the real dependence structure deviates from this assumption. In this work, we overcome this limitation by introducing a flexible copula dependence model to capture complex dependencies. Here, a hierarchical Bayesian multimodel approach is proposed to quantify uncertainty in dependence model-form and model parameters that result from small data sets. This approach begins by identifying, through Bayesian multimodel inference, a set of candidate marginal models and their corresponding model probabilities, and then estimating the uncertainty in the copula-based dependence structure, which is conditional on the marginals and their parameters. The overall uncertainties integrating marginals and copulas are probabilistically represented by an ensemble of multivariate candidate densities. A novel importance sampling reweighting approach is proposed to efficiently propagate the overall uncertainties through a computational model. Through an example studying the influence of constituent properties on the out-of-plane properties of transversely isotropic E-glass fiber composites, we show that the composite property with copula-based dependence model converges to the true estimate as data set size increases, while an independence or arbitrary Gaussian correlation assumption leads to a biased estimate.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Photon detector response function methodology using MCNP and shift hybrid radiation transport code for wide-area contamination assay applications

Here, radiation transport modeling using the Monte Carlo N-Particle (MCNP) radiation transport code and Monte Carlo code, Shift, were employed to model detector responses for a variety of wide-area photon contamination scenarios. In this study, 2" × 2" and 3" × 3" cylindrical NaI(Tl) scintillation detector configurations at source detector-distances of 0.5 cm, 1 cm, 2.54 cm, 10 cm, and 30 cm were modeled. Media of soil, concrete, and steel were evaluated for contamination depths ranging from surface to a depth of an infinite thickness in each medium for photon energies ranging from 20 keV to 3 MeV, which correspond to the energies that current detectors can discern. Monoenergetic photon surface contamination detector responses for each of the media, source–detector distances, and detectors were estimated using MCNP v6.2. Shift was harnessed for improved variance reduction of particle transport in highly attenuating media to obtain average cell fluxes in the two MCNP NaI(Tl) scintillation detector configurations. Average cell flux values in Shift were coupled with detector responses from MCNP to convert average cell flux in a void to energy distribution of pulses in the NaI(Tl) scintillation detector crystal of interest. An optimized detector response function methodology was developed by coupling the MCNP radiation transport method with the Consistent Adjoint Driven Importance Sampling (CADIS) hybrid radiation transport method built into Shift to significantly decrease the runtime of thousands of MCNP pulse height simulations. The methodology may be utilized to quickly and accurately facilitate the assessment of a broad range of wide-area environmental contamination assay and decommissioning cleanup applications.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗