Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Distributed Processing”

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 271 records · Page 15

Cloud response to co-condensation of water and organic vapors over the boreal forest

Abstract. Accounting for the condensation of organic vapors along with water vapor (co-condensation) has been shown in adiabatic cloud parcel model (CPM) simulations to enhance the number of aerosol particles that activate to form cloud droplets. The boreal forest is an important source of biogenic organic vapors, but the role of these vapors in co-condensation has not been systematically investigated. In this work, the environmental conditions under which strong co-condensation-driven cloud droplet number enhancements would be expected over the boreal biome are identified. Recent measurement technology, specifically the Filter Inlet for Gases and AEROsols (FIGAERO) coupled to an iodide-adduct chemical ionization mass spectrometer (I-CIMS), is utilized to construct volatility distributions of the boreal atmospheric organics. Then, a suite of CPM simulations initialized with a comprehensive set of concurrent aerosol observations collected in the boreal forest of Finland during spring 2014 is performed. The degree to which co-condensation impacts droplet formation in the model is shown to be dependent on the initialization of temperature, relative humidity, updraft velocity, aerosol size distribution, organic vapor concentration, and the volatility distribution. The predicted median enhancements in cloud droplet number concentration (CDNC) due to accounting for the co-condensation of water and organics fall on average between 16 % and 22 %. This corresponds to activating particles 10–16 nm smaller in dry diameter that would otherwise remain as interstitial aerosol. The highest CDNC enhancements (ΔCDNC) are predicted in the presence of a nascent ultrafine aerosol mode with a geometric mean diameter of ∼ 40 nm and no clear Hoppel minimum, indicative of pristine environments with a source of ultrafine particles (e.g., via new particle formation processes). Such aerosol size distributions are observed 30 %–40 % of the time in the studied boreal forest environment in spring and fall when new particle formation frequency is the highest. To evaluate the frequencies with which such distributions are experienced by an Earth system model over the whole boreal biome, 5 years of UK Earth System Model (UKESM1) simulations are further used. The frequencies are substantially lower than those observed at the boreal forest measurement site (< 6 % of the time), and the positive values, peaking in spring, are modeled only over Fennoscandia and the western parts of Siberia. Overall, the similarities in the size distributions between observed and modeled (UKESM1) are limited, which would limit the ability of this model, or any model with a similar aerosol representation, to project the climate relevance of co-condensation over the boreal forest. For the critical aerosol size distribution regime, ΔCDNC is shown to be sensitive to the concentrations of semi-volatile and some intermediate-volatility organic compounds (SVOCs and IVOCs), especially when the overall particle surface area is low. The magnitudes of ΔCDNC remain less affected by the more volatile vapors such as formic acid and extremely low- and low-volatility organic compounds (ELVOCs and LVOCs). The reasons for this are that most volatile organic vapors condense inefficiently due to their high volatility below the cloud base, and the concentrations of LVOCs and ELVOCs are too low to gain significant concentrations of soluble mass to reduce the critical supersaturations enough for droplet activation to occur. A reduction in the critical supersaturation caused by organic condensation emerges as the main driver of the modeled ΔCDNC. The results highlight the potential significance of co-condensation in pristine boreal environments close to sources of fresh ultrafine particles. For accurate predictions of co-condensation effects on CDNC, also in larger-scale models, an accurate representation of the aerosol size distribution is critical. Further studies targeted at finding observational evidence and constraints for co-condensation in the field are encouraged.

54 ENVIRONMENTAL SCIENCES↗

Development of high-performance roll-to-roll-coated gas-diffusion-electrode-based fuel cells

This study focuses on determining fabrication conditions to create high-performance roll-to-roll-coated (R2R-coated) gas-diffusion electrodes (GDEs) for proton-exchange-membrane fuel cells (PEMFCs). Here, we examine how process conditions influence the distribution of ionomer in the electrode, which is shown to be critical for high performance. Using a combination of Kelvin probe, X-ray photoelectron spectroscopy, and nano-scale X-ray computed tomography we show that formation of an ionomer-rich surface is promoted by using a higher drying rate. We show that R2R-coated GDEs have higher surface ionomer concentration than spray-coated GDEs, which enables these R2R-coated GDEs to not need an additional ionomer overlayer, as is typically the case for spray-coated GDEs. This will reduce the number of processing steps and lower material costs in a manufacturing setting. This work shows that with the appropriate selection of materials, ink formulation, and processing conditions, direct-coated GDEs are a viable pathway for fuel cell manufacturing.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ducted Assembly Steady-State Heat Transfer Software (DASSH): Theory Manual

The Ducted Assembly Steady-State Heat transfer software (DASSH) is being developed at Argonne National Laboratory to obtain a steady-state thermal fluids result for the coolant flow and temperature distribution in a hexagonally gridded reactor core with ducted assemblies. In DASSH, each assembly is assumed to contain a hexagonal bundle of wire-wrapped pins. DASSH divides the coolant in the pin bundle into subchannels on which it applies energy balances to determine coolant and duct temperatures. DASSH offers multiple models for inter assembly heat transfer and features the ability to calculate pin temperatures. This document describes the models and equations implemented in DASSH. It includes a description of the geometric system to which the models are applied and a derivation of the subchannel energy balance formulation is presented. It also reviews the available correlations and provides an overview of the process by which power distributions are generated. Finally, it covers the code execution procedure, discussing the main code objects and their interactions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Recombination of Autodissociated Water Ions in a Nanoscale Pure Water Droplet

The recombination of water ions has diverse scientific and practical implications, ranging from acid-base chemistry and biological systems to planetary environments and applications in fuel cell and carbon conversion technologies. While spatial confinement affects the physicochemical properties of water dynamics, its impact on the recombination process has rarely been studied. In this work, we investigate the dynamics of water, the water ion distribution, and the ion recombination process in water droplets as a function of droplet size through molecular dynamics simulations and adaptive quantum mechanical/molecular mechanical calculations. We compare the dynamics of recombination in water droplet sizes ranging from 100 to 18 000 waters, both in their interiors and on their surfaces. We found that the self-diffusion of water dramatically decreases in droplets with a diameter below 2.2 nm. Using a classical RexPoN force-field, we found that the ions in 1000 H2O's spend almost 50% of the time on the surface and 0.5 nm beneath it with a slight preference for OH- ion to reside longer on the surface. We estimate that, on average, recombination in these drops occurs at 400 ps in 1000 H2O's and 1 ns in 3000 H2O's. We also found that recombination is not limited by the local structure of the surface or the size of the droplet but can be influenced by the geometry of the water wire connecting the ions as they approach each other, which can often prevent recombination. Our results provide insights to the reaction microenvironments presented by nanoscopic water droplets.

Kwon, Soonho↗

Quantifying Pore Morphology in Spray-Formed Tantalum Using X-ray Micro-computed Tomography

In order to establish quantitative process–structure–property relationships in thermal spray coatings, a robust framework for defining (micro)structural characteristics is needed. Here, we present a quantitative characterization of the three-dimensional morphology of porosity in spray-formed tantalum samples based on high-resolution X-ray micro-computed tomography. Using synchrotron facilities, we acquired dozens of high-resolution scans, enabling a statistically meaningful comparison across multiple samples, different regions within samples, and spray processes. We quantify the spatial distribution, size, and topology of porous inclusions, with a significant focus on variability across samples and different spray processes (plasma and cold sprayed), as well as sensitivity to image segmentation and resolution. Based on a typical segmentation, we report porosities ranging from 0.9 to 1.7 pct for all samples tested, with significant sensitivity due to image segmentation resulting in estimates as low as 0.6 pct and as high as 4.8 pct. For the complex pore space morphology observed in these materials, we argue that a conventional analysis based on identifying individual pores is not well suited, and propose an alternative approach based on morphological metrics with a rich history in porous media literature, such as spatial correlations, local pore thickness, and scale-dependent sub-sampling. Spatial correlations indicate anisotropic splat structures, but only mildly anisotropic pores. Various measures of pore size show a wide distribution of sizes, ranging from sub-micron to 10-micron length scales. Scale-dependent variations in porosity suggest that representative volumes of several hundred microns are required for convergence of morphological metrics, with larger volumes for cold-spray materials. This work provides a robust quantitative basis for describing three-dimensional pore structure in thermal spray coatings.

36 MATERIALS SCIENCE↗

Multi-Fidelity Modeling and Control for Building Temperature Control

The ability to control energy loads such as a building's heating, ventilation, and air conditioning (HVAC) system can help facilitate increased penetration of variable renewable energy sources into the electric grid. To be able to control these HVAC systems more effectively, detailed simulations of the corresponding building physics is becoming increasingly important. These detailed simulations can be complex, nonlinear, and can require immense computational power when used in an advanced control method such as model predictive control (MPC), prompting the need to explore less computationally intensive strategies. In this work, a multi-fidelity approach is proposed to combine samples from a complex, high-fidelity model with a simple, low-fidelity model within the MPC control loop. More specifically, the parameters of a reduced-order, linear building model are periodically updated with knowledge from its high-fidelity counterpart - an EnergyPlus model - in an online fashion using a Gaussian Process surrogate model. Hence, highly accurate predictions of current and future conditions in a building are maintained with a substantially reduced computational burden compared to using the high-fidelity models alone. In other words, this linear parameter varying model preserves the low computational requirements of a low-order linear model while accurately modeling a building's dynamics. This allows a building controller to take highly informed actions without requiring a large computational budget.

building modeling↗

Thermodynamic Implementations of Quantum Processes

Abstract Recent understanding of the thermodynamics of small-scale systems have enabled the characterization of the thermodynamic requirements of implementing quantum processes for fixed input states. Here, we extend these results to construct optimal universal implementations of a given process, that is, implementations that are accurate for any possible input state even after many independent and identically distributed (i.i.d.) repetitions of the process. We find that the optimal work cost rate of such an implementation is given by the thermodynamic capacity of the process, which is a single-letter and additive quantity defined as the maximal difference in relative entropy to the thermal state between the input and the output of the channel. Beyond being a thermodynamic analogue of the reverse Shannon theorem for quantum channels, our results introduce a new notion of quantum typicality and present a thermodynamic application of convex-split methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Physics-Informed Gaussian Process Regression for States Estimation and Forecasting in Power Grids

Real-time state estimation and forecasting are critical for the efficient operation of power grids. In this paper, a physics-informed Gaussian process regression (PhI-GPR) method is presented and used for forecasting and estimating the phase angle, angular speed, and wind mechanical power of a three-generator power grid system using sparse measurements. In standard data-driven Gaussian process regression (GPR), parameterized models for the prior statistics are fit by maximizing the marginal likelihood of observed data. In the PhI-GPR method, we propose to compute the prior statistics offline by solving stochastic differential equations (SDEs) governing the power grid dynamics. The short-term forecast of a power grid system dominated by wind generation is complicated by the stochastic nature of the wind and the resulting uncertainty in wind mechanical power. Here, we assume that the power grid dynamics are governed by swing equations, with the wind mechanical power fluctuating randomly in time. We solve these equations for the mean and covariances of the power grid states using the Monte Carlo simulation method. We demonstrate that the proposed PhI-GPR method can accurately forecast and estimate observed and unobserved states. For the considered problem, PhI-GPR has computational advantages over the ensemble Kalman filter (EnKF) method: In PhI-GPR, ensembles are computed offline and independently of the data acquisition process, whereas for EnFK, ensembles are computed online with data acquisition, rendering real-time forecast more challenging. We also demonstrate that the PhI-GPR forecast is more accurate than the EnKF forecast when the random mechanical wind power is non-Markovian. In contrast, the two methods produce similar forecasts for the Markovian mechanical wind power. For observed states, we show that PhI-GPR provides a forecast comparable to the standard data-driven GPR; both forecasts are significantly more accurate than the autoregressive integrated moving average (ARIMA) forecast. We also show that the ARIMA forecast is more sensitive to observation frequency and measurement errors than the PhI-GPR forecast.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Short-Term Load Forecasting Considering EV Charging Loads with Prediction Interval Evaluation

Short-term load forecasting plays a critical role in power system planning and operation. Along with the electrification of various loads, electricity demands are becoming increasingly hard to predict. Notably, the recent rise in electric vehicles (EVs) has further contributed to this unpredictability. To address this issue, this paper proposes a probabilistic load forecasting strategy utilizing Gaussian process regression, structured in a day-ahead manner. While many works focus on deterministic prediction, probabilistic forecasting offers additional insights into variability and uncertainty, enabling more flexible and reliable operation for power systems. To enhance the accuracy of the load forecasting model, the inputs include features related to EV charging habits as well as commonly used weather information. The load forecasting results are evaluated using various metrics, including conventional ones that assess the accuracy of point forecasts, as well as additional metrics that test the reliability of prediction intervals. The proposed load forecasting method is finally tested on real residential power consumption data and EV charging data sampled from real-world sources. The results prove that the new features can greatly improve the performance of the load forecasting method.

electrical vehicle↗

Process and feedstock driven microstructure for laser powder bed fusion of 316L stainless steel

Here in the pursuit of improving additively manufactured (AM) component quality and reliability, fine-tuning critical process parameters such as laser power and scan speed is a great first step toward limiting defect formation and optimizing the microstructure. However, the synergistic effects between these process parameters, layer thickness, and feedstock attributes (e.g. powder size distribution) on part characteristics such as microstructure, density, hardness, and surface roughness are not as well-studied. In this work, we investigate 316L stainless steel density cubes built via laser powder bed fusion (L-PBF), emphasizing the significant microstructural changes that occur due to altering the volumetric energy density (VED) via laser power, scan speed, and layer thickness changes, coupled with different starting powder size distributions. This study demonstrates that there is not one ideal process set and powder size distribution for each machine. Instead, there are several combinations or feedstock/process parameter ‘recipes’ to achieve similar goals. This study also establishes that for equivalent VEDs, changing powder size can significantly alter part density, GND density, and hardness. Through proper parameter and feedstock control, part attributes such as density, grain size, texture, dislocation density, hardness, and surface roughness can be customized, thereby creating multiple high-performance regions in the AM process space.

316L stainless steel↗

Data-Efficient Strategies for Probabilistic Voltage Envelopes under Network Contingencies

This work presents an efficient data-driven method to construct probabilistic voltage envelopes (PVE) using power flow learning in grids with network contingencies. First, a network-aware Gaussian process (GP) termed Vertex-Degree Kernel (VDK-GP), developed in prior work, is used to estimate voltage–power functions for a few network configurations. The paper introduces a novel multi-task vertex degree kernel (MT-VDK) that amalgamates the learned VDK-GPs to determine power flows for unseen networks, with a significant reduction in the computational complexity and hyperparameter requirements compared to alternate approaches. Simulations on the IEEE 30-Bus network demonstrate the retention and transfer of power flow knowledge in both N-1 and N-2 contingency scenarios. The MT-VDK-GP approach achieves over 50 % reduction in mean prediction error for novel N-1 contingency network configurations in low training data regimes (50–250 samples) over VDK-GP. Additionally, MT-VDK-GP outperforms a hyper-parameter based transfer learning approach in over 75 % of N-2 contingency network structures, even without historical N-2 outage data. Furthermore, the proposed method demonstrates the ability to achieve PVEs using sixteen times fewer power flow solutions compared to Monte-Carlo sampling-based methods.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Limitations of Bin and Bulk Microphysics in Reproducing the Observed Spatial Structure of Light Precipitation

Abstract Coarse-gridded atmospheric models often account for subgrid-scale variability by specifying probability distribution functions (PDFs) of process rate inputs such as cloud and rainwater mixing ratios ( q c and q r , respectively). PDF parameters can be obtained from numerous sources: in situ observations, ground- or space-based remote sensing, or fine-scale modeling such as large-eddy simulation (LES). LES is appealing to constrain PDFs because it generates large sample sizes, can simulate a variety of cloud regimes/case studies, and is not subject to the ambiguities of observations. However, despite the appeal of using model output for parameterization development, it has not been demonstrated that LES satisfactorily reproduces the observed spatial structure of microphysical fields. In this study, the structure of observed and modeled microphysical fields are compared by applying bifractal analysis, an approach that quantifies variability across spatial scales, to simulations of a drizzling stratocumulus field that span a range of domain sizes, drop concentrations (a proxy for mesoscale organization), and microphysics schemes (bulk and bin). Simulated q c closely matches observed estimates of bifractal parameters that measure smoothness and intermittency. There are major discrepancies between observed and simulated q r properties, though, with bulk simulated q r consistently displaying the bifractal properties of observed clouds (smooth, minimally intermittent) rather than rain while bin simulations produce q r that is appropriately intermittent but too smooth. These results suggest fundamental limitations of bulk and bin schemes to realistically represent higher-order statistics of the observed rain structure.

54 ENVIRONMENTAL SCIENCES↗

Multidimensional visualization of the dynamic evolution of Li metal via in situ/operando methods

The growing demands for high-energy density electrical energy storage devices stimulate the coupling of conversion-type cathodes and lithium (Li) metal anodes. While promising, the use of these “Li-free” cathodes brings new challenges to the Li anode interface, as Li needs to be dissolved first during cell operation. In this study, we have achieved a direct visualization and comprehensive analysis of the dynamic evolution of the Li interface. The critical metrics of the interfacial resistance, Li growth, and solid electrolyte interface (SEI) distribution during the initial dissolution/deposition processes were systematically investigated by employing multidimensional analysis methods. They include three-electrode impedance tests, in situ atomic force microscopy, scanning electrochemical microscopy, and cryogenic scanning transmission electron microscopy. The high-resolution imaging and real-time observations show that a loose, diffuse, and unevenly distributed SEI is formed during the initial dissolution process. This leads to the dramatically fast growth of Li during the subsequent deposition, deviating from Fick’s law, which exacerbates the interfacial impedance. The compactness of the interfacial structure and enrichment of electrolyte species at the surface during the initial deposition play critical roles in the long-term stability of Li anodes, as revealed by operando confocal Raman spectroscopic mapping. In conclusion, our observations relate to ion transfer, morphological and structural evolution, and Li (de)solvation at Li interfaces, revealing the underlying pathways influenced by the initial dissolution process, which promotes a reconsideration of anode investigations and effective protection strategies.

25 ENERGY STORAGE↗

The ezSpectra suite: An easy‐to‐use toolkit for spectroscopy modeling

Abstract A molecule's spectrum encodes information about its structure and electronic properties. It is a unique fingerprint that can serve as a molecular ID. Quantum chemistry calculations provide key ingredients for interpreting spectra, but modeling the spectra rarely ends there; it requires additional steps that entail combined treatments of electronic and nuclear degrees of freedom and account for specifics of the experimental setup (light energy, polarization, averaging over molecular orientations, temperature, etc.). This Software Focus article describes the ezSpectra suite, which currently comprises two stand‐alone open‐source codes: ezFCF and ezDyson . ezFCF calculates Franck–Condon factors, which yield vibrational progressions for polyatomic molecules, within the double‐harmonic approximation. ezDyson calculates absolute cross‐sections for photodetachment/photoionization processes and photoelectron angular distributions using Dyson orbitals computed by a quantum chemistry program. This article is categorized under: Electronic Structure Theory > Ab Initio Electronic Structure Methods Theoretical and Physical Chemistry > Spectroscopy Software > Simulation Methods

Gozem, Samer↗

Dark Radiation Constraints on Heavy QCD Axions

The naturalness problem of PQ symmetry motivates study of the heavy QCD axion, with masses m a > 1 MeV generated at scales above the QCD scale, and low values of the PQ symmetry breaking scale, ƒ a . We compute the abundance of such axions in a model-independent way, assuming only that they freeze-out after reheating from inflation, and are not subsequently diluted by new physics. If these axions decay between neutrino decoupling and the last scatter era of the Cosmic Microwave Background (CMB), they dilute the neutrinos and their abundance is constrained by CMB measurements of the energy density in dark radiation, N eff . We accurately compute this bound using a numerical code to evolve the axion momentum distribution, including many key processes and effects previously ignored. We assume that the only relevant axion decays are to final states involving Standard Model particles. We determine regions of (m a , ƒ a ) that will give a signal in N eff at CMB Stage 4 experiments. We similarly compute the N eff bound and CMB Stage 4 signal for heavy axions that can decay to light mirror photons. Finally, we compute the bounds on heavy axions with mass below 1 MeV that decay after the era of CMB last scatter, from their contribution to cold or hot dark matter or N eff at this era.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Zonal-Based Emission Source Term Model for Predicting Particulate Emission Factors in Wildfire Simulations

A physics/chemistry-based numerical model for predicting the emission of fine particles from wildfires is proposed. This model implements the fundamental mechanisms of soot formation in a combustion environment: soot nucleation, surface growth, agglomeration, oxidation, and particle fragmentation. These mechanisms occur on a scale too fine for the discretization of most wildfire models, which need to simulate landscape-scale dynamics. As a result this model implements a zonal approach, where the computed soot particle distribution is partitioned into process zones within a single resolved grid cell. These process zones include: an inception zone (for nucleation), a heating zone (for coagulation, surface growth, and fragmentation), a reaction zone (for oxidation), and a quenched zone (for atmospheric processes). Governing mechanisms are applied to the appropriate zones to predict total particle growth and emission. Additionally, the proposed model is implemented into HIGRAD/FIRETEC, a physics-based wildfire simulation code which couples interactions between fire, fuels, atmosphere, and topography on a landscape scale. Fire simulations among grasslands and conifer forests are performed and compared against experimental data for emission factors.

54 ENVIRONMENTAL SCIENCES↗

Control methods for mitigating flow oscillations in a supercritical CO 2 recompression closed Brayton cycle

A dynamic model of a 10 MWe supercritical CO 2 (sCO 2 ) recompression closed Brayton power generation cycle is used to investigate control methods for mitigating oscillations in flow conditions and power demand load, especially as the load ramp rate is increased and the load setpoint is more closely tracked. The focus is on control of the main compressor inlet temperature (MCIT) since oscillations are prompted by the strong nonlinearities in sCO 2 fluid properties near the critical point and by changes in system inventory for efficient load management. In this study, the oscillatory behavior during load ramping is first shown and then two potential control solutions for substantially reducing oscillations are given: 1) the use of the main compressor inlet guide vanes (IGV) and 2) bypassing a portion of the sCO 2 around the cooler to maintain the MCIT. These control methods reduce the effect of cycle feedback from interactions between the MCIT and load. When implementing either control method, a load setpoint ramp of 7.5%/min was maintainable while achieving closer load setpoint tracking with limited oscillations. Finally, the impact of sensor noise on the control was examined since derivative action was used in the IGV and cooler sCO 2 bypass control. Finally, while some effect was evident, signal noise was not problematic to the control.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Distribution of rare earth elements in the pilot-scale processing of fly ashes derived from eastern Kentucky coals: Comparisons of the feed and processed ashes

Beneficiated fly ash from the combustion of Central Appalachian high volatile bituminous coals was extracted with HNO 3 in a pilot-scale processing plant. Several major oxides (notably CaO and SO 3 , but also including Fe 2 O 3 , MgO, K 2 O, and P 2 O 5 ) and minor elements (Mn, As, Sr, Ba, and Pb) are depleted in the post-HNO 3 -extraction spent ash. The total lanthanides, Y, and Sc concentration is reduced by about 20% in the spent ash, with Gd showing the greatest decrease. Along with Gd, Nd and Dy are also well differentiated between the feed and spent ashes, with La and Sm showing minimal partitioning. The Gd decrease is correlated with the depletion of Fe 2 O 3 . The heavy rare earth elements (REE heavier than Eu) and Y are disproportionately concentrated in the HNO 3 -leachate compared to the light REE. For the ashes studied, Sc did not partition between the feed and spent ashes. Pozzolanicity tests show that the compressive strength and strength activity indices of the spent ash + ordinary Portland cement (OPC) mixes are comparable to 100% OPC, indicating that the spent ashes produced in the pilot-scale runs have the potential to be sold as a Class F fly ash. Finally, the beneficiated ash chemistry influences the chemistry of the post-HNO 3 -extraction spent ash and the HNO 3 -leachate. A 500-ppm-REE fly ash will presumably be a more economically favorable feedstock than an ash with a significantly lesser concentration.

01 COAL, LIGNITE, AND PEAT↗