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 577 records · Page 32

Small-Scale Gravity Waves in ER-2 MMS/MTP Wind and Temperature Measurements during CRYSTAL-FACE

Lower stratospheric wind and temperature measurements made from NASA's high-altitude ER-2 research aircraft during the CRYSTAL-FACE campaign in July 2002 were analyzed to retrieve information on small scale gravity waves (GWs) at the aircraft's flight level (typically approximately 20 km altitude). For a given flight segment, the S-transform (a Gaussian wavelet transform) was used to search for and identify small horizontal scale GW events, and to estimate their apparent horizontal wavelengths. The horizontal propagation directions of the events were determined using the Stokes parameter method combined with the cross S-transform analysis. The vertical temperature gradient was used to determine the vertical wavelengths of the events. GW momentum fluxes were calculated from the cross S-transform. Other wave parameters such as intrinsic frequencies were calculated using the GW dispersion relation. More than 100GW events were identified. They were generally high frequency waves with vertical wavelength of approximately 5 km and horizontal wavelength generally shorter than 20 km. Their intrinsic propagation directions were predominantly toward the east, whereas their ground-based propagation directions were primarily toward the west. Among the events, approximately 20% of them had very short horizontal wavelength, very high intrinsic frequency, and relatively small momentum fluxes, and thus they were likely trapped in the lower stratosphere. Using the estimated GW parameters and the background winds and stabilities from the NCAR/NCEP reanalysis data, we were able to trace the sources of the events using a simple reverse ray-tracing. More than 70% of the events were traced back to convective sources in the troposphere, and the sources were generally located upstream of the locations of the events observed at the aircraft level. Finally, a probability density function of the reversible cooling rate due to GWs was obtained in this study, which may be useful for cirrus cloud models.

Wang, L.↗

Single-Shot Scalar-Triplet Measurements in High-Pressure Swirl-Stabilized Flames for Combustion Code Validation

In support of NASA ARMD's code validation project, we have made significant progress by providing the first quantitative single-shot multi-scalar data from a turbulent elevated-pressure (5 atm), swirl-stabilized, lean direct injection (LDI) type research burner operating on CH4-air using a spatially-resolved pulsed-laser spontaneous Raman diagnostic technique. The Raman diagnostics apparatus and data analysis that we present here were developed over the past 6 years at Glenn Research Center. From the Raman scattering data, we produce spatially-mapped probability density functions (PDFs) of the instantaneous temperature, determined using a newly developed low-resolution effective rotational bandwidth (ERB) technique. The measured 3-scalar (triplet) correlations, between temperature, CH4, and O2 concentrations, as well as their PDF s, also provide a high-level of detail into the nature and extent of the turbulent mixing process and its impact on chemical reactions in a realistic gas turbine injector flame at elevated pressures. The multi-scalar triplet data set presented here provides a good validation case for CFD combustion codes to simulate by providing both average and statistical values for the 3 measured scalars.

Kojima, Jun↗

Simulation of Shallow Cumuli and Their Transition to Deep Convective Clouds by Cloud-resolving Models with Different Third-order Turbulence Closures

The abilities of cloud-resolving models (CRMs) with the double-Gaussian based and the single-Gaussian based third-order closures (TOCs) to simulate the shallow cumuli and their transition to deep convective clouds are compared in this study. The single-Gaussian based TOC is fully prognostic (FP), while the double-Gaussian based TOC is partially prognostic (PP). The latter only predicts three important third-order moments while the former predicts all the thirdorder moments. A shallow cumulus case is simulated by single-column versions of the FP and PP TOC models. The PP TOC improves the simulation of shallow cumulus greatly over the FP TOC by producing more realistic cloud structures. Large differences between the FP and PP TOC simulations appear in the cloud layer of the second- and third-order moments, which are related mainly to the underestimate of the cloud height in the FP TOC simulation. Sensitivity experiments and analysis of probability density functions (PDFs) used in the TOCs show that both the turbulence-scale condensation and higher-order moments are important to realistic simulations of the boundary-layer shallow cumuli. A shallow to deep convective cloud transition case is also simulated by the 2-D versions of the FP and PP TOC models. Both CRMs can capture the transition from the shallow cumuli to deep convective clouds. The PP simulations produce more and deeper shallow cumuli than the FP simulations, but the FP simulations produce larger and wider convective clouds than the PP simulations. The temporal evolutions of cloud and precipitation are closely related to the turbulent transport, the cold pool and the cloud-scale circulation. The large amount of turbulent mixing associated with the shallow cumuli slows down the increase of the convective available potential energy and inhibits the early transition to deep convective clouds in the PP simulation. When the deep convective clouds fully develop and the precipitation is produced, the cold pools produced by the evaporation of the precipitation are not favorable to the formation of shallow cumuli.

Cheng, Anning↗

Non-Gaussian Analysis of Turbulent Boundary Layer Fluctuating Pressure on Aircraft Skin Panels

The purpose of the study is to investigate the probability density function (PDF) of turbulent boundary layer fluctuating pressures measured on the outer sidewall of a supersonic transport aircraft and to approximate these PDFs by analytical models. Experimental flight results show that the fluctuating pressure PDFs differ from the Gaussian distribution even for standard smooth surface conditions. The PDF tails are wider and longer than those of the Gaussian model. For pressure fluctuations in front of forward-facing step discontinuities, deviations from the Gaussian model are more significant and the PDFs become asymmetrical. There is a certain spatial pattern of the skewness and kurtosis behavior depending on the distance upstream from the step. All characteristics related to non-Gaussian behavior are highly dependent upon the distance from the step and the step height, less dependent on aircraft speed, and not dependent on the fuselage location. A Hermite polynomial transform model and a piecewise-Gaussian model fit the flight data well both for the smooth and stepped conditions. The piecewise-Gaussian approximation can be additionally regarded for convenience in usage after the model is constructed.

Rizzi, Stephen A.↗

LSPRAY-III: A Lagrangian Spray Module

LSPRAY-III is a Lagrangian spray solver developed for application with parallel computing and unstructured grids. It is designed to be massively parallel and could easily be coupled with any existing gas-phase flow and/or Monte Carlo Probability Density Function (PDF) solvers. The solver accommodates the use of an unstructured mesh with mixed elements of either triangular, quadrilateral, and/or tetrahedral type for the gas flow grid representation. It is mainly designed to predict the flow, thermal and transport properties of a rapidly vaporizing spray because of its importance in aerospace application. The manual provides the user with an understanding of various models involved in the spray formulation, its code structure and solution algorithm, and various other issues related to parallelization and its coupling with other solvers. With the development of LSPRAY-III, we have advanced the state-of-the-art in spray computations in several important ways.

Raju, M. S.↗

Evaluation of Cloud Physical Properties of ECMWF Analysis and Re-Analysis (ERA-40 and ERA Interim) against CERES Tropical Deep Convective Cloud Object Observations

This study presents an approach that converts the vertical profiles of grid-averaged cloud properties from large-scale models to probability density functions (pdfs) of subgrid-cell cloud physical properties measured at satellite footprints. Cloud physical and radiative properties, rather than just cloud and precipitation occurrences, of assimilated cloud systems by the European Center for Medium-range Weather Forecasts (ECMWF) operational analysis (EOA) and ECMWF Re-Analyses (ERA-40 and ERA Interim) are validated against those obtained from Earth Observing System satellite cloud object data for January-August 1998 and March 2000 periods. These properties include ice water path (IWP), cloud-top height and temperature, cloud optical depth and solar and infrared radiative fluxes. Each cloud object, a contiguous region with similar cloud physical properties, is temporally and spatially matched with EOA and ERA-40 data. Results indicate that most pdfs of EOA and ERA-40 cloud physical and radiative properties agree with those of satellite observations of the tropical deep convective cloud-object type for the January-August 1998 period. There are, however, significant discrepancies in selected ranges of the cloud property pdfs such as the upper range of EOA cloud top height. A major discrepancy is that the dependence of the pdfs on the cloud object size for both EOA and ERA-40 is not as strong as in the observations. Modifications to the cloud parameterization in ECMWF that occurred in October 1999 eliminate the clouds near the tropopause but shift power of the pdf to lower cloud-top heights and greatly reduce the ranges of IWP and cloud optical depth pdfs. These features persist in ERA-40 due to the use of the same cloud parameterizations. The downgrade of data assimilation technique and the lack of snow water content information in ERA-40, not the coarser horizontal grid resolution, are also responsible for the disagreements with observed pdfs of cloud physical properties although the detection rates of cloud object occurrence are improved for small size categories. A possible improvement to the convective parameterization is to introduce a stronger dependence of updraft penetration heights with grid-cell dynamics. These conclusions will be rechecked using the ERA Interim data, due to recent changes in the ECMWF convective parameterization (Bechtold et al. 2004, 2008). Results from the ERA Interim will be presented at the meeting.

Xu, Kuan-Man↗

Evaluation of a New Cloud Droplet Activation Parameterization wtih in Situ Data from CRYSTAL-FACE and CSTRIPE

The accuracy of the 2003 prognostic, physically based aerosol activation parameterization of A. Nenes and J. H. Seinfeld (NS) with modification introduced by C. Fountoukis and A. Nenes in 2005 (modified NS) is evaluated against extensive microphysical data sets collected on board the Center for Interdisciplinary Remotely Piloted Aircraft Studies (CIRPAS) Twin Otter aircraft for cumuliform and stratiform clouds of marine and continental origin. The cumuliform cloud data were collected during NASA's Cirrus Regional Study of Tropical Anvils and Cirrus Layers-Florida Area Cirrus Experiment (CRYSTAL-FACE, Key West, Florida, July 2002), while the stratiform cloud data were gathered during Coastal Stratocumulus Imposed Perturbation Experiment (CSTRIPE, Monterey, California, July 2003). In situ data sets of aerosol size distribution, chemical composition, and updraft velocities are used as input for the NS parameterization, and the evaluation is carried out by comparing predicted cloud droplet number concentrations (CDNC) with observations. This is the first known study in which a prognostic cloud droplet activation parameterization has been evaluated against a wide range of observations. On average, predicted droplet concentration in adiabatic regions is within -20% of observations at the base of cumuliform clouds and -30% of observations at different altitudes throughout the stratiform clouds, all within experimental uncertainty. Furthermore, CDNC is well parameterized using either a single mean updraft velocity w or by weighting droplet nucleation rates with a Gaussian probability density function of w. This study suggests that for nonprecipitating warm clouds of variable microphysics, aerosol composition, and size distribution the modified NS parameterization can accurately predict cloud droplet activation and can be successfully implemented for describing the aerosol activation process in global climate models.

Meskhidze, Nicholas↗

Simulation of Boundary-Layer Cumulus and Stratocumulus Clouds using a Cloud-Resolving Model With Low- and Third-Order Turbulence Closures

The effects of subgrid-scale condensation and transport become more important as the grid spacings increase from those typically used in large-eddy simulation (LES) to those typically used in cloud-resolving models (CRMs). Incorporation of these effects can be achieved by a joint probability density function approach that utilizes higher-order moments of thermodynamic and dynamic variables. This study examines how well shallow cumulus and stratocumulus clouds are simulated by two versions of a CRM that is implemented with low-order and third-order turbulence closures (LOC and TOC) when a typical CRM horizontal resolution is used and what roles the subgrid-scale and resolved-scale processes play as the horizontal grid spacing of the CRM becomes finer. Cumulus clouds were mostly produced through subgrid-scale transport processes while stratocumulus clouds were produced through both subgrid-scale and resolved-scale processes in the TOC version of the CRM when a typical CRM grid spacing is used. The LOC version of the CRM relied upon resolved-scale circulations to produce both cumulus and stratocumulus clouds, due to small subgrid-scale transports. The mean profiles of thermodynamic variables, cloud fraction and liquid water content exhibit significant differences between the two versions of the CRM, with the TOC results agreeing better with the LES than the LOC results. The characteristics, temporal evolution and mean profiles of shallow cumulus and stratocumulus clouds are weakly dependent upon the horizontal grid spacing used in the TOC CRM. However, the ratio of the subgrid-scale to resolved-scale fluxes becomes smaller as the horizontal grid spacing decreases. The subcloud-layer fluxes are mostly due to the resolved scales when a grid spacing less than or equal to 1 km is used. The overall results of the TOC simulations suggest that a 1-km grid spacing is a good choice for CRM simulation of shallow cumulus and stratocumulus.

Xu, Kuan-Man↗

NASA Tech Briefs, December 2009

Topics include: A Deep Space Network Portable Radio Science Receiver; Detecting Phase Boundaries in Hard-Sphere Suspensions; Low-Complexity Lossless and Near-Lossless Data Compression Technique for Multispectral Imagery; Very-Long-Distance Remote Hearing and Vibrometry; Using GPS to Detect Imminent Tsunamis; Stream Flow Prediction by Remote Sensing and Genetic Programming; Pilotless Frame Synchronization Using LDPC Code Constraints; Radiometer on a Chip; Measuring Luminescence Lifetime With Help of a DSP; Modulation Based on Probability Density Functions; Ku Telemetry Modulator for Suborbital Vehicles; Photonic Links for High-Performance Arraying of Antennas; Reconfigurable, Bi-Directional Flexfet Level Shifter for Low-Power, Rad-Hard Integration; Hardware-Efficient Monitoring of I/O Signals; Video System for Viewing From a Remote or Windowless Cockpit; Spacesuit Data Display and Management System; IEEE 1394 Hub With Fault Containment; Compact, Miniature MMIC Receiver Modules for an MMIC Array Spectrograph; Waveguide Transition for Submillimeter-Wave MMICs; Magnetic-Field-Tunable Superconducting Rectifier; Bonded Invar Clip Removal Using Foil Heaters; Fabricating Radial Groove Gratings Using Projection Photolithography; Gratings Fabricated on Flat Surfaces and Reproduced on Non-Flat Substrates; Method for Measuring the Volume-Scattering Function of Water; Method of Heating a Foam-Based Catalyst Bed; Small Deflection Energy Analyzer for Energy and Angular Distributions; Polymeric Bladder for Storing Liquid Oxygen; Pyrotechnic Simulator/Stray-Voltage Detector; Inventions Utilizing Microfluidics and Colloidal Particles; RuO2 Thermometer for Ultra-Low Temperatures; Ultra-Compact, High-Resolution LADAR System for 3D Imaging; Dual-Channel Multi-Purpose Telescope; Objective Lens Optimized for Wavefront Delivery, Pupil Imaging, and Pupil Ghosting; CMOS Camera Array With Onboard Memory; Quickly Approximating the Distance Between Two Objects; Processing Images of Craters for Spacecraft Navigation; Adaptive Morphological Feature-Based Object Classifier for a Color Imaging System; Rover Slip Validation and Prediction Algorithm; Safety and Quality Training Simulator; Supply-Chain Optimization Template; Algorithm for Computing Particle/Surface Interactions; Cryogenic Pupil Alignment Test Architecture for Aberrated Pupil Images; and Thermal Transport Model for Heat Sink Design.

Source record↗

The Principle of Energetic Consistency

A basic result in estimation theory is that the minimum variance estimate of the dynamical state, given the observations, is the conditional mean estimate. This result holds independently of the specifics of any dynamical or observation nonlinearity or stochasticity, requiring only that the probability density function of the state, conditioned on the observations, has two moments. For nonlinear dynamics that conserve a total energy, this general result implies the principle of energetic consistency: if the dynamical variables are taken to be the natural energy variables, then the sum of the total energy of the conditional mean and the trace of the conditional covariance matrix (the total variance) is constant between observations. Ensemble Kalman filtering methods are designed to approximate the evolution of the conditional mean and covariance matrix. For them the principle of energetic consistency holds independently of ensemble size, even with covariance localization. However, full Kalman filter experiments with advection dynamics have shown that a small amount of numerical dissipation can cause a large, state-dependent loss of total variance, to the detriment of filter performance. The principle of energetic consistency offers a simple way to test whether this spurious loss of variance limits ensemble filter performance in full-blown applications. The classical second-moment closure (third-moment discard) equations also satisfy the principle of energetic consistency, independently of the rank of the conditional covariance matrix. Low-rank approximation of these equations offers an energetically consistent, computationally viable alternative to ensemble filtering. Current formulations of long-window, weak-constraint, four-dimensional variational methods are designed to approximate the conditional mode rather than the conditional mean. Thus they neglect the nonlinear bias term in the second-moment closure equation for the conditional mean. The principle of energetic consistency implies that, to precisely the extent that growing modes are important in data assimilation, this term is also important.

Cohn, Stephen E.↗

Broken Ergodicity in MHD Turbulence

Ideal magnetohydrodynamic (MHD) turbulence may be represented by finite Fourier series, where the inherent periodic box serves as a surrogate for a bounded astrophysical plasma. Independent Fourier coefficients form a canonical ensemble described by a Gaussian probability density function containing a Hermitian covariance matrix with positive eigenvalues. The eigenvalues at lowest wave number can be very small, resulting in a large-scale coherent structure: a turbulent dynamo. This is seen in computations and a theoretical explanation in terms of 'broken ergodicity' contains Taylor s theory of force-free states. An important problem for future work is the case of real, i.e., dissipative flows. In real flows, broken ergodicity and coherent structure are still expected to occur in MHD turbulence at the largest scale, as suggested by low resolution simulations. One challenge is to incorporate coherent structure at the largest scale into the theory of turbulent fluctuations at smaller scales.

Shebalin, John V.↗

Displaying Properties of PDFs

PDFVis is a computer program that assists in visualization of uncertainty as represented by a probability density function (PDF) located at each grid cell in a spatial domain. The functions that PDFV performs are listed.

Kao, David↗

Analytic Development of a Reference Profile for the First Entry in a Skip Atmospheric Entry

This note shows that a feasible reference drag profile for the first entry portion of a skip entry can be generated as a polynomial expression of the velocity. The coefficients of that polynomial are found through the resolution of a system composed of m + 1 equations, where m is the degree of the drag polynomial. It has been shown that a minimum of five equations (m = 4) are required to establish the range and the initial and final conditions on velocity and flight path angle. It has been shown that at least one constraint on the trajectory can be imposed through the addition of one extra equation in the system, which must be accompanied by the increase in the degree of the drag polynomial. In order to simplify the resolution of the system of equations, the drag was considered as being a probability density function of the velocity, with the velocity as a distribution function of the drag. Combining this notion with the introduction of empirically derived constants, it has been shown that the system of equations required to generate the drag profile can be successfully reduced to a system of linear algebraic equations. For completeness, the resulting drag profiles have been flown using the feedback linearization method of differential geometric control as a guidance law with the error dynamics of a second order homogeneous equation in the form of a damped oscillator. Satisfactory results were achieved when the gains in the error dynamics were changed at a certain point along the trajectory that is dependent on the velocity and the curvature of the drag as a function of the velocity. Future work should study the capacity to update the drag profile in flight when dispersions are introduced. Also, future studies should attempt to link the first entry, as presented and controlled in this note, with a more standard control concept for the second entry, such as the Apollo entry guidance, to try to assess the overall skip entry performance. A guidance law that includes an integral feedback term, as is the case in the actual Space Shuttle entry guidance and as is proposed in Ref 29, could be tried in future studies to assess whether its use results in an improvement of the tracking performance, and to evaluate the design needs when determining the control gains.

Garcia-Llama, Eduardo↗

Conservational PDF Equations of Turbulence

Recently we have revisited the traditional probability density function (PDF) equations for the velocity and species in turbulent incompressible flows. They are all unclosed due to the appearance of various conditional means which are modeled empirically. However, we have observed that it is possible to establish a closed velocity PDF equation and a closed joint velocity and species PDF equation through conditions derived from the integral form of the Navier-Stokes equations. Although, in theory, the resulted PDF equations are neither general nor unique, they nevertheless lead to the exact transport equations for the first moment as well as all higher order moments. We refer these PDF equations as the conservational PDF equations. This observation is worth further exploration for its validity and CFD application

Shih, Tsan-Hsing↗

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

Characterization of Cloud Water-Content Distribution

The development of realistic cloud parameterizations for climate models requires accurate characterizations of subgrid distributions of thermodynamic variables. To this end, a software tool was developed to characterize cloud water-content distributions in climate-model sub-grid scales. This software characterizes distributions of cloud water content with respect to cloud phase, cloud type, precipitation occurrence, and geo-location using CloudSat radar measurements. It uses a statistical method called maximum likelihood estimation to estimate the probability density function of the cloud water content.

Lee, Seungwon↗

The Vertical Distribution of Thin Features Over the Arctic Analysed from CALIPSO Observations: Aerosols - Part 2

Influx of aerosols from the mid-latitudes has a wide range of impacts on the Arctic atmosphere. In this study, the capability of the CALIPSO-CALIOP instrument to provide accurate observations of aerosol layers is exploited to characterize their vertical distribution, probability density functions (PDFs) of aerosol layer thickness, base and top heights, and optical depths over the Arctic for the 4-yr period from June 2006 to May 2010. It is shown that the bulk of aerosols, from about 65% in winter to 45% in summer, are confined below the lowermost kilometer of the troposphere. In the middle troposphere (3-5 km), spring and autumn seasons show slightly higher aerosol amounts compared to other two seasons. The relative vertical distribution of aerosols shows that clean continental aerosol is the largest contributor in all seasons except in summer, when layers of polluted continental aerosols are almost as large. In winter and spring, polluted continental aerosols are the second largest contributor to the total number of observed aerosol layers, whereas clean marine aerosol is the second largest contributor in summer and autumn. The PDFs of the geometrical thickness of the observed aerosol layers peak about 400-700 m. Polluted continental and smoke aerosols, which are associated with the intrusions from mid-latitudes, have much broader distributions of optical and geometrical thicknesses, suggesting that they appear more often optically thicker and higher up in the troposphere.

Devasthale, Abhya↗

Broken Ergodicity in MHD Turbulence in a Spherical Domain

Broken ergodicity (BE) occurs in Fourier method numerical simulations of ideal, homogeneous, incompressible magnetohydrodynamic (MHD) turbulence. Although naive statistical theory predicts that Fourier coefficients of fluid velocity and magnetic field are zero-mean random variables, numerical simulations clearly show that low-wave-number coefficients have non-zero mean values that can be very large compared to the associated standard deviation. In other words, large-scale coherent structure (i.e., broken ergodicity) in homogeneous MHD turbulence can spontaneously grow out of random initial conditions. Eigenanalysis of the modal covariance matrices in the probability density functions of ideal statistical theory leads to a theoretical explanation of observed BE in homogeneous MHD turbulence. Since dissipation is minimal at the largest scales, BE is also relevant for resistive magnetofluids, as evidenced in numerical simulations. Here, we move beyond model magnetofluids confined by periodic boxes to examine BE in rotating magnetofluids in spherical domains using spherical harmonic expansions along with suitable boundary conditions. We present theoretical results for 3-D and 2-D spherical models and also present computational results from dynamical simulations of 2-D MHD turbulence on a rotating spherical surface. MHD turbulence on a 2-D sphere is affected by Coriolus forces, while MHD turbulence on a 2-D plane is not, so that 2-D spherical models are a useful (and simpler) intermediate stage on the path to understanding the much more complex 3-D spherical case.

Shebalin, John V.↗