SEARCH · Engineering Papers
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.
A model of the largest gust amplitude and gust length is presented which uses the properties of the bivariate gamma distribution. The gust amplitude and length are strongly dependent on the filter function; the amplitude increases with altitude and is larger in winter than in summer.
The Ozone Mapping and Profiler Suite Limb Profiler (OMPS/LP) has been flying on the Suomi National Polar-orbiting Partnership (S-NPP) satellite since October 2011. It is designed to produce ozone and aerosol vertical profiles at 2-kilometer vertical resolution over the entire sunlit globe. Aerosol extinction profiles are computed with Mie theory using radiances measured at 675 nanometers. The operational Version 1.0 (V1.0) aerosol extinction retrieval algorithm assumes a bimodal lognormal aerosol size distribution (ASD) whose parameters were derived by combining an in situ measurement of aerosol microphysics with the Stratospheric Aerosol and Gas Experiment (SAGE II) aerosol extinction climatology. Internal analysis indicates that this bimodal lognormal ASD does not sufficiently explain the spectral dependence of LP-measured radiances. In this paper we describe the derivation of an improved aerosol size distribution, designated Version 1.5 (V1.5), for the LP retrieval algorithm. The new ASD uses a gamma function distribution that is derived from Community Aerosol and Radiation Model for Atmospheres (CARMA)-calculated results. A cumulative distribution fit derived from the gamma function ASD gives better agreement with CARMA results at small particle radii than bimodal or unimodal functions. The new ASD also explains the spectral dependence of LP-measured radiances better than the V1.0 ASD. We find that the impact of our choice of ASD on the retrieved extinctions varies strongly with the underlying reflectivity of the scene. Initial comparisons with collocated extinction profiles retrieved at 676 nanometers from the SAGE III instrument on the International Space Station (ISS) show a significant improvement in agreement for the LP V1.5 retrievals. Zonal mean extinction profiles agree to within 10 percent between 19 and 29 kilometers, and regression fits of collocated samples show improved correlation and reduced scatter compared to the V1.0 product. This improved agreement will motivate development of more sophisticated ASDs from CARMA results that incorporate latitude, altitude and seasonal variations in aerosol properties.
The CMIP5 configurations of the GISS Model-E2 GCM simulated a mid- and high latitude ice IWP that decreased by 50 relative to that simulated for CMIP3 (Jiang et al. 2012; JGR). Tropical IWP increased by 15 in CMIP5. While the tropical IWP was still within the published upper-bounds of IWP uncertainty derived using NASA A-Train satellite observations, it was found that the upper troposphere (200 mb) ice water content (IWC) exceeded the published upper-bound by a factor of 2. This was largely driven by IWC in deep-convecting regions of the tropics.Recent advances in the model-E2 convective parameterization have been found to have a substantial impact on tropical IWC. These advances include the development of both a cold pool parameterization (Del Genio et al. 2015) and new convective ice parameterization. In this presentation, we focus on the new parameterization of convective cloud ice that was developed using data from the NASA TC4 Mission. Ice particle terminal velocity formulations now include information from a number of NASA field campaigns. The new parameterization predicts both an ice water mass weighted-average particle diameter and a particle cross sectional area weighted-average size diameter as a function of temperature and ice water content. By assuming a gamma-distribution functional form for the particle size distribution, these two diameter estimates are all that are needed to explicitly predict the distribution of ice particles as a function of particle diameter.GCM simulations with the improved convective parameterization yield a 50 decrease in upper tropospheric IWC, bringing the tropical and global mean IWP climatologies into even closer agreement with the A-Train satellite observation best estimates.
Ths Special Sensor Microwave/Imager radiometer on board the DMSP satellite measured microwave radiation at 19.35, 22.235, 37.0, and 85.5 GHz with a swath width of 1400 km, providing an opportunity to study global precipitation distributions. A monthly averaged rainfall index was derived using only the 19.35 GHz data. It covers the + or - 50 deg portion of the world with a 5 deg x 5 deg grid. The brightness temperature histogram is analyzed to derive the rainfall distribution. Estimates of the freezing level heights and the water vapor contents are utilized to understand the rainfall distribution better. Log normal and gamma distribution functions are fitted to the data and the derived monthly rainfall index compares favorably with historical precipitation patterns.
This paper examines the relative impacts on grid-averaged longwave flux transmittance (emittance) for Marine Boundary Layer (MBL) cloud fields arising from horizontal variability of optical depth tau and cloud sides, First, using fields of Landsat-inferred tau and a Monte Carlo photon transport algorithm, it is demonstrated that mean all-sky transmittances for 3D variable MBL clouds can be computed accurately by the conventional method of linearly weighting clear and cloudy transmittances by their respective sky fractions. Then, the approximations of decoupling cloud and radiative properties and assuming independent columns are shown to be adequate for computation of mean flux transmittance. Since real clouds have nonzero geometric thicknesses, cloud fractions A'(sub c) presented to isotropic beams usually exceed the more familiar vertically projected cloud fractions A(sub c). It is shown, however, that when A(sub c)less than or equal to 0.9, biases for all-sky transmittance stemming from use of A(sub c) as opposed to A'(sub c) are roughly 2-5 times smaller than, and opposite in sign to, biases due to neglect of horizontal variability of tau. By neglecting variable tau, all-sky transmittances are underestimated often by more than 0.1 for A(sub c) near 0.75 and this translates into relative errors that can exceed 40% (corresponding errors for all-sky emittance are about 20% for most values of A(sub c). Thus, priority should be given to development of General Circulation Model (GCM) parameterizations that account for the effects of horizontal variations in unresolved tau, effects of cloud sides are of secondary importance. On this note, an efficient stochastic model for computing grid-averaged cloudy-sky flux transmittances is furnished that assumes that distributions of tau, for regions comparable in size to GCM grid cells, can be described adequately by gamma distribution functions. While the plane-parallel, homogeneous model underestimates cloud transmittance by about an order of magnitude when 3D variable cloud transmittances are less than or equal to 0.2 and by approx. 20% to 100% otherwise, the stochastic model reduces these biases often by more than 80%.
This research effort produced new methods to analyze the performance of linear predictors that track non-stationary processes. Specifically, prediction methods have been applied to the vibration pattern of rotorcraft drivetrains. This analysis is part or a larger rotorcraft Health and Usage Monitoring System (HUMS) that can diagnose immediate failures of the subsystems, as indicated by abrupt change in the vibration signature, and prognosticate future health, by examining the vibration patterns against long-term trends. This problem is described by a earlier joint paper co-authored by members of the funding agency and the recipient institutions prior to this grant effort. Specific accomplishments under this grant include the following: (1) Definition of a framework for analysis of non-stationary time-series estimation using the coefficients of an adaptive filter. (2) Description of a novel method of combining short-term predictor error and long-term regression error to analyze the performance of a non-stationary predictor. (3) Formulation of a multi-variate probability density function that quantifies the performance of a adaptive predictor by using the short- and long-term error variables in a Gamma function distribution. and (4) Validation of the mathematical formulations with empirical data from NASA flight tests and simulated data to illustrate the utility beyond the domain of vibrating machinery.
Overflow and underflow in sums prevented. Cumulative Poisson Distribution Program, CUMPOIS, one of two computer programs that make calculations involving cumulative Poisson distributions. Both programs, CUMPOIS (NPO-17714) and NEWTPOIS (NPO-17715), used independently of one another. CUMPOIS determines cumulative Poisson distribution, used to evaluate cumulative distribution function (cdf) for gamma distributions with integer shape parameters and cdf for X (sup2) distributions with even degrees of freedom. Used by statisticians and others concerned with probabilities of independent events occurring over specific units of time, area, or volume. Written in C.
We examine the width of the gamma ray burst luminosity function through the distribution of Gamma Ray Burst (GRB) peak fluxes as detected by the Pioneer Venus Orbiter (PVO) and the Burst and Transient Source Experiment (BATSE). The strength of the analysis is greatly enhanced by using a merged catalog of peak fluxes from both instruments with good cross-calibration of their sensitivities. The range of peak fluxes is increased by approximately a factor of 20 relative to the BATSE catalog. Thus, more sensitive investigations of the log N - log P distribution are possible. We place constraints on the width of the luminosity function of gamma-ray bursts brighter than the BATSE completeness limit by comparing the intensity distribution in the merged catalog with those produced by a variety of spatial density and luminosity functions. For the models examined, 90% of the detectable bursts have peak luminosities within a range of 10, indicating that the peak luminosities of gamma-ray bursts span a markedly less wide range of values than many other of their measurable properties. We also discuss for which slopes of a power-law luminosity function the observed width is at the upper end of the constrained range. This is essential in determining the power-law slopes for which luminosity-duration correlations could be important.
An analysis of two days of in situ observations of ice particle size spectra, in convectively generated cirrus, obtained during NASA s Tropical Composition, Cloud, and Climate Coupling (TC4) mission is presented. The observed spectra are examined for their fit to the exponential, gamma, and lognormal function distributions. Characteristic particle size and concentration density scales are determined using two (for the exponential) or three (for the gamma and lognormal functions) moments of the spectra. It is shown that transformed exponential, gamma, and lognormal distributions should collapse onto standard curves. An examination of the transformed spectra, and of deviations of the transformed spectra from the standard curves, shows that the lognormal function provides a better fit to the observed spectra.
Density and distribution functions of the sum of independent variables, each having a truncated gamma density function, were derived for use in the measurement of complex physical phenomena.
A model of multicomponent-liquid (MC-liquid) drop evaporation in a three-dimensional mixing layer is here exercised at larger Reynolds numbers than in a previous study, and transitional states are obtained. The gas phase is followed in an Eulerian frame and the multitude of drops is described in a Lagrangian frame. Complete coupling between phases is included with source terms in the gas conservation equations accounting for the drop/flow interaction in terms of drop drag, drop heating and species evaporation. The liquid composition, initially specified as a single-Gamma (SG) probability distribution function (PDF) depending on the molar mass is allowed to evolve into a linear combination of two SGPDFs, called the double-Gamma PDF (DGPDF). The compositions of liquid and vapor emanating from the drops are calculated through four moments of the DGPDFs, which are drop-specific and location-specific, respectively. The mixing layer is initially excited to promote the double pairing of its four initial spanwise vortices into an ultimate vortex in which small scales proliferate. Simulations are performed for four liquids of different compositions and the effect of the initial mass loading and initial free-stream gas temperature are explored. For reference, Simulations are also performed for gaseous multicomponent mixing layers for which the effect of Reynolds number is investigated. The results encompass examination of the global layer characteristics, flow visualizations and homogeneous-plane statistics at transition. Comparisons are performed with previous pre-transitional MC-liquid simulations and with transitional single-component (SC) liquid studies. It is found that MCC flows at transition, the classical energy cascade is of similar strength, but that the smallest scales contain orders of magnitude less energy than SC flows, which is confirmed by the larger viscous dissipation in the former case. Contrasting to pre-transitional MC flows, the vorticity and drop organization depend on the initial gas temperature, this being due to the drop/turbulence coupling. The vapor-composition mean molar mass and standard deviation distributions strongly correlate with the initial liquid-composition PDF; such a correlation only exists for the magnitude of the mean but not for that of the standard deviation. Unlike in pre-transitional situations, regions of large composition standard deviation no longer necessarily coincide with regions of large mean molar mass. The kinetic energy, rotational and composition characteristics, and dissipation are liquid specific and the variation among liquids is amplified with increasing free-stream gas temperature. Eulerian and Lagrangian statistics of gas-phase quantities show that the different. Observation framework may affect the perception of the flow characteristics. The gas composition, of which the first four moments are calculated, is shown to be close to, but distinct from a SGPDF. The PDF of the scalar dissipation rate is calculated for drop-laden layers and is shown to depart more significantly from the typically assumed Gaussian in gaseous flows than experimentally measured gaseous scalar dissipation rates, this being attributed to the increased heterogeneity due to drop/flow interactions.
A new technique for measuring cloud liquid water, mean droplet radius and droplet number density is outlined. The technique is based on simultaneously measuring Raman and Mie scattering from cloud liquid droplets using a Raman lidar. Laboratory experiments on liquid micro-spheres have shown that the intensity of Raman scattering is proportional to the amount of liquid present in the spheres. This fact is used as a constraint on calculated Mie intensity assuming a gamma function particle size distribution. The resulting retrieval technique is shown to give stable solutions with no false minima. It is tested using Raman lidar data where the liquid water signal was seen as an enhancement to the water vapor signal. The general relationship of retrieved average radius and number density is consistent with traditional cloud physics models. Sensitivity to the assumed maximum cloud liquid water amount and the water vapor mixing ratio calibration are tested. Improvements to the technique are suggested.
Observations by the airborne X-band Doppler radar (EDOP) and the NCAR S-band polarimetric (S-POL) radar from two field experiments are used to evaluate the Surface ref'ercnce technique (SRT) for measuring the path integrated attenuation (PIA) and to study attenuation in deep convective storms. The EDOP, flying at an altitude of 20 km, uses a nadir beam and a forward pointing beam. It is found that over land, the surface scattering cross-section is highly variable at nadir incidence but relatively stable at forward incidence. It is concluded that measurement by the forward beam provides a viable technique for measuring PIA using the SRT. Vertical profiles of peak attenuation coefficient are derived in vxo deep convective storms by the dual-wavelength method. Using the measured Doppler velocity, the reflectivities at. the two wavelengths, the differential reflectivity and the estimated attenuation coefficients, it is shown that: supercooled drops and dry ice particles probably co-existed above the melting level in regions of updraft, that water-coated partially melted ice particles probably contributed to high attenuation below the melting level, and that the data are not readil~ explained in terms of a gamma function raindrop size distribution.
In this study, we evaluate a decadal ozone profile record derived from Suomi NPP Ozone Mapping and Profiler Suite (OMPS) Limb Profiler (LP). The data from OMPS LP were re-processed with the new retrieval algorithm from 2012 to present. The new reprocessing includes adjustments in both Level 1 and Level 2 algorithms. The sensor’s altitude registration, updates in radiometric calibrations, under-sampling and stray light corrections were addressed in Level 1 processing. The new retrieval algorithm (version 2.7) combines measurements from UV and VIS parts of the spectra to retrieve a single vertical ozone profile between 12.5 km (or cloud tops) and 57.5 km, unlike to previous version 2.5 where ozone profiles were retrieved independently from UV and VIS. To improve ozone retrievals in the Upper Troposphere Lower Stratosphere (UTLS) region, we implemented a new aerosol correction that assumes a gamma-function particle size distribution. In addition, ozone and nitrogen dioxide cross sections and climatologies were updated. We discuss sensitivity of the ozone retrievals to implemented changes in Level 1 and Level 2 algorithms. We focus our analysis on evaluation of the stability of the 10-year record from OMPS LP. We show results from the internal methods for characterizing the stability of the LP altitude registration and compare the LP reflectivity and ozone with the similar quantities measured by OMPS nadir sensors. Finally, we will show results from comparisons with independent satellite (MLS and SAGE III) and ozone sonde measurements.
In this presentation we evaluate ozone profile retrievals from Ozone Mapping and Profiling Suite Limb Profiler (OMPS-LP) launched in October 2011 on Suomi NPP satellite. The data from OMPS-LP were re-processed with the new retrieval algorithm that combines measurements from UV and VIS parts of the spectra to retrieve a single ozone profile between 12.5 km (or cloud tops) and 57.5 km. To improve ozone retrievals in the Upper Troposphere Lower Stratosphere (UTLS) region, a new aerosol correction with a gamma-function particle size distribution was implemented. The OMPS-LP altitude registration was also adjusted in the new version to ensure the stability of the ozone record. The continuous ozone profile observations with Suomi NPP OMPS LP exceeds 10 years. We will compare OMPS-LP ozone profiles with coincident stratospheric ozone measurements from NASA Aura MLS, ISS SAGE III, and OMPS NP and with the observations in UTLS and middle stratosphere from a ground-based network of ozone sonde stations to evaluate accuracy and precision of OMPS-LP retrievals and overall stability of the record. We will also investigate the quality of stratospheric ozone columns obtained from OMPS-LP and estimate the contribution of various errors to integrated stratospheric columns from OMPS-LP.
We evaluate a decadal ozone profile record derived from the Suomi National Polar-orbiting Partnership (SNPP) Ozone Mapping and Profiler Suite (OMPS) Limb Profiler (LP) satellite instrument. In 2023, the OMPS LP data were re-processed with the new version 2.6 retrieval algorithm that combines measurements from the ultraviolet (UV) and visible (VIS) parts of the spectra and employs the second order Tikhonov regularization to retrieve a single vertical ozone profile between 12.5 km (or cloud tops) and 57.5 km with the vertical resolution of about 1.9 - 2.5 km between 20-55 km. The algorithm uses radiances measured at six UV ozone-sensitive wavelengths (295, 302, 306, 312, 317 and 322 nm) paired with 353 nm, and one VIS wavelength at 606 nm combined with 510 nm and 675 nm to form a triplet. Each wavelength pair or triplet is used over a limited range of tangent altitudes where the sensitivity to ozone changes are strongest. A new implemented aerosol correction scheme is based on a gamma-function particle size distribution. Numerous calibration changes that affected ozone retrievals were also applied to measured LP radiances, including updates in altitude registration, radiometric calibration, stray light, and spectral registration. The key version 2.6 improvement is the reduction in relative drifts between LP ozone and correlative measurements, linked previously to a drift in the version 2.5 LP altitude registration. We compare LP ozone profiles with those from Aura Microwave Limb Sensor (MLS) to quantify ozone changes in version 2.6.
The nearest BL Lac object, Mrk 421, has a gamma-ray spectrum which is approximately flat in EF-sub E from E less than about 50 MeV to E greater than about 1 TeV. Inverse Compton scattering can explain this smooth spectrum, despite the structure in the Klein-Nishina cross section, if the injected electron distribution function is proportional to gamma exp -2, where gamma is the electron Lorentz factor. When this is the case, the structure imprinted on the steady state electron distribution function by the structure in the Klein-Nishina cross section is almost exactly compensated in the radiated spectrum. Because particle acceleration in strong shocks injects particles with this distribution function, this shape injection function is in fact quite plausible. Other blazars may be explained by the same model if the cutoff below TeV energies observed in other objects is due to pair production on the cosmological IR background, as suggested by Stecker et al. (1992).
Both radiation and increased iron stores can independently increase oxidative damage, resulting in protein, lipid and DNA oxidation. Oxidative stress increases the risk of many health problems including cancer, cataracts, and heart disease. This study, a subset of a larger interdisciplinary investigation of the combined effect of iron overload on sensitivity to radiation injury, monitored immune parameters in the peripheral blood of rats subjected to gamma radiation, high dietary iron or both. Specific immune measures consisted of: (1) peripheral leukocyte distribution, (2) plasma cytokine levels and (3) cytokine production profiles following whole blood mitogenic stimulation