Optimum filters for deriving winds from FPS-16 radar/Jimsphere measurements
FPS-16/Jimsphere wind profiles measurement, discussing effect of data smoothing on accuracy and resolution
SEARCH · Engineering Papers
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.
FPS-16/Jimsphere wind profiles measurement, discussing effect of data smoothing on accuracy and resolution
The processing of data obtained from solar absorption radiation measurements is discussed. The position of the satellite was obtained by numerical integration of the differential equations of motion using initial conditions. The position of the sun was calculated as a function of time, and the tangential elevation was determined approximately from the positions of the satellite and the sun. The coefficients of an approximation formula and of a data smoothing process were determined, and the inversion of an Abel integral equation is solved analytically.
A computer program for studying linear ocean wave refraction is described. The program features random-access modular bathymetry data storage. Three bottom topography approximation techniques are available in the program which provide varying degrees of bathymetry data smoothing. Refraction diagrams are generated automatically and can be displayed graphically in three forms: Ray patterns with specified uniform deepwater ray density, ray patterns with controlled nearshore ray density, or crest patterns constructed by using a cubic polynomial to approximate crest segments between adjacent rays.
An algorithm is presented for generating an m-ary summation tree. The algorithm is completely general and may be applied to any length input string. For an N length sequence summed in groups of m sub l at each level l a maximum of 3L - 2 storage is required. A special case of the general m-ary tree where all m sub l are equal is used to smooth data in a radio frequency interference experiment. The maximum storage required when m l sub l = m for all l reduces to the closed form 3 log m N - 2.
Life prediction curves for rocket engine pump parts were developed from the results of high-cycle fatigue tests run on cast-aluminum specimens. Notched and smooth specimens were cyclically tested at different mean stress levels at -320 F (78 K). The notch size and mean stress enveloped the design operating conditions. Local stress computed in the groove of the notched specimen was used to represent its fatigue strength. The von Mises criterion was used to determine effective cyclic stresses. The Goodman rule was applied to determine equivalent reversed alternating stresses. The procedure permitted the notched and smooth data sets to each be described by a single curve. High-cycle fatigue life curves were provided for the stress state, mean stress, and stress concentration spanned by the data.
Image reduction used to simulate increase in altitude of an acquisition platform is equivalent to data smoothing, and can be achieved either by neighborhood averaging or by filtering techniques. The averaging approach is limited for accurate simulation. A filtering method is described which was based on the hypothesis that all changes due to altitude increase can be represented by a point spread function. Determination of the scale function and factor are discussed as well as implementation of the filtering. Filtering can be either in the spatial or frequency domain. In the spatial domain, filtering consists of the convolution of the image with the weights mask, and then of the declination of the points according to the appropriates scale factor. A simulation of an aircraft day image in the infrared channel is examined.
The utilization of LANDSAT data in forestry inventories in California is summarized. Problems in selecting meaningful classification systems, utilization of supervised or unsupervised classification systems, and data smoothing are discussed.
An examination of limitations, requirements, and precision of the linear multiple-regression technique for quantification of marine environmental parameters is conducted. Both environmental and optical physics conditions have been defined for which an exact solution to the signal response equations is of the same form as the multiple regression equation. Various statistical parameters are examined to define a criteria for selection of an unbiased fit when upwelled radiance values contain error and are correlated with each other. Field experimental data are examined to define data smoothing requirements in order to satisfy the criteria of Daniel and Wood (1971). Recommendations are made concerning improved selection of ground-truth locations to maximize variance and to minimize physical errors associated with the remote sensing experiment.
Following reception and analog to digital conversion (A/D) conversion, atmospheric radar backscatter echoes need to be processed so as to obtain desired information about atmospheric processes and to eliminate or minimize contaminating contributions from other sources. Various signal processing techniques have been implemented at mesosphere-stratosphere-troposphere (MST) radar facilities to estimate parameters of interest from received spectra. Such estimation techniques need to be both accurate and sufficiently efficient to be within the capabilities of the particular data-processing system. The various techniques used to parameterize the spectra of received signals are reviewed herein. Noise estimation, electromagnetic interference, data smoothing, correlation, and the Doppler effect are among the specific points addressed.
Noteworthy properties of one-dimensional splines are presented. Stochastic processes, kernel functions, interpolation, and data smoothing were discussed. Theorems and concomitant proofs are provided.
A multiple high-order derivative analysis algorithm has been developed which can automatically extract absorption band positions from low-quality reflectance spectra with little degredation of accuracy. Overlapping bands with comparable widths and intensities can be resolved whose centers are as close as 0.3-0.5 W, with safer resolution limits of 0.6-1.0 W band center separations suggested for overlapping bands that are dissimilar. The segment length for smoothing is continually adjusted to about 0.5 W to minimize signal distortion, and a spectral pattern recognition algorithm predicts the signal spectrum and calculates approximate W across the spectrum using its second derivative. A single-pass cubic spline is applied to the smoothed data, and a sliding segment sixth-order polynomial is fit to the spectrum, with the length of the segment being continuously locally adjusted to 1.0 W across the spectrum. Good reliability and consistency of the algorithm is demonstrated with application to laboratory and earth-based telescope spectra.
The focus is on the lower temperature, uncoated and notched features of gas turbine blades. Constitutive and fatigue life prediction models applicable to these regions are being developed. Fatigue results are presented which were obtained thus far. Fatigue tests are being conducted on PWA 1480 single crystal material using smooth strain controlled specimens and three different notched specimens. Isothermal fatigue tests were conducted at 1200, 1400, and 1600 F. The bulk of the tests were conducted at 1200 F. The strain controlled tests were conducted at 0.4 percent per second strain rate and the notched tests were cycled at 1.0 cycle per second. A clear orientation dependence is observed in the smooth strain controlled fatigue results. The fatigue lifes of the thin, mild notched specimens agree fairly well with this smooth data when elastic stress range is used as a correlating parameter. Finite element analyses were used to calculate notch stresses. Fatigue testing will continue to further explore the trends observed thus far. Constitutive and life prediction models are being developed.
A tool which generates optimal trajectory/control histories in an integrated manner is generically adapted to the treatment of single-stage-to-orbit air-breathing hypersonic vehicles. The methodology is implemented as a two point boundary value problem solution technique. Its use permits an assessment of an entire near-minimum-fuel trajectory and desired control strategy from takeoff to orbit while satisfying physically derived inequality constraints and while achieving efficient propulsive mode phasing. A simpler analysis strategy that partitions the trajectory into several boundary condition matched segments is also included to construct preliminary trajectory and control history representations with less computational burden than is required for the overall flight profile assessment. A demonstration was accomplished using a tabulated example (winged-cone accelerator) vehicle model that is combined with a newly developed multidimensional cubic spline data smoothing routine. A constrained near-fuel-optimal trajectory, imposing a dynamic pressure limit of 1000 psf, was developed from horizontal takeoff to 20,000 ft/sec relative air speed while aiming for a polar orbit. Previously unspecified propulsive discontinuities were located. Flight regimes demanding rapid attitude changes were identified, dictating control effector and closed-loop controller authority was ascertained after evaluating effector use for vehicle trim. Also, inadequacies in vehicle model representations and specific subsystem models with insufficient fidelity were determined based on unusual control characteristics and/or excessive sensitivity to uncertainty.
We use the first accurate measurements of current densities in the plasma sheet to calculate the half-thickness and position of the current sheet as a function of time. Our technique assumes a Harris current sheet model, which is parameterized by lobe magnetic field B(o), current sheet half-thickness h, and current sheet position z(sub o). Cluster measurements of magnetic field, current density, and plasma pressure are used to infer the three parameters as a function of time. We find that most long timescale (6-12 hours) current sheet crossings observed by Cluster cannot be described by a static Harris current sheet with a single set of parameters B(sub o), h, and z(sub o). Noting the presence of high-frequency fluctuations that appear to be superimposed on lower frequency variations, we average over running 6-min intervals and use the smoothed data to infer the parameters h(t) and z(sub o)(t), constrained by the pressure balance lobe magnetic field B(sub o)(t). Whereas this approach has been used in previous studies, the spatial gnuhen& now provided by the Cluster magnetometers were unavailable or not well constrained in earlier studies. We place the calculated hdf&cknessa in a magnetospheric context by examining the change in thickness with substorm phase for three case study events and 21 events in a superposed epoch analysis. We find that the inferred half-thickness in many cases reflects the nominal changes experienced by the plasma sheet during substorms (i.e., thinning during growth phase, thickening following substorm onset). We conclude with an analysis of the relative contribution of (Delta)B(sub z)/(Delta)X to the cross-tail current density during substorms. We find that (Delta)B(sub z)/(Delta)X can contribute a significant portion of the cross-tail c m n t around substorm onset.
This workshop presentation discusses the design and implementation of numerical methods for the quantification of statistical uncertainty, including a-posteriori error bounds, for output quantities computed using CFD methods. Hydrodynamic realizations often contain numerical error arising from finite-dimensional approximation (e.g. numerical methods using grids, basis functions, particles) and statistical uncertainty arising from incomplete information and/or statistical characterization of model parameters and random fields. The first task at hand is to derive formal error bounds for statistics given realizations containing finite-dimensional numerical error [1]. The error in computed output statistics contains contributions from both realization error and the error resulting from the calculation of statistics integrals using a numerical method. A second task is to devise computable a-posteriori error bounds by numerically approximating all terms arising in the error bound estimates. For the same reason that CFD calculations including error bounds but omitting uncertainty modeling are only of limited value, CFD calculations including uncertainty modeling but omitting error bounds are only of limited value. To gain maximum value from CFD calculations, a general software package for uncertainty quantification with quantified error bounds has been developed at NASA. The package provides implementations for a suite of numerical methods used in uncertainty quantification: Dense tensorization basis methods [3] and a subscale recovery variant [1] for non-smooth data, Sparse tensorization methods[2] utilizing node-nested hierarchies, Sampling methods[4] for high-dimensional random variable spaces.
The compiled near-surface data and satellite crustal magnetic measured data are modeled with a regionally concentrated spherical harmonic presentation technique over Australia and Antarctica. Global crustal magnetic anomaly studies have used a spherical harmonic analysis to represent the Earth's magnetic crustal field. This global approach, however is best applied where the data are uniformly distributed over the entire Earth. Satellite observations generally meet this requirement, but unequally distributed data cannot be easily adapted in global modeling. Even for the satellite observations, due to the errors spread over the globe, data smoothing is inevitable in the global spherical harmonic presentations. In addition, global high-resolution modeling requires a great number of global spherical harmonic coefficients for the regional presentation of crustal magnetic anomalies, whereas a lesser number of localized spherical coefficients will satisfy. We compared methods in both global and regional approaches and for a case where the errors were propagated outside the region of interest. For observations from the upcoming Swarm constellation, the regional modeling will allow the production a lesser number of spherical coefficients that are relevant to the region of interest
Matlab software was written to provide smoothing of radar tracking data to simulate ADS-B (Automatic Dependent Surveillance-Broadcast) data in order to test a tactical conflict probe. The probe, called TSAFE (Tactical Separation-Assured Flight Environment), is designed to handle air-traffic conflicts left undetected or unresolved when loss-of-separation is predicted to occur within approximately two minutes. The data stream that is down-linked from an aircraft equipped with an ADS-B system would include accurate GPS-derived position and velocity information at sample rates of 1 Hz. Nation-wide ADS-B equipage (mandated by 2020) should improve surveillance accuracy and TSAFE performance. Currently, position data are provided by Center radar (nominal 12-sec samples) and Terminal radar (nominal 4.8-sec samples). Aircraft ground speed and ground track are estimated using real-time filtering, causing lags up to 60 sec, compromising performance of a tactical resolution tool. Offline smoothing of radar data reduces wild-point errors, provides a sample rate as high as 1 Hz, and yields more accurate and lag-free estimates of ground speed, ground track, and climb rate. Until full ADS-B implementation is available, smoothed radar data should provide reasonable track estimates for testing TSAFE in an ADS-B-like environment. An example illustrates the smoothing of radar data and shows a comparison of smoothed-radar and ADS-B tracking. This document is intended to serve as a guide for using the smoothing software.
Charged particle calibrations based on differenced range versus integrated Doppler (DRVID) are smoothed and fitted with a polynomial prior to its application in the orbit determination process. A description is given of the results of the tests performed on the computer program performing these calculations to determine its characteristics and to evaluate the acceptability of its computations. This program, called MEDIA, is described.