Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parameter 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.

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

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Trepanning Graphite Samples at a Shutdown Magnox Reactor: A US-UK Joint Exercise for Nuclear Verification

In April 2022, a team sponsored by the US DOE Office of Nuclear Verification (ONV) shipped trepanning equipment to the nuclear power station at Trawsfynydd, Wales, UK to collect samples from the graphite moderator of Reactor 1. The sampling campaign to collect 48 graphite samples was completed in five weeks beginning July 11. Half the samples will be analyzed in the US to determine the total energy production of the facility. The other half will be analyzed by Magnox Ltd. to support decisions for options for the ultimate decommissioning path for the facility. This paper provides a description of the equipment and the operations required to trepan, collect, and package the samples while controlling the small amounts of radioactive graphite dust entrained on and within the tool. The most important sample parameter needed is the precise location of the sample, to within ±6 mm in a fuel channel, while the tool is deployed to between 12 and 20 meters below the operations floor; our means of achieving this precision is given.

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

Comparison of NMC estimates with trans-stilbene, EJ-309, and He-3 detection systems

Neutron multiplicity counting (NMC) is a technique for the assay of fissile material. In this work, three detection systems are utilized for active interrogation assay of shells of the Rocky Flats shells (highly enriched uranium, 93% 235U) stacked from 13.25-54.92 kg assemblies. The singles ($R_1$) and doubles ($R_2$) rates are calculated with each system to estimate two sample parameters: $M_L$- the leakage multiplication and F - the sample fission rate. The estimated mass, m, is found by dividing F by a constant activity per unit mass. Since we are interrogating HEU, the α ratio of (α; n) to fission neutrons is taken to be zero. The system of equations to calculate these quantities was originally derived. The Neutron Multiplicity 3 He Array Detector (NOMAD) consists of 15 3 He tubes inside polyethylene and is the traditional, capture-based detection system for NMC. The polyethylene moderates incoming neutrons for thermal capture in individual tubes. The low gamma background, discrete capture signals, and high efficiency of the NOMAD are beneficial for NMC. However, the time to slow down neutrons to thermal energies leads to a system die away time on the scale of microseconds. Comparatively, the Rossi-alpha Measurements – Rapid Organic (n, γ) Discrimination Detector (RAMRODD) and the Organic Scintillator Array (OSCAR) are scatter-based detection systems. RAMRODD consists of 8, 5.08 by 5.08 cm EJ-309 liquid scintillators in 4 pairs, evenly-spaced with 90 degrees separation about the center of each assembly. OSCAR consists of a single array of 12, 5.08 by 5.08 cm trans-stilbene crystals aligned with the center of the assembly. Scatter-based systems detect fast neutrons without any moderation, leading to system die away times on the scale of tens of nanoseconds. This allows much shorter timing gate widths compared to thermal systems, thus increasing counting statistics of true correlated fission events. However, scatter-based systems are susceptible to neutron cross-talk when an incident neutron scatters off one detector and interacts in an adjacent detector, causing two seemingly correlated detection signals. The equations developed adjust for cross-talk to conduct NMC with scatter-based systems.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Online MCMC Thinning with Kernelized Stein Discrepancy

A fundamental challenge in Bayesian inference is efficient representation of a target distribution. Many nonparametric approaches do so by sampling a large number of points using variants of Markov chain Monte Carlo (MCMC). Here, we propose an MCMC variant that retains only those posterior samples which exceed a kernelized Stein discrepancy (KSD) threshold, which we call KSD thinning. We establish the convergence and complexity trade-offs for several settings of KSD thinning as a function of the KSD threshold parameter, sample size, and other problem parameters. We provide experimental comparisons against other online nonparametric Bayesian methods that generate low-complexity posterior representations. We observe superior consistency/complexity trade-offs across a range of settings including MCMC sampling on two Bayesian inference problems from the biological sciences, and 10 × inference speedup and storage reduction for Bayesian neural networks with no loss of accuracy and no increase in training time. Our code is available at https://github.com/colehawkins/KSD-Thinning.

Bayesian inference↗
DOE OSTI2020

A Possible Radiation-Induced Transition from Monazite-(Ce) to Xenotime-(Y)

This study examines two pegmatitic monazite samples (2a and 4b, these numbers are related to a previous study) to determine their crystal chemistry and effects of internal radiation damage using synchrotron high-resolution powder X-ray diffraction and electron-probe micro-analysis. Both the huttonite and cheralite substitutions are discussed. Rietveld structure refinement of sample 2a shows three different phases [2a = monazite-(Ce), 2b = monazite-(Ce), and 2c = xenotime-(Y)] with distinct structural parameters. The changes among the unit-cell parameters between the two monazite-(Ce) phases is more pronounced in the a followed by the b and c unit-cell parameters. Sample 4a is a single-phase monazite-(Sm) that contains 0.164 apfu Th. Phase 2c with space group I41/amd arises from redistribution of La, Ce, Pr, Nd, Sm, Gd, Dy, Si, and Y atoms from those in monazite (space group P21/n). A possible cause for the phase transition from monazite-(Ce) to xenotime-(Y) is α-radiation events over a long geological time. However, other chemical processes cannot be ruled out as a cause for the transition.

structural variations↗
NASA NTRS1970

FASTER 3: A generalized-geometry Monte Carlo computer program for the transport of neutrons and gamma rays. Volume 2: Users manual

A description of the FASTER-III program for Monte Carlo Carlo calculation of photon and neutron transport in complex geometries is presented. Major revisions include the capability of calculating minimum weight shield configurations for primary and secondary radiation and optimal importance sampling parameters. The program description includes a users manual describing the preparation of input data cards, the printout from a sample problem including the data card images, definitions of Fortran variables, the program logic, and the control cards required to run on the IBM 7094, IBM 360, UNIVAC 1108 and CDC 6600 computers.

Jordan, T. M.↗
DOE OSTI2022

Optimizing the accuracy of viscoelastic characterization with AFM force–distance experiments in the time and frequency domains

Atomic Force Microscopy (AFM) force-distance (FD) experiments have emerged as an attractive alternative to traditional micro-rheology measurement techniques owing to their versatility of use in materials of a wide range of mechanical properties. Here, we show that the range of time dependent behaviour which can reliably be resolved from the typical method of FD inversion (fitting constitutive FD relations to FD data) is inherently restricted by the experimental parameters: sampling frequency, experiment length, and strain rate. Specifically, we demonstrate that violating these restrictions can result in errors in the values of the parameters of the complex modulus. In the case of complex materials, such as cells, whose behaviour is not specifically understood a priori, the physical sensibility of these parameters cannot be assessed and may lead to falsely attributing a physical phenomenon to an artifact of the violation of these restrictions. We use arguments from information theory to understand the nature of these inconsistencies as well as devise limits on the range of mechanical parameters which can be reliably obtained from FD experiments. The results further demonstrate that the nature of these restrictions depends on the domain (time or frequency) used in the inversion process, with the time domain being far more restrictive than the frequency domain. Lastly, we demonstrate how to use these restrictions to better design FD experiments to target specific timescales of a material's behaviour through our analysis of a polydimethylsiloxane (PDMS) polymer sample.

information theory↗
DOE OSTI2022

Reactive burn model calibration using high-throughput initiation experiments at sub-millimeter length scales

We report a first-of-its-kind model calibration was performed using Sandia National Laboratories’ high-throughput initiation (HTI) experiment for two types of vapor-deposited explosive films consisting of hexanitrostilbene (HNS) or pentaerythritol tetranitrate (PETN). These films exhibit prompt initiation, and they reach steady detonation at sub-millimeter length scales. Following prior work on HNS, we test the hypothesis of approximating these explosive films as fine-grained homogeneous solids with simple Arrhenius kinetics burn models. The model calibration process is described herein using a single-step as well as a two-step Arrhenius rate law, and it consists of systematic parameter sampling leading to a reduction in the model degrees of freedom. Multiple local minima are observed; results are given for seven different optimized parameter sets. Each model set is further evaluated in a two-dimensional simulation of the critical failure thickness for a sustained detonation. Overall, the two-step Arrhenius kinetics model captures the observed behavior for HNS; however, neither model produces a good fit to the PETN data. We hypothesize that the HTI results for PETN correspond to a heterogeneous response, owing to the smaller reaction zone of PETN compared to HNS (i.e., it does not homogenize the fine-grained hot spots as well). Future work should consider using the ignition and growth model for PETN, as well as other reactive burn models such as xHVRB, AWSD, PiSURF, and CREST.

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

Molten Salt Sampling Techniques and Analytical Approaches

Recent global interest in pyroprocessing and molten salt reactors has brought salt sampling methods and techniques back to the forefront of nuclear safeguards concerns. Issues with uranium supplies have also encouraged various countries to pursue advanced nuclear fuel cycles. Tracking nuclear material in molten salt has proven to be a challenge and updating molten salt sampling will greatly help in this endeavor. Molten salt is problematic to sample due to salt stratification, lack of homogeneity, solids, and difficulty with hot cell adaptations. Various salt sampling techniques have been used since before the 1960s including surface, spoon/spatula, and bar solidification. Since then, new types of sampling techniques have been developed to improve sampling results. These include rod/dip, pipet, suction, filtered sampling along with devices such as the Valve Core Sampler and the Multi-Level Sampler. These different approaches are being analyzed and improved upon along with developing requirements for an improved salt sampling device. Work continues to develop salt samplers that are more robust, easier to segment, collect at a specific depth, can work with filters, and can collect fines. Sampling parameters are also being narrowed in terms of stirring, settling time, filtration, depth, etc. In the future, we hope to address deficiencies for process control and nuclear material accountancy control by determining the best way to collect samples that minimizes contaminants and is representative. A compilation of salt sampling approaches, analyses techniques, and an evaluation of findings will be presented.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Automated Construction of Artificial Lattice Structures with Designer Electronic States

Manipulating matter with a scanning tunneling microscope (STM) enables the creation of atomically defined artificial structures that host designer quantum states. However, the time-consuming nature of the manipulation process, coupled with the sensitivity of the STM tip, constrains the exploration of diverse configurations and limits the size of the designed features. In this study, we present a reinforcement learning (RL)-based framework for creating artificial structures by spatially manipulating carbon monoxide (CO) molecules on a copper substrate by using the STM tip. The automated workflow combines molecule detection and manipulation, employing deep-learning-based object detection to locate CO molecules and linear assignment algorithms to allocate these molecules to designated target sites. We initially perform molecule maneuvering based on randomized parameter sampling for sample bias, tunneling current set point, and manipulation speed. This data set is then structured into an action trajectory used to train an RL agent. The model is subsequently deployed on the STM for real-time fine-tuning of the manipulation parameters during structure construction. Our approach incorporates path-planning protocols coupled with active drift compensation to enable atomically precise fabrication of structures with significantly reduced human input while realizing larger-scale artificial lattices with the desired electronic properties. Furthermore, using our approach, we demonstrate the automated construction of an extended artificial graphene lattice and confirm the existence of a characteristic Dirac point in its electronic structure. Further challenges regarding the RL-based structural assembly scalability are discussed.

Algorithms↗
NASA NTRS1985

On percolation as a cosmological test

Difficulties in the use of percolation as a complementary statistic for the galaxy clustering pattern are pointed out by studying simple toy models and dynamical N-body models that represent the competing clustering scenarios. The percolation properties are found not to be very sensitive to the presence of pancakes and strings once they are clumpy, and hence they do not distinguish properly between models that are very different. In the case of very smooth pancakes, the ability to percolate depends on sampling parameters, such as the mean number density and the volume, in a way which is unknown a priori because it depends on the same properties that the test ought to measure. This problem could, in principle, be eased by using volume-limited samples of high mean number density (an order of magnitude denser than the CfA redshift survey volume limited at 4000 km/s) and by comparing to models of identical number density and volume. An alternative approach, based on the sampling effects themselves, may provide a qualitative test for pancakes in samples of lower densities.

Dekel, A.↗
NASA NTRS1982

Pilot opinions of sampling effects in lateral-directional control

Flight experiments with a microprocessor control system were conducted to determine the effects of variations in sampling parameters on several pilots' opinions of lateral-directional flying qualities. Princeton's variable-response research aircraft (VRA), which is equipped with a microprocessor based digital flight control system (Micro-DFCS), was the test vehicle. Two U.S. Navy pilots evaluated the effects of sampling rate, quantization, and pure time delay during tracking, approach, and landing. Aircraft carrier approach tasks were conducted using a Navy approach mirror. Acquisition and tracking of fixed objects on the ground provided additional information related to the Navy misson. The longitudinal controls were implemented with analog electronics, while the lateral-directional pilot inputs (stick and rudder) were fed to the Micro-DFCS, which commanded the ailerons and rudder. The conceptual relationship between the evaluation pilot's lateral-directional inputs, the flight computer, and the aircraft are illustrated.

Stengel, R. F.↗
DOE OSTI2021

Turn-key constrained parameter space exploration for particle accelerators using Bayesian active learning

Abstract Particle accelerators are invaluable discovery engines in the chemical, biological and physical sciences. Characterization of the accelerated beam response to accelerator input parameters is often the first step when conducting accelerator-based experiments. Currently used techniques for characterization, such as grid-like parameter sampling scans, become impractical when extended to higher dimensional input spaces, when complicated measurement constraints are present, or prior information known about the beam response is scarce. Here in this work, we describe an adaptation of the popular Bayesian optimization algorithm, which enables a turn-key exploration of input parameter spaces. Our algorithm replaces the need for parameter scans while minimizing prior information needed about the measurement’s behavior and associated measurement constraints. We experimentally demonstrate that our algorithm autonomously conducts an adaptive, multi-parameter exploration of input parameter space, potentially orders of magnitude faster than conventional grid-like parameter scans, while making highly constrained, single-shot beam phase-space measurements and accounts for costs associated with changing input parameters. In addition to applications in accelerator-based scientific experiments, this algorithm addresses challenges shared by many scientific disciplines, and is thus applicable to autonomously conducting experiments over a broad range of research topics.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Relation between sampling, sensitivity and precision in strain mapping using the Geometric Phase Analysis method in Scanning Transmission Electron Microscopy

The sensitivity and the precision of the Geometric Phase Analysis (GPA) method for strain characterization is a topic widely discussed in the literature and is usually difficult to quantify. Indeed, the GPA precision is intricately linked to the resolution of the strain maps defined when masking the periodic reflections in Fourier space. In this study an additional parameter, sampling, is proposed to be analyzed regarding the precision of GPA by developing the concept of a phase noise in the GPA equations. Both experimentally and theoretically, the following article demonstrates how the precision, and the sensitivity of the GPA method is improved when using a larger pixel spacing to record an electron micrograph in Scanning Transmission Electron Microscopy (STEM). In conclusion, the counterintuitive concept of increasing the field of view to improve the GPA precision results in an extension of the application of strain characterization methods in STEM towards low deformation levels.

47 OTHER INSTRUMENTATION↗
NASA NTRS1993

Electrotastic wave excitation in planetary magnetospheres: Application to Neptune

Voyager 2 observations of electrostatic electron and ion harmonic waves in Neptune's magnetosphere are addressed. A model of electron Bernstein modes generated by a loss cone distribution of superthermal electrons is scaled to Neptune parameters and a comparison of theory with the observed electron flux shows good agreement. A model of proton Bernstein modes generated by a ring distribution of Tritonogenic nitrogen ions is also investigated and satisfactory agreement with the data are obtained compatible with known properties of the magnetosphere. The success of the model in accounting for electrostatic emission observed by Voyager over a wide range of sampled parameters recommends its general applicability to planetary magnetospheres.

Barbosa, D. D.↗
NASA NTRS1976

Vector wind and vector wind shear models 0 to 27 km altitude for Cape Kennedy, Florida, and Vandenberg AFB, California

The techniques are presented to derive several statistical wind models. The techniques are from the properties of the multivariate normal probability function. Assuming that the winds can be considered as bivariate normally distributed, then (1) the wind components and conditional wind components are univariate normally distributed, (2) the wind speed is Rayleigh distributed, (3) the conditional distribution of wind speed given a wind direction is Rayleigh distributed, and (4) the frequency of wind direction can be derived. All of these distributions are derived from the 5-sample parameter of wind for the bivariate normal distribution. By further assuming that the winds at two altitudes are quadravariate normally distributed, then the vector wind shear is bivariate normally distributed and the modulus of the vector wind shear is Rayleigh distributed. The conditional probability of wind component shears given a wind component is normally distributed. Examples of these and other properties of the multivariate normal probability distribution function as applied to Cape Kennedy, Florida, and Vandenberg AFB, California, wind data samples are given. A technique to develop a synthetic vector wind profile model of interest to aerospace vehicle applications is presented.

Smith, O. E.↗
DOE OSTI2024

Dark energy survey year 3 results: likelihood-free, simulation-based w CDM inference with neural compression of weak-lensing map statistics

We present simulation-based cosmological wcold dark matter (wCDM) inference using dark energy survey year 3 weak-lensing maps, via neural data compression of weak-lensing map summary statistics: power spectra, peak counts, and direct map-level compression/inference with convolutional neural networks (CNN). Using simulation-based inference, also known as likelihood-free or implicit inference, we use forward-modelled mock data to estimate posterior probability distributions of unknown parameters. This approach allows all statistical assumptions and uncertainties to be propagated through the forward-modelled mock data; these include sky masks, non-Gaussian shape noise, shape measurement bias, source galaxy clustering, photometric redshift uncertainty, intrinsic galaxy alignments, non-Gaussian density fields, neutrinos, and non-linear summary statistics. We include a series of tests to validate our inference results. This paper also describes the Gower Street simulation suite: 791 full-sky pkdgrav3 dark matter simulations, with cosmological model parameters sampled with a mixed active-learning strategy, from which we construct over 3000 mock dark energy survey lensing data sets. For wCDM inference, for which we allow –1 < w < –$\frac{1}{3}$⁠, our most constraining result uses power spectra combined with map-level (CNN) inference. Using gravitational lensing data only, this map-level combination gives Ω m = 0.283$^{+0.020}_{–0.027}$⁠, S 8 = 0.804$^{+0.025}_{–0.017⁠}$, and w < –0.80 (with a 68 per cent credible interval); compared to the power spectrum inference, this is more than a factor of two improvement in dark energy parameter (Ω⁠ DE , w⁠) precision.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Residual-based error correction for neural operator accelerated infinite-dimensional Bayesian inverse problems

We explore using neural operators, or neural network representations of nonlinear maps between function spaces, to accelerate infinite-dimensional Bayesian inverse problems (BIPs) with models governed by nonlinear parametric partial differential equations (PDEs). Neural operators have gained significant attention in recent years for their ability to approximate the parameter-to-solution maps defined by PDEs using as training data solutions of PDEs at a limited number of parameter samples. The computational cost of BIPs can be drastically reduced if the large number of PDE solves required for posterior characterization are replaced with evaluations of trained neural operators. However, reducing error in the resulting BIP solutions via reducing the approximation error of the neural operators in training can be challenging and unreliable. We provide an a priori error bound result that implies certain BIPs can be ill-conditioned to the approximation error of neural operators, thus leading to inaccessible accuracy requirements in training. To reliably deploy neural operators in BIPs, we consider a strategy for enhancing the performance of neural operators: correcting the prediction of a trained neural operator by solving a linear variational problem based on the PDE residual. We show that a trained neural operator with error correction can achieve a quadratic reduction of its approximation error, all while retaining substantial computational speedups of posterior sampling when models are governed by highly nonlinear PDEs. The strategy is applied to two numerical examples of BIPs based on a nonlinear reaction–diffusion problem and deformation of hyperelastic materials. We demonstrate that posterior representations of the two BIPs produced using trained neural operators are greatly and consistently enhanced by error correction.

97 MATHEMATICS AND COMPUTING↗