Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Probability density function”

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 469 records · Page 26

Airborne lidar and radiometric observations of PBL- and low clouds

Boundary layer- and low altitude clouds over open ocean and continent areas have been studied during several field campaigns since mid-1990 using the French airborne backscatter lidar LEANDRE in conjunction with on-board IR and visible radiometers. LEANDRE is an automatic system, and a modification of the instrumental parameters, when airborne, is computer controlled through an operator keyboard. The vertical range squared lidar signals and instrument status are displayed in real time on two dedicated monitors. The lidar is used either down- or up-looking while the aircraft is flying above or below clouds. A switching of the viewing configuration takes about a minute. The lidar measurements provide a high resolution description of cloud morphology and holes in cloud layers. The flights were conducted during various meteorological conditions on single or multilayer stratocumulus and cumulus decks. Analysis on a single shot basis of cloud top (or bottom) altitude and a plot of the corresponding histogram allows one to determine a probability density function (PDF). The preliminary results show the PDFs for cloud top are not Gaussian and symmetric about the mean value. The skewness varies with atmospheric conditions. An example of results recorded over the Atlantic ocean near Biarritz is displayed, showing: (1) the range squared lidar signals as a function of time (here 100 s corresponds to about 8 km, 60 shots are averaged on horizontal); the Planetary Boundary Layer (PBL) - up to 600 m - is observed at the beginning of the leg as well as on surface returns, giving an indication of the porosity; (2) the cloud top altitude variation between 2.4 to 2.8 km during the 150 to 320 s section; and (3) the corresponding PDF. Similar results are obtained on stratocumulus over land. Single shot measurements can be used also to determine an optical porosity at a small scale as well as a fractional cloudiness at a larger scale. A comparison of cloud top altitude retrieved from lidar and narrowbeam IR radiometer is conducted to study the scale integration problem. A good agreement within less than 100 m relies on spatial uniformity and an optically thick layer. In the presence of holes, a discrepancy is observed. This is illustrated in figure 2, displaying as a function of time (1) the lidar signals; (2) the target temperature (either clouds or sea surface) retreived from a narrowbeam IR radiometer, 17 C is the sea surface temperature on that day; and (3) the visible flux, linked to cloud albedo, measured by a pyranometer. In preparation of ASTEX, down- and up-looking measurements where conducted on stratocumulus clouds over the Atlantic Ocean near Quimper in Brittany. Depending on the flight pattern orientation with respect to the wind, the top and bottom cloud morphologies are different. Preliminary results are given on cloud morphology, cloud top PDFs, optical porosity, fractional cloudiness, and comparison of lidar and radiometric measurements.

Flamant, P. H.↗

Towards understanding turbulent scalar mixing

In an effort towards understanding turbulent scalar mixing, we study the effect of molecular mixing, first in isolation and then by accounting for the effects of the velocity field. The chief motivation for this approach stems from the strong resemblance of the scalar probability density function (PDF) obtained from the scalar field evolving from the heat conduction equation that arises in a turbulent velocity field. However, the evolution of the scalar dissipation is different for the two cases. We attempt to account for these differences, which are due to the velocity field, using a Lagrangian frame analysis. After establishing the usefulness of this approach, we use the heat-conduction simulations (HCS), in lieu of the more expensive direct numerical simulations (DNS), to study many of the less understood aspects of turbulent mixing. Comparison between the HCS data and available models are made whenever possible. It is established that the beta PDF characterizes the evolution of the scalar PDF during mixing from all types of non-premixed initial conditions.

Girimaji, Sharath S.↗

Numerical modeling for dilute and dense sprays

We have successfully implemented a numerical model for spray-combustion calculations. In this model, the governing gas-phase equations in Eulerian coordinate are solved by a time-marching multiple pressure correction procedure based on the operator-splitting technique. The droplet-phase equations in Lagrangian coordinate are solved by a stochastic discrete particle technique. In order to simplify the calculation procedure for the circulating droplets, the effective conductivity model is utilized. The k-epsilon models are utilized to characterize the time and length scales of the gas phase in conjunction with turbulent modulation by droplets and droplet dispersion by turbulence. This method entails random sampling of instantaneous gas flow properties and the stochastic process requires a large number of computational parcels to produce the satisfactory dispersion distributions even for rather dilute sprays. Two major improvements in spray combustion modelings were made. Firstly, we have developed a probability density function approach in multidimensional space to represent a specific computational particle. Secondly, we incorporate the Taylor Analogy Breakup (TAB) model for handling the dense spray effects. This breakup model is based on the reasonable assumption that atomization and drop breakup are indistinguishable processes within a dense spray near the nozzle exit. Accordingly, atomization is prescribed by injecting drops which have a characteristic size equal to the nozzle exit diameter. Example problems include the nearly homogeneous and inhomogeneous turbulent particle dispersion, and the non-evaporating, evaporating, and burning dense sprays. Comparison with experimental data will be discussed in detail.

Chen, C. P.↗

Bit error rate performance of pi/4-DQPSK in a frequency-selective fast Rayleigh fading channel

The bit error rate (BER) performance of pi/4-differential quadrature phase shift keying (DQPSK) modems in cellular mobile communication systems is derived and analyzed. The system is modeled as a frequency-selective fast Rayleigh fading channel corrupted by additive white Gaussian noise (AWGN) and co-channel interference (CCI). The probability density function of the phase difference between two consecutive symbols of M-ary differential phase shift keying (DPSK) signals is first derived. In M-ary DPSK systems, the information is completely contained in this phase difference. For pi/4-DQPSK, the BER is derived in a closed form and calculated directly. Numerical results show that for the 24 kBd (48 kb/s) pi/4-DQPSK operated at a carrier frequency of 850 MHz and C/I less than 20 dB, the BER will be dominated by CCI if the vehicular speed is below 100 mi/h. In this derivation, frequency-selective fading is modeled by two independent Rayleigh signal paths. Only one co-channel is assumed in this derivation. The results obtained are also shown to be valid for discriminator detection of M-ary DPSK signals.

Liu, Chia-Liang↗

Mixing in straight shear layers

Planar laser-induced fluorescence measurements were performed in a liquid plane mixing layer to extract the probability density function (pdf) of the mixture fraction of a passive scalar across the layer. Three Reynolds number (Re) cases were studied, 10,000, 33,000 and 90,000, with Re based on velocity difference and visual thickness. The results show that a non-marching pdf (central hump invariant from edge to edge of the layer) exists for Re = 10,000 but that a marching type pdf characterizes the Re = 33,000 and Re = 90,000 cases. For all cases, a broad range of mixture fraction values is found at each location across the layer. Streamwise and spanwise ramps across the layer, and structure-to-structure variation were observed and are believed to be responsible for the above behavior of the composition field. Tripping the boundary layer on the high-speed side of the splitter plate for each of the above three cases resulted in increased three-dimensionality and a change in the composition field. Average and average mixed fluid compositions are reported for all cases.

Karasso, P. S.↗

Probabilistic material degradation model for aerospace materials subjected to high temperature, mechanical and thermal fatigue, and creep

A probabilistic general material strength degradation model has been developed for structural components of aerospace propulsion systems subjected to diverse random effects. The model has been implemented in two FORTRAN programs, PROMISS (Probabilistic Material Strength Simulator) and PROMISC (Probabilistic Material Strength Calibrator). PROMISS calculates the random lifetime strength of an aerospace propulsion component due to as many as eighteen diverse random effects. Results are presented in the form of probability density functions and cumulative distribution functions of lifetime strength. PROMISC calibrates the model by calculating the values of empirical material constants.

Boyce, L.↗

Modeling turbulent/chemistry interactions using assumed pdf methods

Two assumed probability density functions (pdfs) are employed for computing the effect of temperature fluctuations on chemical reaction. The pdfs assumed for this purpose are the Gaussian and the beta densities of the first kind. The pdfs are first used in a parametric study to determine the influence of temperature fluctuations on the mean reaction-rate coefficients. Results indicate that temperature fluctuations significantly affect the magnitude of the mean reaction-rate coefficients of some reactions depending on the mean temperature and the intensity of the fluctuations. The pdfs are then tested on a high-speed turbulent reacting mixing layer. Results clearly show a decrease in the ignition delay time due to increases in the magnitude of most of the mean reaction rate coefficients.

Gaffney, R. L, Jr.↗

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.↗

A PDF closure model for compressible turbulent chemically reacting flows

The objective of the proposed research project was the analysis of single point closures based on probability density function (pdf) and characteristic functions and the development of a prediction method for the joint velocity-scalar pdf in turbulent reacting flows. Turbulent flows of boundary layer type and stagnation point flows with and without chemical reactions were be calculated as principal applications. Pdf methods for compressible reacting flows were developed and tested in comparison with available experimental data. The research work carried in this project was concentrated on the closure of pdf equations for incompressible and compressible turbulent flows with and without chemical reactions.

Kollmann, W.↗

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.↗

PDF turbulence modeling and DNS

The problem of time discontinuity (or jump condition) in the coalescence/dispersion (C/D) mixing model is addressed in probability density function (pdf). A C/D mixing model continuous in time is introduced. With the continuous mixing model, the process of chemical reaction can be fully coupled with mixing. In the case of homogeneous turbulence decay, the new model predicts a pdf very close to a Gaussian distribution, with finite higher moments also close to that of a Gaussian distribution. Results from the continuous mixing model are compared with both experimental data and numerical results from conventional C/D models. The effect of Coriolis forces on compressible homogeneous turbulence is studied using direct numerical simulation (DNS). The numerical method used in this study is an eight order compact difference scheme. Contrary to the conclusions reached by previous DNS studies on incompressible isotropic turbulence, the present results show that the Coriolis force increases the dissipation rate of turbulent kinetic energy, and that anisotropy develops as the Coriolis force increases. The Taylor-Proudman theory does apply since the derivatives in the direction of the rotation axis vanishes rapidly. A closer analysis reveals that the dissipation rate of the incompressible component of the turbulent kinetic energy indeed decreases with a higher rotation rate, consistent with incompressible flow simulations (Bardina), while the dissipation rate of the compressible part increases; the net gain is positive. Inertial waves are observed in the simulation results.

Hsu, A. T.↗

Modeling of turbulent supersonic H2-air combustion with an improved joint beta PDF

Attempts at modeling recent experiments of Cheng et al. indicated that discrepancies between theory and experiment can be a result of the form of assumed probability density function (PDF) and/or the turbulence model employed. Improvements in both the form of the assumed PDF and the turbulence model are presented. The results are again used to compare with measurements. Initial comparisons are encouraging.

Baurle, R. A.↗

Large eddy simulations and direct numerical simulations of high speed turbulent reacting flows

The primary objective of this research is to extend current capabilities of Large Eddy Simulations (LES) and Direct Numerical Simulations (DNS) for the computational analyses of high speed reacting flows. Our efforts in the first two years of this research have been concentrated on a priori investigations of single-point Probability Density Function (PDF) methods for providing subgrid closures in reacting turbulent flows. In the efforts initiated in the third year, our primary focus has been on performing actual LES by means of PDF methods. The approach is based on assumed PDF methods and we have performed extensive analysis of turbulent reacting flows by means of LES. This includes simulations of both three-dimensional (3D) isotropic compressible flows and two-dimensional reacting planar mixing layers. In addition to these LES analyses, some work is in progress to assess the extent of validity of our assumed PDF methods. This assessment is done by making detailed companions with recent laboratory data in predicting the rate of reactant conversion in parallel reacting shear flows. This report provides a summary of our achievements for the first six months of the third year of this program.

Givi, P.↗

Experimental study of cross flow mixing in cylindrical and rectangular ducts

An experimental investigation of non-reacting cross flow jet injection and mixing in cylindrical and rectangular ducts has been conducted with application to a low emissions combustor. Quantitative measurement of injectant concentration distributions perpendicular to the duct axis were obtained by planar digital imaging of the Mie-scattered light from an aerosol seed mixed with the injectant. The flowfield unmixedness was evaluated using (1) a mixing parameter derived from the ratio of the jet concentration fluctuations to the fully mixed concentration, and (2) probability density functions of the concentration distributions. Mixing rate was measured for 45 degree slant slot and round orifice injectors.

Liscinsky, D. S.↗

On the modeling of scalar diffusion in isotropic turbulence

The objectives of the study were to establish the behavior of conditional scalar dissipation and diffusion at the extreme values of mass fraction and to derive and evaluate closure models for these terms using the scalar probability density function. The conditional scalar dissipation, its derivative with respect to mass fraction, and conditional scalar diffusion are all found to be zero at extreme values of mass fraction. A model for conditional scalar dissipation is derived which exhibits the correct behavior at the extreme values of scalar concentration.

Girimaji, Sharath S.↗

Relationship between sea ice freeboard and draft in the Arctic Basin, and implications for ice thickness monitoring

This study confirms the finding of Comiso et al. (1991) that the probability density function (pdf) of the ice freeboard in the Arctic Ocean can be converted to a pdf of ice draft by applying a simple coordinate factor. The coordinate factor, R, which is the ratio of mean draft to mean freeboard pdf is related to the mean material (ice plus snow) density, rho(m), and the near-surface water density rho(w) by the relationship R = rho(m)/(rho(w) - rho(m)). The measured value of R was applied to each of six 50-km sections north of Greenland of a joint airborne laser and submarine sonar profile obtained along nearly coincident tracks from the Arctic Basin north of Greenland and was found to be consistent over all sections tested, despite differences in the ice regime. This indicates that a single value of R might be used for measurements done in this season of the year. The mean value R from all six sections was found to be 7.89.

Wadhams, P.↗

Application of a Reynolds stress turbulence model to a supersonic radiating hydrogen-air diffusion flame

A second-order differential Reynolds Stress turbulence model has been applied to the Favre-averaged Navier-Stokes equations for the study of supersonic flows with finite-rate chemistry and radiation. An assumed Beta Probability Density Function is applied to account for the chemical source terms and the radiative flux terms in the conservation equations. A seven-species, seven-reaction finite rate chemistry mechanism is used to simulate the combustion process. The tangent slab approximation is used in radiative flux formulation. A pseudo-gray gas model is used to represent the absorption-emission characteristics of the participating species. The turbulence/radiation interaction is achieved via a new formulation. The resulting formulation is validated by comparison with experimental data on reacting supersonic axisymmetric jets. Results obtained for specific conditions indicate that the effect of chemical reaction on the turbulence is significant. Also, the radiative heat transfer is enhanced by the turbulence.

Chandrasekhar, R.↗

The isolation limits of stochastic vibration

The vibration isolation problem is formulated as a 1D kinematic problem. The geometry of the stochastic wall trajectories arising from the stroke constraint is defined in terms of their significant extrema. An optimal control solution for the minimum acceleration return path determines a lower bound on platform mean square acceleration. This bound is expressed in terms of the probability density function on the significant maxima and the conditional fourth moment of the first passage time inverse. The first of these is found analytically while the second is found using a Monte Carlo simulation. The rms acceleration lower bound as a function of available space is then determined through numerical quadrature.

Knopse, C. R.↗