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 1,261 records · Page 70

Navigational demands of low-level helicopter flight

The present study was designed to assess the verbal references to map, terrain, direction, and position, that navigators and pilots communicate during a simulated low-level flight. Two-person crews were asked to communicate normally while negotiating six short flight missions that varied widely in regard to map-terrain characteristics. Half of the crews performed the exercises with fixed, north-up maps and the other half used movable maps that were adjusted to maintain a track-up correspondence to the flight route. An analysis was performed to compare differences in crew communication patterns between map orientations, characteristics, and navigation tasks. The results showed differences in the frequency of communication across categories for map characteristics but not for map orientations or tasks. A difference in the proportion of communications between pilot and navigator occurred when crew were lost. Pilots communicated more than navigators when crews were lost and the reverse was true when crews were not lost. The results support and extend findings from a distance error performance analysis performed on the same flight missions.

Delzell, Susanne↗

Uncertainty analysis of signal deconvolution using a measured instrument response function

A common analysis procedure minimizes the ln-likelihood that a set of experimental observables matches a parameterized model of the observation. The model includes a description of the underlying physical process as well as the instrument response function (IRF). Here, we investigate the National Ignition Facility (NIF) neutron time-of-flight (nTOF) spectrometers, the IRF is constructed from measurements and models. IRF measurements have a finite precision that can make significant contributions to the uncertainty estimate of the physical model’s parameters. Finally, we apply a Bayesian analysis to properly account for IRF uncertainties in calculating the ln-likelihood function used to find the optimum physical parameters.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Gaussian Process for Flight Delay Prediction: Learning a Stochastic Process

This paper presents a machine-learning approach to predict flight delays. Whereas neural networks are extensively studied for predictive capabilities, they involve non-intuitive design and extensive analysis, particularly in training and optimization processes. Instead, the proposed framework employs Gaussian Processes as a supervised learning technique for flight delay prediction. This data-driven approach trains the model using prior information, specifically the mean and covariance tied to existing data. The proposed Gaussian Process Regression (GPR) model employs the day of flight as a pivotal feature for delay forecasting. We analyze flights from various routes and gauge the accuracy of the presented learning technique by comparing the predicted delays with the actual ones. Given the inherent challenges in precisely forecasting delays, we predict the delays with a 95 % confidence interval. Also, an error propagation analysis in the prediction horizon is carried out to determine the optimal time frame for prediction. The proposed method for flight delay prediction is important as airlines can strategize flight operations and issue timely advisories.

stochastic↗

The effects of mapping CT images to Monte Carlo materials on GEANT4 proton simulation accuracy

Purpose: Monte Carlo simulations of radiation therapy require conversion from Hounsfield units (HU) in CT images to an exact tissue composition and density. The number of discrete densities (or density bins) used in this mapping affects the simulation accuracy, execution time, and memory usage in GEANT4 and other Monte Carlo code. The relationship between the number of density bins and CT noise was examined in general for all simulations that use HU conversion to density. Additionally, the effect of this on simulation accuracy was examined for proton radiation. Methods: Relative uncertainty from CT noise was compared with uncertainty from density binning to determine an upper limit on the number of density bins required in the presence of CT noise. Error propagation analysis was also performed on continuously slowing down approximation range calculations to determine the proton range uncertainty caused by density binning. These results were verified with Monte Carlo simulations. Results: In the presence of even modest CT noise (5 HU or 0.5%) 450 density bins were found to only cause a 5% increase in the density uncertainty (i.e., 95% of density uncertainty from CT noise, 5% from binning). Larger numbers of density bins are not required as CT noise will prevent increased density accuracy; this applies across all types of Monte Carlo simulations. Examining uncertainty in proton range, only 127 density bins are required for a proton range error of <0.1 mm in most tissue and <0.5 mm in low density tissue (e.g., lung). Conclusions: By considering CT noise and actual range uncertainty, the number of required density bins can be restricted to a very modest 127 depending on the application. Reducing the number of density bins provides large memory and execution time savings in GEANT4 and other Monte Carlo packages.

60 APPLIED LIFE SCIENCES↗

Error propagation in a digital avionic processor: A simulation-based study

An experimental analysis to study error propagation from the gate to the chip level is described. The target system is the CPU in the Bendix BDX-930, an avionic miniprocessor. Error activity data for the study was collected via a gate-level simulation. A family of distributions to characterize the error propagation, both within the chip and at the pins, was then generated. Based on these distributions, measures of error propagation and severity were defined. The analysis quantifies the dependency of the measured error propagation on the location of the fault and the type of instruction/microinstruction executed.

Lomelino, D.↗

THE NANOGRAV NINE-YEAR DATA SET: MONITORING INTERSTELLAR SCATTERING DELAYS

We report on an effort to extract and monitor interstellar scintillation parameters in regular timing observations collected for the North American Nanohertz Observatory for Gravitational Waves pulsar timing array. Scattering delays are measured by creating dynamic spectra for each pulsar and observing epoch of wide-band observations centered near 1500 MHz and carried out at the Green Bank Telescope and the Arecibo Observatory. The ∼800 MHz wide frequency bands imply dramatic changes in scintillation bandwidth across the bandpass, and a stretching routine has been included to account for this scaling. For most of the 10 pulsars for which the scaling has been measured, the bandwidths scale with frequency less steeply than expected for a Kolmogorov medium. We find estimated scattering delay values that vary with time by up to an order of magnitude. The mean measured scattering delays are similar to previously published values and are slightly higher than predicted by interstellar medium models. We investigate the possibility of increasing the timing precision by mitigating timing errors introduced by the scattering delays. For most of the pulsars, the uncertainty in the time of arrival of a single timing point is much larger than the maximum variation of the scattering delay, suggesting that diffractive scintillation remains as only a negligible part of their noise budget.

79 ASTRONOMY AND ASTROPHYSICS↗

Evaluating E. coli genome‐scale metabolic model accuracy with high‐throughput mutant fitness data

Abstract The Escherichia coli genome‐scale metabolic model (GEM) is an exemplar systems biology model for the simulation of cellular metabolism. Experimental validation of model predictions is essential to pinpoint uncertainty and ensure continued development of accurate models. Here, we quantified the accuracy of four subsequent E. coli GEMs using published mutant fitness data across thousands of genes and 25 different carbon sources. This evaluation demonstrated the utility of the area under a precision–recall curve relative to alternative accuracy metrics. An analysis of errors in the latest (iML1515) model identified several vitamins/cofactors that are likely available to mutants despite being absent from the experimental growth medium and highlighted isoenzyme gene‐protein‐reaction mapping as a key source of inaccurate predictions. A machine learning approach further identified metabolic fluxes through hydrogen ion exchange and specific central metabolism branch points as important determinants of model accuracy. This work outlines improved practices for the assessment of GEM accuracy with high‐throughput mutant fitness data and highlights promising areas for future model refinement in E. coli and beyond.

59 BASIC BIOLOGICAL SCIENCES↗

Michigan Orbital DEbris Survey Telescope Observations of the Geosynchronous Orbital Debris Environment. Observing Years: 2007-2009

NASA uses the Michigan Orbital DEbris Survey Telescope (MODEST), the University of Michigan's 0.61-m aperture Curtis-Schmidt telescope at the Cerro Tololo Inter-American Observatory in Chile, to help characterize the debris environment in geosynchronous orbit; this began in February 2001 and continues to the present day. Detected objects that are found to be on the U.S. Space Surveillance Network cataloged objects list are termed correlated targets (CTs), while those not found on the list are called uncorrelated targets (UCTs). This Johnson Space Center report provides details of observational and data-reduction processes for the entire MODEST dataset acquired in calendar years (CYs) 2007, 2008, and 2009. Specifically, this report describes the collection and analysis of 36 nights of data collected in CY 2007, 43 nights of data collected in CY 2008, and 43 nights of data collected in CY 2009. MODEST is equipped with a 2048 x 2048-pixel charged coupled device camera with a 1.3 by 1.3 deg field of view. This system is capable of detecting objects fainter than 18th magnitude (R filter) using a 5-s integration. This corresponds to a 20-cm diameter, 0.175-albedo object at 36,000 km altitude assuming a diffuse Lambertian phase function. The average number of detections each night over all 3 years was 26. The percentage of this number that represented the UCT population ranged from 34% to 18%, depending on the observing strategy and the field center location. Due to the short orbital arc over which observations are made, the eccentricity of the object s orbit is extremely difficult to measure accurately. Therefore, a circular orbit was assumed when calculating the orbital elements. A comparison of the measured inclination (INC), right ascension of ascending node (RAAN), and mean motion to the quantities for CTs from the U.S. Space Surveillance Network shows acceptable errors. This analysis lends credibility to the determination of the UCT orbital distributions. Figure 1 shows the size distribution of 3,143 objects detected in the data processed for CYs 2007, 2008, and 2009. The actual peak of the absolute magnitude distribution for the functional correlated targets is 10th magnitude, whereas the peak was 11th magnitude in 2002 2003 and 10th magnitude for 2004-2006. An absolute magnitude of 10.5 corresponds to objects with average diameters of 6.3 m, assuming an albedo of 0.175 and a diffuse Lambertian phase function. This result generally agrees with the known sizes of intact satellites. The absolute magnitude distribution for the UCTs is broad, but starts to roll off near 25 cm diameter or 17.5 magnitude. This roll off in the distribution reflects the detection capability of MODEST, not the true nature of the population. The true population is believed to continue at the same slope through fainter magnitudes based on comparisons with the LEO break-up law.

Abercromby, K. J.↗

New experiment WAGASCI to measure cross sections of neutrino interactions in water and hydrocarbon using J-PARC beam

The main objective of the long base T2K,Tokai-to-Kamioka, Japan, experiment is the investigation of neutrino oscillations with a high-intensity beam of muon neutrinos and antineutrinos. Neutrino interaction measurements in the near detector complex (ND280) of the Т2К experiment are used for adjustment of the parameters defining event simulation in the far detector (Super-Kamiokande), which is the key idea of the experiment; this allows one to considerably increase the oscillation analysis accuracy. The difference in the target material of the far (water) and near (scintillator, hydrocarbon) detectors results in a substantial systematic error in oscillation analysis. Systematic error can be reduced by direct measurements, a new water grid and scintillator detector WAGASCI was proposed for this purpose. The detector will operate at the J-PARC neutrino beam and measure the ratio of charged-current neutrino cross sections between water and hydrocarbon with a few percent precision. Further physical programs may include precise measurements of various charged channels of neutrino interaction. This paper presents the detector concept and the current plan for its installation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

An Analysis of Gravity-Field Estimation Based on Intersatellite Dual-1-Way Biased Ranging

The GRACE (Gravity Recovery And Climate Experiment) mission is designed to make global, highly accurate measurements of the Earth's gravity field with high spatial resolution. Ancillary GPS occultation measurements are also to be carried out for atmospheric monitoring. In the dual-1-way biased ranging of this mission, the range between two satellites separated by 100 to 200 km in nearly polar, coplanar, circular orbits, is measured to very high precision, to within an additive constant, through the exchange of K- and Ka-band sinusoidal signals. Such biased ranging data, along with GPS L-band range and phase data, can be processed and fit over successive multiday intervals to obtain accurate estimates of the Earth's gravity field. This report approximately models and analyzes this process, from the generation of the RF signals at the two satellites through the extraction of the geopotential. The steps include generation of the transmitted signals, processing the received signals to extract high-rate baseband phase, carrying out a dual-1-way combination of baseband phase to extract high-rate biased range for each band, combining K- and Ka-band ranges to correct for the ionosphere effect, and processing the resulting high-rate biased range values to extract three types of reduced-rate observables: biased range, range rate and range acceleration. The version of dual-1-way biased ranging developed by this report improves upon previous versions in a number of ways: highly accurate satellite-timetag corrections derived from concurrent GPS data, better baseband phase extraction using highly digital processing, highly accurate USO-rate calibration derived from concurrent GPS data, an improved method for extracting high-rate biased range from baseband phase, improved filtering for extracting reduced- rate observables from high-rate biased range, and parallel extraction of three observable types.

SINUSOIDAL SIGNALS↗

Mixed precision s –step Lanczos and conjugate gradient algorithms

Compared to the classical Lanczos algorithm, the s-step Lanczos variant has the potential to improve performance by asymptotically decreasing the synchronization cost per iteration. However, this comes at a price; despite being mathematically equivalent, the s-step variant may behave quite differently in finite precision, potentially exhibiting greater loss of accuracy and slower convergence relative to the classical algorithm. It has previously been shown that the errors in the s-step version follow the same structure as the errors in the classical algorithm, but are amplified by a factor depending on the square of the condition number of the O(s)-dimensional Krylov bases computed in each outer loop. As the condition number of these s-step bases grows (in some cases very quickly) with s, this limits the s values that can be chosen and thus can limit the attainable performance. In this work, we show that if a select few computations in s-step Lanczos are performed in double the working precision, the error terms then depend only linearly on the conditioning of the s-step bases. This has the potential for drastically improving the numerical behavior of the algorithm with little impact on per-iteration performance. Our numerical experiments demonstrate the improved numerical behavior possible with the mixed precision approach, and also show that this improved behavior extends to mixed precision s-step CG. Here, we present preliminary performance results on NVIDIA V100 GPUs that show that the overhead of extra precision is minimal if one uses precisions implemented in hardware.

97 MATHEMATICS AND COMPUTING↗

Data-Driven Multi-Step Demand Prediction for Ride-Hailing Services Using Convolutional Neural Network

Ride-hailing services are growing rapidly and becoming one of the most disruptive technologies in the transportation realm. Accurate prediction of ride-hailing trip demand not only enables cities to better understand people's activity patterns, but also helps ride-hailing companies and drivers make informed decisions to reduce deadheading vehicle miles traveled, traffic congestion, and energy consumption. In this study, a convolutional neural network (CNN)-based deep learning model is proposed for multi-step ride-hailing demand prediction using the trip request data in Chengdu, China, offered by DiDi Chuxing. The CNN model is capable of accurately predicting the ride-hailing pick-up demand at each 1-km by 1-km zone in the city of Chengdu for every 10 min. Compared with another deep learning model based on long short-term memory, the CNN model is 30% faster for the training and predicting process. The proposed model can also be easily extended to make multi-step predictions, which would benefit the on-demand shared autonomous vehicles applications and fleet operators in terms of supply-demand rebalancing. The prediction error attenuation analysis shows that the accuracy stays acceptable as the model predicts more steps.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Computationally inexpensive part-scale thermal history of additive friction-stir deposition

This study presents an analytical model for steady-state power generation and tool heat loss in additive friction-stir deposition (AFSD), developed to enable part-scale thermal simulation while remaining computationally inexpensive. The model predicts total generated power, yielding 3.7–4.7 kW across deposition temperature setpoints of 400–460 °C for the deposition of AA6061 with a Be-Cu tool. This corresponds to 90–95% of the reported spindle power. Tool heat loss is experimentally determined by calibrating a steady-state energy balance between the generated power, the substrate-deposition thermal gradient, and a temperature dependent tool heat loss term: q tool (T) = a + b (T - 400°C) with a = 2.7 x 10 6 Wm -2 and b = 9.5 x 10 3 Wm -2 K -1 . The calibration indicates that about 69% of the generated heat is conducted into the tool for this configuration, which is much higher than previously reported. The calibrated heat-source is implemented in finite element software (Adamantine) to simulate the transient thermal history of a 100 cm 3 representative build in 8 min on a standard desktop (at 0.635 mm build-height resolution). For the first three layers, the substrate temperatures between simulation and experiment are within 10% mean absolute percentage error. Sensitivity analysis indicates that uncertainties in average deposition temperature and deformation localization (stir-zone geometry, depth, and spatial dependance of strain-rate and flow stress) dominate model variance, motivating additional experimental verification.

Additive Friction-Stir Deposition↗

A Mixed Discontinuous Galerkin Method for Linear Elasticity with Strongly Imposed Symmetry

In this paper, we study a mixed discontinuous Galerkin (MDG) method to solve linear elasticity problem with arbitrary order discontinuous finite element spaces in d -dimension ( \(d=2,3\) ). This method uses polynomials of degree \(k+1\) for the stress and of degree k for the displacement ( \(k\ge 0\) ). The mixed DG scheme is proved to be well-posed under proper norms. Specifically, we prove that, for any \(k \ge 0\) , the \(H(\mathrm{div})\) -like error estimate for the stress and \(L^2\) error estimate for the displacement are optimal. We further establish the optimal \(L^2\) error estimate for the stress provided that the \({\mathcal {P}}_{k+2}-{\mathcal {P}}_{k+1}^{-1}\) Stokes pair is stable and \(k \ge d\) . Finally, we also provide numerical results of MDG showing that the orders of convergence are actually sharp.

97 MATHEMATICS AND COMPUTING↗

A $C^1$-Conforming Arbitrary-Order Two-Dimensional Virtual Element Method for the Fourth-Order Phase-Field Equation

We present a two-dimensional conforming virtual element method for the fourth-order phase-field equation. Our proposed numerical approach to the solution of this high-order phase-field (HOPF) equation relies on the design of an arbitrary-order accurate, virtual element space with $C^1$ global regularity. Such regularity is guaranteed by taking the values of the virtual element functions and their full gradient at the mesh vertices as degrees of freedom. Attaining high-order accuracy requires also edge polynomial moments of the trace of the virtual element functions and their normal derivatives. In this work, we detail the scheme construction, and prove its convergence by deriving error estimates in different norms. A set of representative test cases allows us to assess the behavior of the method.

97 MATHEMATICS AND COMPUTING↗

An efficient approach for solving Klein–Gordon equation arising in quantum field theory using wavelets

A new approach for the solution of Klein–Gordon equation using Legendre wavelet-based approximation method is presented. The properties of Legendre wavelets are used to reduce the problem to the solution of system of algebraic equations. Usually 2{sup (k−1)}M{sup 2} connection coefficients are required to solve the Klein–Gordon equation by Legendre wavelets method as mentioned in the literature. But, our proposed method based on Legendre wavelets and algebraic polynomials require only 2{sup (k−1)}M connection coefficients instead of 2{sup (k−1)}M{sup 2} . Also the convergence analysis and error estimation for the proposed function approximation through the truncated series have been discussed and approved with the exact solution. Illustrative examples are discussed to demonstrate the validity and applicability of the technique with lesser computational effort.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Towards efficient backward-in-time adjoint computations using data compression techniques

In the context of a posteriori error estimation for nonlinear time-dependent partial differential equations, the state-of-the-practice is to use adjoint approaches which require the solution of a backward-in-time problem defined by a linearization of the forward problem. One of the major obstacles in the practical application of these approaches, we found, is the need to store, or recompute, the forward solution to define the adjoint problem and to evaluate the error representation. Our study considers the use of data compression techniques to approximate forward solutions employed in the backward-in-time integration. The development derives an error representation that accounts for the difference between the standard-approach and the compressed approximation of the forward solution. This representation is algorithmically similar to the standard representation and only requires the computation of the quantity of interest for the forward solution and the data-compressed reconstructed solution (i.e. scalar quantities that can be evaluated as the forward problem is integrated). This approach is then compared with existing techniques, such as checkpointing and time-averaged adjoints. Lastly, we provide numerical results indicating the potential efficiency of our approach on a transient diffusion–reaction equation and on the Navier–Stokes equations. These results demonstrate memory compression ratios up to 450×450× while maintaining reasonable accuracy in the error-estimates.

42 ENGINEERING↗