Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Statistical 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 901 records · Page 50

The structure of the vorticity field in turbulent channel flow. Part 1: Analysis of instantaneous fields and statistical correlations

An investigation into the existence of hairpin vortices in turbulent channel flow is conducted using a database generated by the large eddy simulation technique. It is shown that away from the wall the distribution of the inclination angle of vorticity vector attains its maximum at about 45 deg to the wall. Two point correlations of velocity and vorticity fluctuations strongly support a flow model consisting of vortical structures inclined at 45 deg to the wall. The instantaneous vorticity vectors plotted in planes inclined at 45 deg show that the flow contains an appreciable number of hairpins. Vortex lines are used to display the three dimensional structure of hairpins, which are shown to be generated from deformation of transverse vortex filaments.

Moin, P.↗

The structure of the vorticity field in turbulent channel flow. I - Analysis of instantaneous fields and statistical correlations

An investigation into the existence of hairpin vortices in turbulent channel flow is conducted using a database generated by the large eddy simulation technique. It is shown that away from the wall the distribution of the inclination angle of vorticity vector attains its maximum at about 45 deg to the wall. Two point correlations of velocity and vorticity fluctuations strongly support a flow model consisting of vortical structures inclined at 45 deg to the wall. The instantaneous vorticity vectors plotted in planes inclined at 45 deg show that the flow contains an appreciable number of hairpins. Vortex lines are used to display the three dimensional structure of hairpins, which are shown to be generated from deformation of transverse vortex filaments.

Moin, P.↗

Simulations of sonic boom ray tube area fluctuations for propagation through atmospheric turbulence including caustics via a Monte Carlo method

A theory which gives statistical predictions for how often sonic booms propagating through the earth's turbulent boundary layer will encounter caustics, given the spectral properties of the atmospheric turbulence, is outlined. The theory is simple but approximately accounts for the variation of ray tube areas along ray paths. This theory predicts that the variation of ray tube areas is determined by the product of two similar area factors, psi (x) and phi (x), each satisfying a generic harmonic oscillator equation. If an area factor increases the peak acoustic pressure decreases, and if the factor decreases the peak acoustic pressure increases. Additionally, if an area factor decreases to zero and becomes negative, the ray has propagated through a caustic, which contributes a phase change of 90 degrees to the wave. Thus, it is clear that the number of times that a sonic boom wave passes through a caustic should be related to the distorted boom waveform received on the ground. Examples are given based on a characterization of atmospheric turbulence due to the structure function of Tatarski as modified by Crow.

Sparrow, Victor W.↗

Jet Noise Physics and Modeling Using First-principles Simulations

An extensive analysis of our jet DNS database has provided for the first time the complex correlations that are the core of many statistical jet noise models, including MGBK. We have also for the first time explicitly computed the noise from different components of a commonly used noise source as proposed in many modeling approaches. Key findings are: (1) While two-point (space and time) velocity statistics are well-fitted by decaying exponentials, even for our low-Reynolds-number jet, spatially integrated fourth-order space/retarded-time correlations, which constitute the noise "source" in MGBK, are instead well-fitted by Gaussians. The width of these Gaussians depends (by a factor of 2) on which components are considered. This is counter to current modeling practice, (2) A standard decomposition of the Lighthill source is shown by direct evaluation to be somewhat artificial since the noise from these nominally separate components is in fact highly correlated. We anticipate that the same will be the case for the Lilley source, and (3) The far-field sound is computed in a way that explicitly includes all quadrupole cancellations, yet evaluating the Lighthill integral for only a small part of the jet yields a far-field noise far louder than that from the whole jet due to missing nonquadrupole cancellations. Details of this study are discussed in a draft of a paper included as appendix A.

Freund, Jonathan B.↗

A time-dependent gyro-kinetic model of thermal ion upflows in the high-latitude F region

Ample evidence supports the significance of the high-latitude ionospheric contribution to magnetospheric plasma. Assuming flux conservation along a flux tube, the upward field-aligned ion flows observed in the magnetosphere require high-latitude ionospheric field-aligned ion upflows of the order of 10(exp 8) to 10(exp 9)/sq cm/s. Since radar and satellite observations of high-latitude F region flows at times exceed this flux requirement by an order of magnitude, the thermal ionospheric upflows are not simply the ionospheric response to a magnetospheric flux requirement. Several ionospheric ion upflow mechanisms have been proposed, but simulations based on fluid theory do not reproduce all the observed features of ionospheric ion upflows. Certain asymmetries in the statistical morphology of high-latitude F region ion upflows suggest that the ion upflows may be generated by ion-neutral frictional heating. We developed a single-component (O(+)), time-dependent gyro-kinetic model of the high-latitude F region response to frictional heating in which the neutral exobase is a discontinuous boundary between fully collisional and collisionless plasmas. The concept of a discontinuous neutreal exobase and the assumption of a constant and uniform polarization electric field reduce the ion velocity distribution function, from which we can compute the ion density, parallel velocity, parallel and perpendicular temperature, and parallel flux. Using our model, we simulated the response of a convecting flux tube between 500 km and 2500 km to various frictional heating inputs; the results were both qualitatively and quantitatively different from fluid model results, which may indicate an inadequacy of the fluid theory approach. The gyro-kinetic frictional heating model responses to the various simulations were qualitatively similar: (1) initial perturbations of all the modeled parameters propagated rapidly up the flux tube, (2) transient values of the ion parallel velocity, temperature, and flux exceeded 3 km/s, 2 x 10(exp 4) K, and 10(exp 9)/sq cm/s, respectively, (3) a second transient regime developed wherein the parallel temperature drops to very low values (a few hundred Kelvins), and (4) well after heating ceased, large parallel temperatures and large downward parallel velocities and fluxes developed as the flux tube slowly returned to diffusive equilibrium. The ion velocity distributions during the simulation are often non-Maxwellian and are sometimes composed of two distinct ion populations.

Loranc, M.↗

Abundance recovery error analysis using simulated AVIRIS data

Measurement noise and imperfect atmospheric correction translate directly into errors in the determination of the surficial abundance of materials from imaging spectrometer data. The effects of errors on abundance recovery were investigated previously using Monte Carlo simulation methods by Sabol et. al. The drawback of the Monte Carlo approach is that thousands of trials are needed to develop good statistics on the probable error in abundance recovery. This computational burden invariably limits the number of scenarios of interest that can practically be investigated. A more efficient approach is based on covariance analysis. The covariance analysis approach expresses errors in abundance as a function of noise in the spectral measurements and provides a closed form result eliminating the need for multiple trials. Monte Carlo simulation and covariance analysis are used to predict confidence limits for abundance recovery for a scenario which is modeled as being derived from Airborne Visible/Infrared Imaging Spectrometer (AVIRIS).

Stoner, William W.↗

Transient structure in the high-energy X-ray light curve of NP 0532

The paper reports the observation of pulsed fractions in the primary and secondary peaks, as well as in the interpulse region, of the high-energy X-ray light curve of NP 0532. A statistical analysis of light-curve data is performed, and a similar analysis is carried out using simulated data. It is concluded that a previously reported third peak in the light curve was transient in nature.

Ryckman, S. G.↗

A study of digital holographic filters generation. Phase 2: Digital data communication system, volume 1

An empirical study of the performance of the Viterbi decoders in bursty channels was carried out and an improved algebraic decoder for nonsystematic codes was developed. The hybrid algorithm was simulated for the (2,1), k = 7 code on a computer using 20 channels having various error statistics, ranging from pure random error to pure bursty channels. The hybrid system outperformed both the algebraic and the Viterbi decoders in every case, except the 1% random error channel where the Viterbi decoder had one bit less decoding error.

Ingels, F. M.↗

Parameter estimation in X-ray astronomy using maximum likelihood

Methods of estimation of parameter values and confidence regions by maximum likelihood and Fisher efficient scores starting from Poisson probabilities are developed for the nonlinear spectral functions commonly encountered in X-ray astronomy. It is argued that these methods offer significant advantages over the commonly used alternatives called minimum chi-squared because they rely on less pervasive statistical approximations and so may be expected to remain valid for data of poorer quality. Extensive numerical simulations of the maximum likelihood method are reported which verify that the best-fit parameter value and confidence region calculations are correct over a wide range of input spectra.

Wachter, K.↗

Remote sensing of atmospheric water content from Bhaskara SAMIR data

The 19.35 GHz and 22.235 GHz passive microwave radiometers (SAMIR) on board the Indian satellite Bhaskara have provided very useful data. From these data has been demonstrated the feasibility of deriving atmospheric and ocean surface parameters such as water vapor content, liquid water content, rainfall rate and ocean surface winds. Different approaches have been tried for deriving the atmospheric water content. The statistical and empirical methods have been used by others for the analysis of the Nimbus data. A simulation technique has been attempted for the first time for 19.35 GHz and 22.235 GHz radiometer data. The results obtained from three different methods are compared with radiosonde data. A case study of a tropical depression has been undertaken to demonstrate the capability of Bhaskara SAMIR data to show the variation of total water vapor and liquid water contents.

Gohil, B. S.↗

Results of the joint NASA/NOAA AMTS/HIRS2 sounding simulation test

The relative accuracies of atmospheric temperature profiles retrieved from HIRS2, the current operational infra-red temperature sounder, and AMTS, a proposed advanced high-spectral resolution infra-red sounder were compared. The sounding simulation test compared retrievals generated by GLAS, using their physical retrieval algorithm, and NESDIS, using their operational statistical regression algorithm, for both instruments under clear and cloudy conditions. Test results for both conditions are summarized in tables showing RMS temperature errors for 18 tropospheric layers, 4 stratospheric layers, and surface skin temperature. The improvement of the AMTS instrument over HIRS2, particularly for cloudy conditions is demonstrated. The GLAS retrievals are more accurate than the operational NESDIS retrieval, especially with increasing cloud cover.

Susskind, J.↗

Computer analysis of X-band radar data

The effectiveness of using currently available computer techniques for interpretation of MSS data to interpret SAR imagery for forest monitoring was assessed. Data were gathered with NASA's airborne APQ-102 dual-polarized, X-band SAR in a flight at 60,000 ft. Microdensitometry was employed to digitize the HH- and HV-polarized imagery. A ground spatial resolution of 15 m was obtained, control points were identified, a second order biquadratic transformation was applied to compensate for orientation, and rms errors were calculated. A second data set was taken with 30 m resolution in order to simulate thematic mapper operation. Classification was performed with pixel-by-pixel and textural classification algorithms. A statistical analysis was also carried out to find any significant differences between classifiers in a data set for a given classifier. Each polarization featured an independent distortion which required appropriate preprocessing to correct. Further studies are recommended with multiple frequencies viewing and multiple polarizations and look angles to define the actual forest classifications that can be made with the SAR imagery.

Knowlton, D. J.↗

Composite load spectra for select space propulsion structural components

The objective of this program is to develop generic load models to simulate the composite load spectra (CLS) that are induced in space propulsion system components representative of the space shuttle main engines (SSME). These models are being developed through describing individual component loads with an appropriate mix of deterministic and state-of-the-art probabilistic models that are related to key generic variables. Combinations of the individual loads are used to synthesize the composite loads spectra. A second approach for developing the composite loads spectra load model simulation, the option portion of the contract will develop coupled models which combine the individual load models. Statistically varying coefficients of the physical models will be used to obtain the composite load spectra.

Newell, J. F.↗

The effect of exercise on venous gas emboli and decompression sickness in human subjects at 4.3 psia

The contribution of upper body exercise to altitude decompression sickness while at 4.3 psia after 3.5 or 4.0 hours of 100% oxygen prebreathing at 14.7 psia was determined by comparing the incidence and patterns of venous gas emboli (VGE), and the incidence of Type 1 decompression sickness (DCS) in 43 exercising male subjects and 9 less active male Doppler Technicians (DT's). Each subject exercised for 4 minutes at each of 3 exercise stations while at 4.3 psia. An additional 4 minutes were spent monitoring for VGE by the DT while the subject was supine on an examination cot. In the combined 3.5 and 4.0 hour oxygen prebreathe data, 13 subjects complained of Type 1 DCS compared to 9 complaints from DT's. VGE were detected in 28 subjects compared to 14 detections from DT's. A chi-square analysis of proportions showed no statistically significantly difference in the incidence of Type 1 DCS or VGE between the two groups; however, the average time to detect VGE and to report Tyep 1 DCS symptoms were statistically different. It was concluded that 4 to 6 hours of upper body exercise at metabolic rates simulating EVA metabolic rates hastens the initial detection of VGE and the time to report Type 1 DCS symptoms as compared to DT's.

Conkin, Johnny↗

Guidelines for the use of programmable display pushbuttons on the Space Station's telerobot control panel

Simulations of a Remote Manipulator System task on microcomputer prototypes were performed to establish guidelines for using programmable display pushbuttons (PDP) on a telerobot control panel. Simulations of control panels with and without PDP were conducted. It is shown that there is a statistically significant increase in the number of commands issued in the non-PDP control panel. Guidelines for using PDP panels are listed.

Stuart, Mark A.↗

On the prediction of the behavior of the velocity, temperature and humidity fields above a thermally stratified air-water interface

A turbulent kinetic energy-dissipation (k-epsilon) model is employed to numerically simulate the behavior of the mean velocity, temperature and humidity fields, along with such second-order statistics as Reynolds stresses, turbulent heat and mass transfer correlation, and bulk transfer coefficients in the boundary layer above wind-generated laboratory waves. The k and epsilon equations have been properly modified to account for thermal stratification, the presence of free surface, and favorable pressure gradient effects. Predictions are compared with available data.

Papadimitrakis, Yiannis A.↗

Transport Of Passive Scalars In A Turbulent Channel Flow

Computer simulation of transport of passive scalars in turbulent channel flow described in report. Shows flow structures and statistical properties. As used here, "passive scalars" means scalar quantities like fluctuations in temperature or concentrations of contaminants that do not disturb flow appreciably. Examples include transport of heat in heat exchangers, gas turbines, and nuclear reactors and dispersal of pollution in atmosphere.

Kim, John↗

Systematic methods for knowledge acquisition and expert system development

Nine cooperating rule-based systems, collectively called AUTOCREW which were designed to automate functions and decisions associated with a combat aircraft's subsystems, are discussed. The organization of tasks within each system is described; performance metrics were developed to evaluate the workload of each rule base and to assess the cooperation between the rule bases. Simulation and comparative workload results for two mission scenarios are given. The scenarios are inbound surface-to-air-missile attack on the aircraft and pilot incapacitation. The methodology used to develop the AUTOCREW knowledge bases is summarized. Issues involved in designing the navigation sensor selection expert in AUTOCREW's NAVIGATOR knowledge base are discussed in detail. The performance of seven navigation systems aiding a medium-accuracy INS was investigated using Kalman filter covariance analyses. A navigation sensor management (NSM) expert system was formulated from covariance simulation data using the analysis of variance (ANOVA) method and the ID3 algorithm. ANOVA results show that statistically different position accuracies are obtained when different navaids are used, the number of navaids aiding the INS is varied, the aircraft's trajectory is varied, and the performance history is varied. The ID3 algorithm determines the NSM expert's classification rules in the form of decision trees. The performance of these decision trees was assessed on two arbitrary trajectories, and the results demonstrate that the NSM expert adapts to new situations and provides reasonable estimates of the expected hybrid performance.

Belkin, Brenda L.↗