Search NASASearch

SEARCH · Search NASA

Results for “process optimization”

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 181 records · Page 10

Cost of Deconstruction Depots for Diversified, Waste-Based Lignocellulosic Sugars Using Distillable Solvents

Transitioning to a bioeconomy that makes use of low-emission and waste feedstocks requires greater flexibility to accommodate seasonal variations and mitigate long-term storage challenges, such as material loss and fire risk. To achieve this goal, biomass deconstruction technologies must efficiently handle diverse feedstocks. Here, we assess the cost of using butylaminea distillable solventto deconstruct 22 different biomass feedstocks: 7 herbaceous, 9 woody, 4 food processing residues, and 2 blends. Lignocellulosic sugar production costs, based on current empirical data, range from $1.3 to 6.1/kg, suggesting that substantial improvements are required to compete with conventional sugars. The high solvent loading (850 g/kg of whole slurry) is a process bottleneck. Lowering the solvent loading to 59 g/kg of whole slurry, demonstrated in an L-scale reactor using poplar biomass, reduces the minimum sugar selling price by 33%. Solvent loading and recovery, solid loading, sugar yield, enzyme use, and delivered biomass cost all play key roles in reaching sugar production costs of $0.45-0.79/kg. Strategic feedstock blending to maximize carbohydrate content, process optimization to improve conversion efficiency, and the selection of low-cost feedstocks are important to advancing feedstock-flexible biorefineries.

TEA

Technoeconomic and Life Cycle Analysis of an Integrated Fermentation and Microbial Electrochemical Process for Volatile Fatty Acid Production from Food Waste

Techno-economic analysis (TEA) and life cycle assessment (LCA) were conducted for an integrated system designed for the production of volatile fatty acid (VFA) from food waste. The TEA estimated a production cost of $\$$3.12/kg VFA, and the LCA predicted negative greenhouse gas (GHG) emissions of -0.4 kg CO 2 e/kg VFA, driven primarily by diverting organic waste from landfills and avoiding methane emissions while producing valuable chemical products. Hotspot analysis showed arrested methanogenesis (AM) fermentation as the largest contributor to costs (37%) and environmental burden (47%), driven by high sodium hydroxide (NaOH) consumption. Distillation and microbial electrosynthesis (MES) units were the next-largest environmental contributors (28% and 18%). Major cost drivers also included residuals management (biosolids and wastewater) and the equipment and operating costs for AM, MES, and sonication pretreatment units. Although the new integrated system is environmentally benign, its costs and environmental impacts can be further reduced by integrating alternative energy sources, minimizing chemical and energy inputs through process optimization, and improving efficiency. In conclusion, this work highlighted the viability of waste-derived VFA production and provided a clear, data-driven strategy to accelerate the commercialization of waste valorization technology.

Carboxylic Acid Production

Visualizing Size-Dependent Dynamics of CeO 2-δ {100}-Supported CoO x Nanoparticles Under CO 2 Hydrogenation Conditions

Carbon dioxide is a major greenhouse gas. In order to optimize processes focused on its chemical valorization, one needs detailed information about the effects of CO 2 and/or CO 2 /H 2 mixtures on the structure and morphology of metal/oxide catalysts. Here, in this study, the evolution of a catalyst with cobalt supported on CeO 2 -cube nanostructures under CO 2 hydrogenation conditions was investigated by using a set of in situ characterization techniques (X-ray absorption fine structure, X-ray diffraction, diffuse reflectance infrared Fourier transform spectroscopy, and environmental transmission electron microscopy (TEM)). The {100} facets of the ceria support displayed an unexpectedly high stability due to strong interactions with the aggregates of cobalt oxide. A significant influence of interfacial bonding between CoO x and CeO 2-δ {100} is evident through a clear preference in the orientation of CoO x nanoparticles (NPs) with respect to the substrate. For initially reduced Co/CeO 2 -cube nanostructures, a kinetically controlled oxidation of cobalt upon the introduction of CO 2 was observed during the early stages of CO 2 hydrogenation. Environmental TEM revealed the size-dependent morphological behavior of cobalt oxide NPs due to strong interactions with the CeO 2 {100} surface. When the environment was switched from H 2 to a mixture of H 2 and CO 2 at 250 °C, small CoO x NPs (in the largest dimension < 2.5 nm) rapidly transformed from a pyramidal three-dimensional (3D) form to a planar, monatomic layer attached to the concurrently oxidized CeO 2-δ {100} surface. This maximizes the number of sites available for the binding of CO 2 or reaction intermediates. The shape transformation reflected the oxophilic character of cobalt and strong metal–support interactions. The removal of CO 2 from the gas phase led to a reduction of the cobalt oxide NPs by hydrogen and a reversible two-dimensional → 3D transformation. In contrast, no significant morphological changes, apart from further oxidation, were observed for big CoO x NPs (in the largest dimension > 3 nm). These trends are not seen for nanoparticles of noble metals. The observed morphological and structural changes in the small CoO x NPs affected the stability of reaction intermediates and modified the selectivity of the CoO x /CeO 2 catalyst system for methane production.

36 MATERIALS SCIENCE

EF-Hand Battle Royale: Hetero-ion Complexation in Lanmodulin

The lanmodulin (LanM) protein has emerged as an effective means for rare earth element (REE) extraction and separation from complex feedstocks without the use of organic solvents. Whereas the binding of LanM to individual REEs has been well characterized, little is known about the thermodynamics of mixed metal binding complexes (i.e., heterogeneous ion complexes), which limits the ability to accurately predict separation performance for a given metal ion mixture. In this paper, we employ the law of mass action to establish a theory of perfect cooperativity for LanM-REE complexation at the two highest-affinity binding sites. The theory is then used to derive an equation that explains the nonintuitive REE binding behavior of LanM, where separation factors for binary pairs of ions vary widely based on the ratio of ions in the aqueous phase, a phenomenon that is distinct from single-ion-binding chemical chelators. We then experimentally validate this theory and perform the first quantitative characterization of LanM complexation with heterogeneous ion pairs using resin-immobilized LanM. Importantly, the resulting homogeneous and heterogeneous constants enable accurate prediction of the equilibrium state of LanM in the presence of mixtures of up to 10 REEs, confirming that the perfect cooperativity model is an accurate mechanistic description of REE complexation by LanM. We further employ the model to simulate separation performance over a range of homogeneous and heterogeneous binding constants, revealing important insights into how mixed binding differentially impacts REE separations based on the relative positioning of the ion pairs within the lanthanide series. In addition to informing REE separation process optimization, these results provide mathematical and experimental insight into competition dynamics in other ubiquitous and medically relevant, cooperative binding proteins, such as calmodulin.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

A holistic platform for accelerating sorbent-based carbon capture

Abstract Reducing carbon dioxide (CO 2 ) emissions urgently requires the large-scale deployment of carbon-capture technologies. These technologies must separate CO 2 from various sources and deliver it to different sinks 1,2 . The quest for optimal solutions for specific source–sink pairs is a complex, multi-objective challenge involving multiple stakeholders and depends on social, economic and regional contexts. Currently, research follows a sequential approach: chemists focus on materials design 3 and engineers on optimizing processes 4,5 , which are then operated at a scale that impacts the economy and the environment. Assessing these impacts, such as the greenhouse gas emissions over the plant’s lifetime, is typically one of the final steps 6 . Here we introduce the PrISMa (Process-Informed design of tailor-made Sorbent Materials) platform, which integrates materials, process design, techno-economics and life-cycle assessment. We compare more than 60 case studies capturing CO 2 from various sources in 5 global regions using different technologies. The platform simultaneously informs various stakeholders about the cost-effectiveness of technologies, process configurations and locations, reveals the molecular characteristics of the top-performing sorbents, and provides insights on environmental impacts, co-benefits and trade-offs. By uniting stakeholders at an early research stage, PrISMa accelerates carbon-capture technology development during this critical period as we aim for a net-zero world.

Science & Technology - Other Topics

Acid–base concentration swing for direct air capture of carbon dioxide

This work demonstrates the first experimental evidence of the acid–base concentration swing (ABCS) for direct air capture of CO 2 . This process is based on the effect that concentrating particular acid–base chemical reactants will strongly acidify solution, through Le Chatelier's principle, and result in outgassing absorbed CO 2 . After collecting the outgassed CO 2 , diluting the solution will result in a reversal of the acid–base reaction, basifying the solution and allowing for atmospheric CO 2 absorption. The experimental study examines a system that includes sodium cation as the alkalinity carrier, boric acid, and a polyol complexing agent that reversibly reacts with boric acid to strongly acidify solution upon concentration. Though the tested experimental system faces absorption rate and water capacity limitations, the ABCS process described here provides a basis for further process optimization. A generalized theoretical ABCS reaction framework is developed and different reaction orders and conditions are studied mathematically. Higher order reactions yield favorable cycle output results, reaching volumetric cycle capacity above 50 mM for third-order and 80 mM for fourth-order reactions. Optimal equilibrium constants are determined in order to guide alternative chemical searches and synthetic chemistry design targets. There is a substantial energetic benefit for reaction orders above the first, with second- and third-order ABCS cycles exhibiting a thermodynamic minimum work for the concentrating and outgassing steps around 150 kJ per mole of CO 2 . A significant advantage of the ABCS is that it can be driven through well-developed and widely-deployed desalination technologies, such as reverse osmosis, with opportunities for energy recovery when recombining the concentrated and diluted streams, and extraction can occur directly from the liquid phase upon vacuum application.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Demystifying group-4 polyolefin hydrogenolysis catalysis: Gaseous propane hydrogenolysis mechanism over the same catalysts

A kinetic/mechanistic investigation of gaseous propane hydrogenolysis over the single-site heterogeneous polyolefin depolymerization catalysts AlS/ZrNp 2 and AlS/HfNp 2 (AlS = sulfated alumina, Np = neopentyl), is use to probe intrinsic catalyst properties without the complexities introduced by time- and viscosity-dependent polymer medium effects. In a polymer-free automated plug-flow catalytic reactor, propane hydrogenolysis turnover frequencies approach 3,000 h −1 at 150 °C. Both catalysts exhibit approximately linear relationships between rate and [H 2 ] at substoichiometric [H 2 ] with rate law orders of 0.66 ± 0.09 and 0.48 ± 0.07 for Hf and Zr, respectively; at higher [H 2 ], the rates approach zero-order in [H 2 ]. Reaction orders in [C 3 H 8 ] and [catalyst] are essentially zero-order under all conditions, with the former implying rapid, irreversible alkane binding/activation. This rate law, activation parameter, and DFT energy span analysis support a scenario in which [H 2 ] is pivotal in one of two plausible and competing rate-determining transition states—bimolecular metal-alkyl bond hydrogenolysis vs. unimolecular β-alkyl elimination. The Zr and Hf catalyst activation parameters, ΔH ‡ = 16.8 ± 0.2 kcal mol −1 and 18.2 ± 0.6 kcal mol −1 , respectively, track the relative turnover frequencies, while ΔS ‡ = −19.1 ± 0.8 and −16.7 ± 1.4 cal mol −1 K −1 , respectively, imply highly organized transition states. These catalysts maintain activity up to 200 °C, while time-on-stream data indicate multiday activities with an extrapolated turnover number ~92,000 at 150 °C for the Zr catalyst. This methodology is attractive for depolymerization catalyst discovery and process optimization.

03 NATURAL GAS

AI-assisted detector design for the EIC (AID(2)E)

Artificial Intelligence is poised to transform the design of complex, large-scale detectors like ePIC at the future Electron Ion Collider. Featuring a central detector with additional detecting systems in the far forward and far backward regions, the ePIC experiment incorporates numerous design parameters and objectives, including performance, physics reach, and cost, constrained by mechanical and geometric limits. This project aims to develop a scalable, distributed AI-assisted detector design for the EIC (AID(2)E), employing state-of-the-art multiobjective optimization to tackle complex designs. Supported by the ePIC software stack and using G EANT 4 simulations, our approach benefits from transparent parameterization and advanced AI features. The workflow leverages the PanDA and iDDS systems, used in major experiments such as ATLAS at CERN LHC, the Rubin Observatory, and sPHENIX at RHIC, to manage the compute intensive demands of ePIC detector simulations. Tailored enhancements to the PanDA system focus on usability, scalability, automation, and monitoring. Ultimately, this project aims to establish a robust design capability, apply a distributed AI-assisted workflow to the ePIC detector, and extend its applications to the design of the second detector (Detector-2) in the EIC, as well as to calibration and alignment tasks. Additionally, we are developing advanced data science tools to efficiently navigate the complex, multidimensional trade-offs identified through this optimization process.

97 MATHEMATICS AND COMPUTING

Microscopic Imprints of Learned Solutions in Tunable Networks

In physical networks trained using supervised learning, physical parameters are adjusted to produce desired responses to inputs. An example is an electrical contrastive local learning network of nodes connected by edges that adjust their conductances during training. When an edge conductance changes, it upsets the current balance of every node. In response, physics adjusts the node voltages to minimize the dissipated power. Learning in these systems is therefore a coupled double-optimization process, in which the network descends both a cost landscape in the high-dimensional space of edge conductances and a physical landscape—the power dissipation—in the high-dimensional space of node voltages. Because of this coupling, the physical landscape of a trained network contains information about the learned task. Here, we derive a structure-function relation for trained tunable networks and demonstrate that all the physical information relevant to the trained input-output relation can be captured by a tuning susceptibility, an experimentally measurable quantity. We supplement our theoretical results with simulations to show that the tuning susceptibility is correlated with functional importance and that we can extract physical insight into how the system performs the task from the conductances of highly susceptible edges. Our analysis is general and can be applied directly to mechanical networks, such as networks trained for protein-inspired function such as allostery.

36 MATERIALS SCIENCE

Optimization Methodology of Polar Direct-Drive Illumination for the National Ignition Facility

Improved laser illumination uniformity drives shocks and implosions to create more extreme high energy density environments. Predominantly, the geometry of experiments that can be performed is dictated by the layout of beams at laser facilities, limiting inter-facility and multiscale investigations. This Letter presents the first automated, algorithmic approach for generating illumination configurations for high energy density experiments. The method is demonstrated in comparison to a polar direct drive solid target experiment at the National Ignition Facility. The new illumination configuration is simulated to create greater than ×3 higher peak pressure and almost ×2 higher density by maintaining better shock uniformity. Furthermore, the optimization process is performed with reduced computational expense and isotropic plasma profiles, while accounting for the impact of cross-beam energy transfer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Bridging length scales in hard materials with ultra-small angle X-ray scattering – a critical review

Owing to their exceptional properties, hard materials such as advanced ceramics, metals and composites have enormous economic and societal value, with applications across numerous industries. Understanding their microstructural characteristics is crucial for enhancing their performance, materials development and unleashing their potential for future innovative applications. However, their microstructures are unambiguously hierarchical and typically span several length scales, from sub-ångstrom to micrometres, posing demanding challenges for their characterization, especially for in situ characterization which is critical to understanding the kinetic processes controlling microstructure formation. This review provides a comprehensive description of the rapidly developing technique of ultra-small angle X-ray scattering (USAXS), a nondestructive method for probing the nano-to-micrometre scale features of hard materials. USAXS and its complementary techniques, when developed for and applied to hard materials, offer valuable insights into their porosity, grain size, phase composition and inhomogeneities. We discuss the fundamental principles, instrumentation, advantages, challenges and global status of USAXS for hard materials. Using selected examples, we demonstrate the potential of this technique for unveiling the microstructural characteristics of hard materials and its relevance to advanced materials development and manufacturing process optimization. We also provide our perspective on the opportunities and challenges for the continued development of USAXS, including multimodal characterization, coherent scattering, time-resolved studies, machine learning and autonomous experiments. Our goal is to stimulate further implementation and exploration of USAXS techniques and inspire their broader adoption across various domains of hard materials science, thereby driving the field toward discoveries and further developments.

36 MATERIALS SCIENCE

Seasonal Reconfiguration of Electrical Distribution Systems to Mitigate the Impact of Electric Vehicle Charging

Power grids face challenges in their infrastructure related to the integration of electric vehicles (EV). In particular, EV charging stations may induce instability in key system parameters such as substantial voltage drops, active power losses, and transformer overload due to high demand during charging periods. This article presents a seasonal reconfiguration strategy based on the differential evolution (DE) algorithm, aimed at enhancing system performance under highly variable and stochastic load profiles, particularly those driven by EV charging. The DEA algorithm is hybridized with the find-union (FU) algorithm to efficiently ensure network radiality throughout the optimization process. The proposed methodology is validated on a hybrid distribution system composed of the IEEE 33-bus network, a modified IEEE 13-bus system, and a specific 13-bus microgrid. Results have demonstrated that seasonal reconfiguration significantly reduces active power losses and mitigates transformer loading during critical demand hours, thereby quantifiably increasing the system’s performance. As an integral component of the proposed approach, an analysis of CO2 emissions associated with energy losses is included, allowing a contextualized assessment of the environmental benefits of seasonal reconfiguration in various geographical areas.

24 POWER TRANSMISSION AND DISTRIBUTION

Goated: goal-oriented tensor decompositions in python

SAND2026-20464O Goated performs goal-oriented tensor decompositions in Python, enabling efficient compression of multi-dimensional simulation data. It extends common tensor decomposition methods by incorporating domain-specific knowledge, such as conservation laws in physics, through a penalty term in the optimization process. This approach improves data compression and modeling accuracy across various applications, including physics simulations, by using specialized algorithms and structure-aware subroutines to accelerate solver performance. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-NA0003525.

SciDAC

A miniaturized feedstocks-to-fuels pipeline for screening the efficiency of deconstruction and microbial conversion of lignocellulosic biomass

Sustainably grown biomass is a promising alternative to produce fuels and chemicals and reduce the dependency on fossil energy sources. However, the efficient conversion of lignocellulosic biomass into biofuels and bioproducts often requires extensive testing of components and reaction conditions used in the pretreatment, saccharification, and bioconversion steps. This restriction can result in a significant and unwieldy number of combinations of biomass types, solvents, microbial strains, and operational parameters that need to be characterized, turning these efforts into a daunting and time-consuming task. Here we developed a high-throughput feedstocks-to-fuels screening platform to address these challenges. The result is a miniaturized semi-automated platform that leverages the capabilities of a solid handling robot, a liquid handling robot, analytical instruments, and a centralized data repository, adapted to operate as an ionic-liquid-based biomass conversion pipeline. The pipeline was tested by using sorghum as feedstock, the biocompatible ionic liquid cholinium phosphate as pretreatment solvent, a “one-pot” process configuration that does not require ionic liquid removal after pretreatment, and an engineered strain of the yeast Rhodosporidium toruloides that produces the jet-fuel precursor bisabolene as a conversion microbe. By the simultaneous processing of 48 samples, we show that this configuration and reaction conditions result in sugar yields (~70%) and bisabolene titers (~1500 mg/L) that are comparable to the efficiencies observed at larger scales but require only a fraction of the time. We expect that this Feedstocks-to-Fuels pipeline will become an effective tool to screen thousands of bioenergy crop and feedstock samples and assist process optimization efforts and the development of predictive deconstruction approaches.

09 BIOMASS FUELS

TRANSFER LEARNING FOR FIELD EMISSION MITIGATION IN CEBAF SRF CAVITIES

The Continuous Electron Beam Accelerator Facility (CEBAF) at Jefferson Lab operates hundreds of super-conducting radio frequency (SRF) cavities in its two linear accelerators (linacs). Field emission (FE) is an ongoing operational challenge in higher gradient SRF cavities. FE generates high levels of neutron and gamma radiation leading to damaged accelerator hardware and a radiation hazard environment. During machine development periods, we performed gradient scans to record data capturing the relationship between cavity gradients and radiation levels measured throughout the linacs. However, the field emission environment at CEBAF varies considerably over time as the configuration of the radio frequency (RF) gradients changes and due to the changing behaviour of field emitters. An artificial intelligence/machine learning (AI/ML) approach with transfer learning could be a valuable tool to mitigate FE and lower the radiation levels. In this work, we mainly focus on leveraging the RF trip data gathered during CEBAF operations. We develop a transfer learning-based surrogate model for radiation detector readings given RF cavity gradients to track the CEBAF?s changing configuration and environment. Then, we could use the developed model as an optimization process for redistributing the RF gradients within a linac to minimize radiation levels.

Ahammed, K.

TRANSFER LEARNING FOR FIELD EMISSION MITIGATION IN CEBAF SRF CAVITIES

The Continuous Electron Beam Accelerator Facility (CEBAF) at Jefferson Lab operates hundreds of super-conducting radio frequency (SRF) cavities in its two linear accelerators (linacs). Field emission (FE) is an ongoing operational challenge in higher gradient SRF cavities. FE generates high levels of neutron and gamma radiation leading to damaged accelerator hardware and a radiation hazard environment. During machine development periods, we performed gradient scans to record data capturing the relationship between cavity gradients and radiation levels measured throughout the linacs. However, the field emission environment at CEBAF varies considerably over time as the configuration of the radio frequency (RF) gradients changes and due to the changing behaviour of field emitters. An artificial intelligence/machine learning (AI/ML) approach with transfer learning could be a valuable tool to mitigate FE and lower the radiation levels. In this work, we mainly focus on leveraging the RF trip data gathered during CEBAF operations. We develop a transfer learning-based surrogate model for radiation detector readings given RF cavity gradients to track the CEBAF?s changing configuration and environment. Then, we could use the developed model as an optimization process for redistributing the RF gradients within a linac to minimize radiation levels.

Ahammed, K.

Preliminary Look at the LBEG & MPEG Beam Transmissions between HPSim and Operation Data

This report summarizes recent work on estimating the transmission of LANSCE H- beams from the end of the present DTL to both the PSR stripper foil (LBEG) and WNR target 4 (MPEG). The LBEG beam might be considered a more typical LINAC beam and lends itself to continuous monitoring of transmission through the various stages of the accelerator. For the MPEG, however, the widely-spaced micropulses and lack of the requisite sensitivity current monitors throughout the accelerator make these measurements extremely difficult and more uncertain. Therefore, to assist in estimating the MPEG beam transmission, beam-dynamics simulations using HPSim were employed. These beam transmission estimates from the end of the Drift Tube Linac (DTL) to Target 4 (MPEG) and the PSR stripper foil (LBEG) are vital to determine the charge requirements for LAMP’s front end. In this technote, we demonstrate three major efforts in determining the transmission: (1) Convert the WNR beamline lattice from TRANSPORT to HPSim for use in the simulation; (2) Simulate the optimization process in the Central Control Room (CCR) that brings down the losses by up to a factor of 4000 between the end of the initial physics-based tuneup phase and production beam operation; (3) Analyze operational data to deduce measured transmissions. Table 1 shows the estimated losses with HPSim and operational analysis. Finally, a better measurement and other improvements to refine the results are also proposed.

43 PARTICLE ACCELERATORS

Integration of LIBS with Machine Learning for Real-Time Monitoring of Feedstock in H 2 Gasification Applications

This project, funded by the U.S. Department of Energy (DOE) – Office of Fossil Energy under Award Number DE-FE0032177, aimed to assess the feasibility of an integrated Laser-Induced Breakdown Spectroscopy (LIBS) system with advanced machine learning (ML) models for real-time characterization and potential control of hydrogen gasifiers running on waste materials as feedstocks. This was a multidisciplinary effort that encompassed the acquisition and standardized analysis of individual and blended feedstocks—comprising biomass, coal waste, and plastic waste, followed by the development of a dynamic LIBS bench system for material sample analysis and development of predictive ML models. Comprehensive laboratory testing enabled the creation of a robust elemental dataset that served as the foundation for ML model training. Techniques such as Random Forest, Gradient Boosting, Support Vector Regression, and Neural Networks were employed to predict key feedstock properties, including higher heating value (HHV), moisture content, thermal conductivity, and ash composition with high accuracy. The results were validated against experimental data and demonstrated strong potential for real-time application in gasifier control systems. The project concluded with a study on the integration of the LIBS+ML approach for gasifier control and a techno-economic analysis of the implementation of the approach into hydrogen (H 2 ) gasification systems. Dissemination of results was carried out at a DOE meeting. This work establishes a scalable framework for automated, in-line feedstock quality assessment, offering significant implications for process optimization and emissions reduction in hydrogen production.

01 COAL, LIGNITE, AND PEAT