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 649 records · Page 36

HRDI Observations of Inertia-Gravity Waves in the Mesosphere and Lower Thermosphere

Vertical profiles of High-resolution Doppler imager (HRDI) mesospheric winds have small-scale structure (vertical wavelengths between 10 and 20 km) that is virtually always present. Fourier analysis of HRDI zonal and meridional wind profiles have been carried out, and the spectral characteristics are sorted by latitude, month and local time. Power spectral density (PSD) exhibits a universal exp(-km) structure in the 10-20km wavelength regime, with K lying between 2 and 3. The observed PSD for wavelengths between 10 and 20 km is a factor of 3 higher than a null spectrum constructed from HRDI reported error bars multiplied by randomly varying numbers between -1 and +1. Stokes parameters were consolidated by month into Northern and Southern hemisphere middle and high latitudes belts (40-72 degrees), tidal belts (32-16 degrees) and a tropical belt (8S-8N). Vertical waves between 10 and 15 km in wavelength are about 10-15% polarized everywhere. The inferred propagation direction in the middle and high latitude Southern hemisphere is predominantly meridional during solstice, and significantly more zonal during equinoxes. In the tropical belt, the wave orientations are nearly North-South during solstices, with a slightly higher east-west component during equinox. In the tidal belts where the background wind includes a strong meridional tidal wind, the preferred wave orientation has a significant zonal component during equinox. These findings are consistent with the interpretation of wave filtering by the background wind.

Lieberman, Ruth S.↗

Scaling Relation for Occulter Manufacturing Errors

An external occulter is a spacecraft own along the line-of-sight of a space telescope to suppress starlight and enable high-contrast direct imaging of exoplanets. The shape of an external occulter must be specially designed to optimally suppress starlight and deviations from the ideal shape due to manufacturing errors can result loss of suppression in the shadow. Due to the long separation distances and large dimensions involved for a space occulter, laboratory testing is conducted with scaled versions of occulters etched on silicon wafers. Using numerical simulations for a flight Fresnel occulter design, we show how the suppression performance of an occulter mask scales with the available propagation distance for expected random manufacturing defects along the edge of the occulter petal. We derive an analytical model for predicting performance due to such manufacturing defects across the petal edges of an occulter mask and compare this with the numerical simulations. We discuss the scaling of an extended occulter test-bed.

Sirbu, Dan↗

Data Fusion of In-Flight Aerothermodynamic Heating Measurements Using Kalman Filtering

On February 18th, 2021, the Mars 2020 entry system successfully delivered the Perserverance rover to the surface of Mars at Jezero Crater. The entry capsule carried a set of instrumentation installed on the heat shield and backshell, named the Mars 2020 Entry, Descent, and Landing Instrumentation. The instruments include pressure transducers, thermocouples, heat flux gauges, and radiometers to measure the aerodynamic and aerothermodynamic performance of the entry vehicle. Three of these sensors, a thermocouple plug, heat flux gauge, and a radiometer, are located in close proximity on the backshell. Each sensor is exposed to roughly the same environment, but measured these environments in different ways, each with its own set of modeling and measurement error complications. This paper develops a method for blending each of these measurements together in a single algorithm to produce estimates of the aerothermodynamic environments at that location on the backshell. The approach makes use of the Kalman filter methodology for solving state estimation problems. The filter has a predictor/corrector structure in which simplified process models are used to propagate the aeroheating states forward through time, which are then updated based on the measurement data. The method has been tested on simulated datasets.

Christopher D. Karlgaard↗

Uncertainty Analysis of Instrument Calibration and Application

Experimental aerodynamic researchers require estimated precision and bias uncertainties of measured physical quantities, typically at 95 percent confidence levels. Uncertainties of final computed aerodynamic parameters are obtained by propagation of individual measurement uncertainties through the defining functional expressions. In this paper, rigorous mathematical techniques are extended to determine precision and bias uncertainties of any instrument-sensor system. Through this analysis, instrument uncertainties determined through calibration are now expressed as functions of the corresponding measurement for linear and nonlinear univariate and multivariate processes. Treatment of correlated measurement precision error is developed. During laboratory calibration, calibration standard uncertainties are assumed to be an order of magnitude less than those of the instrument being calibrated. Often calibration standards do not satisfy this assumption. This paper applies rigorous statistical methods for inclusion of calibration standard uncertainty and covariance due to the order of their application. The effects of mathematical modeling error on calibration bias uncertainty are quantified. The effects of experimental design on uncertainty are analyzed. The importance of replication is emphasized, techniques for estimation of both bias and precision uncertainties using replication are developed. Statistical tests for stationarity of calibration parameters over time are obtained.

Tripp, John S.↗

Reply to Comment on “Coronal Mass Ejections, Interplanetary Ejecta and Geomagnetic Storms” by Gopalswamy Et Al.

The comment of Gopalswamy et al. (thereafter GMY) relates to a letter discussing coronal mass ejections (CMEs), interplanetary ejecta and geomagnetic storms. GMY contend that Cane et al. incorrectly identified ejecta (interplanetary CMEs) and hypothesize that this is because Cane et al. fail to understand how to separate ejecta from "shock sheaths" when interpreting solar wind and energetic particle data sets. They (GMY) are wrong be cause the relevant section of the paper was concerned with the propagation time to 1 AU of any potentially geoeffective structures caused by CMEs, i.e. upstream compression regions with or without shocks, or ejecta. In other words, the travel times used by Cane et al. were purposefully and deliberately distinct from ejecta travel times (except for those slow ejecta, approx. 30% of their events, which generated no upstream features), and no error in identification was involved. The confusion of GMY stems from the description did not characterize the observations sufficiently clearly.

Cane, H. V.↗

Short-Block Protograph-Based LDPC Codes

Short-block low-density parity-check (LDPC) codes of a special type are intended to be especially well suited for potential applications that include transmission of command and control data, cellular telephony, data communications in wireless local area networks, and satellite data communications. [In general, LDPC codes belong to a class of error-correcting codes suitable for use in a variety of wireless data-communication systems that include noisy channels.] The codes of the present special type exhibit low error floors, low bit and frame error rates, and low latency (in comparison with related prior codes). These codes also achieve low maximum rate of undetected errors over all signal-to-noise ratios, without requiring the use of cyclic redundancy checks, which would significantly increase the overhead for short blocks. These codes have protograph representations; this is advantageous in that, for reasons that exceed the scope of this article, the applicability of protograph representations makes it possible to design highspeed iterative decoders that utilize belief- propagation algorithms.

Divsalar, Dariush↗

Assimilation of Satellite-Based Snow Cover and Freeze/Thaw Observations Over the High Mountain Asia

Toward qualifying hydrologic changes in the High Mountain Asia (HMA) region, this study explores the use of a hyper-resolution (1 km) land data assimilation (DA)framework developed within the NASA Land Information System using the Noah Multi-parameterization Land Surface Model (Noah-MP) forced by the meteorological boundary conditions from Modern-Era Retrospective analysis for Research and Applications, Version 2 data. Two different sets of DA experiments are conducted:(1) the assimilation of a satellite-derived snow cover map (MOD10A1) and (2) the assimilation of the NASA MEaSUREs landscape freeze/thaw product from 2007 to 2008. The performance of the snow cover assimilation is evaluated via comparisons with available remote sensing-based snow water equivalent product and ground-based snow depth measurements. For example, in the comparison against ground-based snow depth measurements, the majority of the stations (13 of 14) show slightly improved goodness-of-fit statistics as a result of the snow DA, but only four are statistically significant. In addition, comparisons to the satellite-based land surface temperature products (MOD11A1 and MYD11A1) show that freeze/thaw DA yields improvements (at certain grid cells) of up to 0.58 K in the root-mean-square error (RMSE) and 0.77K in the absolute bias (relative to model-only simulations). In the comparison against three ground-based soil temperature measurements along the Himalayas, the bias and the RMSE in the 0-10 cm soil temperature are reduced (on average) by 10 and 7%,respectively. The improvements in the top layer of soil estimates also propagate through the deeper soil layers, where the bias and the RMSE in the 10-40 cm soil temperature are reduced (on average) by 9 and 6%, respectively. However, no statistically significant skill differences are observed for the freeze/thaw DA system in the comparisons against ground-based surface temperature measurements at mid-to-low altitude. Therefore, the two proposed DA schemes show the potential of improving the predictability of snow mass, surface temperature, and soil temperature states across HMA, but more ground-based measurements are still required, especially at high-altitudes, in order to document a more statistically significant improvement as a result of the two DA schemes.

High Mountain Asia↗

Dependence of the interstellar propagation of ultraheavy cosmic-ray nuclei on various parameters

A stepwise program of interstellar-cosmic-ray propagation suitable for implementation on a small work-station has been used to study the sensitivity of propagation models of ultraheavy cosmic-ray nuclei (UH nuclei) on the input parameters. In this program an initial set of relevant parameters describing the astrophysics of the propagations process and the nuclear physics representing the interactions was assumed. These parameters were then varied within the estimated uncertainties in order to evaluate the sensitivity of the predictions on the assumed values. It is shown that current uncertainties in the estimates of the UH nuclei nuclear cross sections for the production of lighter fragments present the most significant sources of error and are sufficient to preclude reaching detailed conclusions about the source composition. Other uncertainties, such as those in the escape probabilities, or the mean interstellar density, while still significant, are less serious.

Clinton, Robert R.↗

Variable Step Integration Coupled with the Method of Characteristics Solution for Water-Hammer Analysis, A Case Study

One-dimensional water-hammer modeling involves the solution of two coupled non-linear hyperbolic partial differential equations (PDEs). These equations result from applying the principles of conservation of mass and momentum to flow through a pipe, and usually the assumption that the speed at which pressure waves propagate through the pipe is constant. In order to solve these equations for the interested quantities (i.e. pressures and flow rates), they must first be converted to a system of ordinary differential equations (ODEs) by either approximating the spatial derivative terms with numerical techniques or using the Method of Characteristics (MOC). The MOC approach is ideal in that no numerical approximation errors are introduced in converting the original system of PDEs into an equivalent system of ODEs. Unfortunately this resulting system of ODEs is bound by a time step constraint so that when integrating the equations the solution can only be obtained at fixed time intervals. If the fluid system to be modeled also contains dynamic components (i.e. components that are best modeled by a system of ODEs), it may be necessary to take extremely small time steps during certain points of the model simulation in order to achieve stability and/or accuracy in the solution. Coupled together, the fixed time step constraint invoked by the MOC, and the occasional need for extremely small time steps in order to obtain stability and/or accuracy, can greatly increase simulation run times. As one solution to this problem, a method for combining variable step integration (VSI) algorithms with the MOC was developed for modeling water-hammer in systems with highly dynamic components. A case study is presented in which reverse flow through a dual-flapper check valve introduces a water-hammer event. The predicted pressure responses upstream of the check-valve are compared with test data.

Turpin, Jason B.↗

Advantages of estimating rate corrections during dynamic propagation of spacecraft rates: Applications to real-time attitude determination of SAMPEX

This paper describes real-time attitude determination results for the Solar, Anomalous, and Magnetospheric Particle Explorer (SAMPEX), a gyroless spacecraft, using a Kalman filter/Euler equation approach denoted the real-time sequential filter (RTSF). The RTSF is an extended Kalman filter whose state vector includes the attitude quaternion and corrections to the rates, which are modeled as Markov processes with small time constants. The rate corrections impart a significant robustness to the RTSF against errors in modeling the environmental and control torques, as well as errors in the initial attitude and rates, while maintaining a small state vector. SAMPLEX flight data from various mission phases are used to demonstrate the robustness of the RTSF against a priori attitude and rate errors of up to 90 deg and 0.5 deg/sec, respectively, as well as a sensitivity of 0.0003 deg/sec in estimating rate corrections in torque computations. In contrast, it is shown that the RTSF attitude estimates without the rate corrections can degrade rapidly. RTSF advantages over single-frame attitude determination algorithms are also demonstrated through (1) substantial improvements in attitude solutions during sun-magnetic field coalignment and (2) magnetic-field-only attitude and rate estimation during the spacecraft's sun-acquisition mode. A robust magnetometer-only attitude-and-rate determination method is also developed to provide for the contingency when both sun data as well as a priori knowledge of the spacecraft state are unavailable. This method includes a deterministic algorithm used to initialize the RTSF with coarse estimates of the spacecraft attitude and rates. The combined algorithm has been found effective, yielding accuracies of 1.5 deg in attitude and 0.01 deg/sec in the rates and convergence times as little as 400 sec.

Challa, M. S.↗

Some aspects of high-order numerical solutions of the linear convection equation with forced boundary conditions

A six-stage low-storage Runge-Kutta time-marching method is presented and shown to be an efficient method for use with high-accuracy spatial difference operators for wave propagation problems. The accuracy of the method for inhomogeneous ordinary differential equations is demonstrated through numerical solutions of the linear convection equation with forced boundary conditions. Numerical experiments are presented simulating a sine wave and a Gaussian pulse propagating into and through the domain. For practical levels of mesh refinement corresponding to roughly ten points per wavelength, the six-stage Runge-Kutta method is more accurate than the popular fourth-order Runge-Kutta method. Further numerical experiments are presented which show that the numerical boundary scheme at an inflow boundary can be a significant source of error when high-accuracy spatial discretizations are used.

Zingg, D. W.↗

Demand Capacity Balancing at Vertiports for Initial Strategic Conflict Management of Urban Air Mobility Operations

Urban Air Mobility (UAM) is a new transportation concept that enables highly automated, cooperative, passenger or cargo-carrying air transportation services in and around urban areas. To achieve the high level of operational density and complexity desired by the UAM community, an airspace system that allows UAM operators to readily access and operate safely and efficiently in the airspace is needed. This airspace system will require air traffic management designed to reduce the risk of conflicts and loss of separation between UAM flights. In general, strategic conflict management is considered as the first layer of conflict management for safe flight operations to condition the traffic to reduce the need for airborne separation provision, the second layer of conflict management. Demand Capacity Balancing (DCB) is one of the concept components to achieve strategic conflict management. DCB strategically evaluates traffic demand and resource capacities to allow UAM operators to determine when, where and how they operate, while mitigating conflicting needs for airspace and vertiport capacity. DCB can be applied whenever UAM demand exceeds the capacity in airspace or at vertiports. As the UAM ecosystem evolves with advanced technologies and matured operational procedures, more complicated conflict management will likely be needed. In the current UAM ‘Concept of Operation (ConOps) 1.0’ operational stage defined by FAA, however, it will be meaningful to explore the demand capacity balancing at vertiports only, as an initial strategic conflict management approach for UAM operations because vertiport capacity seems to be a bottleneck of UAM traffic. For this research, we developed a demand-capacity imbalance detection and resolution service for UAM. This DCB service identifies the demand from operators and compares the demand to a given capacity at the shared resources (i.e., vertiports) over the upcoming time horizon which is divided into time bins having a constant interval. When a new flight plan is submitted, the algorithm embedded in the DCB service checks the available time bins based on the desired departure time and estimated arrival time at origin and destination vertiports, respectively. If the time bins for the originally desired times are already occupied by other flights (i.e., demand is at or above capacity), the algorithm finds the next available time bins for takeoff and landing and shifts the conflicting departure time to the earliest time that satisfies the capacity constraints at both origin and destination vertiports. The details of the algorithm will be described in the final manuscript. Figure 1 shows that the proposed DCB algorithm works well for a sample traffic scenario. In this example, a total of 144 flights, split between two operators, are planned over 2 hours, traveling 10 routes between five vertiports. In the heatmaps, the horizontal axis shows 12 time bins where each bin represents a 12-minute interval, and the vertical axis shows five vertiports. The number in each cell shows the number of operations, counting both departures and arrivals, at a specific vertiport in each time bin. For the given capacity of 2 operations/vertiport/bin, Figure 1 shows that the original demand sometimes exceeds the capacity, but the modified demand is reduced to the given capacity after resolving demand-capacity imbalances. When UAM flights are operated, it is expected that many practical issues would arise in the federated system architecture with multiple operators. UAM operators may experience a time synchronization issue due to communication delay between operator and vehicle. UAM vehicles would fly at different flight speeds, depending on vehicle models. Actual departure and arrival times can have large variations, compared to the schedule. The lead time from flight plan submission to desired departure time can vary by service type (e.g., regular shuttle service vs. on-demand service). Using the proposed DCB algorithm, we also investigated how the actual flight schedule and DCB performance are affected by these uncertainties such as unsynchronized times between operators, flight speed differences, lead time differences, and departure time errors. The final manuscript will include the background of this research work, the description of the DCB algorithm and its use cases with traffic scenarios. It will also provide the analytical results about the impact of various uncertainties that can occur in actual UAM operations on the DCB at vertiports, in terms of demand distribution changes, number of simultaneous operations, and delay propagation.

Urban Air Mobility↗

Evaluation of fatigue-crack growth rates by polynomial curve fitting

Fundamental characterization of the constant-amplitude fatigue crack propagation is achieved by an analysis of the rate of change of crack length with change in number of applied loading cycles, defining the rate values such that they are consistent with the basic assumption of smoothness and continuity in the fatigue crack growth process. The technique used to satisfy the analytical conditions and minimize the effects of local material anomalies and experimental errors is that of fitting a smooth curve to the entire set of basic data by least square regression. This yields a well-behaved function relating the number of cycles to the crack length. By taking the first derivative of the function, the crack growth rate is obtained for each point. The class of curve fitting functions used in the analysis is the polynomial of degree n.

Davies, K. B.↗

Quantification of Error Associated with Using Misaligned Meshes in Continuum Damage Mechanics Material Models for Matrix Crack Growth Predictions in Composites

The ability of a material model to capture in-plane matrix mode I and mode II crack growth is an essential component for modeling ply level damage evolution in composite structures. Previous studies using a continuum damage mechanics (CDM) approach have shown success in satisfying benchmark solutions for mode I and II crack growth. However, success was shown using a fiber-aligned meshing strategy, which encourages matrix cracks to propagate in a single band of elements, along the fiber direction. Generating a fiber-aligned mesh becomes a highly involved process for laminates including off-axis (non 0° or 90°) plies. The objective of this study is to quantify the effect of non-fiber aligned mesh discretization on predictions of in-plane matrix crack propagation. The approach taken incrementally varies the mesh orientation angle relative to the fiber orientation; more specifically, misaligned meshes are used to quantify the effect of element angle orientation relative to the initial crack orientation on the energy released during matrix crack propagation simulations using a CDM method. CDM solutions obtained with the misaligned meshes are evaluated against known benchmarks for mode I and II matrix crack growth. The CDM solutions reveal a near-polynomial trend of increased predicted failure stress with increased mesh misalignment angle; hence implying a potential relationship between element orientation angle and apparent fracture toughness.

Justusson, Brian↗

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.↗

Influence of Rayleigh-Doppler broadening on the selection of H2O dial system parameters

Computer simulations have enabled the performance of a H2O Differential Absorption Lidar (DIAL) system to be studied by spectrally analyzing the forward propagating and backscattered laser energy. The simulations were done for a high altitude (21 km) DIAL system operating in a nadir-viewing mode. The influence of Rayleigh Doppler broadening on DIAL measurement accuracies were evaluated and show that the Rayleigh broadening influence, which can be corrected to first order in regions free of large aerosol gradients, reduces the sensitivity of DIAL H2O measurement errors in the upper tropospheric region. The ability to correct the Rayleigh broadening and the selection of H2O DIAL parameters when all the systematic effects are combined, were discussed.

Ismail, S.↗

Land-mobile-satellite fade measurements in Australia

Attenuation measurements were implemented at L-band (1.5 GHz) in southeastern Australia during an 11-day period in October 1988 as part of a continuing examination of the propagation effects due to roadside trees and terrain for mobile-satellite service. Beacon transmissions from the geostationary ETS-V and IPORS satellites were observed. The Australian campaign expanded to another continent our Mobile Satellite Service data base of measurements executed in the eastern and southwestern United States regions. An empirical fade distribution model based on U.S. data predicted the Australian results with errors generally less than 1 dB in the 1-20 percent probability region. Directive antennas are shown to suffer deeper fades under severe shadowing conditions (3 dB excess at 4 percent), the equal-probability isolation between co- and cross-polarized transmissions deteriorated to 10 dB at the 5 dB fade level, and antenna diversity reception may reduce unavailability of the system by a factor of 2-8.

Vogel, Wolfhard J.↗

The tropospheric calibration model for Mariner Mars 1971

The effects of the troposphere on arrival time, range, and range rate measurement of radio signals propagating from spacecraft are examined. A tropospheric calibration model which takes into account seasonal variations, various elevations, and apparent location changes was developed. A comparison was made of the effects before and after calibration. Without calibration, a 16-m change in the distance off the spin axis, is induced for a 10-deg minimum elevation symmetrical pass. Even after calibration with the old model, the error in spin axis is still greater than 1 m. Calibration results and estimated uncertainties of the model are given in tabular form. Test results also show that the model met the required accuracy for the Mariner Mars 1971 mission.

C. C. Chao↗