SEARCH · Engineering Papers

Results for “Graphite modeling”

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.

SourceResearch connectionFrom yearThrough year

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

At least 451 records · Page 25

All sources NASA NTRS DOE OSTI

NASA NTRS1991

Simulation studies for surfaces and materials strength

Investigations were carried out in the following three areas. First, a simulation study was conducted to investigate the early stages of a nucleation and growth mechanism taking place on diamond surfaces. The energetically most favorable binding sites were calculated along with the binding energies for up to three carbon atoms deposited on the planes of a diamond crystal. Secondly, calculations were carried out to analyze energies and structural properties for restructured patterns of the surface of a diamond. Two different model functions, developed recently for carbon, were employed in the calculations. Lastly, a crack propagation process for the graphitic basal plane, was investigated using a molecular dynamics technique. Interactions among carbon atoms in the system were calculated using the Tersoff potential energy function.

Halicioglu, Timur↗

NASA NTRS1996

The infrared spectrum of the Galactic center and the composition of interstellar dust

We have obtained 5-8 micrometers spectra of the Galactic center from the Kuiper Airborne Observatory at resolving powers of approximately 50, approximately 150, and approximately 300. These spectra show absorption features at 5.5, 5.8, 6.1, and 6.8 micrometers. Together with previously observed features in the 3 micrometers region, these features are compared with laboratory spectra of candidate materials. The 3.0 and 6.1 micrometers features are due to the OH stretching and bending variations of H2O and are well fitted by water of hydration in silicates (e.g., talc). The 3.0 micrometer band is equally well fitted by ice mixtures containing 30% H2O, but such mixtures do not provide a good fit to the observed 6.1 micrometer band. The 3.4 and 6.8 micrometers features are identified with the CH stretching and deformation modes in CH2 and CH3 groups in saturated aliphatic hydrocarbons. The 6.1 micrometer band shows a short wavelength shoulder centered on 5.8 micrometer, attributed to carbonyl (C double bond O) groups in this interstellar hydrocarbon dust component. Finally, the narrow 5.5 micrometer feature is also attributed to carbonyl groups, but in the form of metal carbonyls [e.g., Fe(CO)4]. We have derived column densities and abundances along the line of sight toward the Galactic center for the various identified dust components. This analysis shows that hydrocarbon grains contain only 0.08 of the elemental abundance of C and contribute only a relatively minor fraction (0.1) of the total dust volume. Most of the interstellar dust volume is made up of silicates (approximately 0.6). Small graphite grains, responsible for the 2200 angstroms bump, account for 0.07 of the total dust volume. The remaining one-quarter of the interstellar dust volume consists of a material(s) without strong IR absorption features. Likely candidates include large graphite grains, diamonds, or amorphous carbon grains, which all have weak or no IR active modes. Finally, various models for the origin of the hydrocarbon dust component of the interstellar dust are discussed. All of them face some problems in explaining the observations, in particular, the absence of the spectroscopic signature of hydrocrbon grains in sources associated with molecular clouds.

NASA Discipline Exobiology↗

NASA NTRS2006

Experimental Constraints on Ureilite Petrogenesis

This experimental study explores the petrogenesis of ureilites by a partial melting/smelting process. Experiments have been performed over temperature (1150-1280 C), pressure (5-12.5 MPa), and low oxygen fugacity (graphite-CO gas) conditions appropriate for a hypothetical ureilite parent body approximately 200 km in size. Experimental and modeling results indicate that a partial melting/smelting model of ureilite petrogenesis can explain many of the unique characteristics displayed by this meteorite group. Compositional information preserved in the pigeonite-olivine ureilites was used to estimate the composition of melts in equilibrium with the ureilites. The results of 20 experiments saturated with olivine, pyroxene, metal, and liquid with appropriate ureilite compositions are used to calibrate the phase coefficients and pressure-temperature dependence of the smelting reaction. The calibrated coefficients are used to model the behavior of a hypothetical residue that is experiencing fractional smelting. The residue is initially olivine-rich and smelting progressively depletes the olivine content and enriches the pyroxene and metal contents of the residues. The modeled residue composition at 1260 C best reproduces the trend of ureilite bulk compositions. The model results also indicate that as a ureilite residue undergoes isothermal decompression smelting over a range of temperatures, Ca/Al values and Cr203 contents are enriched at lower temperatures (below about 1240 C) and tend to decrease at higher temperatures. Therefore, fractional smelting can account for the high Ca/A1 and Cr203 wt% values observed in ureilites. We propose that ureilites were generated from an olivine-rich, cpx-bearing residue. Smelting began when the residue was partially melted and contained liquid, olivine, and carbon. These residues experienced varying degrees of fractional smelting to produce the compositional variability observed within the pigeonite-bearing ureilites. Variations in mineral composition, modal proportions, and isotopic signatures are best described by heterogeneous accretion of the ureilite parent body followed by minimal and variable degrees of igneous processing.

Singletary, Steven↗

NASA NTRS1995

Graphite fiber textile preform/copper matrix composites

Graphite fiber reinforced/copper matrix composites have sufficiently high thermal conduction to make them candidate materials for critical heat transmitting and rejection components. The term textile composites arises because the preform is braided from fiber tows, conferring three-dimensional reinforcement and near net shape. The principal issues investigated in the past two years have centered on developing methods to characterize the preform and fabricated composite and on braidability. It is necessary to have an analytic structural description for both processing and final property modeling. The structure of the true 3-D braids used is complex and has required considerable effort to model. A structural mapping has been developed as a foundation for analytic models for thermal conduction and mechanical properties. The conductivity has contributions both from the copper and the reinforcement. The latter is accomplished by graphitization of the fibers, the higher the amount of graphitization the greater the conduction. This is accompanied by an increase in the fiber modulus, which is desirable from a stiffness point of view but decreases the braidability; the highest conductivity fibers are simply too brittle to be braided. Considerable effort has been expended on determining the optimal braidability--conductivity region. While a number of preforms have been fabricated, one other complication intervenes; graphite and copper are immiscible, resulting in a poor mechanical bond and difficulties in infiltration by molten copper. The approach taken is to utilize a proprietary fiber coating process developed by TRA, of Salt Lake City, Utah, which forms an itermediary bond. A number of preforms have been fabricated from a variety of fiber types and two sets of these have been infiltrated with OFHC copper, one with the TRA coating and one without. Mechanical tests have been performed using a small-scale specimen method and show the coated specimens to have superior mechanical properties. Final batches of preforms, including a finned, near net shape tube, are being fabricated and will be infiltrated before summer.

Gilatovs, G. J.↗

DOE OSTI2026

SAM Code Enhancements for Fission Product Tracking of Noble Gases and Metals in MSRs

This report documents fiscal year 2026 enhancements to the System Analysis Module (SAM) for modeling fission product transport in liquid-fueled molten salt reactors (MSRs). The work advances three principal areas: noble gas transport, noble metal deposition, and user interface improvements. The noble gas transport capability integrates drift-flux gas transport, Henry’s law two-film interphase mass transfer with pressure-based nucleation suppression, Knudsen-regime pore diffusion into porous graphite with a conjugate salt-graphite interface constraint, built-in material properties, five Sherwood-number mass transfer correlations including three derived from high-fidelity NekRS simulations, and xenon-135 reactivity feedback through SAM’s point-kinetics model. This work also presents a comprehensive verification test suite, including new analytically verified cases for pressure-dependent onset of interphase gas transfer in a stagnant vertical pipe, a postulated FLiBe-graphite Xe extraction permeator, a gravity riser with a fission-product source, and a descending pipe with gas redissolution driven by hydrostatic pressure. A machine learning framework for bubble rise velocity prediction in molten salt systems is developed and benchmarked on molten-salt and diverse aqueous bubble datasets. The best-performing fine-tuned transfer-learning networks achieve an 82% reduction in RMSE relative to the Clift correlation, and is implemented directly in SAM. The noble metal transport capability is developed, including a liquid-wall deposition model and a gas-surface flotation mechanism that transfers insoluble particles entrained by sparging gas to wetted structures. Verification tests and demonstration cases cover the surface deposition, flotation efflux, and flotation shedding. Finally, a new [SpeciesTransport] input structure replaces positional global vectors with selfcontained, order-independent, named species blocks, simplifies the specification of multiphase species and decay chains, and remains fully compatible with existing SAM input files. Together, these developments improve the physical fidelity, verification basis, and usability of SAM for system-level analyses of fissionproduct behavior in MSRs.

Mui, Travis (ORCID:0000000303736470)↗

DOE OSTI2023

Informing air–carbon ablation modeling with theoretical calculations of atomic oxygen and nitrogen interacting with carbon surfaces

To understand the gas–surface chemistry above the thermal protection system of a hypersonic vehicle, it is necessary to map out the kinetics of key elementary reaction steps. In this work, extensive periodic density functional theory (DFT) calculations are performed to elucidate the interaction of atomic oxygen and nitrogen with both the basal plane and edge sites of highly oriented pyrolytic graphite (HOPG). Reaction energies and barriers are determined for adsorption, desorption, diffusion, recombination, and several reactions. These DFT results are compared with the most recent finite-rate model for air–carbon ablation. Our DFT results corroborated some of the parameters used in the model but suggest that further refinement may be necessary for others. The calculations reported here will help to establish a predictive kinetic model for the complex reaction network present under hypersonic flight conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Nano-engineering the material structure of preferentially oriented nano-graphitic carbon for making high-performance electrochemical micro-sensors

Direct synthesis of thin-film carbon nanomaterials on oxide-coated silicon substrates provides a viable pathway for building a dense array of miniaturized (micron-scale) electrochemical sensors with high performance. However, material synthesis generally involves many parameters, making material engineering based on trial and error highly inefficient. Here, we report a two-pronged strategy for producing engineered thin-film carbon nanomaterials that have a nano-graphitic structure. First, we introduce a variant of the metal-induced graphitization technique that generates micron-scale islands of nano-graphitic carbon materials directly on oxide-coated silicon substrates. A novel feature of our material synthesis is that, through substrate engineering, the orientation of graphitic planes within the film aligns preferentially with the silicon substrate. This feature allows us to use the Raman spectroscopy for quantifying structural properties of the sensor surface, where the electrochemical processes occur. Second, we find phenomenological models for predicting the amplitudes of the redox current and the sensor capacitance from the material structure, quantified by Raman. Our results indicate that the key to achieving high-performance micro-sensors from nano-graphitic carbon is to increase both the density of point defects and the size of the graphitic crystallites. Our study offers a viable strategy for building planar electrochemical micro-sensors with high-performance.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

NASA NTRS1976

Chemistry of primitive solar material

The paper reviews chemical processes that occurred in the cooler outer regions of the primitive solar nebula (PSN) at the time of intimate chemical contact between the preplanetary condensate and the nebular gas. The elemental composition of the PSN is discussed, the 15 most abundant elements in it are listed, and numerical models of it are examined. Various condensation models are described and tested against observed properties of the planets, their satellites, and the asteroids. The chemistry of abundant volatile elements in the PSN is investigated along with stability limits of graphite in a solar-composition gas, regions of dominance of the most abundant carbon-containing gas species in the same gas, and implications of the moon's composition for its origin. Some theories that have been proposed as alternatives to the condensation models are noted.

Barshay, S. S.↗

DOE OSTI2022

Computational investigation of incident ion angles and material erosion at rough graphite and silicon carbide divertor surfaces

We present a computational investigation of the dependence of material erosion on the incident ion angle at rough graphite and silicon carbide divertor surfaces. Ion angle distributions (IADs) for D plasmas at NSTX-U and DIII-D divertors were calculated by an equation-of-motion model that traces the ion trajectories in the sheath. Then, the effective sputtering yields and ion shadowed area fractions were calculated by a Monte Carlo micro-patterning and roughness code that applied the calculated IADs to surface topographic data that were obtained from optical confocal microscopy of rough graphite and SiC divertor surfaces from NSTX-U and DIII-D experiments. The calculations found that the effective sputtering yields, the sputtering pattern, and the shadowed area are determined by the detailed surface topology rather than the root mean square roughness R RMS , which represents deviations from a flat surface. The suppression of the effective sputtering yields for rough surfaces compared to the yield for a smooth surface was accounted for by the change of the mean local incident ion angle (LIIA) (θ'). The mean surface inclination angle distribution (SIAD) (δ) was found to be a useful parameter to estimate the LIIA from the calculated IADs. We report global empirical formulas for the mean LIIA and fraction of the area shadowed from the main ions for D plasmas for rough surfaces with B-field incident angles α = 85°–89° as a function of the mean SIAD (δ). We propose the use of the mean LIIA (θ') to estimate the sputtering yield for rough surfaces from the angular dependence of the sputtering yield.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRSDate not supplied

Macro-scale Testing and Micromechanics Modeling of Fracture Behaviors for Boron Carbide Composites with Hierarchical Microstructures

With a unique combination of properties including high hardness, low density, wear and corrosion resistance, thermal stability, high neutron absorption, and semi-conductivity, boron carbide (B4C) is a candidate material for various engineering applications that involve extreme environment. The current applications of boron carbide, however, are limited by its intrinsic brittleness due to its strong covalent bonding. To toughen boron carbide, in this work hierarchical microstructure designs was used to provide multiple toughening mechanisms including crack deflection/bridging, micro-crack toughening, etc. Using field assisted sintering, B4C composites with hierarchical microstructure features including graphite platelets, micron and sub-micron sized TiB2 reinforcements were fabricated. The fracture toughness of fabricated B4C composites were previously measured at micro-scale using micro-indentation followed by post-testing microstructure inspection. However, questions including whether the fracture toughness enhancement measured at micro-scale can translate to macro-scale mechanical properties, and what are the fundamental mechanisms behind observed fracture toughness enhancement, remain to be answered. In this study, the fabricated B4C composites were tested using standardized four-point bending method to obtain fracture toughness at macro-scale. In addition, micromechanics modeling was conducted using MAC/GMC code and crack-band model to study the effect of residual stress and weak interphases on fracture behaviors of B4C composites reinforced with TiB2 particles. Through standardized four-point bending tests, fracture toughness enhancements up to 2.85, 3.32, and 3.65 MPa∙m(exp 1/2) (from 2.38 MPa∙m(exp 1/2)) were achieved for B4C composites with graphite platelets addition (micro/nano B4C), with TiB2 formation (micro B4C-TiB2), and with both graphite and TiB2 addition (micro/nano B4C-TiB2) respectively. Micromechanics modeling indicated that introduction of thermal residual stress and weak interphases caused enhanced micro-cracking behavior and resulted in the observed fracture toughness enhancement. These results furthered understanding about the mechanical behaviors at macro-scale and the mechanisms behind observed fracture toughness enhancement for B4C composites with hierarchical microstructures and can provide reference data for the future design of B4C composites with optimized microstructures for further fracture toughness enhancement.

Boron Carbide↗

NASA NTRS1997

Formation of Carbon Nanotube Based Gears: Quantum Chemistry and Molecular Mechanics Study of the Electrophilic Addition of o-Benzyne to Fullerenes, Graphene, and Nanotubes

Considerable progress has been made in recent years in chemical functionalization of fullerene molecules. In some cases, the predominant reaction products are different from those obtained (using the same reactants) from polycyclic aromatic hydrocarbons (PAHs). One such example is the cycloaddition of o-benzyne to C60. It is well established that benzyne adds across one of the rings in naphthalene, anthracene and other PAHs forming the [2+4] cycloaddition product (benzobicyclo[2.2.2.]-octatriene with naphthalene and triptycene with anthracene). However, Hoke et al demonstrated that the only reaction path for o-benzyne with C60 leads to the [2+2] cycloaddition product in which benzyne adds across one of the interpentagonal bonds (forming a cyclobutene ring in the process). Either reaction product results in a loss of aromaticity and distortion of the PAH or fullerene substrate, and in a loss of strain in the benzyne. It is not clear, however, why different products are preferred in these cases. In the current paper, we consider the stability of benzyne-nanotube adducts and the ability of Brenner's potential energy model to describe the structure and stability of these adducts. The Brenner potential has been widely used for describing diamondoid and graphitic carbon. Recently it has also been used for molecular mechanics and molecular dynamics simulations of fullerenes and nanotubes. However, it has not been tested for the case of functionalized fullerenes (especially with highly strained geometries). We use the Brenner potential for our companion nanogear simulations and believe that it should be calibrated to insure that those simulations are physically reasonable. In the present work, Density Functional theory (DFT) calculations are used to determine the preferred geometric structures and energetics for this calibration. The DFT method is a kind of ab initio quantum chemistry method for determining the electronic structure of molecules. For a given basis set expansion, it is comparable in accuracy to the MP2 method (better than Hartree Fock, but less accurate than more extensive electron correlation methods such as MP4 or CCSD). However, for systems with large numbers of basis functions it more efficient than any other methods that include electron correlation effects. In this presentation we show the results of DFT calculations for the reaction of benzyne with naphthalene, C60, and nanotube models. We compare energies for [2+2] and [2+4] cycloaddition products. The preferred products for the naphthalene and C60 reactions have been determined by experiment and, thus, these cases serve as a validation of our quantum chemical approach. We also compare the DFT and Brenner potential results. Finally we can predict the likelihood of reaction between benzyne and nanotubes.

Jaffe, Richard↗

DOE OSTI2024

Molecular simulation using transfer-learned potentials for the disordered nanoscale structure of nitrogen-doped nanoporous carbons

Machine learning (ML)-based molecular dynamics (MD) simulations of the formation of a class of N-doped nanoporous carbons are performed to assess their disordered partially graphitized nanoscale structure. The study is motivated by the effectiveness of so-called nitrogen assembly carbons (NACs) for catalysis applications. Benchmark simulations for pure-C disordered graphitic systems reveal the importance of reliably capturing the vdW component of the potentials in order to accurately describe the tendency for layering of disordered graphene-like sheets. In our modeling, this is achieved by a transfer learning strategy incorporating features of the energetics from the optB88-vdW DFT functional into potentials initially trained with a less expensive functional, thereby providing a superior description of the pure-C systems. Generation from MD simulations of realistic partially graphitized structures is significantly more challenging for N-doped versus for pure C systems. However, such structures are achieved by a tailored MD simulation protocol mimicking the experimental synthesis process and in particular incorporating an annealing and subsequent quenching stages. Simulated PXRD patterns effectively reproduce the features of experimental observations for NACs, including the appearance of a prominent but broad (002) peak at around 25, and the development of another weaker feature associated with in-layer ordering of mixed C-N graphene-like sheets.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2003

Molecular modeling of polymer composite-analyte interactions in electronic nose sensors

We report a molecular modeling study to investigate the polymer-carbon black (CB) composite-analyte interactions in resistive sensors. These sensors comprise the JPL electronic nose (ENose) sensing array developed for monitoring breathing air in human habitats. The polymer in the composite is modeled based on its stereoisomerism and sequence isomerism, while the CB is modeled as uncharged naphthalene rings with no hydrogens. The Dreiding 2.21 force field is used for the polymer, solvent molecules and graphite parameters are assigned to the carbon black atoms. A combination of molecular mechanics (MM) and molecular dynamics (NPT-MD and NVT-MD) techniques are used to obtain the equilibrium composite structure by inserting naphthalene rings in the polymer matrix. Polymers considered for this work include poly(4-vinylphenol), polyethylene oxide, and ethyl cellulose. Analytes studied are representative of both inorganic and organic compounds. The results are analyzed for the composite microstructure by calculating the radial distribution profiles as well as for the sensor response by predicting the interaction energies of the analytes with the composites. c2003 Elsevier Science B.V. All rights reserved.

NASA Center JPL↗

NASA NTRS2007

Frequency Response of an Aircraft Wing with Discrete Source Damage Using Equivalent Plate Analysis

An equivalent plate procedure is developed to provide a computationally efficient means of matching the stiffness and frequencies of flight vehicle wing structures for prescribed loading conditions. Several new approaches are proposed and studied to match the stiffness and first five natural frequencies of the two reference models with and without damage. One approach divides the candidate reference plate into multiple zones in which stiffness and mass can be varied using a variety of materials including aluminum, graphite-epoxy, and foam-core graphite-epoxy sandwiches. Another approach places point masses along the edge of the stiffness-matched plate to tune the natural frequencies. Both approaches are successful at matching the stiffness and natural frequencies of the reference plates and provide useful insight into determination of crucial features in equivalent plate models of aircraft wing structures.

Krishnamurthy, T.↗

DOE OSTI2021

Modeling Interface Debonding in Coated Fuel Particles with BISON

Here we introduce a meshed debonding model intended for use with tristructural isotropic (TRISO) fuel particles. The ability to better examine the potential effects of interface debonding between the coating layers will enable a host of investigations of the stress state when complex loading and interactions between coated fuel particles are taking place. In this work, the debonding model is discussed in detail. Then, simulations are presented to show the capabilities and limitations of the model. First, a traditional high-temperature gas reactor TRISO particle is considered. Then, a TRISO particle in a graphitic matrix is analyzed, followed by a TRISO particle in a silicon carbide matrix.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

Engineering in 5 “Deimos, an Advanced Reactor Criticality Experiment” [Slides]

Problem: Advanced nuclear reactors come in many sizes, materials, and designs; Criticality experiments are necessary to understand, quantify, and validate these reactors. Approach: The National Criticality Experiments Research Center (NCERC) operated by Los Alamos National Laboratory (LANL) conducts criticality experiments; Deimos is a graphite-moderated, beryllium-reflected, high assay low-enriched uranium (HALEU) tri-structural isotropic (TRISO) fueled experiment; Modeling and simulation is needed to design Deimos.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

NASA NTRS1974

Measurement of the X-ray emission of the solar atmosphere during a period of low activity

A large-area high-sensitivity X-ray spectrometer has been constructed and used to measure the 1.8- to 5.3-A X-ray emission of the sun under quiescent conditions. The instrument utilizes Bragg reflection from mosaic graphite crystals. The data indicate that the X-ray emission can best be accounted for by a multitemperature model of the solar atmosphere in which both the overall corona and active regions contribute to the X-ray spectrum. Theoretical calculations of the X-ray flux of a hot, optically thin plasma have been used to estimate the solar conditions at the time when the measurements were made.

Wolff, R. S.↗

NASA NTRS1980

Experimental determination of unsteady blade element aerodynamics in cascades. Volume 2: Translation mode cascade

A two dimensional cascade of harmonically oscillating airfoils was designed to model a near tip section from a rotor which was known to have experienced supersonic translational model flutter. This five bladed cascade had a solidity of 1.52 and a setting angle of 0.90 rad. Unique graphite epoxy airfoils were fabricated to achieve the realistic high reduced frequency level of 0.15. The cascade was tested over a range of static pressure ratios approximating the blade element operating conditions of the rotor along a constant speed line which penetrated the flutter boundary. The time steady and time unsteady flow field surrounding the center cascade airfoil were investigated.

Riffel, R. E.↗