Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “secondary phases”

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 397 records · Page 22

Performing spectroscopic and specific heat studies of improper ferroelectrics

The results of infrared measurements on Ni-Br, Cu-Cl, and Fe-I boracite improper ferroelectrics and far infrared measurements of Ni-Br boracite are presented. The boracites have the general formula X3B7O3Y, where X = divalent metal and Y = halogen. They undergo a first order phase transition from a high temperature paraelectric phase with cubic symmetry to a ferroelectric phase with orthorhombic symmetry. The boracites are "improper ferroelectrics" since the spontaneous polarization is not the primary order parameter in the cubic-orthorhombic phase transition. Current understanding of these materials is that the primary order parameter is associated with a doubly degenerate zone-boundary phonon in the cubic phase. The degenerate critical modes become homogeneous and split into the A sub 1 and A sub 2 modes in the orthorhombic phase, doubling the volume of the primitive cell. An harmonic coupling between the softing A sub 1 and a low frequency A sub 1 optic mode induces a spontaneous polarization as a secondary effect in the ferroelectric phase. This secondary non-critical nature of the ferroelectric mode earns these materials the "improper" title and is responsible for their unique properties and high figure of merit in detector use.

Coleman, L. B.↗

Discovery of charge order above room-temperature in the prototypical kagome superconductor La(Ru1−xFex)3Si2

Abstract The kagome lattice is an intriguing and rich platform for discovering, tuning and understanding the diverse phases of quantum matter, crucial for advancing modern and future electronics. Despite considerable efforts, accessing correlated phases at room temperature has been challenging. Using single-crystal X-ray diffraction, we discovered charge order above room temperature in La(Ru 1−x Fe x ) 3 Si 2 (x= 0, 0.01, 0.05), where charge order related to out-of-plane Ru atom displacements appears belowT CO,I ≃ 400 K. The secondary charge ordered phase emerges belowT CO,II ≃ 80–170 K. Furthermore, first principles calculations reveal both the kagome flat band and the van Hove point near the Fermi energy in LaRu 3 Si 2 , driven by Ru-dz 2 orbitals. Our results identify LaRu 3 Si 2 as the kagome superconductor with the highest known charge ordering temperature, offering a promising avenue for researching room temperature quantum phases and developing related technologies.

Physics↗

Implications of the Secondary Eclipse of Exoplanet HAT-P-11b

We observed exoplanet HAT-P-11b and have successfully detected its secondary eclipse. We conducted observations using the Spitzer Space Telescope in the post-cryo mission at 3.6 microns for a period of 22 hours centered on the anticipated secondary eclipse time, to detect the eclipse and determine its phase. Having detected the secondary eclipse, we are at present making a more focused series of observations in both the 3.6 and 4.5 micron bands to fully characterize it. HAT-P-11b is one of only two known exo-Neptunes and has a period of 4.8878 days, radius of 0.422 RJ, mass of 0.081 MJ and semi-major axis 0.053 AU. Measurements of the secondary eclipse will serve to clarify two key issues; 1) the planetary brightness temperature and the nature of its atmosphere, and 2) the eccentricity of its orbit, with implications for its dynamical evolution. We discuss implications of these observations.

Barry, Richard K.↗

Investigating Weathering of Basaltic Materials in Gale Crater, Mars: A Combined Laboratory, Modeling and Terrestrial Field Approach

Recent observations from Gale Crater, Mars document past aqueous alteration both in the formation of the Stimson sandstone unit, as well as in the formation of altered fractures within that unit. Geochemical and mineralogical data from Curiosity also suggest Fe-rich amorphous weathering products are present in most samples measured to date. Here we interpret conditions of possible past weathering in Gale Crater using a combination of field, laboratory, and modeling work. In order to better understand secondary Fe-rich phases on Mars, we are examining formation of weathering products in high Fe and Mg and low Al serpentine soils in the Klamath Mountains, CA. We have isolated potential weathering products from these soils, and are analyzing them using synchrotron μXRF and μXRD as well as FullPat for a direct comparison to analyses from Gale Crater. In order to interpret the implications of the persistence of potential secondary Fe-containing phases on Mars, we are also measuring the dissolution rates of the secondary weathering products allophane, Fe-rich allophane, and hisingerite. Ongoing dissolution experiments of these materials suggest that they dissolve significantly more rapidly than more crystalline secondary minerals with similar chemical compositions. Finally, to quantify the specific conditions of past aqueous alteration in Gale Crater we are performing reactive transport modeling of a range of possible past environmental conditions. Specifically, we are testing the conditions under which a Stimson unit-like material forms from a parent material similar to Rocknest or Bagnold eolian deposits, and the conditions under which observed altered fracture zones form from a Stimson unit-like parent material. Our modeling results indicate that the formation of the Stimson unit is consistent with leaching of an eolian deposit with a solution of pH = 6-8, and that formation of the altered fracture zones is consistent with leaching with a very acidic (pH = 2-3) high sulfate solution containing Ca. These results suggest circumneutral pH conditions during authigenesis or early diagenesis in the Stimson formation sediments followed by diagenetic alteration by very acidic solutions along fracture zones.

Hausrath, Elisabeth↗

DNS of Laminar-Turbulent Transition in Swept-Wing Boundary Layers

Direct numerical simulation (DNS) is performed to examine laminar to turbulent transition due to high-frequency secondary instability of stationary crossflow vortices in a subsonic swept-wing boundary layer for a realistic natural-laminar-flow airfoil configuration. The secondary instability is introduced via inflow forcing and the mode selected for forcing corresponds to the most amplified secondary instability mode that, in this case, derives a majority of its growth from energy production mechanisms associated with the wall-normal shear of the stationary basic state. An inlet boundary condition is carefully designed to allow for accurate injection of instability wave modes and minimize acoustic reflections at numerical boundaries. Nonlinear parabolized stability equation (PSE) predictions compare well with the DNS in terms of modal amplitudes and modal shape during the strongly nonlinear phase of the secondary instability mode. During the transition process, the skin friction coefficient rises rather rapidly and the wall-shear distribution shows a sawtooth pattern that is analogous to the previously documented surface flow visualizations of transition due to stationary crossflow instability. Fully turbulent features are observed in the downstream region of the flow.

Duan, L.↗

Faster-than-real-time Simulation with Demonstration for Resilient DER Integration

The US electric grid is facing operational, stability, and security challenges. Transmission system operators need some measure of visibility into distribution system renewable generation. Distribution system generation needs to support transmission system voltage. The grid is experiencing an expansion in measurement systems. How to take full advantage of this expansion and defend against attacks, both cyber and physical, poses additional challenges. The Faster-than-real-time Simulation with demonstration for Resilient DER Integration project set out to do the following: a. Flatten the voltage profile through the feeders and system for cost saving and voltage stabilization needs. b. Increase the amount of intermittent distributed energy resources (IDERs) that could be deployed on a utility feeder and provide 100% or more energy needed for the demands on that feeder, and c. based on an accurate model (Digital Twin) of the utilities system, be able to detect any abnormalities on the utilities distribution system. To manage the voltage and increase IDER penetration (a,b), Graph Trace Analysis is employed in a time-series, optimal power flow to coordinate the time-varying feedback control setpoints of a distribution feeder’s utility control devices. Under the coordinated control are a Load Tap Changing Transformer, a voltage regulator, and five switched capacitor banks. The feeder serves over 2000 customers, the feeder secondaries are modeled, and the feeder has 2.3 MW of PV generation, corresponding to a 17.4% penetration of PV generation. The feeder model has over 12,000 components, where every customer load bus and PV generator are modeled. The accuracy of the power flow solution is compared against historical meter voltage measurements, the improvement in conservation voltage reduction energy savings as a function of the coordinated control desired voltage profile is investigated, and the increase in PV penetration of the coordinated control over the existing control is presented. To achieve improved control performance while observing system operation constraints, bellwether Advanced Metering Infrastructure (AMI) voltage measurements are used to adjust the desired voltage profile used by the optimal power flow analysis. To detect and alleviate or negate attacks or failures on the distribution and transmission utility grids (c) the grid needs to be resilient and self-healing. In this project software was designed to do just that. At the center of the software is an Integrated System Model (ISM) that spans from transmission to secondary distribution. The ISM is employed in real-time abnormality detection, voltage stability forecasting, and multi-mode control. Testing results are presented for: 1—attacks on utility infrastructure; 2—energy savings from optimal control; 3—distribution system control response during a low voltage transmission system event; 4—cyber-attacks on PV inverters, where physical inverters are used in hard-ware-in-the-simulation-loop studies. Contributions of this work include real-time analysis that spans from three-phase transmission through secondary distribution; an approach for detecting abnormalities that employs measurements from three independent measurement systems; and a multi-mode distribution system control that responds to cyber-attacks, physical attacks, equipment failures, and transmission system needs.

Integrated System Model, Graph Trace Analysis, Adv↗

The effect of original and initial saturation on residual nonwetting phase capillary trapping efficiency

Injection of supercritical carbon dioxide (CO 2 ) into geological formations is a strategy for both atmospheric greenhouse gas reduction (climate change mitigation) and enhanced oil recovery. To understand CO 2 trapping efficiency, the capillary trapping behaviors that immobilize subsurface fluids were analyzed at the pore-scale using pairs of proxy fluids representing the range of in situ nonwetting and wetting fluid properties encountered in geologic storage reservoirs. The pairs of fluids were cycled through imbibition and drainage processes using a flow cell apparatus containing a sintered glass bead column. Computed x-ray microtomography (microCT) was used to identify immobilized nonwetting fluid after imbibition and drainage events. From microCT images, the trapped residual (post-secondary imbibition) nonwetting phase was spatially correlated to both the original (post-primary imbibition) and the initial (post-primary drainage) nonwetting phase; this relationship is referred to here as the original saturation dependence (SO-dependence) and initial saturation dependence (SI-dependence), respectively. Significant trends of decreasing SO- and SI-dependence with increasing wetting and nonwetting fluid phase viscosities were observed. This finding implies that the amount of CO 2 injected and ultimately trapped is dependent on the nonwetting phase (e.g. oil or gas) already present in the formation, as well as on the manner in which supercritical CO 2 is initially injected, which are factors not considered in current trapping models. To explore the potential effect that varying viscosity and interfacial tension (IFT) might have on this process, we report on a variety of fluid pairs with different viscosities and IFT.

42 ENGINEERING↗

Uranium nitride (UN) pellets with controllable microstructure and phase – fabrication by spark plasma sintering and their thermal-mechanical and oxidation properties

Dense uranium mononitride (UN) pellets with controlled microstructures and tailored grain size from large-grained to a few microns are synthesized by spark plasma sintering (SPS) combined with high energy ball milling. The impacts of the sintering conditions on fuel microstructure, grain size, physical density, and phase behavior are systematically investigated, and the thermal-mechanical properties and oxidation behavior of the SPS densified UN pellets are characterized. Higher sintering temperatures and longer ball milling durations and thus finer starting UN powders promote sintering and densification, and dense UN pellets above 95% theoretical density can be achieved by SPS at 1873 K for 10 min. UN phase purity is maintained in the SPS-densified pellets sintered at a lower temperature and short duration. A phase heterogeneity with secondary UO 2 or uranium sesquinitride (U 2 N 3 ) occurs for the UN pellets sintered at higher temperatures using finer UN powders. The hardness and fracture toughness of the SPS-densified UN pellets increase with smaller grain sizes and higher densities to 7.9 GPa and 3.5 MPa m 1/2 , respectively. Furthermore, both small (1–2 μm) and large grain-sized (30–50 µm) UN pellets exhibit good thermal conductivity. Dynamic oxidation testing by a thermogravimetric analyzer in air shows that the onset temperature for oxidation varies with microstructure and phase heterogeneity of the SPS densified UN pellets. Particularly, the smaller-grained (micron-sized) UN pellets containing uranium oxides and U 2 N 3 display lower weight gain and significantly-reduced oxidation kinetics, and full oxidation completes at a temperature above 1173 K when tested with a ramp rate of 10 K/min.

36 MATERIALS SCIENCE↗

Transmission and Distribution Real-Time Analysis Software for Monitoring and Control: Design and Simulation Testing

The US electric grid is facing operational, stability, and security challenges. Transmission system operators need some measure of visibility into distribution system renewable generation. Distribution system generation needs to support transmission system voltage. The grid is experiencing an expansion in measurement systems. How to take full advantage of this expansion and defend against attacks, both cyber and physical, poses additional challenges. This paper introduces software designed to meet these challenges. At the center of the software is an Integrated System Model (ISM) that spans from transmission to secondary distribution. The ISM is employed in real-time abnormality detection, voltage stability forecasting, and multi-mode control. The software architecture along with selected analysis modules is presented. Testing results are presented for: 1—attacks on utility infrastructure; 2—energy savings from optimal control; 3—distribution system control response during a low voltage transmission system event; 4—cyber-attacks on PV inverters, where physical inverters are used in hardware-in-the-simulation-loop studies. Contributions of this work include real-time analysis that spans from three-phase transmission through secondary distribution; an approach for detecting abnormalities that employs measurements from three independent measurement systems; and a multi-mode distribution system control that responds to cyber-attacks, physical attacks, equipment failures, and transmission system needs.

14 SOLAR ENERGY↗

Impact of phase state and non-ideal mixing on equilibration timescales of secondary organic aerosol partitioning

Abstract. Evidence has accumulated that secondary organic aerosols (SOAs) exhibit complex morphologies with multiple phases that can adopt amorphous semisolid or glassy phase states. However, experimental analysis and numerical modeling on the formation and evolution of SOA still often employ equilibrium partitioning with an ideal mixing assumption in the particle phase. Here we apply the kinetic multilayer model of gas–particle partitioning (KM-GAP) to simulate condensation of semi-volatile species into a core–shell phase-separated particle to evaluate equilibration timescales of SOA partitioning. By varying bulk diffusivity and the activity coefficient of the condensing species in the shell, we probe the complex interplay of mass transfer kinetics and the thermodynamics of partitioning. We found that the interplay of non-ideality and phase state can impact SOA partitioning kinetics significantly. The effect of non-ideality on SOA partitioning is slight for liquid particles but becomes prominent in semisolid or solid particles. If the condensing species is miscible with a low activity coefficient in the viscous shell phase, the particle can reach equilibrium with the gas phase long before the dissolution of concentration gradients in the particle bulk. For the condensation of immiscible species with a high activity coefficient in the semisolid shell, the mass concentration in the shell may become higher or overshoot its equilibrium concentration due to slow bulk diffusion through the viscous shell for excess mass to be transferred to the core phase. Equilibration timescales are shorter for the condensation of lower-volatility species into semisolid shell; as the volatility increases, re-evaporation becomes significant as desorption is faster for volatile species than bulk diffusion in a semisolid matrix, leading to an increase in equilibration timescale. We also show that the equilibration timescale is longer in an open system relative to a closed system especially for partitioning of miscible species; hence, caution should be exercised when interpreting and extrapolating closed-system chamber experimental results to atmosphere conditions. Our results provide a possible explanation for discrepancies between experimental observations of fast particle–particle mixing and predictions of long mixing timescales in viscous particles and provide useful insights into description and treatment of SOA in aerosol models.

54 ENVIRONMENTAL SCIENCES↗

Determining the extent of potential fugitive fluid migration from geologic carbon storage in hydrocarbon-bearing reservoirs: Insights from one-dimensional numerical modeling

Numerical modeling of Geologic Carbon Sequestration in permeable reservoirs initially containing hydrocarbons is conducted using the multi-phase, multi-component thermohydrologic simulator TOGA (TOUGH Oil, Gas, Aqueous; TOUGH stands for Transport Of Unsaturated Groundwater and Heat), to determine how phase and composition of the original fluids influence the extent of the zone where upward fugitive fluid migration could potentially occur, denoted R f . The area within R f comprises regions of substantially elevated pressure and free-phase CO 2 saturation, where a breach in reservoir sealing capacity would lead to upward fugitive fluid migration. The model examines the conditions within the storage reservoir that could lead to fugitive flow, but does not model the fugitive flow itself. A one-dimensional radial model of the storage reservoir is used, and three initial phase conditions are considered: single-phase aqueous, two-phase gas-aqueous, and three-phase oil-gas-aqueous. Components that may be present are H 2 O, CO 2 , CH 4 , C 4 H 10 , and C 10 H 22 . The most important factors controlling Rf are (1) the initial gas-phase saturation within the reservoir, and (2) the lateral extent of multi-phase initial conditions, particularly CO 2 . The composition of liquid and gas phases has a secondary effect. The impact of reservoir depth, thickness, injection rate, and hydrologic properties are also briefly examined, with thickness (or equivalently injection rate) having the biggest effect. These results can help to understand important trends in potential response of CO 2 -EOR fields being considered for dedicated CO 2 storage.

CO₂ plume migration↗

Investigating rapid alpha-decay induced aging of 238 PuO 2 by Raman Spectroscopy

Investigating alpha-decay induced aging in PuO 2 is useful in nuclear forensics; helping to determine the time since last calcination. This was previously investigate by monitoring the damage to 240 PuO 2 over the course of a few years. The rate of alpha-decay induced aging has long been assumed to scale only with the decay rate of other isotopes without direct verification. This article reports the first alpha-decay aging study of 238 PuO 2 by Raman spectroscopy. Contrary to the expected 10 days for 238 PuO 2 to reach steady state based on 240 PuO 2 studies, the alpha aging curve reached an approximate steady state ∼24–30 h after laser annealing. While the cause of the order of magnitude decrease is unknown, it is speculated that dynamic annealing could contribute to the differences in rate of induced aging. The Raman spectra of annealed and aged 238 PuO 2 matched the expected features in 239 PuO 2 and 240 PuO 2 spectra; showing no formation of new stable chemical species in the material such as secondary Pu oxide phases (e.g., Pu 4 O 9 ). In conclusion, results indicate that 238 Pu could be leveraged for rapid alpha-decay aging studies to characterize alpha-decay induced features and better understand matrix temperatures and the annealing of Frankel pair defects.

Actinide↗

Recovering high-purity uranyl nitrate from simulated used nuclear fuel dissolver solutions by crystallization: rejecting technetium

The separation of U from Tc and other problematic fission product elements like Mo and Ru, along with Sr, Zr, Cs, and Nd, has been achieved via the crystallization of uranyl nitrate hexahydrate (UNH). Rejection of technetium as pertechnetate anion ( 99 TcO 4 – ) is an especially important feature of this system, as it otherwise tends to follow U (VI) in extractive separations. It also raises the salient question regarding why this oxoanion cannot replace nitrate within the crystalline lattice of UNH. Results showed high-yield (>90 %), high-purity (>99 %) recovery of U as UNH from solutions containing 99 TcO 4 – by simple reduction of temperature from 60°C to 20°C. There was no observable interaction of 99 TcO 4 – with UO 2 2+ . The addition of other cations like, Sr 2+ , Zr 4+ , Cs + , and Nd 3+ , also did not form secondary, contaminant solid phases, leaving the > 99 % of the fission product elements in the mother liquor, while the U was recovered at > 90 %. Similarly, Mo and Ru, when added to the mixture, were shown to behave as the other fission-product elements, remaining in the mother liquor during crystallization. As a result, DFT calculations showed that, despite the higher binding strength of TcO 4 – , HMoO 4 – , and BiO 3 – with the UO 2 2+ cation compared to NO 3 – , the hydrogen-bonding network of the two coordinated ions and four waters of hydration in the UNH crystal structure is the driving force for the high specificity of this separation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Hydration structure and energetics of iron-exchanged montmorillonite

The interaction between iron released from corroded steel canisters and bentonite is a key process influencing the long-term performance of nuclear waste repositories. In particular, the migration of Fe²⁺ into montmorillonite (Mnt) interlayers may alter its hydration, swelling, and ion-transport properties. In the present study, molecular dynamics simulations were performed to investigate the hydration behavior, structural response, and transport properties of Fe-exchanged montmorillonite (Fe-Mnt) under varying hydration states. The simulations focus on short- to intermediate-time-scale Fe 2+ and Fe 3+ interlayer exchange and hydration effects, and do not consider long-term structural substitution, Fe-bearing clay phase stabilization, or secondary iron mineral precipitation. Systems containing Na + -, Fe 2+ -, and Fe 3+ -Mnt were examined using both periodic and edge-exposed configurations to evaluate interlayer structure, ion exchange, and free energy of Fe intercalation. The results show that Fe ions influence the interlayer spacing primarily at low water contents (<1 bilayer), where Fe-Mnt exhibits a d-spacing 1–2 Å larger than Na-Mnt due to stronger hydration. The calculated hydration energies follow the order Fe 2+ $<$Fe 3+ $<$Na + . Both water and ion diffusion coefficients decrease upon Fe ion intercalation, with Fe 2+ ions diffusing an order of magnitude more slowly than those of Na + . Free energy profiles further confirm that Fe 2+ and Fe 3+ ions are thermodynamically favored in the interlayer, with Fe 3+ being the most stable. In conclusion, these findings provide molecular-scale insights into the mechanisms of Fe–Na exchange and their implications for bentonite alteration in repository environments.

clay↗

Thermophysical modeling of raw glaze liquidus temperature and viscosity

Raw glazes, which do not contain pre-melted glass frits, are widely used due to their low cost. The natural variability of the mineral compositions used in raw glazes affects their viscosity during firing. Molten glaze viscosity plays an important role in determining the final surface quality; therefore, it is important to determine how raw glaze composition and firing conditions affect viscosity. Thermophysical modeling provides a way to evaluate the liquidus temperature and viscosity of molten glazes as a function of composition and oxygen activity. Equilibrium phase volume fraction and composition of 242 raw glaze recipes were evaluated under oxidizing and reducing conditions. Undissolved or precipitation of solid phases explain the secondary flux transition from anti-fluxing to fluxing near the liquidus temperature. The liquidus temperature and viscosity of iron glazes decrease as a function of increasing iron content under a reducing atmosphere. The empirical probability distribution of molten glaze viscosity follows a lognormal distribution, with lead glazes having significantly lower viscosity compared with glazes without lead. Here, the peak firing temperature viscosity of lead-free glazes is near the working point of glass. Multiple linear regression analysis shows that the peak firing temperature viscosity and liquidus temperature significantly predict the viscosity factor.

36 MATERIALS SCIENCE↗

Fate of Listeria monocytogenes Serotypes on Frozen Mixed Vegetables During Consumer‐Simulated Thawing and Storage

ABSTRACT Recent outbreaks and recalls associated with frozen vegetables in the United States and Europe have been linked to Listeria monocytogenes . This study aims to understand the extent to which frozen vegetables support the growth of L. monocytogenes once thawed and held at different temperatures. Six L. monocytogenes strains, two of each from serotypes 1/2a, 1/2b, and 4b, were individually inoculated onto frozen vegetables and stored at −18°C for 7 days. After 7 days, the vegetables were thawed and stored at 5°C or 10°C for up to 14 days or at 25°C for up to 7 days. L. monocytogenes was enumerated from the thawed vegetables throughout the storage period. Population data were fitted to the primary Baranyi model to estimate growth rates and lag phase durations; the secondary Ratkowsky square root model was used to model the relationship of the growth rates with storage temperature. Five of the L. monocytogenes strains survived and grew on the thawed vegetables (population increases of > 1 log CFU/g) stored at 5°C, and all six of the strains proliferated at 10°C and 25°C (population increases of > 3 log CFU/g after 14 days and > 4 log CFU/g after 7 days, respectively). A secondary model was successfully generated based on the growth rates of the six L. monocytogenes strains on the thawed vegetables ( r 2 = 0.8888, RMSE = 0.2057). Results from this study fill a data gap associated with L. monocytogenes survival on thawed vegetables and can be used to determine safe handling and storage practices for these products to protect public health.

Salazar, Joelle K. [Division of Food Processing Sc↗

Sub-pilot-scale Production of High-Value Products from U.S. Coals

Investigators from the University of Utah, University of Wyoming and Marshall University pursued a program to study the conversion of raw coal to high-value products of carbon fiber and silicon carbide. Team members also developed an initial framework for a data portal that can incorporate laboratory data on coal processing and product quality, and also work with tools for machine learning for data analysis, data visualization and economic assessment. Experimental R&D efforts focused on the conversion of raw coal to coal tar and other byproducts, and the resulting tar intermediates were upgraded to form anisotropic and isotropic pitch materials. These pitch materials were produced from coal using both thermal (pyrolysis) and chemical (mild solvolysis liquefaction) decomposition of raw coal. Four different coals were studied: Utah bituminous coal (Sufco), Wyoming PRB coal (Black Thunder), Illinois bituminous coal (Illinois #6), and West Virginia bituminous coal (Flying Eagle). Both metallurgical-grade coking coals and lower-grade steam coals were investigated, and controlled secondary gas-phase reactions were used during a two-stage pyrolysis process to induce cracking and condensation reactions among the pyrolytic tar species. This approach successfully improved the performance of the lower grade coals for yielding pitch materials, with properties more consistent with a commercial-grade pitch that had previously demonstrated success for quality carbon fiber production. The use of waste plastic materials was also studied, to help improve physical and chemical characteristics of the intermediate tars and final pitch product; in particular, for lowering the pitch softening point to an acceptable level for melt spinning carbon fiber. Mild solvolysis liquefaction was also used as a method for producing pitch for carbon fiber production. As expected, significantly higher pitch yields were obtained using this approach, and waste plastic materials were also successfully used to reduce pitch softening point to an acceptable level. The plastic materials were also utilized to create a solvent for the mild solvolysis process, and this plastic-derived solvent was shown to provide results consistent with more expensive commercial chemical solvents, and could thus avoid the need for costly recovery and recycle of a liquefaction solvent. Additional experimental R&D focused on the production of silicon carbide (β-SiC) from the residual char byproduct from pitch production, and also on the production of carbon fiber from the anisotropic pitch. SiC was successfully synthesized using a mixture of residual char and sandstone at a ratio of 1:1. Reaction temperature and residence time were optimized and yielded a product purity of 81%. For carbon fiber production, the most successful pitch samples were obtained from the mild solvolysis liquefaction approach, combined with the use of a plastic (HDPE)-derived solvent. Fiber properties improved over time as laboratory fiber production methodologies improved, and final yields of carbon fiber were obtained with a diameter of 12.14 ± 1.10 um, Modulus of 173.73 ± 15.25 GPa, and Tensile Strength of 1.04 ± 0.10 GPa. A proof-of-concept Modern Community Research Data Portal (MCRDP) was developed and deployed for coal and coal-derived pitch characterization, with the full support of (i) remote web-based access, (ii) distributed analysis, (iii) interactive visualization and exploration, (iv) shared and long-term data access, (v) advanced query capabilities and (vi) real-time collaboration. The Coal to Products Data Portal “coaltoproducts.org” provides researchers with space to store and share data within a project, tools for analyzing and understanding data for scientific investigation, and the ability to publish data to the broader community for reproducibility. The portal leverages the Material Commons 2.0 (MC) platform developed by the Center for PRedictive Integrated Structural Materials Science (PRISMS) of the University of Michigan, to achieve long-term longevity of data collections and, more importantly, collaborative science. A number of data visualization tools were also assessed and implemented for interrogating the experimental and modeling data. The machine learning portion of this project analyzed datasets from two different coal conversion processes performed on a diverse set of coal samples from both the coal pyrolysis experiments and the solvent liquefaction experiments. The work was initiated by exploring standard regression models on the pyrolysis data, aiming to understand the impact of sample characteristics and processing conditions on key product metrics. Over the course of the project, the focus expanded to include a variety of machine learning tools, delving into both supervised and unsupervised learning methods. Models tested on the pyrolysis data included linear, ridge, lasso, elastic-net, Gaussian process, random forest regression, and AutoSklearn, and the approach was continually refined to enhance predictive accuracy and model interpretability. Similar techniques were applied to the liquefaction data with an additional focus on feature engineering. Along with mesophase content, additional outputs of interest were the pitch yield, softening point, and QI content. Insights derived from these analyses are crucial in determining the factors influencing the quality and yield of coal-derived products. As the work progressed, the research evolved from foundational model comparisons to analyses of random forests, decision paths, and feature importance scores. A thorough market analysis was performed to examine the prospects of coal-based carbon fibers. The best opportunities for coal come from its lower and more stable price relative to petroleum, particularly for subbituminous coals, which is the primary advantage that a coal refinery may have over a petroleum refinery. Before a commercial CTP production facility can be modeled, however, several things need to be understood regarding the nature of the would-be coal refinery. These include the technology to be deployed, the size of facility, the volume(s) of co-product(s), and the waste and emissions profile of the plant. The volume of co-products and waste may be substantial and will require separate market analysis to ensure viability. In the near-term, the importance of coal tar pitch, in the form of carbon pitch, to the aluminum and steel industries is likely to overshadow the alternative use of this material as an input for carbon fiber. The importance of steel and aluminum in building materials, and the need for carbon materials in their manufacturing, will ensure that demand for these products remains for the long run. In addition, carbon fiber may also be the best substitute for steel and aluminum well into the future. While society will eventually be able to shift production of much of its electricity needs to renewables, it will not be able to shift away from fossil fuels for production of high-strength construction and vehicular materials. Demand for carbon fiber is expected to increase quickly, but the volume of carbon fiber and the amount of coal that would be needed to produce even a sizeable share of this market may still be relatively small compared to current coal production. Thus, other coal-based products like graphene, graphite, carbon foams, resins, and carbon-based building products will play important roles in sustaining coal production as coal-fired power generation continues to decline.

01 COAL, LIGNITE, AND PEAT↗