Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “PDF”

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 199 records · Page 11

Modeling of the reactant conversion rate in a turbulent shear flow

Results are presented of direct numerical simulations (DNS) of spatially developing shear flows under the influence of infinitely fast chemical reactions of the type A + B yields Products. The simulation results are used to construct the compositional structure of the scalar field in a statistical manner. The results of this statistical analysis indicate that the use of a Beta density for the probability density function (PDF) of an appropriate Shvab-Zeldovich mixture fraction provides a very good estimate of the limiting bounds of the reactant conversion rate within the shear layer. This provides a strong justification for the implementation of this density in practical modeling of non-homogeneous turbulent reacting flows. However, the validity of the model cannot be generalized for predictions of higher order statistical quantities. A closed form analytical expression is presented for predicting the maximum rate of reactant conversion in non-homogeneous reacting turbulence.

Frankel, S. H.↗

Pre-mixed flame simulations for non-unity Lewis numbers

A principal effect of turbulence on premixed flames in the flamelet region is to wrinkle the flame fronts. For non-unity Lewis numbers (Le), the local flame structure is altered in curved regions. This effect is examined using direct numerical simulations of the three dimensional, constant density, decaying isotropic turbulence with a single step, finite rate chemical reaction. Simulations of Lewis numbers 0.8, 1.0, and 1.2 are compared. The turbulent flame speed, S(sub T), increases as Le decreases. The correlation between S(sub T) and u prime found in previous Le = 1 simulations has a strong Lewis number dependency. The variance of the pdf of the flame curvature increases as Le decreases, indicating that the flames become more wrinkled. A strong correlation between local flame speed and curvature was found. For Le greater than 1, the flame speed increases in regions concave towards the products and decreases in convex regions. The opposite correlation was found for Le less than 1. The mean temperature of the products was also found to vary with Lewis number. For Le = 0.8, it is less than the adiabatic flame temperature and for Le = 1.2 it is greater.

Rutland, C. J.↗

A subgrid mixing model for LES of nonpremixed turbulent combustion

A new sub-grid mixing model for use in large eddy simulations of turbulent combustion is presented and applied to a hydrogen-air diffusion flame. The sub-grid model is based on Kerstein's Linear Eddy Model (1988). The model is first used to predict the mixing of a conserved scalar in a turbulent shear flow. The model correctly predicts the behavior of the pdf of the scalar field and displays a nonmarching peak at the preferred mixture fraction as the shear layer is traversed. It is then illustrated how a reduced chemical mechanism can be implemented within the linear eddy subgrid model formulation. The model is used to predict NO formation in hydrogen-air diffusion flame using a reduced chemical mechanism involving nine reactive scalars.

Mcmurtry, P. A.↗

Turbulent reacting flow computations including turbulence-chemistry interactions

A two-equation (k-epsilon) turbulence model has been extended to be applicable for compressible reacting flows. A compressibility correction model based on modeling the dilatational terms in the Reynolds stress equations has been used. A turbulence-chemistry interaction model is outlined. In this model, the effects of temperature and species mass concentrations fluctuations on the species mass production rates are decoupled. The effect of temperature fluctuations is modeled via a moment model, and the effect of concentration fluctuations is included using an assumed beta-pdf model. Preliminary results obtained using this model are presented. A two-dimensional reacting mixing layer has been used as a test case. Computations are carried out using the Navier-Stokes solver SPARK using a finite rate chemistry model for hydrogen-air combustion.

Narayan, J. R.↗

Conceptual aspects of a statistical-dynamical approach to represent landscape subgrid-scale heterogeneities in atmospheric models

A parameterization of land surfaces based on a statistical-dynamical approach is presented. With this approach, the most important characteristics of the soil-plant-atmosphere system that affect the partition of energy (e.g., plant stomatal conductance, soil humidity, and surface roughness) are represented by a pdf rather than by a single 'representative' value. A primary simplified version of this parameterization is used to estimate the land-surface energy fluxes that are produced at the grid scale by various distributions of stomatal conductance under a broad range of environmental conditions. To demonstrate the approach's potential, results are compared with the same fluxes calculated with a big leaf model using the mean stomatal conductance that corresponds to the distributions. Large absolute and relative differences are obtained between the two schemes for many combinations of stomatal conductance pdfs and environmental conditions.

Avissar, Roni↗

Comparison of laser transit anemometry data analysis techniques

Two techniques for the extraction of two-dimensional flow information from laser transit anemometry (LTA) data sets are presented and compared via a simulation study and experimental investigation. The methods are a probability density function (PDF) estimation technique and a marginal distribution analysis technique. The simulation study builds on the results of previous work and provides a quantification of the accuracy of both techniques for various LTA data acquisition scenarios. The experimental comparison consists of using an LTA system to survey the flow downstream of a turbulence generator in a small low-speed wind tunnel. The collected data sets are analyzed and compared.

Humphreys, William M., Jr.↗

Simulated meteorite impacts and volcanic explosions: Ejecta analyses and planetary implications

Past cratering studies have focused primarily on crater morphology. However, important questions remain about the nature of crater deposits. Phenomena that need to be studied include the distribution of shock effects in crater deposits and crater walls; the origin of mono- and polymict breccia; differences between local and distal ejecta; deformation induced by explosive volcanism; and the production of unshocked, high-speed ejecta that could form the lunar and martian meteorites found on the Earth. To study these phenomena, one must characterize ejecta and crater wall materials from impacts produced under controlled conditions. New efforts at LLNL simulate impacts and volcanism and study resultant deformation. All experiments use the two-stage light-gas gun facility at LLNL to accelerate projectiles to velocities of 0.2 to 4.3 km/s, including shock pressures of 0.9 to 50 GPa. We use granite targets and novel experimental geometries to unravel cratering processes in crystalline rocks. We have thus far conducted three types of simulations: soft recovery of ejecta, 'frozen crater' experiments, and an 'artificial volcano. Our ejecta recovery experiments produced a useful separation of impactites. Material originally below the projectile remained trapped there, embedded in the soft metal of the flyer plate. In contrast, material directly adjacent to the projectile was jetted away from the impact, producing an ejecta cone that was trapped in the foam recovery fixture. We find that a significant component of crater ejecta shows no signs of strong shock; this material comes from the near-surface 'interference zone' surrounding the impact site. This phenomenon explains the existence of unshocked meteorites on the Earth of lunar and martian origin. Impact of a large bolide on neighboring planets will produce high-speed, weakly shocked ejecta, which may be trapped by the Earth's gravitational field. 'Frozen crater' experiments show that the interference zone is highly localized; indeed, disaggregation does not extend beyond approx. 1.5 crater radii. A cone-shaped region extending downward from the impact site is completely disaggregated, including powdered rock that escaped into the projectile tube. Petrographic analysis of crater ejecta and wall material will be presented. Finally, study of ejecta from 0.9- and 1.3-GPa simulations of volcanic explosions reveal a complete lack of shock metamorphism. The ejecta shows no evidence of PDF's, amorphization, high-pressure phases, or mosaicism. Instead, all deformation was brittle, with fractures irregular (not planar) and most intergranular. The extent of fracturing was remarkable, with the entire sample reduced to fragments of gravel size and smaller.

Gratz, A. J.↗

Where's the Beaverhead beef?

Only rare quartz grains with single-set planar (1013) deformation features (PDF's) are present in breccia dikes found in association with uniformly oriented shatter cones that occur over an area 8 x 25 km. This suggests that the Beaverhead shocked rocks come from only the outer part of the central uplift of what must have been a large (greater than 100 km diameter) complex impact structure. An impact event of this magnitude on continental crust (thought to have occurred in late Precambrian or ealy Paleozoic time) could be expected to punctuate local geologic history. Furthermore, although it may now be covered, its scar should remain despite all the considerable subsequent erosion/deposition and tectonism since the impact. The following are three large-scale singularities or anomalies that may reflect the event and mark its source. (1) The Lemhi Arch is a major structural uplift that occurred in late Proterozoic-early Paleozoic time in East Central Idaho and caused the erosion of at least 4 km of sedimentary cover. This may be directly related to the impact. (2) Of the many thrust sheets comprising the Cordilleran belt, the Cabin plate that carries the shocked rocks is unique in that it alone intersected the crystalline basement. It also now marks the apex of the Southwest Montana Recess in the Sevier front. The basement uplift remaining from the impact may have constituted a mechanical obstacle to the advancing thrust sheets in Cretaceous time, causing the recess. (3) What could be interpreted as a roughly circular aeromagnetic anomaly approx. 70 km in diameter can be discerned in the state aeromagnetic map centered about 20 km southeast of Challis, Idaho, in the Lost River range. It is in approximately the right place, and ignoring the possibility that the anomalies have diverse causes and the circular pattern is coincidental, it may mark what remains of the buried central uplift structure.

Hargraves, R. B.↗

Center for Modeling of Turbulence and Transition (CMOTT): Research Briefs, 1992

The progress is reported of the Center for Modeling of Turbulence and Transition (CMOTT). The main objective of the CMOTT is to develop, validate and implement the turbulence and transition models for practical engineering flows. The flows of interest are three-dimensional, incompressible and compressible flows with chemical reaction. The research covers two-equation (e.g., k-e) and algebraic Reynolds-stress models, second moment closure models, probability density function (pdf) models, Renormalization Group Theory (RNG), Large Eddy Simulation (LES) and Direct Numerical Simulation (DNS).

Liou, William W.↗

A mapping closure for turbulent scalar mixing using a time-evolving reference field

A general mapping-closure approach for modeling scalar mixing in homogeneous turbulence is developed. This approach is different from the previous methods in that the reference field also evolves according to the same equations as the physical scalar field. The use of a time-evolving Gaussian reference field results in a model that is similar to the mapping closure model of Pope (1991), which is based on the methodology of Chen et al. (1989). Both models yield identical relationships between the scalar variance and higher-order moments, which are in good agreement with heat conduction simulation data and can be consistent with any type of epsilon(phi) evolution. The present methodology can be extended to any reference field whose behavior is known. The possibility of a beta-pdf reference field is explored. The shortcomings of the mapping closure methods are discussed, and the limit at which the mapping becomes invalid is identified.

Girimaji, Sharath S.↗

General probability-matched relations between radar reflectivity and rain rate

An improved method for transforming radar-observed reflectivities Ze into rain rate R is presented. The method is based on a formulation of a Ze-R function constrained such that (1) the radar-retrieved pdf of R and all of its moments are identical to those determined from the gauges over a sufficiently large domain, and (2) the fraction of the time that it is raining above a low but still has an accurately measurable rain intensity is identical for both the radar and for simultaneous measurements of collocated gauges on average. Data measured by a 1.65-deg beamwidth C-band radar and 22 gauges located in the vicinity of Darwin, Australia, are used. The resultant Ze-R functions show a strong range dependence, especially for the rain regimes characterized by strong reflectivity gradients and substantial attenuation. The application of these novel Ze-R functions to the radar data produces excellent matches to the gauge measurements without any systematic bias.

Rosenfeld, Daniel↗

Coherent anti-Stokes Raman spectroscopy temperature measurements in a hydrogen-fueled supersonic combustor

Coherent anti-Stokes Raman spectroscopy (CARS) thermometry has been used to obtain static temperature cross sections in a three-dimensional supersonic combustor flowfield. Data were obtained in three spanwise planes downstream of a single normal fuel injector which was located downstream of a rearward-facing step. The freestream flow was nominally Mach 2 and was combustion heated to a total temperature of 1440 K (yielding a static temperature of about 800 K in the freestream) to simulate the inflow to a combustor operating at a flight Mach number of about 5.4. Since a broadband probe laser was used an instantaneous temperature sample was obtained with each laser shot at a repetition rate of 10 Hz. Thus root-mean-square (rms) temperatures and temperature probability density functions (pdf's) were obtained in addition to mean temperatures.

Smith, Michael W.↗

Reactant conversion in homogeneous turbulence - Mathematical modeling, computational validations, and practical applications

The presently obtained closed-form analytical expressions, which predict the limiting rate of mean reactant conversion in homogeneous turbulent flows under the influence of a binary reaction, are derived via the single-point pdf method based on amplitude mapping closure. With this model, the maximum rate of the mean reactant's decay can be conveniently expressed in terms of definite integrals of the parabolic cylinder functions. The results obtained are shown to be in good agreement with data generated by direct numerical simulations.

Madnia, C. K.↗

Computational methods for efficient structural reliability and reliability sensitivity analysis

This paper presents recent developments in efficient structural reliability analysis methods. The paper proposes an efficient, adaptive importance sampling (AIS) method that can be used to compute reliability and reliability sensitivities. The AIS approach uses a sampling density that is proportional to the joint PDF of the random variables. Starting from an initial approximate failure domain, sampling proceeds adaptively and incrementally with the goal of reaching a sampling domain that is slightly greater than the failure domain to minimize over-sampling in the safe region. Several reliability sensitivity coefficients are proposed that can be computed directly and easily from the above AIS-based failure points. These probability sensitivities can be used for identifying key random variables and for adjusting design to achieve reliability-based objectives. The proposed AIS methodology is demonstrated using a turbine blade reliability analysis problem.

Wu, Y.-T.↗

Issues related to waveform computations for radar altimeter applications

An algorithm has been developed to model the average return power waveforms available from general radar altimeter systems, such as the Multimode Airborne Radar Altimeter (MARA) system operated at NASA Goddard Space Flight Center - Wallops Flight Facility. The algorithm is based on a convolutional model comprised of three functions: the average flat surface impulse response (FSIR), the radar system point target response (PTR), and the height pdf of the specular points on the sea surface. The FSIR is modified to account for the asymmetric antenna beam used by the MARA system, and then certain properties of this modified SIR are exploited to obtain closed-form expressions that can be rapidly evaluated. An FFT convolution routine is used to further speed up the computations. The result is an algorithm that can be used to study the effects of pointing errors in surface measurements.

Newkirk, Michael H.↗

Probability-matched reflectivity-rainfall relations for a hurricane from aircraft observations

The relation between the distribution of equivalent reflectivity, Z(sub e), measured by an airborne C-band radar and that for rain rate, R, measured concurrently by a disdrometer on the same aircraft in the eyewall and outer bands of Hurricane Anita (1977) was investigated using the probability-matching method (PMM). It was found that the pdf of R is much narrower in the outer rainbands of Hurricane Anita than in the eyewall, and that the mean rain rates are also much smaller. It is shown that the use of PMM provides more realistic relations than the conventional power-law Z-R regression relations based upon scatter-plots of disdrometer measurements, showing higher eyewall reflectivities for the same rain rates at rates in excess of 6.3 mm/h.

Marks, Frank D., Jr.↗

CFD validation of subsonic turbulent planar shear layers

The primary objective of the present study is to assess the limitations and capabilities of RPLUS in predicting the entrainment, mixing, and burning characteristics of a high subsonic planar shear layer; these computations are performed in conjunction with an experiment being conducted at NASA Lewis Research Center. Turbulence is modeled by a two-equation k-epsilon closure modified for compressibility effects. Combustion is modeled by two different models: one in which a finite-rate laminar chemistry model is used for H2-O2 oxidation, and the other based on a composition joint pdf approach to account for the turbulence-chemistry interaction effects. Between upwind and centered differencing schemes, minimal differences are observed for the flow properties of a nonreacting shear layer case studied. In comparison with experimental data, computed growth rates are underpredicted in both reacting and nonreacting cases of a shear layer, but more severely for the reacting flows. However, the qualitative agreement between computation and experiment is reasonably good for the other observable characteristics, and indicative of potential for further improvements.

Lai, H. T.↗