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.

286 records · Page 16

Toward a NASA Deep Space Optical Communications System

As discussed at SpaceOps in 2016, we expect the data rates from deep space missions to increase approximately one order of magnitude per decade for the next 50 years. The first order of magnitude improvement will come from existing plans for radio frequency (RF) communications including enhancements to both spacecraft and Deep Space Network (DSN) facilities. The next two orders of magnitude are predicted to come from the introduction of deep space optical communications. Studies indicate that optical receive apertures of between 8m12m are desired. The large cost of dedicated receive telescopes makes this method unrealistic – at least in the near-term. The cost of large optical ground terminals is driven primarily by the cost of the optics and by the cost of a stable structure for the telescope. We propose a novel hybrid design in which existing DSN 34m beam waveguide (BWG) radio antennas can be modified to include an 8m equivalent optical primary. By utilizing a low-cost segmented spherical mirror optical design, pioneered by the optical astronomical community, and by exploiting the already existing extremely stable large radio aperture structures in the DSN, we can minimize both of these cost drivers for implementing large optical communications ground terminals. Two collocated hybrid RF/optical antennas could be arrayed to synthesize the performance of an 11.3m receive aperture to support more capable or more distant space missions or used separately to communicate with two optical spacecraft simultaneously. NASA is in the midst of building six new 34m BWG antennas in the DSN. The final two are planned to be built at the DSN Goldstone, California and Canberra complexes. We are now investigating building these last two antennas as RF/optical hybrids. By delaying their operational dates by two years, we would be able to add the 8m optical receive capability for these two antennas while fitting within existing budgetary constraints. This paper describes the hybrid antenna design, the technical challenges being addressed, and plan for using this concept, together with ongoing work on optical flight terminals, to infuse operation optical communications into deep space missions.

Cornwell, Donald M.↗

High-degree Gravity Models from GRAIL Primary Mission Data

We have analyzed Ka‒band range rate (KBRR) and Deep Space Network (DSN) data from the Gravity Recovery and Interior Laboratory (GRAIL) primary mission (1 March to 29 May 2012) to derive gravity models of the Moon to degree 420, 540, and 660 in spherical harmonics. For these models, GRGM420A, GRGM540A, and GRGM660PRIM, a Kaula constraint was applied only beyond degree 330. Variance‒component estimation (VCE) was used to adjust the a priori weights and obtain a calibrated error covariance. The global root‒mean‒square error in the gravity anomalies computed from the error covariance to 320×320 is 0.77 mGal, compared to 29.0 mGal with the pre‒GRAIL model derived with the SELENE mission data, SGM150J, only to 140×140. The global correlations with the Lunar Orbiter Laser Altimeter‒derived topography are larger than 0.985 between l = 120 and 330. The free‒air gravity anomalies, especially over the lunar farside, display a dramatic increase in detail compared to the pre‒GRAIL models (SGM150J and LP150Q) and, through degree 320, are free of the orbit‒track‒related artifacts present in the earlier models. For GRAIL, we obtain an a posteriori fit to the S‒band DSN data of 0.13 mm/s. The a posteriori fits to the KBRR data range from 0.08 to 1.5 micrometers/s for GRGM420A and from 0.03 to 0.06 micrometers/s for GRGM660PRIM. Using the GRAIL data, we obtain solutions for the degree 2 Love numbers, k20=0.024615+/-0.0000914, k21=0.023915+/-0.0000132, and k22=0.024852+/-0.0000167, and a preliminary solution for the k30 Love number of k30=0.00734+/-0.0015, where the Love number error sigmas are those obtained with VCE.

Primary 1↗

Crustal Magnetic Spectra from Correlated Sources on Mars

The spectral method for distinguishing crustal from core-source magnetic fields has been re-examined, modified and applied to both a comprehensive geomagnetic field model and an altitude normalized magnetic map of Mars. These observational spectra are fairly fitted by theoretical forms expected from certain elementary classes of magnetic sources. For Earth we found fields from a core radius 3512 +/- 64 km and a crust represented by a shell of random dipolar sources at radius 6367 +/- 14 km. For Mars we found only a field from a crust represented in the same way, but 46 +/- 10 km below the planetary menu radius. More realistic theoretical spectra, allowing for crustal thickness, oblateness and magnetization by a planet centered dipole, were derived and discussed, as were spectral effects of laterally correlated sources. The main effect of laterally correlated sources is to soften the spectrum at high degrees. We tend to over-estimate source shell depth when this is omitted. To include this effect simply size and magnetization distribution functions for extended sources are recast as a characteristic diameter and mean square magnetization amplitude for an ensemble of vertically magnetized spherical caps on a shell. For small caps, and at moderated degrees, the practical derivatives of the log-theoretical spectrum with respect to amplitude, shell radius, and cap diameter are approximately proportional to l, n, and -n **2, respectively. Separation of diameter from amplitude from layer thickness. Results from applications to observational spectra are discussed, noting that there are now several fine field models for Mars; moreover, the terrestrial magnetic spectrum at high degrees, as revised, updated with high precision Oerstad data, and upgrade with high resolution Champ data, appears softer than before.

Voorhies, Coerte V.↗

Modelling CLPX IOP3 Radiometric Data by Means of the Dense Media Theory: Preliminary Results for the LSOS Test Site

The capabilities of the Dense Media Radiative Transfer model using the Quasi Crystalline Approximation with Coherent Potential (QCA-CP) to reproduce measured radiometric data were tested using the University of Tokyo Ground Based Microwave Radiometer (GBMR-7) during the third Intensive Observation Period (IOP3) of the NASA Cold-land Processes Field Experiment (CLPX). The data were collected at the Local-Scale Observation Site (LSOS), a 0.8-ha study site consisting of two open meadows separated by trees. Intensive measurements were also made of snow depth and temperature, density, and grain size profiles. A DMRT model is needed to describe radiative transfer in a medium such as snow because the assumption of independent scattering used in classical radiative transfer theory (CRT) is not valid. Validation of the DMRT approach requires a relationship between measured snow grain size and the DMRT approximation of snow grain radius as spherical particles with a mean radius of the log-normal particle-size distribution. This relationship is very important for a better understanding of snow modelling and for practical applications. DMRT simulations were compared with observations of microwave brightnesses at 18.7, 36.5 and 89-GHz (V and H polarizations) collected on February-1 9-25, 2003. Observation angles ranged from 30\deg to 70\deg. Model inputs included measured snow parameters except mean grain size. The average snow temperature, fractional volume and depth were held constant, together with the ice and soil permittivities. The minimum and maximum measured mean grain sizes were used to test the capabilities of the DMRT to reproduce the brightnesses as upper and lower limits. The sensitivity to the largest and smallest measured grain size in the three classes of minimum, medium and maximum observed grain sizes was also investigated. DMRT particle sizes yielding a best-fit to the experimental data for each date were computed. Results show that the measured brightnesses fall within the range of simulated brightnesses using the smallest and largest measured grain size values. The DMRT best-fit radii are comparable to the average radii for the medium observed grain sizes.

Tedesco, Marco↗

Effect of surface roughness on characteristics of spherical shock waves

An investigation has been conducted on a small-scale test layout in which direct observation of the shock wave movement with time could be made in order to determine the effects of surface roughness on the characteristics of spherical shock waves. Data were obtained with 15-gram pentolite charges at four heights of burst, both for a smooth surface and for a surface completely covered with pyramid-shaped roughness elements. The observations resulted in determinations of shock peak overpressure and Mach stem height as a function of distance for each test. Comparison of the smooth-surface data with those obtained for the extremely rough condition showed a small net effort of roughness on the shock peak overpressures at the surface for all burst heights, the effect being to lower the overpressures. The effect of surface roughness on the Mach stem formation and growth was to delay the formation at the greatest charge height and to lower the height of the Mach stem for all heights.Comparison of the free-air shock peak overpressures with larger scale data showed good similarity of the overpressure-distance relationships. The data did not fit a geometrical similarity parameter for the path of the triple point at different heights of burst suggested by other investigators. A simple similarity parameter (relating the horizontal distance to the theoretical point of Mach formation) was found which showed only a small influence of burst height on the path of the triple point. While the data presented provide knowledge of the effect of many surface-roughness elements on the overall shock characteristics, the data do not provide insight into the details of the air-flow characteristics along the surface, nor the relative contribution of individual roughness elements to the results obtained.

Huber, Paul W↗

Lunar Gravity Field Determination Using SELENE Same-Beam Differential VLBI Tracking Data

A lunar gravity field model up to degree and order 100 in spherical harmonics, named SGM 100i, has been determined from SELENE and historical tracking data, with an emphasis on using same-beam S-band differential VLBI data obtained in the SELENE mission between January 2008 and February 2009. Orbit consistency throughout the entire mission period of SELENE as determined from orbit overlaps for the two sub-satellites of SELENE involved in the VLBI tracking improved consistently from several hundreds of metres to several tens of metres by including differential VLBI data. Through orbits that are better determined, the gravity field model is also improved by including these data. Orbit determination performance for the new model shows improvements over earlier 100th degree and order models, especially for edge-on orbits over the deep far side. Lunar Prospector orbit determination shows an improvement of orbit consistency from I-day predictions for 2-day arcs of 6 m in a total sense, with most improvement in the along and cross-track directions. Data fit for the types and satellites involved is also improved. Formal errors for the lower degrees are smaller, and the new model also shows increased correlations with topography over the far side. The estimated value for the lunar GM for this model equals 4902.80080 +/- 0.0009 cu km/sq s (10 sigma). The lunar degree 2 potential Love number k2 was also estimated, and has a value of 0.0255 +/- 0.0016 (10 sigma as well).

Goossens, S.↗

X-ray spectrum and variability of 3C 120

The X-ray luminosity of 3C 120 is variable by a factor of 2.5 on time scales of days to months. The spectral slope of the 2-10 keV X-rays changes systematically in the sense that the higher intensity states are steeper. Simple power laws fitted to the data are consistent with an energy index alpha of 0.51 + or - 0.13 at minimum intensity and 0.83 + or - 0.03 at maximum, with a continuous variation in between. A more complicated spectral shape cannot be ruled out, but changes in absorbing column density alone cannot explain the data. The results are interpreted in terms of beamed synchrotron emission from a relativistic jet, a model which has been invoked previously for the spectra of BL Lac objects. The spectral variability is in marked contrast to the behavior of 3C 273, another superluminal source in which the X-ray spectral index is constant to a high degree of accuracy despite variability in time scales and amplitudes which are quite similar to those of 3C 120. The difference may be due to the degree of beaming of the respective core components. Another model which is consistent with the data involves thermal emission from isothermal shocks in a hot (approximately 10 to the 8th K) medium. Thermal Compton and spherical accretion X-ray source models do not reproduce the correct behavior of alpha versus flux.

Halpern, J. P.↗

Graphite-fiber-reinforced polyimide liners of various compositions in plain spherical bearings

A plain spherical bearing design with a ball diameter of 28.6 mm, a race length of 12.7 mm, and a 1.7-mm-thick, molded composite liner was evaluated. The liner material is a self-lubricating composite of graphite-fiber-reinforced polyimide resin (GFRPI). The liner is prepared by transfer molding a mixture of one part chopped graphite fiber and one part partially polymerized resin into the space between the bearing ball and the outer race and then completing the polymerization under heat and pressure. Several liner compositions were evaluated: two types of polyimide, condensation and addition; two types of graphite fiber, low and high modulus; and four powder additives - cadmium oxide, cadmium iodide, graphite fluoride, and molybdenum disulfide. The bearings were oscillated + or - 15 deg at 1 Hz for 20 kilocycles under a radial unit load of 29 MN sq m (4200 psi) in dry air at 25, 200, or 315 C. Both types of fiber and polyimide gave low friction and wear. A simple equation was developed to fit the wear-time data and adequately predicted wear to 100 kilocycles.

Sliney, H. E.↗

Reanalysis of x ray emission from M87 II: Multiphase models

In a previous paper by the authors (hereafter TB), it was shown that the hot intracluster gas around M87 could not be adequately explained in terms of a spherically symmetric, single phase model. It was found that although data from the Einstein satellite High Resolution Imager (HRI), Imaging Proportional Counter (IPC), and Focal Point Crystal Spectrometer (FPCS) could indeed be simultaneously explained by a single phase model, data from the Solid State Spectrometer (SSS) and optically determined mass estimates could not be similarly explained. A qualitative discussion was then given to indicate how the adoption of a specific multiphase model could plausibly improve upon the single phase model, but a detailed discussion was delayed until now. The x ray data is reanalyzed for the gas around M87 in the same spirit as the work in TB. That is, by assuming several cooling flow models, the shortcomings are improved upon. The procedure of TB was adopted and the data reanalyzed in the following manner: a multiphase model is assumed which can compute the surface brightness as seen by HRI and IPC. Also, the line fluxes can be computed of the lines seen by the FPCS and the mass profile after assuming hydrostatic equilibrium. The parameters of the model are then adjusted to fit the data. A check is then made to see whether the resulting model is consistent with the SSS spectrum and the equivalent width of the 7 keV complex of Fe lines as seen by large field of view instruments.

Tsai, John C.↗

Long-term modulation of galactic cosmic radiation and its model for space exploration

As the human exploration of space has received new attention in the United States, studies find that exposure to space radiation could adversely impact the mission design. Galactic Cosmic Radiation (GCR), with its very wide range of charges and energies, is particularly important for a mission to Mars, because it imposes a stiff mass penalty for spacecraft shielding. Dose equivalent versus shielding thickness calculations, show a rapid initial drop on exposure with thickness, but an asymptotic behavior at a higher shielding thickness. Uncertainties in the radiobiology are largely unknown. For a fixed radiation risk, this leads to large uncertainties in shielding thickness for small uncertainties in estimated dose. In this paper we investigate the application of steady-state, spherically-symmetric diffusion-convection theory of solar modulation to individual measurements of differential energy spectra from 1954 to 1989 in order to estimate the diffusion coefficient, kappa(r,t), as a function of time. We have correlated the diffusion coefficient to the Climax neutron monitor rates and show that, if the diffusion coefficient can be separated into independent functions of space and time: kappa(r,t)=K(t) k(sub 0) beta Pk(sub 1) (r), where beta is the particle velocity and P the rigidity, then (i) The time dependent quantity 1/K(t), which is proportional to the deceleration potential, phi(r,t), is linearly related to the Climax neutron monitor counting rate. (ii) The coefficients obtained from hydrogen or helium intensity measurements are the same. (iii) There are different correlation functions for odd and even solar cycles. (iv) The correlation function for the Climax neutron monitor counting rate for given time, t, can be used to estimate mean deceleration parameter phi(t) to within +/- 15% with 90% confidence. We have shown that k(r,t) determined from hydrogen and/or helium data, can be used to fit the oxygen and iron differential energy spectra with a root mean square error of about +/- 10%, and essentially independent of the particle charge or energy. We have also examined the ion chamber and C-14 measurements which allow the analysis to be extended from the year 1906 to 1990. Using this model we have defined reference GCR spectra at solar minimum and solar maximum. These can be used for space exploration studies and provide a quantitative estimate of the error in dose due to changes in GCR intensities.

Badhwan, G. D.↗

Flow Regime Transitions and Pressure Drop in Gas-Liquid Flows through Packed Beds in Microgravity

The Packed Bed Reactor Experiment (PBRE) focused on the hydrodynamics of gas-liquid flow in porous media, specifically on flow regime transition and pressure drop correlation. The reactor included a cylindrical test section of 60 cm in length, 5 cm in diameter and filled with 2 mm spherical glass beads. Liquid water and Nitrogen gas were flowed into a centrifugal phase separator after exiting the test section. Five absolute pressure transducers and two high-speed video cameras spanned the test section to measure the pressure drop and capture flow regime behavior. Video observations, analysis of pressure traces, and pressure drop data were used to identify flow patterns. It was found that the slope of pressure gradient data when plotted versus the gas or liquid Reynolds numbers changed when the flow regime changed. This observation was used to identify the flow regime boundaries. At high liquid flow rates (high interaction regime), dispersed bubble and pulse flow regimes were detected. At low liquid flow rates (low interaction regime), large bubble and gas channeling regimes were observed. The transition from the big bubble regime to the gas channeling regime occurred at a critical gas flow rate independent of the liquid flow rate. Similarly, the transition from the low interaction to the high interaction regime occurred at a critical liquid flow rate. For these different flow regimes, different two-phase flow friction factors and pressure gradient correlations were fitted as a function of the gas and liquid modified Reynolds numbers and Suratman number. It was found that in all the flow regimes the capillary contribution to the pressure drop is significant and can be as high as 90% of the total pressure drop.

Brian Motil↗

New Features of the NEQAIR Radiation Code

The longest-lived code for predicting shock layer radiation, NEQAIR, is now in its 5th decade of service. Substantial changes to the code have been made over the previous decade, the most recent report of which was at the 5th Workshop on Radiation in High Temperature Gases in 2014, for the version referred to as NEQAIR14. This paper will review some of the improvements made to the NEQAIR code since then, which is now at v15.2. Some of these features are discussed briefly below. NEQAIR15 and subsequent versions have enabled parallel evaluation of multiple lines of sight. This is accomplished by utilizing the HDF5 file format and placing multiple lines into a single file, LOS.h5, which is used for both input and output. This approach enables straightforward parallel execution both over the number of lines of sight and the number of points per line. For large problems, runtime reduces linearly with the number of nodes deployed since each line is processed independently by a subset of MPI ranks. Three applications of the multi-line solver are discussed. The first has to do with performing loosely coupled radiation-flowfield solutions. In this case the computed absorption and emission coefficients are used to evaluate the total energy absorbed or emitted at each point, allowing evaluation of the volumetric source term in the flowfield. The second computation is for obtaining heat flux from nonuniform flows, which require integration over spherical co-ordinates. These are of particular interest for evaluating radiation on the vehicle backshell. This 3D option improves the angular integration scheme and allows adaptive line selection that together reduce the number of lines required by about an order of magnitude. The final application is for remote observation, which is essentially the 3D integration problem over a small solid angle. For all three of these computations, data can be stored in the HDF5 file which allows a NEQAIR run to be restarted when it times out, or to add atmospheric absorption or instrument scan functions. An additional level of parallelism is enabled in NEQAIR15.2 using GPU routines. The GPU parallelism has realized up to 8x speed-up when running on a single core but diminishes as CPU parallelism is increased. For running multi-line simulations, it may be easier to reserve a large number of CPU nodes than to obtain the number of GPU nodes required for similar performance. A GUI, known as NEQTPY, allows for reading and creating input files, running NEQAIR, and displaying results. A significant feature of NEQTPY is the ability to perform spectral fits to data. The fits can operate on a single line spectrum (radiance vs. wavelength) or a 3D input file with multiple columns of data. Other new features include improved constants, additional species, more detailed non-Boltzmann modelling, advanced user controls, the ability to read and calculate spectra from HITRAN datafiles, photodissociation and photoionization cross-sections. A “fast” automatic grid option may reduce the size and time of spectral calculations while still maintaining good accuracy for total heat flux.

Brett A Cruden↗

Using FU Orionis outbursts to constrain self-regulated protostellar disk models

One-dimensional, convective, vertical structure models and one dimensional time-dependent, radial diffusion models are combined to create a self-consistent picture in which FU Orionis outbursts occur in young stellar objects (YSOs) as the result of a large-scale, self-regulated, thermal ionization instability in the surrounding protostellar accretion disk. Although active accretion disks have long been postulated to be ubiqitous among low-mass YSOs, few constraints have until now been imposed on physical conditions in these disks. By fitting the results of time-dependent disk models to observed timescales of FU Orionis events, we estimate the magnitude of the effective viscous stress in the inner disk (r approximately less than 1 AU) to be, in accordance with an ad hoc 'alpha' prescription, the product of the local sound speed, pressure scale height, and an efficiency factor alpha of 10(exp -4) where hydrogen is neutral and 10(exp 3) where hydrogen is ionized. We hypothesize that all YSOs receive infall onto their outer disks which is steady (or slowly declining with time) and that FU Orionis outbursts are self-regulated, disk outbursts which occur only in systems which transport matter inward at a rate sufficiently high to cause hydrogen to be ionized in the inner disk. We estimate a critical mass flux of dm(sub crit)/dt = 5 x 10(exp 7) solar mass/yr independent of the magnitude of alpha for systems with one solar mass, three solar radius central objects. Infall accretion rates in the range of dm(sub in)/dt = 1-10) x 10(exp -6) solar mass/yr produce observed FU Orionis timescales consistent with estimates of spherical molecular cloud core collapse rates. Modeled ionization fronts are typically initiated near the inner edge of the disk and propogate out to a distance of several tens of stellar radii. Beyond this region, the disk transports mass steadily inward at the supplied constant infall rate. Mass flowing through the innermost disk annulus is equal to dm(sub in)/dt only in a time-averaged sense and is regulated by the ionization of hydrogen in the inner disk such that long intervals (approximately 1000 yr) of low-mass flux: (1-30) x 10(exp -8) solar mass/yr are punctuated by short intervals (approximately 100 yr) of high-mass flux: (1-30) x 10(exp -5) solar mass/yr. Timescales and mass fluxes derived for quiescent and outburst stages are consistent with estimates from observations of T Tauri and FU Orionis systems, respectively.

Bell, K. R.↗

Impact experiments in viscous fluid media

Available phase and group velocity data are inverted by a procedure which includes the effects of transverse anisotropy, anelastic dispersion, sphericity, and gravity. The resulting models, for average Earth, average ocean, and oceanic regions divided according to the age of the ocean floor, are quite different from previous results which ignore the above effects. The models show a low-veocity zone with age dependent anisotropy and velocities higher than derived in previous surface wave studied. The correspondence between the anisotropy variation with age and a physical model based on flow aligned olivine is suggestive. For most of the Earth SHSV in the vicinity of the low-velocity zone. Near the East Pacific Rise, however, SVSH at depth, consistent with ascending flow. Anisotropy is as important as temperature in causing radial and lateral variations in velocity. The models have a high velocity nearly isotropic layer at the top of the mantle that thickens with age. This layer defines the LID, or seismic lithosphere. In the Pacific, the LID thickens with age to a maximum thickness of about 50 km. This thickness is comparable to the thickness of the elastic lithosphere. The LID thickness is thinner than derived using isotropic or pseudo-isotropic procedures A new model for Average Earth is obtained which includes a thin LID. This model extends the fit of a P.R.E.M. type model to shorter period surface waves.

Greeley, R.↗

A Model of Thermal Conductivity for Planetary Soils: Theory for Cemented Soils - 2

A numerical model of heat conduction through particulate media made of spherical grains cemented by various bonding agents is presented. The pore-filling gas conductivity, volume fraction, and thermal conductivity of the cementing phase are tunable parameters. Cement fractions <0.001-0.01% in volume have small effects on the soil bulk thermal conductivity. A significant conductivity increase (factor 3-8) is observed for bond fractions of 0.01 to 1% in volume. In the 1 to 15% bond fraction domain, the conductivity increases continuously but less intensely (25-100% conductivity increase compared to a 1% bond system). Beyond 15% of cements, the conductivity increases vigorously and the bulk conductivity rapidly approaches that of bedrock. The composition of the cements (i.e. conductivity) has little influence on the bulk thermal inertia of the soil, especially if the volume of bond <10%. These results indicate that temperature measurements are sufficient to detect cemented soils and quantify the amount of cementing phase, but the mineralogical nature of the bonds and the typical grain size are unlikely to be determined from orbit. On Mars, a widespread surface unit characterized by a medium albedo (0.19-0.26) and medium/high thermal inertia (200-600 J s(0.5)/sq m/K) has long been hypothesized to be associated with a duricrust. The fraction of cement required to fit the thermal data is less than approx.1-5% by volume. This small amount of material is consistent with orbital observations, confirming that soil cementation is an important factor controlling the thermal inertia of the Martian surface

Piqueux, S.↗

Velocity-Resolved Far-Infrared Spectra of [Fe 2]: Evidence for Mixing and Clumping in SN 1987A

We present approx. 400 km/s resolution profiles of the 17.94 and 25.99 micron [Fe II] transitions from SN 1987A at t approx. 400 days after core collapse. These observations used the facility cooled grating spectrometer aboard NASA's Kuiper Airborne Observatory. The two profiles are similar and have FWHM line widths of approx. 2700 km/s. The higher signal-to-noise 18 micron profile is somewhat asymmetric, falling off more steeply on the redshifted side than on the blue. Gaussian fits to the profiles yield an average centroid velocity of 280 +/- 140 km/s relative to the Large Magellanic Cloud. The wings of the profiles extend to velocities is approx. greater than 3000 km/s. This shows that a significant fraction of the iron has been mixed outward into the hydrogen-rich envelope, which has a minimum expansion velocity of 2100-2400 km/s. Both profiles also contain an unresolved 3-5 sigma emission feature on the redshifted wing at nu(LSR) approx. + 3900 km/s. We interpret this feature as emission from a high-velocity clump of material containing approx. 3% of the total iron mass. The total line flux of the 26 micron ground-state transition yields an optically thin, singly ionized iron mass of 0.026 solar mass, relatively independent of the assumed temperature. This is significantly less than the 0.06 Me of Fe+ determined from the decline of the optical light curve and the ionization of measured nickel lines, implying that the iron transitions still have appreciable optical depth. However, because of the small change in the 26 micron line flux from our measurement at 250 days, and the similarity of our profiles to the 1.26 micron [Fe II] profile, most of the emission is believed to originate from optically thin material with a temperature of 4406 +/- 400 K. A comparison of the data with spherically symmetric models indicates a power-law density exponent of -3.2 +/- 1.1 and a minimum expansion velocity of 650 +/- 650 km/s for this optically thin component. The [Fe II] line fluxes and profiles also imply that the remainder of the material has high optical depth and is distributed in clumps throughout the ejecta, rather than being concentrated at low velocities in the center of a smooth density distribution.

Haas, Michael R.↗