Search NASA⌕ Search

SEARCH · Search NASA

Results for “Relative agreement model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Revisiting Reentry Roll Resonance

Here, this paper revisits the subject of reentry roll resonance, which occurs when the roll rate and pitch frequency become equal in a spinning reentry vehicle. A first-principles model is developed and analyzed numerically and by using a reduced-order pendulum model for the relative phase of the roll and pitch/yaw axes. The results show good qualitative agreement with recent flight data and provide a natural physical interpretation. Based on the model it is suggested that sustained roll resonance depends primarily on the initial phase relationship between the pitch and roll axes, in addition to the vehicle moments of inertia and the vehicle symmetry. Finally the model is used to propose a mechanism for sustained resonance in vehicles spinning about a major axis, a condition previously thought to be forbidden.

33 ADVANCED PROPULSION SYSTEMS↗
NASA NTRS2007

Global Water Cycle Agreement in the Climate Models Assessed in the IPCC AR4

This study examines the fidelity of the global water cycle in the climate model simulations assessed in the IPCC Fourth Assessment Report. The results demonstrate good model agreement in quantities that have had a robust global observational basis and that are physically unambiguous. The worst agreement occurs for quantities that have both poor observational constraints and whose model representations can be physically ambiguous. In addition, components involving water vapor (frozen water) typically exhibit the best (worst) agreement, and fluxes typically exhibit better agreement than reservoirs. These results are discussed in relation to the importance of obtaining accurate model representation of the water cycle and its role in climate change. Recommendations are also given for facilitating the needed model improvements.

climate change↗
DOE OSTI2023

Quantitative Systems Modeling for Critical Infrastructure Predictions in Climate Change: A National Defense Framework

Any national defense is dependent on the efficacy of the available physical infrastructure. Whatever degrades infrastructure, structurally, physically, architecturally, or esthetically, is of interest to the nation. Climate change is now a major significant factor of interest impinging on national critical infrastructure. The devastating effects of climate change have increasing pervasiveness throughout the world. The impact on critical infrastructure is of particular interest to researchers. In consonance with ASCE-ASME's special issue on critical infrastructure protection and resilience, this paper presents a systems-modeling approach for critical infrastructure and predictions in relation to climate change agreements in COP26. The paper addresses high-level critical systems-based assessment of the social, legal, economic, and technical nuances impinging on the viability of COP26 agreements. The paper was written using a multinational collaboration approach. In conclusion, the specific focus of national defense is used as the backdrop for the methodology of the paper.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Challenges of standard halo models in constraining galaxy properties from cosmic infrared background anisotropies

The halo model, combined with halo occupation distribution (HOD) prescriptions, is widely used to interpret cosmic infrared background (CIB) anisotropies and extract physical information about star-forming galaxies and their connection to large-scale structures. Recent CIB-specific implementations of the halo model have adopted more physical parameterizations. However, the extent to which these models can reliably recover meaningful physical parameters remains uncertain. We assessed whether the current parameterization of CIB halo models is sufficient to recover astrophysical quantities, such as star formation efficiency, η(M h , z), and halo mass at which the peak of star formation efficiency occurs, M max , when fit to mock data. We also assessed whether discrepancies arise from assumptions about galaxy emission (the HOD ingredients) or from more fundamental components in the halo model, such as bias and matter clustering. We fit the M21 CIB HOD model, implemented within the halo model framework, to mock CIB power spectra and star formation rate density (SFRD) data generated from the SIDES-Uchuu simulation, and compared the best-fit parameters to the known simulation inputs. We then repeated the analysis using a simplified version of the simulation (SSU), explicitly designed to match the HOD assumptions. A detailed comparison of model and simulation outputs was carried out to trace the origin of observed discrepancies. While the M21 HOD model provides a good fit to the mock data, it failed to recover the intrinsic parameters accurately, particularly the halo mass at which star formation efficiency peaks. This mismatch persists even when fitting data generated with the same model assumptions. We find strong agreement (within 5%) in the emission-related components (SFRD, emissivity), but observe a scale- and redshift-dependent offset exceeding 20% in the two-halo term of the CIB power spectrum. This likely arises from limitations in the treatment of halo bias and matter clustering within the linear approximation. Additionally, incorporating scatter in the SFR–halo mass relation and the spectral energy distribution (SED) templates significantly affects the shot noise (∼50%), but has only a modest impact (less than 10%) on the clustered component. These results suggest that recovering physical parameters from CIB clustering requires improvements to the cosmological ingredients of the halo model framework, such as adopting scale-dependent halo bias and nonlinear matter power spectra in addition to careful modeling of emission physics.

cosmic background radiation↗
NASA NTRS1986

Probing the inner regions of the coma by UV observations of CS in comets

New vibrational bands of CS have been identified in several comets. The relative intensities of these bands are in agreement with a simple model, based upon most of the CS radicals being in the lowest vibrational level. The present observational accuracy is not sufficient to decide between this simple model and an earlier radiative equilibrium model. The observed intensity ratios appear to be similar in different comets, which means that they cannot depend strongly upon the total gas production rates.

Butterworth, P. S.↗
NASA NTRS1994

Estimated SAGE II ozone mixing ratios in early 1993 and comparisons with Stratospheric Photochemistry, Aerosols and Dynamic Expedition measurements

An empirical time-series model for estimating ozone mixing ratios based on Stratospheric Aerosols and Gas Experiment II (SAGE II) monthly mean ozone data for the period October 1984 through June 1991 has been developed. The modeling results for ozone mixing ratios in the 10- to 30- km region in early months of 1993 are presented. In situ ozone profiles obtained by a dual-beam UV-absorption ozone photometer during the Stratospheric Photochemistry, Aerosols and Dynamics Expedition (SPADE) campaign, May 1-14, 1993, are compared with the model results. With the exception of two profiles at altitudes below 16 km, ozone mixing ratios derived by the model and measured by the ozone photometer are in relatively good agreement within their individual uncertainties. The identified discrepancies in the two profiles are discussed.

Yue, G. K.↗
NASA NTRS2013

Modeling X-Ray Binary Evolution in Normal Galaxies: Insights from SINGS

We present the largest-scale comparison to date between observed extragalactic X-ray binary (XRB) populations and theoretical models of their production. We construct observational X-ray luminosity functions (oXLFs) using Chandra observations of 12 late-type galaxies from the Spitzer Infrared Nearby Galaxy Survey (SINGS). For each galaxy, we obtain theoretical XLFs (tXLFs) by combining XRB synthetic models, constructed with the population synthesis code StarTrack, with observational star formation histories (SFHs). We identify highest-likelihood models both for individual galaxies and globally, averaged over the full galaxy sample. Individual tXLFs successfully reproduce about half of oXLFs, but for some galaxies we are unable to find underlying source populations, indicating that galaxy SFHs and metallicities are not well matched and/or XRB modeling requires calibration on larger observational samples. Given these limitations, we find that best models are consistent with a product of common envelope ejection efficiency and central donor concentration approx.. = 0.1, and a 50% uniform - 50% "twins" initial mass-ratio distribution. We present and discuss constituent subpopulations of tXLFs according to donor, accretor and stellar population characteristics. The galaxy-wide X-ray luminosity due to low-mass and high-mass XRBs, estimated via our best global model tXLF, follows the general trend expected from the L(sub X) - star formation rate and L(sub X) - stellar mass relations of Lehmer et al. Our best models are also in agreement with modeling of the evolution both of XRBs over cosmic time and of the galaxy X-ray luminosity with redshift.

Normal Galaxies↗
NASA NTRS2019

The Effects of High-Temperature Fractional Crystallization on Calcium Isotopic Composition

We report Ca isotope fractionation among a co-genetic suite of samples from the Mauritanian Guelb el Azib ultramafic-mafic-anorthosite complex (GAC), which represents the fractional crystallization sequence of an igneous magma chamber. We investigate how the composition of an evolving crystallizing liquid and the resultant mineralogy of co-crystallizing assemblages affects inter-phase Ca isotope fractionation. Because the GAC is an Archean metamorphosed complex, we first investigate the resistance of the Ca isotope signature to secondary hydrothermal alteration and to metamorphism. While we cannot conclude here that Ca isotopes are undisturbed at the mineral scale, we show that they are preserved at the bulk rock scale. This is adequate for our study where several samples are almost monomineralic at the bulk scale level. The δ(exp 44/40)Ca (relative to NIST standard SRM 915a) of GAC layers range from −1.53 to 1.61‰, with the earliest-formed ultramafic cumulate layer being isotopically heaviest and the later, more felsic layers being isotopically lighter. For the first time, we show that plagioclase is much more enriched in light Ca isotopes compared to olivine and (Ca-) pyroxene. Monte Carlo simulations suggest that Ca isotopes are fractionated among co-existing silicate minerals during fractional crystallization, possibly in relation with the residual liquid composition. In qualitative agreement with computational models based on first principles lattice dynamics, we observe that Ca isotope fractionation is mineralogically controlled and importantly, that the degree of fractionation can vary according to the CaO composition of the residual liquid. While previous studies have aimed to understand partial melting as a source of Ca isotope fractionation, our results suggest that fractional crystallization is also a source of Ca isotope variability among co-genetic samples.

Maria C. Valdes↗
NASA NTRS1994

TOPEX/POSEIDON tides estimated using a global inverse model

Altimetric data from the TOPEX/POSEIDON mission will be used for studies of global ocean circulation and marine geophysics. However, it is first necessary to remove the ocean tides, which are aliased in the raw data. The tides are constrained by the two distinct types of information: the hydrodynamic equations which the tidal fields of elevations and velocities must satisfy, and direct observational data from tide gauges and satellite altimetry. Here we develop and apply a generalized inverse method, which allows us to combine rationally all of this information into global tidal fields best fitting both the data and the dynamics, in a least squares sense. The resulting inverse solution is a sum of the direct solution to the astronomically forced Laplace tidal equations and a linear combination of the representers for the data functionals. The representer functions (one for each datum) are determined by the dynamical equations, and by our prior estimates of the statistics or errors in these equations. Our major task is a direct numerical calculation of these representers. This task is computationally intensive, but well suited to massively parallel processing. By calculating the representers we reduce the full (infinite dimensional) problem to a relatively low-dimensional problem at the outset, allowing full control over the conditioning and hence the stability of the inverse solution. With the representers calculated we can easily update our model as additional TOPEX/POSEIDON data become available. As an initial illustration we invert harmonic constants from a set of 80 open-ocean tide gauges. We then present a practical scheme for direct inversion of TOPEX/POSEIDON crossover data. We apply this method to 38 cycles of geophysical data records (GDR) data, computing preliminary global estimates of the four principal tidal constituents, M(sub 2), S(sub 2), K(sub 1) and O(sub 1). The inverse solution yields tidal fields which are simultaneously smoother, and in better agreement with altimetric and ground truth data, than previously proposed tidal models. Relative to the 'default' tidal corrections provided with the TOPEX/POSEIDON GDR, the inverse solution reduces crossover difference variances significantly (approximately 20-30%), even though only a small number of free parameters (approximately equal to 1000) are actually fit to the crossover data.

Egbert, Gary D.↗
NASA NTRS2014

Implementing Badhwar-O'Neill Galactic Cosmic Ray Model for the Analysis of Space Radiation Exposure

For the analysis of radiation risks to astronauts and planning exploratory space missions, accurate energy spectrum of galactic cosmic radiation (GCR) is necessary. Characterization of the ionizing radiation environment is challenging because the interplanetary plasma and radiation fields are modulated by solar disturbances and the radiation doses received by astronauts in interplanetary space are likewise influenced. A model of the Badhwar‐O'Neill 2011 (BO11) GCR environment, which is represented by GCR deceleration potential theta, has been derived by utilizing all of the GCR measurements from balloons, satellites, and the newer NASA Advanced Composition Explorer (ACE). In the BO11 model, the solar modulation level is derived from the mean international sunspot numbers with time‐delay, which has been calibrated with actual flight instrument measurements to produce better GCR flux data fit during solar minima. GCR fluxes provided by the BO11 model were compared with various spacecraft measurements at 1 AU, and further comparisons were made for the tissue equivalent proportional counters measurements at low Earth orbits using the high‐charge and energy transport (HZETRN) code and various GCR models. For the comparison of the absorbed dose and dose equivalent calculations with the measurements by Radiation Assessment Detector (RAD) at Gale crater on Mars, the intensities and energies of GCR entering the heliosphere were calculated by using the BO11 model, which accounts for time‐dependent attenuation of the local interstellar spectrum of each element. The BO11 model, which has emphasized for the last 24 solar minima, showed in relatively good agreement with the RAD data for the first 200 sols, but it was resulted in to be less well during near the solar maximum of solar cycle 24 due to subtleties in the changing heliospheric conditions. By performing the error analysis of the BO11 model and the optimization in reducing overall uncertainty, the resultant BO13 model corrects the fit at solar maxima as well as being accurate at solar minima. The BO13 model is implemented to the NASA Space Cancer Risk model for the assessment of radiation risks. Overall cumulative probability distribution of solar modulation parameters represents the percentile rank of the average interplanetary GCR environment, and the probabilistic radiation risks can be assessed for various levels of GCR environment to support mission design and operational planning for future manned space exploration missions.

Kim, Myung-Hee Y.↗
NASA NTRS2005

The QBO as Potential Amplifier and Conduit to Lower Altitudes of Solar Cycle Influence?

The solar cycle (SC) effect in the lower atmosphere has been linked observationally to the Quasi-biennial Oscillation (QBO), which is generated primarily by small-scale gravity waves (GW). Following up on a 2D study with our Numerical Spectral Model (NSM), we discuss here a 3D study in which we simulated the QBO under the influence of the SC. For a SC period of 10 years, the relative amplitude of radiative forcing is taken to vary exponentially with height, i.e., 0.2% at the surface, 2% at 50 km, 20% at 100 km and above. Applying spectral analysis to filter out and identify the SC signature, this model produces a relatively large modulation in the QBO of the lower stratospheric circulation, which is in qualitative agreement with the results obtained by Salby and Callaghan (2000) who analyzed zonal wind observations covering more than 40 years. The modeled SC modulation of the QBO extends to high latitudes where it produces temperature variations of about 1 K in the troposphere below 10 km. We report that the SC also generates in the model a hemispherically symmetric Equatorial Annual Oscillation (EAO, with 12-month period) that is largely confined to low latitudes. Under the influence of the GWs, the SC modulated EA0 propagates down into the lower stratosphere like the QBO. As is the case for the QBO, the energy of this EA0 is partially redistributed by the meridional circulation and planetary waves, presumably, to generate measurable SC signatures in the tropospheric temperature of the polar regions, which may be related to the so called Arctic Oscillation (Thompson and Wallace, 1998). The larger SC influence at higher altitudes is apparently transferred to the lower and denser regions of the atmosphere by tapping the momentum from the upward propagating GWs that drive the oscillations. The SC modulation of the QBO period could prove to be very effective in this process, as our earlier 2D study indicated. Further studies are needed, (1) to make sure that the SC effects are real and the numerical results are robust, and (2) to explore more fully the mechanism(s) that may conspire to amplify the SC effect. Quasi-decadal oscillations, generated internally by the QBO interacting with the seasonal cycles, may interfere with, or aid, the SC influence.

Mayr, Hans G.↗
DOE OSTI2025

Sensitivities of time-dependent temperature profile predictions for NSTX with the multi-mode model

The Multi-Mode Model (MMM) for turbulent transport was applied to a large set of well-analyzed discharges from the National Spherical Torus Experiment (NSTX) in order to evaluate its sensitivities to a wide range of plasma conditions. MMM calculations were performed for hundreds of milliseconds in each discharge by performing time-dependent predictive simulations with the 1.5D tokamak integrated modeling code TRANSP. A closely related study (Lestz et al 2025 Plasma Phys. Control. Fusion 67 105029) concluded that MMM predicted electron and ion temperature profiles that were in reasonable agreement with NSTX observations, generally outperforming a different reduced transport model, TGLF. This finding motivates the more thorough investigation of the characteristics of the MMM predictions conducted in this work. The simulations with MMM have electron energy transport dominated by electron temperature gradient modes in the examined discharges with relatively low plasma β (ratio of kinetic plasma pressure to magnetic field pressure) and high collisionality, transitioning to a mixture of different modes for higher β and lower collisionality. The thermal ion diffusivity predicted by MMM is much smaller than the neoclassical contribution, in line with previous experimental analysis of NSTX. Nonetheless, the electron and ion temperature profiles are coupled via collisional energy exchange and thus sensitive to which transport channels are predicted. The time-dependent simulations with MMM are robust to the simulation start time, converging to remarkably similar temperature profiles later during the discharge. MMM typically overpredicts confinement relative to NSTX observations, leading to the prediction of overly steep temperature profiles. Plasmas with spatially broader temperature profiles, higher plasma β, and longer energy confinement times tend to be predicted by MMM with better agreement with the experiment. As a result, these findings provide useful context for understanding the regime-dependent tendencies of MMM in anticipation of self-consistent, time-dependent predictive simulations of NSTX-U discharges with these same modeling tools.

MMM↗
DOE OSTI2022

Charge state distributions in dense plasmas

Charge state distributions in hot, dense plasmas are a key ingredient in the calculation of spectral quantities like the opacity. However, they are challenging to calculate, as models like Saha–Boltzmann become unreliable for dense, quantum plasmas. Here, we present a new variational model for the charge state distribution, along with a simple model for the energy of the configurations that includes the orbital relaxation effect. Comparison with other methods reveals generally good agreement with average atom-based calculations, the breakdown of the Saha–Boltzmann method, and mixed agreement with a chemical model. We conclude that the new model gives a relatively inexpensive, but reasonably high fidelity method of calculating the charge state distribution in hot dense plasmas, in local thermodynamic equilibrium.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Expansion-history preferences of DESI DR2 and external data

We explore the origin of the preference of Dark Energy Spectroscopic Instrument (DESI) Data Release 2 (DR2) baryon acoustic oscillation measurements and external data from cosmic microwave background (CMB) and type Ia supernovae (SNIa) that dark energy behavior departs from that expected in the standard cosmological model with vacuum energy (Λ ⁢CDM). In our analysis, we allow a flexible scaling of the expansion rate with redshift that nevertheless allows reasonably tight constraints on the quantities of interest, and adopt and validate a simple yet accurate compression of the CMB data that allows us to constrain our phenomenological model of the expansion history. We find that data consistently show a preference for a 3%–4% increase in the expansion rate at 𝑧 ≃ 0.7 relative to that predicted by the standard Λ⁢ CDM model, in excellent agreement with results from the less flexible (𝑤 0 ,𝑤 𝑎 ) parametrization which was used in previous analyses. Even though our model allows a departure from the best-fit Λ⁢ CDM model at zero redshift, we find no evidence for such a signal. We also find no evidence (at greater than 1⁢𝜎 significance) for a departure of the expansion rate from the Λ ⁢CDM predictions at higher redshifts for any of the data combinations that we consider. Altogether, our results strengthen the robustness of the findings using the combination of DESI, CMB, and SNIa data to dark-energy modeling assumptions.

Cosmological parameters↗
DOE OSTI2021

Lithium-ion battery physics and statistics-based state of health model

A pseudo-2d model using COMSOL Multiphysics® software is developed to simulate performance and performance degradation of Li-ion batteries consisting of layered and olivine cathodes with graphite anode when subjected to peak shaving grid service. Multiple degradation pathways are considered, including solid electrolyte interphase (SEI) formation and breakdown at the anode, cathode dissolution and its synergistic effect on SEI formation at the anode. The model is validated by simulating commercial cylindrical cell performance. A global model is developed to simulate performance across all chemistries, along with individual chemistry models using global model parameters as initial values. There is good agreement between these models for various optimization parameters such as SEI equilibrium potential, cathode dissolution exchange current density, solvent diffusivity in the SEI and SEI ionic conductivity. To circumvent time constraints related to the COMSOL model, a 0d global model is developed which fits data well and provides more clarity on differences in cathode dissolution exchange current density. Again, good agreement for various optimization parameters is obtained among the COMSOL global & individual chemistry models and the 0-d model. The lessons learned from the physics-based model is used to develop a top down statistics-based model using current, voltage and anode volumetric change per mole lithium intercalated, along with their interactions as degradation predictors. This model predicts out of sample degradation for multiple grid services and electric vehicle drive cycle with high accuracy and provides the pathway to develop an efficient battery management system combining machine learning and findings from physics-based computationally intensive algorithms.

Crawford, Aladsair J.↗
DOE OSTI2024

Ne viii in the Warm-hot Circumgalactic Medium of FIRE Simulations and in Observations

Abstract The properties of warm-hot gas around ∼ L * galaxies can be studied with absorption lines from highly ionized metals. We predict Ne viii column densities from cosmological zoom-in simulations of halos with masses in ∼10 12 and ∼10 13 M ☉ from the Feedback in Realistic Environments (FIRE) project. Ne viii traces the volume-filling, virial-temperature gas in ∼10 12 M ☉ halos. In ∼10 13 M ☉ halos the Ne viii gas is clumpier, and biased toward the cooler part of the warm-hot phase. We compare the simulations to observations from the COS Absorption Survey of Baryon Harbors (or CASBaH) and COS Ultraviolet Baryon Survey (or CUBS). We show that when inferring halo masses from stellar masses to compare simulated and observed halos, it is important to account for the scatter in the stellar-mass–halo-mass relation, especially at M ⋆ ≳ 10 10.5 M ☉ . Median Ne viii columns in the fiducial FIRE-2 model are about as high as observed upper limits allow, while the simulations analyzed do not reproduce the highest observed columns. This suggests that the median Ne viii profiles predicted by the simulations are consistent with observations, but that the simulations may underpredict the scatter. We find similar agreement with analytical models that assume a product of the halo gas fraction and metallicity (relative to solar) ∼0.1, indicating that observations are consistent with plausible circumgalactic medium temperatures, metallicities, and gas masses. Variants of the FIRE simulations with a modified supernova feedback model and/or active galactic nuclei feedback included (as well as some other cosmological simulations from the literature) more systematically underpredict Ne viii columns. The circumgalactic Ne viii observations therefore provide valuable constraints on simulations that otherwise predict realistic galaxy properties.

Wijers, Nastasha A. (ORCID:0000000163747185)↗
DOE OSTI2021

Divergent controls of soil organic carbon between observations and process-based models

Abstract The storage and cycling of soil organic carbon (SOC) are governed by multiple co-varying factors, including climate, plant productivity, edaphic properties, and disturbance history. Yet, it remains unclear which of these factors are the dominant predictors of observed SOC stocks, globally and within biomes, and how the role of these predictors varies between observations and process-based models. Here we use global observations and an ensemble of soil biogeochemical models to quantify the emergent importance of key state factors – namely, mean annual temperature, net primary productivity, and soil mineralogy – in explaining biome- to global-scale variation in SOC stocks. We use a machine-learning approach to disentangle the role of covariates and elucidate individual relationships with SOC, without imposing expected relationships a priori . While we observe qualitatively similar relationships between SOC and covariates in observations and models, the magnitude and degree of non-linearity vary substantially among the models and observations. Models appear to overemphasize the importance of temperature and primary productivity (especially in forests and herbaceous biomes, respectively), while observations suggest a greater relative importance of soil minerals. This mismatch is also evident globally. However, we observe agreement between observations and model outputs in select individual biomes – namely, temperate deciduous forests and grasslands, which both show stronger relationships of SOC stocks with temperature and productivity, respectively. This approach highlights biomes with the largest uncertainty and mismatch with observations for targeted model improvements. Understanding the role of dominant SOC controls, and the discrepancies between models and observations, globally and across biomes, is essential for improving and validating process representations in soil and ecosystem models for projections under novel future conditions.

58 GEOSCIENCES↗
NASA NTRS2011

Dust Processing in Supernova Remnants: Spitzer MIPS SED and IRS Observations

We present Spitzer MIPS SED and IRS observations of 14 Galactic Supernova Remnants previously identified in the GLIMPSE survey. We find evidence for SNR/molecular cloud interaction through detection of [OI] emission, ionic lines, and emission from molecular hydrogen. Through black-body fitting of the MIPS SEDs we find the large grains to be warm, 29-66 K. The dust emission is modeled using the DUSTEM code and a three component dust model composed of populations of big grains, very small grains, and polycyclic aromatic hydrocarbons. We find the dust to be moderately heated, typically by 30-100 times the interstellar radiation field. The source of the radiation is likely hydrogen recombination, where the excitation of hydrogen occurred in the shock front. The ratio of very small grains to big grains is found for most of the molecular interacting SNRs to be higher than that found in the plane of the Milky Way, typically by a factor of 2--3. We suggest that dust shattering is responsible for the relative over-abundance of small grains, in agreement with prediction from dust destruction models. However, two of the SNRs are best fit with a very low abundance of carbon grains to silicate grains and with a very high radiation field. A likely reason for the low abundance of small carbon grains is sputtering. We find evidence for silicate emission at 20 $\mu$m in their SEDs, indicating that they are young SNRs based on the strong radiation field necessary to reproduce the observed SEDs.

Hewitt, John W.↗