Theoretical analysis and engineering application of controllable shock wave technology for enhancing coalbed methane in soft and low‑permeability coal seams

Coalbed methane(CBM)is a significant factor in triggering coal and gas outburst disaster,while also serving as a clean fuel.With the increasing depth of mining operations,coal seams that exhibit high levels of gas content and low permeability have become increasingly prevalent.While controllable shockwave(CSW)technology has proven effective in enhancing CBM in laboratory settings,there is a lack of reports on its field applications in soft and low-permeability coal seams.This study establishes the governing equations for stress waves induced by CSW.Laplace numerical inversion was employed to analyse the dynamic response of the coal seam during CSW antireflection.Additionally,quantitative calculations were performed for the crushed zone,fracture zone,and effective CSW influence range,which guided the selection of field test parameters.The results of the field test unveiled a substantial improvement in the gas permeability coefficient,the average rate of pure methane flowrate,and the mean gas flowrate within a 10 m radius of the antireflection borehole.These enhancements were notable,showing increases of 3 times,13.72 times,and 11.48 times,respectively.Furthermore,the field test performed on the CSW antireflection gas extraction hole cluster demonstrated a noticeable improvement in CBM extraction.After antireflection,the maximum peak gas concentration and maximum peak pure methane flow reached 71.2%and 2.59 m^(3)/min,respectively.These findings will offer valuable guidance for the application of CSW antireflection technology in soft and low-permeability coal seams.

A stable implicit nodal integration-based particle finite element method(N-PFEM)for modelling saturated soil dynamics

In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a generalised Hellinger-Reissner(HR)variational principle,creating an implicit PFEM formulation.To mitigate the volumetric locking issue in low-order elements,we employ a node-based strain smoothing technique.By discretising field variables at the centre of smoothing cells,we achieve nodal integration over cells,eliminating the need for sophisticated mapping operations after re-meshing in the PFEM.We express the discretised governing equations as a min-max optimisation problem,which is further reformulated as a standard second-order cone programming(SOCP)problem.Stresses,pore water pressure,and displacements are simultaneously determined using the advanced primal-dual interior point method.Consequently,our numerical model offers improved accuracy for stresses and pore water pressure compared to the displacement-based PFEM formulation.Numerical experiments demonstrate that the N-PFEM efficiently captures both transient and long-term hydro-mechanical behaviour of saturated soils with high accuracy,obviating the need for stabilisation or regularisation techniques commonly employed in other nodal integration-based PFEM approaches.This work holds significant implications for the development of robust and accurate numerical tools for studying saturated soil dynamics.

Numerical parametric study on the influence of location and inclination of large-scale asperities on the shear strength of concreterock interfaces of small buttress dams

When assessing the sliding stability of a concrete dam,the influence of large-scale asperities in the sliding plane is often ignored due to limitations of the analytical rigid body assessment methods provided by current dam assessment guidelines.However,these asperities can potentially improve the load capacity of a concrete dam in terms of sliding stability.Although their influence in a sliding plane has been thoroughly studied for direct shear,their influence under eccentric loading,as in the case of dams,is unknown.This paper presents the results of a parametric study that used finite element analysis(FEA)to investigate the influence of large-scale asperities on the load capacity of small buttress dams.By varying the inclination and location of an asperity located in the concrete-rock interface along with the strength of the rock foundation material,transitions between different failure modes and correlations between the load capacity and the varied parameters were observed.The results indicated that the inclination of the asperity had a significant impact on the failure mode.When the inclinationwas 30and greater,interlocking occurred between the dam and foundation and the governing failure modes were either rupture of the dam body or asperity.When the asperity inclination was significant enough to provide interlocking,the load capacity of the dam was impacted by the strength of the rock in the foundation through influencing the load capacity of the asperity.The location of the asperity along the concrete-rock interface did not affect the failure mode,except for when the asperity was located at the toe of the dam,but had an influence on the load capacity when the failure occurred by rupture of the buttress or by sliding.By accounting for a single large-scale asperity in the concrete-rock interface of the analysed dam,a horizontal load capacity increase of 30%e160%was obtained,depending on the inclination and location of the asperity and the strength of the foundation material.

Numerical analysis of downward progressive landslides in long natural slopes with sensitive clay

Landslides occurring in sensitive clay often result in widespread destruction,posing a significant risk to human lives and property due to the substantial decrease in undrained shear strength during deformation.Assessing the consequences of these landslides is challenging and necessitates robust numerical methods to comprehensively investigate their failure mechanisms.While studies have extensively explored upward progressive landslides in sensitive clays,understanding downward progressive cases remains limited.In this study,we utilised the nodal integration-based particle finite element method(NPFEM)with a nonlinear strain-softening model to analyse downward progressive landslides in sensitive clay on elongated slopes,induced by surcharge loads near the crest.We focused on elucidating the underlying failure mechanisms and evaluating the effects of different soil parameters and strainsoftening characteristics.The simulation results revealed the typical pattern for downward landslides,which typically start with a localised failure in proximity to the surcharge loads,followed by a combination of different types of failure mechanisms,including single flow slides,translational progressive landslides,progressive flow slides,and spread failures.Additionally,inclined shear bands occur within spread failures,often adopting distinctive ploughing patterns characterised by triangular shapes.The sensitive clay thickness at the base,the clay strength gradient,the sensitivity,and the softening rate significantly influence the failure mechanisms and the extent of diffused displacement.Remarkably,some of these effects mirror those observed in upward progressive landslides,underscoring the interconnectedness of these phenomena.This study contributes valuable insights into the complex dynamics of sensitive clay landslides,shedding light on the intricate interplay of factors governing their behaviour and progression.

Evolution of nonmetallic inclusions in 80-t 9CrMoCoB large-scale ingots during electroslag remelting process

In combination with theoretical calculations,experiments were conducted to investigate the evolution behavior of nonmetallic inclusions(NMIs)during the manufacture of large-scale heat-resistant steel ingots using 9CrMoCoB heat-resistant steel and CaF_(2)–CaO–Al_(2)O_(3)–SiO_(2)–B_(2)O_(3)electroslag remelting(ESR)-type slag in an 80-t industrial ESR furnace.The main types of NMI in the consumable electrode comprised pure alumina,a multiphase oxide consisting of an Al_(2)O_(3)core and liquid CaO–Al_(2)O_(3)–SiO_(2)–MnO shell,and M_(23)C_(6)carbides with an MnS core.The Al_(2)O_(3)and MnS inclusions had higher precipitation temperatures than the M_(23)C_(6)-type carbide under equilibrium and nonequilibrium solidification processes.Therefore,inclusions can act as nucleation sites for carbide layer precipitation.The ESR process completely removed the liquid CaO–Al_(2)O_(3)–SiO_(2)–MnO oxide and MnS inclusion with a carbide shell,and only the Al_(2)O_(3)inclusions and Al_(2)O_(3)core with a carbide shell occupied the remelted ingot.The M_(23)C_(6)-type carbides in steel were determined as Cr_(23)C_(6)based on the analysis of transmission electron microscopy results.The substitution of Cr with W,Fe,or/and Mo in the Cr_(23)C_(6)lattice caused slight changes in the lattice parameter of the Cr_(23)C_(6)carbide.Therefore,Cr_(21.34)Fe_(1.66)C_(6),(Cr_(19)W_(4)C_(6),Cr_(18.4)Mo_(4.6)C_(6),and Cr_(16)Fe_(5)Mo_(2)C_(6)can match the fraction pattern of Cr_(23)C_(6)carbide.The Al_(2)O_(3)inclusions in the remelted ingot formed due to the reduction of CaO,SiO_(2),and MnO components in the liquid inclusion.The increased Al content in liquid steel or the higher supersaturation degree of Al_(2)O_(3)precipitation in the remelted ingot than that in the electrode can be attributed to the evaporation of CaF_(2)and the increase in CaO content in the ESR-type slag.

Effect of hafnium and molybdenum addition on inclusion characteristics in Co-based dual-phase high-entropy alloys

Specific grades of high-entropy alloys(HEAs)can provide opportunities for optimizing properties toward high-temperature applications.In this work,the Co-based HEA with a chemical composition of Co_(47.5)Cr_(30)Fe_(7.5)Mn_(7.5)Ni_(7.5)(at%)was chosen.The refractory metallic elements hafnium(Hf)and molybdenum(Mo)were added in small amounts(1.5at%)because of their well-known positive effects on high-temperature properties.Inclusion characteristics were comprehensively explored by using a two-dimensional cross-sectional method and extracted by using a three-dimensional electrolytic extraction method.The results revealed that the addition of Hf can reduce Al_(2)O_(3)inclusions and lead to the formation of more stable Hf-rich inclusions as the main phase.Mo addition cannot influence the inclusion type but could influence the inclusion characteristics by affecting the physical parameters of the HEA melt.The calculated coagulation coefficient and collision rate of Al_(2)O_(3)inclusions were higher than those of HfO_(2)inclusions,but the inclusion amount played a larger role in the agglomeration behavior of HfO_(2)and Al_(2)O_(3)inclusions.The impurity level and active elements in HEAs were the crucial factors affecting inclusion formation.

Color and Gloss Changes of a Lignin-Based Polyurethane Coating under Accelerated Weathering

The purpose of this research study was to investigate the properties of polyurethane coatings based on lignin nano-particles.For this purpose,the prepared coatings were applied to pine wood surfaces and weathered artificially.Subsequently,color and gloss of the coatings were measured before and after the weathering test.Field emission scanning electron microscopy(FE-SEM)micrographs prepared from the coatings showed that the average size of nano-particles in the polyurethane substrate was approximately 500 nm.Nuclear magnetic resonance(13C-NMR)spectroscopy showed that strong urethane bonds were formed in the nano-lignin-based polyurethane.Differential calorimetric analysis(DSC)test revealed that the glass-transition temperature(Tg)of lignin nanoparticles modified with diethylenetriamine(DETA)was 112.8℃ and Tg of lignin nano-particles modified with ethylenediamine(EDA)was 102.5℃,which is lower than the Tg of un-modified lignin(114.6℃)and lignin modified with DETA(126.8℃)and lignin modified with EDA(131.3℃).The coatings modified with lignin nano-particles had a greater change in gloss.The lignin nano-particles in the modified coating are trapping hydroxyl radicals which reduces photoactivity and yellowing of the polyurethane by about 3 times compared to unmodified polyurethane coatings.After weathering test,the nano-lignin-based coating had a rougher surface with a lower contact angle(0.78°)compared to the unmodified polyurethane coating(0.85°).

Advances and challenges of electrolyzers for large-scale CO_(2) electroreduction

CO_(2) electroreduction(CO_(2) ER)to high value-added chemicals is considered as a promising technology to achieve sustainable carbon neutralization.By virtue of the progressive research in recent years aiming at design and understanding of catalytic materials and electrolyte systems,the CO_(2) ER performance(such as current density,selectivity,stability,CO_(2) conversion,etc.)has been continually increased.Unfortunately,there has been relatively little attention paid to the large-scale CO 2 electrolyzers,which stand just as one obstacle,alongside series-parallel integration,challenging the practical application of this infant technology.In this review,the latest progress on the structures of low-temperature CO_(2) electrolyzers and scale-up studies was systematically overviewed.The influence of the CO_(2) electrolyzer configurations,such as the flow channel design,gas diffusion electrode(GDE)and ion exchange membrane(IEM),on the CO_(2) ER performance was further discussed.The review could provide inspiration for the design of large-scale CO_(2) electrolyzers so as to accelerate the industrial application of CO_(2) ER technology.

Viscous behavior of 1-hexyl-methylimidazolium bis(trifluoromethylsulfonyl)imide/titanium dioxide/polyethylene glycol

Viscous behavior is important for the process design, especially for the non-Newtonian fluid. In this study, the viscous behaviors of slurry, i.e., 1-hexyl-3-methylimidazolium bis(trifluoromethylsulfonyl)im ide([Hmim][NTf_(2)])/titanium dioxide(TiO_(2))-polyethylene glycol(PEG200), were determined experimentally and systematically. The pressure drop was estimated when [Hmim][NTf_(2)]/TiO_(2)-PEG200 was used as the solvent in the absorption/desorption towers. The results show that the slurry belongs to the non-Newtonian fluid with shear-thinning behavior. High temperature and low solid content are beneficial to reduce the viscosity of [Hmim][NTf_(2)]/TiO_(2)-PEG200, and the presence of [Hmim][NTf_(2)] can effectively reduce the viscosity of the slurry. In addition, high temperature is preferable for reducing the pressure drop, and the pressure drop of slurry with the solid content value of 8.0%(mass) can reduce by 28.0%when the temperature increases from 313 to 333 K.

Modelling erosion of a single rock block using a coupled CFD-DEM approach

Rock block removal is the prevalent physical mechanism for rock erosion and could affect the stability of dam foundations and spillways.Despite this,understanding of block removal is still inadequate because of the complex interactions among block characteristics,hydraulic forces,and erosive processes acting on the block.Herein,based on a previously conducted physical experiment of erosion of a single rock block,the removal processes of two different protruding blocks are represented by a coupled computational fluid dynamics-discrete element model(CFD-DEM)approach under varied flow conditions.Additionally,the blocks could be rotated with respect to the flow direction to consider the effect of the discontinuity orientation on the block removal process.Simulation results visualize the entire block removal process.The simulations reproduce the effects of the discontinuity orientation on the critical flow velocity inducing block incipient motion and the trajectory of the block motion observed in the physical experiments.The numerical results present a similar tendency of the critical velocities at different discontinuity orientations but have slightly lower values.The trajectory of the block in the simulations fits well with the experimental measurements.The relationship between the dimensionless critical shear stress and discontinuity orientation observed from the simulations shows that the effect of block protrusion becomes more dominant on the block incipient motion with the increase of relative protrusion height.To our knowledge,this present study is the first attempt to use the coupled finite volume method(FVM)-DEM approach for modelling the interaction behavior between the block and the flowing water so that the block removal process can be reproduced and analyzed.

Efficient SO_(2)removal using aqueous ionic liquid at low partial pressure

The development of novel absorbents is essential for SO_(2)removal.In this study,a novel ionic liquid(IL,[BHEP][HSO_(4)])was prepared,and water was selected as the co-solvent.The density and viscosity of aqueous[BHEP][HSO_(4)]were measured and the SO_(2)absorption performance was systematically investigated.Furthermore,the thermodynamic properties of SO_(2)in aqueous[BHEP][HSO_(4)]were calculated.Additionally,the mechanism of SO_(2)absorption in aqueous[BHEP][HSO_(4)]was confirmed using Fouriertransform infrared and nuclear magnetic resonance spectroscopy.It showed that[BHEP][HSO_(4)]absorbed0.302 g·g^(-1)(g SO_(2)/g IL)at an SO_(2)partial pressure of 2000μl·L^(-1)at 303.2 K,and the SO_(2)desorption enthalpy was-39.63 k J·mol^(-1).The mechanistic study confirmed the chemical absorption of SO_(2)in aqueous[BHEP][HSO_(4)].

Polysaccharides-based pyrite depressants for green flotation separation:An overview

Froth flotation is an essential processing technique for upgrading low-grade ores.Flotation separation would not be efficient without chemical surfactants(collectors,depressants,frothers,etc.).Depressants play a critical role in the selective separation of minerals in that they deactivate unfavorable mineral surfaces and hinder them from floating into the flotation concentration zone.Pyrite is the most common and challenging sulfide gangue,and its conventional depressants could be highly harmful to nature and humans.Therefore,using available,affordable,eco-friendly polymers to assist or replace hazardous reagents is mandatory for a green transition.Polysaccharide-based(starch,dextrin,carboxymethyl cellulose,guar gum,etc.)polymers are one of the most used biodegradable depressant groups for pyrite depression.Despite the satisfactory flotation results obtained using these eco-friendly depressants,several gaps still need to be addressed,specifically in investigating surface interactions,adsorption mechanisms,and parameters affecting their depression performance.As a unique approach,this review comprehensively discussed previously conducted studies on pyrite depression with polysaccharide-based reagents.Additionally,practical suggestions have been provided for future assessments and developments of polysaccharide-based depressants,which pave the way to green flotation.This robust review also explored the depression efficiency and various adsorption aspects of naturally derived depressants on the pyrite surface to create a possible universal trend for each biodegradable depressant derivative.

Evaluation of Actual Evapotranspiration and Crop Coefficient in Carrot by Remote Sensing Methodology Using Drainage and River Water to Overcome Reduced Water Availability

Searching for alternative methods for traditional irrigation is World trend at days due to a reduction in water and increased of drought due to climate changes therefore farmers need use modern methods of scheduling water and minimizing water losses while also increasing yield. To meet the future increasing demands water and food there is a need to utilize alternative methods to reduce evaporation, transpiration and deep percolation of water. Any countries use recycled water (drain and sewage) and desalination water from the sea or drains to irrigate crops plus computing actual crop evapotranspiration (ETc) so as to calculate the amount of water to apply to a crop. The paper aims to assess the actual evaporation and evaporation coefficient of carrots, by planting carrots in a field and the crop was exposed to several sources of water (DW and RW) and comparing ETc, Kc and production among plots of three sites (A, B and C). The study used two types of irrigation water (drain water (DW) and river water (RW)). The results were to monthly rate and accumulated actual evapotranspiration to C (irrigation by RW only) more than A (67% RW and 33% DW) and B (17% RW and 83% DW) via 7% and 58%, respectively. The yield to C more than A and B by 17% and 75%, respectively. In conclusion the use of DW can cause a reduction in crop consumptive of carrot crops also causes a reduction in yield, crop length, root length, root size, canopy of crop, number of leaves and biomass of the plant therefore, the drainage water needs to treated before irrigating crops And making use of it to irrigate the fields and fill the shortfall in the amount of water from the river. The drain water helped on filling the water shortage due to climate changes and giving production of carrot crop but less than river water.

Quantitative investigation on micro-parameters of cemented paste backfill and its sensitivity analysis

The mechanical properties of cemented paste backfill(CPB) depend heavily on its pore structural characteristics and micro-structural changes. In order to explore the variation mechanisms of macro-mechanical characteristics and micro-structure of CPB. CPB specimens with different mass concentrations prepared from the full tailings of Xianglushan Tungsten Ore were micro-tests. Moreover, acquired pore digital images were processed by using the pores(particles) and cracks analysis system(PCAS), and a sensitivity analysis was performed. The results show that as the mass concentration of CPB increases from 70% to 78%, the porosity, the average pore area and the number of pores drop overall, leading to a decline in the pores opening degree and enhancing the mechanical characteristics. As the mass concentration of CPB increases, the trend of fractal dimension, probability entropy and roundness is reduced, constant and increased, which can result in an enhancement of the uniformity, an unchanged directionality and more round pores. According to the definition of sensitivity, the sensitivities of various micro-parameters were calculated and can be ranked as porosity > average pore area > number of pores > roundness > fractal dimension > probability entropy.

Interpretable modeling of metallurgical responses for an industrial coal column flotation circuit by XGBoost and SHAP-A “conscious-lab” development

Surprisingly,no investigation has been explored relationships between operating variables and metallurgical responses of coal column flotation(CF) circuits based on industrial databases for under operation plants.As a novel approach,this study implemented a conscious-lab "CL" for filling this gap.In this approach,for developing the CL dedicated to an industrial CF circuit,SHapley Additive explanations(SHAP) and extreme gradient boosting(XGBoost) were powerful unique machine learning systems for the first time considered.These explainable artificial intelligence models could effectively convert the dataset to a basis that improves human capabilities for better understanding,reasoning,and planning the unit.SHAP could provide precise multivariable correlation assessments between the CF dataset by using the Tabas Parvadeh coal plant(Kerman,Iran),and showed the importance of solid percentage and washing water on the metallurgical responses of the coal CF circuit.XGBoost could predict metallurgical responses(R-square> 0.88) based on operating variables that showed quite higher accuracy than typical modeling methods(Random Forest and support vector regression).

Blast-induced ground vibration prediction in granite quarries:An application of gene expression programming,ANFIS,and sine cosine algorithm optimized ANN

Blasting of rocks has intrinsic environmental impacts such as ground vibration,which can interfere with the safety of lives and property.Hence,accurate prediction of the environmental impacts of blasting is imperative as the empirical models are not accurate as evident in the literature.Therefore,there is need to consider some robust predictive models for accurate prediction results.Gene expression programming(GEP),adaptive neuro-fuzzy inference system(ANFIS),and sine cosine algorithm optimized artificial neural network(SCA-ANN)models are proposed for predicting the blast-initiated ground vibration in five granite quarries.The input parameters into the models are the distance from the point of blasting to the point of measurement(D),the weight of charge per delay(W),rock density(q),and the Schmidt rebound hardness(SRH)value while peak particle velocity(PPV)is the targeted output.100 datasets were used in developing the proposed models.The performance of the proposed models was examined using the coefficient of determination(R2)and error analysis.The R2 values obtained for the GEP,ANFIS,and SCA-ANN models are 0.989,0.997,and 0.999,respectively,while their errors are close to zero.The proposed models are compared with an empirical model and are found to outperform the empirical model.

Three-dimensional upper bound limit analysis of underground cavities using nonlinear Baker failure criterion

A generalized nonlinear Baker failure criterion is employed with the upper bound limit analysis to study the surrounding rock stability of underground cavities. A three-dimensional(3D) failure mode is established by extending the two-dimensional(2D) failure mode, which offers an upper bound expression of the surrounding rock pressure. This method is validated with a series of examples before the influence of four parameters of scale parameter, curvature parameter, shift parameter and lateral pressure coefficient, on the surrounding rock pressure is analyzed. According to these results, failure ranges of the underground cavities are determined. The following conclusions are reached:(1) the proposed approach is more accurate to predict surrounding rock pressure than the Mohr-Coulomb failure criterion;(2) the surrounding rock with large scale parameter, curvature parameter, shift parameter, and lateral pressure coefficient can lead to a more stable underground cavity;(3) the failure range in 3D mode can be predicted according to the upper bound solutions.

Ionic liquids/deep eutectic solvents for CO_(2) capture:Reviewing and evaluating

The CO_(2)solubilities(including CO_(2)Henry’s constant)in physical-and chemical-based ILs/DESs and the COSMO-RS models describing these properties were comprehensively collected and summarized.The summarized results indicate that chemical-based ILs/DESs are superior to physical-based ILs/DESs for CO_(2)capture,especially those ILs have functionalized cation and anion,and superbase DESs;some of the superbase DESs have higher CO_(2)solubilities than those of ILs;the best physical-and chemical-based ILs,as well as physical-and chemicalbased DESs are[BMIM][BF4](4.20 mol kg^(-1)),[DETAH][Im](11.91 mol kg^(-1)),[L-Arg]-Gly 1:6(4.92 mol kg^(-1))and TBD-EG 1:4(12.90 mol kg^(-1)),respectively.Besides the original COSMO-RS mainly providing qualitative predictions,six corrected COSMO-RS models have been proposed to improve the prediction performance based on the experimental data,but only one model is with universal parameters.The newly determined experimental results were further used to verify the perditions of original and corrected COSMO-RS models.The comparison indicates that the original COSMO-RS qualitatively predicts CO_(2)solubility for some but not all ILs/DESs,while the quantitative prediction is incapable at all.The original COSMO-RS is capable to predict CO_(2)Henry’s constant qualitatively for both physical-based ILs and DESs,and quantitative prediction is only available for DESs.For the corrected COSMO-RS models,only the model with universal parameters provides quantitative predictions for CO_(2)solubility in physical-based DESs,while other corrected models always show large deviations(>83%)compared with the experimental CO_(2)Henry’s constants.

Ionic liquids for CO_(2) electrochemical reduction

Electrochemical reduction of CO_(2) is a novel research field towards a CO_(2)-neutral global economy and combating fast accelerating and disastrous climate changes while finding new solutions to store renewable energy in value-added chemicals and fuels.Ionic liquids(ILs),as medium and catalysts(or supporting part of catalysts)have been given wide attention in the electrochemical CO_(2) reduction reaction(CO_(2) RR)due to their unique advantages in lowering overpotential and improving the product selectivity,as well as their designable and tunable properties.In this review,we have summarized the recent progress of CO_(2) electro-reduction in IL-based electrolytes to produce higher-value chemicals.We then have highlighted the unique enhancing effect of ILs on CO_(2) RR as templates,precursors,and surface functional moieties of electrocatalytic materials.Finally,computational chemistry tools utilized to understand how the ILs facilitate the CO_(2) RR or to propose the reaction mechanisms,generated intermediates and products have been discussed.

Prediction and verification of heat capacities for pure ionic liquids

The heat capacity of ionic liquids is an important physical property,and experimental measuring is usually used as a common method to obtain them.Owing to the huge number of ionic liquids that can be potentially synthesized,it is desirable to acquire theoretical predictions.In this work,the Conductor-like Screening Model for Real Solvents(COSMO-RS)was used to predict the heat capacity of pure ionic liquids,and an intensive literature survey was conducted for providing a database to verify the prediction of COSMO-RS.The survey shows that the heat capacity is available for 117 ionic liquids at temperatures ranging 77.66-520 K since 2004,and the 4025 data points in total with the values from 76.37 to 1484 J·mol^(-1)·K^(-1) have been reported.The prediction of heat capacity with COSMO-RS can only be conducted at two temperatures(298 and 323 K).The comparison with the experimental data proves the prediction reliability of COSMO-RS,and the average relative deviation(ARD)is 8.54%.Based on the predictions at two temperatures,a linear equation was obtained for each ionic liquid,and the heat capacities at other temperatures were then estimated via interpolation and extrapolation.The acquired heat capacities at other temperatures were then compared with the experimental data,and the ARD is only 9.50%.This evidences that the heat capacity of a pure ionic liquid follows a linear equation within the temperature range of study,and COSMO-RS can be used to predict the heat capacity of ionic liquids reliably.