Search NASA⌕ Search

SEARCH · Search NASA

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.

468 records · Page 26
All sourcesNASA NTRSDOE OSTI
NASA NTRS2023

Effects of Heating Configuration and Operating Parameters on Heat Transfer and Interfacial Physics of Microgravity Flow Boiling With Subcooled Inlet Conditions –Experiments Onboard the International Space Station

This study is part of the Flow Boiling and Condensation Experiment (FBCE), a collaborative effort between the Purdue University Boiling and Two-Phase Flow Laboratory (PU-BTPFL) and the NASA Glenn Research Center. The FBCE fitted with the Flow Boiling Module (FBM) was launched to the International Space Station (ISS) in August 2021 and experiments were successfully performed from February to July 2022 to amass a large microgravity-flow-boiling database. This study is focused on heat transfer and flow visualization of microgravity flow boiling of n-Perfluorohexane in a rectangular channel of 5.0 mm height, 2.5 mm width (heated), and 114.6 mm length, with subcooled inlet conditions. High-speed-video photography is utilized to present flow patterns and temporal interfacial behavior. Heat transfer results are presented in the form of flow boiling curves and both parametric curves and streamwise profiles of wall temperature and heat transfer coefficient. Firstly, the parametric effects of mass velocity (199.4 – 3200.0 kg/m 2 s), inlet subcooling (0.2 – 46.0°C), and inlet pressure (124.2 – 176.7 kPa), on the aforementioned aspects are assessed for double-sided heating to establish them for a microgravity environment. Of these three parameters, mass velocity and inlet subcooling mostly determine the microgravity flow boiling behavior, while inlet pressure plays an insignificant role. Flow patterns for double-sided heating are more complex than those for single-sided heating due to interaction between the two vapor layers. Vapor interaction is minimized at high subcoolings and high mass velocities due to strong condensation offered by the subcooled bulk liquid layer separating them. Despite the different flow patterns, both single- and double-sided heating generally result in similar parametric trends and local heat transfer coefficients for similar operating conditions. Flow instabilities manifest as temporal flow anomalies and temperature oscillations, and their severity increases with increasing boiling number. Secondly, the effects of heating configuration are analyzed by comparing and contrasting several aspects of single- and double-sided heating data. The heat fluxes at which onset of nucleate boiling degradation (ONBD) and critical heat flux (CHF) occur are distinctly different for single- and double-sided heating. There exists a threshold inlet subcooling demarcating the dominance of flow acceleration and condensation effects in vapor removal from the near-wall region and replenishment of fresh liquid for boiling. Above the threshold, condensation from the near-wall region is dominant and single-sided heating yields higher heat fluxes, and below it, acceleration is dominant and double-sided yields higher heat fluxes. At mass velocity in the range of 200 – 2400 kg/m 2 s, the threshold inlet subcooling lies in the approximate range of 20 – 30°C (corresponding inlet quality of roughly -0.40 – -0.20).

Microgravity↗
NASA NTRS2000

Hierarchical Statistical 3D ' Atomistic' Simulation of Decanano MOSFETs: Drift-Diffusion, Hydrodynamic and Quantum Mechanical Approaches

When MOSFETs are scaled to deep submicron dimensions the discreteness and randomness of the dopant charges in the channel region introduces significant fluctuations in the device characteristics. This effect, predicted 20 year ago, has been confirmed experimentally and in simulation studies. The impact of the fluctuations on the functionality, yield, and reliability of the corresponding systems shifts the paradigm of the numerical device simulation. It becomes insufficient to simulate only one device representing one macroscopical design in a continuous charge approximation. An ensemble of macroscopically identical but microscopically different devices has to be characterized by simulation of statistically significant samples. The aims of the numerical simulations shift from predicting the characteristics of a single device with continuous doping towards estimating the mean values and the standard deviations of basic design parameters such as threshold voltage, subthreshold slope, transconductance, drive current, etc. for the whole ensemble of 'atomistically' different devices in the system. It has to be pointed out that even the mean values obtained from 'atomistic' simulations are not identical to the values obtained from continuous doping simulations. In this paper we present a hierarchical approach to the 'atomistic' simulation of aggressively scaled decanano MOSFETs. A full scale 3D drift-diffusion'atomostic' simulation approach is first described and used for verification of the more economical, but also more restricted, options. To reduce the processor time and memory requirements at high drain voltage we have developed a self-consistent option based on a thin slab solution of the current continuity equation only in the channel region. This is coupled to the Poisson's equation solution in the whole simulation domain in the Gummel iteration cycles. The accuracy of this approach is investigated in comparison with the full self-consistent solution. At low drain voltage only single solution of the nonlinear Poisson equation is sufficient to extract the current with satisfactory accuracy. A pilot version of a hydrodynamic 'atomistic' simulator has been developed in order to study the effect of the nonequilibrium, non local transport in decanano MOSFETs on the random dopant induced current fluctuations. For the first time we have also applied the density gradient approach in 3D to investigate the effect of the quantum confinement on the threshold voltage fluctuations. The developed 'atomistic' simulation techniques have been applied to study various fluctuation resistant MOSFET architectures including epitaxial and delta doped devices.

Asenov, Asen↗
NASA NTRS1994

Steps toward determination of the size and structure of the broad-line region in active galactic nuclei. 5: Variability of the ultraviolet continuum and emission lines of NGC 3783

We report on the results of intensive ultraviolet spectral monitoring of the Seyfert 1 galaxy NGC 3783. The nucleus of NGC 3783 was observed with the International Ultraviolet Explorer satellite on a regular basis for a total of 7 months, once every 4 days for the first 172 days and once every other day for the final 50 days. Significant variability was observed in both continuum and emission-line fluxes. The light curves for the continuum fluxes exhibited two well-defined local minima or 'dips,' the first lasting is less than or approximately 20 days and the second is less than or approximately 4 days, with additional episodes of relatively rapid flickering of approximately the same amplitude. As in the case of NGC 5548 (the only other Seyfert galaxy that has been the subject of such an intensive, sustained monitoring effort), the largest continuum variations were seen at the shortest wavelengths, so that the continuum became 'harder' when brighter. The variations in the continuum occurred simultaneously at all wavelengths (delta(t) is less than 2 days). Generally, the amplitude of variability of the emission lines was lower than (or comparable to) that of the continuum. Apart from Mg II (which varied little) and N V (which is relatively weak and badly blended with Ly(alpha), the light curves of the emission lines are very similar to the continuum light curves, in each case with a small systematic delay or 'lag.' As for NGC 5548, the highest ionization lines seem to respond with shorter lags than the lower ionization lines. The lags found for NGC 3783 are considerably shorter than those obtained for NGC 5548, with values of (formally) approximately 0 days for He II + O III), and approximately 4 days for Ly(alpha) and C IV. The data further suggest lags of approximately 4 days for Si IV + O IV) and 8-30 days for Si III + C III). Mg II lagged the 1460 A continuum by approximately 9 days, although this result depends on the method of measuring the line flux and may in fact be due to variability of the underlying Fe II lines. Correlation analysis further shows that the power density spectrum contains substantial unresolved power over timescales of is less than or approximately 2 days, and that the character of the continuum variability may change with time.

Reichert, G. A.↗
NASA NTRS2002

Simulations of Ground and Space-Based Oxygen Atom Experiments

Fast, pulsed atomic oxygen sources are a key tool in ground-based investigations of spacecraft contamination and surface erosion effects. These technically challenging ground-based studies provide a before and after picture of materials under low-earth-orbit (LEO) conditions. It would be of great interest to track in real time the pulsed flux from the source to the surface sample target and beyond in order to characterize the population of atoms and molecules that actually impact the surface and those that make it downstream to any coincident detectors. We have performed simulations in order to provide such detailed descriptions of these ground-based measurements and to provide an assessment of their correspondence to the actual LEO environment. Where possible we also make comparisons to measured fluxes and erosion yields. To perform the calculations we use a detailed description of a measurement beam and surface geometry based on the W, pulsed apparatus at Montana State University. In this system, a short pulse (on the order of 10 microseconds) of an O/O2 beam impacts a flat sample about 40 cm downstream and slightly displaced &om the beam s central axis. Past this target, at the end of the beam axis is a quadrupole mass spectrometer that measures the relative in situ flux of 0102 to give an overall normalized erosion yield. In our simulations we use the Direct Simulation Monte Carlo (DSMC) method, and track individual atoms within the atomic oxygen pulse. DSMC techniques are typically used to model rarefied (few collision) gas-flows which occur at altitudes above approximately 110 kilometers. These techniques are well suited for the conditions here, and multi-collision effects that can only be treated by this or a similar technique are included. This simulation includes collisions with the surface and among gas atoms that have scattered from the surface. The simulation also includes descriptions of the velocity spread and spatial profiles of the O/O2 beam obtained from separate measurements. These computations use basic engineering models for the gas-gas and gas-surface scattering and focus on the influence of multi-collision effects. These simulations characterize many important quantities of interest including the actual flux of atoms that reach the surface, the energy distribution of this flux, as well as the direction of the velocity of the flux that strikes the surface. These quantities are important in characterizing the conditions which give rise to measured surface erosion. The calculations also yield time- snapshots of the pulse as it impacts and flows around the surface. These snapshots reveal the local environment of gas near the surface for the duration of the pulse. We are also able to compute the flux of molecules that travel downstream and reach the spectrometer, and we characterize their velocity distribution. The number of atoms that reach the spectrometer can in fact be influenced by the presence of the surface due to gas-gas collisions from atoms scattered h m the surface, and it will generally be less than that with the surface absent. This amounts to an overall normalization factor in computing erosion yields. We discuss these quantities and their relationship to the gas-surf$ce interaction parameters. We have also performed similar calculations corresponding to conditions (number densities, temperatures, and velocities) of low-earth orbit. The steady-state nature and lower overall flux of the actual space environment give rise to differences in the nature of the gas-impacts on the surface from those of the ground-based measurements using a pulsed source.

Minton, T. K.↗
NASA NTRS2013

NIR Imaging Spectroscopy of the Inner Few Arcseconds of NGC 4151 with OSIRIS at Keck

We present H- and K-band data from the inner arcsecond of the Seyfert 1.5 galaxy NGC 4151 obtained with the adaptive optics assisted near-infrared imaging field spectrograph OSIRIS at the Keck Observatory. The angular resolution is about a few parsecs on-site and thus competes easily with optical images taken previously with the Hubble Space Telescope. We present the morphology and dynamics of most species detected but focus on the morphology and dynamics of the narrow line region (as traced by emission of [FeII]λ1.644 μm), the interplay between plasma ejected from the nucleus (as traced by 21 cm continuum radio data) and hot H2 gas and characterize the detected nuclear HeIλ2.058 μm absorption feature as a narrow absorption line (NAL) phenomenon. Emission from the narrow line region (NLR) as traced by [FeII] reveals a biconical morphology and we compare the measured dynamics in the [FeII] emission line with models proposing acceleration of gas in the NLR and simple ejection of gas into the NLR. In the inner 2.5 arcseconds the acceleration model reveals a better fit to our data than the ejection model.We also see evidence that the jet very locally enhances emission in [FeII] at certain positions in our field-of-view such that we were able to distinct the kinematics of these clouds from clouds generally accelerated in the NLR. Further, the radio jet is aligned with the bicone surface rather than the bicone axis such that we assume that the jet is not the dominant mechanism responsible for driving the kinematics of clouds in the NLR. The hot H2 gas is thermal with a temperature of about 1700 K. We observe a remarkable correlation between individual H2 clouds at systemic velocity with the 21 cm continuum radio jet. We propose that the radio jet is at least partially embedded in the galactic disk of NGC 4151 such that deviations from a linear radio structure are invoked by interactions of jet plasma with H2 clouds that are moving into the path of the jet because of rotation of the galactic disk of NGC 4151. Additionally, we observe a correlation of the jet as traced by the radio data, with gas as traced in Brγ and H2, at velocities between systemic and +/- 200 km/s at several locations along the path of the jet. The HeIλ2.058 μm line in NGC 4151 appears in emission with a blueshifted absorption component from an outflow. The emission (absorption) component has a velocity offset of 10 km/s (-280 km/s) with a Gaussian (Lorentzian) full-width (half-width) at half maximum of 160 km/s (440 km/s). The absorption component remains spatially unresolved and its kinematic measures differ from that of UV resonance absorption lines. From the amount of absorption we derive a lower limit of the HeI 2S column density of 1 × 10(exp 14) cm−2 with a covering factor along the line-of-sight of C(sub los) approximately equal to 0.1.

K-band data↗
NASA NTRSDate not supplied

Effect of Vacuum on Force Response of an Ultrasonic Penetrator

Introduction: The Apollo astronauts encountered higher than expected resistances when interacting with the lunar soil via the Apollo Lunar Surface Drill (ALSD) and the trenching tool. Reducing the force required to move tools or other mechanical components through regolith will impact many steps of the resource extraction process. Force reduction has been achieved in soil materials by imparting vibration to tooling interfaces such as a vibratory farming cultivator, a percussive scoop, and ultrasonically resonant penetrators. Vibration-assisted tools in granular media reduce interaction forces by fluidizing a volume around the tool, allowing the tool to progress through a dynamic (fluid) medium instead of a static (solid) medium. This work seeks to quantify ultrasonic vibration’s effect on the force response of a penetrator in lunar soil simulant in vacuum sufficient to be within the molecular flow regime of any disturbed gases. Methods: A custom vacuum chamber setup, CUBEvac, was designed and built to facilitate penetration testing in a high vacuum environment, for comparison to penetration behavior in ambient terrestrial environment. A two-stage pumping system (Agilent Triscroll 600 roughing pump, Agilent VHS-6 oil diffusion pump) reached chamber pressures of about 5x10-6 Torr with regolith simulant in place. Figure 1 is a schematic of the heart of the assembly (note the penetration drive mechanisms above the chamber feedthrough and the regolith simulant sample in the bottom are not shown). The penetration actuation stack was comprised of a stepper motor driving a lead screw to move the ultra-sonic probe vertically inside the chamber. Motion was coordinated with an Arduino Uno. GRC-3 lunar simulant was used for this set of experiments. Samples were prepared in a four-liter stain-less steel, cylindrical pot with an internal diameter of approximately 15.56 cm (6.125 in) and a depth of 19.37 cm (7 5/8 in) for testing. The soil was baked out prior to compaction preparation as a measure to reduce soil moisture which interfered with pump down capacity. The soil was not baked again if it was removed from the vacuum chamber, prepped, and immediately returned to the vacuum chamber for pump down. The soil was compacted using a 60 Hz vibration table with a surcharge of 34 kg place on top of the soil in the container. Prepared soil samples weighed approximately 6.5 kg (bulk density 1.895-1.934 g/cm3). Two probe end effectors were tested: A cone penetrometer (static only) with a nominal diameter of 12.7 mm (0.5 in) and a nominal height of 28.6 mm (1.125 in); and a vibrating cylindrical probe measuring 12.7 mm in diameter and 50.8 mm in effective length from the tip (Figure 2). The cylindrical probe vibrated resonantly at 20 kHz with an amplitude 23 μm. The cone penetration tests were conducted to assess potential soil behavior differences in vacuum. The cylinder probe tests were conducted as the primary subject of this investigation to assess force response in vacuum. For each test, a regolith simulant sample was load-ed and compacted in the chamber, which was then evacuated for roughly 18 hours to reach the lowest possible pressure (approximately 5x10-6 Torr for most tests). The probe was then moved to about 10 mm above the soil surface before being pushed to a depth of 50 mm for the cylinder probe tests and to a depth of 100 mm for the cone penetrometer tests, both at a speed of 2 mm/s. The simulant samples were prepared the same for all tests. Ideally, they would respond consistently to probe penetration under ambient and vacuum conditions. This was evaluated by measuring the resistance of representative prepared simulant beds with a standard cone penetrometer in both environments. Results and Discussion: The resistance of the simulant samples in the vacuum tests was consistently lower than in the ambient tests as determined by the cone penetration tests. Thus, the ambient and vacuum results cannot be compared directly; work is underway to de-confound and better correlate the data. Still, figures 3-6 show that probe penetration forces are lower overall in the vacuum environment. In both environments, resonant vibration of the probe provides two useful effects: It reduces the probe penetration force and smooths the force-depth curve, significantly reducing local maxima. These effects have implications for various potential applications, such as astronaut hand-tools, where benefits (reducing astronaut effort) outweigh the cost of the additional energy re-quired to generate vibration. These results demonstrate that resonantly vibrating tools can meaningfully reduce the penetration force required for excavation, probing, and drilling tools in simulated lunar regolith deposits under vacuum levels approaching those that will be experience on the Moon’s surface. Lunar-gravity, ambient environment tests are scheduled soon. The effects of realistic temperatures and temperature gradients and deeper vacuum remain to be tested.

E Rezich↗
NASA NTRSDate not supplied

Large Primtive Asteroids: Thermal and Dynamical Context

Primitiveasteroids, most of which are located in the outer belt and Jupiter’s Trojan clouds, provide information related to the origin and evolution of the solar system and the conditions in which the solar nebula was formed. These asteroids are widely thought to be the origin of the least-altered car-bonaceous chondrite that allow us to put crucial con-straints on the current dynamical and thermal theories of the formation and evolution of the early solar sys-tem. The nature of surface composition of large and low-albedo asteroids, like (1) Ceres, (10) Hygiea and (52) Europa, is still under intense debate and different interpretations have beenput forth to explain the ab-sorption features in these objects (e.g.,[1, 2, 3,4,5,6]). Laboratory and spectroscopic experiments on me-teorites that represent all nine carbonaceous chondrite types also found no spectral matches for these large asteroids (e.g.,[7]). Previous studies of asteroid Ceres (the largest primitive asteroid in the solar system) have been conducted to constrain and estimate its surface composition (e.g., [2, 5]). Using linear mixing,[2]found hydroxide brucite, serpentines, and carbonates, to be consistent with Ceres’ ground-based spectra.[5]estimated the surface composition of Ceres and found evidence of widespread NH3-phyllosilicates across its surface using best-fit solutions to Dawn’s NIR spectra. The presence of NH3-phyllosilicates implies that mate-rial from the outer solar system was incorporated into large primitive asteroids, either during their formation at great heliocentric distance or by incorporation of material transported into the Main Belt region. Here we presentnew largeand primitive asteroidsthat share the same spectral similaritieswith the largest asteroid in the solar system, Ceres.We also present the context of these new observations in terms oftheirthermal and dynamicalevolution.Large Primitive Asteroids: Dynamical and Thermal Context:Orbital Distribution of LargePrimitiveAsteroids: Constraining the mineralogy and surface composition of large primitive asteroids will place crucial con-straints on existing dynamical and thermal theories of the formation of the early solar system. We identified several additional asteroids in the Europa-like group in addition to asteroids Europa, Euphrosyne, and Patien-tia, that wereidentified in [8]using NASA IRTFtele-scope (e.g.,Figure 1). These new asteroids along with the already-observed Ceres-and Europa-like group members are localized in the 2.8 < a < 3.4 AU region and characterized by larger sizes, showing an interest-ing orbital distribution (Figure 2). Ceres-likeand Eu-ropa-like groups, which include the largest asteroids in the solar system, show an interesting orbital distribu-tion. These groups are located in the 2.6 < a < 3.6 AU region that contains the snow-line. The snow-line’s location may have been driftedinwards due to the disk’s cooling and evolution[9,10]. Recent dynamical models[11,12]suggested that a substantial fraction of primitive asteroids originated between or beyond the giant planets (a > 5 AU), where water ice would have been stable, and then implanted in the outer Main Belt region because of the giant planets’ growth.Figure 1. Two asteroids, 94 Aurora and 423 Diotima, showing spectra similar to asteroids (1)Ceres and (52)Europa.Thermal modeling and evolution of primitive aster-oids:Primitive water-rich asteroids are thought to be originally composed of mixtures of anhydrous materi-als and water ice that waslater melted by heating sources such as the decay of 26Al, reacting with anhy-drous materials to form H2O/OH-rich minerals. Calcu-lations of the evolution of the temperature and struc-ture of icy planetesimals were performedusing a 1D finite differences thermal evolution model[13, 14] for 26Al-heated planetesimals. In particular, thermally activated compaction due to hot pressing of bodies with an initially unconsolidated porous structure is included. An ice-rich initial composition that leads to a material dominated by phyllosilicates upon aqueous alteration (with 25 vol% H2O and a rock fraction that contains 85 vol% phyllosilicates and 15 vol% olivine upon aque-ous alteration, similar to CIand CM chondrites) was assumed. Atypical initial porosity of 40%[15] is re-duced following the change of the strain rate that is calculated as Voigt approximation from the strain rates of components[16]. Material properties (thermal con-ductivity, density, heat capacity, etc.) correspond to the composition assumed and are adjusted with tempera-ture and porosity. Melting of the water ice as well as water-rock separation are included[14]. Both short-and long-lived radionuclides are consideredas heat sources. Figure 3 shows the maximum temperature calculated as a function of radius and accretion time. A variety of internal structures is obtained, ranging from primordial (no melting of water ice) over partially melted or partially differentiated (melting of water ice, hydration, formation of a rocky core and water ocean below an undifferentiated layer) to completely differ-entiated ones (rocky core, water mantle, Enceladus-like case). The heating and differentiation of planetesi-mals is determined by the availability of 26Al, i.e., by the accretion time t0 relative to the formation of the calcium-aluminum-rich inclusions (CAIs), such that maximum temperatures and structures vary stronglyfor t0< 6 Myr rel. to CAIs. However, for a later accre-tion only the size of the body determines its maximum tempera-ture and structure due to the nearly constant heating by long-lived radionuclides. Average densities of Ceres-and Europa-like group membersimply high-ly porous interiors and, consequently, relatively late accretion at t0> 3 Myr rel. to CAIs with a maximum temperature of < 600 K (Figure3)

D Takir↗
NASA NTRS2001

Atmospheric Constraints on Landing Site Selection

The Martian atmosphere is a significant part of the environment that the Mars Exploration Rovers (MER) will encounter. As such, it imposes important constraints on where the rovers can and cannot land. Unfortunately, as there are no meteorological instruments on the rovers, there is little atmospheric science that can be accomplished, and no scientific preference for landing sites. The atmosphere constrains landing site selection in two main areas, the entry descent and landing (EDL) process and the survivability of the rovers on the surface. EDL is influenced by the density profile and boundary layer winds (up to altitudes of 5 to 10 km). Surface survivability involves atmospheric dust, temperatures and winds. During EDL, the atmosphere is used to slow the lander down, both ballistically and on the parachute. This limits the maximum elevation of the landing site to -1.3 km below the MOLA reference aeroid. The landers need to encounter a sufficiently dense atmosphere to be able to stop, and the deeper the landing site, the more column integrated atmosphere the lander can pass through before reaching the surface. The current limit was determined both by a desire to be able to reach the hematite region and by a set of atmosphere models we developed for EDL simulations. These are based on Thermal Emission Spectrometer (TES) atmospheric profile measurements, Ames Mars General Circulation Model (MGCM) results, and the 1-D Ames GCM radiative/convective model by J. Murphy. The latter is used for the near surface diurnal cycle. The current version of our model encompasses representative latitude bands, but we intend to make specific models for the final candidate landing sites to insure that they fall within the general envelope. The second constraint imposed on potential landing sites through the EDL process is the near surface wind. The wind in the lower approximately 5 km determines the horizontal velocity that the landers have when they land. Due to the mechanics of the landing process, the total velocity (including both the horizontal and vertical components) determines whether or not the landers are successful. Unfortunately, the landing system has no easy way to nullify any horizontal velocity imparted by the wind, so the landing sites selected need to have as little wind as possible. In addition to the mean wind velocity, the landing system is sensitive to vertical wind shear in the lowest kilometer or so. Wind shear can deflect the retro rockets (RADs) from their nominal vertical orientation producing unwanted horizontal spacecraft velocities. Both mean velocity and wind shear are dominated by the the local topography and other surface properties (in particular albedo and thermal inertia which control the surface temperature). This is seen even in simplified 2-D mesoscale models. The effects in a fully 3-D model are expected to he even more topographically dependent. In particular there is potential for wind channeling in canyons and other terrain features. Boundary layer winds and wind shear are currently being modeled based on terrestrial data and boundary layer scaling laws modified for Martian conditions. We hope to supplement this with mesoscale model results (from several sources) once the number of landing sites is reduced to a manageable number.

Kass, David M.↗
DOE OSTI2024

A static quantum embedding scheme based on coupled cluster theory

Here, we develop a static quantum embedding scheme that utilizes different levels of approximations to coupled cluster (CC) theory for an active fragment region and its environment. To reduce the computational cost, we solve the local fragment problem using a high-level CC method and address the environment problem with a lower-level Møller–Plesset (MP) perturbative method. This embedding approach inherits many conceptual developments from the hybrid second-order Møller–Plesset (MP2) and CC works by Nooijen [J. Chem. Phys. 111, 10815 (1999)] and Bochevarov and Sherrill [J. Chem. Phys. 122, 234110 (2005)]. We go beyond those works here by primarily targeting a specific localized fragment of a molecule and also introducing an alternative mechanism to relax the environment within this framework. We will call this approach MP-CC. We demonstrate the effectiveness of MP-CC on several potential energy curves and a set of thermochemical reaction energies, using CC with singles and doubles as the fragment solver, and MP2-like treatments of the environment. The results are substantially improved by the inclusion of orbital relaxation in the environment. Using localized bonds as the active fragment, we also report results for N=N bond breaking in azomethane and for the central C–C bond torsion in butadiene. We find that when the fragment Hilbert space size remains fixed (e.g., when determined by an intrinsic atomic orbital approach), the method achieves comparable accuracy with both a small and a large basis set. Additionally, our results indicate that increasing the fragment Hilbert space size systematically enhances the accuracy of observables, approaching the precision of the full CC solver.

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

First-Generation Jet Propulsion Laboratory "Hockey-Puck" Free-Flying Magnetometers for Distributed In-Situ Multiprobe Measurement of Current Density Filamentation in the Northern Auroral Zone: Enstrophy Mission

The sub-orbital rocket mission was a collaborative project between the University of New Hampshire, Cornell University, and the Jet Propulsion Laboratory (JPL) to study filamentation phenomena in the northern Auroral zone. The Enstrophy mission test flies the JPL Free-Flying Magnetometer (FFM) concept. The FFM technology development task has been funded by NASA develop miniaturized, low-power, integrated "sensorcrafts". JPL's role was to design, integrate, test, and deliver four FFMs for deployment from the sounding rocket, allowing a unique determination of curl-B. This provides a direct measurement of magnetic-field-aligned current density along the rocket trajectory. A miniaturized three-axis fluxgate magnetometer was integrated with a 4-channel 22-bit sigma-delta Analog to Digital Converter (ADC), four temperature sensors, digital control electronics, seven (Li-SOCl2) batteries, two (4 deg x 170 deg field of view) sun-sensors, a fan-shaped-beam laser diode beacon, a (16 MHz) stable Temperature Compensated Crystal Oscillator (TCXO) clock, Radio Frequency (RF) communication subsystem, and an antenna for approximately 15 minutes of operation where data was collected continuously and transmitted in three (3) bursts (approximately 26 seconds each) to ground station antennas at Poker Flat, Alaska. FFMs were stowed within two trays onboard the rocket during the rocket launch and were released simultaneously using the spinning action of the rocket at approximately 300 km altitude (approximately 100 sec. into the flight). FFMs were deployed with spin rate of approximately 17 Hz and approximately 3 m/sec linear velocity with respect to the rocket. For testing purposes while the rocket was in the launch pad and during flight prior to release of FFMs from the rocket, commands (such as "power on", "test", "flight", "power off', and clock "Reset" signal) were transmitted via a infrared Light Emitting Diode to an infrared detector in the FFM. Special attention was paid to low magnetic signature electronic design and choice of materials in packaging. The miniaturized fluxgate magnetometers had a range of 1-60000 nT with 0.1% full-scale linearity. The frequency range of interest for magnetic measurement was 10 mHz - 50 Hz. Digital data from the magnetometer's three axes were placed in a 4MB Static Random Access Memory (SRAM) in data packages (frames) formatted together with time tags and frame ID. After a specified time was elapsed, the data were Viterbi encoded and transmitted at a rate of 100 kbps (BPSK). Each of the four FFMs transmitted at different frequency. These carrier frequencies were in the range of 2200-2300 MHz. The antenna was a single patch on a high dielectric constant substrate covering one end-plate of the hockey-puck-sized unit. The local clocks aboard the FFMs were reset at the start of the mission and stayed synchronized within 3 msec during the mission. Position of each FFM with respect to the rocket is calculated by the knowledge of its release velocity (measured at exit point of the FFM launcher tract) providing an accuracy of 1 m over the maximum range of 3 km. Spatial and temporal nature of observants can be separated to within 3 m in space or 3 msec time interval.

Javadi, H.↗
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↗
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.↗
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.↗