Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Gaussian simulation”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Head-on collision of viscous vortex rings

The head-on collision of two identical axisymmetric viscous vortex rings is studied through direct simulations of the incompressible Navier-Stokes equations. The initial vorticity distributions considered are those of Hill's spherical vortex and of rings with circular Gaussian cores, each at Reynolds numbers of about 350 and 1000. The Reynolds number is defined by Gamma/Nu, the ratio of circulation to viscosity. As the vortices approach each other by self-induction, the radii increase by mutual induction, and vorticy cancels through viscous cross-diffusion across the collision plane. Following contact, the vorticity distribution in the core forms a head-tail structure (for the cases considered). The characteristic time of vorticity annihilation is compared with that of a 3D collision experiment and 3D numerical simulations. It is found that the annihilation time is somewhat longer in the axisymmetric case than it is in the symmetry plane of the experiment and 3D numerical simulation. By comparing the annihilatiom time with a viscous timescale and a circulation timescale, it is deduced that both the strain rate due to local effects and to 3D vorticity realignment are important.

Stanaway, S.↗
NASA NTRS1987

Roundoff noise analysis for digital signal power processors using Welch's power spectrum estimation

The noise due to finite-word-length effects is analyzed for digital-signal power processors using Welch's power-spectrum estimation technique to measure the power of Gaussian random signals over a frequency band of interest. The input of the digital signal processor contains a finite-length time interval in which the true Gaussian signal is contaminated by Gaussian noise. The roundoff noise-to-signal ratio in the measurement of the signal power is derived, and computer simulations which validate the analytical results are presented. These results can be used in tradeoff studies of hardware design, such as the number of bits required at each processing stage. The results presented in this paper are currently being used in the design of a digital Doppler processor (Chi et al., 1986) for a radar scatterometer.

Chi, Chong-Yung↗
NASA NTRS1988

Analysis of the light scattering interferometry for spheres larger than the light wavelength

A theoretical model, based on the geometrical optics approach, has been developed to simulate various aspects of the phase Doppler particle analyzer (PDPA). The model has taken into consideration the non-uniform (Gaussian) illumination of the particles as they pass through the probe volume. Instrument response curves have been generated for various scattering angles by performing spatial and temporal integration of the scattered intensity distribution over the receiver surface. Experimental and theoretical investigations have established the applicability of this instrument to both forward and backscatter angles.

Bachalo, W. D.↗
NASA NTRS1994

Probing the large-scale velocity field with clusters of galaxies

What is the role of clusters of galaxies in probing the large-scale velocity field of the universe? We investigate the distribution of peculiar velocities of clusters of galaxies in the popular low-density (omega = 0.3) flat cold dark matter (CDM) cosmological model, which best fits many large-scale structure observations. An omega = 1 CDM model is also studied for comparison. We find that clusters of galaxies are efficient tracers of the large-scale velocity field. The clusters exhibit a Maxwellian distribution of peculiar velocities, as expected from Gaussian initial density fluctuations. The cluster three-dimensional velocity distribution for the omega = 0.3 model peaks at nu approximately greater than 400 km/s and extends to high velocities of nu approximately 1200 km/s. The rms peculiar velocity of the clusters is 440 km/s. Approximately 10% of all model clusters move with high peculiar velocities nu greater or equal to 700 km/s. The observed velocity distribution of clusters of galaxies is compared with the predictions from cosmological models. The observed data exhibit a larger velocity tail than seen in the model simulations; however, due to the large observational uncertainties, the data are consistent at approximately equal to 3 sigma level with the odel predictions, and with a Gaussian initial density field. The large peculiar velocities reported for some clusters of galaxies (nu approximately greater than 3000 km/s) are likely to be overestimated, if the current model is viable.

Bahcall, Neta A.↗
NASA NTRS1999

Large Eddy Simulation of Aircraft Wake Vortices in a Homogeneous Atmospheric Turbulence: Vortex Decay and Descent

The effects of ambient turbulence on decay and descent of aircraft wake vortices are studied using a validated, three-dimensional: large-eddy simulation model. Numerical simulations are performed in order to isolate the effect of ambient turbulence on the wake vortex decay rate within a neutrally-stratified atmosphere. Simulations are conducted for a range of turbulence intensities, by injecting wake vortex pairs into an approximately homogeneous and isotropic turbulence field. The decay rate of the vortex circulation increases clearly with increasing ambient turbulence level, which is consistent with field observations. Based on the results from the numerical simulations, simple decay models are proposed as functions of dimensionless ambient turbulence intensity (eta) and dimensionless time (T) for the circulation averaged over a range of radial distances. With good agreement with the numerical results, a Gaussian type of vortex decay model is proposed for weak turbulence: while an exponential type of Tortex decay model can be applied for strong turbulence. A relationship for the vortex descent based on above vortex decay model is also proposed. Although the proposed models are based on simulations assuming neutral stratification, the model predictions are compared to Lidar vortex measurements observed during stable, neutral, and unstable atmospheric conditions. In the neutral and unstable atmosphere, the model predictions appear to be in reasonable agreement with the observational data, while in the stably-stratified atmosphere, they largely underestimate the observed circulation decay with consistent overestimation of the observed vortex descent. The underestimation of vortex decay during stably-stratified conditions suggests that stratification has an important influence on vortex decay when ambient levels of turbulence are weak.

Han, Jongil↗
NASA NTRS1994

Resolution requirements for velocity gradients in turbulence

Since high resolution numerical simulations of turbulent flows, or experiments at high Reynolds numbers, represent a substantial investment in resources, the estimation of the minimum resolution required for the study of a given property has been the subject of continued interest. Different properties require, in general, different resolutions, and the present paper is dedicated to the requirements for the measurement of the probability distribution functions of the velocity gradients and, in particular, of their low order moments. The deviation of these quantities from the values corresponding to a Gaussian distribution was one of the first indications of the presence of Reynolds number-dependent intermittency and has been the object of recent interest as numerical simulations have become able to explore the distribution of gradients in the low Reynolds number range, while new experiments have extended the range to increasingly high Reynolds numbers. We will use progressive filtering of the results of numerical simulations of isotropic turbulence as a model for the effect of a sensor of finite size. The numerical issues will be addressed first to ensure that the simulations are fully resolved from the point of view of the velocity gradients. This will also give us an estimate for the numerical resolution required for the different quantities.

Jimenez, Javier↗
NASA NTRS1972

Digital simulation of continuous error models with application to an instrument landing system error

A digital simulation of the continuous error of the localized beam of a conventional instrument landing system is discussed. The digital simulation was developed during the analysis of space shuttle navigation capabilities. A discrete mathematical model for use on a digital computer is described. The model generates an output random sequence which is equivalent, for simulation purposes, to the desired random process. The model is a system of difference equations driven by a zero-mean Gaussian random sequence.

Merrick, R. B.↗
NASA NTRS1977

A theoretical study of heterojunction and graded band gap type solar cells

A computer program was designed for the analysis of variable composition solar cells and applied to several proposed solar cell structures using appropriate semiconductor materials. The program simulates solar cells made of a ternary alloy of two binary semiconductors with an arbitrary composition profile, and an abrupt or Gaussian doping profile of polarity n-on-p or p-on-n with arbitrary doping levels. Once the device structure is specified, the program numerically solves a complete set of differential equations and calculates electrostatic potential, quasi-Fermi levels, carrier concentrations and current densities, total current density and efficiency as functions of terminal voltage and position within the cell. These results are then recorded by computer in tabulated or plotted form for interpretation by the user.

Sutherland, J. E.↗
NASA NTRS1991

Coding gains and error rates from the Big Viterbi Decoder

A prototype hardware Big Viterbi Decoder (BVD) was completed for an experiment with the Galileo Spacecraft. Searches for new convolutional codes, studies of Viterbi decoder hardware designs and architectures, mathematical formulations, and decompositions of the deBruijn graph into identical and hierarchical subgraphs, and very large scale integration (VLSI) chip design are just a few examples of tasks completed for this project. The BVD bit error rates (BER), measured from hardware and software simulations, are plotted as a function of bit signal to noise ratio E sub b/N sub 0 on the additive white Gaussian noise channel. Using the constraint length 15, rate 1/4, experimental convolutional code for the Galileo mission, the BVD gains 1.5 dB over the NASA standard (7,1/2) Maximum Likelihood Convolution Decoder (MCD) at a BER of 0.005. At this BER, the same gain results when the (255,233) NASA standard Reed-Solomon decoder is used, which yields a word error rate of 2.1 x 10(exp -8) and a BER of 1.4 x 10(exp -9). The (15, 1/6) code to be used by the Cometary Rendezvous Asteroid Flyby (CRAF)/Cassini Missions yields 1.7 dB of coding gain. These gains are measured with respect to symbols input to the BVD and increase with decreasing BER. Also, 8-bit input symbol quantization makes the BVD resistant to demodulated signal-level variations which may cause higher bandwidth than the NASA (7,1/2) code, these gains are offset by about 0.1 dB of expected additional receiver losses. Coding gains of several decibels are possible by compressing all spacecraft data.

Onyszchuk, I. M.↗
NASA NTRS1992

Application of piezoelectric devices to vibration suppression - From modeling and controller designs to implementation

Embedded piezoelectric devices may be ideally suited for vibration control of space structures, which lack an inertial ground. When subjected to an input voltage, an embedded piezoelectric actuator changes its dimensions, which in turn generates a pair of forces exerted on adjacent structural members. From the direct piezoelectric effect, an embedded piezoelectric transducer generates an electric charge proportional to the structural dynamic response. In this paper, the implementation, testing and modeling of an active truss structure consisting of piezoelectric sensors and actuators are described. Linear quadratic Gaussian (LQG), second-order, and direct rate feedback control schemes are designed to suppress the vibrations of the active structure. Simulation and test results are presented. It is shown that special model reduction considerations are required to achieve good correlation between test and analysis.

Won, Chin C.↗
NASA NTRS1993

Neural control of magnetic suspension systems

The purpose of this research program is to design, build and test (in cooperation with NASA personnel from the NASA Langley Research Center) neural controllers for two different small air-gap magnetic suspension systems. The general objective of the program is to study neural network architectures for the purpose of control in an experimental setting and to demonstrate the feasibility of the concept. The specific objectives of the research program are: (1) to demonstrate through simulation and experimentation the feasibility of using neural controllers to stabilize a nonlinear magnetic suspension system; (2) to investigate through simulation and experimentation the performance of neural controllers designs under various types of parametric and nonparametric uncertainty; (3) to investigate through simulation and experimentation various types of neural architectures for real-time control with respect to performance and complexity; and (4) to benchmark in an experimental setting the performance of neural controllers against other types of existing linear and nonlinear compensator designs. To date, the first one-dimensional, small air-gap magnetic suspension system has been built, tested and delivered to the NASA Langley Research Center. The device is currently being stabilized with a digital linear phase-lead controller. The neural controller hardware is under construction. Two different neural network paradigms are under consideration, one based on hidden layer feedforward networks trained via back propagation and one based on using Gaussian radial basis functions trained by analytical methods related to stability conditions. Some advanced nonlinear control algorithms using feedback linearization and sliding mode control are in simulation studies.

Gray, W. Steven↗
NASA NTRSDate not supplied

Structured Covariance Gaussian Networks for Orion Crew Module Aerodynamic Uncertainty Quantification

In this paper we propose a new approach for nonlinear regression and uncertainty quantification. The method is based on a pair of neural networks which parameterize mean and dense covariance functions of a multivariate Gaussian process, trained together to maximize the log-likelihood of observing the given data. The covariance matrix is made positive definite at every input by construction. We also propose a sampling approach that produces viable surrogate function realizations from the Gaussian process. We call the proposed model a Structured Covariance Gaussian Network (SCGN). We illustrate the use of SCGNs for learning an aerodynamic response surface with built-in uncertainty for the Orion crew module. We find that SCGN provides an efficient and systematic way to learn nonlinear functional relationships and dense covariances. We compare results to a baseline Gaussian process regressor and observe that the SCGN provides comparable uncertainty descriptions with improved scalability to dataset size. The sample functions generated by SCGN are fast to evaluate online and are therefore convenient for use in trajectory simulations. These results suggest that SCGN may be a viable computational method for aerodynamic uncertainty quantification.

machine learning↗
NASA NTRSDate not supplied

Structured Covariance Gaussian Networks for Orion Crew Module Aerodynamic Uncertainty Quantification

In this paper we propose a new approach for nonlinear regression and uncertainty quantification. The method is based on a pair of neural networks which parameterize mean and dense covariance functions of a multivariate Gaussian process, trained together to maximize the log-likelihood of observing the given data. The covariance matrix is made positive definite at every input by construction. We also propose a sampling approach that produces viable surrogate function realizations from the Gaussian process. We call the proposed model a Structured Covariance Gaussian Network (SCGN). We illustrate the use of SCGNs for learning an aerodynamic response surface with built-in uncertainty for the Orion crew module. We find that SCGN provides an efficient and systematic way to learn nonlinear functional relationships and dense covariances. We compare results to a baseline Gaussian process regressor and observe that the SCGN provides comparable uncertainty descriptions with improved scalability to dataset size. The sample functions generated by SCGN are fast to evaluate online and are therefore convenient for use in trajectory simulations. These results suggest that SCGN may be a viable computational method for aerodynamic uncertainty quantification.

machine learning↗
NASA NTRS1988

The topology of large-scale structure. II - Nonlinear evolution of Gaussian models

The evolution of non-Gaussian behavior in the large-scale universe from Gaussian initial conditions is studied. Topology measures developed in previous papers are applied to the smoothed initial, final, and biased matter distributions of cold dark matter, white noise, and massive neutrino simulations. When the smoothing length is approximately twice the mass correlation length or larger, the evolved models look like the initial conditions, suggesting that random phase hypotheses in cosmology can be tested with adequate data sets. When a smaller smoothing length is used, nonlinear effects are recovered, so nonlinear effects on topology can be detected in redshift surveys after smoothing at the mean intergalaxy separation. Hot dark matter models develop manifestly non-Gaussian behavior attributable to phase correlations, with a topology reminiscent of bubble or sheet distributions. Cold dark matter models remain Gaussian, and biasing does not disguise this.

Melott, Adrian L.↗
NASA NTRS1975

Performance of quasi-optimum digital FM demodulators for fading channels

This paper deals with the problem of digital demodulation of FM signals transmitted over Rayleigh and Rician fading channels. The Rayleigh and Rician fading channels are represented by two quadrature multiplicative nonzero mean white Gaussian processes in addition to an additive zero-mean white Gaussian noise. Quasi-optimum digital baseband demodulation algorithms using various nonlinear estimation techniques are derived. The digital demodulator structures are then simulated on a digital computer for an FM system with first order message spectrum for various values of the parameters for Rayleigh and Rician channels.

Dharamsi, M. T.↗
NASA NTRS1992

An investigation of error characteristics and coding performance

The performance of forward error correcting coding schemes on errors anticipated for the Earth Observation System (EOS) Ku-band downlink are studied. The EOS transmits picture frame data to the ground via the Telemetry Data Relay Satellite System (TDRSS) to a ground-based receiver at White Sands. Due to unintentional RF interference from other systems operating in the Ku band, the noise at the receiver is non-Gaussian which may result in non-random errors output by the demodulator. That is, the downlink channel cannot be modeled by a simple memoryless Gaussian-noise channel. From previous experience, it is believed that those errors are bursty. The research proceeded by developing a computer based simulation, called Communication Link Error ANalysis (CLEAN), to model the downlink errors, forward error correcting schemes, and interleavers used with TDRSS. To date, the bulk of CLEAN was written, documented, debugged, and verified. The procedures for utilizing CLEAN to investigate code performance were established and are discussed.

Ebel, William J.↗
NASA NTRS2005

Initial Kernel Timing Using a Simple PIM Performance Model

This presentation will describe some initial results of paper-and-pencil studies of 4 or 5 application kernels applied to a processor-in-memory (PIM) system roughly similar to the Cascade Lightweight Processor (LWP). The application kernels are: * Linked list traversal * Sun of leaf nodes on a tree * Bitonic sort * Vector sum * Gaussian elimination The intent of this work is to guide and validate work on the Cascade project in the areas of compilers, simulators, and languages. We will first discuss the generic PIM structure. Then, we will explain the concepts needed to program a parallel PIM system (locality, threads, parcels). Next, we will present a simple PIM performance model that will be used in the remainder of the presentation. For each kernel, we will then present a set of codes, including codes for a single PIM node, and codes for multiple PIM nodes that move data to threads and move threads to data. These codes are written at a fairly low level, between assembly and C, but much closer to C than to assembly. For each code, we will present some hand-drafted timing forecasts, based on the simple PIM performance model. Finally, we will conclude by discussing what we have learned from this work, including what programming styles seem to work best, from the point-of-view of both expressiveness and performance.

BRIEFING CHARTS↗
NASA NTRS1994

Comparison of dynamical approximation schemes for non-linear gravitational clustering

We have recently conducted a controlled comparison of a number of approximations for gravitational clustering against the same n-body simulations. These include ordinary linear perturbation theory (Eulerian), the adhesion approximation, the frozen-flow approximation, the Zel'dovich approximation (describable as first-order Lagrangian perturbation theory), and its second-order generalization. In the last two cases we also created new versions of approximation by truncation, i.e., smoothing the initial conditions by various smoothing window shapes and varying their sizes. The primary tool for comparing simulations to approximation schemes was crosscorrelation of the evolved mass density fields, testing the extent to which mass was moved to the right place. The Zel'dovich approximation, with initial convolution with a Gaussian e(exp -k(exp 2)/k(exp 2, sub G)) where k(sub G) is adjusted to be just into the nonlinear regime of the evolved model (details in text) worked extremely well. Its second-order generalization worked slightly better. All other schemes, including those proposed as generalizations of the Zel'dovich approximation created by adding forces, were in fact generally worse by this measure. By explicitly checking, we verified that the success of our best-choice was a result of the best treatment of the phases of nonlinear Fourier components. Of all schemes tested, the adhesion approximation produced the most accurate nonlinear power spectrum and density distribution, but its phase errors suggest mass condensations were moved to slightly the wrong location. Due to its better reproduction of the mass density distribution function and power spectrum, it might be preferred for some uses. We recommend either n-body simulations or our modified versions of the Zel'dovich approximation, depending upon the purpose. The theoretical implication is that pancaking is implicit in all cosmological gravitational clustering, at least from Gaussian initial conditions, even when subcondensations are present.

Melott, Adrian L.↗