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 361 records · Page 20

Sensitivity Analysis and Uncertainty Quantification of a Mars Ascent Vehicle Concept

The design of a conceptual Mars ascent vehicle is a challenging problem. In order to aid the vehicle and mission concept design it is important to understand the driving design parameters and the expected performance in the presence of model errors and uncertainties. An existing six degree of freedom simulation model is analyzed on a statistical basis using the methods available in the Design Analysis Kit for Optimization and Terascale Applications toolkit. The methods utilized include conventional Monte Carlo techniques, metamodeling via polynomial chaos expansions, and global variance-based sensitivity analyses. Two additional analysis methods referred to as “Monte Carlo filtering” and ”Classification trees” are used to determine which uncertain parameters are driving the performance of the vehicle. Monte Carlo filtering provides a methodology to determine which parameters cause qualitatively different behavior while the classification trees use heuristics to partition the input space and assign probabilities to each partition. These methods serve as qualitative descriptors of model sensitivity while variance-based global sensitivity analysis seeks a quantitative mapping from total output variance to the variance of individual inputs. Application of these techniques to several outputs of a Mars ascent vehicle concept simulation indicates that only a select few input factors dominate their variance.

Noyes, Connor↗

Nonstationary atmospheric boundary layer turbulence simulation

Report on a new and general technique for simulating atmospheric turbulence-like random processes which are statistically homogeneous along the horizontal and nonhomogeneous along the vertical. This technique is general in the sense that it can be used for a broad class of similar problems. Like the other presently available schemes, the techniques presented are based on the Dryden hypothesis and Taylor's frozen eddy hypothesis; however, they go a step further by utilizing certain self-similarity properties of the Dryden spectral density function which permits the development of height invariant filters. These filters are in turn used to generate vertically homogeneous (statistically) random processes from which turbulence at any specified level in the boundary layer can be simulated, thus facilitating the simulation of a nonstationary turbulence process along the flight path of an aircraft during take-off or landing.

Fichtl, G. H.↗

On the error statistics of Viterbi decoding and the performance of concatenated codes

Computer simulation results are presented on the performance of convolutional codes of constraint lengths 7 and 10 concatenated with the (255, 223) Reed-Solomon code (a proposed NASA standard). These results indicate that as much as 0.8 dB can be gained by concatenating this Reed-Solomon code with a (10, 1/3) convolutional code, instead of the (7, 1/2) code currently used by the DSN. A mathematical model of Viterbi decoder burst-error statistics is developed and is validated through additional computer simulations.

Miller, R. L.↗

Digital simulation of a communication link for Pioneer Saturn Uranus atmospheric entry probe, part 1

A digital simulation study is presented for a candidate modulator/demodulator design in an atmospheric scintillation environment with Doppler, Doppler rate, and signal attenuation typical of the conditions of an outer planet atmospheric probe. The simulation results indicate that the mean channel error rate with and without scintillation are similar to theoretical characterizations of the link. The simulation gives information for calculating other channel statistics and generates a quantized symbol stream on magnetic tape from which error correction decoding is analyzed. Results from the magnetic tape data analyses are also included. The receiver and bit synchronizer are modeled in the simulation at the level of hardware component parameters rather than at the loop equation level and individual hardware parameters are identified. The atmospheric scintillation amplitude and phase are modeled independently. Normal and log normal amplitude processes are studied. In each case the scintillations are low pass filtered. The receiver performance is given for a range of signal to noise ratios with and without the effects of scintillation. The performance is reviewed for critical reciever parameter variations.

Hinrichs, C. A.↗

POLARRIS: A POLArimetric Radar Retrieval and Instrument Simulator

This paper introduces a synthetic polarimetric radar simulator and retrieval package, POLArimetric Radar Retrieval and Instrument Simulator (POLARRIS), for evaluating cloud‐resolving models (CRMs). POLARRIS is composed of forward (POLARRIS‐f) and inverse (retrieval and diagnostic) components (iPOLARRIS) to generate not only polarimetric radar observables (Zh, Zdr, Kdp, ρhv) but also radar‐consistent geophysical parameters such as hydrometeor identification, vertical velocity, and rainfall rates retrieved from CRM data. To demonstrate its application and uncertainties, POLARRIS is applied to simulations of a mesoscale convective system over the Southern Great Plains on 23 May 2011, using the Weather Research and Forecasting model with both spectral bin microphysics (SBM) and the Goddard single‐moment bulk 4ICE microphysics. Statistical composites reveal a significant dependence of simulated polarimetric observables (Zdr, Kdp) on the assumptions of the particle axis ratio (oblateness) and orientation angle distributions. The simulated polarimetric variables differ considerably between the SBM and 4ICE microphysics in part due to the differences in their ice particle size distributions as revealed by comparisons with aircraft measurements. Regardless of these uncertainties, simulated hydrometeor identification distributions overestimate graupel and hail fractions, especially from the simulation with SBM. To minimize uncertainties in forward model, the particle shape and orientation angle distributions of frozen particles should be predicted in a microphysics scheme in addition to the size distributions and particle densities.

Toshi Matsui↗

Sampling errors for satellite-derived tropical rainfall - Monte Carlo study using a space-time stochastic model

Estimates of monthly average rainfall based on satellite observations from a low earth orbit will differ from the true monthly average because the satellite observes a given area only intermittently. This sampling error inherent in satellite monitoring of rainfall would occur even if the satellite instruments could measure rainfall perfectly. The size of this error is estimated for a satellite system being studied at NASA, the Tropical Rainfall Measuring Mission (TRMM). First, the statistical description of rainfall on scales from 1 to 1000 km is examined in detail, based on rainfall data from the Global Atmospheric Research Project Atlantic Tropical Experiment (GATE). A TRMM-like satellite is flown over a two-dimensional time-evolving simulation of rainfall using a stochastic model with statistics tuned to agree with GATE statistics. The distribution of sampling errors found from many months of simulated observations is found to be nearly normal, even though the distribution of area-averaged rainfall is far from normal. For a range of orbits likely to be employed in TRMM, sampling error is found to be less than 10 percent of the mean for rainfall averaged over a 500 x 500 sq km area.

Bell, Thomas L.↗

Use of High-resolution WRF Simulations to Forecast Lightning Threat

Recent observational studies have confirmed the existence of a robust statistical relationship between lightning flash rates and the amount of large precipitating ice hydrometeors in storms. This relationship is exploited, in conjunction with the capabilities of recent forecast models such as WRF, to forecast the threat of lightning from convective storms using the output fields from the model forecasts. The simulated vertical flux of graupel at -15C is used in this study as a proxy for charge separation processes and their associated lightning risk. Six-h simulations are conducted for a number of case studies for which three-dimensional lightning validation data from the North Alabama Lightning Mapping Array are available. Experiments indicate that initialization of the WRF model on a 2 km grid using Eta boundary conditions, Doppler radar radial velocity and reflectivity fields, and METAR and ACARS data yield the most realistic simulations. An array of subjective and objective statistical metrics are employed to document the utility of the WRF forecasts. The simulation results are also compared to other more traditional means of forecasting convective storms, such as those based on inspection of the convective available potential energy field.

McCaul, William E.↗

Digital simulation of hybrid loop operation in RFI backgrounds.

A digital computer model for Monte-Carlo simulation of an imperfect second-order hybrid phase-locked loop (PLL) operating in radio-frequency interference (RFI) and Gaussian noise backgrounds has been developed. Characterization of hybrid loop performance in terms of cycle slipping statistics and phase error variance, through computer simulation, indicates that the hybrid loop has desirable performance characteristics in RFI backgrounds over the conventional PLL or the costas loop.

Ziemer, R. E.↗

Three Dimensional CFD Analysis of the GTX Combustor

The annular combustor geometry of a combined-cycle engine has been analyzed with three-dimensional computational fluid dynamics. Both subsonic combustion and supersonic combustion flowfields have been simulated. The subsonic combustion analysis was executed in conjunction with a direct-connect test rig. Two cold-flow and one hot-flow results are presented. The simulations compare favorably with the test data for the two cold flow calculations; the hot-flow data was not yet available. The hot-flow simulation indicates that the conventional ejector-ramjet cycle would not provide adequate mixing at the conditions tested. The supersonic combustion ramjet flowfield was simulated with frozen chemistry model. A five-parameter test matrix was specified, according to statistical design-of-experiments theory. Twenty-seven separate simulations were used to assemble surrogate models for combustor mixing efficiency and total pressure recovery. ScramJet injector design parameters (injector angle, location, and fuel split) as well as mission variables (total fuel massflow and freestream Mach number) were included in the analysis. A promising injector design has been identified that provides good mixing characteristics with low total pressure losses. The surrogate models can be used to develop performance maps of different injector designs. Several complex three-way variable interactions appear within the dataset that are not adequately resolved with the current statistical analysis.

Steffen, C. J., Jr.↗

Three Dimensional CFD Analysis of the GTX Combustor

The annular combustor geometry of a combined-cycle engine has been analyzed with three-dimensional computational fluid dynamics. Both subsonic combustion and supersonic combustion flowfields have been simulated. The subsonic combustion analysis was executed in conjunction with a direct-connect test rig. Two cold-flow and one hot-flow results are presented. The simulations compare favorably with the test data for the two cold flow calculations; the hot-flow data was not yet available. The hot-flow simulation Indicates that the conventional ejector-ramjet cycle would not provide adequate mixing at the conditions tested. The supersonic combustion ramjet flowfield was simulated with frozen chemistry model. A five-parameter test matrix was specified, according to statistical design-of-experiments theory. Twenty-seven separate simulations were used to assemble surrogate models for combustor mixing efficiency and total pressure recovery. Scramjet injector design parameters (injector angle, location, and fuel split) as well as mission variables (total fuel mass flow and freestream Mach number) were included in the analysis. A promising injector design has been identified that provides good mixing characteristics with low total pressure losses. The surrogate models can be used to develop performance maps of different injector designs. Several complex three-way variable interactions appear within the dataset that are not adequately resolved with the current statistical analysis.

Steffen, C. J., Jr.↗

The Latitude Dependence of the Effect of Pinatubo on Stratospheric Ozone

Statistical analysis of TOMS and SBLT total ozone data indicate that the eruption of Pinatubo in 1991 led to a significant decrease in ozone at northern midlatitudes with little or no effect at southern midlatitudes. We argue that this puzzling absence of a southern hemisphere effect may be an artifact of the statistical analysis. We have run a 3D CTM simulation of the past 30 years of stratospheric photochemistry with variable forcing due to chlorine/bromine compounds, solar ultraviolet radiation, and volcanic aerosols. This integration used winds from the FVGCM, which has similar interannual variability to the atmosphere. When this CTM output was examined with a standard time-series analysis, we found an effect of Pinatubo in the southern hemisphere, but not in the northern hemisphere. We then reran the CTM without volcanic aerosols. The subtraction of the two simulations indicated that, as expected, that Pinatubo affected both hemispheres in the model. This means that the northern hemisphere effect was in the model but did not show up in the statistical analysis. We also had an on-line parameterized chemical ozone tracer with seasonally repeating production and loss over the simulation. We used this as a dynamical surrogate to remove interannual variability from the original model output. The residual time series was then analyzed for the Pinatubo effect and we were able to find it in both hemispheres. We suggest that the combination of the two volcanoes, El Chichon and Pinatubo, with the solar cycle and interannual variability led to this problem of analysis in the northern hemisphere of our model. We furthermore suggest that a similar think may be occurring in the southern hemisphere of the data. An analysis of the atmosphere's southern hemisphere with a good dynamical surrogate may solve the mystery of the missing southern hemisphere effect of Pinatubo on ozone.

Stolarski, Richard S.↗

A Model Assessment of Satellite Observed Trends in Polar Sea Ice Extents

For more than three decades now, satellite passive microwave observations have been used to monitor polar sea ice. Here we utilize sea ice extent trends determined from primarily satellite data for both the Northern and Southern Hemispheres for the period 1972(73)-2004 and compare them with results from simulations by eleven climate models. In the Northern Hemisphere, observations show a statistically significant decrease of sea ice extent and an acceleration of sea ice retreat during the past three decades. However, from the modeled natural variability of sea ice extents in control simulations, we conclude that the acceleration is not statistically significant and should not be extrapolated into the future. Observations and model simulations show that the time scale of climate variability in sea ice extent in the Southern Hemisphere is much larger than in the Northern Hemisphere and that the Southern Hemisphere sea ice extent trends are not statistically significant.

Vinnikov, Konstantin Y.↗

Mission Analysis Program for Solar Electric Propulsion (MAPSEP). Volume 2: User's manual for earth orbital MAPSEP

The trajectory simulation mode (SIMSEP) requires the namelist SIMSEP to follow TRAJ. The SIMSEP contains parameters which describe the scope of the simulation, expected dynamic errors, and cumulative statistics from previous SIMSEP runs. Following SIMSEP are a set of GUID namelists, one for each guidance correction maneuver. The GUID describes the strategy, knowledge or estimation uncertainties and cumulative statistics for that particular maneuver. The trajectory display mode (REFSEP) requires only the namelist TRAJ followed by scheduling cards, similar to those used in GODSEP. The fixed field schedule cards define: types of data displayed, span of interest, and frequency of printout. For those users who can vary the amount of blank common storage in their runs, a guideline to estimate the total MAPSEP core requirements is given. Blank common length is related directly to the dimension of the dynamic state (NDIM) used in transition matrix (STM) computation, and, the total augmented (knowledge) state (NAUG). The values of program and blank common must be added to compute the total decimal core for a CDC 6500. Other operating systems must scale these requirements appropriately.

Source record↗

Subgrid or Reynolds stress-modeling for three-dimensional turbulence computations

A review is given of recent advances in two distinct computational methods for evaluating turbulence fields, namely, statistical Reynolds stress modeling and turbulence simulation, where large eddies are followed in time. It is shown that evaluation of the mean Reynolds stresses, rather than use of a scalar eddy viscosity, permits an explanation of streamline curvature effects found in several experiments. Turbulence simulation, with a new volume averaging technique and third-order accurate finite-difference computing is shown to predict the decay of isotropic turbulence in incompressible flow with rather modest computer storage requirements, even at Reynolds numbers of aerodynamic interest.

Rubesin, M. W.↗

Satellite disintegration dynamics

The subject of satellite disintegration is examined in detail. Elements of the orbits of individual fragments, determined by DOD space surveillance systems, are used to accurately predict the time and place of fragmentation. Dual time independent and time dependent analyses are performed for simulated and real breakups. Methods of statistical mechanics are used to study the evolution of the fragment clouds. The fragments are treated as an ensemble of non-interacting particles. A solution of Liouville's equation is obtained which enables the spatial density to be calculated as a function of position, time and initial velocity distribution.

Dasenbrock, R. R.↗

Compressor seal rub energetics study

The rub mechanics of compressor abradable blade tip seals at simulated engine conditions were investigated. Twelve statistically planned, instrumented rub tests were conducted with titanium blades and Feltmetal fibermetal rubstrips. The tests were conducted with single stationary blades rubbing against seal material bonded to rotating test disks. The instantaneous rub torque, speed, incursion rate and blade temperatures were continuously measured and recorded. Basic rub parameters (incursion rate, rub depth, abradable density, blade thickness and rub velocity) were varied to determine the effects on rub energy and heat split between the blade, rubstrip surface and rub debris. The test data was reduced, energies were determined and statistical analyses were completed to determine the primary and interactive effects. Wear surface morphology, profile measurements and metallographic analysis were used to determine wear, glazing, melting and material transfer. The rub energies for these tests were most significantly affected by the incursion rate while rub velocity and blade thickness were of secondary importance. The ratios of blade wear to seal wear were representative of those experienced in engine operation of these seal system materials.

Laverty, W. F.↗

Geological terrain models

The initial phase of a program to determine the best interpretation strategy and sensor configuration for a radar remote sensing system for geologic applications is discussed. In this phase, terrain modeling and radar image simulation were used to perform parametric sensitivity studies. A relatively simple computer-generated terrain model is presented, and the data base, backscatter file, and transfer function for digital image simulation are described. Sets of images are presented that simulate the results obtained with an X-band radar from an altitude of 800 km and at three different terrain-illumination angles. The simulations include power maps, slant-range images, ground-range images, and ground-range images with statistical noise incorporated. It is concluded that digital image simulation and computer modeling provide cost-effective methods for evaluating terrain variations and sensor parameter changes, for predicting results, and for defining optimum sensor parameters.

Kaupp, V. H.↗