SEARCH · Engineering Papers
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.
Structure of the program Graphite material issues of interest How do things fit together Status of the program Baseline AGC Experiment Update AGC-4 and HDG-1 Irradiation Oxidation Modeling ASME Code development Discussion of how the data can be used Modeling behavior of core components Licensing/code development Interaction with vendors Collaborations (domestic and international)
A model has been developed for reflection nebulae around luminous IR sources embedded in dense dust clouds. The shape of the IR spectrum is shown to be the result of a combination of the scattering properties of the dust, the spectrum of the illuminating source, and foreground extinction, while geometry plays a minor role. Comparison of the model results with IR observations of the reflection nebula surrounding OMC-2/IRS 1 shows that either a grain size distribution like that found in the diffuse ISM, or consisting of larger grains, can explain the observed shape of the spectrum. However, the absolute intensity level of the scattered light, as well as the observed polarization, requires large grains. By adding water-ice mantles to the silicate and graphite cores, the 3.08 micron ice-band feature observed in the spectra of several IR reflection nebulae has been modeled. It is shown that this ice band arises naturally in optically thick reflection nebulae containing ice-coated grains.
Electrocatalysts for the oxygen reduction reaction (ORR) based on complexes of iron and nitrogen in a carbon matrix (Fe-N-C) are a promising alternative to platinum group metal (PGM) based catalysts in polymer electrolyte membrane (PEM) fuel cells. Further improvements of Fe-N-C catalysts would benefit from model thin film studies of activity and stability of catalytic sites, but synthesis of Fe-N-C model thin films is challenging. Here we report on synthesis and characterization of Fe-N-C thin films produced by co-sputtering iron and carbon in a reactive nitrogen atmosphere onto removable glassy carbon rotating disk electrode (RDE) tips. Scanning electron microscopy (SEM) measurements indicate that the Fe-N-C films deposited at high temperature are smoother than the films annealed at high temperature. Electrocatalytic activity measured on the thin Fe-N-C films is greater for both high-temperature samples than for the room-temperature sample. From the analysis of X-ray photoelectron spectroscopy (XPS) data, exposure of the films to high temperatures results in increased graphitization of the carbon with the Fe-N-C films, and increased relative amount of graphitic and hydrogenated nitrogen species. Overall the results of this study demonstrate the feasibility of a thin film model system approach for studying active sites in PGM-free catalysts.
Particle fuel forms generally consist of a dispersion of fuel, such as tristructural isotropic (TRISO) particles, within a refractory matrix (e.g., graphite or silicon carbide). The density of matrix materials for particle fuel forms is of interest for modeling fuel form strength and thermal properties and may be specified as a quality control parameter, depending on reactor design. Some of the uncertainty associated with traditional, manual approaches can be eliminated by performing x-ray computed tomography (XCT) on the fuel forms and applying image processing methods to generate a precise count of the number of particles. This also removes the need to include determination of particle count within each individual fuel form during fabrication. Unfortunately, reconstruction artifacts from high-Z uranium-bearing kernels prevent accurate measurement of individual particle volumes using this approach, so the use of mean particle mass and volume are still necessary for computation of average fuel form matrix density. This method of using XCT to count particles in individual fuel form for determination of average matrix density was applied to three archived compacts from the Advanced Gas Reactor Fuel Development and Qualification (AGR)-1 campaign, four archived compacts with uranium carbide/uranium oxide (UCO) TRISO from the AGR-2 campaign, and three archived UO 2 -TRISO compacts from the AGR-2 campaign. The resulting density values were compared with those previously reported, showing slight changes due to uncertainties in the previously used number of particles in each of these cylindrical, graphite matrix compacts.
An orthotropic double cantilever beam (DCB) model is used to study dynamic crack propagation and arrest in 90 deg unidirectional Hercules AS/3501-6 graphite fiber epoxy composites. The dynamic fracture toughness of the composite is determined from tests performed on the long-strip specimen and DCB crack arrest experiments are conducted. By using the dynamic fracture toughness in a finite-difference solution of the DCB governing partial differential equations, a numerical solution of the crack propagation and arrest events is computed. Excellent agreement between the experimental and numerical crack arrest results are obtained.
The predictions of the model of Puget et al. (1985) for the emission from Very Small Grains (VSGs) including both graphitic and silicate components are compared with published 8-13-micron observations of astronomical sources. The VSGs are found to be mainly graphitic and an upper limit is placed on the relative mass of silicates based on lack of the 9.7-micron silicate emission feature on M 82 and NGC 2023. This dissymetry in the composition of VSGs supports the suggestion that they are formed in grain-grain collisions where the behaviors of graphite and silicate grains are expected to be quite different.
Directional solidification experiments in a Bridgman-type furnace were used to study particle behavior at the liquid/solid interface in aluminum metal matrix composites. Graphite or silicon-carbide particles were first dispersed in aluminum-base alloys via a mechanically stirred vortex. Then, 100-mm-diameter and 120-mm-long samples were cast in steel dies and used for directional solidification. The processing variables controlled were the direction and velocity of solidification and the temperature gradient at the interface. The material variables monitored were the interface energy, the liquid/particle density difference, the particle/liquid thermal conductivity ratio, and the volume fraction of particles. These properties were changed by selecting combinations of particles (graphite or silicon carbide) and alloys (Al-Cu, Al-Mg, Al-Ni). A model which consideres process thermodynamics, process kinetics (including the role of buoyant forces), and thermophysical properties was developed. Based on solidification direction and velocity, and on materials properties, four types of behavior were predicted. Sessile drop experiments were also used to determine some of the interface energies required in calculation with the proposed model. Experimental results compared favorably with model predictions.
A method of predicting the crack-related energy-absorption capability of composite tubes is presented. The method is based upon a phenomenological model of the crushing process exhibited by continuous-fiber-reinforced tubes. A finite element method is used to model the crushing process. The analysis is compared with experiments on Kevlar-epoxy and graphite-epoxy tubes. Reasonable agreement is obtained between the analysis and experiment.
Interstellar dust is characterized by strong absorption in the ultraviolet and the mid-IR. Current models of interstellar dust are based on three chemically distinct components: a form of carbon (usually graphite), a silicate, and a blend of polycyclic aromatic hydrocarbons or other carbonaceous material. Previous work using effective medium theories to understand the optical properties of cometary dust suggested that an amalgam of materials could reproduce the observed interstellar and cometary dust features. Recently, Lawler and Brownlee (1992) re-analyzed the PIA and PUMA-1 data sets from the Giotto flyby of P/Halley and discovered that the so-called "CHON" particles were actually composed of a blend of carbon-bearing and silicon-bearing materials. Based on effective medium theories, the absorption spectrum of such a material would display the spectral features of each of the components - strong UV absorption from the carbonaceous component and strong absorption in the IR from the silicate component. To test this idea, vapor-deposited samples were created using two different deposition techniques: sputtering with an argon RF magnetron and deposition from an argon plasma torch. Two different compositions were tested: a blend of graphite and silica in a 7:1 ratio and an amalgam of materials whose approximate composition matches the "CHON"-silicate abundances for the uncompressed PIA data set of Lawler and Brownlee: graphite, iron oxide, magnesium oxide, ammonium sulfate, calcium carbonate, and silica in mass ratios of 6:4.3:4:2.2:1:9. The samples were finely ground and pressed into 2" diameter disks using a 40 ton press. In all, four different experiments were performed: one with each of the compositions (C:SiO and "CHON") in both the RF magnetron and the plasma torch chambers. The RF magnetron created a uniform dark thin film on the substrate surface, and the plasma torch created a coating of small (<100 micron) diameter grey particles. The spectra of all four samples show a strong, broad absorption feature at around 220 nm as well as a strong but narrower absorption peak near 10 microns. The RF magnetron sputtered samples showed some sub-structure in the UV, and the peak of the absorption was shifted toward longer wavelengths. The UV absorption in the plasma torch deposited samples have no sub-structure, and the peak absorption is very near 220 nm. Strong absorption near 9 microns is seen in the spectra from both sample preparation techniques, and is consistent with the IR spectra of some terrestrial silicates. Other features, particularly at 6.2 and 8.6 microns, are seen in the interstellar medium. A strong feature near 2 microns is due to absorbed water in the sample. Based on the results of these experiments, there is evidence that a material with a composition similar to that detected in "CHON" particles in the coma of P/Halley have a spectral signature which reproduces the main absorption features of interstellar dust. This suggests that the "CHON" particles could be the interstellar component of cometary dust.
We investigate how systematically increasing the accuracy of various molecular dynamics modeling techniques influences the structure and capacitance of ionic liquid electric double layers (EDLs). The techniques probed concern long-range electrostatic interactions, electrode charging (constant charge versus constant potential conditions), and electrolyte polarizability. Our simulations are performed on a quasi-two-dimensional, or slab-like, model capacitor, which is composed of a polarizable ionic liquid electrolyte, [EMIM][BF4], interfaced between two graphite electrodes. To ensure an accurate representation of EDL differential capacitance, we derive new fluctuation formulas that resolve the differential capacitance as a function of electrode charge or electrode potential. The magnitude of differential capacitance shows sensitivity to different long-range electrostatic summation techniques, while the shape of differential capacitance is affected by charging technique and the polarizability of the electrolyte. For long-range summation techniques, errors in magnitude can be mitigated by employing two-dimensional or corrected three dimensional electrostatic summations, which lead to electric fields that conform to those of a classical electrostatic parallel plate capacitor. With respect to charging, the changes in shape are a result of ions in the Stern layer (i.e. ions at the electrode surface) having a higher electrostatic affinity to constant potential electrodes than to constant charge electrodes. For electrolyte polarizability, shape changes originate from induced dipoles that soften the interaction of Stern layer ions with the electrode. The softening is traced to ion correlations vertical to the electrode surface that induce dipoles that oppose double layer formation. In general, our analysis indicates an accuracy dependent differential capacitance profile that transitions from the characteristic camel shape with coarser representations to a more diffuse profile with finer representations.
An explanation is offered for the extreme sensitivity that has been observed in the power law parameters of the T300/934 graphite epoxy material systems during experiments to evaluate the system's viscoelastic response. It is shown that the singularity associated with the power law can explain the sensitivity as well as the observed variability in the calculated parameters. Techniques for minimizing errors are suggested.
A bimodal model of interstellar grains consisting of silicate cores of size approximately 0.06 micron with modified ice mantles of approximately 0.1 micron and very small particles of silicate and/or graphite of size approximately 0.005 micron is followed through various stages of evolution. Starting with a distribution of core-mantle cylinders producing average wavelength dependence of polarization and extinction, changes in physical and optical characteristics of the model are studied going first into regions of dense cloud condensation and coming out at the other end of star formation into young H2 regions. Physical justification is presented to show that not only in dense clouds must the core mantle grains be larger than normal, but also in young H2 regions.
Silicon carbide (SiC) represents a promising but largely untested plasma-facing material (PFM) for next-step fusion devices. In this work, an analytic mixed-material erosion model is developed by calculating the physical (via SDTrimSP) and chemical (via empirical scalings) sputtering yield from SiC, Si, and C. The Si content in the near-surface SiC layer is predicted to increase during D plasma bombardment due to more efficient physical and chemical sputtering of C relative to Si. Silicon erosion from SiC thereby occurs primarily from sputtering of the enriched Si layer, rather than directly from the SiC itself. SiC coatings on ATJ graphite, manufactured via chemical vapor deposition, were exposed to repeated H-mode plasma discharges in the DIII-D tokamak to test this model. The qualitative trends from analytic modeling are reproduced by the experimental measurements, obtained via spectroscopic inference using the S/XB method. Quantitatively the model slightly under-predicts measured erosion rates, which is attributed to uncertainties in the ion impact angle distribution, as well as the effect of edge-localized modes. After exposure, minimal changes to the macroscopic or microscopic surface morphology of the SiC coatings were observed. Compositional analysis reveals Si enrichment of about 10%, in line with expectations from the erosion model. Extrapolating to a DEMO-type device, an order-of-magnitude decrease in impurity sourcing, and up to a factor of 2 decrease in impurity radiation, is expected with SiC walls, relative to graphite, if low C plasma impurity content can be achieved. Furthermore, these favorable erosion properties motivate further investigations of SiC as a low-Z, non-metallic PFM.
The third Advanced Gas Reactor (AGR) irradiation experiment, AGR-3/4, was designed to investigate the migration of fission products in fuel compact graphitic matrix and reactor graphite components. Using destructive methods, radial fission product concentration profiles were measured for gamma-emitting fission products (e.g., Ag-110m, Cs-134, and Eu-154) and beta-emitting Sr-90 in irradiated graphitic and graphite components from six different AGR-3/4 irradiation capsules. These new measured concentration profiles can now be compared to non-destructive measurements and fission product transport simulations and will be used to derive new diffusivities and sorptivities to support refinement of fission product transport models and high-temperature gas-cooled reactor (HTGR) source-term analyses. Each capsule in the AGR-3/4 experiment had four fuel compacts in the middle of two concentric rings of graphitic matrix material, PCEA graphite, or IG-110 graphite. In addition to the approximately 1898 tristructural isotropic (TRISO) coated particles in each compact, there were 20 designed to fail (DTF) particles coated only in pyrocarbon so that they released fission products into the surrounding cylindrical rings of carbonaceous materials. Destructive sampling of the rings involved machining/milling material from around the circumference of the rings, collecting that material, and performing radiochemical analyses on it. Milling operations were performed in multiple steps or segments, and each segment was generally 0.508 mm (0.020 in) thick. Knowing the radial position at which each segment was milled, the volume of the milled material at each segment, and the fission product content in each segment, the radial fission product concentrations were constructed for select isotopes in each ring. Ag-110m profiles had the most variation. Some profiles were peaked at an inner or outer surface. Some were peaked at the middle of the ring wall thickness. Some increased radially outward, and some decreased radially outward. These types of variations and the fact that the measured profiles do not generally compare favorably with the transport model employed for AGR-3/4 may adversely impact the ability to extract reasonable transport parameters for this isotope. In many cases, the Cs-134 profiles decreased somewhat linearly in the outward radial direction, and in cursory comparisons, the shapes of these profiles appeared similar to those from model predictions. The step changes in concentration across the inner-outer ring gap were generally consistent with the model predictions as well. In some cases, there were local maxima in concentration at the outer surface of the rings. This suggests that fission products could have transported in the small gaps between the inner ring and the outer ring and between the outer ring and the sink ring such that some portion of a given fission product can bypass diffusion through the ring itself. The analysis of the small nubs on the outer surfaces of some of the outer rings revealed fission product concentrations in the nubs that were often higher than in the outermost segments of the rings. This further supports the hypothesis that short-circuit, gap transport occurred, causing relatively high surface concentrations on the outer surfaces of the rings. Eu-154 and Sr-90 profiles tend to have very similar shapes, suggesting that they transport via the same mechanisms. The observed profiles were indicative of a transport process where the isotopes are sorbed on the inner surface of the ring, but diffusion into the ring from that surface is quite slow. Some elevated concentrations of Sr-90 (relative to Eu-154) on the outer surface of a ring suggested that rapid, gas-gap transport of gaseous precursor Kr-90 and volatile Rb-90 could have occurred prior to their decaying to Sr-90. Overall, the Eu 154 and Sr-90 profiles were still very similar, which indicates that the transport of short-lived Sr-90 precursors is not a major effect. In some capsules, the qualitative Sr-90 behavior across the ring gaps was consistent with the model (using the available legacy Sr-90 transport parameters), but in other capsules the model was inconsistent with the measurements and seems to underestimate the amount of Sr-90 in the outer rings. The total ring Sr-90 inventories were estimated for all the rings that were subject to physical sampling. These results will be used to adjust the predicted particle and/or compact releases used in the AGR-3/4 fission product transport model. Given the different irradiation temperatures among the capsules and the rings, it was not possible to discern fundamental differences in the transport of isotopes within the different carbon materials, i.e., graphitic matrix, IG-110, or PCEA. It may be possible to do this in the course of determining transport from the concentration profiles in future work.
Two micromechanical models were developed to investigate the thermal expansion of graphite/copper (Gr/Cu) composites. The models incorporate the effects of temperature-dependent material properties, matrix inelasticity, initial residual stresses due to processing history, and nonuniform fiber distribution. The first model is based on the multiple concentric cylinder geometry, with each cylinder treated as a two-phase composite with a characteristic fiber volume fractions. By altering the fiber volume fraction of the individual cylinders, unidirectional composites with radially nonuniform fiber distributions can be investigated using this model. The second model is based on the inelastic lamination theory. By varying the fiber content in the individual laminae, composites with nonuniform fiber distribution in the thickness direction can be investigated. In both models, the properties of the individual regions (cylinders or laminae) are calculated using the method of cells micromechanical model. Classical incremental plasticity theory is used to model the inelastic response of the copper matrix at the microlevel. The models were used to characterize the effects of nonuniform fiber distribution on the thermal expansion of Gr/Cu. These effects were compared to the effects of matrix plasticity, choice of stress-free temperature, and slight fiber misalignment. It was found that the radially nonuniform fiber distribution has little effect on the thermal expansion of Gr/Cu but could become significant for composites with large fiber-matrix transverse CTE and Young's modulus mismatch. The effect of nonuniform fiber distribution in the through-thickness direction of a laminate was more significant, but only approached that of the stress-free temperature for the most extreme cases that include large amounts of bending. Subsequent comparison with experimental thermal expansion data indicated the need for more accurate characterization of the graphite fiber thermomechanical properties. Correlation with cyclic data revealed the presence of a mechanism not considered in the developed models. The predicted response did, however, exhibit ratcheting behavior that has been observed experimentally in Gr/Cu. Finally, simulation of the actual fiber distribution of particular specimens had little effect on the predicted thermal expansion.
Radiation damage estimation is an important component of the post irradiation analysis of the Advanced Graphite Creep (AGC) experiment. It depends primarily on the fast fluence, which is determined using well established methods of spectral adjustment. These are based on best estimates from models such as Monte Carlo N-Particle (MCNP), input cross-sections, and measured activities from flux wires in the experiment. Each of these parameters can propagate uncertainties which will affect the uncertainty in the calculated dose levels for AGC, or any experiment irradiated within a reactor. While the methods of propagating uncertainty are well-established, the final uncertainty estimates they provide are only as good as the estimates of uncertainty in the inputs on which they are based. The purpose of this work is to outline some deficiencies in the ways these input uncertainties are presently estimated, and to outline a methodology by which they can be improved. The fast fluence and radiation damage received by graphite specimens irradiated in the Advanced Graphite Creep (AGC) experiments is presently estimated using spectral adjustment methods that are based on both flux wire activity measurements, and MCNP model predictions. This work describes an ongoing effort to quantify and propagate uncertainties in inputs to the spectral adjustment process, and thereby quantify the resultant error in radiation damage (dpa) estimates. The effort is multi-faceted, and we consider the impacts of both the set of flux wires selected, and the counting process. An expanded set of flux wires is identified that provides a more comprehensive data set on the fast spectrum. To address the counting process itself, a series of round-robin measurements in several reactor metrology laboratories across the Department of Energy (DOE) complex and nuclear industry are being undertaken to refine the American Society for Testing and Materials (ASTM) standards for flux wire measurements. To address the contribution of uncertainty in the MCNP model predictions, an uncertainty quantification (UQ) tool has been developed that statistically samples the model input parameters, runs a series of cases, and assimilates the results to provide an overall uncertainty. The impact of the MCNP UQ tool results is demonstrated by re-analyzing previous AGC flux wire and irradiation data. While the expanded flux wire set obviously cannot be added to these experiments retroactively, plans for future graphite irradiations are outlined.
Nuclear graphite is a critical material in high-temperature nuclear reactors due to its superior thermal and mechanical properties. The manufacturing process leaves multi-scale ‘pre-existing’ defects that can affect thermal transport characteristics. Because these defects are remnant of graphitization temperature, they cannot be thermally annealed. This study employs a non-thermal electron wind force (EWF) annealing technique to avoid this obstacle. 2 min of EWF treatment of the as-received graphite IG-110 at temperatures below 100 °C led up to 67% increase in thermal diffusivity and ~ 35% decrease in electrical resistivity in average. Differential scanning calorimetry also showed similar outcome for specific heat. X-ray diffraction characterization was performed by fitting a bi-modal crystallite size distribution model to reveal the enhancement in crystallinity after the EWF treatment. The findings emphasize the potential of EWF annealing for optimizing thermal performance in nuclear graphite and its implications for reactor efficiency and safety.
A rational methodology of modelling both nonlinear and elastic dissipative response of transversely isotropic fibrous composites is developed and illustrated with the aid of the observed response of graphite-polyimide off-axis coupons. The methodology is based on the internal variable formalism employed within the text of classical irreversible thermodynamics and entails extension of Valanis' endochronic theory to transversely isotropic media. Applicability of the theory to prediction of various response characteristics of fibrous composites is illustrated by accurately modelling such often observed phenomena as: stiffening reversible behavior along fiber direction; dissipative response in shear and transverse tension characterized by power-laws with different hardening exponents; permanent strain accumulation; nonlinear unloading and reloading; and stress-interaction effects.