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 1,621 records · Page 90
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Cloud optical parameters as derived from the multispectral cloud radiometer

Simultaneous measurements of the liquid water content and particle size have assumed an important role in cloud physics as they help elucidate the mechanism of cloud particle formation and the mechanism of air mass-mixing in stratus clouds. Such measurements can reveal the modification of cloud air masses by anthropogenic aerosol particles (Coakley et al. 1987, Durkee 1989). Studies of the climatic impact of these modification processes on cloud microphysics seems to be urgent for understanding mechanisms of climate change. General Circulation Model (GCM) simulations can be improved by introducing a parameterization of cloud optical properties in terms of integrated liquid water content (liquid water path) and particle size (Slingo 1989). Motivated by the above mentioned circumstances, remote sensing techniques were developed for simultaneously retrieving the cloud optical thickness and effective particle radius, from which the liquid water path can be inferred. Statistical features of the cloud optical thickness (or liquid water path) and effective particle size for marine stratocumulus clouds are presented. These results were obtained during 4 days (7, 10, 13, and 16 July 1987) of observations with the Multispectral Cloud Radiometer (ER-2) and the Thematic Mapper (LANDSAT-5) during the First ISCCP Regional Experiment (FIRE).

Nakajima, Teruyuki↗
NASA NTRS1990

Structure of three-dimensional turbulent boundary layers

The changes that occur in the Reynolds-stress-producing motion when a cross-stream pressure gradient is applied to an initially two-dimensional turbulent flow are discussed. Two examples are used: (1) a temporal simulation of a channel flow with crossflow applied by a spanwise pressure gradient for t is greater than 0; and (2) a spatial simulation of the boundary layer on an infinite swept wing. Evidence examined to date suggests that the structural changes in the two cases are similar, but the mechanisms may be significantly different, even if effects peculiar to the viscous wall region are ignored. The results from (2) are provisional, based on too short a time series for accurate statistical averages to be obtained. We treat turbulence 'statistics' (solutions of the Navier-Stokes equations in three space dimensions and time) in the same way as experiments: both have limitations of accuracy but both are acceptable representations of real fluid flows.

Bradshaw, P.↗
NASA NTRS1999

Dynamical Signatures of Living Systems

One of the main challenges in modeling living systems is to distinguish a random walk of physical origin (for instance, Brownian motions) from those of biological origin and that will constitute the starting point of the proposed approach. As conjectured, the biological random walk must be nonlinear. Indeed, any stochastic Markov process can be described by linear Fokker-Planck equation (or its discretized version), only that type of process has been observed in the inanimate world. However, all such processes always converge to a stable (ergodic or periodic) state, i.e., to the states of a lower complexity and high entropy. At the same time, the evolution of living systems directed toward a higher level of complexity if complexity is associated with a number of structural variations. The simplest way to mimic such a tendency is to incorporate a nonlinearity into the random walk; then the probability evolution will attain the features of diffusion equation: the formation and dissipation of shock waves initiated by small shallow wave disturbances. As a result, the evolution never "dies:" it produces new different configurations which are accompanied by an increase or decrease of entropy (the decrease takes place during formation of shock waves, the increase-during their dissipation). In other words, the evolution can be directed "against the second law of thermodynamics" by forming patterns outside of equilibrium in the probability space. Due to that, a specie is not locked up in a certain pattern of behavior: it still can perform a variety of motions, and only the statistics of these motions is constrained by this pattern. It should be emphasized that such a "twist" is based upon the concept of reflection, i.e., the existence of the self-image (adopted from psychology). The model consists of a generator of stochastic processes which represents the motor dynamics in the form of nonlinear random walks, and a simulator of the nonlinear version of the diffusion equation which represents the mental dynamics. It has been demonstrated that coupled mental-motor dynamics can simulate emerging self-organization, prey-predator games, collaboration and competition, "collective brain," etc.

Zak, M.↗
NASA NTRS2008

Vibroacoustic Response of the NASA ACTS Spacecraft Antenna to Launch Acoustic Excitation

The Advanced Communications Technology Satellite was an experimental NASA satellite launched from the Space Shuttle Discovery. As part of the ground test program, the satellite s large, parabolic reflector antennas were exposed to a reverberant acoustic loading to simulate the launch acoustics in the Shuttle payload bay. This paper describes the modelling and analysis of the dynamic response of these large, composite spacecraft antenna structure subjected to a diffuse acoustic field excitation. Due to the broad frequency range of the excitation, different models were created to make predictions in the various frequency regimes of interest: a statistical energy analysis (SEA) model to capture the high frequency response and a hybrid finite element-statistical energy (hybrid FE-SEA) model for the low to mid-frequency responses. The strengths and limitations of each of the analytical techniques are discussed. The predictions are then compared to the measured acoustic test data and to a boundary element (BEM) model to evaluate the performance of the hybrid techniques.

Larko, Jeffrey M.↗
NASA NTRS2011

Recent Advances on INSAR Temporal Decorrelation: Theory and Observations Using UAVSAR

We review our recent advances in understanding the role of temporal decorrelation in SAR interferometry and polarimetric SAR interferometry. We developed a physical model of temporal decorrelation based on Gaussian-statistic motion that varies along the vertical direction in forest canopies. Temporal decorrelation depends on structural parameters such as forest height, is sensitive to polarization and affects coherence amplitude and phase. A model of temporal-volume decorrelation valid for arbitrary spatial baseline is discussed. We tested the inversion of this model to estimate forest height from model simulations supported by JPL/UAVSAR data and lidar LVIS data. We found a general good agreement between forest height estimated from radar data and forest height estimated from lidar data.

polarimetric SAR interferometry↗
NASA NTRS2017

The Substorm Cycle as Reproduced by Global MHD Models

Recently, Gordeev et al. (2015) suggested a method to test global MHD models against statistical empirical data. They showed that four community-available global MHD models supported by the Community Coordinated Modeling Center (CCMC) produce a reasonable agreement with reality for those key parameters (the magnetospheric size, magnetic field, and pressure) that are directly related to the large-scale equilibria in the outer magnetosphere. Based on the same set of simulation runs, here we investigate how the models reproduce the global loading-unloading cycle. We found that in terms of global magnetic flux transport, three examined CCMC models display systematically different response to idealized2 h north then 2 h south interplanetary magnetic field (IMF) Bz variation. The LFM model shows a depressed return convection and high loading rate during the growth phase as well as enhanced return convection and high unloading rate during the expansion phase, with the amount of loaded unloaded magnetotail flux and the growth phase duration being the closest to their observed empirical values during isolated substorms. Two other models exhibit drastically different behavior. In the BATS-R-US model the plasma sheet convection shows a smooth transition to the steady convection regime after the IMF southward turning. In the Open GGCM a weak plasma sheet convection has comparable intensities during both the growth phase and the following slow unloading phase. We also demonstrate potential technical problem in the publicly available simulations which is related to post processing interpolation and could affect the accuracy of magnetic field tracing and of other related procedures.

Gordeev, E.↗
NASA NTRS2019

Modeling Mg II During Solar Flares. II. Nonequilibrium Effects

To extract the information that the Mg II NUV spectra (observed by the Interface Region Imaging Spectrograph) carry about the chromosphere during solar flares, and to validate models of energy transport via model–data comparison, forward modeling is required. The assumption of statistical equilibrium (SE) is typically used to obtain the atomic level populations from snapshots of flare atmospheres, due to computational necessity. However, it is possible that relying on SE could lead to spurious results. We compare solving the atomic level populations via SE versus a non-equilibrium (NEQ) time-dependent approach. This was achieved using flare simulations from RADYN alongside the minority species version MS_RADYN from which the time-dependent Mg II atomic level populations and radiation transfer were computed in complete frequency redistribution. The impacts on the emergent profiles, light curves, line ratios, and formation heights are discussed. In summary we note that NEQ effects during flares are typically important only in the initial stages and for a short period following the cessation of the energy injection. An analysis of the timescales of ionization equilibrium reveals that for most of the duration of the flare, when the temperatures and densities are sufficiently enhanced, the relaxation timescales are short (τrelax < 0.1 s), so that the equilibrium solution is an adequate approximation. These effects vary with the size of the flare, however. In weaker flares, effects can be more pronounced. We recommend that NEQ effects be considered when possible but that SE is sufficient at most stages of the flare.

Kerr, Graham S.↗
NASA NTRSDate not supplied

The OSSE Framework at the NASA Global Modeling and Assimilation Office (GMAO)

The Global Modeling and Assimilation Office (GMAO) has been developing an Observing System Simulation Experiment (OSSE) framework. The OSSE system is currently based on the Nature Run (NR) developed by GMAO called G5NR. The G5NR is currently being tested for use in the GMAO OSSE framework. Synthetic observations have been generated based on the G5NR fields, including conventional observations, GPS, and satellite radiances. These synthetic observations are ingested using the Gridpoint Statistical Interpolation data assimilation system, with forecasts performed by the GEOS-5 model at 55 km/72L.

OSSE Framework↗
NASA NTRS2013

Monte Carlo Bayesian Inference on a Statistical Model of Sub-gridcolumn Moisture Variability Using High-resolution Cloud Observations : Sensitivity Tests and Results - Part II

Part I presented a Monte Carlo Bayesian method for constraining a complex statistical model of GCM sub-gridcolumn moisture variability using high-resolution MODIS cloud data, thereby permitting large-scale model parameter estimation and cloud data assimilation. This part performs some basic testing of this new approach, verifying that it does indeed significantly reduce mean and standard deviation biases with respect to the assimilated MODIS cloud optical depth, brightness temperature and cloud top pressure, and that it also improves the simulated rotational-Ramman scattering cloud optical centroid pressure (OCP) against independent (non-assimilated) retrievals from the OMI instrument. Of particular interest, the Monte Carlo method does show skill in the especially difficult case where the background state is clear but cloudy observations exist. In traditional linearized data assimilation methods, a subsaturated background cannot produce clouds via any infinitesimal equilibrium perturbation, but the Monte Carlo approach allows finite jumps into regions of non-zero cloud probability. In the example provided, the method is able to restore marine stratocumulus near the Californian coast where the background state has a clear swath. This paper also examines a number of algorithmic and physical sensitivities of the new method and provides guidance for its cost-effective implementation. One obvious difficulty for the method, and other cloud data assimilation methods as well, is the lack of information content in the cloud observables on cloud vertical structure, beyond cloud top pressure and optical thickness, thus necessitating strong dependence on the background vertical moisture structure. It is found that a simple flow-dependent correlation modification due to Riishojgaard (1998) provides some help in this respect, by better honoring inversion structures in the background state.

Monte Carlo Bayesian inference↗
NASA NTRS2022

Premature Deaths in Africa Due to Particulate Matter Under High and Low Warming Scenarios

Sustainable development and climate change mitigation can provide enormous public health benefits via improved air quality, especially in polluted areas. We use the latest state-of-the-art composition-climate model simulations to contrast human exposure to fine particulate matter in Africa under a “baseline” scenario with high material consumption, population growth, and warming to that projected under a sustainability scenario with lower consumption, population growth, and warming. Evaluating the mortality impacts of these exposures, we find that under the low warming scenario annual premature deaths due to PM2.5 are reduced by roughly 515,000 by 2050 relative to the high warming scenario (100,000, 175,000, 55,000, 140,000, and 45,000 in Northern, West, Central, East, and Southern Africa, respectively). This reduction rises to ∼800,000 by the 2090s, though by that time much of the difference is attributable to the projected differences in population. By contrast, during the first half of the century benefits are driven predominantly by emissions changes. Depending on the region, we find large intermodel spreads of ∼25%–50% in projected future exposures owing to different physics across the ensemble of 6 global models. The spread of projected deaths attributable to exposure to fine particulate matter, including uncertainty in the exposure-response function, are reduced in every region to ∼20%–35% by the non-linear exposure-response function. Differences between the scenarios have an even narrower spread of ∼5%–25% and are highly statistically significant in all regions for all models. These results provide valuable information for policy-makers to consider when working toward climate change mitigation and sustainable development goals.

PM2.5↗
NASA NTRS1979

Simulation of the pressure field near a jet by randomly distributed vortex rings

Fluctuations of the pressure field in the vicinity of a jet are simulated numerically by a flow model consisting of axially symmetric vortex rings with viscous cores submerged in a uniform stream. The time interval between the shedding of successive vortices is taken to be a random variable with a probability distribution chosen to match that from experiments. It is found that up to 5 diameters downstream of the jet exit, statistics of the computed pressure field are in good agreement with experimental results. Statistical comparisons are provided for the overall sound pressure level, the peak amplitude, and the Strouhal number based on the peak frequency of the pressure signals.

Fung, Y. T.↗
NASA NTRS1981

Waveband evaluation of proposed thematic mapper in forest cover classification

This study involved the evaluation of the characteristics of multispectral scanner data relative to forest cover type mapping, using NASA's NS-001 multispectral scanner to simulate the proposed Thematic Mapper (TM). The objectives were to determine: (1) the optimum number of wavebands to utilize in computer classifications of TM data; (2) which channel combinations provide the highest expected classification accuracy; and (3) the relative merit of each channel in the context of the cover classes examined. Transformed divergence was used as a measure of statistical distance between spectral class densities associated with each of twelve cover classes. The maximum overall mean pair-wise transformed divergence was used as the basis for evaluating all possible waveband combinations available for use in computer-assisted forest cover classifications.

Latty, R. S.↗
NASA NTRS1990

The influence of Lewis number and nonhomogeneous mixture on premixed turbulent flame structure

The structure of a premixed flame front in two-dimensional turbulence is investigated using full numerical simulation including heat release, variable fluid properties, and one-step Arrhenius chemistry. Non-unity Lewis number (Le) effects are described by comparing the local instantaneous turbulent flame structure to the steady one dimensional laminar flame structure for the same thermochemical parameters. Flame surface area, mean reactant consumption rate per unit area of flame, turbulent flame speed, and statistical descriptions of the flame geometry are also reported. Principal findings are that the local flame structure correlates strongly with the local flame curvature, while global properties depend both on strain and on curvature. Preliminary results for cases with nonhomogeneous reactant mixture strength are reported. Here, the emphasis is on the ability of a propagating flame to recover after encountering a fuel-lean pocket in an otherwise homogeneous mixture. Both one dimensional and two dimensional cases are described.

Haworth, D. C.↗
NASA NTRS1993

Concepts for conformal and body-axis attitude information for spatial awareness presented in a helmet-mounted display

A piloted simulation study has been conducted to evaluate two methods of presenting attitude information in a helmet-mounted display (HMD) for spatial awareness in a fighter airplane. One method, the body-axis concept, displayed the information relative to the body axis of the airplane. The quantitative results of this study favored the body-axis concept. Although no statistically significant differences were noted for either the pilots' understanding of roll attitude or target position, the pilots made pitch judgment errors three times more often with the conformal display. The subjective results showed the body-axis display did not cause attitude confusion, a prior concern with this display. In the posttest comments, the pilots overwhelmingly selected the body-axis display as the display of choice.

Jones, Denise R.↗
NASA NTRS1992

Passive microwave precipitation mapping and retrieval simulation

A numerical simulation of the satellite-based multispectral passive microwave mapping and precipitation retrieval process is described. The simulation contains three major components: (1) the forward radiative transfer calculation; (2) the sensor observation simulation; and (3) the precipitation retrieval algorithm. Secondary components include a Karhunen-Loeve transformation module and antenna pattern deconvolution module. The simulation facilitates the investigation of candidate precipitation retrieval algorithms using sub-Nyquist brightness temperature imagery, and the comparison of the relative merits of proposed radiometric systems such as the EOS multispectral imaging microwave radiometer (MIMR), the tropical rainfall measurement mission (TRMM) microwave imager (TRMM TMI), and the geosynchronous satellite sensors. A procedure for retrieving rain rates using a statistical iterative approach is outlined.

Skofronick, G. M.↗
NASA NTRS2004

Interaction of a Boundary Layer with a Turbulent Wake

The objective of this grant was to study the transition mechanisms on a flat-plate boundary layer interacting with the wake of a bluff body. This is a simplified configuration presented and designed to exemplify the phenomena that occur in multi-element airfoils, in which the wake of an upstream element impinges on a downstream one. Some experimental data is available for this configuration at various Reynolds numbers. The first task carried out was the implementation and validation of the immersed-boundary method. This was achieved by performing calculations of the flow over a cylinder at low and moderate Reynolds numbers. The low-Reynolds number results are discussed, which is enclosed as Appendix A. The high-Reynolds number results are presented in a paper in preparation for the Journal of Fluid Mechanics. We performed calculations of the wake-boundary-layer interaction at two Reynolds numbers, Re approximately equal to 385 and 1155. The first case is discussed and a comparison of the two calculations is reported. The simulations indicate that at the lower Reynolds number the boundary layer is buffeted by the unsteady Karman vortex street shed by the cylinder. This is shown: long streaky structures appear in the boundary layer in correspondence of the three-dimensionalities in the rollers. The fluctuations, however, cannot be self-sustained due to the low Reynolds-number, and the flow does not reach a turbulent state within the computational domain. In contrast, in the higher Reynolds-number case, boundary-layer fluctuations persist after the wake has decayed (due, in part, to the higher values of the local Reynolds number Re achieved in this case); some evidence could be observed that a self-sustaining turbulence generation cycle was beginning to be established. A third simulation was subsequently carried out at a higher Reynolds number, Re=3900. This calculation gave results similar to those of the Re=l155 case. Turbulence was established at fairly low Reynolds number, as a consequence of the high level of the free-stream perturbation. An instantaneous flow visualization for that case is shown. A detailed examination of flow statistics in the transitional and turbulent regions, including the evolution of the turbulent kinetic energy (TKE) budget and frequency spectra showed the formation and evolution of turbulent spots characteristic of the bypass transition mechanism. It was also observed that the turbulent eddies achieved an equilibrium, fully developed turbulent states first, as evidenced by the early agreement achieved by the terms in the TKE budget with those observed in turbulent flows. Once a turbulent Reynolds stress profile had been established, the velocity profile began to resemble a turbulent one, first in the inner region and later in the outer region of the wall layer. An extensive comparison of the three cases, including budgets, mean velocity and Reynolds stress profiles and flow visualization, is included. The results obtained are also presented.

Piomelli, Ugo↗
NASA NTRS2004

Spatial Distribution of Large Cloud Drops

The analysis of aircraft measurements of individual drop sizes in clouds suggests that for sufficiently small volumes the mean number of cloud drops with a given radius is proportional to volume powered by a drop-size dependent exponent. For abundant small drops present, the exponent is 1 as assumed in conventional approach. However, for rarer large drops, the exponents fall below unity. We show striking examples of the spatial distribution of large cloud drops using models that simulate the observed power laws. In contrast to currently used models that assume homogeneity and therefore a Poisson distribution of cloud drops, these models show strong drop clustering, the more so the larger the drops. The degree of clustering is determined by the observed exponents. The strong clustering of large drops arises naturally from the observed power-law statistics. This clustering has vital consequences for rain physics explaining how rain can form so fast and also helps explain why remotely sensed cloud drop size is generally biased.

Marshak, Alexander↗
NASA NTRS2010

Measurement of Insulation Compaction in the Cryogenic Fuel Tanks at Kennedy Space Center by Fast/Thermal Neutron Techniques

The liquid hydrogen and oxygen cryogenic storage tanks at John F. Kennedy Space Center (KSC) use expanded perlite as thermal insulation. Th ere is evidence that some of the perlite has compacted over time, com promising the thermal performance and possibly also structural integr ity of the tanks. Therefore an Non-destructive Testing (NDT) method for measuring the perlite density or void fraction is urgently needed. Methods based on neutrons are good candidates because they can readil y penetrate through the 1.75 cm outer steel shell and through the ent ire 120 cm thickness of the perlite zone. Neutrons interact with the nuclei of materials to produce characteristic gamma rays which are the n detected. The gamma ray signal strength is proportional to the atom ic number density. Consequently, if the perlite is compacted then the count rates in the individual peaks in the gamma ray spectrum will i ncrease. Perlite is a feldspathic volcanic rock made up of the major elements Si, AI, Na, K and 0 along with some water. With commercially available portable neutron generators it is possible to produce simul taneously fluxes of neutrons in two energy ranges: fast (14 MeV) and thermal (25 meV). Fast neutrons produce gamma rays by inelastic scatt ering which is sensitive to Fe and O. Thermal neutrons produce gamma rays by radiative capture in prompt gamma neutron activation (PGNA) and this is sensitive to Si, AI, Na, Kand H. Thus the two energy ranges produce complementary information. The R&D program has three phases: numerical simulations of neutron and gamma ray transport with MCNP s oftware, evaluation of the system in the laboratory on test articles and finally mapping of the perlite density in the cryogenic tanks at KSC. The preliminary MCNP calculations have shown that the fast/therma l neutron NDT method is capable of distinguishing between expanded an d compacted perlite with excellent statistics.

Livingston, R. A.↗