Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Average Accuracy”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Estimates of Long Term Surface Soil Moisture in the Midwestern U.S. Derived from Satellite Microwave Observations

Soil moisture is a key component of the water and energy balances of the Earth's surface, and has been identified as a parameter of significant potential for improving the accuracy of large-scale land surface-atmosphere interaction models. However, soil moisture is often somewhat difficult to measure accurately in both space and time, especially at large spatial scales. Soil moisture is highly variable, and while point measurements are typically quite accurate, subsequent areal averaging of these measurements often leads to large errors. Since remotely sensed land surface observations are already a spatially averaged or areally integrated value, they are a logical input parameter to regional or larger scale land process models. A database of long-term soil moisture was compared to satellite microwave observations over test sites in the Midwestern United States. Ground measurements of average volumetric surface soil moisture in the top ten cm were made bimonthly at 19 locations throughout the state of Illinois. Nighttime microwave brightness temperatures were observed at a frequency of 6.6 GHz, by the Scanning Multichannel Microwave Radiometer (SMMR), onboard the Nimbus 7 satellite. The life of the SMMR instrument spanned from Nov. 1978 to Aug. 1987. At 6.6 GHz, the instrument provided a spatial resolution of approximately 150 km, and a temporal frequency over the test area of about 3 nighttime orbits per week. Vegetation radiative transfer characteristics, such as the canopy transmissivity, were estimated from vegetation indices such as the Normalized Difference Vegetation Index (NDVI) and the 37 GHz Microwave Polarization Difference Index (MPDI). Because the time of satellite coverage does not always coincide with the ground measurements of soil moisture, the existing ground data were used to calibrate a water balance for the top IO cm surface layer in order to interpolate daily surface moisture values. Such a climate-based approach is often more appropriate for estimating large-area average soil moisture because meteorological data are generally more spatially representative than isolated point measurements of soil moisture, Passive microwave remote sensing presents the greatest potential for providing regular spatially representative estimates of surface soil moisture at global scales. Real time estimates should improve weather and climate modelling efforts, while the development of historical data sets will provide necessary information for simulation and validation of long-term climate and global change studies.

Owe, M.↗
NASA NTRS1999

Practical Aerodynamic Design Optimization Based on the Navier-Stokes Equations and a Discrete Adjoint Method

Compressible and incompressible versions of a three-dimensional unstructured mesh Reynolds-averaged Navier-Stokes flow solver have been differentiated and resulting derivatives have been verified by comparisons with finite differences and a complex-variable approach. In this implementation, the turbulence model is fully coupled with the flow equations in order to achieve this consistency. The accuracy demonstrated in the current work represents the first time that such an approach has been successfully implemented. The accuracy of a number of simplifying approximations to the linearizations of the residual have been examined. A first-order approximation to the dependent variables in both the adjoint and design equations has been investigated. The effects of a "frozen" eddy viscosity and the ramifications of neglecting some mesh sensitivity terms were also examined. It has been found that none of the approximations yielded derivatives of acceptable accuracy and were often of incorrect sign. However, numerical experiments indicate that an incomplete convergence of the adjoint system often yield sufficiently accurate derivatives, thereby significantly lowering the time required for computing sensitivity information. The convergence rate of the adjoint solver relative to the flow solver has been examined. Inviscid adjoint solutions typically require one to four times the cost of a flow solution, while for turbulent adjoint computations, this ratio can reach as high as eight to ten. Numerical experiments have shown that the adjoint solver can stall before converging the solution to machine accuracy, particularly for viscous cases. A possible remedy for this phenomenon would be to include the complete higher-order linearization in the preconditioning step, or to employ a simple form of mesh sequencing to obtain better approximations to the solution through the use of coarser meshes. An efficient surface parameterization based on a free-form deformation technique has been utilized and the resulting codes have been integrated with an optimization package. Lastly, sample optimizations have been shown for inviscid and turbulent flow over an ONERA M6 wing. Drag reductions have been demonstrated by reducing shock strengths across the span of the wing. In order for large scale optimization to become routine, the benefits of parallel architectures should be exploited. Although the flow solver has been parallelized using compiler directives. The parallel efficiency is under 50 percent. Clearly, parallel versions of the codes will have an immediate impact on the ability to design realistic configurations on fine meshes, and this effort is currently underway.

Grossman, Bernard↗
NASA NTRS1999

Practical Aerodynamic Design Optimization Based on the Navier-Stokes Equations and a Discrete Adjoint Method

The technical details are summarized below: Compressible and incompressible versions of a three-dimensional unstructured mesh Reynolds-averaged Navier-Stokes flow solver have been differentiated and resulting derivatives have been verified by comparisons with finite differences and a complex-variable approach. In this implementation, the turbulence model is fully coupled with the flow equations in order to achieve this consistency. The accuracy demonstrated in the current work represents the first time that such an approach has been successfully implemented. The accuracy of a number of simplifying approximations to the linearizations of the residual have been examined. A first-order approximation to the dependent variables in both the adjoint and design equations has been investigated. The effects of a "frozen" eddy viscosity and the ramifications of neglecting some mesh sensitivity terms were also examined. It has been found that none of the approximations yielded derivatives of acceptable accuracy and were often of incorrect sign. However, numerical experiments indicate that an incomplete convergence of the adjoint system often yield sufficiently accurate derivatives, thereby significantly lowering the time required for computing sensitivity information. The convergence rate of the adjoint solver relative to the flow solver has been examined. Inviscid adjoint solutions typically require one to four times the cost of a flow solution, while for turbulent adjoint computations, this ratio can reach as high as eight to ten. Numerical experiments have shown that the adjoint solver can stall before converging the solution to machine accuracy, particularly for viscous cases. A possible remedy for this phenomenon would be to include the complete higher-order linearization in the preconditioning step, or to employ a simple form of mesh sequencing to obtain better approximations to the solution through the use of coarser meshes. . An efficient surface parameterization based on a free-form deformation technique has been utilized and the resulting codes have been integrated with an optimization package. Lastly, sample optimizations have been shown for inviscid and turbulent flow over an ONERA M6 wing. Drag reductions have been demonstrated by reducing shock strengths across the span of the wing.

Grossman, Bernard↗
NASA NTRS2000

SHDOM

SHDOM is a general purpose, publicly available, three-dimensional atmospheric radiative transfer model. SHDOM is an explicit method, which means it solves for the whole radiation field, as distinct from Monte Carlo methods which solve for particular radiative outputs. SHDOM is particularly well suited for remote sensing applications, where it can compute outgoing radiances at many angles from a cloud field at virtually no extra cost. SHDOM is not appropriate for calculating domain average quantities for which Monte Carlo methods excel. The I3RC intercomparison offers an opportunity to explore the pros and cons of SHDOM and Monte Carlo models on some real world inhomogeneous cloud fields. Specifically, we wish to determine the computer resources required to achieve a particular accuracy for a certain number of outputs using SHDOM and Monte Carlo models. This will help guide modelers on the appropriate choice of SHDOM or Monte Carlo for their applications. To emphasize the importance of this accuracy versus CPU time tradeoff, we are submitting two SHDOM entries (low and high resolution) in the I3RC.

Evans, K. Franklin↗
NASA NTRS2010

A New Approach for Checking and Complementing CALIPSO Lidar Calibration

We have been studying the backscatter ratio of the two CALIPSO wavelengths for 3 different targets. We are showing the ratio of integrate attenuated backscatter coefficient for cirrus clouds, ocean surface and liquid. Water clouds for one month of nightime data (left:July,right:December), Only opaque cirrus classified as randomly oriented ice[1] are used. For ocean and water clouds, only the clearest shots, determined by a threshold on integrated attenuated backscatter are used. Two things can be immediately observed: 1. A similar trend (black dotted line) is visible using all targets, the color ratio shows a tendency to be higher north and lower south for those two months. 2. The water clouds average value is around 15% lower than ocean surface and cirrus clouds. This is due to the different multiple scattering at 532 nm and 1064 nm [2] which strongly impact the water cloud retrieval. Conclusion: Different targets can be used to improve CALIPSO 1064 nm calibration accuracy. All of them show the signature of an instrumental calibration shift. Multiple scattering introduce a bias in liquid water cloud signal but it still compares very well with all other methods and should not be overlooked. The effect of multiple scattering in liquid and ice clouds will be the subject of future research. If there really is a sampling issue. Combining all methods to increase the sampling, mapping the calibration coefficient or trying to reach an orbit per orbit calibration seems an appropriate way.

Josset, Damien B.↗
NASA NTRS2018

Convergence Analysis of Turbulent Flow Solutions

Data from the "Turbulence Modeling Resource" website for turbulent flow over an NACA-0012 airfoil is analyzed to determine the convergence behavior of three second-order CFD (Computational Fluid Dynamics) codes: CFL3D (Computational Fluids Lab 3 Dimensional flow solver), FUN3D (Fully Unstructured Navier-stokes flow solver), and TAU (German Aerospace Center (DLR) 2 dimensional code for unstructured hybrid grids solving the Reynolds-Averaged Navier-Stokes equations or the Euler equations). The convergence of both integrated properties and pointwise data are examined. Several different methods for estimating errors and computing convergence rates are compared. A high-order extension to the Richardson extrapolation is developed that improves the accuracy of the mesh limit values and provides a quantitative estimate of the threshold of the asymptotic regime. The coefficient of total drag exhibits second-order convergence for all three codes, and convergence is monotone over a sequence of 7 grids. Other force coefficients are not so well behaved. The convergence rates of the viscous component of drag on the three nest grids ranges from 3:0 for CFL3D to 1:0 for FUN3D. The three codes are converging to similar but not identical solutions. The largest differences between the codes are in the coefficient of lift for which the difference between CFL3D and FUN3D is greater than 10 (sup minus 4). The best agreement occurs in the viscous component of drag, which is the only force component for which all three codes are converging towards each other at a rate of second-order. The agreement between the two unstructured grid codes is good with all properties except lift converging towards common values at a rate of second-order. No one code was universally better than the other. The TAU code has the lowest error in total drag, FUN3D has the lowest error in lift, and CFL3D has the lowest error in the viscous component of drag.

Atkins, Harold L.↗
NASA NTRS2019

Global Validation of Columnar Water Vapor Derived from EOS MODIS-MAIAC Algorithm Against the Ground-Based AERONET Observations

The water vapor is a relevant greenhouse gas in the Earth's climate system, and satellite products become one of the most effective way to characterize and monitor the columnar water vapor (CWV) content at global scale. Recently, a new product (MCD19) was released as part of MODIS (Moderate Resolution Imaging Spectroradiometer) Collection 6 (C6). This operational product from the Multi-Angle Implementation for Atmospheric Correction (MAIAC) algorithm includes a high 1-kilometer resolution CWV retrievals. This study presents the first global validation of MAIAC C6 CWV obtained from MODIS MCD19A2 product. This evaluation was performed using Aerosol Robotic Network (AERONET) observations at 265 sites (2000-2017). Overall, the results show a good agreement between MAIAC/AERONET CWV retrievals, with correlation coefficient higher than 0.95 and RMS (Root Mean Square) error lower than 0.250 centimeters. The binned error analysis revealed an underestimation (approximately 10 percent) of Aqua CWV retrievals with negative bias for CWV higher than 3.0 centimeters. In contrast, Terra CWV retrievals show a slope of regression close to unity and a low mean bias of 0.075 centimeters. While the accuracy is relatively similar between 1.0 and 5.0 centimeters for both sensor products, Terra dataset is more reliable for applications in humid tropical areas (less than 5.0 centimeters). The expected error was defined as plus or minus 15 percent, with less than 68 percent of retrievals falling within this envelope. However, the accuracy is regionally dependent, and lower error should be expected in some regions, such as South America and Oceania. Since MODIS instruments have exceeded their design lifetime, time series analysis was also presented for both sensor products. The temporal analysis revealed a systematic offset of global average between Terra and Aqua CWV records. We also found an upward trend (approximately 0.2 centimeters per decade) in Terra CWV retrievals, while Aqua CWV retrievals remain stable over time. The sensor degradation influences the ability to detect climate signals, and this study indicates the need for revisiting calibration of the MODIS bands 17-19, mainly for Terra instrument, to assure the quality of the MODIS water vapor product. Finally, this study presents a comprehensive validation analysis of MAIAC CWV over land, raising the understanding of its overall quality.

MAIAC↗
NASA NTRS2018

Application of Lattice Boltzmann and Navier-Stokes Methods to NASA's Wall Mounted Hump

Lattice Boltzmann (LB) based Large Eddy Simulation (LES), Reynolds-averaged Navier-Stokes (RANS) as well as hybrid RANS/LES methods within the Launch Ascent and Vehicle Aerodynamics (LAVA) solver framework are applied to NASA's wall-mounted hump. Computational results are compared with experiments performed by Greenblatt et al. A detailed comparison between the accuracy and resolution requirements of the two approaches for turbulence resolving simulations, as well as the suitability of different grid paradigms (body-fitted curvilinear and block structured Cartesian) are presented. This test case is part of NASA's Revolutionary Computational Aerosciences (RCA) sub-project which addresses the technical challenge of predicting flow separation and reattachment accurately. Improvements in predictive accuracy by as much as 90% are demonstrated using LB as well as hybrid RANS/LES approaches compared to state-of-the-art steady state RANS simulations.

Lattice Boltzmann↗
NASA NTRS2021

Florida Transportation and Infrastructure: Monitoring Water Quality Along Southern Florida Seaports to Assess Impact on Coral Reef Tracts from Harbor Deepening Projects

The U.S. Army Corps of Engineers (USACE) and National Oceanic and Atmospheric Administration (NOAA) National Marine Fisheries Service (NMFS) will be supervising a harbor deepening project in Port Everglades, Florida. The project raises concerns about potential impacts on the nearby Florida reef tract through increased turbidity and sediment from the dredging. To better understand these potential impacts, the NASA DEVELOP team created an interactive Google Earth Engine tool to help establish a historical baseline of water quality parameters and assist monitoring these parameters more frequently than traditional sampling. This Seaport & Harbor Area Resource Quality (SHARQ) tool incorporates remotely sensed data from Sentinel-2 Multispectral Instrument, Landsat 5 Thematic Mapper, Landsat 7 Enhanced Thematic Mapper+, Landsat 8 Operational Land Imager, and Aqua Moderate Resolution Imaging Spectroradiometer. It allows users to view true color images and calculate water quality parameters, like turbidity and chlorophyll-a, for any given study area and time period from 1984 onward. The SHARQ tool also generate time series charts, allowing users to interpret changes in water quality over a given time range. The accuracy of the remotely sensed water quality parameter algorithms was determined using in situ data to calculate percent difference and root mean square error values (RMSE), which ranged from 0.32 to 0.58 error between sites. Using the SHARQ tool’s time series analysis feature, a baseline average turbidity metric of ~6.8 FNUs provides a historical baseline average for turbidity between September 2000 and 2020 and can assist in future decision-making for determining thresholds for turbidity.

Benjamin Rocha↗
NASA NTRS1979

Present obliquity oscillations of Mars - Fourth-order accuracy in orbital e and I

A long period analysis of solar system orbital evolution, correct to fourth order in orbital eccentricities and inclinations (Bretagnon, 1974), and an improved value of the planet's moment of inertia (Reasenberg, 1977) have been incorporated in a recalculation of the obliquity oscillations of Mars. A linearized solution predicts a maximum oscillation amplitude of 13.6 deg centered on a long-term average value of 24.4 deg. A numerical integration of the obliquity for the past 10,000,000 years is also presented. Epochs of minimal oscillation like the present occur at intervals of the order of 4 m.y.

Ward, W. R.↗
NASA NTRS1980

Initial geomagnetic field model from Magsat vector data

Magsat data from the magnetically quiet days of November 5-6, 1979, were used to derive a thirteenth degree and order spherical harmonic geomagnetic field model, MGST(6/80). The model utilized both scalar and high-accuracy vector data and fit that data with root-mean-square deviations of 8.2, 6.9, 7.6 and 7.4 nT for the scalar magnitude, B(r), B(theta), and B(phi), respectively. The model includes the three first-order coefficients of the external field. Comparison with averaged Dst indicates that zero Dst corresponds with 25 nT of horizontal field from external sources. When compared with earlier models, the earth's dipole moment continues to decrease at a rate of about 26 nT/yr. Evaluation of earlier models with Magsat data shows that the scalar field at the Magsat epoch is best predicted by the POGO(2/72) model but that the WC80, AWC/75 and IGS/75 are better for predicting vector fields.

Langel, R. A.↗
NASA NTRS1984

Determination of monthly mean humidity in the atmospheric surface layer over oceans from satellite data

A simple statistical technique is described to determine monthly mean marine surface-layer humidity, which is essential in the specification of surface latent heat flux, from total water vapor in the atmospheric column measured by space-borne sensors. Good correlation between the two quantities was found in examining the humidity soundings from radiosonde reports of mid-ocean island stations and weather ships. The relation agrees with that obtained from satellite (Seasat) data and ship reports averaged over 2 deg areas and a 92-day period in the North Atlantic and in the tropical Pacific. The results demonstrate that, by using a local regression in the tropical Pacific, total water vapor can be used to determine monthly mean surface layer humidity to an accuracy of 0.4 g/kg. With a global regression, determination to an accuracy of 0.8 g/kg is possible. These accuracies correspond to approximately 10 to 20 W/sq m in the determination of latent heat flux with the bulk parameterization method, provided that other required parameters are known.

Liu, W. T.↗
NASA NTRS1990

Multiprocessing on supercomputers for computational aerodynamics

Very little use is made of multiple processors available on current supercomputers (computers with a theoretical peak performance capability equal to 100 MFLOPs or more) in computational aerodynamics to significantly improve turnaround time. The productivity of a computer user is directly related to this turnaround time. In a time-sharing environment, the improvement in this speed is achieved when multiple processors are used efficiently to execute an algorithm. The concept of multiple instructions and multiple data (MIMD) through multi-tasking is applied via a strategy which requires relatively minor modifications to an existing code for a single processor. Essentially, this approach maps the available memory to multiple processors, exploiting the C-FORTRAN-Unix interface. The existing single processor code is mapped without the need for developing a new algorithm. The procedure for building a code utilizing this approach is automated with the Unix stream editor. As a demonstration of this approach, a Multiple Processor Multiple Grid (MPMG) code is developed. It is capable of using nine processors, and can be easily extended to a larger number of processors. This code solves the three-dimensional, Reynolds averaged, thin-layer and slender-layer Navier-Stokes equations with an implicit, approximately factored and diagonalized method. The solver is applied to generic oblique-wing aircraft problem on a four processor Cray-2 computer. A tricubic interpolation scheme is developed to increase the accuracy of coupling of overlapped grids. For the oblique-wing aircraft problem, a speedup of two in elapsed (turnaround) time is observed in a saturated time-sharing environment.

Yarrow, Maurice↗
NASA NTRS2010

Satellite Sampling and Retrieval Errors in Regional Monthly Rain Estimates from TMI AMSR-E, SSM/I, AMSU-B and the TRMM PR

Passive and active microwave rain sensors onboard earth-orbiting satellites estimate monthly rainfall from the instantaneous rain statistics collected during satellite overpasses. It is well known that climate-scale rain estimates from meteorological satellites incur sampling errors resulting from the process of discrete temporal sampling and statistical averaging. Sampling and retrieval errors ultimately become entangled in the estimation of the mean monthly rain rate. The sampling component of the error budget effectively introduces statistical noise into climate-scale rain estimates that obscure the error component associated with the instantaneous rain retrieval. Estimating the accuracy of the retrievals on monthly scales therefore necessitates a decomposition of the total error budget into sampling and retrieval error quantities. This paper presents results from a statistical evaluation of the sampling and retrieval errors for five different space-borne rain sensors on board nine orbiting satellites. Using an error decomposition methodology developed by one of the authors, sampling and retrieval errors were estimated at 0.25 resolution within 150 km of ground-based weather radars located at Kwajalein, Marshall Islands and Melbourne, Florida. Error and bias statistics were calculated according to the land, ocean and coast classifications of the surface terrain mask developed for the Goddard Profiling (GPROF) rain algorithm. Variations in the comparative error statistics are attributed to various factors related to differences in the swath geometry of each rain sensor, the orbital and instrument characteristics of the satellite and the regional climatology. The most significant result from this study found that each of the satellites incurred negative longterm oceanic retrieval biases of 10 to 30%.

Fisher, Brad↗
NASA NTRS2020

Tuning Neural Network Models for Improved Prediction of Boundary Layer Transition

Boundary layer transition can strongly impact flight vehicle performance as it influences surface skin friction and aerodynamic heating, making accurate transition prediction a key to designing next generation aircraft. Artificial neural networks (ANNs) have shown promise toward predicting laminar-turbulent transition based on linear stability correlations. The computational efficiency of ANNs and the substantially reduced user involvement in relation to direct computations based on the linear stability theory (LST) makes them an attractive methodology for integrating the LST based correlations in computational fluid dynamics codes. Tollmien-Schlichting (TS) waves correspond to the dominant transition mechanism in 2D or weakly 3D subsonic boundary layers, such as those encountered in general aviation applications. Improvements to neural network model accuracy in predicting the amplification rates of TS instability waves have been investigated by leveraging recent machine learning developments in conjunction with surrogate optimization techniques and via suitable augmentation of the data used to train the networks. The optimized models trained on the modified dataset reduced the average transition location errors on different airfoils at several flow conditions by 51% of the original manually-tuned network’s errors on the same flow cases. The actual transition locations were derived from the Langley Stability and Transition Analysis Code (LASTRAC).

Machine Learning↗
NASA NTRS1994

An improved method of measuring tropospheric NO2, NO3, HO2, and RO2 by Matrix Isolation and Electronic Spin Resonance (MIESR)

The MIESR method consists of two steps (1) collection of the radicals present in the ambient air at 77K in a polycrystalline D2O matrix and (2) identification and quantification of the different radicals in the laboratory by Electron Spin Resonance spectroscopy. In step (1), the sampling efficiency for sampling NO2 and RO2 was determined to be greater than or equal to 95 percent, with a measured accuracy of plus or minus 5 percent. In step (2), after collection, the samples are maintained at 77 K and spectra are recorded in the laboratory using a standard 9.5 GHz ESR system (Varian E-line). About 50 individual scans of each spectrum are recorded and digitally averaged in order to improve the signal-to-noise ratio. The ESR-spectra are analyzed with a recently developed numerical procedure which was demonstrated to allow speciation of NO2, NO3, HO2, CH3C(O)O2, and the sum of the alkylperoxy radicals. The detection limit is 5ppt for HO2, RO2, and NO2 and 3ppt for NO3 due to its narrower ESR-linewidth.

Source record↗
NASA NTRS1998

Application of Pressure Sensitive Paint to Confined Flow at Mach Number 2.5

Pressure sensitive paint (PSP) is a novel technology that is being used frequently in external aerodynamics. For internal flows in narrow channels, and applications at elevated nonuniform temperatures, however, there are still unresolved problems that complicate the procedures for calibrating PSP signals. To address some of these problems, investigations were carried out in a narrow channel with supersonic flows of Mach 2.5. The first set of tests focused on the distribution of the wall pressure in the diverging section of the test channel downstream of the nozzle throat. The second set dealt with the distribution of wall static pressure due to the shock/wall interaction caused by a 25 deg. wedge in the constant Mach number part of the test section. In addition, the total temperature of the flow was varied to assess the effects of temperature on the PSP signal. Finally, contamination of the pressure field data, caused by internal reflection of the PSP signal in a narrow channel, was demonstrated. The local wall pressures were measured with static taps, and the wall pressure distributions were acquired by using PSP. The PSP results gave excellent qualitative impressions of the pressure field investigated. However, the quantitative results, specifically the accuracy of the PSP data in narrow channels, show that improvements need to be made in the calibration procedures, particularly for heated flows. In the cases investigated, the experimental error had a standard deviation of +/- 8.0% for the unheated flow, and +/- 16.0% for the heated flow, at an average pressure of 11 kpa.

Lepicovsky, J.↗
NASA NTRS1998

Improved Modeling of Finite-Rate Turbulent Combustion Processes in Research Combustors

The objective of this thesis is to further develop and test a stochastic model of turbulent combustion in recirculating flows. There is a requirement to increase the accuracy of multi-dimensional combustion predictions. As turbulence affects reaction rates, this interaction must be more accurately evaluated. In this work a more physically correct way of handling the interaction of turbulence on combustion is further developed and tested. As turbulence involves randomness, stochastic modeling is used. Averaged values such as temperature and species concentration are found by integrating the probability density function (pdf) over the range of the scalar. The model in this work does not assume the pdf type, but solves for the evolution of the pdf using the Monte Carlo solution technique. The model is further developed by including a more robust reaction solver, by using accurate thermodynamics and by more accurate transport elements. The stochastic method is used with Semi-Implicit Method for Pressure-Linked Equations. The SIMPLE method is used to solve for velocity, pressure, turbulent kinetic energy and dissipation. The pdf solver solves for temperature and species concentration. Thus, the method is partially familiar to combustor engineers. The method is compared to benchmark experimental data and baseline calculations. The baseline method was tested on isothermal flows, evaporating sprays and combusting sprays. Pdf and baseline predictions were performed for three diffusion flames and one premixed flame. The pdf method predicted lower combustion rates than the baseline method in agreement with the data, except for the premixed flame. The baseline and stochastic predictions bounded the experimental data for the premixed flame. The use of a continuous mixing model or relax to mean mixing model had little effect on the prediction of average temperature. Two grids were used in a hydrogen diffusion flame simulation. Grid density did not effect the predictions except for peak temperature and tangential velocity. The hybrid pdf method did take longer and required more memory, but has a theoretical basis to extend to many reaction steps which cannot be said of current turbulent combustion models.

VanOverbeke, Thomas J.↗