Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Clay”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

The evolution of soil microstructure and micromineralogy along a soil age gradient on the Galápagos Islands (Ecuador)

The islands of the Galápagos archipelago are of different ages due to the combination of a geostationary volcanic hotspot and the eastward movement of the Nazca tectonic plate. Taking advantage of this, a soil chronosequence (1.5–1070 ka) has recently been established covering four of the islands. This was used to investigate the temporal development of soil microstructure and micromineralogy using micromorphological methods in combination with synchrotron-based mineralogical and geochemical microanalyses. Along the chronosequence, the microstructure changed from intergrain microaggregates in the young soils to angular blocky microaggregates in the older soils. In the young soils, a high amount of vitric scoria and weatherable minerals occur, but they quickly disappear with age, and micromass increased fast. The micromorphological studies as well as micro-XRD and micro-XRF revealed increasing heterogeneity of soil microstructure and micromineralogy due to weathering and pedogenic processes until intermediate soil age. In the oldest soils, the microstructure became more uniform again, de-mixing pedofeatures like clay coatings have completely disappeared, and kaolinite, hematite and gibbsite dominate the mineralogical composition. Our study shows that the development of soil microstructure mirrors macroscopic soil profile development showing increasing differentiation until intermediate developmental stages and receding heterogeneity in highly weathered soils.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Implications of rootless geothermal models: Missing processes, parameter compensation, and imposter convection

Numerical models of geothermal systems commonly capture only the top of a reservoir. Deeper areas of the reservoir are simplified to a boundary at the base of the model domain. Commonly, the basal boundary is given either a heat source or a source of mass and enthalpy. Here we developed and present simple numerical experiments which demonstrate that these approaches do not produce the correct model behavior in comparison to a model that captures the entire convecting domain with a heat flux only. We describe a variety of incorrect types of model behavior that arise directly from the choice of boundary condition, independent of the specific parameterization of the model. Heat sources are sensitive to the thickness of the domain and parameters take unphysical values to compensate for the reduced height. The combined mass/heat boundary can produce temperatures that show similarities to circulating geothermal systems, but with incorrect fluid flow and a strong boundary layer focused at the base of the simulated clay cap. These errors likely cause parameters to adjust their values to compensate for the incorrect physics. We highlight these issues and show an example from a developed reservoir model. Initial calibrations to natural state temperature were unsuitable for history matching. The 3D model required parameter adjustments to achieve a more realistic production model. Appropriate mitigation measures should be considered to reduce parameter compensation and improve decision-support models.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Stability study of aqueous foams under high-temperature and high-pressure conditions relevant to Enhanced Geothermal Systems (EGS)

In Enhanced Geothermal System (EGS) hydraulic fracturing is carried out by injecting cold water into deep Hot Dry Rocks (HDR) under carefully controlled conditions to create new or reopen existing fractures. Water-based fracturing fluids demonstrate some challenges including immense quantity of water usage, water sensitivity of the formations, water blocking, and lack of proppant carrying capacity and transportation. Thus, an alternative is to use foam-based fracturing fluid which offers potential advantage over conventional water-based fracturing fluid such as minimum water usage, reduced wellbore damage, high proppant carrying capacity, and less environmental damage. However, foams are complex mixture of gaseous phase and liquid phase which are thermodynamically unstable at downhole conditions, and their stability can decrease over time due to liquid drainage, bubble coarsening, and coalescence. This paper shows laboratory experiments executed to study foam stability at high temperature (200°C) and high pressure (6.9MPa) conditions which simulates the geothermal environment. Foam stability was characterized by half-life of foam, which is defined as the time taken by the foam to decreases by 50% of its original height due to drainage. In this paper, two types of gaseous phases, nitrogen (N 2 ) and carbon dioxide (CO 2 ) were investigated. Also, based on successful practice of foam-based fracturing fluid in oil and gas industries, four surfactants, including Alpha olefin sulfonate (AOS), Sodium dodecyl sulfonate (SDS), Tergitol™ (NP – 40), and Cetyltrimethylammonium chloride (CTAC) at optimum concentration of 1wt.% were tested for best stability performance. In addition, different stabilizing agents including guar gum, bentonite clay, crosslinker, silicon dioxide nanoparticles (SiO 2 ), graphene oxide (GO) were also studied. Experimental results showed that N 2 foams were more stable than CO 2 foams. It was observed that foam half-life decreased with the increase in temperature. Among all the surfactants, AOS foams showed the most promising thermal stability at high temperatures. Moreover, with the addition of stabilizing agents, foam's half-life was enhanced. Stabilizing agents such as crosslinker and GO dispersion showed the most stable foams with half-life recorded at 20min and 17min, respectively, at 200°C and 6.9MPa. Finally, pressure also showed a positive effect on foam stability; with increased pressure, foam half-life was increased. Based on the experimental data, analytical models for the effect of temperature and pressure were developed, considering foam degradation is a first-order kinetic reaction that linearly depends on the foam drainage mechanism. The effect of temperature on foam half-life was studied as an exponential decay model. In this model, foam half-life is a function of drainage rate constant (D A ) and activation energy (E a ) of the foam system. The effect of pressure on foam half-life was found to obey a power-law model where an increase in pressure showed an increase in foam half-life. Furthermore, a linear relation was studied for the effect of pressure on foam activation energy and drainage rate. Then the combined effects of temperature and pressure were studied, which yielded an analytical model to predict the foam stabilities in terms of half-life for different foam compositions. In conclusion, this research indicates that with an appropriate selection of surfactants and stabilizing agents, it is possible to obtain stable foams, which could replace conventional water fracturing fluid under EGS conditions.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

Geothermal play fairway analysis, part 1: Example from the Snake River Plain, Idaho

The Snake River Plain (SRP) volcanic province overlies the track of the Yellowstone hotspot, a thermal anomaly that extends deep into the mantle. Most of the area is underlain by a basaltic volcanic province that overlies a mid-crustal intrusive complex, which in turn provides the long-term heat flux needed to sustain geothermal systems. Previous studies have identified several known geothermal resource areas within the SRP. For the geothermal study presented herein, our goals were to: (1) adapt the methodology of Play Fairway Analysis (PFA) for geothermal exploration to create a formal basis for its application to geothermal systems, (2) assemble relevant data for the SRP from publicly available and private sources, and (3) build a geothermal PFA model for the SRP and identify the most promising plays, using GIS-based software tools that are standard in the petroleum industry. The study focused on identifying three critical resource parameters for exploitable hydrothermal systems in the SRP: heat source, reservoir and recharge permeability, and cap or seal. Data included in the compilation for heat source were heat flow, distribution and ages of volcanic vents, groundwater temperatures, thermal springs and wells, helium isotope anomalies, and reservoir temperatures estimated using geothermometry. Reservoir and recharge permeability was inferred from the analysis of stress orientations and magnitudes, post-Miocene faults, and subsurface structural lineaments based on magnetics and gravity data. Data for cap or seal included the distribution of impermeable lake sediments and clay-seal associated with hydrothermal alteration below the regional aquifer. These data were used to compile Common Risk Segment maps for heat, permeability, and seal, which were combined to create a Composite Common Risk Segment map for all southern Idaho that reflects the risk associated with geothermal resource exploration and identifies favorable resource tracks. Our regional assessment indicated that undiscovered geothermal resources may be located in several areas of the SRP. Two of these areas, the western SRP and Camas Prairie, were selected for more detailed assessment, during which heat, permeability, and seal were evaluated using newly collected field data and smaller grid parameters to refine the location of potential resources. These higher resolution assessments illustrate the flexibility of our approach over a range of scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

A note on the instability and pattern formation of shrinkage cracks in viscoplastic soils

In this work we present a theoretical study on the conditions for the onset of cracks, as well as the corresponding pattern formation, in saturated viscoplastic soils under isotropic loading (extension). The type of stress applied is left unspecified, to cover a variety of loadings including shrinkage due to desiccation, isotropic thermal expansion, mechanical loading and so forth. By treating the saturated soil as rigid viscoplastic, we obtain a 2D extension of the Cnoidal Waves equations (Veveakis and Regenauer-Lieb, 2015). By numerically solving the corresponding boundary value problem, we retrieve conditions for the onset of cracking instability in 2D loading, and identify the characteristic spacing between cracks to be a length scale combining all the hydro-mechanical parameters of the problem. Finally, we show that in a rectangular slab of clay under isotropic extension, patterns of triangular, rectangular and hexagonal cracks can tessellate the domain, with the hexagonal pattern being the energetically favored, as it minimizes the free energy of the system.

58 GEOSCIENCES↗
DOE OSTI2023

Advective gas flow in bentonite: Development and comparison of enhanced multi-phase numerical approaches

Understanding the impact of repository gas, generated from degradation of waste and its interaction with the host rock, is essential when assessing the performance and safety function of long-term disposal systems for radioactive waste. Numerical models based on conventional multi-phase flow theory have historically been applied to predict the outcome and impact of gas flow on different repository components. However, they remain unable to describe the full complexity of the physical processes observed in water-saturated experiments (e.g., creation of dilatant pathways) and thus, the development of novel representations for their description is required when assessing fully saturated clay-based systems. This was the primary focus of Task A within the international cooperative project DECOVALEX-2019 (D-2019) and refinement of these approaches is the primary focus of this study (Task B in the current phase of DECOVALEX-2023). This paper summarises development of enhanced numerical representations of key processes and compares the performance of each model against high-quality laboratory test data. Experimental data reveals that gas percolation in water-saturated compacted bentonite is characterised by four key features: (i) a quiescence phase, followed by (ii) the gas breakthrough, which leads to a (iii) peak value, which is then followed by (iv) a negative decay. Three models based on the multiphase flow theory have been developed. These models can provide good initial values and reasonable responses for gas breakthrough (although some of them still predict a too-smooth response). Peak gas pressure values are in general reasonably well captured, although maximum radial stress differences are observed at 48 mm from the base of the sample. Here, numerical peak values of 12.8 MPa are predicted, whereas experimental values are about 11 MPa. These models are also capable of providing a reasonable representation of the negative pressure decay following peak pressure. However, other key specific features (such as the timing of gas breakthrough) still require a better representation. The model simulations and their comparison with experimental data show that these models need to be further improved with respect to model parameter calibration, the numerical representation of spatial heterogeneities in material properties and flow localisation, and the upscaling of the related physical processes and parameters. To further understand gas flow localisation, a new conceptual model has been developed, which shows that discrete channels can possibly be induced through the instability of gas-bentonite interface during gas injection, thus providing a new perspective for modeling gas percolation in low-permeability deformable media.

58 GEOSCIENCES↗
DOE OSTI2024

Soft interface instability and gas flow channeling in low-permeability deformable media

Understanding gas percolation through a clay layer or a shale formation is of great importance for the development of a geologic repository for nuclear waste disposal, a subsurface system for gas storage, and an engineering approach for hydrocarbon extraction from unconventional reservoirs. Gas injection experiments have revealed complex dynamic behaviours of gas percolation through water saturated compacted bentonite, characterized by a high breakthrough pressure, rapid breakthrough, a pressure/stress decay after the breakthrough, a relatively high migration rate, high-frequency periodic/nonperiodic variations in flow rate, stepwise rate reductions during relaxation, and low gas saturation over the whole process, all indicating channelling nature of the processes. Using linear stability analyses, we show that this channelling can autonomously emerge from the instability of the deformable interface between the injected gas and the compacted bentonite matrix driven by local stress concentration, pore dilation, and hydrologic gradient. Channel patterns formed would possess a fractal geometry. We further show that, once a percolating channel is established, the gas injected would percolate through the channel in a chain of gas bubbles, also due to the interface instability, resulting in periodic/chaotic variations in gas flow rate. Our work provides a unified explanation for key features observed for gas percolation in low-permeability deformable media. The work also suggests a possibility of designing an engineered barrier system for a nuclear waste repository that can have controllable gas release while limit water transport.

Bentonite↗
DOE OSTI2023

Outgassing behavior and heat treatment optimization of JSC-1A lunar regolith simulant

As NASA strives towards a long duration presence on the Moon, it has become increasingly important to learn how to better utilize resources from the lunar surface for everything from habitats, vehicle infrastructure, and chemical extraction. To that end, a variety of lunar simulants have been sourced from terrestrially available volcanic minerals and glass as Apollo regolith is unavailable for experimentation needing large masses. However, while mineralogy and chemical composition can approach that of lunar material in these simulants, there are still distinct non-lunar phases such as hydrates, carbonates, sulfates, and clays that can cause simulants to behave distinctly non-lunar in a variety of processing conditions that maybe applied in-situ to lunar material. Notably, severe glassy bubbling has been documented in a variety of vacuum sintering experiments on JSC-1A lunar mare simulant heated via microwaves. The origins of this outgassing have not been well understood but are normally attributed to the decomposition of non-lunar contaminates intrinsic to virtually all terrestrially sourced simulants. As such, a series of controlled environmental tests were performed to ascertain the origins of the high temperature outgassing and to develop heat treatments that can drive JSC-1A closer to lunar composition and behavior. It was found that in JSC-1A at elevated temperatures distinct gas evolutions of water, carbon dioxide, and sulfur dioxide occur in both inert gas and vacuum. Additionally, the presence of hydrogen during heat treatments was shown to dramatically change gas evolutions, leading to distinctly more lunar-like composition and behavior from JSC-1A simulant.

58 GEOSCIENCES↗
DOE OSTI2025

Elevated-Mn ChemCam targets illuminating Mn redox cycling and diagenesis in the Bradbury Rise, Gale Crater, Mars

Manganese plays a crucial role as a paleo-environmental and geological indicator due to its sensitivity to redox potential and pH variations in the environment. On Earth, the association between the rise of atmospheric oxygen during the Great Oxidation Event and the presence of Mn in the sedimentary rock record underscores its significance. Here, in this study, we reexamined ChemCam targets from the first 600 sols of the Mars Science Laboratory mission, focusing on identifying instances of above-average Mn within these targets. These elevated-Mn targets were categorized into distinct geologic classes, revealing a pattern linking heightened Mn levels with diagenetically altered materials, such as calcium-sulfate veins and concretions, as well as clay minerals within the same targets, indicating a compelling relationship between Mn enrichment and diagenetic processes. High concentrations of Mn were observed in chemically altered targets, suggesting the occurrence of multiple fluid events: the first to alter the material and the second to deposit Mn. The observed patterns suggest multiple diagenetic events and redox cycling that facilitated the deposition and transport of Mn subsequent to the initial dissolution of basaltic materials. This research sheds light on the complexity of martian diagenetic processes and their implications for the planet’s environmental evolution.

58 GEOSCIENCES↗
DOE OSTI2020

Influence of geochemical reactions on the creep behavior of Mt. Simon sandstone

The capture and subsurface storage of carbon dioxide is a sustainable option that is currently pursued worldwide to mitigate greenhouse gas effect. However, predicting the long-term mechanical integrity of CO 2 underground storage systems remains a challenge. To address that question, it is essential to understand the influence of fluid-rock chemo-mechanical interactions on the long-term and on the time-dependent mechanical properties. In turn, the long-term mechanical response and the time-dependent mechanical behavior can be represented by the creep response. We investigate the impact of CO 2 -induced geochemical reaction on the creep response of Mt. Simon sandstone with a 50–400 micron grain size. We perform static and dynamic flow experiments on Mt. Simon sandstone specimens under geological conditions, at a temperature of T = 50 - 53 °C and for CO 2 pressures of P = 8:62; 17:2 MPa under both static flow and dynamic flow-through conditions. After aging, we employ creep indentation testing, high-resolution SEM-EDS, computer vision, machine learning, and micromechanics modeling to probe changes on the microstructure and mechanical properties. Following both static and dynamic flow-through experiments, we observe a 10–22% decrease in quartz volume fraction and an increase in both the microporosity (7–28%) and nanoporosity (60–65%). Additional CO 2 -induced microstructural changes include an enlargement of pore throats and the formation of channels. These observations point to the presence of CO 2 -induced K-feldspar dissolution and clay dissolution reactions. Here, the macroscopic creep behavior is logarithmic and the macroscopic creep modulus varies depending on the microporosity and the relative quartz and feldspar content. As a result of these geochemical reactions and of the related microstructural changes, a 55-60% decrease in the macroscopic logarithmic creep modulus is predicted.

58 GEOSCIENCES↗
DOE OSTI2024

Carbon sequestration into lime-stabilized soils: Engineering and mineralogical characterization

Stabilizing weak clayey soils with lime is an effective method for improving the mechanical properties of soil. However, lime production is an energy-intensive process producing significant CO 2 emissions in lime-stabilized soils, which can be counteracted through accelerated carbonation that enhances its engineering performance. The present study evaluates accelerated carbonation of lime-treated soils by adding gaseous (CO 2 -rich gas), liquid (water-CO 2 mixture), and solid (sodium bicarbonate) CO 2 sources. Results indicated that samples carbonated with gaseous CO 2 exhibited 100% lime carbonation, while samples treated with solid and aqueous sources of CO 2 had a mean lime carbonation of 60% and 40%, respectively. Further, all lime-treated-carbonated samples exhibited a mean 50% increase in unconfined compressive strength compared to the untreated samples after a 7-day curing period. Durability evaluation through cyclic wetting and drying indicated that the carbonated samples had higher durability than the untreated samples. X-ray computed tomography showed that adding solid and liquid sources of CO 2 facilitated the flocculation of montmorillonite, reducing the porosity. However, a higher dosage of solid CO 2 induced clay dispersion, increasing the porosity. X-ray diffraction and thermogravimetric analysis verified CO 2 sequestration through the formation of calcite, a thermodynamically stable polymorph of calcium carbonate.

58 GEOSCIENCES↗
DOE OSTI2020

In situ precipitation of hydrous ferric oxide (HFO) for remediation of subsurface iodine contamination

A practical approach for in-situ hydrous ferric oxide (HFO) precipitation was developed for iodine immobilization under field-scale conditions at the Hanford Site. A series of 1D meter-long bench top column experiments packed with Hanford sediments was conducted with a single acidic ferric solution (0.1 M, pH = 1.5) injection. Because carbonate and clay minerals are widely present in sediments, self-pH buffering of the injected acidic ferric solution occurred due to mineral dissolution, leading to HFO precipitation under a neutral condition. Up to ~12 mg/g HFO successfully precipitated and evenly distributed in the column sediments, and the remobilization of the neo-formed HFO precipitates was limited (< 4 wt.% after 100 pore volumes of flushing). The transport of IO3- in the HFO-amended sediments was strongly retarded through both adsorption and co-precipitation processes. However, reversible adsorption of iodine on HFO was observed, which might limit its application to slow moving groundwater systems.

Wang, Guohui↗
DOE OSTI2022

Logistics simulation of a remediation effort for a hypothetical radiological contamination scenario

To mitigate the effects following a large-scale nuclear or radiological material release in an urban environment and to expedite recovery, the Integrated Wash-Aid Treatment Emergency Reuse System (IWATERS) was developed. IWATERS consists of three operations: washing contaminated surfaces with an ionic wash solution, collecting, and treating the contaminated wash solution on-site to remove contaminants, and reusing the treated solution throughout operations to preserve the clean water resource. This study develops a framework to simulate the logistics of IWATERS deployment, thereby gaining an understanding of the timeline for decontamination operations. For this purpose, the Analysis of Mobility Platform and GoldSim were leveraged for a hypothetical contamination scenario covering 65,200 m 2 of an urban center. The framework reveals that remediation progress is limited by several resources, notably the availability of vermiculite, a reactive clay that is required to treat the contaminated wash solution. Further, this study also presents how the simulation approach can be used to characterize alternatives to reduce the influence of limited resources on operational progress. Overall, this work lays the foundation for evaluating different decontamination methods through detailed logistics simulation, i.e., by refining simulation assumptions and expanding the range of scenarios the simulation can depict.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2023

Evolving radon diffusion through earthen barriers at uranium waste disposal sites

Field measurements of Rn-222 fluxes from the tops and bottoms of compacted clay radon barriers were used to calculate effective Rn diffusion coefficients (D Rn ) at four uranium waste disposal sites in the western United States to assess cover performance after more than 20 years of service. Values of D Rn ranged from 7.4 × 10 -7 to 6.0 × 10 -9 m 2 /s, averaging 1.42 × 10 -7 . Water saturation (S W ) from soil cores indicated that there was relatively little control of D Rn by S W , especially at higher moisture levels, in contrast to estimates from most steady-state diffusion models. Further, this is attributed to preferential pathways intrinsic to construction of the barriers or to natural process that have developed over time including desiccation cracks, root channels, and insect burrows in the engineered earthen barriers. A modification to some models in which fast and slow pathway D Rn values are partitioned appears to give a good representation of the data; 4% of the fast pathway was needed to fit the data regression. For locations with high S w and highest D Rn (and fluxes) at each site, the proportion of fast pathway ranged from 1.7% to 34%, but for many locations with lower fluxes, little if any fast pathway was needed.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Analysis of time-dependent strain heterogeneity in shales using X-ray microscopy and digital volume correlation

In this paper, we report the results of a 4-dimensional (3 spatial dimensions plus time) in-situ creep experiment on a milli-core shale plug from Barnett formation in Texas. An outstanding challenge in shale mechanics is connecting observable characteristics to predictable mechanical and hydraulic behavior. We employ novel imaging and analysis techniques to statistically link sample mineralogy and structure to observed deformation. In particular, we focus on mechanisms of porosity reduction during transient response to step changes in stress, analogous to the aseismic deformation that occurs after an episode of fracturing. Our main goal in this study is to take the first steps to visualize and capture the small-scale deformation mechanisms involved in creep of shales, using imaging techniques. To do so, we conducted a uniaxial creep experiment inside an Xradia Versa 520 micro-CT system for 12 h at room temperature under a constant uniaxial stress of 25 MPa. We employed a Digital Volume Correlation (DVC) method to measure micro-scale deforspmation of the sample by comparing the resulting volumetric images with a high-resolution scan of the sample obtained prior to the creep experiment. The results indicate that the time-dependent deformation in the compliant clay plus kerogen matrix constituents facilitates rigid grain rotation and compaction of intergranular porosity at the boundary between rigid grains and matrix.

58 GEOSCIENCES↗
DOE OSTI2022

Effect of zeolite type, temperature, and pH on Stage III glass alteration behavior for two nuclear waste glasses

We report long-term laboratory scale static alteration tests (up to 1500 d) were conducted on two low-activity waste glasses (LAWA44 and IDF18-A161) at S/V = 2000 m –1 . Tests were seeded with either analcime, clinoptilolite, zeolite P1, or zeolite P2 to induce Stage III glass alteration behavior. The effect of several parameters on glass Stage III behavior was studied with consideration of several variables, including temperature (22 °C – 90 °C) and pH (floating or initially set using KOH at values from 9.5 to 11.5). Temporal solution data and solids characterization at the end of the experiments are reported. Solution data demonstrated the occurrence of a sustained rate acceleration at temperatures down to 40 °C and the activation energy for these Stage III rates was determined for both glasses. Three major secondary phases were identified at the end of the experiments for both glasses: zeolite P1, zeolite P2, and/or analcime (along with minor phases: zeolites, clays, carbonates). The identification of the crystalline phases at the end of the experiment often differed from the zeolite seed phase. In addition, the use of the Avrami equation showed unimpeded zeolite growth (3-dimensional) in certain conditions, suggesting that in certain conditions the glass dissolution rate is controlled by zeolite formation kinetics whereas in other conditions the glass corrosion process impeded zeolite formation kinetics. The morphological evolution of the samples revealed the growth of the zeolites which completely fill the intergranular space in the powder bed. These results are discussed with regards to previous seeded and unseeded Stage III dissolution rate experiments performed on low- and high-activity nuclear waste glasses.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Leaching of lanthanide and yttrium from a Central Appalachian coal and the ashes obtained at 550–950 °C

Here in this work, we investigated leaching of lanthanide and yttrium (REY) from a Central Appalachian coal and its ashes obtained at 550–950 °C with the main purpose of understanding the impact of ashing temperature on REY leachability in water, ammonium sulfate, and hydrochloric acid. It is found that the coal contains a negligible amount of water-soluble REY, less than 1% ion-exchangeable REY, and about 28% of HCl-soluble REY. Ashing leads to dramatic changes in REY leachability in both ammonium sulfate and hydrochloric acid solutions, which is believed to be related to transformation and redistribution of organically-associated REY in coal during the ashing process. Ashing temperature significantly affects REY leaching from coal ashes; higher ashing temperature results in lower REY leachability in both solutions. Clay minerals may play a significant role in changing the leachability of REY after ashing. In addition, the results also suggest that the organic matter in the coal is relatively enriched in heavy REY.

01 COAL, LIGNITE, AND PEAT↗