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 721 records · Page 40

A numerical study of viscous vortex rings using a spectral method

Viscous, axisymmetric vortex rings are investigated numerically by solving the incompressible Navier-Stokes equations using a spectral method designed for this type of flow. The results presented are axisymmetric, but the method is developed to be naturally extended to three dimensions. The spectral method relies on divergence-free basis functions. The basis functions are formed in spherical coordinates using Vector Spherical Harmonics in the angular directions, and Jacobi polynomials together with a mapping in the radial direction. Simulations are performed of a single ring over a wide range of Reynolds numbers (Re approximately equal gamma/nu), 0.001 less than or equal to 1000, and of two interacting rings. At large times, regardless of the early history of the vortex ring, it is observed that the flow approaches a Stokes solution that depends only on the total hydrodynamic impulse, which is conserved for all time. At small times, from an infinitely thin ring, the propagation speeds of vortex rings of varying Re are computed and comparisons are made with the asymptotic theory by Saffman. The results are in agreement with the theory; furthermore, the error is found to be smaller than Saffman's own estimate by a factor square root ((nu x t)/R squared) (at least for Re=0). The error also decreases with increasing Re at fixed core-to-ring radius ratio, and appears to be independent of Re as Re approaches infinity). Following a single ring, with Re=500, the vorticity contours indicate shedding of vorticity into the wake and a settling of an initially circular core to a more elliptical shape, similar to Norbury's steady inviscid vortices. Finally, we consider the case of leapfrogging vortex rings with Re=1000. The results show severe straining of the inner vortex core in the first pass and merging of the two cores during the second pass.

Stanaway, S. K.↗

Isotopic Dependence of GCR Fluence behind Shielding

In this paper we consider the effects of the isotopic composition of the primary galactic cosmic rays (GCR), nuclear fragmentation cross-sections, and isotopic-grid on the solution to transport models used for shielding studies. Satellite measurements are used to describe the isotopic composition of the GCR. For the nuclear interaction data-base and transport solution, we use the quantum multiple-scattering theory of nuclear fragmentation (QMSFRG) and high-charge and energy (HZETRN) transport code, respectively. The QMSFRG model is shown to accurately describe existing fragmentation data including proper description of the odd-even effects as function of the iso-spin dependence on the projectile nucleus. The principle finding of this study is that large errors (+/-100%) will occur in the mass-fluence spectra when comparing transport models that use a complete isotopic-grid (approx.170 ions) to ones that use a reduced isotopic-grid, for example the 59 ion-grid used in the HZETRN code in the past, however less significant errors (<+/-20%) occur in the elemental-fluence spectra. Because a complete isotopic-grid is readily handled on small computer workstations and is needed for several applications studying GCR propagation and scattering, it is recommended that they be used for future GCR studies.

Cucinotta, Francis A.↗

Numerical Study of Buoyancy and Different Diffusion Effects on the Structure and Dynamics of Triple Flames

Numerical simulations of 2-D triple flames under gravity force have been implemented to identify the effects of gravity on triple flame structure and propagation properties and to understand the mechanisms of instabilities resulting from both heat release and buoyancy effects. A wide range of gravity conditions, heat release, and mixing widths for a scalar mixing layer are computed for downward-propagating (in the same direction with the gravity vector) and upward-propagating (in the opposite direction of the gravity vector) triple flames. Results of numerical simulations show that gravity strongly affects the triple flame speed through its contribution to the overall flow field. A simple analytical model for the triple flame speed, which accounts for both buoyancy and heat release, is developed. Comparisons of the proposed model with the numerical results for a wide range of gravity, heat release and mixing width conditions, yield very good agreement. The analysis shows that under neutral diffusion, downward propagation reduces the triple flame speed, while upward propagation enhances it. For the former condition, a critical Froude number may be evaluated, which corresponds to a vanishing triple flame speed. Downward-propagating triple flames at relatively strong gravity effects have exhibited instabilities. These instabilities are generated without any artificial forcing of the flow. Instead disturbances are initiated by minute round-off errors in the numerical simulations, and subsequently amplified by instabilities. A linear stability analysis on mean profiles of stable triple flame configurations have been performed to identify the most amplified frequency in spatially developed flows. The eigenfunction equations obtained from the linearized disturbance equations are solved using the shooting method. The linear stability analysis yields reasonably good agreements with the observed frequencies of the unstable triple flames. The frequencies and amplitudes of disturbances increase with the magnitude of the gravity vector. Moreover, disturbances appear to be most amplified just downstream of the premixed branches. The effects of mixing width and differential diffusion are investigated and their roles on the flame stability are studied.

Chen, Jyh-Yuan↗

Physics of the Isotopic Dependence of Galactic Cosmic Ray Fluence Behind Shielding

For over 25 years, NASA has supported the development of space radiation transport models for shielding applications. The NASA space radiation transport model now predicts dose and dose equivalent in Earth and Mars orbit to an accuracy of plus or minus 20%. However, because larger errors may occur in particle fluence predictions, there is interest in further assessments and improvements in NASA's space radiation transport model. In this paper, we consider the effects of the isotopic composition of the primary galactic cosmic rays (GCR) and the isotopic dependence of nuclear fragmentation cross-sections on the solution to transport models used for shielding studies. Satellite measurements are used to describe the isotopic composition of the GCR. Using NASA's quantum multiple-scattering theory of nuclear fragmentation (QMSFRG) and high-charge and energy (HZETRN) transport code, we study the effect of the isotopic dependence of the primary GCR composition and secondary nuclei on shielding calculations. The QMSFRG is shown to accurately describe the iso-spin dependence of nuclear fragmentation. The principal finding of this study is that large errors (plus or minus 100%) will occur in the mass-fluence spectra when comparing transport models that use a complete isotope grid (approximately 170 ions) to ones that use a reduced isotope grid, for example the 59 ion-grid used in the HZETRN code in the past, however less significant errors (less than 20%) occur in the elemental-fluence spectra. Because a complete isotope grid is readily handled on small computer workstations and is needed for several applications studying GCR propagation and scattering, it is recommended that they be used for future GCR studies.

Cucinotta, Francis A.↗

Space-Time Reduced-Order Modeling for Uncertainty Quantification

This work focuses on the space-time reduced-order modeling (ROM) method for solving large-scale uncertainty quantification (UQ) problems with multiple random coefficients. In contrast with the traditional space ROM approach, which performs dimension reduction in the spatial dimension, the space-time ROM approach performs dimension reduction on both the spatial and temporal domains, and thus enables accurate approximate solutions at a low cost. We incorporate the space-time ROM strategy with various classical stochastic UQ propagation methods such as stochastic Galerkin and Monte Carlo. Numerical results demonstrate that our methodology has significant computational advantages compared to state-of-the-art ROM approaches. By testing the approximation errors, we show that there is no obvious loss of simulation accuracy for space-time ROM given its high computational efficiency.

97 MATHEMATICS AND COMPUTING↗

Balloon atmospheric propagation experiment measurements

High altitude balloon measurements on laser beam fading during propagation through turbulent atmosphere show that a correlation between fading strength and stellar scintillation magnitudes exists. Graphs for stellar scintillation as a function of receiver aperture are used to predict fading bit error rates for neodymium-yag laser communication system.

Minott, P. O.↗

Evaluation of a scale-model experiment to investigate long-range acoustic propagation

Tests were conducted to evaluate the feasibility of using a scale-model experiment situated in an anechoic facility to investigate long-range sound propagation over ground terrain. For a nominal scale factor of 100:1, attenuations along a linear array of six microphones colinear with a continuous-wave type of sound source were measured over a wavelength range from 10 to 160 for a nominal test frequency of 10 kHz. Most tests were made for a hard model surface (plywood), but limited tests were also made for a soft model surface (plywood with felt). For grazing-incidence propagation over the hard surface, measured and predicted attenuation trends were consistent for microphone locations out to between 40 and 80 wavelengths. Beyond 80 wavelengths, significant variability was observed that was caused by disturbances in the propagation medium. Also, there was evidence of extraneous propagation-path contributions to data irregularities at more remote microphones. Sensitivity studies for the hard-surface and microphone indicated a 2.5 dB change in the relative excess attenuation for a systematic error in source and microphone elevations on the order of 1 mm. For the soft-surface model, no comparable sensitivity was found.

Parrott, Tony L.↗

The use of holographic interferometry for measurements of temperature in a rectangular heat pipe

Holographic interferometry is a nonintrusive method and as such possesses considerable advantages such as not disturbing the velocity and temperature field by creating obstacles which would alter the flow field. These optical methods have disadvantages as well. Holography, as one of the interferometry methods, retains the accuracy of older methods, and at the same time eliminates the system error of participating components. The holographic interferometry consists of comparing the objective beam with the reference beam and observing the difference in lengths of optical paths, which can be observed during the propagation of the light through a medium with locally varying refractive index. Thus, change in refractive index can be observed as a family of nonintersecting surfaces in space (wave fronts). The object of the investigation was a rectangular heat pipe. The goal was to measure temperatures in the heat pipe, which yields data for computer code or model assessment. The results were obtained by calculating the temperatures by means of finite fringes.

Marn, Jure↗

Assessing the Impact of SMAP Soil Moisture Data Assimilation on the Simulation and Prediction of Tropical Cyclone Idai

The role of soil moisture in the evolution of tropical cyclones (TCs) approaching land has long been recognized. Dry land surface conditions can lead to faster dissipation of a TC over land (often associated with changes in precipitation structure), whereas very wet conditions can help sustain or even re-intensify a TC. The ability to forecast post-landfall TC evolution and mitigate the associated socio-economic impact thus hinges on accurate knowledge of land surface conditions prior to landfall. The NASA Soil Moisture Active Passive (SMAP) mission provides accurate observations of soil moisture globally and at high revisit times of 2-3 days. It has been shown that the assimilation of SMAPbrightness temperatures (Tb) significantly improves modeled land surface states and thus has the potential to constrain land surface initial conditions in TC forecasts. In this presentation, we investigate this potential through an extensive set of Observing System Experiments that systematically assess the impact of assimilating SMAP Tbs on TC forecast skill in the Goddard Earth Observing System. Focusing on the case of TC Idai, we show that the assimilation of SMAP generally leads to a drier soil moisture analysis for the land surface underneath much of the storm’s circulation. An exception is the direct TC track, where SMAP assimilation increases soil moisture. These changes are reflected in the surface fluxes, which propagate the land surface state changes to the atmosphere. This results a TC that is overall weaker with lower winds, but with a better-defined eye and a slightly smaller track error than a simulation without SMAP DA. The changes in soil moisture resulting from the assimilation of SMAP also impact the total precipitation amounts as well as the precipitation structure associated with TC Idai.

SMAP↗

Evaluation of Two Computational Techniques of Calculating Multipath Using Global Positioning System Carrier Phase Measurements

Two computational techniques are used to calculate differential phase errors on Global Positioning System (GPS) carrier war phase measurements due to certain multipath-producing objects. The two computational techniques are a rigorous computati electromagnetics technique called Geometric Theory of Diffraction (GTD) and the other is a simple ray tracing method. The GTD technique has been used successfully to predict microwave propagation characteristics by taking into account the dominant multipath components due to reflections and diffractions from scattering structures. The ray tracing technique only solves for reflected signals. The results from the two techniques are compared to GPS differential carrier phase ns taken on the ground using a GPS receiver in the presence of typical International Space Station (ISS) interference structures. The calculations produced using the GTD code compared to the measured results better than the ray tracing technique. The agreement was good, demonstrating that the phase errors due to multipath can be modeled and characterized using the GTD technique and characterized to a lesser fidelity using the DECAT technique. However, some discrepancies were observed. Most of the discrepancies occurred at lower devations and were either due to phase center deviations of the antenna, the background multipath environment, or the receiver itself. Selected measured and predicted differential carrier phase error results are presented and compared. Results indicate that reflections and diffractions caused by the multipath producers, located near the GPS antennas, can produce phase shifts of greater than 10 mm, and as high as 95 mm. It should be noted tl the field test configuration was meant to simulate typical ISS structures, but the two environments are not identical. The GZ and DECAT techniques have been used to calculate phase errors due to multipath o the ISS configuration to quantify the expected attitude determination errors.

Gomez, Susan F.↗

Exploring experimental conditions to reduce uncertainties in the optical potential

Uncertainty quantification for nuclear theories has gained a more prominent role in the field, with more groups attempting to understand the uncertainties in their calculations. However, recent studies have shown that the uncertainties on the optical potentials are too large for the theory to be useful. The goal of this work is to explore possible experimental conditions that may reduce the uncertainties on elastic scattering and single-nucleon transfer cross sections that come from the fitting of the optical model parameters to experimental data. Using Bayesian methods, we explore the effect of the uncertainties of optical model parameters on the angular grid of the differential cross section, including cross-section data at nearby energies, and changes in the experimental error bars. We also study the effect on the resulting uncertainty when other observables are included in the fitting procedure, particularly the total (reaction) cross sections. We study proton and neutron elastic scattering on 48 Ca and 208 Pb. We explore the parameter space with the Markov-chain Monte Carlo method, produce posterior distributions for the optical model parameters, and construct the corresponding 95% confidence intervals on the elastic-scattering cross sections. We also propagate the uncertainties on the optical potentials to the 48 Ca(d,p) 49 Ca(g.s.) and 208 Pb(d,p) 209 Pb(g.s.) cross sections. We find little sensitivity to the angular grid and an improvement of up to a factor of 2 on the uncertainties by including data at a nearby energy. Although reducing the error bars on the data does reduce the uncertainty, the gain is often considerably smaller than one would naively expect. We also find that the inclusion of total reaction cross section can improve the uncertainty although the magnitude of the effect depends strongly on the cases considered.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Development of improved gravity field estimates using existing altimeter data

The formulation of a Fourier series and geographic representation expressions for orbit error were reported. Addition activities include the study of undulation errors in the Lagrangian formulation, and the development and testing of all the necessary software for the above expression. Additionally, the matrix formulations that led to the formation of observed equations and covariance propagation were developed.

Source record↗

Hubble Space Telescope characterized by using phase-retrieval algorithms

We describe several results characterizing the Hubble Space Telescope from measured point spread functions by using phase-retrieval algorithms. The Cramer-Rao lower bounds show that point spread functions taken well out of focus result in smaller errors when aberrations are estimated and that, for those images, photon noise is not a limiting factor. Reconstruction experiments with both simulated and real data show that the calculation of wave-front propagation by the retrieval algorithms must be performed with a multiple-plane propagation rather than a simple fast Fourier transform to ensure the high accuracy required. Pupil reconstruction was performed and indicates a misalignment of the optical axis of a camera relay telescope relative to the main telescope. After we accounted for measured spherical aberration in the relay telescope, our estimate of the conic constant of the primary mirror of the HST was - 1.0144.

Fienup, J. R.↗

Neural Network Prediction of Failure of Damaged Composite Pressure Vessels from Strain Field Data Acquired by a Computer Vision Method

This effort used a new and novel method of acquiring strains called Sub-pixel Digital Video Image Correlation (SDVIC) on impact damaged Kevlar/epoxy filament wound pressure vessels during a proof test. To predict the burst pressure, the hoop strain field distribution around the impact location from three vessels was used to train a neural network. The network was then tested on additional pressure vessels. Several variations on the network were tried. The best results were obtained using a single hidden layer. SDVIC is a fill-field non-contact computer vision technique which provides in-plane deformation and strain data over a load differential. This method was used to determine hoop and axial displacements, hoop and axial linear strains, the in-plane shear strains and rotations in the regions surrounding impact sites in filament wound pressure vessels (FWPV) during proof loading by internal pressurization. The relationship between these deformation measurement values and the remaining life of the pressure vessels, however, requires a complex theoretical model or numerical simulation. Both of these techniques are time consuming and complicated. Previous results using neural network methods had been successful in predicting the burst pressure for graphite/epoxy pressure vessels based upon acoustic emission (AE) measurements in similar tests. The neural network associates the character of the AE amplitude distribution, which depends upon the extent of impact damage, with the burst pressure. Similarly, higher amounts of impact damage are theorized to cause a higher amount of strain concentration in the damage effected zone at a given pressure and result in lower burst pressures. This relationship suggests that a neural network might be able to find an empirical relationship between the SDVIC strain field data and the burst pressure, analogous to the AE method, with greater speed and simplicity than theoretical or finite element modeling. The process of testing SDVIC neural network analysis and some encouraging preliminary results are presented in this paper. Details are given concerning the processing of SDVIC output data such that it may be used as back propagation neural network (BPNN) input data. The software written to perform this processing and the BPNN algorithm are also discussed. It will be shown that, with limited training, test results indicate an average error in burst pressure prediction of approximately six percent,

Russell, Samuel S.↗

The Radial Speed-Expansion Speed Relation for Earth-Directed CMEs

Earth-directed coronal mass ejections (CMEs) are the main drivers of major geomagnetic storms. Therefore, a good estimate of the disturbance arrival time at Earth is required for space weather predictions. The STEREO and SOHO spacecraft were viewing the Sun in near quadrature during January 2010 to September 2012, providing a unique opportunity to study the radial speed (V (sub rad)) to expansion speed(V (sub exp)) relationship of Earth-directed CMEs. This relationship is useful in estimating the V (sub rad) of Earth-directed CMEs, when they are observed from Earth view only. We selected 19 Earth-directed CMEs observed by the Large Angle and Spectrometric Coronagraph (LASCO)/C3 coronagraph on SOHO and the Sun Earth Connection Coronal and Heliospheric Investigation (SECCHI)/COR2 coronagraph on STEREO during January 2010 to September 2012. We found that of the three tested geometric CME models the full ice-cream cone model of the CME describes best the V (sub rad) to V (sub exp) relationship, as suggested by earlier investigations. We also tested the prediction accuracy of the empirical shock arrival (ESA) model proposed by Gopalswamy et al.(2005a), while estimating the CME propagation speeds from the CME expansion speeds. If we use STEREO observations to estimate the CME width required to calculate the V (sub rad) from the V (sub exp) measurements, the mean absolute error (MAE) of the shock arrival times of the ESA model is 8.4 hours. If the LASCO measurements are used to estimate the CME width, the MAE still remains below 17 hours. Therefore, by using the simple V (sub rad) to V (sub exp) relationship to estimate the V (sub rad) of the Earth-directed CMEs, the ESA model is able to predict the shock arrival times with accuracy comparable to most other more complex models.

Makela, P.↗

WFIRST coronagraph optical modeling

End-to-end numerical optical modeling of the WFIRST coronagraph incorporating wavefront sensing and control is used to determine the performance of the coronagraph with realistic errors, including pointing jitter and polarization. We present the performance estimates of the current flight designs as predicted by modeling. We also describe the release of a new version of the PROPER optical propagation library, our primary modeling tool, which is now available for Python and Matlab in addition to IDL.

Zhou, Hanying↗

Calibration of Predictor Models Using Multiple Validation Experiments

This paper presents a framework for calibrating computational models using data from several and possibly dissimilar validation experiments. The offset between model predictions and observations, which might be caused by measurement noise, model-form uncertainty, and numerical error, drives the process by which uncertainty in the models parameters is characterized. The resulting description of uncertainty along with the computational model constitute a predictor model. Two types of predictor models are studied: Interval Predictor Models (IPMs) and Random Predictor Models (RPMs). IPMs use sets to characterize uncertainty, whereas RPMs use random vectors. The propagation of a set through a model makes the response an interval valued function of the state, whereas the propagation of a random vector yields a random process. Optimization-based strategies for calculating both types of predictor models are proposed. Whereas the formulations used to calculate IPMs target solutions leading to the interval value function of minimal spread containing all observations, those for RPMs seek to maximize the models' ability to reproduce the distribution of observations. Regarding RPMs, we choose a structure for the random vector (i.e., the assignment of probability to points in the parameter space) solely dependent on the prediction error. As such, the probabilistic description of uncertainty is not a subjective assignment of belief, nor is it expected to asymptotically converge to a fixed value, but instead it casts the model's ability to reproduce the experimental data. This framework enables evaluating the spread and distribution of the predicted response of target applications depending on the same parameters beyond the validation domain.

Crespo, Luis G.↗

Filter Tuning Using the Chi-Squared Statistic

The Goddard Space Flight Center (GSFC) Flight Dynamics Facility (FDF) performs orbit determination (OD) for the Aqua and Aura satellites. Both satellites are located in low Earth orbit (LEO), and are part of what is considered the A-Train satellite constellation. Both spacecraft are currently in the science phase of their respective missions. The FDF has recently been tasked with delivering definitive covariance for each satellite.The main source of orbit determination used for these missions is the Orbit Determination Toolkit developed by Analytical Graphics Inc. (AGI). This software uses an Extended Kalman Filter (EKF) to estimate the states of both spacecraft. The filter incorporates force modelling, ground station and space network measurements to determine spacecraft states. It also generates a covariance at each measurement. This covariance can be useful for evaluating the overall performance of the tracking data measurements and the filter itself. An accurate covariance is also useful for covariance propagation which is utilized in collision avoidance operations. It is also valuable when attempting to determine if the current orbital solution will meet mission requirements in the future.This paper examines the use of the Chi-square statistic as a means of evaluating filter performance. The Chi-square statistic is calculated to determine the realism of a covariance based on the prediction accuracy and the covariance values at a given point in time. Once calculated, it is the distribution of this statistic that provides insight on the accuracy of the covariance.For the EKF to correctly calculate the covariance, error models associated with tracking data measurements must be accurately tuned. Over estimating or under estimating these error values can have detrimental effects on the overall filter performance. The filter incorporates ground station measurements, which can be tuned based on the accuracy of the individual ground stations. It also includes measurements from the NASA space network (SN), which can be affected by the assumed accuracy of the TDRS satellite state at the time of the measurement.The force modelling in the EKF is also an important factor that affects the propagation accuracy and covariance sizing. The dominant force in the LEO orbit regime is the drag force caused by atmospheric drag. Accurate accounting of the drag force is especially important for the accuracy of the propagated state. The implementation of a box and wing model to improve drag estimation accuracy, and its overall effect on the covariance state is explored.The process of tuning the EKF for Aqua and Aura support is described, including examination of the measurement errors of available observation types (Doppler and range), and methods of dealing with potentially volatile atmospheric drag modeling. Predictive accuracy and the distribution of the Chi-square statistic, calculated based of the ODTK EKF solutions, are assessed versus accepted norms for the orbit regime.

Covariance Analysis↗