Search NASA⌕ Search

SEARCH · Search NASA

Results for “Statistical homogeneity”

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 37 records · Page 2

LRAT: Lightning Radiative Transfer

In this report, we extend to cloud physics the work done for single and multiple scattering of electromagnetic waves. We consider the scattering of light, visible or infrared, by a spherical cloud represented by a statistically homogeneous ensemble of configurations of N identical spherical water droplets whose centers are uniformly distributed in its volume V. The ensemble is specified by the average number rho of scatterers in unit volume and by rho f(R) with f(R) as the distribution function for separations R of pairs. The incident light, vector-phi(sub 0) a plane electromagnetic wave with harmonic time dependence, is from outside the cloud. The propagation parameter kappa(sub 0) and the index of refraction eta(sub 0) determine physically the medium outside the distribution of scatterers. We solve the interior problem separately to obtain the bulk parameters for the scatterer equivalent to the ensemble of spherical droplets. With the interior solution or the equivalent medium approach, the multiple scattering problem is reduced to that of an equivalent single scatterer excited from outside illumination. A dispersion relation which determines the bulk propagation parameter K and the bulk index of refraction eta of the cloud is given in terms of the vector equivalent scattering amplitude vector-G and the dyadic scattering amplitude tilde-g of the single object in isolation. Based on this transfer model we will have the ability to consider clouds composed of inhomogeneous distribution of water and/or ice particles and we will be able to take into account particle size distributions within the cloud. We will also be able to study the effects of cloud composition (i.e., particle shape, size, composition, orientation, location) on the polarization of the single or the multiple scattered waves. Finally, this study will provide a new starting point for studying the problem of lightning radiative transfer.

Phanord, Dieudonne D.↗

Two space scatterer formalism calculation of bulk parameters of thunderclouds

In a previous study, we used a modified two-space scatterer formalism of Twersky to establish for a cloud modeled as a statistically homogeneous distribution of spherical water droplets, the dispersion relations that determine its bulk propagation numbers and bulk indexes of refraction in terms of the vector equivalent scattering amplitude and the dyadic scattering amplitude of the single water droplet in isolation. The results were specialized to the forward direction of scattering while demanding that the scatterers preserve the incident polarization. We apply this approach to obtain specific numerical values for the macroscopic parameters of the cloud. We work with a cloud of density rho = 100 cm(exp -3), a wavelength lambda = 0.7774 microns, and with spherical water droplets of common radius alpha = 10 microns. In addition, the scattering medium is divided into three parts, the medium outside the cloud, moist air (the medium inside the cloud but outside the droplets), and the medium inside the spherical water droplets. The results of this report are applicable to a cloud of any geometry since the boundary does not interfere with the calculations. Also, it is important to notice the plane wave nature of the incidence wave in the moist atmosphere.

Phanord, Dieudonne D.↗

A Lagrangian dynamic subgrid-scale model turbulence

A new formulation of the dynamic subgrid-scale model is tested in which the error associated with the Germano identity is minimized over flow pathlines rather than over directions of statistical homogeneity. This procedure allows the application of the dynamic model with averaging to flows in complex geometries that do not possess homogeneous directions. The characteristic Lagrangian time scale over which the averaging is performed is chosen such that the model is purely dissipative, guaranteeing numerical stability when coupled with the Smagorinsky model. The formulation is tested successfully in forced and decaying isotropic turbulence and in fully developed and transitional channel flow. In homogeneous flows, the results are similar to those of the volume-averaged dynamic model, while in channel flow, the predictions are superior to those of the plane-averaged dynamic model. The relationship between the averaged terms in the model and vortical structures (worms) that appear in the LES is investigated. Computational overhead is kept small (about 10 percent above the CPU requirements of the volume or plane-averaged dynamic model) by using an approximate scheme to advance the Lagrangian tracking through first-order Euler time integration and linear interpolation in space.

Meneveau, C.↗

Constraints on galactic distributions of gamma-ray burst sources from BATSE observations

The paradigm that gamma-ray bursts originate from Galactic sources is studied in detail using the angular and intensity distributions observed by the Burst and Transient Source Experiment (BATSE) on NASA's Compton Gamma Ray Observatory (CGRO). Monte Carlo models of gamma-ray burst spatial distributions and luminosity functions are used to simulate bursts, which are then folded through mathematical models of BATSE selection effects. The observed and computed angular intensity distributions are analyzed using modifications of standard statistical homogeneity and isotropy studies. Analysis of the BATSE angular and intensity distributions greatly constrains the origins and luminosities of burst sources. In particular, it appears that no single population of sources confined to a Galactic disk, halo, or localized spiral arm satisfactorily explains BATSE observations and that effects of the burst luminosity function are secondary when considering such models. One family of models that still satisfies BATSE observations comprises sources located in an extended spherical Galactic corona. Coronal models are limited to small ranges of burst luminosity and core radius, and the allowed parameter space for such models shrinks with each new burst BATSE observes. Multiple-population models of bursts are found to work only if (1) the primary population accounts for the general isotropy and inhomogeneity seen in the BATSE observations and (2) secondary populations either have characteristics similar to the primary population or contain numbers that are small relative to the primary population.

Hakkila, Jon↗

LES on unstructured deforming meshes: Towards reciprocating IC engines

A variable explicit/implicit characteristics-based advection scheme that is second-order accurate in space and time has been developed recently for unstructured deforming meshes (O'Rourke & Sahota 1996a). To explore the suitability of this methodology for Large-Eddy Simulation (LES), three subgrid-scale turbulence models have been implemented in the CHAD CFD code (O'Rourke & Sahota 1996b): a constant-coefficient Smagorinsky model, a dynamic Smagorinsky model for flows having one or more directions of statistical homogeneity, and a Lagrangian dynamic Smagorinsky model for flows having no spatial or temporal homogeneity (Meneveau et al. 1996). Computations have been made for three canonical flows, progressing towards the intended application of in-cylinder flow in a reciprocating engine. Grid sizes were selected to be comparable to the coarsest meshes used in earlier spectral LES studies. Quantitative results are reported for decaying homogeneous isotropic turbulence, and for a planar channel flow. Computations are compared to experimental measurements, to Direct-Numerical Simulation (DNS) data, and to Rapid-Distortion Theory (RDT) where appropriate. Generally satisfactory evolution of first and second moments is found on these coarse meshes; deviations are attributed to insufficient mesh resolution. Issues include mesh resolution and computational requirements for a specified level of accuracy, analytic characterization of the filtering implied by the numerical method, wall treatment, and inflow boundary conditions. To resolve these issues, finer-mesh simulations and computations of a simplified axisymmetric reciprocating piston-cylinder assembly are in progress.

Haworth, D. C.↗

Diffusion of Sound Waves in a Turbulent Atmosphere

The directional and frequency diffusion of a plane monochromatic 2 sound wave in statistically homogeneous, isotropic, and stationary turbulence is analyzed theoretically. The treatment is based on the diffusion equation for the energy density of sound waves, using the scattering cross section derived by Kraichnan for the type of turbulence assumed here. A form for the frequency-wave number spectrum of the turbulence is adopted which contains the pertinent parameters of the flow and is adapted to ease of calculation. A new approach to the evaluation of the characteristic period of the flow is suggested. This spectrum is then related to the scattering cross section. Finally, a diffusion equation is derived as a small-angle scattering approximation to the rigorous transport equation. The rate of spread of the incident wave in frequency and direction is calculated, as well as the power spectrum and autocorrelation for the wave.

Lyon, Richard H.↗

Ensemble Averaged Probability Density Function (APDF) for Compressible Turbulent Reacting Flows

In this paper, we present a concept of the averaged probability density function (APDF) for studying compressible turbulent reacting flows. The APDF is defined as an ensemble average of the fine grained probability density function (FG-PDF) with a mass density weighting. It can be used to exactly deduce the mass density weighted, ensemble averaged turbulent mean variables. The transport equation for APDF can be derived in two ways. One is the traditional way that starts from the transport equation of FG-PDF, in which the compressible Navier- Stokes equations are embedded. The resulting transport equation of APDF is then in a traditional form that contains conditional means of all terms from the right hand side of the Navier-Stokes equations except for the chemical reaction term. These conditional means are new unknown quantities that need to be modeled. Another way of deriving the transport equation of APDF is to start directly from the ensemble averaged Navier-Stokes equations. The resulting transport equation of APDF derived from this approach appears in a closed form without any need for additional modeling. The methodology of ensemble averaging presented in this paper can be extended to other averaging procedures: for example, the Reynolds time averaging for statistically steady flow and the Reynolds spatial averaging for statistically homogeneous flow. It can also be extended to a time or spatial filtering procedure to construct the filtered density function (FDF) for the large eddy simulation (LES) of compressible turbulent reacting flows.

Shih, Tsan-Hsing↗

Statistics of pressure and pressure gradient in homogeneous isotropic turbulence

The statistics of pressure and pressure gradient in stationary isotropic turbulence are measured within direct numerical simulations at low to moderate Reynolds numbers. It is found that the one-point pdf of the pressure is highly skewed and that the pdf of the pressure gradient is of stretched exponential form. The power spectrum of the pressure P(k) is found to be larger than the corresponding spectrum P(sub G)(k) computed from a Gaussian velocity field having the same energy spectrum as that of the DNS field. The ratio P(k)/P(sub G)(k), a measure of the pressure-field intermittence, grows with wavenumber and Reynolds number as -R(sub lambda)(exp 1/2)log(k/k(sub d)) for k less than k(sub d)/2 where k(sub d) is the Kolmogorov wavenumber. The Lagrangian correlations of pressure gradient and velocity are compared and the Lagrangian time scale of the pressure gradient is observed to be much shorter than that of the velocity.

Gotoh, T.↗

Homogeneity of spatial correlation statistics of tropical oceanic rainfall

The possibility of uniform horizontal correlation scales for tropical oceanic rainfall has been examined by a study of satellite-observed microwave data as a proxy measure of rain rates. From the brightness temperatures from the electrically scanning microwave radiometer on Nimbus 5 near nadir during the year 1974, the mean spatial autocorrelation function as a function of simultaneous pixel separation was calculated in each 5 by 5 deg grid box over the tropical Pacific and Atlantic for each season. The equal-time spatially lagged correlations were compared for geographical dependence to investigate the hypothesis of homogeneity. A simple model of the spatial statistics of the microwave brightness temperatures was used, consisting of a mixture of uncorrelated spatial white noise incoherently superimposed on a spatially correlated field (spatial red noise). The red noise signals are presumed to be generated by convective activity in the tropical atmosphere. The parameters of the red noise are consistent with this scheme over the tropical oceans, yielding a uniform spatial scale of about 50 km throughout.

Shin, Kyung-Sup↗

The Statistical Mechanics of Ideal Homogeneous Turbulence

Plasmas, such as those found in the space environment or in plasma confinement devices, are often modeled as electrically conducting fluids. When fluids and plasmas are energetically stirred, regions of highly nonlinear, chaotic behavior known as turbulence arise. Understanding the fundamental nature of turbulence is a long-standing theoretical challenge. The present work describes a statistical theory concerning a certain class of nonlinear, finite dimensional, dynamical models of turbulence. These models arise when the partial differential equations describing incompressible, ideal (i.e., nondissipative) homogeneous fluid and magnetofluid (i.e., plasma) turbulence are Fourier transformed into a very large set of ordinary differential equations. These equations define a divergenceless flow in a high-dimensional phase space, which allows for the existence of a Liouville theorem, guaranteeing a distribution function based on constants of the motion (integral invariants). The novelty of these particular dynamical systems is that there are integral invariants other than the energy, and that some of these invariants behave like pseudoscalars under two of the discrete symmetry transformations of physics, parity, and charge conjugation. In this work the 'rugged invariants' of ideal homogeneous turbulence are shown to be the only significant scalar and pseudoscalar invariants. The discovery that pseudoscalar invariants cause symmetries of the original equations to be dynamically broken and induce a nonergodic structure on the associated phase space is the primary result presented here. Applicability of this result to dissipative turbulence is also discussed.

Shebalin, John V.↗

Coherent Eigenmodes in Homogeneous MHD Turbulence

The statistical mechanics of Fourier models of ideal, homogeneous, incompressible magnetohydrodynamic (MHD) turbulence is discussed, along with their relevance for dissipative magnetofluids. Although statistical theory predicts that Fourier coefficients of fluid velocity and magnetic field are zero-mean random variables, numerical simulations clearly show that certain coefficients have a non-zero mean value that can be very large compared to the associated standard deviation, i.e., we have coherent structure. We use eigenanalysis of the modal covariance matrices in the probability density function to explain this phenomena in terms of `broken ergodicity', which is defined to occur when dynamical behavior does not match ensemble predictions on very long time-scales. We provide examples from 2-D and 3-D magnetohydrodynamic simulations of homogeneous turbulence, and show new results from long-time simulations of MHD turbulence with and without a mean magnetic field

Shebalin, John V.↗

Problems in the simulation of atmospheric boundary layer flows

The realistic simulation of flow in the atmospheric boundary layers at heights greater than two kilometers is discussed. Information concerning horizontally homogeneous and statistically stationary atmospheric boundary layer flows is presented. The problems related to the incorporation of the information into atmospheric wind simulation programs are analyzed. The information which the meteorologist must acquire in order to provide a basis for improving the simulation of atmospheric boundary flows is explained.

Fichtl, G. H.↗

A new image enhancement algorithm with applications to forestry stand mapping

The theory and applications are presented of a new image enhancement algorithm which refines computer classification maps of multispectral data. The refinement eliminates connected sets smaller than a prespecified size and merges them to the surrounding area. Conventional practices in forestry timber stand mapping requires small geographic areas to be absorbed by surrounding large areas to form homogeneous stands. This homogeneity is often incompatible with the statistical formulation of homogeneity. Elements within a timber stand which should be labeled as one feature often correspond to more than one class mapped by existing computer classification techniques. The new algorithm is designed to postprocess classification maps to result in more usable timber stand maps. The new image enhancement technique is compared with an accepted neighbor-checking postprocessing technique, demonstrating the superiority of the new technique for forestry stand mapping.

Kan, E. P.↗

Direct Numerical Simulation of Transitional Multicomponent-Species Gaseous and Multicomponent-Liquid Drop-Laden Mixing

A model of multicomponent-liquid (MC-liquid) drop evaporation in a three-dimensional mixing layer is here exercised at larger Reynolds numbers than in a previous study, and transitional states are obtained. The gas phase is followed in an Eulerian frame and the multitude of drops is described in a Lagrangian frame. Complete coupling between phases is included with source terms in the gas conservation equations accounting for the drop/flow interaction in terms of drop drag, drop heating and species evaporation. The liquid composition, initially specified as a single-Gamma (SG) probability distribution function (PDF) depending on the molar mass is allowed to evolve into a linear combination of two SGPDFs, called the double-Gamma PDF (DGPDF). The compositions of liquid and vapor emanating from the drops are calculated through four moments of the DGPDFs, which are drop-specific and location-specific, respectively. The mixing layer is initially excited to promote the double pairing of its four initial spanwise vortices into an ultimate vortex in which small scales proliferate. Simulations are performed for four liquids of different compositions and the effect of the initial mass loading and initial free-stream gas temperature are explored. For reference, Simulations are also performed for gaseous multicomponent mixing layers for which the effect of Reynolds number is investigated. The results encompass examination of the global layer characteristics, flow visualizations and homogeneous-plane statistics at transition. Comparisons are performed with previous pre-transitional MC-liquid simulations and with transitional single-component (SC) liquid studies. It is found that MCC flows at transition, the classical energy cascade is of similar strength, but that the smallest scales contain orders of magnitude less energy than SC flows, which is confirmed by the larger viscous dissipation in the former case. Contrasting to pre-transitional MC flows, the vorticity and drop organization depend on the initial gas temperature, this being due to the drop/turbulence coupling. The vapor-composition mean molar mass and standard deviation distributions strongly correlate with the initial liquid-composition PDF; such a correlation only exists for the magnitude of the mean but not for that of the standard deviation. Unlike in pre-transitional situations, regions of large composition standard deviation no longer necessarily coincide with regions of large mean molar mass. The kinetic energy, rotational and composition characteristics, and dissipation are liquid specific and the variation among liquids is amplified with increasing free-stream gas temperature. Eulerian and Lagrangian statistics of gas-phase quantities show that the different. Observation framework may affect the perception of the flow characteristics. The gas composition, of which the first four moments are calculated, is shown to be close to, but distinct from a SGPDF. The PDF of the scalar dissipation rate is calculated for drop-laden layers and is shown to depart more significantly from the typically assumed Gaussian in gaseous flows than experimentally measured gaseous scalar dissipation rates, this being attributed to the increased heterogeneity due to drop/flow interactions.

tubulent multicompenent-species↗

Statistical Calibration and Validation of a Homogeneous Ventilated Wall-Interference Correction Method for the National Transonic Facility

Wind tunnel experiments will continue to be a primary source of validation data for many types of mathematical and computational models in the aerospace industry. The increased emphasis on accuracy of data acquired from these facilities requires understanding of the uncertainty of not only the measurement data but also any correction applied to the data. One of the largest and most critical corrections made to these data is due to wall interference. In an effort to understand the accuracy and suitability of these corrections, a statistical validation process for wall interference correction methods has been developed. This process is based on the use of independent cases which, after correction, are expected to produce the same result. Comparison of these independent cases with respect to the uncertainty in the correction process establishes a domain of applicability based on the capability of the method to provide reasonable corrections with respect to customer accuracy requirements. The statistical validation method was applied to the version of the Transonic Wall Interference Correction System (TWICS) recently implemented in the National Transonic Facility at NASA Langley Research Center. The TWICS code generates corrections for solid and slotted wall interference in the model pitch plane based on boundary pressure measurements. Before validation could be performed on this method, it was necessary to calibrate the ventilated wall boundary condition parameters. Discrimination comparisons are used to determine the most representative of three linear boundary condition models which have historically been used to represent longitudinally slotted test section walls. Of the three linear boundary condition models implemented for ventilated walls, the general slotted wall model was the most representative of the data. The TWICS code using the calibrated general slotted wall model was found to be valid to within the process uncertainty for test section Mach numbers less than or equal to 0.60. The scatter among the mean corrected results of the bodies of revolution validation cases was within one count of drag on a typical transport aircraft configuration for Mach numbers at or below 0.80 and two counts of drag for Mach numbers at or below 0.90.

Walker, Eric Lee↗