Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “error propagation”

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 379 records · Page 21

Performance of a 2-micrometer coherent Doppler lidar for wind measurements

Measurements of boundary layer winds are presented using a 2-micrometer coherent Doppler lidar and the optimal performance of the maximum likelihood estimator. The systematic error for single-shot estimates was estimated as 3.6 cm/s using measurements from a stationary hard target. The estimation error for measurements of the radial component of the wind field was determined, as well as the fraction of the estimates that are randomly distributed over the velocity search space, when the signal power is low and speckle fading is important. The results from actual data are compared with the results from ideal simulations. The first direct estimation of the spatial structure function of the radial wind field and of the energy dissipation rate is presented for both horizontal and vertical directions of propagation. The rms estimation error of the velocity estimates is found to be within 30% of ideal performance based on simulation.

Frehlich, Rod↗

Effect of ephemeris errors on the accuracy of the computation of the tangent point altitude of a solar scanning ray as measured by the SAGE 1 and 2 instruments

An analysis was made of the error in the minimum altitude of a geometric ray from an orbiting spacecraft to the Sun. The sunrise and sunset errors are highly correlated and are opposite in sign. With the ephemeris generated for the SAGE 1 instrument data reduction, these errors can be as large as 200 to 350 meters (1 sigma) after 7 days of orbit propagation. The bulk of this error results from errors in the position of the orbiting spacecraft rather than errors in computing the position of the Sun. These errors, in turn, result from the discontinuities in the ephemeris tapes resulting from the orbital determination process. Data taken from the end of the definitive ephemeris tape are used to generate the predict data for the time interval covered by the next arc of the orbit determination process. The predicted data are then updated by using the tracking data. The growth of these errors is very nearly linear, with a slight nonlinearity caused by the beta angle. An approximate analytic method is given, which predicts the magnitude of the errors and their growth in time with reasonable fidelity.

Buglia, James J.↗

Uncertainty Quantification of Expected Time-of-Arrival in UAV Flight Trajectory

One of the foremost requirements for accurate in-flight safety monitoring of autonomous unmanned aerial vehicles (UAVs) is tracking of their flight trajectory. Existing UAVs leverage autonomous flight functionalities based on trajectory generation algorithms developed in robotic applications such as polynomial or spline curves in order to facilitate kinematic smoothness, minimum vibrations and fuel efficiency. However in practice, the actual path may be subjected to unexpected local weather conditions, unexpected obstacles along the path or abrupt traffic changes in the low-altitude airspace resulting in large errors of the predicted time-of-arrival at way-points. In this study, an approach to quantify and propagate uncertainty in 4D trajectories is proposed. The paper presents a simple error interval propagation method based on the expected cruise speed of the UAV and its associated uncertainty. The uncertainty is then propagated in time to estimate reasonable confidence intervals on the times-of-arrival of the vehicle at each way-point as well as along the entire flight-path. The uncertainty propagation is demonstrated on a state-of-the-art trajectory generation algorithm based on non-uniform rational B-spline (NURBS) curves. Further, the effect of a stationary wind field is incorporated in the uncertainty propagation approach. The proposed method is implemented on synthetic and real data obtained from flight experiments with a small UAV.

Uncertainty Quantification↗

Advanced Finite-Volume Numerics and Source Term Assumptions for Kernel and G-Equation Modelling of Propane/Air Flames

Here G-Equation models represent propagating flame fronts with an implicit two-dimensional surface representation (level-set). Level-set methods are fast, as transport source terms for the implicit surface can be solved with finite-volume operators on the finite-volume domain, without having to build the actual surface. However, they include approximations whose practical effects are not properly understood. In this study, we improved the numerics of the FRESCO CFD code’s G-Equation solver and developed a new method to simulate kernel growth using signed distance functions and the analytical sphere-mesh overlap. We analyzed their role for simulating propane/air flames, using three well-established constant-volume configurations: a one-dimensional, freely propagating laminar flame; a disc-shaped, constant-volume swirl combustor; and torch-jet flame development through an orifice from a two-chamber device. We tested the explicit (sub-cycled) vs. implicit formulation for the standard transport operators (advection, diffusion, compressibility). In addition to the accurate flame swept-volume method for chemistry and species source term, we developed a more accurate estimator for the burnt/unburnt split cell composition. Then, we developed a signed-distance-function (SDF) based method which provides a more stable reinitialization of the level-set field at every time-step. We found that simplifying assumptions common to several G-Equation implementations, for straightforward terms such as compressibility and advection, lead to large errors in predicting the propagation of even laminar flames, with deviations up to ~300% in simulated vs. formulated flame speed. Conversely, the enhanced numerics enabled through the SDF field reinitialization and improved chemistry source term improve simulation stability and smooth flame propagation even with significantly larger solver time-steps.

42 ENGINEERING↗

Efficient data reduction for time-of-flight neutron scattering experiments on single crystals

Event-mode data collection presents remarkable new opportunities for time-of-flight neutron scattering studies of collective excitations, diffuse scattering from short-range atomic and magnetic structures, and neutron crystallography. In these experiments, large volumes of the reciprocal space are surveyed, often using different wavelengths and counting times. These data then have to be added together, with accurate propagation of the counting errors. This paper presents a statistically correct way of adding and histogramming the data for single-crystal time-of-flight neutron scattering measurements. In order to gain a broader community acceptance, particular attention is given to improving the efficiency of calculations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The effect of initial flaw configuration on leak critical pressure vessel design

The effect of initial flaw shape and depth on the cycles required to produce a pressure leak condition in parent and welded 2219-T87 aluminum was examined for thicknesses from 0.04 to 0.3 inch. Initial flaw aspect ratios of 0.5, 0.3, and 0.15 were examined for tests conducted at a frequency of 4 Hz at room temperature and at 300 F. A computer program developed by Forman was used to predict the cycles to leak curve as a function of initial crack aspect ratio, material thickness, and specimen type and the results compared with the experimental data. Sources of error in the crack propagation analysis and the interface with NDI capabilities are discussed.

Pettit, D. E.↗

On the effect of quarter-point element size on fracture criteria

The effect of the quarter-point element size on the solution parameters of fracture criteria is investigated in view of its effect on the computed stress intensity factors. An analytical relationship between the error in calculating crack propagation increment and stress intensity for single mode problems is obtained while, for mixed mode problems, the investigation is based on numerical experimentation. The singular element size is found to have a substantial effect on both crack propagation angle and increment for the geometries considered.

Yehia, N. A. B.↗

Survey of Primary Flow Measurement Parameters at the NASA Langley Transonic Dynamics Tunnel

An assessment of the methods and locations used to measure the primary flow conditions in the NASA Langley Transonic Dynamics Tunnel was conducted during calibration activities following the facility conversion from a Freon-12 heavy-gas test medium to R-134a. A survey of stagnation pressure, plenum static pressure, and stagnation temperature was undertaken at many pertinent locations in the settling chamber, plenum, and contraction section of the wind tunnel and these measurements were compared to those of the existing primary flow measurement systems. Local flow velocities were measured in the settling chamber using a pitot probe. Results illustrate that small discrepancies exist between measured primary tunnel flow conditions and the survey measurements. These discrepancies in tunnel stagnation pressure, plenum pressure, and stagnation temperature were found to be approximately +/- 1-3 psf and 2-3 degrees Fahrenheit. The propagation of known instrument errors in measured primary flow conditions and its impact on tunnel Mach number, dynamic pressure, flow velocity, and Reynolds number have been investigated analytically and shown to require careful attention when considering the uncertainty in measured test section conditions.

Piatak, David J.↗

Shuttle Debris Impact Tool Assessment Using the Modern Design of Experiments

Computational tools have been developed to estimate thermal and mechanical reentry loads experienced by the Space Shuttle Orbiter as the result of cavities in the Thermal Protection System (TPS). Such cavities can be caused by impact from ice or insulating foam debris shed from the External Tank (ET) on liftoff. The reentry loads depend on cavity geometry and certain Shuttle state variables, among other factors. Certain simplifying assumptions have been made in the tool development about the cavity geometry variables. For example, the cavities are all modeled as shoeboxes , with rectangular cross-sections and planar walls. So an actual cavity is typically approximated with an idealized cavity described in terms of its length, width, and depth, as well as its entry angle, exit angle, and side angles (assumed to be the same for both sides). As part of a comprehensive assessment of the uncertainty in reentry loads estimated by the debris impact assessment tools, an effort has been initiated to quantify the component of the uncertainty that is due to imperfect geometry specifications for the debris impact cavities. The approach is to compute predicted loads for a set of geometry factor combinations sufficient to develop polynomial approximations to the complex, nonparametric underlying computational models. Such polynomial models are continuous and feature estimable, continuous derivatives, conditions that facilitate the propagation of independent variable errors. As an additional benefit, once the polynomial models have been developed, they require fewer computational resources to execute than the underlying finite element and computational fluid dynamics codes, and can generate reentry loads estimates in significantly less time. This provides a practical screening capability, in which a large number of debris impact cavities can be quickly classified either as harmless, or subject to additional analysis with the more comprehensive underlying computational tools. The polynomial models also provide useful insights into the sensitivity of reentry loads to various cavity geometry variables, and reveal complex interactions among those variables that indicate how the sensitivity of one variable depends on the level of one or more other variables. For example, the effect of cavity length on certain reentry loads depends on the depth of the cavity. Such interactions are clearly displayed in the polynomial response models.

DeLoach, R.↗

A Critical Look at Deriving Monthly Aerosol Optical Depth from Satellite Data

Satellite-derived aerosol data sets, such as those provided by NASA's Moderate Resolution Imaging Spectroradiometer (MODIS) instruments, are greatly improving our understanding of global aerosol optical depth (AOD). Yet, there are sampling issues. MODIS specific orbital geometry, convolved with the need to avoid bright surfaces (glint, desert, clouds, etc.), means that AOD can be under- or over-sampled in places. When deriving downstream products, such as daily or monthly gridded AOD, one must consider the spatial and temporal density of the measurements relative to the gradients of the true AOD. Additionally, retrieval confidence criteria should be considered. Averaged products are highly dependent on choices made for data aggregation and weighting, and sampling errors can be further propagated when deriving regional or global mean AOD. Different choices for aggregation and weighting result in estimates of regional and global means varying by 30% or more. The impacts of a particular averaging algorithm vary by region and surface type and can be shown to represent different tolerance for clouds and retrieval confidence.

Levy, R. C.↗

An Error Analysis of the Phased Array Antenna Pointing Algorithm for STARS Flight Demonstration No. 2

STARS is a multicenter NASA project to determine the feasibility of using space-based assets, such as the Tracking and Data Relay Satellite System (TDRSS) and Global Positioning System (GPS), to increase flexibility (e.g. increase the number of possible launch locations and manage simultaneous operations) and to reduce operational costs by decreasing the need for ground-based range assets and infrastructure. The STARS project includes two major systems: the Range Safety and Range User systems. The latter system uses broadband communications (125 kbps to 500 kbps) for voice, video, and vehicle/payload data. Flight Demonstration #1 revealed the need to increase the data rate of the Range User system. During Flight Demo #2, a Ku-band antenna will generate a higher data rate and will be designed with an embedded pointing algorithm to guarantee that the antenna is pointed directly at TDRS. This algorithm will utilize the onboard position and attitude data to point the antenna to TDRS within a 2-degree full-angle beamwidth. This report investigates how errors in aircraft position and attitude, along with errors in satellite position, propagate into the overall pointing vector.

Carney, Michael P.↗

Assessing Uncertainties in Gridded Emissions: A Case Study for Fossil Fuel Carbon Dioxide (FFCO2) Emission Data

Fossil fuel carbon dioxide (CO2) emissions (FFCO2) are the largest input to the global carbon cycle on a decadal time scale. Because total emissions are assumed to be reasonably well constrained by fuel statistics, FFCO2 often serves as a reference in order to deduce carbon uptake by poorly understood terrestrial and ocean sinks. Conventional atmospheric CO2 flux inversions solve for spatially explicit regional sources and sinks and estimate land and ocean fluxes by subtracting FFCO2. Thus, errors in FFCO2 can propagate into the final inferred flux estimates. Gridded emissions are often based on disaggregation of emissions estimated at national or regional level. Although national and regional total FFCO2 are well known, gridded emission fields are subject to additional uncertainties due to the emission disaggregation. Assessing such uncertainties is often challenging because of the lack of physical measurements for evaluation. We first review difficulties in assessing uncertainties associated with gridded FFCO2 emission data and present several approaches for evaluation of such uncertainties at multiple scales. Given known limitations, inter-emission data differences are often used as a proxy for the uncertainty. The popular approach allows us to characterize differences in emissions, but does not allow us to fully quantify emission disaggregation biases. Our work aims to vicariously evaluate FFCO2 emission data using atmospheric models and measurements. We show a global simulation experiment where uncertainty estimates are propagated as an atmospheric tracer (uncertainty tracer) alongside CO2 in NASA's GEOS model and discuss implications of FFCO2 uncertainties in the context of flux inversions. We also demonstrate the use of high resolution urban CO2 simulations as a tool for objectively evaluating FFCO2 data over intense emission regions. Though this study focuses on FFCO2 emission data, the outcome of this study could also help improve the knowledge of similar gridded emissions data for non-CO2 compounds with similar emission characteristics.

FFCO↗

Experiment Concepts in The Trans-Pacific HDR Satcom Experiment - Phase 2

This paper describes experiment concepts for the Phase-2 of the Trans-Pacific HDR Satcom Experiment (Astronomy). Broadband Satellites serve as a key link in establising a worldwide network of networks, while imposing upon us the particular charcteristics of propagation delay, bandwidth, and error rates which differ from those of terrestrial networks.

Trans-Pacific↗

Multiphenomenology explosion monitoring (MultiPEM): a general framework for data interpretation and yield estimation

SUMMARY An underground nuclear explosion (UNE) couples mechanical energy into crustal rock, which propagates as seismic and acoustic waves. These different physical phenomena transport, by different pathways, to standoff detectors at varying distances. The transport pathways attenuate the original signal but in different ways. Enabled by correct statistical weighting, signal attenuation models can be used to combine these disparate sensor data to estimate the yield of an UNE. Contemporaneous statistical models, used in yield estimation, can be improved with an advanced partition of error for these physical signal propagation models. We present an advanced multivariate approach to error modelling of multiphenomenology physical signatures. In addition to measurement error, our error model represents physical model biases as random with a physics-based covariance structure. To illustrate this proposed framework, we demonstrate the estimation of explosion yield using openly available seismic and acoustic data from chemical single-point explosions.

Williams, Brian J.↗

An evaluation of soundings, analyses and model forecasts derived from TIROS-N and NOAA-6 satellite data

TIROS-N and NOAA-6 temperature soundings over North America during three days in January 1980, and synoptic analyses and numerical-model forecasts derived from them, are compared with conventional data and analyses from NMC's limited-area fine-mesh model (LFM). The collocated sounding comparison revealed significant errors, especially near the surface and the tropopause. Satellite-derived thermal gradients were found to be weak, and thickness-analysis difference fields to propagate eastward, suggesting that sounding errors are correlated with synoptic patterns. The same pattern of anomalies is seen in the model forecasts. More detailed determinations of the correlation detected here could be used to optimize the assimilation of satellite soundings to conventional data.

Koehler, T. L.↗

Two-stage formation-energy correction (NbZr, TaZr, VZr)

This bundle contains the scripts, the raw and corrected per-structure data, and the manuscript plots for the NbZr / TaZr / VZr BCC binary formation energies and the associated RMSDs. Why a two-stage correction is necessary: The "raw" formation energy of every relaxed VASP configuration is computed in the usual way, FE_raw(c) = E_alloy(c) - sum_i x_i * E_pure_i , where E_pure_i are the per-atom total energies of the pure-element reference structures (Nb, Ta, V, Zr in the same BCC supercell, with identical INCAR / KPOINTS / PAW choices). With perfectly consistent reference runs the raw FE should vanish at the two pure-element endpoints (x = 0 and x = 1) by construction. In practice this does not hold for two reasons that are present in our dataset: 1. Reference-energy inconsistency (composition-dependent bias). Even with identical input parameters, the pure-element runs (stored in `corrected_DFT_pure_element_runs/`) differ slightly from the values that would be implied by the alloy runs at near-pure compositions (a few meV/atom). This bias is approximately linear in concentration, because the residual error in E_pure_Nb (or E_pure_Ta / E_pure_V) propagates into FE_raw(c) as (1 - x) * dE_pure_1, and the corresponding error in E_pure_Zr propagates as x * dE_pure_2. Left uncorrected, this produces a non-physical "tilt" of FE_raw(x) and shifts the entire FE-vs-x cloud away from zero at the endpoints. 2. Endpoint anchoring against the audited true endpoints. The strict endpoint values (FE_x0_meVatom, FE_x1_meVatom in `corrected_fe_strict_endpoints_20260518/strict_endpoint_check_20260518.csv`) were re-derived from an independent cross-check of the pure-element runs. After stage 1 removes the linear bias, the near-pure compositions in the alloy dataset still extrapolate to values that differ slightly from these audited endpoints — because stage 1 is fit from a few near-end alloy bins, not from the audited pure-element references themselves. The README.txt file discusses how these issues are addressed by the two-stage correction, and describes folder layout, pipeline summary, and how to re-run.

36 MATERIALS SCIENCE↗

Errors in reduction methods

A mathematical basis is given for comparing the relative merits of various techniques used to reduce the order of large linear and nonlinear dynamics problems during their numerical integration. In such techniques as Guyan-Irons, path derivatives, selected eigenvectors, Ritz vectors, etc., the nth order initial value problem of /y(dot) = f(y) for t greater than 0, y(0) given/ is typically reduced to the mth order (m is much less than n) problem of /z(dot) = g(z) for t greater than 0, z(0) given/ by the transformation y = Pz where P changes from technique to technique. This paper gives an explicit approximate expression for the reduction error e-i in terms of P and the Jacobian of f. It is shown that: (a) reduction techniques are more accurate when the time rate of change of the response y is relatively small; (b) the change in response between two successive stations contributes to the errors at future stations after the change in response is transformed by a filtering matrix H, defined in terms of P; (c) the error committed at a station propagates to future stations by a mixing and scaling matrix G, defined in terms of P, Jacobian and of f, and time increment h. The paper discusses the conditions under which the reduction errors may be minimized and gives guidelines for selecting the reduction basis vector, i.e., the columns of P.

Utku, S.↗