Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “exponential probability density”

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.

30 records · Page 2

Eulerian Mapping Closure Approach for Probability Density Function of Concentration in Shear Flows

The Eulerian mapping closure approach is developed for uncertainty propagation in computational fluid mechanics. The approach is used to study the Probability Density Function (PDF) for the concentration of species advected by a random shear flow. An analytical argument shows that fluctuation of the concentration field at one point in space is non-Gaussian and exhibits stretched exponential form. An Eulerian mapping approach provides an appropriate approximation to both convection and diffusion terms and leads to a closed mapping equation. The results obtained describe the evolution of the initial Gaussian field, which is in agreement with direct numerical simulations.

He, Guowei↗

Highly ionized interstellar gas located in the Galactic disk and halo

High-resolution IUE absorption line spectra have been obtained for 40 distant stars in order to study the distribution of interstellar H I, Si IV, C IV, and N V in the Galactic disk and lower halo. Respective midplane densities of 2 x 10 to the -9th, 7 x 10 to the -9th, and 3 x 10 to the -9th are found for Si IV, C IV, and Ni V. Both column density and velocity data indicate that the highly ionized gas (HIG) is considerably more extended in directions away from the Galactic plane than is H I or Si II. The absorption-line velocities for the halo HIG are consistent with the notion that halo gas in the inner Galaxy rotates more slowly than gas in the underlying disk. The derived column densities suggest an exponential scale height for the HIG of about 3 kpc; however, a simple exponential distribution is a poor representation of the distribution of the gas. It is concluded that a full explanation of the origin of the halo HIG will probably require a blending of ideas from the Galactic fountain and the photoionized halo models.

Savage, Blair D.↗

Dispersion and Lifetime of the SO2 Cloud from the August 2008 Kasatochi Eruption

Hemispherical dispersion of the SO2 cloud from the August 2008 Kasatochi eruption is analyzed using satellite data from the Ozone Monitoring Instrument (OMI) and the Goddard Trajectory Model (GTM). The operational OMI retrievals underestimate the total SO2 mass by 20-30% on 8-11 August, as compared with more accurate offline Extended Iterative Spectral Fit (EISF) retrievals, but the error decreases with time due to plume dispersion and a drop in peak SO2 column densities. The GTM runs were initialized with and compared to the operational OMI SO2 data during early plume dispersion to constrain SO2 plume heights and eruption times. The most probable SO2 heights during initial dispersion are estimated to be 10-12 km, in agreement with direct height retrievals using EISF algorithm and IR measurements. Using these height constraints a forward GTM run was initialized on 11 August to compare with the month-long Kasatochi SO2 cloud dispersion patterns. Predicted volcanic cloud locations generally agree with OMI observations, although some discrepancies were observed. Operational OMI SO2 burdens were refined using GTM-predicted mass-weighted probability density height distributions. The total refined SO2 mass was integrated over the Northern Hemisphere to place empirical constraints on the SO2 chemical decay rate. The resulting lower limit of the Kasatochi SO2 e-folding time is approx.8-9 days. Extrapolation of the exponential decay back in time yields an initial erupted SO2 mass of approx.2.2 Tg on 8 August, twice as much as the measured mass on that day.

Krotkov, N. A.↗

Results of critical velocity experiments with barium, strontium, and calcium releases from CRRES satellite

As part of the NASA Combined Release and Radiation Effects Satellite (CRRES) chemical release program in September 1990, two Ba and also one each Sr and Ca canisters of a boron-titanium thermite mixture, which vaporizes the element on ignition, were released near perigee after dusk in the South Pacific to study the critical velocity effect proposed by Alfven. The critical velocities of these three elements are 2.7, 3.5, and 5.4 km/s respectively, all well below the orbital velocity of 9.4 km/s. On September 10, 1990, a Sr and Ba pair (G-13, or critical ionization velocity (CIV) I) was released near Rarotonga at approximately 515 km altitude in a background electron density of 3.4 x 10(exp 6)/cu cm. On September 14, 1990, G-14 or CIV II released a Ca and Ba pair west of New Caledonia near 595 km at an electron density of 1.5 x 10(exp 6)/cu cm. Ions of all three elements were observed with low-light level imagers from two aircraft after they had transited up the magnetic field lines into the sunlight. Emissions from the spherically expanding neutral gas shells below the solar terminator, observed with cameras filtered for the Ba(+) ion line at 4554 A and also in unfiltered imagers for approximately 15 s after release, are probably due to excitation by hot electrons created in the CIV process. The ions created clearly lost much of their energy, which we now show can be explained by elastic collisions: Ba(+) + O. Inventories of the observed ions indicate yields of 0.15% and 1.84% for Ba in the first and second experiments, 0.02% for Sr and 0.27% for Ca. Ionization from all the releases continued along the satellite trajectory much longer (greater than 45 s) than expected for a CIV process. The ion production along the satellite track versus time typically shows a rapid rise to a peak in a few seconds followed by an exponential decrease to a level essentially constant rate. The characteristic distances for CIV I and II are 47 and 62 km, respectively. We interpret the early time rise and exponential fall to be due to CIV ionization, of 0.014% (CIV I) and 0.40% (CIV II) for the Ba releases. The later ions produced at a constant rate probably have origins from other such processes as stripping and associative ionization collisions with atmospheric constituents primarily O, and charge exchange with O(+), He(+), and H(+). We suggest that the much larger Ba ionization rate in CIV II than CIV I is due to the fact that the release occurred in the peak Ca density where hot electrons were already present.

Wescott, E. M.↗

Coronal gas in the Galactic halo: ORFEUS observations of NGC 346 #1

We present high resolution (lambda/delta lambda = 3000) observations of the Small Magellanic Cloud (SMC) star NGC 346 #1 near the O VI resonance doublet (1031.9, 1037.6 A). The data were collected with the Berkeley extreme/far ultraviolet (EUV/FUV) spectrometer in the ORFEUS telescope aboard the space platform Astro-SPAS during the missiom of space shuttle Discovery flown in September 1993. The O VI features shows broad absorption centered near nu(sub LSR) is approximately equal to +140 km/sec. This is almost certainly associated with the OB cluster NGC 346 itself, which also shows deep absorption at Si IV and C IV at comparable velocities. An absorption feature corresponding to O VI near nu(sub LSR) = -230 km/sec is present in our data, but this feature is caused by contamination from a second star, probably an early B-type main sequence star at lower standard of rest (LSR) velocity, which is also present in the 20 min circular aperture. Between these velocity extrema lies one broad feature (b is approximately equal 70 km/sec) or multiple narrower features, with a combined column density of about 1.4 x 10(exp 14)/sq cm. Assuming an exponential distribution of O VI with the in-plane density of Jenkins, we infer a scale height, H(sub 0), of no more than approximately equal 1.2 kpc toward NGC 346 #1. Our O VI column density is about a factor of 4 below the predictions of cooling fountain models, but our low scale height is consistent with more recent predictions.

Hurwitz, M.↗

Hidden Statistics Approach to Quantum Simulations

Recent advances in quantum information theory have inspired an explosion of interest in new quantum algorithms for solving hard computational (quantum and non-quantum) problems. The basic principle of quantum computation is that the quantum properties can be used to represent structure data, and that quantum mechanisms can be devised and built to perform operations with this data. Three basic non-classical properties of quantum mechanics superposition, entanglement, and direct-product decomposability were main reasons for optimism about capabilities of quantum computers that promised simultaneous processing of large massifs of highly correlated data. Unfortunately, these advantages of quantum mechanics came with a high price. One major problem is keeping the components of the computer in a coherent state, as the slightest interaction with the external world would cause the system to decohere. That is why the hardware implementation of a quantum computer is still unsolved. The basic idea of this work is to create a new kind of dynamical system that would preserve the main three properties of quantum physics superposition, entanglement, and direct-product decomposability while allowing one to measure its state variables using classical methods. In other words, such a system would reinforce the advantages and minimize limitations of both quantum and classical aspects. Based upon a concept of hidden statistics, a new kind of dynamical system for simulation of Schroedinger equation is proposed. The system represents a modified Madelung version of Schroedinger equation. It preserves superposition, entanglement, and direct-product decomposability while allowing one to measure its state variables using classical methods. Such an optimal combination of characteristics is a perfect match for simulating quantum systems. The model includes a transitional component of quantum potential (that has been overlooked in previous treatment of the Madelung equation). The role of the transitional potential is to provide a jump from a deterministic state to a random state with prescribed probability density. This jump is triggered by blowup instability due to violation of Lipschitz condition generated by the quantum potential. As a result, the dynamics attains quantum properties on a classical scale. The model can be implemented physically as an analog VLSI-based (very-large-scale integration-based) computer, or numerically on a digital computer. This work opens a way of developing fundamentally new algorithms for quantum simulations of exponentially complex problems that expand NASA capabilities in conducting space activities. It has been illustrated that the complexity of simulations of particle interaction can be reduced from an exponential one to a polynomial one.

Zak, Michail↗

Analysis of the photospheric epoch spectra of type 1a supernovae SN 1990N and SN 1991T

This paper presents an LTE analysis of the photospheric epoch optical and IUE UV spectra of the recent Type Ia supernovae SN 1990N and SN 1991T. It is found that SN 1990N has matter moving as fast as 40,000 km/s, while SN 1991T has matter moving as fast as 20,000 km/s. The exponential density profile with e-folding velocity of 3160 km/s suggested by deflation and delayed/late-detonation explosion models is adequate for calculating fits to the observed spectra. Findings on the iron peak elements suggest that in both SNe some nuclear burning continues into the outer ejecta or the newly synthesized elements are mixed into the outer ejecta. The observed spectra of SN 1990N are consistent with a composition of the inner envelope like that of model W7. Intermediate-mass elements with abundances higher than solar are probably in the outer ejecta of both SNe, while silicon, sulfur, and calcium are underabundance in SN 1991T relative to SN 1990N by factors of order 3, 3, and 120, respectively.

Jeffery, David J.↗

An IUE survey of interstellar H I Ly alpha absorption. 2: Interpretations

We present an analysis of interstellar neutral hydrogen column densities obtained from archival Ly alpha absorption line data toward the complete sample of B2 and hotter stars observed at high spectral resolution with the IUE satellite. The full sample includes 554 stars and more than doubles the number of lines of sight previously observed. Many of the B2 and B1.5 stars exhibit contamination from stellar Ly alpha absorption, and were excluded from the subsequent analysis. The final working sample includes 393 stars. We present statistical averages for a number of quantities derived for the sample. The stars range in distance from 0.12 to 11 kpc with an average distance of 2.1 kpc. Values for average sight-line density range from 0.017 to 8.62 atoms cm(exp -3), with an overall average of 0.23 atoms cm(exp -3). Neutral hydrogen and dust are well correlated, with N(H I)/E(B - V) = 4.93 x 10(exp 21) cm(exp -2) mag(exp -1), and N(H I)/E(Bump) = 2.17 x 10(exp 21) cm(exp -2) mag(exp -1). There is an increase in the H I to dust ratios for the densest sight lines. Evidently, the H I to dust ratios in dense clouds are higher than in the lower density medium between the clouds and in interarm directions. The effect is probably an indication of dust modification in the dense clouds. The H I to dust ratios are approx. 17% smaller for the sight lines to O stars compared to the ratios found toward B stars. We have identified sight lines with anomalously large and small values of N(H I)/E(B - V) and H(H I)/E(Bump). Some of these sight lines are well-known examples of peculiar extinction (i.e., HD 147933 and HD 37061 (NU Ori)). Others are less well known and may provide important new examples of regions with highly modified dust. The number of stars with anomalously low H I to dust ratios is disproportionately larger than the number of stars with anomalously high H I to dust ratios. The distribution of the gas away form the Galactic plane for the sample of objects is complex. For the subset of 375 stars whose lines of sight have a statistically low probablility of intersecting a large cloud, the distribution may be roughly approximated by an exponential density distribution, with a midplane density of 0.366 atoms cm(exp -3), a scale height of 195 pc, and a random logarithmic scattering parameter sigma (sub p) of 0.159 dex. For the same sample of stars, interstellar dust is found to be confined closer to the plane than the H I, with (E(B - V)/r)(sub 0) = 0.257 mag kpc(exp -1), scale height h = 152 pc, and sigma(sub p) = 0.263 dex. The dust distribution is more inhomogeneous than the gas distribution. The observed distribution of H I away from the Galactic plane can also be fitted with a two-component exponential model including the effects of patchiness. Unfortunately, the foreground contamination produced by the cloudy and more confined component of H I makes it difficult to estimate the scale height of the extended component without invoking an independent estimate of its midplane density. Assuming n(sub 2)(0) = 0.16 cm(exp -3) from the Copernicus satellite H I/H(sub 2) survey, we obtain n(sub 1)(0) = 0.247 cm(exp -3), h(sub 1) = 73 PC, and h(sub 2) + 357 pc.

Diplas, Athanassios↗

Day 640 infrared line and continuum measurements: Dust formation in SN 1987A

We have measured day 640-645 line and continuum spectra of (Ni II) 6.6 micrometer (Ne II) 12.8 micrometer (line emission was not detected), and (Fe II) 17.9 and 26.0 micrometer from SN 1987A. The high velocity feature at v(sub HVF) approximately 3900 km/sec found in both of our day 410 (Fe II) spectra is again detected in the day 640 (Ni II) spectrum, although the signal-to-noise of the day 640 (Fe II) spectra is insufficient to show this feature. The continuum fluxes provide clear evidence for the formation of dust between day 410 and day 640 and are best fitted by a graybody spectrum with a temperature of 342 +/- 17 K at day 640 and a surface area corresponding to a minimum dust velocity v(sub dust) = 1910 +/- 170 km/sec. Optically thin dust emissivity laws proportional to lambda(exp -1) or lambda(exp -2) are inconsistent with the data. Either the dust grains are large (radius a much greater than 4 micrometer and radiate like individual blackbodies, or else they are located in clumps optically thick in the 6-26 micrometer range. The (Ni II) 6.6 micrometer line flux yields a minimum Ni(+) mass of 5.8 +/- 1.6 x 10(exp -4) solar mass and a Ni/Fe abundance ratio of 0.06 +/- 0.02, equal to the solar value. The ratio of the two (Fe II) line profiles implies a gas temperature 2600 +/- 700 K, a drop of 1800 +/- 800 K from our day 410 measurement. The (Fe II) 26.0 micrometer line flux has decreased by a factor of 2 and the day 640 (Ni II) profile is blueshifted by -440 +/- 270 km/sec, relative to observations before day 500. We show that the decrease in the (Fe II) flux and the blueshift are not produced by a decrease in electron scattering optical depth, electron density, or temperature, but rather are probably due to obscuration by the same dust which produces the infrared continuum. This supports the interpretation that the dust spectrum is produced by optically thick clumps. We discuss possible explanations for the discrepancy between the mass of Fe(+) detected and the total iron mass required to power the light curve. The decrease in the (Fe II) fluxes relative to the decrease required to account for the blueshifts of optical lines from non-iron-group elements and the similarity between v(sub dust) and the Ni(+) expansion velocity imply a spatial association between the dust clumps and the iron-group elements. In addition, the larger blueshift observed for the near and far-infrared, heavy metal transitions relative to non-iron-group lines suggests that the iron-group elements are somewhat segregated from lighter elements such as the Mg(sup 0) and O(sup 0) responsible for shorter wavelength lines. We speculate that FeS may be an important constituent of the dust. A comparison of our line profiles with radiative transfer models shows that while power law and exponential density distributions yield reasonable fits to the data, polytrope distributions provided significantly worse agreement. The best fits require a substantial fraction of the iron to be undetectable, and are consistent with maximum expansion velocities of v(sub max) approximately 3000 km/sec.

Colgan, Sean W. J.↗

The Impact of Pixel Size on the Characterization of Deep Convective Clouds for Calibration

The NASA CERES project provides the scientific community the observed TOA SW and LW fluxes for climate monitoring and climate model validation. CERES utilizes hourly geostationary imager derived broadband fluxes, which rely on the channel radiances and associated cloud retrievals, are used to estimate the broadband fluxes between CERES observations. This requires stable and consistent cross-platform imager visible channel calibration. The CERES project utilizes deep convective clouds (DCC) as an invariant Earth target to both monitor the stability of sensors and for radiometric scaling. GSICS, an international collaboration, is also evaluating and implementing the DCC invariant target calibration methodology to provide consistent calibration coefficients across geostationary imagers anchored to the AquaMODIS or the NOAA-20 VIIRS calibration reference. Tropical DCC are the brightest, coldest, most Lambertian, top of the atmosphere Earth targets. The DCC invariant target calibration methodology relies on a large ensemble of tropical D CC-identified pixel-level reflectances, which are aggregated as probability density functions (PDF). By assuming the monthly PDF shape is otherwise consistent in time excepting shifts in reflectance caused by changes in the sensor calibration, the imager stability is monitored. Radiometric scaling is accomplished by ratioing the sensor pair DCC PDF reflectance values. The success of the DCC methodology relies on consistent PDF distributions. The goal of this study is to determine the impact of pixel resolution on the DCC reflectance distribution. Single SNPP-VIIRS 750-m and Landsat 8 OLI 30-m granules are aggregated to degrade the pixel resolution from the native level. The DCC pixels are identified using a BT threshold. Most of the brightest DCC pixels are also the coldest, although there are exceptions. It was found that increasing the BT threshold exponentially increased the number of darker pixels. The pixel resolution did not seem to impact the DCC reflectance PDF distribution for pixel resolutions less than 3 km, which suggests that imagers of varying pixel resolutions may be radiometrically scaled to each other using DCC targets.

DCC↗

Dynamics of flux tubes in accretion disks

The study of magnetized plasmas in astrophysics is complicated by a number of factors, not the least of which is that in considering magnetic fields in stars or accretion disks, we are considering plasmas with densities well above those we can study in the laboratory. In particular, whereas laboratory plasmas are dominated by the confining magnetic field pressure, stars, and probably accretion disks, have magnetic fields whose beta (ratio of gas pressure to magnetic field pressure) is much greater than 1. Observations of the Sun suggest that under such circumstances the magnetic field breaks apart into discrete flux tubes with a small filling factor. On the other hand, theoretical treatments of MHD turbulence in high-beta plasmas tend to assume that the field is more or less homogeneously distributed throughout the plasma. Here we consider a simple model for the distribution of magnetic flux tubes in a turbulent medium. We discuss the mechanism by which small inhomogeneities evolve into discrete flux tubes and the size and distribution of such flux tubes. We then apply the model to accretion disks. We find that the fibrilation of the magnetic field does not enhance magnetic buoyancy. We also note that the evolution of an initially diffuse field in a turbulent medium, e.g., any uniform field in a shearing flow, will initially show exponential growth as the flux tubes form. This growth saturates when the flux tube formation is complete and cannot be used as the basis for a self-sustaining dynamo effect. Since the typical state of the magnetic field is a collection of intense flux tubes, this effect is of limited interest. However, it may be important early in the evolution of the galactic magnetic field, and it will play a large role in numerical simulations. Finally, we note that the formation of flux tubes is an essential ingredient in any successful dynamo model for stars or accretion disks.

Vishniac, E. T.↗

Heating of the solar middle chromosphere by large-scale electric currents

A global resistive, two-dimensional, time-dependent magnetohydrodynamic (MHD) model is used to introduce and support the hypothesis that the quiet solar middle chromosphere is heated by resistive dissipation of large-scale electric currents which fill most of its volume. The scale height and maximum magnitude of the current density are 400 km and 31.3 m/sq m, respectively. The associated magnetic field is almost horizontal, has the same scale height as the current density, and has a maximum magnitude of 153 G. The current is carried by electrons flowing across magnetic field lines at 1 m/s. The resistivity is the electron contribution to the Pedersen resitivity for a weakly ionized, strongly magnetized, hydrogen gas. The model does not include a driving mechanism. Most of the physical quantities in the model decrease exponentially with time on a resistive timescale of 41.3 minutes. However, the initial values and spatial; dependence of these quantities are expected to be essentially the same as they would be if the correct driving mechanism were included in a more general model. The heating rate per unit mass is found to be 4.5 x 10(exp 9) ergs/g/s, independent of height and latitude. The electron density scale height is found to be 800 km. The model predicts that 90% of the thermal energy required to heat the middle chromosphere is deposited in the height range 300-760 km above the temperature minimum. It is shown to be consistent to assume that the radiation rate per unit volume is proportional to the magnetic energy density, and then it follows that the heating rate per unit volume is also proportional to the energy from the photosphere into the overlying chromosphere are briefly discussed as possible driving mechanisms for establishing and maintaining the current system. The case in which part of or all of the current is carried by protons and metal ions, and the contribution of electron-proton scattering to the current are also considered, with the conclusion that these effects do not change the qualitative prediction of the model, but probably change the quantitative predictions slightly, mainly by increasing the maximum magntiude of the current density and magnetic field to at most approximately 100 mA/m and approximately 484 G, respectively. The heating rate per unit mass, current density scale height, magnetic field scale height, temperatures, and pressures are unchanged or are only slightly changed by including these additional effects due to protons and ions.

Goodman, M. L.↗