Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “local density approximations”

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.

729 records · Page 41
All sourcesNASA NTRSDOE OSTI
NASA NTRS2017

New Global Electron Density Observations from GPS-RO in the D- and E-Region Ionosphere

A novel retrieval technique is developed for electron density (N(sub e)) in the D- and E-region (80-120 km) using the high-quality 50-Hz GPS radio occultation (GPS-RO) phase measurements. The new algorithm assumes a slow, linear variation in the F-region background when the GPS-RO passes through the D- and E-region, and extracts the N(sub e) profiles at 80-130 km from the phase advance signal caused by N(sub e). Unlike the conventional Abel function, the new approach produces a sharp N(sub e) weighting function in the lower ionosphere, and the N(sub e) retrievals are in good agreement with the IRI (International Reference Ionosphere) model in terms of monthly maps, zonal means and diurnal variations. The daytime GPS-RO N(sub e) profiles can be well characterized by the alpha-Chapman function of three parameters (N(sub mE), h(sub mE) and H), showing that the bottom of E-region is deepening and sharpening towards the summer pole. At high latitudes the monthly GPS-RO N(sub e) maps at 80-120 km reveal clear enhancement in the auroral zones, more prominent at night, as a result of energetic electron precipitation (EEP) from the outer radiation belt. The D-/E-region auroral N(sub e) is strongly correlated with K(sub p) on a daily basis. The new N(sub e) data allow further comprehensive analyses of the sporadic E (E(sub s)) phenomena in connection with the background N(sub e) in the E-region. The layered (2-10 km) and fluctuated (less than 2 km) E(sub s) components, namely N(sub e_Layer) than N(sub e_Pert), are extracted with respect to the background N( sub e_Region) on a profile-by-profile basis. The N(sub e_Layer) component has a strong but highly-refined peak at approximately 105 km, with an amplitude smaller than N(sub e_Region) approximately by an order of magnitude. The N(sub e_Pert) component, which was studied extensively in the past, is approximately 2 orders of magnitude weaker than N(sub e_Layer). Both N(sub e_Layer) and N(sub e_Pert) are subject to significant diurnal and semidiurnal variations, showing downward progression with local time in amplitude. The 11-year solar cycle dominates the N(sub e) interannual variations, showing larger N(sub e_Region) and N(sub e_Layer) but smaller N(sub e_Pert) amplitudes in the solar maximum years. Enhanced Ne profiles are often observed in the polar winter, showing good correlation with solar proton events (SPEs) and geomagnetic activity. The new methodology offers great potential for retrieving low N(sub e) in the D-region, where radio propagation and communication blackouts can occur due to enhanced ionization. For space weather applications it is recommended for GPSRO operations to raise the top of high-rate data acquisition to approximately 140 km in the future.

GSP-R↗
DOE OSTI2024

How well do one-electron self-interaction-correction methods perform for systems with fractional electrons?

Recently developed locally scaled self-interaction correction (LSIC) is a one-electron SIC method that, when used with a ratio of kinetic energy densities (z σ ) as iso-orbital indicator, performs remarkably well for both thermochemical properties as well as for barrier heights overcoming the paradoxical behavior of the well-known Perdew–Zunger self-interaction correction (PZSIC) method. In this work, we examine how well the LSIC method performs for the delocalization error. Our results show that both LSIC and PZSIC methods correctly describe the dissociation of $H$$^{+}_{2}$ and $H$$^{+}_{2}$ but LSIC is overall more accurate than the PZSIC method. Likewise, in the case of the vertical ionization energy of an ensemble of isolated He atoms, the LSIC and PZSIC methods do not exhibit delocalization errors. For the fractional charges, both LSIC and PZSIC significantly reduce the deviation from linearity in the energy vs number of electrons curve, with PZSIC performing superior for C, Ne, and Ar atoms while for Kr they perform similarly. The LSIC performs well at the endpoints (integer occupations) while substantially reducing the deviation. The dissociation of LiF shows both LSIC and PZSIC dissociate into neutral Li and F but only LSIC exhibits charge transfer from Li + to F – at the expected distance from the experimental data and accurate ab initio data. Overall, both the PZSIC and LSIC methods reduce the delocalization errors substantially.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2010

Ocean Turbulence, III: New GISS Vertical Mixing Scheme

We have found a new way to express the solutions of the RSM (Reynolds Stress Model) equations that allows us to present the turbulent diffusivities for heat, salt and momentum in a way that is considerably simpler and thus easier to implement than in previous work. The RSM provides the dimensionless mixing efficiencies Gamma-alpha (alpha stands for heat, salt and momentum). However, to compute the diffusivities, one needs additional information, specifically, the dissipation Epsilon. Since a dynamic equation for the latter that includes the physical processes relevant to the ocean is still not available, one must resort to different sources of information outside the RSM to obtain a complete Mixing Scheme usable in OGCMs. As for the RSM results, we show that the Gamma-alpha s are functions of both Ri and Rq (Richardson number and density ratio representing double diffusion, DD); the Gamma-alpha are different for heat, salt and momentum; in the case of heat, the traditional value Gamma-h = 0.2 is valid only in the presence of strong shear (when DD is inoperative) while when shear subsides, NATRE data show that Gamma-h can be three times as large, a result that we reproduce. The salt Gamma-s is given in terms of Gamma-h. The momentum Gamma-m has thus far been guessed with different prescriptions while the RSM provides a well defined expression for Gamma-m(Ri,R-rho). Having tested Gamma-h, we then test the momentum Gamma-m by showing that the turbulent Prandtl number Gamma-m/Gamma-h vs. Ri reproduces the available data quite well. As for the dissipation epsilon, we use different representations, one for the mixed layer (ML), one for the thermocline and one for the ocean;s bottom. For the ML, we adopt a procedure analogous to the one successfully used in PB (planetary boundary layer) studies; for the thermocline, we employ an expression for the variable epsilon/N(exp 2) from studies of the internal gravity waves spectra which includes a latitude dependence; for the ocean bottom, we adopt the enhanced bottom diffusivity expression used by previous authors but with a state of the art internal tidal energy formulation and replace the fixed Gamma-alpha = 0.2 with the RSM result that brings into the problem the Ri, R-rho dependence of the Gamma-alpha; the unresolved bottom drag, which has thus far been either ignored or modeled with heuristic relations, is modeled using a formalism we previously developed and tested in PBL studies. We carried out several tests without an OGCM. Prandtl and flux Richardson numbers vs. Ri. The RSM model reproduces both types of data satisfactorily. DD and Mixing efficiency Gamma-h(Ri,Rq). The RSM model reproduces well the NATRE data. Bimodal epsilon-distribution. NATRE data show that epsilon (Ri < 1) approximately equals 10epsilon(Ri > 1), which our model reproduces. Heat to salt flux ratio. In the Ri much greater than 1 regime, the RSM predictions reproduce the data satisfactorily. NATRE mass diffusivity. The z-profile of the mass diffusivity reproduces well the measurements at NATRE. The local form of the mixing scheme is algebraic with one cubic equation to solve.

Canuto, V. M.↗
DOE OSTI2020

Best Practices for Shielding Analyses of Activated Metals and Spent Resins from Reactor Operation

This report investigated best practices for performing shielding evaluations of Type B waste packages, as defined in 10 CFR Part 71 regulations on packaging and transportation of radioactive material, or packages for which the contents are not defined before loading and may include a broad range of nuclides, geometries, and non-fuel materials. The following non-fuel waste streams were analyzed: Activated metals from decommissioned commercial power reactors, including Type 304 stainless steel, reactor vessel steel, and Inconel, Control blades from boiling water reactors, Neutron-activated corrosion products on surfaces of activated metals, and Spent resins from power plant operations. Measured elemental compositions, including major constituents and impurities, for steel and Inconel samples from commercial power reactors were used in activation calculations to determine radionuclide inventories in activated metals. For a simplified cask model, 60 Co contribution to the total external package dose rate at 30 days after shutdown varied from approximately 60% to 95%, depending on the activated metal, initial cobalt impurity concentration in the metal, and the thickness of the overpack gamma shield. Its maximum contribution to the total external package dose rate of approximately 100% was reached within the time interval of 2 to 5 years after shutdown and was maintained for up to 45 to 60 years after shutdown, depending on material, initial cobalt impurity concentration, and shield thickness. Thereafter, the 60 Co contribution to external package dose rate decreased with increasing decay time. Cobalt-60 is primarily produced by neutron reactions with the cobalt impurity in steel and Inconel. Other important radionuclides in activated metals contributing to package external dose rate are radionuclides with relatively short decay times, including 51 Cr, 59 Fe, 58 Co, and 54 Mn. These radionuclides may be represented as an equivalent 60 Co activity/source because 60 Co gamma ray emissions are bounding in terms of source strength and energy to other important radionuclides identified in the analyzed activated metals. Approaches for modeling the neutron-activated corrosion products that may be attached to activated reactor components were analyzed in this report. It was demonstrated that a surface source is more conservative than a uniform volumetric source for the treatment of neutron-activated corrosion products with respect to external package dose rates. An analysis of the maximum radionuclide loadings reported on spent resins identified the radionuclides 137 Cs, 60 Co, 134 Cs, 65 Zn, and 58 Co as the primary contributors to external package dose rate. For a resin cooled for 3.08 years, the external package dose rate was entirely produced by the reported 137 Cs and 60 Co inventory. The neutron sources from actinides found on spent resins or activated metals produced negligible dose rates and may be ignored in dose rate analyses. Effects of idealized waste material, source geometry, and spatial material/source distributions on external package dose rates were determined based on dose rate results for a simplified cask model under normal conditions of transport. Type 304 stainless steel, zirconium, and aluminum with adjusted mass densities based on a maximum content weight were analyzed for material modeling effects on external package dose rate. These materials produced identical external package dose rates within the statistical uncertainties of the dose rate estimates. Among four different source geometry configurations with homogeneous material of different mass densities, uniform volumetric source distribution, and the same source strength, the geometry configuration with lowest mass density (i.e., minimum self-shielding effect among the four cases) was most conservative. Spatial source distributions that better represent localized peak 60 Co activity values were more conservative than a uniform volumetric source distribution, assuming the same weight and total source strength per package. The increase in external dose rate caused by localized source peaks can be as much as the ratio between source peak activity density to average activity density, depending on the location of the activated metal with peak activity density. Therefore, the shielding analysis may be simplified if localized peak activities and the average activity per package can be measured/determined and documented at the time of cask loading. For simplicity, external package dose rates may be determined based on average source activity and a uniform volumetric source distribution. The dose rate results from that calculation model multiplied by the ratio of peak activity density to average activity density will produce maximum dose rate values for conservative estimates.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1999

UV Spectral Variability and the Lyman-Alpha Forest in the Lensed Quasar Q0957+561

Far-ultraviolet spectra of the gravitational lens components Q0957+561 A and B were obtained with the Hubble Space Telescope Faint Object Spectrograph (FOS) at five equally spaced epochs, one every two weeks. We confirm the flux variability of the quasar's Lyman-alpha and O VI lambda 1037 emission lines reported in IUE (International Ultraviolet Explorer) spectra. The fluxes in these lines vary on a time scale of weeks in the local rest frame (LRF), independently of each other and of the surrounding continuum. The individual spectra of each image were co-added to investigate the properties of the Lyman-alpha forest along the two lines of sight to the quasar. Absorption lines having equivalent width W (sub lambda) greater than or equal to 0.3 Angstroms in the LRF not previously identified as interstellar lines, metal lines, or higher order Lyman lines were taken to be LY-alpha forest lines. The existence of each line in this consistently selected set was then verified by its presence in two archival FOS spectra with approx. 1.5 times higher signal to noise than our co-added spectra. Ly-alpha forest lines with W (sub lambda) greater than or equal to 0.3 Angstroms appear at 42 distinct wavelengths in the spectra of the two images. Two absorption lines in the spectrum of image A have no counterpart at that wavelength in the spectrum of image B, and two lines in image B have no counterpart in image A. Based on the separation of the lines of sight at the redshift of the absorption lines appearing in only one spectrum, the density of the absorbing clouds in the direction of Q0957+561 must change significantly over a distance of 79 (+34, -26) h (sub 50) (sup -1) kpc in the simplified model where the absorbers are treated as spherical clouds and the characteristic dimension is the radius. (We adopt H (sub 0)= 50 h (sub 50) km s (sup -1) kpc (sup -1), q (sub 0) = 1/2, and LAMBDA = 0 throughout the paper.) The two limits define the 68% confidence interval on the characteristic dimension, equivalent to the 1 sigma confidence interval for a Gaussian distribution. The 95% confidence interval extends from (32 - 250) h (sub 50) (sup -1) kpc. We show in the Appendix that the fraction of Ly-alpha forest lines that appear in only one spectrum can be expressed as a rapidly converging power series in 1/r, where r the ratio of the radius of the cloud to the separation of the two LOS at the redshift of the cloud. This power series can be rewritten to give r in terms of the fraction of Ly-alpha forest wavelengths that appear in the spectrum of only one image. A simple linear approximation to the solution which everywhere agrees with the power series solution to better than 0.8% for r greater than or equal to 2 is derived in the Appendix.

Dolan, J. F.↗
NASA NTRS2016

Photospheric Current Spikes and Their Possible Association with Flares - Results from an HMI Data Driven Model

A data driven, near photospheric magnetohydrodynamic model predicts spikes in the horizontal current density, and associated resistive heating rate per unit volume Q. The spikes appear as increases by orders of magnitude above background values in neutral line regions (NLRs) of active regions (ARs). The largest spikes typically occur a few hours to a few days prior to M or X flares. The spikes correspond to large vertical derivatives of the horizontal magnetic field. The model takes as input the photospheric magnetic field observed by the Helioseismic & Magnetic Imager (HMI) on the Solar Dynamics Observatory (SDO) satellite. This 2.5 D field is used to determine an analytic expression for a 3 D magnetic field, from which the current density, vector potential, and electric field are computed in every AR pixel for 14 ARs. The field is not assumed to be force-free. The spurious 6, 12, and 24 hour Doppler periods due to SDO orbital motion are filtered out of the time series of the HMI magnetic field for each pixel using a band pass filter. The subset of spikes analyzed at the pixel level are found to occur on HMI and granulation scales of 1 arcsec and 12 minutes. Spikes are found in ARs with and without M or X flares, and outside as well as inside NLRs, but the largest spikes are localized in the NLRs of ARs with M or X flares. The energy to drive the heating associated with the largest current spikes comes from bulk flow kinetic energy, not the electromagnetic field, and the current density is highly non-force free. The results suggest that, in combination with the model, HMI is revealing strong, convection driven, non-force free heating events on granulation scales, and that it is plausible these events are correlated with subsequent M or X flares. More and longer time series need to be analyzed to determine if such a correlation exists. Above an AR dependent threshold value of Q, the number of events N(Q) with heating rates greater than or equal to Q obeys a scale invariant power law distribution for each AR given by N(Q) varies Q(sup -s), where 0.40 less than or equal to S less than or equal to 0.53, with a mean and standard deviation across the 14 ARs of 0.47 and 0.045, showing there is little variation of S from one AR to another. These properties of N(Q) are in close agreement with those of the distribution N(E) for the total energy E of solar flares, determined from observations to be N(E) = constant x E(sup -alpha). From observations of nanoflares in the 0.7 to 4 MK range, and from observations of flares in hard X-rays, it is found that 0.51 less than or equal to alpha less than or equal to 0.57, and 0.4 less than or equal to alpha less than or equal to 0.6, respectively (Crosby et al. 1993, Sol. Phys., 143, 275; Aschwanden & Parnell 2002, ApJ, 572, 1048). Observations also show that, as is found here for the exponent S, there is little variation of alpha with AR (Wheatland 2000, ApJ, 532, 1209), indicating N(E) and N(Q) are largely independent of individual properties of ARs such as area, total magnetic flux, and distribution of current density (i.e. non-potentiality). Therefore the power law scaling of the photospheric heating rate Q computed here on granulation scales is essentially identical to that found for coronal observations of flare energies on scales 1-2 orders of magnitude larger. This suggests the physical mechanisms that cause Q and coronal flares are closely related. It seems likely that Q is the signature of a magnetic reconnection process in an energy range and volume orders of magnitude smaller than those of flares. In this context, at least the larger spikes in Q might be signatures of UV photospheric or lower chromospheric bombs in which plasma is heated to temperatures approximately 10(exp -5) K (Peter et al. 2014, Science 346, 1255726; Judge 2015, ApJ, 808, 116). In addition, lattice based avalanche simulations of flare energy release predict 0.4 less than or equal to alpha less than or equal to 0.5, while analytic, fractal-diffusive self-organized criticality models predict 0.4 less than or equal to alpha less than or equal to 0.67, in excellent agreement with observations, and the results presented here (Aschwanden & Parnell 2002, ApJ, 572, 1048; Aschwanden 2012, A&A, 539, A2; Aschwanden 2013, in "Self Organized Criticality Systems"; Aschwanden et al. 2016, SSR, 198, 47).

Goodman, Michael↗
NASA NTRS2005

Summary of Results from Analyses of Deposits of the Deep-Ocean Impact of the Eltanin Asteroid

Deposits of the late Pliocene (2.5 Ma) Eltanin impact are unique in the known geological record. The only known example of a km-sized asteroid to impact a deep-ocean (5 km) basin, is the most meterorite-rich locality known. This was discovered as an Ir anomaly in sediments from three cores collected in 1965 by the USNS Eltanin. These cores contained mm-sized shock-melted asteroid materials and unmelted meteorite fragments. Mineral chemistry of meteorite fragments, and siderophole concentrations in melt rocks, indicate that the parent asteroid was a low-metal (4\%) mesosiderite. A geological exploration of the impact in 1995 by Polarstern expedition ANT-XIV4 was near the Freeden Seamounts (57.3S, 90.5 W), and successfully collected three cores with impact deposits. Analyses showed that sediments as old as Eocene were eroded by the impact disturbance and redeposited in three distinct units. The lowermost is a chaotic assemblage of sediment fragments up to 50 cm in size. Above this is a laminated sand-rich unit that deposited as a turbulent flow, and this is overlain by a more fine-grained deposit of silts and clays that settled from a cloud of sediment suspended in the water column. Meteoritic ejecta particles were concentrated near the base of the uppermost unit, where coarse ejecta caught up with the disturbed sediment. Here we will present results from a new suite of cores collected on Polarstern expedition ANT-XVIIU5a. In 2001, the Polarstern returned to the impact area and explored a region of 80,000 sq-km., collecting at least 16 sediment cores with meteoritic ejecta. The known strewn field extends over a region 660 by 200 km. The meteoritic ejecta is most concentrated in cores on the Freeden seamounts, and in the basins to the north, where the amount of meteoritic material deposited on the ocean floor was as much as 3 g/sq-cm. These concentrations drop off to the north and the east to levels as low as approximately 0.1 g/sq-cm. We were unable to sample the impact south and west of the seamounts, as the deposit was buried beyond the reach of our 25 m piston corer. We estimate that ground zero was in the region just north, or northwest, of the seamounts. There is no evidence that the impactor penetrated the ocean floor or formed a crater. The composition of the melted ejecta is inconsistent with mixing between projectile and terrestrial materials other than seawater salts. X-ray radiographs of sediments reveal details not seen in earlier cores. The uppermost impact unit is well-preserved in several cores, found as much as 50 km from the seamounts to the east, north, and west of the seamounts, where at least 25 cm of this unit is preserved. At greater distances burrowing organisms have mixed the sediments so if this unit did exist, it was too thin to survive bioturbation. These finegrained sediments are clearly laminated, and show alternating layers of low- and high-density (meteoritic) sediments, consistent with ripple formation in an energetic flow regime. We have extracted 35 g of meteoritic melt rock and 3 g of meteorite fragments from sieved sediments. Additionally a 9 g, 2.2 cm meteorite was recovered during opening of one core. The fact that 9\% of the coarse ejecta is unmelted meteorites may be characteristic of deep-ocean impacts. This may have significance for delivery of organic matter to the early Earth by small impacts into primordial oceans, where actual meteorite fragments can survive in significant amounts. However, a large portion of the meteoritic debris is buried rapidly by the sediments disturbed by the impact.

Kyte, Frank T.↗
DOE OSTI2025

Modeling and Simulation of Fuel Dispersal During the Loss-of-Coolant Accident

This document is the compilation of the milestone portion to a larger end of project NEUP report. The executive summary of the modeling portion is provided below: In the event of cladding rupture during a postulated LOCA in a pressurized water reactor, fuel particles, along with fission gases, can be expelled into the reactor core from the fractured fuel rod, a phenomenon referred to as fuel dispersal. The initial stage of fuel dispersal is strongly influenced by the high-pressure ejection of fuel fragments, the size and geometry of the ruptured cladding, and the depressurization history of the fuel rod during the postulated LOCA transient. Depending on the location of the burst orifice relative to the quench front, the dispersal event represents an intricate three-phase flow and heat transfer phenomenon, where high-temperature fuel particles carried by the fission gases interact with the coolant within the narrow subchannels of the fuel assemblies, inducing localized phase change. Given the unique multiphysics nature of this phenomena, the current study develops a dedicated computational framework to predict the mass distribution and cooling of dispersing fuel particles, facilitating post-accident assessment and management of the fuel assemblies. Considering the scale of nuclear reactor applications, a continuum three-fluid model is proposed for simulating the transport of solids within the reactor core. With high-temperature fuel fragments within the liquid media, nucleation sites inducing phase changes are dispersed within the flow domain. Coupled with the fact that the transient dispersal event occurs on different time scales than other three-phase flow applications, this study derives a time-averaged three-fluid flow model without losing generality. The assumptions regarding the continuum treatment of the solid phase and the modeling of fuel dispersal behavior are incorporated to simplify the governing equations and derive applicable closure relations. The computational validation of the model was conducted using adiabatic experimental results obtained from ongoing research at Oregon State University, focusing on characterizing fuel dispersal behavior during simulated LOCA conditions. Settlement characteristics of the solids, quantified by the probability distribution of equivalent particles, closely matched the probability density functions reported in experimental studies. The transport of fuel particles within a scaled 5 × 5 lattice of a pressurized-water reactor rod bundle geometry was modeled through a two-fluid Eulerian framework. The required boundary conditions were evaluated from the fuel performance code BISON in a postulated large-break LOCA scenario. The modeling framework considered solid fuel particles as granular matter, interacting with the gaseous dry steam phase and fission gases through the governing interfacial momentum exchange between the participating fluids. The simulation results provided the volume fraction of the solids obtained at the bottom surface of the enclosing tank geometry. Postulated LOCA leading to fuel dispersal phenomena involves the strong coupling between fuel thermomechanics, cladding deformation, thermal-hydraulics, and fuel particle transport. Incorporation of such a strong coupling in numerical simulation is performed by coupling the multiphysics solvers. In the case of fuel dispersal, a strong coupled simulation can be performed by coupling the BISON code for fuel performance, the TRACE code for system-level thermal hydraulics, and fuel particle transport in Multiphysics Object-Oriented Simulation Environment (MOOSE). For such intricate infrastructure, the MOOSE Framework eases the data transfer between codes. The recent version of MOOSE has incorporated the Navier-Stokes module for the fluid flow. An exploratory exercise was done to gain familiarity with finite volume capabilities in the MOOSE framework to incorporate the Spalart-Allmaras (SA) turbulence model. New finite-volume and auxiliary kernels were introduced to assemble the SA transport equation, compute turbulent viscosity, and evaluate wall distance and diagnostic turbulence terms, fully integrated with existing Navier-Stokes modules. A turbulent lid-driven cavity at a Reynolds number of approximately 10,000 is used for verification. MOOSE shows the robust solver convergence and produces the turbulent features. But it underpredicts the velocity profile and turbulent quantities, emphasizing the need to develop improved SA near-wall treatments (e.g., low-Re corrections or wall functions) as a key direction for future work.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1998

Spatially Resolved Spectroscopy of the SNR IC443

IC 443 is a supernova remnant of intermediate age, i.e. a few thousand years. It is especially interesting because part of its periphery is expanding into a molecular cloud while other sections are expanding into a typical interstellar medium of much lower density. Since the evolution of a supernova remnant through its various phases is affected by the density of the medium it expands into with the reasonable assumption that the supernova explosion was approximately symmetric we have an opportunity to observe a single object in two phases simultaneously. It was observed by ASCA in April, 1993 for a short period during the PV phase and more thoroughly in a 42 ksec exposure in March, 1994. The latter measurement provides most of the results that have been reported. Most of the analysis took place after the grant ended but is included here for completeness. The data was sent simultaneously to US and Japanese Pls. We worked independently. The software set of FTOOLs was used to construct images and spectra. They were judged to be rather unintuitive and not at all user friendly. I found I was using one FTOOL to read the header to obtain information that would only be provided to another FTOOL. The Japanese investigators were more successful. They analyzed the data and published results more rapidly. The scientific results summarized below are based primarily on their publications. Since IC 443 is an interesting example of a middle aged SNR in which a variety of processes are occurring it is one of a class. IC 443 exhibits shell-like emission in hard X-rays and extended soft X-rays with thin thermal spectra. It resembles SN 1006 in these respects. IC 443 contains hard X-rays in a semi-circular shell surrounding the thermal component. The total hard X-ray flux in the ASCA FOV is only a half of the Ginga hard component; which suggests that the hard X-rays are not confined only in the shell but some are extended larger than the ASCA FOV of eq 1 degree diameter. Japanese investigators examined the spatial structure of the thermal component and analyzed the GIS spectra with a non-equilibrium plasma model, and found no systematic variation of the interstellar absorption across the remnant. Evidence for shock acceleration of cosmic rays to high energies (10 TeV) was found by Keohane. X-ray imaging spectroscopy with ASCA reveals two regions of particularly hard emission: an unresolved source embedded in an extended emission region, and a ridge of emission coincident with the southeastern rim. Both features are located on part of the radio shell where the shock wave is interacting with molecular gas, and together they account for a majority of the emission at 7 keV. Though we would not have noticed it a priori, the unresolved feature is coincident with one resolved by the ROSAT HRI. The ASCA measurements were combined with higher energy data from the XTE and GRO missions and with radio and TeV gamma-ray data to produce a nonthermal multiwavelength spectrum for IC 443 which was fit with a cosmic ray interaction model. This model calculates the cynchrotron, bremsstrahlung, invers Compton, and neutral pion decay emission produced by locally accelerated cosmic ray interacting with ambient matter, soft photon fields, and magnetic fields.

Gorenstein, P.↗