Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Spherical Fitting”

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 145 records · Page 8

Precision Effects for Solar Image Coordinates Within the FITS World Coordinate System

The FITS world coordinate system (WCS) provides a number of tools for precisely specifying the spatial coordinates of an image. Many of the finer details that the WCS addresses have not historically been taken into account in solar image processing. This paper examines various effects which can affect the expression of coordinates in FITS headers, to determine under what conditions such effects need to be taken into account in data analysis, and under what conditions they can be safely ignored. Effects which are examined include perspective, parallax, spherical projection, optical axis determination, speed-of-light effects, stellar aberration, gravitational deflection, and scattering and refraction at radio wavelengths. Purely instrumental effects, such as misalignment or untreated optical aberrations, are not considered. Since the value of the solar radius is an experimental quantity, the effect of adopting a specific radius value is also examined. These effects are examined in the context of a previous paper outlining a WCS standard for encoding solar coordinates in FITS files. Aspects of that previous paper are clarified and extended in the present work.

Thompson, W. T.↗

X-Ray Gas Temperatures in the Arc Clusters MS0440+204 and MS0302+1658

The cluster of galaxies MS0440+02, originally discovered through its X-ray emission, was part of an optical observational program to search for arcs and arclets in a complete sample of X-ray luminous, medium-distant clusters of galaxies. Mauna Kea CCD images of MS0440+02 showed a remarkable optical morphology. The core of the cluster contains 6 bright galaxies and numerous fainter ones embedded in a low surface brightness halo. Besides, MS0440+02 is the most spectacular example that we have found of an arc system in a compact condensed cluster, with arcs symmetrically distributed to draw almost perfect circles around the cluster center. Giant arcs are magnified images of distant galaxies, gravitationally distorted by massive foreground clusters. It is of great importance to compare the results of the lensing studies with those derived from X-ray observations, as the two are independent methods of studying the mass distribution. Thus MS0440+02 was the ideal target to obtain temperature measurement with ASCA and good spatial resolution X-ray observations with ROSAT. The X-ray data have been used in conjunction with Hubble Space Telescope observations to put more stringent constrains on the mass estimates. Most of the different wavelength datasets have been reduced and analyzed. Mass determinations have been separately obtained from galaxy virial motions and X-ray profile fitting using the cluster gas temperature as measured by the ASCA satellite. Assuming that the hot gas is in hydrostatic equilibrium and in a spherical potential, we find from the X-ray data a mass distribution profile that is well described by a Beta model. From the multiple images formed by gravitational lensing (HST data) using the modelling of the gravitational lensed arcs, we have derived Beta model. To reconcile the mass estimates we have explored the possibility of having a supercluster surrounding the MOS0440 cluster, that is a model with two isothermal spheres, one embedded inside the other. These results have been published or are in press.

Gioia, Isabella M.↗

Radio emission from nova shells

This paper presents data from multifrequency radio observations of the novae HR Del 1967 and FH Ser 1970 made from June 1970 to December 1978 and of V1500 Cyg 1975 made from September 1975 to December 1978. The data for FH Ser and V1500 Cyg constitute a complete history of the important phases of their radio 'light' curves, while all but the rising portion of the optically thick phase of an expanding nova shell is present for HR Del. A simple solution of the spherically symmetric equations of mass and momentum conservation for an isothermal gas, for which pressure gradients are effectively decoupled from the dynamics, is shown to fit the radio data for all three novae quite well. The main features of this model are strong density and velocity gradients that lead to long periods when the 'observable' shell evolves from an optically thick 'disk' to a completely optically thin shell.

Hjellming, R. M.↗

Simultaneous stochastic inversion for geomagnetic main field and secular variation. I - A large-scale inverse problem

The method of stochastic inversion is extended to the simultaneous inversion of both main field and secular variation. In the present method, the time dependency is represented by an expansion in Legendre polynomials, resulting in a simple diagonal form for the a priori covariance matrix. The efficient preconditioned Broyden-Fletcher-Goldfarb-Shanno algorithm is used to solve the large system of equations resulting from expansion of the field spatially to spherical harmonic degree 14 and temporally to degree 8. Application of the method to observatory data spanning the 1900-1980 period results in a data fit of better than 30 nT, while providing temporally and spatially smoothly varying models of the magnetic field at the core-mantle boundary.

Bloxham, Jeremy↗

Dynamic Finite Element Predictions for Mars Sample Return Cellular Impact Test #4

The nonlinear finite element program MSC.Dytran was used to predict the impact pulse for (he drop test of an energy absorbing cellular structure. This pre-test simulation was performed to aid in the design of an energy absorbing concept for a highly reliable passive Earth Entry Vehicle (EEV) that will directly impact the Earth without a parachute. In addition, a goal of the simulation was to bound the acceleration pulse produced and delivered to the simulated space cargo container. EEV's are designed to return materials from asteroids, comets, or planets for laboratory analysis on Earth. The EEV concept uses an energy absorbing cellular structure designed to contain and limit the acceleration of space exploration samples during Earth impact. The spherical shaped cellular structure is composed of solid hexagonal and pentagonal foam-filled cells with hybrid graphite-epoxy/Kevlar cell walls. Space samples fit inside a smaller sphere at the enter of the EEV's cellular structure. The material models and failure criteria were varied to determine their effect on the resulting acceleration pulse. Pre-test analytical predictions using MSC.Dytran were compared with the test results obtained from impact test #4 using bungee accelerator located at the NASA Langley Research Center Impact Dynamics Research Facility. The material model used to represent the foam and the proper failure criteria for the cell walls were critical in predicting the impact loads of the cellular structure. It was determined that a FOAMI model for the foam and a 20% failure strain criteria for the cell walls gave an accurate prediction of the acceleration pulse for drop test #4.

Fasanella, Edwin L.↗

Dynamic Finite Element Predictions for Mars Sample Return Cellular Impact Test #4

The nonlinear, transient dynamic finite element code, MSC.Dytran, was used to simulate an impact test of an energy absorbing Earth Entry Vehicle (EEV) that will impact without a parachute. EEVOs are designed to return materials from asteroids, comets, or planets for laboratory analysis on Earth. The EEV concept uses an energy absorbing cellular structure designed to contain and limit the acceleration of space exploration samples during Earth impact. The spherical shaped cellular structure is composed of solid hexagonal and pentagonal foam-filled cells with hybrid graphite-epoxy/Kevlar cell walls. Space samples fit inside a smaller sphere at the center of the EEVOs cellular structure. Pre-test analytical predictions were compared with the test results from a bungee accelerator. The model used to represent the foam and the proper failure criteria for the cell walls were critical in predicting the impact loads of the cellular structure. It was determined that a FOAM1 model for the foam and a 20% failure strain criteria for the cell walls gave an accurate prediction of the acceleration pulse for cellular impact.

Fasanella, Edwin L.↗

Axisymmetric model of the ionized gas in the Orion Nebula

New ionization and thermal equilibrium models for the ionized gas in the Orion Nebula with an axisymmetric two-dimensional 'blister' geometry/density distribution are presented. The HII region is represented more realistically than in previous models, while the physical detail of the microphysics and radiative transfer of the earlier spherical modeling is maintained. The predicted surface brightnesses are compared with observations for a large set of lines at different positions to determine the best-fitting physical parameters. The model explains the strong singly ionized line emission along the lines of sight near the Trapezium.

Rubin, R. H.↗

Flexural models of trench/outer rise topography of coronae on Venus with axisymmetric spherical shell elastic plates

Magellan altimetry has revealed that many coronae on Venus have trenches or moats around their peripheries and rises outboard of the trenches. This trench/outer rise topographic signature is generally associated with the tectonic annulus of the corona. Sandwell and Schubert have interpreted the trench/outer rise topography and the associated tectonic annulus around coronae to be the result of elastic bending of the Venus lithosphere (though the tectonic structures are consequences of inelastic deformation of the lithosphere). They used two-dimensional elastic plate flexure theory to fit topographic profiles across a number of large coronae and inferred elastic lithosphere thicknesses between about 15 and 40 km, similar to inferred values of elastic thickness for the Earth's lithosphere at subduction zones around the Pacific Ocean. Here, we report the results of using axisymmetric elastic flexure theory for the deformation of thin spherical shell plates to interpret the trench/outer rise topography of the large coronae modeled by Sandwell and Schubert and of coronae as small as 250 km in diameter. In the case of a corona only a few hundred kilometers in diameter, the model accounts for the small planform radius of the moat and the nonradial orientation of altimetric traces across the corona. By fitting the flexural topography of coronae we determine the elastic thickness and loading necessary to account for the observed flexure. We calculate the associated bending moment and determine whether the corona interior topographic load can provide the required moment. We also calculate surface stresses and compare the stress distribution with the location of annular tectonic features.

Moore, W.↗

Comparing the jerk with other global models of the geomagnetic field from 1960 to 1978

About 3300 satellite values of geomagnetic intensity and about 700 observatory values of annual mean magnetic vector components from 1960 to 1978 were fitted by three global models of the geomagnetic field B. Each model includes a spatially constant external field whose time dependence is a constant plus another constant times the Dst index, and each model accepts a time-independent station correction at each observatory. The time dependence of the internal Gauss coefficients is either cubic, quintic, or biquadratic (two independent quadratics, one before and one after January 1, 1970); and g1(0) also has an induced term proportional to the Dst index. The rms residual of the data fit is the same for the cubic and biquadratic models and insignificantly smaller for the quintic model. The quintic and biquadratic models have 1164 adjustable parameters, and the cubic has 1038. At a high level of significance the parameters of the best fitting biquadratic rule out a physical model for the magnetic impulse of 1969 in which the level surfaces of electrical conductivity in the lower mantle are approximately spherical, and the radial magnetic field at the core-mantle boundary goes from one quadratic time dependence to another in a year or less.

Backus, G. E.↗

Stellar Occultation Probe of Triton's Atmosphere

The goals of this research were (i) to better characterize Triton's atmospheric structure by probing a region not well investigated by Voyager and (ii) to begin acquiring baseline data for an investigation of the time evolution of the atmosphere which will set limits on the thermal conductivity of the surface and the total mass of N2 in the atmosphere. Our approach was to use observations (with the Kuiper Airborne Observatory) of a stellar occultation by Triton that was predicted to occur on 1993 July 10. As described in the attached reprint, we achieved these objectives through observation of this occultation and a subsequent one with the KAO in 1995. We found new results about Triton's atmospheric structure from the analysis of the two occultations observed with the KAO and ground-based data. These stellar occultation observations made both in the visible and infrared, have good spatial coverage of Triton including the first Triton central-flash observations, and are the first data to probe the 20-100 km altitude level on Triton. The small-planet light curve model of Elliot and Young (AJ 103, 991-1015) was generalized to include stellar flux refracted by the far limb, and then fitted to the data. Values of the pressure, derived from separate immersion and emersion chords, show no significant trends with latitude indicating that Triton's atmosphere is spherically symmetric at approximately 50 km altitude to within the error of the measurements. However, asymmetry observed in the central flash indicates the atmosphere is not homogeneous at the lowest levels probed (approximately 20 km altitude). From the average of the 1995 occultation data, the equivalent-isothermal temperature of the atmosphere is 47 +/- 1 K and the atmospheric pressure at 1400 km radius (approximately 50 km altitude) is 1.4 +/- 0.1 microbar. Both of these are not consistent with a model based on Voyager UVS and RSS observations in 1989 (Strobel et al, Icarus 120, 266-289). The atmospheric temperature from the occultation is 5 K colder than that predicted by the model and the observed pressure is a factor of 1.8 greater than the model.

Elliot, James L.↗

The Atmospheric Structure of Triton and Pluto

The goal of this research was to better determine the atmospheric structures of Triton and Pluto through further analysis of three occultation data sets obtained with the Kuiper Airborne Observatory (KAO.) As the research progressed, we concentrated our efforts on the Triton data, as this appeared to be the most fruitful. Three papers have been prepared as a result of this research. The first paper presents new results about Triton's atmospheric structure from the analysis of all ground-based stellar occultation data recorded to date, including one single-chord occultation recorded on 1993 July 10 and nine occultation lightcurves from the double-star event on 1995 August 14. These stellar occultation observations made both in the visible and in the infrared have good spatial coverage of Triton, including the first Triton central-flash observations, and are the first data to probe the altitude level 20-100 km on Triton. The small-planet lightcurve model of J. L. Elliot and L. A. Young was generalized to include stellar flux refracted by the far limb, and then fitted to the data. Values of the pressure, derived from separate immersion and emersion chords, show no significant trends with latitude, indicating that Triton's atmosphere is spherically symmetric at approximately 50 km altitude to within the error of the measurements; however, asymmetry observed in the central flash indicates the atmosphere is not homogenous at the lowest levels probed (approximately 20 km altitude). From the average of the 1995 occultation data, the equivalent isothermal temperature of the atmosphere is 47 plus or minus 1 K and the atmospheric pressure at 1400 km radius (approximately 50 km altitude) is 1.4 plus or minus 0.1 microbar. Both of these are not consistent with a model based on Voyager UVS and RSS observations in 1989. The atmospheric temperature from the occultation is 5 K colder than that predicted by the model and the observed pressure is a factor of 1.8 greater than the model. In our opinion, the disagreement in temperature and pressure is probably due to modeling problems at the microbar level, since measurements at this level have not previously been made. Alternatively, the difference could be due to seasonal change in Triton's atmospheric structure. The second paper reports observations of a recent stellar occultation by Triton which, when combined with earlier results, show that Triton has undergone a period of global warming since 1989. The most conservative estimates of the rate of temperature and surface-pressure increase during this period imply that the atmosphere is doubling in bulk every 10 years -- significantly faster than predicted by published frost model for Triton. Our results suggests that permanent polar caps on Triton play a dominant role in regulating seasonal atmospheric changes. Similar processes should also be active on Pluto. A third paper 'Global Warming on Triton' will appear in a the January 1999 issue of Sky and Telescope.

Elliot, James L.↗

Modeling the Volcanic Source at Long Valley, CA, Using a Genetic Algorithm Technique

In this project, we attempted to model the deformation pattern due to the magmatic source at Long Valley caldera using a real-value coded genetic algorithm (GA) inversion similar to that found in Michalewicz, 1992. The project has been both successful and rewarding. The genetic algorithm, coded in the C programming language, performs stable inversions over repeated trials, with varying initial and boundary conditions. The original model used a GA in which the geophysical information was coded into the fitness function through the computation of surface displacements for a Mogi point source in an elastic half-space. The program was designed to invert for a spherical magmatic source - its depth, horizontal location and volume - using the known surface deformations. It also included the capability of inverting for multiple sources.

Tiampo, Kristy F.↗

DDA Computations of Porous Aggregates with Forsterite Crystals: Effects of Crystal Shape and Crystal Mass Fraction

Porous aggregate grains are commonly found in cometary dust samples and are needed to model cometary IR spectral energy distributions (SEDs). Models for thermal emissions from comets require two forms of silicates: amorphous and crystalline. The dominant crystal resonances observed in comet SEDs are from Forsterite (Mg2SiO4). The mass fractions that are crystalline span a large range from 0.0 < or = fcrystal < or = 0.74. Radial transport models that predict the enrichment of the outer disk (>25 AU at 1E6 yr) by inner disk materials (crystals) are challenged to yield the highend-range of cometary crystal mass fractions. However, in current thermal models, Forsterite crystals are not incorporated into larger aggregate grains but instead only are considered as discrete crystals. A complicating factor is that Forsterite crystals with rectangular shapes better fit the observed spectral resonances in wavelength (11.0-11.15 microns, 16, 19, 23.5, 27, and 33 microns), feature asymmetry and relative height (Lindley et al. 2013) than spherically or elliptically shaped crystals. We present DDA-DDSCAT computations of IR absorptivities (Qabs) of 3 micron-radii porous aggregates with 0.13 < or = fcrystal < or = 0.35 and with polyhedral-shaped Forsterite crystals. We can produce crystal resonances with similar appearance to the observed resonances of comet Hale- Bopp. Also, a lower mass fraction of crystals in aggregates can produce the same spectral contrast as a higher mass fraction of discrete crystals; the 11micron and 23 micron crystalline resonances appear amplified when crystals are incorporated into aggregates composed otherwise of spherically shaped amorphous Fe-Mg olivines and pyroxenes. We show that the optical properties of a porous aggregate is not linear combination of its monomers, so aggregates need to be computed. We discuss the consequence of lowering comet crystal mass fractions by modeling IR SEDs with aggregates with crystals, and the implications for radial transport models of our protoplanetary disk.

Crystal Shape and Crystal Mass Fraction↗

Models of the upper atmosphere

Several computer models of the thermosphere are discussed. J70MM, the Jacchia 1970 model with matrix and matrix mean output, and J703X, the Jacchia 1970, 1971, and 1977 models, are of primary interest. The subroutines in these programs were studied in detail, and several mistakes were found and corrected. It is proposed that a weighted average of the three-hour geomagnetic index be used in the models rather than a single index value. The densities (for a given date and time) generated by the J703X program (all three models), using both a single value of the index and weighted averages for three different lag times, are tabulated and discussed. Most of the equations used in the models are the results of empirical curve fits. An attempt was made to generate a theoretical prediction of the thermospheric temperature profile based on a simple physical model of atmospheric heat conduction in a spherically-symmetric shell. The exponential heating theory and the temperature dependence of k are discussed.

Davis, T. M.↗

NASA Galactic Cosmic Radiation Environment Model: Badhwar-O'Neill (2014)

The Badhwar‐O'Neill (BON) Galactic Cosmic Ray (GCR) flux model is used by NASA to certify microelectronic systems and in the analysis of radiation health risks for human space flight missions. Of special interest to NASA is the kinetic energy region below 4.0 GeV/n due to the fact that exposure from GCR behind shielding (e.g., inside a space vehicle) is heavily influenced by the GCR particles from this energy domain. The BON model numerically solves the Fokker‐Planck differential equation to account for particle transport in the heliosphere due to diffusion, convection, and adiabatic deceleration under the assumption of a spherically symmetric heliosphere. The model utilizes a GCR measurements database from various particle detectors to determine the boundary conditions. By using an updated GCR database and improved model fit parameters, the new BON model (BON14) is significantly improved over the previous BON models for describing the GCR radiation environment of interest to human space flight.

O'Neill, P. M.↗

Line asymmetry in the Seyfert Galaxy NGC 3783

We have reanalyzed the 900 ks Chandra X-ray spectrum of NGC 3783, finding evidence on the asymmetry of the spectral absorption lines. The lines are fitted with a parametric expression that results from an analytical treatment of radiatively driven winds. The line asymmetry distribution derived from the spectrum is consistent with a non-spherical outflow with a finite optical depth. Within this scenario, our model explains the observed correlations between the line velocity shifts and the ionization parameter and between the line velocity shift and the line asymmetry. The present results may provide a framework for detailed testing of models for the dynamic and physical properties of warm absorber in Seyfert galaxies.

Ramirez, J. M.↗

NASA Galactic Cosmic Radiation Environment Model: Badhwar - O'Neill (2014)

The Badhwar-O'Neill (BON) Galactic Cosmic Ray (GCR) flux model has been used by NASA to certify microelectronic systems and in the analysis of radiation health risks for human space flight missions. Of special interest to NASA is the kinetic energy region below 4.0 GeV/n due to the fact that exposure from GCR behind shielding (e.g., inside a space vehicle) is heavily influenced by the GCR particles from this energy domain. The BON model numerically solves the Fokker-Planck differential equation to account for particle transport in the heliosphere due to diffusion, convection, and adiabatic deceleration under the assumption of a spherically symmetric heliosphere. The model utilizes a comprehensive database of GCR measurements from various particle detectors to determine boundary conditions. By using an updated GCR database and improved model fit parameters, the new BON model (BON14) is significantly improved over the previous BON models for describing the GCR radiation environment of interest to human space flight.

Golge, S.↗

The Broad-Band X-Ray Spectrum of IC 4329A from a Joint NuSTAR/Suzaku Observation

We have obtained a deep, simultaneous observation of the bright, nearby Seyfert galaxy IC 4329A with Suzaku and NuSTAR (Nuclear Spectroscopic Telescope Array). Through a detailed spectral analysis, we are able to robustly separate the continuum, absorption, and distant reflection components in the spectrum. The absorbing column is found to be modest (approximately 6 times 10 (sup 21) per square centimeter), and does not introduce any significant curvature in the Fe K band. We are able to place a strong constraint on the presence of a broadened Fe K alpha line (E (sub rest) equals 6.46 (from plus 0.08 to minus 0.07) kiloelectronvolts with sigma equals 0.33 (from plus 0.08 to minus 0.07) kiloelectronvolts and EW equals 34 (from plus 8 to minus 7) electronvolts), though we are not able to constrain any of the parameters of a relativistic reflection model. These results highlight the range in broad Fe K line strengths observed in nearby, bright, active galactic nuclei (roughly an order of magnitude), and imply a corresponding range in the physical properties of the inner accretion disk in these sources. We have also updated our previously reported measurement of the high-energy cutoff of the hard X-ray emission using both observatories rather than just NuSTAR alone: E (sub cut) equals 186 plus or minus 14 kiloelectronvolts. This high-energy cutoff acts as a proxy for the temperature of the coronal electron plasma, enabling us to further separate this parameter from the plasmas optical depth and to update our results for these parameters as well. We derive kT 50 (from plus 6 to minus 3) kiloelectronvolts with tau equals 2.34 (from plus +0.16 to minus 0.11) using a spherical geometry, kT equals 61 plus or minus1 kiloelectronvolt with tau equals 0.68 plus or minus 0.02 for a slab geometry, with both having an equivalent goodness-of-fit.

Dynamics - stars↗