Search NASA⌕ Search

SEARCH · Search NASA

Results for “Error Vector Magnitude”

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.

36 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Deep learning surrogate model for kinetic Landau-fluid closure with collision

In this work, the kinetic Landau-fluid (LF) closure with collision and periodic boundary condition is used in the development of the deep learning (DL) surrogate model. A classical neural network, namely, feedforward neural network or sometimes termed multilayer perceptron, is constructed and trained to learn the kinetic LF closure in the static limit and arbitrary mean free path in configuration space. The preliminary relation between best hyperparameters and critical parameters for data generation is found. Compared with the numerical approach (non-Fourier method) of the LF closure, the deep learning surrogate model shows an order of magnitude of improvement in terms of accuracy. Perhaps most importantly, the surrogate model closure has been integrated for the first time with fluid simulations. Our DL-enabled fluid simulations, for the first time, give the correct Landau damping rate for a wide range of wave vectors, while the Hammett–Perkins closure cannot produce the correct damping rate. We correctly connect the collisionless Hammett–Perkins closure and collisional Braginskii closure to reproduce the intrinsic nonlocal feature of the heat flux with DL techniques. We address the most concerning error accumulation problem and find that simulations with the deep learning surrogate model are as good as, if not better than, simulations with the analytic closure in terms of long-term numerical stability in the linear Landau damping test.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Single-qubit multi-party transmission using universal symmetric quantum cloning

This study considers the hypothetical quantum network case where Alice wishes to transmit one qubit of information (specifically a pure quantum state) to M parties, where M is some large number. The remote receivers locally perform single-qubit quantum state tomography on the transmitted qubits in order to compute the quantum state within some error rate (dependent on the tomography technique and the number of transmitted qubits). We show that with the use of an intermediate optimal symmetric universal quantum cloning machine (between Alice and the remote receivers) as a repeater-type node in a hypothetical quantum network, Alice can send significantly fewer qubits compared to direct transmission of the message qubits to each of the M remote receivers. This is possible due to two properties of quantum cloning. The first is that single qubit quantum clones retain the same Bloch angle as the initial quantum state. This means that if the mixed state of the quantum clone can be computed to high enough accuracy, the original pure quantum state can be inferred by extrapolating that vector to the surface of the Bloch sphere. The second property is that the state overlap of approximate quantum clones, with respect to the original pure quantum state, quickly converges (specifically for 1 → M , the limit of the fidelity as M goes to infinity is $\frac{2}{3}$). This means that Alice can prepare a constant number of qubits (which are then passed through the quantum cloning machine) in order to achieve a desired error rate if M is large enough. Combined, these two properties mean that for a large M , Alice can prepare many orders of magnitude fewer qubits in order to achieve the same single qubit transmission accuracy compared to the naive direct qubit transmission approach.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Resolution of a few problems in the application of quasilinear theory to calculating diffusion coefficients in heliophysics [Slides]

Three problems with the historical literature were discussed. Kennel and Engelmann (1966): ‘notational error’ that implies that diffusion coefficients depend on the distribution of wave power with respect to ψ. This dependence is unphysical and is not present in Lyons+ 1971. Kennel and Engelmann (1966) assume that only a single frequency can be associated with each wave vector: not necessarily true for EMIC, e.g. Lyons (1974a, 1974b) incorrectly used the Jacobian to relate the power spectral density in (ω,χ) space to that into (κ ⊥ , κ ∥ space. The correct use of the Jacobian: (1) Eliminates undesirable properties associated with N(ω) : lack of integrability, lack of superposition; and (2) Can produce orders of magnitude changes to the diffusion coefficients, but depends on the wave distribution: bigger effect for more oblique waves, and depends on the plasma parameter: bigger effect for smaller $ω_{pe}/ω_{ce}$.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Extension of Complex Refractive Index Measurements to the Near-Infrared for Liquids: Methodology and Uncertainty Analysis

Optical identification of liquid droplets, aerosols, or thin films is important for many applications. While reference spectra are sometimes available for such measurements, they are not always applicable to the observed spectrum or the given sample morphology. Reference spectra for many forms can be modeled, however, if the n/k vectors (real and imaginary refractive indices) are available. In previous work we have reported protocols to determine the n/k vectors for dozens of liquids, primarily in the mid-infrared (MIR) spectral range from 7500 to 400 cm –1 . In this work we extend the spectral range into the near-infrared (NIR) region, demonstrating a method to measure and merge the data sets to create composite n/k data ranging from 10 000 to 400 cm –1 (1.0 to 25 µm) with absorbance fidelity spanning over four orders of magnitude, and vastly improved signal-to-noise in the NIR. The precision of the composite data is evaluated for three different liquids, focusing primarily on the steps for converting the raw absorbance spectra to k values. The variability in both MIR and NIR data as well as in the final n/k vectors is also investigated for several liquids. For typical liquids, the overall variability (reported as 2σ) in the final n and k-vectors is determined to be ∼0.4% and 3%, respectively. Finally, the derived n/k data are used to calculate absorbance spectra for aerosol droplets, showing marginal variability due to the typical measurement errors in the final n/k vectors.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Accelerating the Discovery of New, Single Phase High Entropy Ceramics via Active Learning

High-entropy ceramics have garnered interest due to their remarkable hardness, compressive strength, thermal stability, and fracture toughness; yet the discovery of new high-entropy ceramics (out of a tremendous number of possible elemental permutations) still largely requires costly, inefficient, trial-and-error experimental and computational approaches. The entropy forming ability (EFA) factor was recently proposed as a computational descriptor that positively correlates with the likelihood that a 5-metal high-entropy carbide (HECs) will form the desired single phase, homogeneous solid solution; however, discovery of new compositions is computationally expensive. If you consider 8 candidate metals, the HEC EFA approach uses 49 optimizations for each of the 56 unique 5-metal carbides, requiring a total of 2744 costly density functional theory calculations. Here, we describe an orders-of-magnitude more efficient active learning (AL) approach for identifying novel HECs. To begin, we compared numerous methods for generating composition-based feature vectors (e.g., magpie and mat2vec), deployed an ensemble of machine learning (ML) models to generate an average and distribution of predictions, and then utilized the distribution as an uncertainty. Here we then deployed an AL approach to extract new training data points where the ensemble of ML models predicted a high EFA value or was uncertain of the prediction. Our approach has the combined benefit of decreasing the amount of training data required to reach acceptable prediction qualities and biases the predictions toward identifying HECs with the desired high EFA values, which are tentatively correlated with the formation of single phase HECs. Using this approach, we increased the number of 5-metal carbides screened from 56 to 15,504, revealing 4 compositions with record-high EFA values that were previously unreported in the literature. Our AL framework is also generalizable and could be modified to rationally predict optimized candidate materials/combinations with a wide range of desired properties (e.g., mechanical stability, thermal conductivity).

36 MATERIALS SCIENCE↗
DOE OSTI2017

Deflation as a method of variance reduction for estimating the trace of a matrix inverse

Many fields require computing the trace of the inverse of a large, sparse matrix. The typical method used for such computations is the Hutchinson method which is a Monte Carlo (MC) averaging over matrix quadratures. To improve its convergence, several variance reductions techniques have been proposed. In this paper, we study the effects of deflating the near null singular value space. We make two main contributions. First, we analyze the variance of the Hutchinson method as a function of the deflated singular values and vectors. Although this provides good intuition in general, by assuming additionally that the singular vectors are random unitary matrices, we arrive at concise formulas for the deflated variance that include only the variance and mean of the singular values. We make the remarkable observation that deflation may increase variance for Hermitian matrices but not for non-Hermitian ones. This is a rare, if not unique, property where non-Hermitian matrices outperform Hermitian ones. The theory can be used as a model for predicting the benefits of deflation. Second, we use deflation in the context of a large scale application of "disconnected diagrams" in Lattice QCD. On lattices, Hierarchical Probing (HP) has previously provided an order of magnitude of variance reduction over MC by removing "error" from neighboring nodes of increasing distance in the lattice. Although deflation used directly on MC yields a limited improvement of 30% in our problem, when combined with HP they reduce variance by a factor of over 150 compared to MC. For this, we pre-computated 1000 smallest singular values of an ill-conditioned matrix of size 25 million. Furthermore, using PRIMME and a domain-specific Algebraic Multigrid preconditioner, we perform one of the largest eigenvalue computations in Lattice QCD at a fraction of the cost of our trace computation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Fusion of Multiple Models for Improving Gross Primary Production Estimation With Eddy Covariance Data Based on Machine Learning

Abstract Terrestrial gross primary production (GPP) represents the magnitude of CO 2 uptake through vegetation photosynthesis, and is a key variable for carbon cycles between the biosphere and atmosphere. Light use efficiency (LUE) models have been widely used to estimate GPP for its physiological mechanisms and availability of data acquisition and implementation, yet each individual GPP model has exhibited large uncertainties due to input errors and model structure, and further studies of systematic validation, comparison, and fusion of those models with eddy covariance (EC) site data across diverse ecosystem types are still needed in order to further improve GPP estimation. We here compared and fused five GPP models (VPM, EC‐LUE, GOL‐PEM, CHJ, and C‐Fix) across eight ecosystems based on FLUXNET2015 data set using the ensemble methods of Bayesian Model Averaging (BMA), Support Vector Machine (SVM), and Random Forest (RF) separately. Our results showed that for individual models, EC‐LUE gave a better performance to capture interannual variability of GPP than other models, followed by VPM and GLO‐PEM, while CHJ and C‐Fix were more limited in their estimation performance. We found RF and SVM were superior to BMA on merging individual models at various plant functional types (PFTs) and at the scale of individual sites. On the basis of individual models, the fusion methods of BMA, SVM, and RF were examined by a five‐fold cross validation for each ecosystem type, and each method successfully improved the average accuracy of estimation by 8%, 18%, and 19%, respectively.

Environmental Sciences & Ecology↗
DOE OSTI2024

Distorted dislocation cores and asymmetric glide resistances in titanium

The determination of Critical Resolved Shear Stress (CRSS) in titanium for basal, prismatic, and pyramidal slip-planes without empirical constants is presented by combining Density Functional Theory (DFT) and anisotropic elasticity. A new mechanism leading to tension-compression (T-C) asymmetry in the CRSS levels has been revealed for the first time. The conditions for this asymmetry are established, involving a complex interplay between the dislocation core-structure and core-advance behavior. The three conditions for T-C asymmetry that need to be simultaneously satisfied can be summarized as: (1) a medium stacking fault width, d, (3 < d/b F < 10, where b F is the magnitude of the Burgers vector of the full dislocation), (2) an asymmetric core-structure of the extended dislocation (ξ 1 ≠ ξ 2 , where ξ 1 and ξ 2 are the core-widths of the first and second partials, respectively), and (3) intermittent motion of the partials (Δd/b F ≠ 0, where Δd is the magnitude of fluctuation in stacking fault width during intermittent motion). Pyramidal-slip in titanium satisfies all three conditions, resulting in significant T-C asymmetry. The CRSS theory considers a Wigner-Seitz (W-S) cell based domain area assigned to each lattice site for the calculations of core-energies accurately capturing the slip-plane lattice. The W-S based approach is essential due to the lower symmetry of the HCP crystal circumventing potential errors associated with the one-dimensional atomic-row approximation. Dislocation core structures are obtained for all the slip-systems in titanium showing significant distortion of the disregistry profile governing the core-advance behavior. The CRSS values predicted from the theory show agreement with the experimental CRSS levels reported in the literature.

Dislocation core↗
DOE OSTI2023

Multifidelity Neural Network Formulations for Prediction of Reactive Molecular Potential Energy Surfaces

Here, this paper focuses on the development of multifidelity modeling approaches using neural network surrogates, where training data arising from multiple model forms and resolutions are integrated to predict high-fidelity response quantities of interest at lower cost. We focus on the context of quantum chemistry and the integration of information from multiple levels of theory. Important foundations include the use of symmetry function-based atomic energy vector constructions as feature vectors for representing structures across families of molecules and single-fidelity neural network training capabilities that learn the relationships needed to map feature vectors to potential energy predictions. These foundations are embedded within several multifidelity topologies that decompose the high-fidelity mapping into model-based components, including sequential formulations that admit a general nonlinear mapping across fidelities and discrepancy-based formulations that presume an additive decomposition. Methodologies are first explored and demonstrated on a pair of simple analytical test problems and then deployed for potential energy prediction for C 5 H 5 using B2PLYP-D3/6-311++G(d,p) for high-fidelity simulation data and Hartree–Fock 6-31G for low-fidelity data. For the common case of limited access to high-fidelity data, our computational results demonstrate that multifidelity neural network potential energy surface constructions achieve roughly an order of magnitude improvement, either in terms of test error reduction for equivalent total simulation cost or reduction in total cost for equivalent error.

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

SU-E-T-291: Dosimetric Accuracy of Multitarget Single Isocenter Radiosurgery

Purpose: To evaluate the accuracy of single-isocenter multiple-target VMAT radiosurgery (SIMT-VMAT-SRS) by analysis of pre-treatment verification measurements. Methods: Our QA procedure used a phantom having a coronal plane for EDR2 film and a 0.125 cm3 ionization chamber. Film measurements were obtained for the largest and smallest targets for each plan. An ionization chamber measurement (ICM) was obtained for sufficiently large targets. Films were converted to dose using a patient-specific calibration curve and compared to treatment planning system calculations. Alignment error was estimated using image registration. The gamma index was calculated for 3%/3 and 3%/1 mm criteria. The median dose in the target region and, for plans having an ICM, the average dose in the central 5 mm was calculated. Results: The average equivalent target diameter of the 48 targets was 15 mm (3–43 mm). Twenty of the 24 plans had an ICM for the plan corresponding to the largest target (diameter 11–43 mm) with a mean ratio of chamber reading to expected dose (ED) and the mean ratio of film to ED (averaged over the central 5 mm) was 1.001 (0.025 SD) and 1.000 (0.029 SD), respectively. For all plans, the mean film to ED (from the median dose in the target region) was 0.997 (0.027 SD). The mean registration vector was (0.15,0.29) mm, with an average magnitude of 0.96 mm. Before (after) registration, the average fraction of pixels having gamma < 1 was 99.3% (99.6%) and 89.1% (97.6%) for 3%/3mm and 3%/1mm, respectively. Conclusion: Our results demonstrate dosimetric accuracy of SIMT-VMAT-SRS for targets as small as 3 mm. Film dosimetry provides accurate assessment of the absolute dose delivered to targets too small for an ionization chamber measurement; however, the relatively large registration vector indicates that image-guidance should replace laser-based setup for patient-specific evaluation of geometric accuracy.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-T-422: Fast Analytical Beamlet Optimization for Volumetric Intensity-Modulated Arc Therapy

Purpose: To implement a fast optimization algorithm on CPU/GPU heterogeneous computing platform and to obtain an optimal fluence for a given target dose distribution from the pre-calculated beamlets in an analytical approach. Methods: The 2D target dose distribution was modeled as an n-dimensional vector and estimated by a linear combination of independent basis vectors. The basis set was composed of the pre-calculated beamlet dose distributions at every 6 degrees of gantry angle and the cost function was set as the magnitude square of the vector difference between the target and the estimated dose distribution. The optimal weighting of the basis, which corresponds to the optimal fluence, was obtained analytically by the least square method. Those basis vectors with a positive weighting were selected for entering into the next level of optimization. Totally, 7 levels of optimization were implemented in the study.Ten head-and-neck and ten prostate carcinoma cases were selected for the study and mapped to a round water phantom with a diameter of 20cm. The Matlab computation was performed in a heterogeneous programming environment with Intel i7 CPU and NVIDIA Geforce 840M GPU. Results: In all selected cases, the estimated dose distribution was in a good agreement with the given target dose distribution and their correlation coefficients were found to be in the range of 0.9992 to 0.9997. Their root-mean-square error was monotonically decreasing and converging after 7 cycles of optimization. The computation took only about 10 seconds and the optimal fluence maps at each gantry angle throughout an arc were quickly obtained. Conclusion: An analytical approach is derived for finding the optimal fluence for a given target dose distribution and a fast optimization algorithm implemented on the CPU/GPU heterogeneous computing environment greatly reduces the optimization time.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

On approaches to analyze the sensitivity of simulated hydrologic fluxes to model parameters in the community land model

Here, effective sensitivity analysis approaches are needed to identify important parameters or factors and their uncertainties in complex Earth system models composed of multi-phase multi-component phenomena and multiple biogeophysical-biogeochemical processes. In this study, the impacts of 10 hydrologic parameters in the Community Land Model on simulations of runoff and latent heat flux are evaluated using data from a watershed. Different metrics, including residual statistics, the Nash-Sutcliffe coefficient, and log mean square error, are used as alternative measures of the deviations between the simulated and field observed values. Four sensitivity analysis (SA) approaches, including analysis of variance based on the generalized linear model, generalized cross validation based on the multivariate adaptive regression splines model, standardized regression coefficients based on a linear regression model, and analysis of variance based on support vector machine, are investigated. Results suggest that these approaches show consistent measurement of the impacts of major hydrologic parameters on response variables, but with differences in the relative contributions, particularly for the secondary parameters. The convergence behaviors of the SA with respect to the number of sampling points are also examined with different combinations of input parameter sets and output response variables and their alternative metrics. This study helps identify the optimal SA approach, provides guidance for the calibration of the Community Land Model parameters to improve the model simulations of land surface fluxes, and approximates the magnitudes to be adjusted in the parameter values during parametric model optimization.

58 GEOSCIENCES↗
DOE OSTI2022

J-PLUS: Support vector regression to measure stellar parameters

Stellar parameters are among the most important characteristics in studies of stars which, in traditional methods, are based on atmosphere models. However, time, cost, and brightness limits restrain the efficiency of spectral observations. The Javalambre Photometric Local Universe Survey (J-PLUS) is an observational campaign that aims to obtain photometry in 12 bands. Owing to its characteristics, J-PLUS data have become a valuable resource for studies of stars. Machine learning provides powerful tools for efficiently analyzing large data sets, such as the one from J-PLUS, and enables us to expand the research domain to stellar parameters. The main goal of this study is to construct a support vector regression (SVR) algorithm to estimate stellar parameters of the stars in the first data release of the J-PLUS observational campaign. The training data for the parameter's regressions are featured with 12-waveband photometry from J-PLUS and are crossidentified with spectrum-based catalogs. These catalogs are from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope, the Apache Point Observatory Galactic Evolution Experiment, and the Sloan Extension for Galactic Understanding and Exploration. We then label them with the stellar effective temperature, the surface gravity, and the metallicity. Ten percent of the sample is held out to apply a blind test. We develop a new method, a multi-model approach, in order to fully take into account, the uncertainties of both the magnitudes and the stellar parameters. The method utilizes more than 200 models to apply the uncertainty analysis. We present a catalog of 2 493 424 stars with the root mean square error of 160 K in the effective temperature regression, 0.35 in the surface gravity regression, and 0.25 in the metallicity regression. We also discuss the advantages of this multi-model approach and compare it to other machine-learning methods.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Exploiting Kronecker structure in exponential integrators: Fast approximation of the action of φ $-$functions of matrices via quadrature

Here, in this article, we propose an algorithm for approximating the action of φ $-$ functions of matrices against vectors, which is a key operation in exponential time integrators. In particular, we consider matrices with Kronecker sum structure, which arise from problems admitting a tensor product representation. The method is based on quadrature approximations of the integral form of the φ $-$ functions combined with a scaling and modified squaring method. Owing to the Kronecker sum representation, only actions of 1D matrix exponentials are needed at each quadrature node and assembly of the full matrix can be avoided. Additionally, we derive a priori bounds for the quadrature error, which show that, as expected by classical theory, the rate of convergence of our method is supergeometric. Guided by our analysis, we construct a fast and robust method for estimating the optimal scaling factor and number of quadrature nodes that minimizes the total cost for a prescribed error tolerance. We investigate the performance of our algorithm by solving several linear and semilinear time-dependent problems in 2D and 3D. The results show that our method is accurate and orders of magnitude faster than the current state-of-the-art.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

First-Principles Quantum Transport Modeling of Spin-Transfer and Spin-Orbit Torques in Magnetic Multilayers

A unified approach for computing (i) spin-transfer torque in magnetic trilayers like spin valve and magnetic tunnel junction, where injected charge current flows perpendicularly to interfaces, and (ii) spin-orbit torque in magnetic bilayers of the type ferromagnet/spin-orbit-coupled material, where injected charge current flows parallel to the interface, was reviewed. The experimentally explored and technologically relevant spin-orbit-coupled materials include 5d heavy metals, topological insulators, Weyl semimetals, and transition metal dichalcogenides. This approach requires to construct the torque operator for a given Hamiltonian of the device and the steady-state nonequilibrium density matrix, where the latter is expressed in terms of the nonequilibrium Green’s functions and split into three contributions. Tracing these contributions with the torque operator automatically yields field-like and damping-like components of spin-transfer torque or spin-orbit torque vector, which is particularly advantageous for spin-orbit torque where the direction of these components depends on the unknown-in-advance orientation of the current-driven nonequilibrium spin density in the presence of spin-orbit coupling. Illustrative examples are provided by computing spin-transfer torque in a one-dimensional toy model of a magnetic tunnel junction and realistic Co/Cu/Co spin valve, both of which are described by first-principles Hamiltonians obtained from noncollinear density functional theory calculations, as well as by computing spin-orbit torque in a ferromagnetic layer described by a tight-binding Hamiltonian which includes spin-orbit proximity effect within ferromagnetic monolayers assumed to be generated by the adjacent monolayer transition metal dichalcogenide. In addition, it is shown here how spin-orbit proximity effect, quantified by computing (via first-principles retarded Green’s function) spectral functions and spin textures on monolayers of realistic ferromagnetic material like Co in contact with heavy metal or monolayer transition metal dichalcogenide, can be tailored to enhance the magnitude of spin-orbit torque. Furthermore, errors made in the calculation of spin-transfer torque are quantified when using Hamiltonian from collinear density functional theory, with rigidly rotated magnetic moments to create noncollinear magnetization configurations, instead of proper (but computationally more expensive) self-consistent Hamiltonian obtained from noncollinear density functional theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Experimental and CFD Studies of Coolant Flow Mixing within Scaled Models of the Upper and Lower Plenums of NGNP Gas-Cooled Reactors

A 1/16th scaled VHTR experimental model was constructed and the preliminary test was performed in this study. To produce benchmark data for CFD validation in the future, the facility was first run at partial operation with five pipes being heated. PIV was performed to extract the vector velocity field for three adjacent naturally convective jets at statistically steady state. A small recirculation zone was found between the pipes, and the jets entered the merging zone at 3 cm from the pipe outlet but diverged as the flow approached the top of the test geometry. Turbulence analysis shows the turbulence intensity peaked at 41-45% as the jets mixed. A sensitivity analysis confirmed that 1000 frames were sufficient to measure statistically steady state. The results were then validated by extracting the flow rate from the PIV jet velocity profile, and comparing it with an analytic flow rate and ultrasonic flowmeter; all flow rates lie within the uncertainty of the other two methods for Tests 1 and 2. This test facility can be used for further analysis of naturally convective mixing, and eventually produce benchmark data for CFD validation for the VHTR during a PCC or DCC accident scenario. Next, a PTV study of 3000 images (1500 image pairs) were used to quantify the velocity field in the upper plenum. A sensitivity analysis confirmed that 1500 frames were sufficient to precisely estimate the flow. Subsequently, three (3, 9, and 15 cm) Y-lines from the pipe output were extracted to consider the output differences between 50 to 1500 frames. The average velocity field and standard deviation error that accrued in the three different tests were calculated to assess repeatability. The error was varied, from 1 to 14%, depending on Y-elevation. The error decreased as the flow moved farther from the output pipe. In addition, turbulent intensity was calculated and found to be high near the output. Reynolds stresses and turbulent intensity were used to validate the data by comparing it with benchmark data. The experimental data gave the same pattern as the benchmark data. A turbulent single buoyant jet study was performed for the case of LOFC in the upper plenum of scaled VHTR. Time-averaged profiles show that 3,000 frames of images were sufficient for the study up to second-order statistics. Self-similarity is an important feature of jets since the behavior of jets is independent of Reynolds number and a sole function of geometry. Self-similarity profiles were well observed in the axial velocity and velocity magnitude profile regardless of z/D where the radial velocity did not show any similarity pattern. The normal components of Reynolds stresses have self-similarity within the expected range. The study shows that large vortices were observed close to the dome wall, indicating that the geometry of the VHTR has a significant impact on its safety and performance. Near the dome surface, large vortices were shown to inhibit the flows, resulting in reduced axial jet velocity. The vortices that develop subsequently reduce the Reynolds stresses that develop and the impact on the integrity of the VHTR upper plenum surface. Multiple jets study, including two, three and five jets, were investigated.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2014

Harmonic analysis for the characterization and correction of geometric distortion in MRI

Purpose: Magnetic resonance imaging (MRI) is gaining widespread use in radiation therapy planning, patient setup verification, and real-time guidance of radiation delivery. Successful implementation of these technologies relies on the development of simple and efficient methods to characterize and monitor the geometric distortions arising due to system imperfections and gradient nonlinearities. To this end, the authors present the theory and validation of a novel harmonic approach to the quantification of system-related distortions in MRI. Methods: The theory of spatial encoding in MRI is applied to demonstrate that the 3D distortion vector field (DVF) is given by the solution of a second-order boundary value problem (BVP). This BVP is comprised of Laplace’s equation and a limited measurement of the distortion on the boundary of a specified region of interest (ROI). An analytical series expansion solving this BVP within a spherical ROI is obtained, and a statistical uncertainty analysis is performed to determine how random errors in the boundary measurements propagate to the ROI interior. This series expansion is then evaluated to obtain volumetric DVF mappings that are compared to reference data obtained on a 3 T full-body scanner. This validation is performed within two spheres of 20 cm diameter (one centered at the scanner origin and the other offset +3 cm along each of the transverse directions). Initially, a high-order mapping requiring measurements at 5810 boundary points is used. Then, after exploring the impact of the boundary sampling density and the effect of series truncation, a reduced-order mapping requiring measurements at 302 boundary points is evaluated. Results: The volumetric DVF mappings obtained from the harmonic analysis are in good agreement with the reference data. Following distortion correction using the high-order mapping, the authors estimate a reduction in the mean distortion magnitude from 0.86 to 0.42 mm and from 0.93 to 0.39 mm within the central and offset ROIs, respectively. In addition, the fraction of points with a distortion magnitude greater than 1 mm is reduced from 35.6% to 2.8% and from 40.4% to 1.5%, respectively. Similarly, following correction using the reduced-order mapping, the mean distortion magnitude reduces to 0.45–0.42 mm within the central and offset ROIs, and the fraction of points with a distortion magnitude greater than 1 mm is reduced to 2.8% and 1.5%, respectively. Conclusions: A novel harmonic approach to the characterization of system-related distortions in MRI is presented. This method permits a complete and accurate mapping of the DVF within a specified ROI using a limited measurement of the distortion on the ROI boundary. This technique eliminates the requirement to exhaustively sample the DVF at a dense 3D array of points, thereby permitting the design of simple, inexpensive phantoms that may incorporate additional modules for auxiliary quality assurance objectives.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2021

Nuclear Data and Cross Section Testing Using ENDF/B-VIII.0

With the release of the Evaluated Nuclear Data File (ENDF)/B-VIII.0 library, nuclear criticality safety practitioners and engineers have access to the latest cross section sets available for their analyses. However, these cross sections must be rigorously tested and validated to ensure that the nuclear data are responsive to the needs of the individuals responsible for developing, implementing, and maintaining computational tools for criticality safety applications. Thus, the ENDF/B-VIII.0 library is tested and validated with a large collection of experiments that were vetted by the International Criticality Safety Benchmark Evaluation Project and made available in the International Handbook of Evaluated Criticality Safety Benchmark Experiments. A selection of benchmark experiments for use within the criticality safety community were prepared and reviewed within the Verified, Archived Library of Inputs and Data (VALID), which is maintained by the Nuclear Energy and Fuel Cycle Division at Oak Ridge National Laboratory. The performance of the ENDF/B-VIII.0 library is assessed by using VALID models of benchmark experiments with the beta 12 version of SCALE 6.3 KENO V.a and KENO-VI Monte Carlo codes. The performance is compared with the results obtained from with the ENDF/B-VII.1 library. This report considers multigroup (MG) and continuous energy (CE) formats of the ENDF/B-VIII.0 and -VII.1 libraries. The benchmark experiments within VALID that validate the ENDF/B-VIII.0 library cover 15 broad system categories by using a range of fissile materials, uranium enrichments, plutonium isotopic vectors, and mixed uranium/plutonium systems. These forms are represented as metals, solutions, or various arrays of rods or plates that cover a variety of neutron energy spectra: thermal, fast, mixed, and intermediate. Over 600 cases were considered for use with the KENO V.a and KENO-VI codes with the ENDF/B-VIII.0 library. The results of the Monte Carlo comparison of ENDF/B-VIII.0 to ENDF/B-VII.1 with both KENO V.a and KENO-VI indicate that there is a less than 0.53% Δk difference between the bias of calculated k eff from the expected values. The CE ENDF/B-VIII.0 library results in smaller magnitude biases than the ENDF/B-VII.1 data for HEU-MET-FAST, HEU-SOL-THERM, IEU-MET-FAST, LEU-SOL-THERM, PU-SOL-THERM, and U233-MET-FAST systems, while the MG results yielded smaller magnitude biases for HEU-MET-FAST, HEU-SOL-THERM, IEU-MET-FAST, LEU-COMP-THERM, LEU-SOL THERM, MIX-COMP-FAST, and U233-MET-FAST systems. Most notable are the adjustments to the plutonium and 233 U cross section data, which has resulted in noticeably lower biases in the ENDF/B VIII.0 results for the mixed, plutonium, and 233 U systems. Results of the sensitivity data file comparison generated from TSUNAMI-3D for selected VALID cases for the ENDF/B-VIII.0 library indicate a very high level of agreement with correlation coefficients of the effect of nuclear data uncertainty on k eff (the c k integral parameter) all above 0.99. This indicates that cases with the ENDF/B-VIII.0 library would see very similar responses to any nuclear data errors or change as those with the ENDF/B-VII.1 library.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗