Search NASA⌕ Search

SEARCH · Search NASA

Results for “data skew”

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 145 records · Page 8

On-orbit control system performance of the HEAO-2 observatory

The on-orbit performance of the High Energy Astronomy Observatory is described. The control system utilizes precision gyros for attitude reference and skewed reaction wheels controlled by a reprogrammable computer. The observatory points at selected targets, maneuvers automatically and acquires guide stars for updating at each target. Performance data indicates maneuver accuracies less than one arc minute and pointing accuracies of 2-5 arc seconds. Use of 7th-9th magnitude guide stars has resulted in several improper updates caused by 'flat-field' stray light effects in the star trackers. The method used to discriminate true stars from the flat-field is also presented.

Rose, R. E.↗

Statistical characteristics of velocity, concentration, mass transport, and momentum transport for coaxial jet mixing in a confined duct

An experimental study of mixing downstream of coaxial jets discharging into an expanded circular duct was conducted to obtain data for the evaluation and improvement of turbulent transport models currently used for combustor flow modeling. A combination of turbulent momentum transport rate and two velocity component data was obtained from simultaneous measurements with a two-color LV system. A combination of turbulent mass transport rate, concentration and velocity data was obtained from simultaneous measurements with laser velocimeter (LV) and laser induced fluorescence (LIF) systems. These measurements were used to obtain mean, second central moment, skewness and kurtosis values for three velocity components and the concentration. These measurements showed the existence of countergradient turbulent axial mass transport where the annular jet fluid was accelerating the inner jet fluid. Results from the study are related to the assumptions employed in the current mass and momentum turbulent transport models.

Johnson, B. V.↗

Calibration of Direct Velocimetry Using Hot Wire Anemometry

We studied the functionality of the IFA300 Hot Wire Anemometry system and its feasibility for the calibration of Direct Estimator Velocimetry (DEV) using Laser Induced Florescence (LIF). The experimental setup consisted of a single-wire hot wire probe measuring one component of velocity, the anemometry computational hardware and software, an HP oscilloscope, and a table fan to produce a simple flow with fluctuations. Measurements were taken at several points in the stream wise and transverse directions in the flow and various parameters recorded such as mean velocities temperature, turbulence intensities, skewness coefficients and flatness coefficients. The IFA300 software also allowed us to perform statistical manipulations such as spectrum analysis on velocities samples and correlation. Utilization of data files, also produced by the anemometry software, and post analysis were performed to produce graphical representations of turbulent intensity versus probe position and a flow field velocity profile. We concluded that the IFA300 Hot Wire Anemometry system is a reliable and functional method for calibration of DEV using LIF. Our future intentions are to set up a test chamber such that both velocity measurement techniques can be applied simultaneously, thus the calibration.

Holland, Michael L.↗

On-Orbit Calibration of Redundant Spacecraft Gyros by Optimal Reduction to Three Axes

The Aqua spacecraft will carry four single-axis gyros configured with three orthogonal axes and one skew axis. This redundancy presents a challenge for batch methods of on-orbit gyro calibration that use a spacecraft rotation model deterministically related to gyro data, in that sensor data can respond to at most three angular velocity components. When the number of gyros, N, is greater than 3, the 3xN matrix, G, that reduces the N gyro measurements to three body-frame angular-velocity components cannot be fully determined by such methods; there are many such matrices that produce essentially the same angular velocity history. In such a case, spacecraft operators require information about the Nx3 gyro linear response matrix, R, that relates gyro outputs to the body-frame angular velocities causing them. This matrix provides sufficient information to determine multiple reduced-dimension G-matrices for use in case of failure or degradation of one or more gyros, as well as to determine an optimal 3xN G for the fully-functional configuration. A viable proposal is to apply a 3xN pre-filter matrix, F, to the N gyro outputs before carrying out a conventional gyro calibration procedure. The angular-velocity history emerging from conventional calibration may then be used as input data, together with the same gyro data that generated it, to fit the alignment, scale-factor, and bias parameters of each gyro axis in turn. A difficulty of such a proposal is the arbitrariness in the choice of F. Due to gyro noise, different pre-filter matrices produce different calibrations. This paper presents a method of choosing F that is based on optimizing gyro consistency in the limit of infinite weight on gyro data, as compared to sensor data. The choice of F is independent of a priori alignment and is based on the gyro data alone. The method is applicable to any N of three or more, but reduces to conventional batch-estimation methodologies when N = 3. Results of computational comparison among calibration simulations using various choices of F will be presented for the Aqua gyro configuration with N = 4.

Radomski, M. S.↗

A comparison of CMG steering laws for High Energy Astronomy Observatories (HEAOs)

A comparison of six selected control moment gyro steering laws for use on the HEAO spacecraft is reported. Basic equations are developed to project the momentum and torque of four skewed, single gimbal CMGs into vehicle coordinates. In response to the spacecraft attitude error signal, six algorithms are derived for controlling the CMG gimbal movements. HEAO performance data are obtained using each steering law and compared on the basis of such factors as accuracy, complexity, singularities, gyro hang-up and failure adaption. Moreover, each law is simulated with and without a magnetic momentum management system. The performance of any steering law is enhanced by the magnetic system. Without magnetics, the gimbal angles get large and there are significant differences in steering law performances due to cross coupling and nonlinearities. The performance of the pseudo inverse law is recommended for HEAO.

Davis, B. G.↗

Random Variables with Moment-Matching Staircase Density Functions

This paper proposes a family of random variables for uncertainty modeling. The variables of interest have a bounded support set, and prescribed values for the first four moments. We present the feasibility conditions for the existence of any of such variables, and propose a class of variables that conforms to such constraints. This class is called staircase because the density of its members is a piecewise constant function. Convex optimization is used to calculate their distributions according to several optimality criteria, including maximal entropy and maximal log-likelihood. The flexibility and efficiency of staircases enable modeling phenomena having a possibly skewed and/or multimodal response at a low computational cost. Furthermore, we provide a means to account for the uncertainty in the distribution caused by estimating staircases from data. These ideas are illustrated by generating empirical staircase predictor models. We consider the case in which the predictor matches the sample moments exactly (a setting applicable to large datasets), as well as the case in which the predictor accounts for the sampling error in such moments (a setting applicable to sparse datasets). A predictor model for the dynamics of an aeroelastic airfoil subject to flutter instability is used as an example. The resulting predictor not only describes the system's response accurately, but also enables carrying out a risk analysis for safe flight.

Luis G. Crespo↗

Bayesian Rules of Thumb: Robust Uncertainty Quantification in Early Project Cost Estimation

Systems engineers often make use of cost Rules ofThumb in order to estimate cost during early phases of projectformulation. These Rules of Thumb typically take the form ofa sequence of percentages over which a total cost is allocatedacross NASA WBS elements. Rules of Thumb can then be usedto extrapolate cost from one or more known WBS elements tothe remaining unknown WBS elements, assisting early projectformulation architecture studies (such as those in JPL’s Team Xand A Team).A number of issues can arise when generating and using costRules of Thumb. For example, many records of project costsconsist of incomplete data. Typical methods of dealing withincomplete cost allocation data include (a) ignoring missionswith incomplete data, or (b) taking averages of the non-zero percentagesacross missions, but both of these methods can result inbiased estimates if the existence of incomplete data correlateswith total mission cost or any particular WBS element. Anothercommon example is cost reported in one or more incorrect WBSelements. This is especially prevalent in smaller missions whereit is more common for engineers to perform tasks that fall underthe purview of multiple WBS elements.Furthermore, a Rule of Thumb estimate is typically reported asa point estimate; there is no reported uncertainty around thepercentages used to generate an allocation. Even in the rarecase in which confidence intervals around mean percentages areprovided, there may be positive or negative correlations betweenWBS elements which can skew estimates.Here we attempt to address these problems by formulatingprobabilistic Rules of Thumb in which a distribution of allocationschemes, rather than a single allocation scheme, is generated.We use a bootstrap imputation method to simultaneouslyaccount for uncertainty in the missing data while using allavailable information contained in the dataset. The imputeddatasets are then input into a multivariate Bayesian modelwhich accounts for correlations between WBS elements andproperly accounts for uncertainty in the final Rule of Thumbpercentages and predictions. We describe the mathematicalmodel and provides snippets of R code utilizing the brms(Bayesian Regression Models using Stan) package. To illustratethis model, we generate a Bayesian Level 2 WBS Cost Rule ofThumb for MIDEX (Medium-Class Explorers) missions withdata extracted from NASA’s CADRe. We then compare thismethod’s performance with the classical Rule of Thumb method.

Hooke, Melissa A↗

Hydrodynamic impeller stiffness, damping, and inertia in the rotordynamics of centrifugal flow pumps

The lateral hydrodynamic forces experienced by a centrifugal pump impeller performing circular whirl motions within several volute geometries were measured. The lateral forces were decomposed into: (1) time averaged lateral forces and (2) hydrodynamic force matrices representing the variation of the lateral forces with position of the impeller center. It is found that these force matrices essentially consist of equal diagonal terms and skew symmetric off diagonal terms. One consequence of this is that during its whirl motion the impeller experiences forces acting normal and tangential to the locus of whirl. Data on these normal and tangential forces are presented; it is shown that there exists a region of positive reduced whirl frequencies, within which the hydrodynamic forces can be destablizing with respect to whirl.

Jery, S.↗

Parallel VLSI architecture emulation and the organization of APSA/MPP

The Applicative Programming System Architecture (APSA) combines an applicative language interpreter with a novel parallel computer architecture that is well suited for Very Large Scale Integration (VLSI) implementation. The Massively Parallel Processor (MPP) can simulate VLSI circuits by allocating one processing element in its square array to an area on a square VLSI chip. As long as there are not too many long data paths, the MPP can simulate a VLSI clock cycle very rapidly. The APSA circuit contains a binary tree with a few long paths and many short ones. A skewed H-tree layout allows every processing element to simulate a leaf cell and up to four tree nodes, with no loss in parallelism. Emulation of a key APSA algorithm on the MPP resulted in performance 16,000 times faster than a Vax. This speed will make it possible for the APSA language interpreter to run fast enough to support research in parallel list processing algorithms.

Odonnell, John T.↗

Decay of Isotropic Turbulence at Low Reynolds Number

Decay of isotropic turbulence is computed using direct numerical simulations. Comparisons with experimental spectra at moderate and low Reynolds numbers (R(sub lambda) less than 70) show good agreement. At moderate to high Reynolds numbers (R(sub lambda) greater 50), the spectra are found to collapse with Kolmogorov scaling at high wave numbers. However, at low Reynolds numbers (R(sub lambda) less than 50) the shape of the spectra at the Kolmogorov length scales is Reynolds number dependent. Direct simulation data from flowfields of decaying isotropic turbulence are used to compute the terms in the equation for the dissipation rate of the turbulent kinetic energy. The development of the skewness and the net destruction of the turbulence dissipation rate in the limit of low Reynolds numbers are presented. The nonlinear terms are found to remain active at surprisingly low Reynolds numbers.

Mansour, N. N.↗

Measurement and Modeling of Steep Ocean Wave Slopes

Our study emphasizes the importance of identifying and quantifying the distribution variance, skewness and kurtosis from optical and microwave scattering observations. Recent field measurements of the sea slope distribution for intermediate-to-long scale gravity waves will be presented. These data were collected using an airborne laser range system designed to estimate the surface slope vector at horizontal scales of 1-2 m. The observed slope distribution tail indicates that the occurrence of steep waves substantially exceeds a Gaussian prediction. This measured peakedness is present over the wide range of sea state and wind speed conditions encountered. Data are further evaluated within the context of Cox and Munk's well-known sea slope investigations. Based on a re-evaluation of the Cox and Munk's reported parameters, we find a consistent picture develops wherein data are shown to consistently indicate non-Gaussian statistics. One fundamental application of such a non-Gaussian slope observation is its place in modifying predicted wave breaking probability to help to better quantify gas transfer processes at the sea surface.

Chapron, B.↗

Magnetic Field Strength Fluctuations in the Heliosheath: Voyager 1 Observations during 2009

We analyze the "microscale fluctuations" of the magnetic field strength Beta on a scale of several hours observed by Voyager I (VI) in the heliosheath during 2009. The microscale fluctuations of Beta range from coherent to stochastic structures. The amplitude of microscale fluctuations of Beta during 1 day is measured by the standard deviation (SD) of 48 s averages of B. The distribution of the daily values of SD is lognormal. SD(t) from day of year (DOY) I to 331, 2009, is very intermittent. SD(t) has a 1/f or "pink noise" spectrum on scales from I to 100 days, and it has a broad multi fractal spectrum f(alpha) with 0.57 much < alpha much < 1.39. The time series of increments SD(t + r) - SD(t) has a pink noise spectrum with alpha' = 0.88 +/- 0.14 on scales from 1 to 100 days. The increments have a Tsallis (q-Gaussian) distribution on scales from 1 to 165 days, with an average q = 1.75 +/- 0.12. The skewness S and kurtosis K have Gaussian and lognormal distributions, respectively. The largest spikes in K(t) and Set) are often associated with a change in Beta across a data gap and with identifiable physical structures. The "turbulence" observed by VI during 2009 was weakly compressible on average but still very intermittent, highly variable, and highly compressible at times. The turbulence observed just behind the termination shock by Voyager 2 was twice as strong. These observations place strong constraints on any model of "turbulence" in the heliosheath.

Burlaga, L. F.↗

Magnetic Field Strength Fluctuations in the Heliosheath: Voyager 1 Observations During 2009

We analyze the microscale fluctuations of the magnetic field strength B on a scale of several hours observed by Voyager1 (V1) in the heliosheath during 2009. The microscale fluctuations of B range from coherent to stochastic structures. The amplitude of microscale fluctuations of B during 1 day is measured by the standard deviation (SD) of 48 s averages of B. The distribution of the daily values of SD is lognormal. SD(t) from day of year (DOY) 1 to 331, 2009, is very intermittent. SD(t) has a 1/f or "pink noise" spectrum on scales from 1 to 100 days, and it has a broad multifractal spectrum f(alpha) with 0.57 less than or equal to alpha less than or equal to 1.39. The time series of increments SD(t + tau) -- SD(t) has a pink noise spectrum with alpha(1) = 0.88 plus or minus 0.14 on scales from 1 to 100 days. The increments have a Tsallis (q-Gaussian) distribution on scales from 1 to 165 days, with an average q = 1.75 plus or minus 0.12. The skewness S and kurtosis K have Gaussian and lognormal distributions, respectively. The largest spikes in K(t) and S(t) are often associated with a change in B across a data gap and with identifiable physical structures. The "turbulence" observed by V1 during 2009 was weakly compressible on average but still very intermittent, highly variable, and highly compressible at times. The turbulence observed just behind the termination shock by Voyager 2 was twice as strong. These observations place strong constraints on any model of turbulence in the heliosheath.

Brulaga, L. F.↗

Shuttle active thermal control system development testing. Volume 2: Modular radiator system tests

Tests were designed to investigate the validity of the "modular" approach to space radiator system design for space shuttle and future applications by gathering performance data on various systems comprised of different numbers of identical panels, subject to nominal and extreme heat loads and environments. Both one-sided and two-sided radiation was tested, and engineering data was gathered on simulated low a/e coatings and system response to changes in outlet temperature control point. The results of the testing showed system stability throughout nominal orbital transients, unrealistically skewed environments, freeze-thaw transients, and rapid changes in outlet temperature control point. Various alternative panel plumbing arrangements were tested with no significant changes in performance being observed. With the MRS panels arranged to represent the shuttle baseline system, a maximum heat rejection of 76,600 Btu/hr was obtained in segmented tests under the expected worst case design environments. Testing of an alternate smaller two-sided radiation configuration yielded a maximum heat rejection of 52,931 Btu/hr under the maximum design environments.

Scheps, P. B.↗

Simulation of spatially evolving turbulence and the applicability of Taylor's hypothesis in compressible flow

For the numerical simulation of inhomogeneous turbulent flows, a method is developed for generating stochastic inflow boundary conditions with a prescribed power spectrum. Turbulence statistics from spatial simulations using this method with a low fluctuation Mach number are in excellent agreement with the experimental data, which validates the procedure. Turbulence statistics from spatial simulations are also compared to those from temporal simulations using Taylor's hypothesis. Statistics such as turbulence intensity, vorticity, and velocity derivative skewness compare favorably with the temporal simulation. However, the statistics of dilatation show a significant departure from those obtained in the temporal simulation. To directly check the applicability of Taylor's hypothesis, space-time correlations of fluctuations in velocity, vorticity, and dilatation are investigated. Convection velocities based on vorticity and velocity fluctuations are computed as functions of the spatial and temporal separations. The profile of the space-time correlation of dilatation fluctuations is explained via a wave propagation model.

Lee, Sangsan↗

Interference Lattice-based Loop Nest Tilings for Stencil Computations

A common method for improving performance of stencil operations on structured multi-dimensional discretization grids is loop tiling. Tile shapes and sizes are usually determined heuristically, based on the size of the primary data cache. We provide a lower bound on the numbers of cache misses that must be incurred by any tiling, and a close achievable bound using a particular tiling based on the grid interference lattice. The latter tiling is used to derive highly efficient loop orderings. The total number of cache misses of a code is the sum of (necessary) cold misses and misses caused by elements being dropped from the cache between successive loads (replacement misses). Maximizing temporal locality is equivalent to minimizing replacement misses. Temporal locality of loop nests implementing stencil operations is optimized by tilings that avoid data conflicts. We divide the loop nest iteration space into conflict-free tiles, derived from the cache miss equation. The tiling involves the definition of the grid interference lattice an equivalence class of grid points whose images in main memory map to the same location in the cache-and the construction of a special basis for the lattice. Conflicts only occur on the boundaries of the tiles, unless the tiles are too thin. We show that the surface area of the tiles is bounded for grids of any dimensionality, and for caches of any associativity, provided the eccentricity of the fundamental parallelepiped (the tile spanned by the basis) of the lattice is bounded. Eccentricity is determined by two factors, aspect ratio and skewness. The aspect ratio of the parallelepiped can be bounded by appropriate array padding. The skewness can be bounded by the choice of a proper basis. Combining these two strategies ensures that pathologically thin tiles are avoided. They do not, however, minimize replacement misses per se. The reason is that tile visitation order influences the number of data conflicts on the tile boundaries. If two adjacent tiles are visited successively, there will be no replacement misses on the shared boundary. The iteration space may be covered with pencils larger than the size of the cache while avoiding data conflicts if the pencils are traversed by a scanning-face method. Replacement misses are incurred only on the boundaries of the pencils, and the number of misses is minimized by maximizing the volume of the scanning face, not the volume of the tile. We present an algorithm for constructing the most efficient scanning face for a given grid and stencil operator. In two dimensions it is based on a continued fraction algorithm. In three dimensions it follows Voronoi's successive minima algorithm. We show experimental results of using the scanning face, and compare with canonical loop orderings.

VanderWijngaart, Rob F.↗

Two-stage fan. 3: Data and performance with rotor tip casing treatment, uniform and distorted inlet flows

A two stage fan with a 1st-stage rotor design tip speed of 1450 ft/sec, a design pressure ratio of 2.8, and corrected flow of 184.2 lbm/sec was tested with axial skewed slots in the casings over the tips of both rotors. The variable stagger stators were set in the nominal positions. Casing treatment improved stall margin by nine percentage points at 70 percent speed but decreased stall margin, efficiency, and flow by small amounts at design speed. Treatment improved first stage performance at low speed only and decreased second stage performance at all operating conditions. Casing treatment did not affect the stall line with tip radially distorted flow but improved stall margin with circumferentially distorted flow. Casing treatment increased the attenuation for both types of inlet flow distortion.

Burger, G. D.↗

A statistical model for interpreting computerized dynamic posturography data

Computerized dynamic posturography (CDP) is widely used for assessment of altered balance control. CDP trials are quantified using the equilibrium score (ES), which ranges from zero to 100, as a decreasing function of peak sway angle. The problem of how best to model and analyze ESs from a controlled study is considered. The ES often exhibits a skewed distribution in repeated trials, which can lead to incorrect inference when applying standard regression or analysis of variance models. Furthermore, CDP trials are terminated when a patient loses balance. In these situations, the ES is not observable, but is assigned the lowest possible score--zero. As a result, the response variable has a mixed discrete-continuous distribution, further compromising inference obtained by standard statistical methods. Here, we develop alternative methodology for analyzing ESs under a stochastic model extending the ES to a continuous latent random variable that always exists, but is unobserved in the event of a fall. Loss of balance occurs conditionally, with probability depending on the realized latent ES. After fitting the model by a form of quasi-maximum-likelihood, one may perform statistical inference to assess the effects of explanatory variables. An example is provided, using data from the NIH/NIA Baltimore Longitudinal Study on Aging.

NASA Discipline Neuroscience↗