Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sigma model”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
NASA NTRS2016

Vibration Penalty Estimates for Indoor Annoyance Caused by Sonic Boom

Commercial supersonic flight is currently forbidden over land because sonic booms have historically caused unacceptable annoyance levels in overflown communities. NASA is providing data and expertise to noise regulators as they consider relaxing the ban for future quiet supersonic aircraft. One key objective is a predictive model for indoor annoyance based on factors such as noise and indoor vibration levels. The current study quantified the increment in indoor sonic boom annoyance when sonic booms can be felt directly through structural vibrations in addition to being heard. A shaker mounted below each chair in the sonic boom simulator emulated vibrations transmitting through the structure to that chair. The vibration amplitudes were determined from numeric models of a large range of residential structures excited by the same sonic boom waveforms used in the experiment. The analysis yielded vibration penalties, which are the increments in sound level needed to increase annoyance as much as the vibration does. For sonic booms at acoustic levels from 75 to 84 dB Perceived Level, vibration signals with lower amplitudes (+1 sigma) yielded penalties from 0 to 5 dB, and vibration signals with higher amplitudes (+3 sigma) yielded penalties from 6 to 10 dB.

Rathsam, Jonathan↗
NASA NTRS1984

Changes in the High-latitude Ocean as Possible Causes of Atmospheric CO2 Variations

Measurements on air enclosed in old polar ice have indicated that the atmospheric CO2 concentration was ca. 50 to 70 ppm lower in late glacial times than during the Holocene. Similar measurements performed on samples from a Greenland ice core, dating ca. 30,000 to 40,000 B.P., and have yielded evidence of several CO2 oscillations with an amplitude of ca. 50 ppm. Each change lasted on the order of a few centuries. A mechanism by which circulation changes in the high-latitude ocean could lead to rapid variations in atmospheric CO2 is proposed. In the Antarctic Ocean a slowing down of the vertical mixing would imply a smaller upward flux of sigma CO2 and nutrients. Assuming constant productivity, sigma CO2 and nutrients would be more completely used which would imply lower CO2 in these high-latitude surface waters. In areas with a warm surface, a slowing down of the circulation would not have a direct impact on CO2 because productivity would automatically decrease by the same factor as the upwelling rate of nutrients. Studies with a simple box model of the ocean-atmosphere system suggest that a suddent decrease by a factor of 2 of the water exchange between the surface and deep sea in high latitudes could lead to a CO2 decrease of ca. 40 to 50 ppm with a time constant of ca. 200 years. Deep-sea sediment studies indicate rapid changes in the high-latitude surface conditions of the North Atlantic and the Antarctic Oceans at the end of the last glaciation. Studies of carbon isotope ratios should help ascertain whether this proposed mechanism was indeed responsible for the CO2 variation.

Siegenthaler, U.↗
NASA NTRS1990

Navigation systems requirement analysis for the Aeroassist Flight Experiment (AFE)

Navigation requirements for Aeroassist Flight Experiment (AFE) spacecraft passing through the earth's atmosphere have been studied using a 6-DOF dynamics model, an Inertial Measurement Unit model, a baseline AFE aeropass flight guidance logic, and a baseline AFE aeropass control model. The goal of this study is to determine, in a statistical sense, how much flight path angle error can be tolerated at Entry Interface (EI) and still have acceptable delta-V requirements at exit to position the AFE spacecraft for recovery. Assuming there is fuel available to produce 370 ft/sec of delta-V at atmospheric exit, a 3-sigma standard deviation in flight path angle error of 0.04 degree at EI would result in a 98 percent probability of mission success. In addition to the required delta-V at exit, other aeropass parameters such as maximum aeroheating rate, fuel consumption, and the science requirements affecting mission success are also investigated.

Huang, K. H.↗
NASA NTRS2013

Enigma Version 12

Enigma Version 12 software combines model building, animation, and engineering visualization into one concise software package. Enigma employs a versatile user interface to allow average users access to even the most complex pieces of the application. Using Enigma eliminates the need to buy and learn several software packages to create an engineering visualization. Models can be created and/or modified within Enigma down to the polygon level. Textures and materials can be applied for additional realism. Within Enigma, these models can be combined to create systems of models that have a hierarchical relationship to one another, such as a robotic arm. Then these systems can be animated within the program or controlled by an external application programming interface (API). In addition, Enigma provides the ability to use plug-ins. Plugins allow the user to create custom code for a specific application and access the Enigma model and system data, but still use the Enigma drawing functionality. CAD files can be imported into Enigma and combined to create systems of computer graphics models that can be manipulated with constraints. An API is available so that an engineer can write a simulation and drive the computer graphics models with no knowledge of computer graphics. An animation editor allows an engineer to set up sequences of animations generated by simulations or by conceptual trajectories in order to record these to highquality media for presentation. Enigma Version 12 Lyndon B. Johnson Space Center, Houston, Texas 28 NASA Tech Briefs, September 2013 Planetary Protection Bioburden Analysis Program NASA's Jet Propulsion Laboratory, Pasadena, California This program is a Microsoft Access program that performed statistical analysis of the colony counts from assays performed on the Mars Science Laboratory (MSL) spacecraft to determine the bioburden density, 3-sigma biodensity, and the total bioburdens required for the MSL prelaunch reports. It also contains numerous tools that report the data in various ways to simplify the reports required. The program performs all the calculations directly in the MS Access program. Prior to this development, the data was exported to large Excel files that had to be cut and pasted to provide the desired results. The program contains a main menu and a number of submenus. Analyses can be performed by using either all the assays, or only the accountable assays that will be used in the final analysis. There are three options on the first menu: either calculate using (1) the old MER (Mars Exploration Rover) statistics, (2) the MSL statistics for all the assays, or This software implements penetration limit equations for common micrometeoroid and orbital debris (MMOD) shield configurations, windows, and thermal protection systems. Allowable MMOD risk is formulated in terms of the probability of penetration (PNP) of the spacecraft pressure hull. For calculating the risk, spacecraft geometry models, mission profiles, debris environment models, and penetration limit equations for installed shielding configurations are required. Risk assessment software such as NASA's BUMPERII is used to calculate mission PNP; however, they are unsuitable for use in shield design and preliminary analysis studies. The software defines a single equation for the design and performance evaluation of common MMOD shielding configurations, windows, and thermal protection systems, along with a description of their validity range and guidelines for their application. Recommendations are based on preliminary reviews of fundamental assumptions, and accuracy in predicting experimental impact test results. The software is programmed in Visual Basic for Applications for installation as a simple add-in for Microsoft Excel. The user is directed to a graphical user interface (GUI) that requires user inputs and provides solutions directly in Microsoft Excel workbooks. This work was done by Shannon Ryan of the USRA Lunar and Planetary Institute for Johnson Space Center. Further information is contained in a TSP (see page 1). MSC- 24582-1 Micrometeoroid and Orbital Debris (MMOD) Shield Ballistic Limit Analysis Program Lyndon B. Johnson Space Center, Houston, Texas Commercially, because it is so generic, Enigma can be used for almost any project that requires engineering visualization, model building, or animation. Models in Enigma can be exported to many other formats for use in other applications as well. Educationally, Enigma is being used to allow university students to visualize robotic algorithms in a simulation mode before using them with actual hardware.

Shores, David↗
NASA NTRS1982

A solar flare X-ray polarimeter for the Space Shuttle

An instrument has recently been built and tested which is designed to measure the polarization of the hard (5-30 keV) X-ray emission from solar flares, and thereby to investigate the energy release mechanism and constrain flare models. In particular, these measurements will help to determine whether hard X-ray bursts are produced by nonthermal or by thermal electrons. The polarimeter makes use of the angular dependence of Thomson scattering from targets of metallic lithium. It has an energy resolution of a few keV, a time resolution of 5 s, and sufficient sensitivity to measure polarization levels (3 sigma) of a few percent in about 10 s for a moderate strength solar flare. The instrumental polarization has been directly measured and found to be within the design goal of approximately 1%. This polarimeter is scheduled to be flown as part of the OSS-1 pallet on an early Space Shuttle mission.

Lemen, J. R.↗
NASA NTRS1994

Statistical association of QSO's with foreground galaxy clusters

We report a statistically significant overdensity of high redshift quasi-stellar objects (QSO's) in the directions of foreground galaxy clusters. QSO's are taken from the Large Bright QSO Survey (LBQS) between 1.4 less than or equal z less than or equal 2.2 with a limiting magnitude of m(sub B) = 18.5. Foreground clusters are regions within 6 Zwicky radii of small Zwicky clusters at a characteristic redshift of about z approximately = 0.2, covering about 40% of the total area surveyed (304 sq. deg). The overdensity, defined as the ratio of the number density of QSO's in the directions of clusters ('association QSO's) to that in the remainder of the fields ('background QSO's), is equal to 1.7, and formally differs from unity at 4.7 sigma significance. The observed overdensity probably is not due to statistical variation in QSO density, intrinsic QSO-QSO and/or cluster-cluster autocorrelations, or patchy Galactic obscuration. We thus interpret this observation as being due to statistical gravitational lensing of background QSO's by galaxy clusters. However, this amplitude of overdensity behind clusters cannot be accounted for in any cluster lensing model if the background QSO number-magnitude counts are similar to the intrinsic (unlensed) counts, and is implausible in any conventional model of cosmic mass distribution.

Rodrigues-Williams, Liliya L.↗
NASA NTRS1994

Measuring the Hubble constant and our Virgo-infall velocity independently

A sample of spiral galaxies with B(sub T) less than 14.5 located in two local volumes, one in the direction of, but behind, the Virgo Cluster (behind-Virgo volume (BV)) and the other in the opposite direction (anti-Virgo volume (AV)), were used via a Tully-Fisher (TF) relation to derive the following two parameters: H(sub AB), the mean Hubble ratio between AV and BV, and delta v(sub parallel), the peculiar velocity of the Local Group in the direction of the Virgo Cluster (VC) with respect to a uniformly expanding reference system defined by our AV and BV sub-samples. The two sampled volumes, separated by a velocity interval of 5600 km/s, form an antipodal pair. This particular geometry not only allows us to derive the two parameters independently but also reduces the dynamical effect of the Local Supercluster on H(sub AB) without increasing the Malmquist bias. By limiting our sample to spiral galaxies having large velocity widths W(sub R), we effectively reduce the TF scatter and Malmquist bias in our sample. The TF zero point and dispersion were then determined by further correcting for the small residual Malmquist bias. An additional sample of fainter galaxies was used to test for a non-Gaussian tail to the TF disperison. We found no evidence for such a tail and formally give an upper limit of about 18% for the fractional contribution of an unseen tail. The average intrinsic TF dispersion for the dominant Gaussian component is sigma(sub TF)(sup 0) approximately 0.33 mag for W(sub R) approximately equal to or greater than 180 km/s. Our numerical results are delta v(sub parallel) approximately equals 414 +/- 82 km/s and H(sub AB) approximately equals (84.0 +/- 2.4)(1 + epsilon) km/s Mpc, where (1 + epsilon) accounts for any systematic error between the calibrators and the sample galaxies. Various dynamical models were tested to explore the effect on H(sub AB) of the uncertainties in the local velocity field. Constrained by our observed delta v(sub parallel) as well as other observational quantities, we found that the rms deviation from unity of H(sub AB)/H(sub 0) (where H(sub 0) is the Hubble constant for each model) is 5%, making H(sub AB) a good indicator for H(sub 0). Taking this variation as an additional error, our formal estimate for the Hubble constant is H(sub 0) approximately equals (84 +/- 5)(1 + epsilon) km/s Mpc.

Lu, Nanyao Y.↗
NASA NTRS2005

Analysis of variability in the burst oscillations of the accreting millisecond pulsar XTE J1814-338

The accreting millisecond pulsar XTE J1814-338 exhibits oscillations at the known spin frequency during Type I X-ray bursts. The properties of the burst oscillations reflect the nature of the thermal asymmetry on the stellar surface. We present an analysis of the variability of the burst oscillations of this source, focusing on three characteristics: fractional amplitude, harmonic content and frequency. Fractional amplitude and harmonic content constrain the size, shape and position of the emitting region, whilst variations in frequency indicate motion of the emitting region on the neutron star surface. We examine both long-term variability over the course of the outburst, and short-term variability during the bursts. For most of the bursts, fractional amplitude is consistent with that of the accretion pulsations, implying a low degree of fuel spread. There is however a population of bursts whose fractional amplitudes are substantially lower, implying a higher degree of fuel spread, possibly forced by the explosive burning front of a precursor burst. For the first harmonic, substantial differences between the burst and accretion pulsations suggest that hotspot geometry is not the only mechanism giving rise to harmonic content in the latter. Fractional amplitude variability during the bursts is low; we can only rule out the hypothesis that the fractional amplitude remains constant at the l(sigma) level for bursts that do not exhibit photospheric radius expansion (PRE). There are no significant variations in frequency in any of the bursts except for the one burst that exhibits PRE. This burst exhibits a highly significant but small (= 0.1Hz) drop in frequency in the burst rise. The timescale of the frequency shift is slower than simple burning layer expansion models predict, suggesting that other mechanisms may be at work.

Watts, Anna L.↗
NASA NTRS2013

Application of Statistically Derived CPAS Parachute Parameters

The Capsule Parachute Assembly System (CPAS) Analysis Team is responsible for determining parachute inflation parameters and dispersions that are ultimately used in verifying system requirements. A model memo is internally released semi-annually documenting parachute inflation and other key parameters reconstructed from flight test data. Dispersion probability distributions published in previous versions of the model memo were uniform because insufficient data were available for determination of statistical based distributions. Uniform distributions do not accurately represent the expected distributions since extreme parameter values are just as likely to occur as the nominal value. CPAS has taken incremental steps to move away from uniform distributions. Model Memo version 9 (MMv9) made the first use of non-uniform dispersions, but only for the reefing cutter timing, for which a large number of sample was available. In order to maximize the utility of the available flight test data, clusters of parachutes were reconstructed individually starting with Model Memo version 10. This allowed for statistical assessment for steady-state drag area (CDS) and parachute inflation parameters such as the canopy fill distance (n), profile shape exponent (expopen), over-inflation factor (C(sub k)), and ramp-down time (t(sub k)) distributions. Built-in MATLAB distributions were applied to the histograms, and parameters such as scale (sigma) and location (mu) were output. Engineering judgment was used to determine the "best fit" distribution based on the test data. Results include normal, log normal, and uniform (where available data remains insufficient) fits of nominal and failure (loss of parachute and skipped stage) cases for all CPAS parachutes. This paper discusses the uniform methodology that was previously used, the process and result of the statistical assessment, how the dispersions were incorporated into Monte Carlo analyses, and the application of the distributions in trajectory benchmark testing assessments with parachute inflation parameters, drag area, and reefing cutter timing used by CPAS.

Romero, Leah M.↗
NASA NTRS1983

HEAO A-2 soft X-ray source catalog

The HEAO A-2 low energy X-ray detectors have surveyed over 95% of the sky in the spectral bands 0.18-0.44 and 0.44-2.8 keV, down to a typical limiting sensitivity of 1 x 10 to the -11th and 3 x 10 to the -11th ergs per sq cm per sec in each band, respectively. Using a significance criterion of 6 sigma for existence, 114 sources are cataloged, of which 54 were previously undiscovered and 32 remain unidentified. The catalog contains a listing of all counterpart identifications and a cross-reference to all HEAO 1 A-2 low energy detector team publications on the catalog sources complete through the end of 1981. Supplementary material allows the estimation of spectral parameters for the spectral models from the observed intensities. Simple graphs allow the extraction of energy fluxes at the earth, both with and without the effects of interstellar absorption. The angular distribution of the cataloged sources in the 0.44-2.8 keV band is concentrated toward the galactic plane and resembles that seen at higher energy. The distribution in the 0.18-0.44 keV band is more uniform and contains many more newly discovered and unidentified sources.

Nugent, J. J.↗
NASA NTRS1994

On determining the spectrum of primordial inhomogeneity from the COBE DMR sky maps: Results of two-year data analysis

A new technique of Fourier analysis on a cut sky has been applied to the two-year Cosmic Background Explorer (COBE) Differential Microwave Radiometer (DMR) 53 and 90 GHz sky maps. The Bayesian power spectrum estimation results are consistent with the Harrison-Zel'dovich n = 1 model. The maximum likelihood estimates of the usual parameters defining the power spectrum of primordial perturbations are n = 1.22 (1.02) and Q(sub rms-PS) = 17 (20) microK including (excluding) the quadrupole. A spectral-index-independent normalization is naturally expressed for the two-year maps in terms of the multipole amplitude a(sub 9) = 8.2 (8.3) microK (to approximately 12 sigma significance). The marginal likelihood function on n obtained by intergration with respect to a(sub 9) renders n = 1.17 +/- 0.31 (0.96 +/- 0.36).

Gorski, K. M.↗
NASA NTRS2014

Atmospheric Characterization of Five Hot Jupiters with the Wide Field Camera 3 on the Hubble Space Telescope

We probe the structure and composition of the atmospheres of five hot Jupiter exoplanets using the Hubble Space Telescope Wide Field Camera 3 (WFC3) instrument. We use the G141 grism (1.1-1.7 micrometers) to study TrES-2b, TrES-4b, and CoRoT-1b in transit; TrES-3b in secondary eclipse; and WASP-4b in both. This wavelength region includes a predicted absorption feature from water at 1.4 micrometers, which we expect to be nondegenerate with the other molecules that are likely to be abundant for hydrocarbon-poor (e.g., solar composition) hot Jupiter atmospheres. We divide our wavelength regions into 10 bins. For each bin we produce a spectrophotometric light curve spanning the time of transit or eclipse. We correct these light curves for instrumental systematics without reference to an instrument model. For our transmission spectra, our mean 1σ precision per bin corresponds to variations of 2.1, 2.8, and 3.0 atmospheric scale heights for TrES-2b, TrES-4b, and CoRoT-1b, respectively. We find featureless spectra for these three planets. We are unable to extract a robust transmission spectrum for WASP-4b. For our dayside emission spectra, our mean 1 sigma precision per bin corresponds to a planet-to-star flux ratio of 1.5 x 10(exp −4) and 2.1 x 10(exp −4) for WASP-4b and TrES-3b, respectively. We combine these estimates with previous broadband measurements and conclude that for both planets isothermal atmospheres are disfavored. We find no signs of features due to water. We confirm that WFC3 is suitable for studies of transiting exoplanets, but in staring mode multivisit campaigns are necessary to place strong constraints on water abundance.

eclipses↗
NASA NTRS2017

Direct Estimation of the Global Distribution of Vertical Velocity Within Cirrus Clouds

Cirrus clouds determine the radiative balance of the upper troposphere and the transport of water vapor across the tropopause.The representation of vertical wind velocity, W, in atmospheric models constitutes the largest source of uncertainty in the calculation of the cirrus formation rate. Using global atmospheric simulations with a spatial resolution of 7 kilometers we obtain for the first time a direct estimate of the distribution of W at the scale relevant for cirrus formation, validated against long-term observations at two different ground sites. The standard deviation in W, sigma (sub W), varies widely over the globe with the highest values resulting from orographic uplift and convection, and the lowest occurring in the Arctic. Globally about 90 of the simulated sigma (sub W) values are below 0.1 meters per second and about one in 10 (sup 4) cloud formation events occur in environments with sigma (sub W) greater than 0.8 meters per second. Combining our estimate with reanalysis products and an advanced cloud formation scheme results in lower homogeneous ice nucleation frequency than previously reported, and a decreasing average ice crystal concentration with decreasing temperature. These features are in agreement with observations and suggest that the correct parameterization of sigma (sub W) is critical to simulate realistic cirrus properties.

Radiative Balance↗
NASA NTRS1976

Recent Uhuru results on Centaurus X-3

The current status of the analysis of Cen X-3 data from Uhuru concerning pulsations, orbital period and eccentricity, and extended lows, are reviewed. The pulse period decreases irregularly. The pulsed fraction (2-7 keV) is 70%-90% for single pulses but significantly less for superpositions of pulses, due to variability in shape. The pulses are narrower at higher energies with a correlated increase in fraction pulsed. The orbital period is found to both decrease and increase with p/p on the order of a few times 0.00001/year. A three sigma upper limit on the eccentricity of 0.003 is obtained; if no significant periastron motion is allowed over two years, the upper limit becomes 0.0016. The orbital period is found to be detectable during some extended lows but with a significantly decreased ratio of eclipsed to non-eclipsed intensity. Two transitions between normal high states and extended lows are studied, and a consistent model is obtained in which extended lows are caused by both burying the source in an increased stellar wind from the companion, and starving the source by decreasing the stellar wind. Changes in fraction pulsed during transitions and systematic differences in the harmonic content of the pulses are also found.

Schreier, E. J.↗
NASA NTRS1983

The parameterization of the planetary boundary layer in the UCLA general circulation model - Formulation and results

A planetary boundary layer (PBL) parameterization for general circulation models (GCMs) is presented. It uses a mixed-layer approach in which the PBL is assumed to be capped by discontinuities in the mean vertical profiles. Both clear and cloud-topped boundary layers are parameterized. Particular emphasis is placed on the formulation of the coupling between the PBL and both the free atmosphere and cumulus convection. For this purpose a modified sigma-coordinate is introduced in which the PBL top and the lower boundary are both coordinate surfaces. The use of a bulk PBL formulation with this coordinate is extensively discussed. Results are presented from a July simulation produced by the UCLA GCM. PBL-related variables are shown, to illustrate the various regimes the parameterization is capable of simulating.

Suarez, M. J.↗
NASA NTRS1990

Preliminary plume characteristics of an arcjet thruster

An experimental investigation of a low power arcjet plume was conducted using emission spectroscopy. A laboratory model arcjet incorporating a segmented anode was run on simulated hydrazine at a flow rate of 5 x 10(exp -5) kg/s. The complete visible spectrum measured in the exit plane of the arcjet showed the presence of N2, N2(+), NH, and H. Radial intensity profiles for the H alpha, H sub beta, and the NH A(sup 3)Pi yields X(sup 3)Sigma(0,0) transitions at four different axial locations were measured. These line of sight intensity measurements, spaced 0.05 mm apart, were deconvoluted to give the radial intensity distribution using an inverse Abel transformation. The ratio between the intensities from the H sub alpha and H sub beta transitions indicated a non-Boltzmann energy distribution between excited states in the plume. Axial intensity profiles taken on center line indicated the decay rate of excited states in the plume. An electron number density of 2 x 10(exp 13)/cu cm at the exit plane was determined based on Stark broadening of the H sub beta line. Rotational temperatures of 750 K, 1750 K, and 2500 K were determined for N2, N2(+), and NH respectively. The results demonstrate that the location of the current attachment on the anode has a measurable effect on the electronically excited species in the plume and that dissociation is the dominant frozen flow loss mechanism in low power arcjets.

Manzella, David H.↗
NASA NTRS1993

A high-sensitivity survey of radio continuum emission from Herbig Ae/Be stars

Results of a high-sensitivity VLA/Australia Telescope survey of radio continuum emission from the 57 Herbig Ae/Be stars and candidates in the 1984 catalog of Finkenzeller and Mundt are presented. Twelve stars were detected at the primary observing wavelength of 3.6 cm, on the basis that not less than 4 sigma radio sources lie within 1 arcsec of the optical positions. It is suggested that the radio emission is predominantly thermal and in many cases wind-related. The unusual eclipsing binary TY CrA is an exception and is classified as a nonthermal radio source on the basis of its decidedly negative spectral index (alpha = -1.2). A simple spherically symmetric free-fall accretion model is used to show that the predicted radio fluxes due to accretion at rates, estimated in the literature, of about 10 exp -6 to 10 exp -5 solar mass/yr are one to four orders of magnitude larger than observed.

Skinner, Stephen L.↗
NASA NTRS1996

Hard Gamma Ray Emission from the Starburst Galaxy NGC 253

We have completed the study to search for hard gamma ray emission from the starburst galaxy NGC 253. Since supernovae are thought to provide the hard gamma ray emission from the Milky Way, starburst galaxies, with their extraordinarily high supernova rates, are prime targets to search for hard gamma ray emission. We conducted a careful search for hard gamma ray emission from NGC 253 using the archival data from the EGRET experiment aboard the CGRO. Because this starburst galaxy happens to lie near the South Galactic Pole, the Galactic gamma ray background is minimal. We found no significant hard gamma ray signal toward NGC 253, although a marginal signal of about 1.5 sigma was found. Because of the low Galactic background, we obtained a very sensitive upper limit to the emission of greater than 100 MeV gamma-rays of 8 x 10(exp -8) photons/sq cm s. Since we expected to detect hard gamma ray emission, we investigated the theory of gamma ray production in a dense molecular medium. We used a leaky-box model to simulate diffusive transport in a starburst region. Since starburst galaxies have high infrared radiation fields, we included the effects of self-Compton scattering, which are usually ignored. By modelling the expected gamma-ray and synchrotron spectra from NGC 253, we find that roughly 5 - 15% of the energy from supernovae is transferred to cosmic rays in the starburst. This result is consistent with supernova acceleration models, and is somewhat larger than the value derived for the Galaxy (3 - 10%). Our calculations match the EGRET and radio data very well with a supernova rate of 0.08/ yr, a magnetic field B approx. greater than 5 x 10(exp -5) G, a density n approx. less than 100/sq cm, a photon density U(sub ph) approx. 200 eV/sq cm, and an escape time scale tau(sub 0) approx. less than 10 Myr. The models also suggest that NGC 253 should be detectable with only a factor of 2 - 3 improvement in sensitivity. Our results are consistent with the standard picture of gamma-ray acceleration by supernovae.

Jackson, James M.↗