Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical errors”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,423 records · Page 79
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Shadow Distillation: Quantum Error Mitigation with Classical Shadows for Near-Term Quantum Processors

Mitigating errors in quantum information processing devices is especially important in the absence of fault tolerance. An effective method in suppressing state-preparation errors is using multiple copies to distill the ideal component from a noisy quantum state. Here, we use classical shadows and randomized measurements to circumvent the need for coherent access to multiple copies at an exponential cost. We study the scaling of resources using numerical simulations and find that the overhead is still favorable compared to full state tomography. We optimize measurement resources under realistic experimental constraints and apply our method to an experiment preparing a Greenberger-Horne-Zeilinger state with trapped ions. In addition to improving stabilizer measurements, the analysis of the improved results reveals the nature of errors affecting the experiment. Hence, our results provide a directly applicable method for mitigating errors in near-term quantum computers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1981

Influence of boundary approximations and conditions on finite difference solutions

Numerical representations of boundary approximations and conditions for three problems are investigated to determine the resulting global accuracy of the steady state solution. Numerical accuracy with various boundary approximations is determined for quasi-one-dimensional inviscid flow in a duct with the interior grid points evaluated using the MacCormack scheme. When an extrapolation approximation with first order local truncation error is used, the global second order accuracy of the difference scheme can be destroyed. For one dimensional flow in a porous medium, an implicit midpoint difference scheme which is consistent with the boundary conditions is developed without the need of boundary approximations. A dissipative model problem is solved with the boundary conditions discretized with first and second order accuracy. The overall second order accuracy of the difference scheme is destroyed if first order numerical representation of one of the boundary conditions is used. With a boundary approximation, the second order global accuracy of the model problem is retained if either second order extrapolation or first order representation of the governing equation is used.

Blottner, F. G.↗
NASA NTRS2014

Neural Network Machine Learning and Dimension Reduction for Data Visualization

Neural network machine learning in computer science is a continuously developing field of study. Although neural network models have been developed which can accurately predict a numeric value or nominal classification, a general purpose method for constructing neural network architecture has yet to be developed. Computer scientists are often forced to rely on a trial-and-error process of developing and improving accurate neural network models. In many cases, models are constructed from a large number of input parameters. Understanding which input parameters have the greatest impact on the prediction of the model is often difficult to surmise, especially when the number of input variables is very high. This challenge is often labeled the "curse of dimensionality" in scientific fields. However, techniques exist for reducing the dimensionality of problems to just two dimensions. Once a problem's dimensions have been mapped to two dimensions, it can be easily plotted and understood by humans. The ability to visualize a multi-dimensional dataset can provide a means of identifying which input variables have the highest effect on determining a nominal or numeric output. Identifying these variables can provide a better means of training neural network models; models can be more easily and quickly trained using only input variables which appear to affect the outcome variable. The purpose of this project is to explore varying means of training neural networks and to utilize dimensional reduction for visualizing and understanding complex datasets.

Liles, Charles A.↗
DOE OSTI2023

Deep nonparametric estimation of intrinsic data structures by chart autoencoders: Generalization error and robustness

Autoencoders have demonstrated remarkable success in learning low-dimensional latent features of high-dimensional data across various applications. Assuming that data are sampled near a low-dimensional manifold, we employ chart autoencoders, which encode data into low-dimensional latent features on a collection of charts, preserving the topology and geometry of the data manifold. Our paper establishes statistical guarantees on the generalization error of chart autoencoders, and we demonstrate their denoising capabilities by considering n noisy training samples, along with their noise-free counterparts, on a d-dimensional manifold. By training autoencoders, we show that chart autoencoders can effectively denoise the input data with normal noise. We prove that, under proper network architectures, chart autoencoders achieve a squared generalization error in the order of n–$\frac{2}{d+2}$log 4 n, which depends on the intrinsic dimension of the manifold and only weakly depends on the ambient dimension and noise level. We further extend our theory on data with noise containing both normal and tangential components, where chart autoencoders still exhibit a denoising effect for the normal component. As a special case, our theory also applies to classical autoencoders, as long as the data manifold has a global parametrization. Furthermore, our results provide a solid theoretical foundation for the effectiveness of autoencoders, which is further validated through several numerical experiments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Thermomechanics coupling to Monte Carlo particle transport on unstructured mesh geometries using Cardinal

Geometry deformation due to thermal expansion influences neutron transport in many systems. Studying this phenomenon involves coupling models for neutronics, thermal hydraulics, and solid mechanics. To enable high fidelity modeling of these coupled physics, new capabilities were introduced in Cardinal, coupling OpenMC Monte Carlo particle transport models with MOOSE thermomechanical physics on unstructured moving-mesh geometries. In this work, we present a fully open-source capability leveraging on-the-fly mesh skinning to automatically regenerate OpenMC geometry, which allows multiphysics feedback from temperature, density, and geometry changes. The new capability is verified using an analytic benchmark slab problem, which couples S 2 neutron transport with thermal conduction, convective boundary conditions, Doppler-broadened cross sections, and nonlinear thermal expansion effects along the heated slab. Cardinal reproduces the analytic solutions for the neutron flux, heating, k eff , and temperature with demonstrated convergence in various error terms including mesh resolution and cross section temperature library spacing. For the nominal benchmark conditions and with a fine mesh, maximum relative errors for neutron flux, temperature, and heating are lower than 1%, while errors in integral quantities such as k eff and slab length are within 1 pcm and 48 µm, respectively. This work (i) presents a new numerical approach to thermomechanics coupling with OpenMC models, (ii) is the first (to our knowledge) to utilize a mechanical partial differential equation (PDE) solution to solve the (Griesheimer and Kooreman, 2022) analytic benchmark, and (iii) develops this verified capability within an open-source package.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Limitations of traditional tools for beyond design basis external hazard PRA

Probabilistic risk assessment (PRA) is being used increasingly by the nuclear industry for safety during normal operations as well as for the protection against external hazards. Computation of total risk in an external hazard PRA is dependent on hazard assessment, fragility assessment, and systems analysis. A systems analysis for propagation of component fragilities is conducted using event and fault trees. The event and fault trees for an actual power plant can be fairly large in size, which imposes computational challenges. Hence, certain assumptions are employed for computational efficiency. These assumptions typically represent the conditions imposed during the design basis (DB) scenario. The traditional PRA tools based on these assumptions are also widely applied to perform risk assessment in the context of beyond design basis (BDB) scenarios. However, some of these assumptions may not be valid for certain BDB scenarios. In addition, the probability of dependent failures also increases in BDB scenarios due to common cause failures (CCF) which usually results from design modifications, human errors, etc. In this manuscript, a simple and a relatively more complex illustrative examples are used to show the limitation of these assumptions in the numerical quantification of risk for the case of BDB conditions. Case studies with CCF events across multiple fault trees are also presented to illustrate the effect of these assumptions when traditional approach is used in BDB risk assessment. It is shown that the assumptions are valid for the case of DB conditions but may lead to excessively conservative risk estimates in the case of BDB conditions. Finally, a Bayesian network based top-down algorithm is proposed as an alternative tool for accurate numerical quantification of total risk in systems analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

A divergence-free constrained magnetic field interpolation method for scattered data

An interpolation method to evaluate magnetic fields, given its unstructured and scattered magnetic data, is presented. The method is based on the reconstruction of the global magnetic field using a superposition of orthogonal functions. The coefficients of the expansion are obtained by minimizing a cost function defined as the L2 norm of the difference between the ground truth and the reconstructed magnetic field evaluated on the training data. The divergence-free condition is incorporated as a constraint in the cost function, allowing the method to achieve arbitrarily small errors in the magnetic field divergence. An exponential decay of the approximation error is observed and compared with the less favorable algebraic decay of local splines. Compared to local methods involving computationally expensive search algorithms, the proposed method exhibits a significant reduction of the computational complexity of the field evaluation, while maintaining a small error in the divergence even in the presence of magnetic islands and stochasticity. Finally, applications to the computation of Poincaré sections using data obtained from numerical solutions of the magnetohydrodynamic equations in toroidal geometry are presented and compared with local methods currently in use.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Local step-flow dynamics in thin film growth with desorption

Desorption of deposited species plays a role in determining the evolution of surface morphology during crystal growth when the desorption time constant is short compared with the time to diffuse to a defect site, step edge, or kink. However, experiments to directly test the predictions of these effects are lacking. Novel techniques such as in situ coherent x-ray scattering can provide significant new information. Herein we present x-ray photon correlation spectroscopy (XPCS) measurements during diindenoperylene (DIP) vapor deposition on thermally oxidized silicon surfaces. DIP forms a nearly complete two-dimensional first layer over the range of temperatures studied (40–120 ⁢°C), followed by mounded growth during subsequent deposition. Local step flow within mounds was observed, and we find that there was a terrace-length-dependent behavior of the step edge dynamics. This led to unstable growth with rapid roughening (β>0.5) and deviation from a symmetric error-function-like height profile. At high temperatures, the grooves between the mounds tend to close up leading to nearly flat polycrystalline films. In conclusion, numerical analysis based on a (1+1)-dimensional model suggests that terrace-length dependent desorption of deposited ad-molecules is an essential cause of the step dynamics, and it influences the morphology evolution.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Model Predictive Control-Based Trajectory Shaper for Safe and Efficient Adaptive Cruise Control

Recent studies show that commercially-available adaptive cruise control (ACC) systems are string-unstable, indicating that ACC-driven vehicles amplify speed fluctuations from downstream traffic and induce stop-and-go waves. Moreover, it is challenging to revise the original control algorithm of an ACC system to achieve string stability due to its internal complexity and powertrain uncertainties. To achieve desired control performance given a string-unstable ACC system and circumvent revising the original control algorithm, this study proposes a model predictive control-based trajectory shaper (MPC-TS), which only modifies the sensor-measured trajectory information (i.e., position and speed) of the preceding vehicle. The proposed MPC-TS leverages the input shaping technique to generate reference trajectory to improve string stability, while incorporating tracking errors and vehicle acceleration/deceleration magnitude in the MPC cost function and constraining fluctuations of vehicle speed and spacing to ensure desired car-following performance. Numerical experiments validate the control performance of ACC with the proposed MPC-TS in terms of string stability, safety, traffic efficiency, and comfort.

Zhou, Anye↗
DOE OSTI2026

Demand Response Under Stochastic, Price-Dependent User Behavior

This letter focuses on price-based demand response (DR) implemented through dynamic adjustments of electricity prices. It extends existing DR models to a stochastic framework in which customer response is represented by price-dependent random variables, leveraging models and tools from the theory of stochastic optimization with decision-dependent distributions. The inherent epistemic uncertainty in the customers' responses renders open-loop, model-based DR strategies impractical. We propose a stochastic, feedback-based pricing strategy to compensate for estimation errors and uncertainty in customer response, establish theoretical results demonstrating the stability and near-optimality of the proposed approach, and validate its effectiveness through numerical simulations.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
NASA NTRS1971

Computational investigations of low-discrepancy point sets

The quasi-Monte Carlo method of integration offers an attractive solution to the problem of evaluating integrals in a large number of dimensions; however, the associated error bounds are difficult to obtain theoretically. Since these bounds are associated with the L2 discrepancy of the set of points used in the integration. Numerical calculations of the L2 discrepancy for several types of quasi-Monte Carlo formulae are presented.

Warnock, T. T.↗
NASA NTRS1972

Quantum-mechanical treatment of an electron undergoing synchrotron radiation.

The problem of an electron moving perpendicular to an intense magnetic field is approached from the framework of quantum mechanics. A numerical solution to the related rate equations describing the probabilities of occupation of the electron's energy states is put forth along with the expected errors involved. The quantum-mechanical approach is found to predict a significant amount of energy broadening with time for an initially monoenergetic electron beam entering a region of an intense magnetic field as long as the product of initial energy and magnetic field is of order 50 MG BeV or larger.

White, D.↗
NASA NTRS1972

Sound propagation in duct shear layers.

Consideration of the propagation of sound in a two-dimensional inviscid shear layer for given initial sound pressure profiles. The eigenvalue problem resulting from an assumption of separable solutions in a form first obtained by Pridmore-Brown (1958) is solved numerically. It is shown that although the resulting eigenfunctions cannot be proven to be orthogonal or complete, they can be combined by a least total error squared method to give a good representation of the initial pressure profile. The acoustic pressure in the duct is then easily calculated. The results verify all the predictions made in an earlier perturbation calculation. Moreover, they show that the refraction effect gets saturated at high subsonic Mach numbers and at high frequencies. The proposed technique may be used, if necessary, with impedance boundary conditions at the duct walls.

Shankar, P. N.↗
NASA NTRS1973

Computer checking of rotational line intensity factors for diatomic transitions.

Previously published analytical formulae for the rotational line intensity factors of diatomic molecules have been compared with results from a comprehensive computer program, which determines numerical intensity factors for both spin-allowed and spin-forbidden electric and magnetic dipole transitions in diatomic molecules. The comparison uncovered several typographical errors and a few algebraic errors in the published formulae. The changes required in the formulae to give agreement with the results from the computer program are tabulated.

Whiting, E. E.↗
NASA NTRS1975

Editing of EIA coded, numerically controlled, machine tool tapes

Editing of numerically controlled (N/C) machine tool tapes (8-level paper tape) using an interactive graphic display processor is described. A rapid technique required for correcting production errors in N/C tapes was developed using the interactive text editor on the IMLAC PDS-ID graphic display system and two special programs resident on disk. The correction technique and special programs for processing N/C tapes coded to EIA specifications are discussed.

Weiner, J. M.↗
NASA NTRS1977

Orbit determination accuracies using satellite-to-satellite tracking

The uncertainty in relay satellite sate is a significant error source which cannot be ignored in the reduction of satellite-to-satellite tracking data. Based on simulations and real data reductions, it is numerically impractical to use simultaneous unconstrained solutions to determine both relay and user satellite epoch states. A Bayesian or least squares estimation technique with an a priori procedure is presented which permits the adjustment of relay satellite epoch state in the reduction of satellite-to-satellite tracking data without the numerical difficulties introduced by an ill-conditioned normal matrix.

Vonbun, F. O.↗
NASA NTRS1979

HEAO 1 observations of 4U1232+07

An analysis of the NRL HEAO 1 observations of the Virgo cluster and 3C 273 has revealed X-ray emission from a region between these two objects. The most satisfactory explanation is the presence of two weak X-ray sources, for which error boxes have been obtained. One source is likely to be the peculiar galaxy IC 3576, suggested by Margon et al. (1972) as the source 2U 1231+07=4U 1232+07. The error box resulting from combining the present observations with the 4U data strengthen the X-ray identification with this extremely interesting galaxy. There are numerous possibilities for the other X-ray source, including the interacting pair of galaxies NGC 4410a and NGC 4410b (NGC 4410a/b).

Bhavsar, S. P.↗
NASA NTRS1981

Simulation of a numerical filter for enhancing earth radiation budget measurements

The Earth Radiation Budget Experiment has the objective to collect the radiation budget data which are needed to determine the radiation budget at the top of the atmosphere (TOA) on a regional scale. A second objective is to determine the accuracy of the results. Three satellites will carry wide and medium field of view radiometers which measure the longwave and shortwave components of radiation. Scanning radiometers will be included to detect small spatial features. A proposal has been made to employ for the nonscanning radiometers a one-dimensional numerical filter which reduces satellite measurements to TOA radiant excitances. The numerical filter was initially formulated by House (1980). It enhances the resolution of the radiation budget along the satellite groundtrack. The accuracy of the numerical filter estimate is studied by simulating the data gathering and measurement inversion process. The results of the study are discussed, taking into account two error sources.

Green, R. N.↗