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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 541 records · Page 30

All sources NASA NTRS DOE OSTI

NASA NTRS1972

A Management Information System Model for Program Management

The development of a model to simulate the information system of a program management type of organization is reported. The model statistically determines the following parameters: type of messages, destinations, delivery durations, type processing, processing durations, communication channels, outgoing messages, and priorites. The total management information system of the program management organization is considered, including formal and informal information flows and both facilities and equipment. The model is written in General Purpose System Simulation 2 computer programming language for use on the Univac 1108, Executive 8 computer. The model is simulated on a daily basis and collects queue and resource utilization statistics for each decision point. The statistics are then used by management to evaluate proposed resource allocations, to evaluate proposed changes to the system, and to identify potential problem areas. The model employs both empirical and theoretical distributions which are adjusted to simulate the information flow being studied.

Shipman, D. L.↗

NASA NTRS2000

Computational Analysis for Rocket-Based Combined-Cycle Systems During Rocket-Only Operation

A series of Reynolds-averaged Navier-Stokes calculations were employed to study the performance of rocket-based combined-cycle systems operating in an all-rocket mode. This parametric series of calculations were executed within a statistical framework, commonly known as design of experiments. The parametric design space included four geometric and two flowfield variables set at three levels each, for a total of 729 possible combinations. A D-optimal design strategy was selected. It required that only 36 separate computational fluid dynamics (CFD) solutions be performed to develop a full response surface model, which quantified the linear, bilinear, and curvilinear effects of the six experimental variables. The axisymmetric, Reynolds-averaged Navier-Stokes simulations were executed with the NPARC v3.0 code. The response used in the statistical analysis was created from Isp efficiency data integrated from the 36 CFD simulations. The influence of turbulence modeling was analyzed by using both one- and two-equation models. Careful attention was also given to quantify the influence of mesh dependence, iterative convergence, and artificial viscosity upon the resulting statistical model. Thirteen statistically significant effects were observed to have an influence on rocket-based combined-cycle nozzle performance. It was apparent that the free-expansion process, directly downstream of the rocket nozzle, can influence the Isp efficiency. Numerical schlieren images and particle traces have been used to further understand the physical phenomena behind several of the statistically significant results.

Steffen, C. J., Jr.↗

NASA NTRS2008

Mars Phoenix Entry, Descent, and Landing Simulation Design and Modelling Analysis

The 2007 Mars Phoenix Lander was launched in August of 2007 on a ten month cruise to reach the northern plains of Mars in May 2008. Its mission continues NASA s pursuit to find evidence of water on Mars. Phoenix carries upon it a slew of science instruments to study soil and ice samples from the northern region of the planet, an area previously undiscovered by robotic landers. In order for these science instruments to be useful, it was necessary for Phoenix to perform a safe entry, descent, and landing (EDL) onto the surface of Mars. The EDL design was defined through simulation and analysis of the various phases of the descent. An overview of the simulation and various models developed to characterize the EDL performance is provided. Monte Carlo statistical analysis was performed to assess the performance and robustness of the Phoenix EDL system and are presented in this paper. Using these simulation and modelling tools throughout the design and into the operations phase, the Mars Phoenix EDL was a success on May 25, 2008.

Prince, Jill L.↗

NASA NTRS2011

CUDA Simulation of Homogeneous, Incompressible Turbulence

We discuss very fast Compute Unified Device Architecture (CUDA) simulations of ideal homogeneous incompressible turbulence based on Fourier models. These models have associated statistical theories that predict that Fourier coefficients of fluid velocity and magnetic fields (if present) are zero-mean random variables. Prior numerical simulations have shown that certain coefficients have a non-zero mean value that can be very large compared to the associated standard deviation. We review the theoretical basis of this "broken ergodicity" as applied to 2-D and 3-D fluid and magnetohydrodynamic simulations of homogeneous turbulence. Our new simulations examine the phenomenon of broken ergodicity through very long time and large grid size runs performed on a state-of-the-art CUDA platform. Results comparing various CUDA hardware configurations and grid sizes are discussed. NS and MHD results are compared.

Morin, Lee↗

NASA NTRS1995

Development of a Probabilistic Component Mode Synthesis Method for the Analysis of Non-Deterministic Substructures

Standard methods of structural dynamic analysis assume that the structural characteristics are deterministic. Recognizing that these characteristics are actually statistical in nature, researchers have recently developed a variety of methods that use this information to determine probabilities of a desired response characteristic, such as natural frequency, without using expensive Monte Carlo simulations. One of the problems in these methods is correctly identifying the statistical properties of primitive variables such as geometry, stiffness, and mass. This paper presents a method where the measured dynamic properties of substructures are used instead as the random variables. The residual flexibility method of component mode synthesis is combined with the probabilistic methods to determine the cumulative distribution function of the system eigenvalues. A simple cantilever beam test problem is presented that illustrates the theory.

Brown, Andrew M.↗

NASA NTRS1996

Probabilistic Component Mode Synthesis of Nondeterministic Substructures

Standard methods of structural dynamic analysis assume that the structural characteristics are deterministic. Recognizing that these characteristics are actually statistical in nature researchers have recently developed a variety of methods that use this information to determine probabilities of a desired response characteristic, such as natural frequency, without using expensive Monte Carlo simulations. One of the problems in these methods is correctly identifying the statistical properties of primitive variables such as geometry, stiffness, and mass. We present a method where the measured dynamic properties of substructures are used instead as the random variables. The residual flexibility method of component mode synthesis is combined with the probabilistic methods to determine the cumulative distribution function of the system eigenvalues. A simple cantilever beam test problem is presented that illustrates the theory.

Brown, Andrew M.↗

NASA NTRS1995

Comparison of 2- and 10-micron coherent Doppler lidar performance

The performance of 2- and 10-micron coherent Doppler lidar is presented in terms of the statistical distribution of the maximum-likelihood velocity estimator from simulations for fixed range resolution and fixed velocity search space as a function of the number of coherent photoelectrons per estimate. The wavelength dependence of the aerosol backscatter coefficient, the detector quantum efficiency, and the atmospheric extinction produce a simple shift of the performance curves. Results are presented for a typical boundary layer measurement and a space-based measurement for two regimes: the pulse-dominated regime where the signal statistics are determined by the transmitted pulse, and the atmospheric-dominated regime where the signal statistics are determined by the velocity fluctuations over the range gate. The optimal choice of wavelength depends on the problem under consideration.

Frehlich, Rod↗

NASA NTRS1993

The compositional structure and the effects of exothermicity in a nonpremixed planar jet flame

Results are presented of direct numerical simulation (DNS) of a randomly perturbed compressible, spatially developing two-dimensional planar jet under the influence of a finite rate chemical reaction of the type F + O yields product. The objectives of the simulations are to assess the compositional structure of the flame and to determine the influence of reaction exothermicity by means of statistical sampling of the DNS generated data. It is shown that even with this idealized kinetics model the simulated results exhibit features in accord with experimental data. These results indicate that the Damkohler number is an important parameter in determining the statistical composition of the reacting field and that the results are not very sensitive to the mechanism by which this parameter is varied. It is demonstrated that as the intensity of mixing is increased and the effect of finite rate chemistry is more pronounced, the magnitudes of the ensemble mean and variance of the product mass fraction decrease and those of the reactants' mass fraction increase. Also, at higher mixing rates the joint probability density functions of the reactants' mass fractions shift towards higher values within the composition domain, indicating a lower reactedness.

Steinberger, C. J.↗

NASA NTRS2008

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 aloft in storms. This relationship is exploited, in conjunction with the capabilities of cloud-resolving forecast models such as WRF, to forecast explicitly the threat of lightning from convective storms using selected output fields from the model forecasts. The simulated vertical flux of graupel at -15C and the shape of the simulated reflectivity profile are tested in this study as proxies for charge separation processes and their associated lightning risk. Our lightning forecast method differs from others in that it is entirely based on high-resolution simulation output, without reliance on any climatological data. short [6-8 h) simulations are conducted for a number of case studies for which three-dmmensional 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 fields, and METAR and ACARS data y&eld satisfactory simulations. __nalyses of the lightning threat fields suggests that both the graupel flux and reflectivity profile approaches, when properly calibrated, can yield reasonable lightning threat forecasts, although an ensemble approach is probably desirable in order to reduce the tendency for misplacement of modeled storms to hurt the accuracy of the forecasts. Our lightning threat forecasts are also compared to other more traditional means of forecasting thunderstorms, such as those based on inspection of the convective available potential energy field.

McCaul, E. W., Jr.↗

NASA NTRS2005

Formation of a Tropopause Cirrus Layer Observed over Florida during CRYSTAL-FACE

On July 13, 2002 a widespread, subvisible tropopause cirrus layer occurred over the Florida region. This cloud was observed in great detail with the NASA Cirrus Regional Study of Tropical Anvils and Cirrus Layers-Florida Area Cirrus Experiment (CRYSTAL-FACE) instrumentation, including in situ measurements with the WB-57 aircraft. In this paper, we use the 13 July cloud as a case study to evaluate the physical processes controlling the formation and evolution of tropopause cirrus layers. Microphysics measurements indicate that ice crystal diameters in the cloud layer ranged from about 7 to 50 microns, and the peak number mode was about 10-25 microns. In situ water vapor and temperature measurements in the cloud indicated supersaturation with respect to ice throughout, with ice saturation ratios as large as 1.8. Even when the ice surface area density was as high as about 500 sq microns/cu cm, ice supersaturations of 20-30% were observed. Trajectory analysis shows that the air sampled near the tropopause on this day generally came from the north and cooled considerably during the previous few days. Examination of infrared satellite imagery along air parcel back trajectories from the WB-57 flight track indicates that the tropopause cloud layer formation was, in general, not simply left over ice from recently generated anvil cirrus. Simulations of cloud formation using time-height curtains of temperature along the trajectory paths show that the cloud could have formed in situ near the tropopause as the air was advected into the south Florida region and cooled to unusually low temperatures. If we assume a high threshold for ice nucleation via homogeneous freezing of aqueous sulfate aerosols, the model reproduces the observed cloud structure, ice crystal size distributions, and ice supersaturation statistics. Inclusion of observed gravity wave temperature perturbations in the simulations is essential to reproduce the observed cloud properties. Without waves, crystal number densities are too low, crystal sizes are too large, and the crystals fall out too fast, leaving very little cloud persisting at the end of the simulations. In the cloud simulations, coincidence of high supersaturations and high surface areas can be produced by either recent nucleation or sedimentation of crystals into supersaturated layers. The agreement between model results and observed supersaturations is improved somewhat if we assume that the steady state relative humidity within cirrus at T<200 K is enhanced by about 30%. The WB-57 measurements and the model results suggest that the cloud layer irreversibly dehydrated air near the tropopause.

Jensen, Eric↗

NASA NTRS2003

Comparison of mode estimation methods and application in molecular clock analysis

BACKGROUND: Distributions of time estimates in molecular clock studies are sometimes skewed or contain outliers. In those cases, the mode is a better estimator of the overall time of divergence than the mean or median. However, different methods are available for estimating the mode. We compared these methods in simulations to determine their strengths and weaknesses and further assessed their performance when applied to real data sets from a molecular clock study. RESULTS: We found that the half-range mode and robust parametric mode methods have a lower bias than other mode methods under a diversity of conditions. However, the half-range mode suffers from a relatively high variance and the robust parametric mode is more susceptible to bias by outliers. We determined that bootstrapping reduces the variance of both mode estimators. Application of the different methods to real data sets yielded results that were concordant with the simulations. CONCLUSION: Because the half-range mode is a simple and fast method, and produced less bias overall in our simulations, we recommend the bootstrapped version of it as a general-purpose mode estimator and suggest a bootstrap method for obtaining the standard error and 95% confidence interval of the mode.

Evolution, Molecular↗

NASA NTRS1977

Large eddy simulation of homogeneous isotropic turbulence

This paper presents results from a comprehensive investigation of large-eddy simulations of homogeneous isotropic turbulence. Calculations have been made using grid meshes, a number of flowfield filters, and two simple subgrid scale turbulence models. Particular attention has been paid to the degree of isotropy and the ability of the approach to predict higher-order statistical quantities (skewness and flatness). Directions for needed improvement in the simulation approach are indicated, and a plea is made to experimenters to process their data in a way which will facilitate comparison with large eddy simulations.

Ferziger, J. H.↗

NASA NTRS1981

Simulation of atmospheric turbulent gusts and gust gradients

An improved three-dimensional, nonrecursive model for atmospheric turbulence has been developed which provides for simulation of both instantaneous gusts and gust gradients along the flight path of the Space Shuttle. The one-dimensional gust and gust gradient spectral models, which form the basis for the simulation scheme, are developed from three-dimensional, von Karman spectra, integrated over finite limits based on the characteristic dimensions (length, width, and thickness) of the flight vehicle. By means of the simulation process nondimensional time series for both gusts and gust gradients have been generated and stored on a series of magnetic tapes for four altitude bands ranging from 0 to 10,000 meters. These Shuttle Simulated Turbulence Tapes (SSTT) have been validated by means of spectral and statistical analyses with satisfactory results. The dimensionless form of the time series coupled with the use of a von Karman spectral model (as opposed to a Dryden model) distinguish the SSTT from earlier turbulence simulation concepts.

Tatom, F. B.↗

NASA NTRS2008

The Integrated Medical Model: Statistical Forecasting of Risks to Crew Health and Mission Success

The Integrated Medical Model (IMM) helps capture and use organizational knowledge across the space medicine, training, operations, engineering, and research domains. The IMM uses this domain knowledge in the context of a mission and crew profile to forecast crew health and mission success risks. The IMM is most helpful in comparing the risk of two or more mission profiles, not as a tool for predicting absolute risk. The process of building the IMM adheres to Probability Risk Assessment (PRA) techniques described in NASA Procedural Requirement (NPR) 8705.5, and uses current evidence-based information to establish a defensible position for making decisions that help ensure crew health and mission success. The IMM quantitatively describes the following input parameters: 1) medical conditions and likelihood, 2) mission duration, 3) vehicle environment, 4) crew attributes (e.g. age, sex), 5) crew activities (e.g. EVA's, Lunar excursions), 6) diagnosis and treatment protocols (e.g. medical equipment, consumables pharmaceuticals), and 7) Crew Medical Officer (CMO) training effectiveness. It is worth reiterating that the IMM uses the data sets above as inputs. Many other risk management efforts stop at determining only likelihood. The IMM is unique in that it models not only likelihood, but risk mitigations, as well as subsequent clinical outcomes based on those mitigations. Once the mathematical relationships among the above parameters are established, the IMM uses a Monte Carlo simulation technique (a random sampling of the inputs as described by their statistical distribution) to determine the probable outcomes. Because the IMM is a stochastic model (i.e. the input parameters are represented by various statistical distributions depending on the data type), when the mission is simulated 10-50,000 times with a given set of medical capabilities (risk mitigations), a prediction of the most probable outcomes can be generated. For each mission, the IMM tracks which conditions occurred and decrements the pharmaceuticals and supplies required to diagnose and treat these medical conditions. If supplies are depleted, then the medical condition goes untreated, and crew and mission risk increase. The IMM currently models approximately 30 medical conditions. By the end of FY2008, the IMM will be modeling over 100 medical conditions, approximately 60 of which have been recorded to have occurred during short and long space missions.

Fitts, M. A.↗

NASA NTRS1980

Numerical simulation of plasma insulator interactions in space. Part 2: Dielectric effects

Any process which may be modeled statistically for a simple plasma particle may be modeled by the particle in cell technique. The success of the calculation is then dependent on having a large enough number of particles that the statistical treatment is meaningful. Thus it is possible to include the effects of secondary emission, backscattering, charge sticking and possibly dielectric breakdown, photoemission, and spallation. The first plasma dielectric interaction included in the computer code for simulating plasma insulator interactions is secondary electron emission. A calculated current density vs. voltage curve is presented and compared to an experimental curve.

Chaky, R. C.↗

NASA NTRSDate not supplied

Overview of Additively Manufactured TPS Proposed Flight Test and Earth Re-Entry Capsule Design

A flight mechanics overview is presented of an Earth flight test designed to investigate a novel, 3D printed thermal protection system (TPS) that is currently in development at NASA as part of the Additive Manufacturing of Thermal Protection Systems project. The project is pioneering a method to print a thermal protection system onto an entry vehicle forebody one layer at a time. This method reduces labor and complexity as compared to traditional manufacturing methods while increasing mission-dependent customization of through-depth materials properties. The flight test has three objectives. First, subject the forebody stagnation point of a capsule equipped with additively manufactured TPS (AMTPS) material to peak heat fluxes in excess of 100 $W/cm^2$. Second, capture in-flight data to enable flight reconstruction and AMTPS material thermal response model improvement. Third, recover the capsule with data storage and forebody AMTPS intact to enable post-flight inspection and analysis of AMTPS performance. The flight test trajectory is designed to achieve a peak stagnation point, cold-wall, entry heat flux of 135~$W/cm^2$. Flight mechanics simulations are performed using the Program to Optimize Simulated Trajectories II (POST2) and Monte-Carlo analysis yields statistical percentiles on vehicle performance at key points along the trajectory. Based on the flight mechanics analysis presented in this paper, a prototype capsule was designed, partially fabricated, and underwent preliminary component stress testing in preparation for fabrication of the flight unit capsule. The capsule outer mold line is a modified version of the heritage Mars Microprobe geometry. The capsule has a 0.356~m diameter, a 30~kg mass, and a hypersonic ballistic coefficient of 300~$kg/m^2$. Sensor selection is guided by flight dynamics simulations with the goal of resolving the re-entry heating pulse. On-board instrumentation include forebody and aftbody pressure sensors and thermocouples, a 9-axis IMU, a GPS receiver, and an Iridium satellite modem, all of which collect and store data throughout flight via on-board avionics systems. A two-stage parachute system is designed to decelerate the capsule to touchdown velocities that will not result in significant fracture or deformation of the charred AMTPS material at ground impact.

Flight Mechanics↗

NASA NTRS1999

Development of a Computational Chemical Vapor Deposition Model: Applications to Indium Nitride and Dicyanovinylaniline

A computational chemical vapor deposition (CVD) model is presented, that couples chemical reaction mechanisms with fluid dynamic simulations for vapor deposition experiments. The chemical properties of the systems under investigation are evaluated using quantum, molecular and statistical mechanics models. The fluid dynamic computations are performed using the CFD-ACE program, which can simulate multispecies transport, heat and mass transfer, gas phase chemistry, chemistry of adsorbed species, pulsed reactant flow and variable gravity conditions. Two experimental setups are being studied, in order to fabricate films of: (a) indium nitride (InN) from the gas or surface phase reaction of trimethylindium and ammonia; and (b) 4-(1,1)dicyanovinyl-dimethylaminoaniline (DCVA) by vapor deposition. Modeling of these setups requires knowledge of three groups of properties: thermodynamic properties (heat capacity), transport properties (diffusion, viscosity, and thermal conductivity), and kinetic properties (rate constants for all possible elementary chemical reactions). These properties are evaluated using computational methods whenever experimental data is not available for the species or for the elementary reactions. The chemical vapor deposition model is applied to InN and DCVA. Several possible InN mechanisms are proposed and analyzed. The CVD model simulations of InN show that the deposition rate of InN is more efficient when pulsing chemistry is used under conditions of high pressure and microgravity. An analysis of the chemical properties of DCVA show that DCVA dimers may form under certain conditions of physical vapor transport. CVD simulations of the DCVA system suggest that deposition of the DCVA dimer may play a small role in the film and crystal growth processes.

Cardelino, Carlos↗

NASA NTRS1992

The compositional structure and the effects of exothermicity in an unpremixed planar jet flame

Results are presented of direct numerical simulations (DNS) of a randomly perturbed compressible, spatially developing planar jet under the influence of a finite rate Arrhenius chemical reaction of the type F + O yields Product + Heat. The objectives of the simulations are to assess the compositional structure of the flame and to determine the influence of reaction exothermicity by means of statistical sampling of the data generated by DNS. It is shown that even with this idealized kinetics model, the simulated results exhibit features in accordance with experimental data. These results indicate that the Damkoehler number is an important parameter in determining the statistical composition of the reacting field, and that the results are insensitive to the mechanism by which this parameter is varied. It is demonstrated that as the intensity of mixing is increased and the effect of finite rate chemistry is more pronounced, the magnitudes of the ensemble mean and rms of the product mass fraction decrease, and those of the reactants mass fraction increase. Also, at higher mixing rates the joint probability density functions of the reactants' mass fractions shift towards higher values within the composition domain indicating a lower reactedness. These trends are consistent with those observed experimentally and are very useful in portraying the statistical structure of non-equilibrium diffusion flames. The DNS generated data are also utilized to examine the applicability of the 'laminar diffusion flamelet model' in predicting the rate of the reactant conversion with finite rate chemistry. This examination indicates that the bounds of the product formation scale reasonably well with those obtained by the flamelet model. However, in the range of the relatively low values of the Damkoehler numbers considered, the scatter of the results is substantially more than that to be completely determined by the closure.

Steinberger, C. J.↗