Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “error analysis”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

The study of frequency-scan photothermal reflectance technique for thermal diffusivity measurement

A frequency scan photothermal reflectance technique to measure thermal diffusivity of bulk samples is studied in this manuscript. Similar to general photothermal reflectance methods, an intensity-modulated heating laser and a constant intensity probe laser are used to determine the surface temperature response under sinusoidal heating. The approach involves fixing the distance between the heating and probe laser spots, recording the phase lag of reflected probe laser intensity with respect to the heating laser frequency modulation, and extracting thermal diffusivity using the phase lag – (frequency) 1/2 relation. The experimental validation is performed on three samples (SiO 2 , CaF 2 and Ge), which have a wide range of thermal diffusivities. The measured thermal diffusivity values agree closely with literature values. Lastly, compared to the commonly used spatial scan method, the experimental setup and operation of the frequency scan method are simplified, and the uncertainty level is equal to or smaller than that of the spatial scan method.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2015

Comparative study of 0° X-cut and Y + 36°-cut lithium niobate high-voltage sensing

A comparison study between Y + 36° and 0° X-cut lithium niobate (LiNbO{sub 3}) was performed to evaluate the influence of crystal cut on the acoustic propagation to realize a piezoelectric high-voltage sensor. The acoustic time-of-flight for each crystal cut was measured when applying direct current (DC), alternating current (AC), and pulsed voltages. Results show that the voltage-induced shift in the acoustic wave propagation time scaled quadratically with voltage for DC and AC voltages applied to X-cut crystals. For the Y + 36° crystal, the voltage-induced shift scales linearly with DC voltages and quadratically with AC voltages. When applying 5 μs voltage pulses to both crystals, the voltage-induced shift scaled linearly with voltage. For the Y + 36° cut, the voltage-induced shift from applying DC voltages ranged from 10 to 54 ps and 35 to 778 ps for AC voltages at 640 V over the frequency range of 100 Hz–100 kHz. Using the same conditions as the Y + 36° cut, the 0° X-cut crystal sensed a shift of 10–273 ps for DC voltages and 189–813 ps for AC voltage application. For 5 μs voltage pulses, the 0° X-cut crystal sensed a voltage induced shift of 0.250–2 ns and the Y + 36°-cut crystal sensed a time shift of 0.115–1.6 ns. This suggests a frequency sensitive response to voltage where the influence of the crystal cut was not a significant contributor under DC, AC, or pulsed voltage conditions. The measured DC data were compared to a 1-D impedance matrix model where the predicted incremental length changed as a function of voltage. When the voltage source error was eliminated through physical modeling from the uncertainty budget, the combined uncertainty of the sensor (within a 95% confidence interval) decreased to 0.0033% using a Y + 36°-cut crystal and 0.0032% using an X-cut crystal for all the voltage conditions used in this experiment.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2016

Validation metrics for turbulent plasma transport

Developing accurate models of plasma dynamics is essential for confident predictive modeling of current and future fusion devices. In modern computer science and engineering, formal verification and validation processes are used to assess model accuracy and establish confidence in the predictive capabilities of a given model. This paper provides an overview of the key guiding principles and best practices for the development of validation metrics, illustrated using examples from investigations of turbulent transport in magnetically confined plasmas. Particular emphasis is given to the importance of uncertainty quantification and its inclusion within the metrics, and the need for utilizing synthetic diagnostics to enable quantitatively meaningful comparisons between simulation and experiment. As a starting point, the structure of commonly used global transport model metrics and their limitations is reviewed. An alternate approach is then presented, which focuses upon comparisons of predicted local fluxes, fluctuations, and equilibrium gradients against observation. The utility of metrics based upon these comparisons is demonstrated by applying them to gyrokinetic predictions of turbulent transport in a variety of discharges performed on the DIII-D tokamak [J. L. Luxon, Nucl. Fusion 42, 614 (2002)], as part of a multi-year transport model validation activity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Estimating the performance of lithium beam measurements of current density and electron density in an H-mode pedestal

The lithium beam is an effective diagnostic tool for investigation of stability and particle transport in the pedestal. It was used successfully to measure edge current density on DIII-D, achieving qualitative agreement with neoclassical models. Electron density profiles were also measured. Proposed upgrades will continue these measurements with higher reliability as well as explore edge current measurements using spectroscopy. The optics will be redesigned to optimize throughput and aperture broadening and to replace the photomultiplier tubes with avalanche photodiodes. Here, the new system will yield detailed measurements of the pedestal, complementing existing diagnostics for investigating pedestal stability, edge localized mode cycle, and particle transport through the pedestal.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Implementation of a Wollaston interferometry diagnostic on OMEGA EP

Here, a Wollaston interferometer is presented for use in measuring the electron density of plasma plumes created in experiments on the OMEGA EP laser system. The diagnostic is installed as an additional arm on the 4ω probe system, a suite of diagnostics that share a 10 ps pulse of 263 nm laser light captured by an imaging system at f /4. The interferometer utilizes a Wollaston prism to create two angularly separated beams from a single input probe beam, split at any angle between 0° and 90°. This configuration is implemented uniquely such that fringe spacing may be altered independently of field of view, magnification, and imaging resolution, from a range of 17 to 76 μm/fringe. The region of overlap between the two beams forms a total field of view of approximately 1.2 × 1.6 mm at the target chamber center with an imaging resolution of 5 μm. Using this configuration, here it is shown that plasma density may be accurately characterized over a range of 3 × 10 18 –1 × 10 20 cm -3 .

47 OTHER INSTRUMENTATION↗
DOE OSTI2015

A STATISTICAL RECONSTRUCTION OF THE PLANET POPULATION AROUND KEPLER SOLAR-TYPE STARS

Using the cumulative catalog of planets detected by the NASA Kepler mission, we reconstruct the intrinsic occurrence of Earth- to Neptune-size (1-4 R {sub ⊕}) planets and their distributions with radius and orbital period. We analyze 76,711 solar-type (0.8 < R {sub *}/R {sub ☉} < 1.2) stars with 430 planets on 20-200 day orbits, excluding close-in planets that may have been affected by the proximity to the host star. Our analysis considers errors in planet radii and includes an ''iterative simulation'' technique that does not bin the data. We find a radius distribution that peaks at 2-2.8 Earth radii, with lower numbers of smaller and larger planets. These planets are uniformly distributed with logarithmic period, and the mean number of such planets per star is 0.46 ± 0.03. The occurrence is ∼0.66 if planets interior to 20 days are included. We estimate the occurrence of Earth-size planets in the ''habitable zone'' (defined as 1-2 R {sub ⊕}, 0.99-1.7 AU for solar-twin stars) as 6.4{sub −1.1}{sup +3.4}%. Our results largely agree with those of Petigura et al., although we find a higher occurrence of 2.8-4 Earth-radii planets. The reasons for this excess are the inclusion of errors in planet radius, updated Huber et al. stellar parameters, and also the exclusion of planets that may have been affected by proximity to the host star.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Three-Phase Power Converter-Based Real-Time Synchronous Generator Emulation

This paper develops a synchronous generator emulator by using a three-phase voltage source converter for transmission level power system testing. Different interface algorithms are compared, and the voltage type ideal transformer model is selected considering accuracy and stability. At the same time, closed-loop voltage control with current feed-forward is proposed to decrease the emulation error. The emulation is then verified through two different ways. First, the output waveforms of the emulator in experiments are compared with the simulation under the same condition. Second, a transfer function perturbation-based error model is obtained and redefined as the relative error for the amplitude and phase between the emulated and the target system over the frequency range of interest. Finally, the major cause of the error is investigated through a quantitative analysis of the error with varying parameters.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Absolute timing calibration of the Resolve microcalorimeter spectrometer on the X-ray Imaging and Spectroscopy Mission

The Resolve microcalorimeter spectrometer on the X-ray Imaging and Spectroscopy Mission (XRISM) is designed to have a good timing capability with the mission-level requirement of 1 ms as the absolute time tagging accuracy to suffice the needs for observatory science. In the ground calibration campaign, the absolute and relative timing offsets were measured using pulsed X-rays from the modulated X-ray sources. These offsets were used to determine calibration parameters, the Resolve timing coefficients, which are used in the offline correction of event times. In the orbit, reevaluation of the timing coefficients was carried out using the Crab pulsar, which was observed in two periods, one in the performance verification phase and the other in the guest observation phase. We report the absolute timing calibration of XRISM/Resolve using the ground and in-orbit data. Although the requirement is likely satisfied with the ground parameters, the timing coefficients have been refined using the in-orbit calibration to improve the timing accuracy. In addition, for the first time, we present the interpretation of the absolute timing offset originating from the analog and digital processing of X-ray events unique to a microcalorimeter spectrometer and present a complete view of the timing error of the Resolve instrument.

Astronomy and AstroPhysics↗
DOE OSTI2014

SU-E-T-310: Targeting Safety Improvements Through Analysis of Near-Miss Error Detection Points in An Incident Learning Database

Purpose: Radiation treatment planning involves a complex workflow that can make safety improvement efforts challenging. This study utilizes an incident reporting system to identify detection points of near-miss errors, in order to guide our departmental safety improvement efforts. Previous studies have examined where errors arise, but not where they are detected or their patterns. Methods: 1377 incidents were analyzed from a departmental nearmiss error reporting system from 3/2012–10/2013. All incidents were prospectively reviewed weekly by a multi-disciplinary team, and assigned a near-miss severity score ranging from 0–4 reflecting potential harm (no harm to critical). A 98-step consensus workflow was used to determine origination and detection points of near-miss errors, categorized into 7 major steps (patient assessment/orders, simulation, contouring/treatment planning, pre-treatment plan checks, therapist/on-treatment review, post-treatment checks, and equipment issues). Categories were compared using ANOVA. Results: In the 7-step workflow, 23% of near-miss errors were detected within the same step in the workflow, while an additional 37% were detected by the next step in the workflow, and 23% were detected two steps downstream. Errors detected further from origination were more severe (p<.001; Figure 1). The most common source of near-miss errors was treatment planning/contouring, with 476 near misses (35%). Of those 476, only 72(15%) were found before leaving treatment planning, 213(45%) were found at physics plan checks, and 191(40%) were caught at the therapist pre-treatment chart review or on portal imaging. Errors that passed through physics plan checks and were detected by therapists were more severe than other errors originating in contouring/treatment planning (1.81 vs 1.33, p<0.001). Conclusion: Errors caught by radiation treatment therapists tend to be more severe than errors caught earlier in the workflow, highlighting the importance of safety checks in dosimetry and physics. We are utilizing our findings to improve manual and automated checklists for dosimetry and physics.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

An Optimization-based Atomistic-to-Continuum Coupling Method

In this paper, we present a new optimization-based method for atomistic-to-continuum (AtC) coupling. The main idea is to cast the latter as a constrained optimization problem with virtual Dirichlet controls on the interfaces between the atomistic and continuum subdomains. The optimization objective is to minimize the error between the atomistic and continuum solutions on the overlap between the two subdomains, while the atomistic and continuum force balance equations provide the constraints. Separation, rather then blending of the atomistic and continuum problems, and their subsequent use as constraints in the optimization problem distinguishes our approach from the existing AtC formulations. Finally, we present and analyze the method in the context of a one-dimensional chain of atoms modeled using a linearized two-body potential with next-nearest neighbor interactions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

fPINNs: Fractional Physics-Informed Neural Networks

Physics-informed neural networks (PINNs), introduced, are effective in solving integer-order partial differential equations (PDEs) based on scattered and noisy data. PINNs employ standard feedforward neural networks (NNs) with the PDEs explicitly encoded into the NN using automatic differentiation, while the sum of the mean-squared PDE residuals and the mean-squared error in initial-boundary conditions is minimized with respect to the NN parameters. Here we extend PINNs to fractional PINNs (fPINNs) to solve space-time fractional advection-diffusion equations (fractional ADEs), and we study systematically their convergence, hence explaining both fPINNs and PINNs for the first time. Specifically, we demonstrate their accuracy and effectiveness in solving multidimensional forward and inverse problems with forcing terms whose values are only known at randomly scattered spatio-temporal coordinates (black-box (BB) forcing terms). A novel element of the fPINNs is the hybrid approach that we introduce for constructing the residual in the loss function using both automatic differentiation for the integer-order operators and numerical discretization for the fractional operators. This approach bypasses the difficulties stemming from the fact that automatic differentiation is not applicable to fractional operators because the standard chain rule in integer calculus is not valid in fractional calculus. To discretize the fractional operators, we employ the Grünwald--Letnikov (GL) formula in one-dimensional fractional ADEs and the vector GL formula in conjunction with the directional fractional Laplacian in two- and three-dimensional fractional ADEs. We first consider the one-dimensional fractional Poisson equation and compare the convergence of the fPINNs against the finite difference method (FDM). We present the solution convergence using both the mean L 2 error as well as the standard deviation due to sensitivity to NN parameter initializations. Using different GL formulas we observe first-, second-, and third-order convergence rates for small size training sets but the error saturates for larger training sets. We explain these results by analyzing the four sources of numerical errors due to discretization, sampling, NN approximation, and optimization. The total error decays monotonically (below 10 –5 for a third-order GL formula) but it saturates beyond that point due to the optimization error. We also analyze the relative balance between discretization and sampling errors and observe that the sampling size and the number of discretization points (auxiliary points) should be comparable to achieve the highest accuracy. As we increase the depth of the NN up to certain value, the mean error decreases and the standard deviation increases, whereas the width has essentially no effect unless its value is either too small or too large. We next consider time-dependent fractional ADEs and compare white-box (WB) and BB forcing. We observe that for the WB forcing, our results are similar to the aforementioned cases; however, for the BB forcing fPINNs outperform FDM. Subsequently, we consider multidimensional time-, space-, and space-time-fractional ADEs using the directional fractional Laplacian and we observe relative errors of 10 –3 ~ 10 –4 . Lastly, we solve several inverse problems in one, two, and three dimensions to identify the fractional orders, diffusion coefficients, and transport velocities and obtain accurate results given proper initializations even in the presence of significant noise.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A Fast Solver for the Fractional Helmholtz Equation

The purpose of this paper is to study a Helmholtz problem with a spectral fractional Laplacian, instead of the standard Laplacian. Recently, it has been established that such a fractional Helmholtz problem better captures the underlying behavior in Geophysical Electromagnetics. In this work, we establish the well-posedness and regularity of this problem. We introduce a hybrid spectral-finite element approach to discretize it and show well-posedness of the discrete system. In addition, we derive a priori discretization error estimates. Finally, we introduce an efficient solver that scales aswell as the best possible solver for the classical integer-order Helmholtz equation. We conclude withseveral illustrative examples that confirm our theoretical findings.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

On the Training and Generalization of Deep Operator Networks

Here, we present a novel training method for deep operator networks (DeepONets), one of the most popular neural network models for operators. DeepONets are constructed by two subnetworks, namely the branch and trunk networks. Typically, the two subnetworks are trained simultaneously, which amounts to solving a complex optimization problem in a high dimensional space. In addition, the nonconvex and nonlinear nature makes training very challenging. To tackle such a challenge, we propose a two-step training method that trains the trunk network first and then sequentially trains the branch network. The core mechanism is motivated by the divide-and-conquer paradigm and is the decomposition of the entire complex training task into two subtasks with reduced complexity. Therein the Gram–Schmidt orthonormalization process is introduced which significantly improves stability and generalization ability. On the theoretical side, we establish a generalization error estimate in terms of the number of training data, the width of DeepONets, and the number of input and output sensors. Numerical examples are presented to demonstrate the effectiveness of the two-step training method, including Darcy flow in heterogeneous porous media.

deep operator networks↗
DOE OSTI2025

Comparison of CERES SYN1deg Radiative Fluxes with Those Derived from Observations at the ARM ENA Site

Profiles of radiative fluxes simulated from thermodynamic and cloud observations made at the Atmospheric Radiation Measurement (ARM) eastern North Atlantic (ENA) site for a 6-yr period are termed as ENARad. ENARad radiative fluxes are compared to those from the Clouds and the Earth’s Radiant Energy System (CERES) 1°-resolution synoptic product (SYN1deg)-simulated radiative flux profiles as well as the CERES instrument observed top-of-the-atmosphere (TOA) fluxes and ground site broadband radiometer measurements. Monthly average differences between ENARad and surface radiometer reported fluxes and differences between ENARad, SYN1deg, and observed fluxes at TOA were statistically insignificant. SYN1deg significantly overestimated surface downwelling shortwave flux by 12 ± 52 W m −2 and surface downwelling longwave flux by 5 ± 20 W m −2 on monthly time scales. Such overestimations were traced to a moister and warmer subcloud layer, a drier cloud layer, and a moister and colder above-cloud-free troposphere in the ancillary thermodynamic and cloud properties used by SYN1deg than observed. Similarly, low-cloud coverage, boundaries, and liquid water paths utilized by SYN1deg were also significantly higher than observed. Intramodel differences in the hourly values of shortwave fluxes exceeded 100 W m−2 at the TOA and the surface. These differences were also due to inaccuracies in the representation of low-cloud properties within the SYN1deg product relative to those determined by ENA ARM instrumentation and used as ENARad ancillary data. Results presented are relevant to investigations employing the CERES SYN1deg data product, studies that estimate radiative fluxes from surface-based or satellite-borne observations, and comparative analyses of radiative fluxes derived using different methodological approaches.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Visual inspection reliability for precision manufactured parts

Sandia National Laboratories conducted an experiment for the National Nuclear Security Administration to determine the reliability of visual inspection of precision manufactured parts used in nuclear weapons. In addition visual inspection has been extensively researched since the early 20th century; however, the reliability of visual inspection for nuclear weapons parts has not been addressed. In addition, the efficacy of using inspector confidence ratings to guide multiple inspections in an effort to improve overall performance accuracy is unknown. Further, the workload associated with inspection has not been documented, and newer measures of stress have not been applied.

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

Variationally Derived Discontinuity Factors for the Asymptotic Homogenized Diffusion Equation

Here in this work, we derive and test variational discontinuity factors (DFs) for the asymptotic homogenized diffusion equation. We begin with a functional for optimally estimating the reactor multiplication factor, then introduce asymptotic expressions for the forward and adjoint angular fluxes, and finally require that all first-order error terms vanish. In this way, the reactor multiplication factor can be calculated with second-order error. The analysis leads to (1) an alternate derivation of the asymptotic homogenized diffusion equation, (2) variational boundary conditions for large periodic systems, and (3) variational DFs to be applied between adjacent periodic regions (e.g., fuel assemblies). Numerical tests show that applying the variational DFs to the asymptotic homogenized diffusion equation yields the most accurate estimates of the reactor multiplication factor compared to other DFs for a wide range of problems. However, the resulting assembly powers are less accurate than those obtained using other DFs for many realistic problems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2018

Utilizing Adjoint-Based Error Estimates for Surrogate Models to Accurately Predict Probabilities of Events

In thist study, we develop a procedure to utilize error estimates for samples of a surrogate model to compute robust upper and lower bounds on estimates of probabilities of events. We show that these error estimates can also be used in an adaptive algorithm to simultaneously reduce the computational cost and increase the accuracy in estimating probabilities of events using computationally expensive high-fidelity models. Specifically, we introduce the notion of reliability of a sample of a surrogate model, and we prove that utilizing the surrogate model for the reliable samples and the high-fidelity model for the unreliable samples gives precisely the same estimate of the probability of the output event as would be obtained by evaluation of the original model for each sample. The adaptive algorithm uses the additional evaluations of the high-fidelity model for the unreliable samples to locally improve the surrogate model near the limit state, which significantly reduces the number of high-fidelity model evaluations as the limit state is resolved. Numerical results based on a recently developed adjoint-based approach for estimating the error in samples of a surrogate are provided to demonstrate (1) the robustness of the bounds on the probability of an event, and (2) that the adaptive enhancement algorithm provides a more accurate estimate of the probability of the QoI event than standard response surface approximation methods at a lower computational cost.

97 MATHEMATICS AND COMPUTING↗