Carnot battery (CB) is a novel electricity storage technology which could realize electricity storage and waste heat recovery simultaneously. The thermal integrated CB (TI-CB) only utilizes the waste heat energy during the charging process in previous study. The waste heat energy utilization is insufficient. A CB with thermal integration during charging and discharging process (DTI-CB) is proposed in this paper. The waste heat energy is effectively integrated with DTI-CB in all periods through being coupled to the Organic Rankine cycle (ORC). A thermo-economic evaluation of the proposed DTI-CB is conducted. Compared with TI-CB, the power output capacity could be increased 95.67% and the Levelized cost of storage could be reduced 30.90%. The thermoeconomic performance of the proposed DTI-CB performs better than those of the single ORC and CB under the conditions with relatively higher heat utilization load. 

The thermal conductivity of soil is a vital parameter in describing its heat transfer properties. Accurate prediction and sensitivity analysis of this parameter can aid in assessing the thermal response in geotechnical engineering and prevent deformation and damage in projects. Based on thermal conductivity experiments by Kersten’s team, we analyzed the factors influencing this parameter. we considered introducing a temperature variable into the traditional empirical formula and conducted validation, resulting in an improved formula with good applicability to clay. Using artificial intelligence algorithms, we established a prediction model for thermal conductivity. The model uses soil type, dry density, water content, and temperature as input variables. Our analysis showed that the Random Forest model, Radial Basis Function Neural Network (RBFNN), and Whale Optimization Algorithm Backpropagation Neural Network (WOA-BP) could all accurately predict thermal conductivity. Among these, the WOA-BP model demonstrated the best performance, followed by Random Forest and RBFNN.We tested the prediction model using a new sample set and found that the model still performed well, indicating a certain level of generalization ability. We employed a Monte Carlo simulation for parameter sensitivity analysis of the improved empirical formula. With the Random Forest model, we ranked feature importance to evaluate the impact of different input variables on model output. Finally, we calculated the sensitivity of influencing factors using a weighted product method combined with WOA-BP. The results from all three methods were consistent. They indicated that the sensitivity of thermal conductivity to changes decreases in the order of water content, dry density, temperature, and soil type.

Chemical looping can achieve high product selectivity using lattice oxygen in oxygen carriers for partial oxidation of methane. Oxygen carrier La_1−xSr_xFe_0.8Al_0.2O₃ was prepared by sol-gel method for chemical looping dry reforming of methane. The reaction performance of the oxygen carrier doped with different proportions of Sr was tested by thermogravimetric and fixed bed reactor. The experimental results showed that the oxygen capacity of x=0.4 oxygen carrier in La_1−xSr_xFe_0.8Al_0.2O₃ is as high as 1.88 mmol·g⁻¹, with excellent reaction performance and less carbon deposition. The stability of the oxygen carrier La_0.6Sr_0.4Fe_0.8Al_0.2O₃ was further tested for 20 redox cycles. The oxygen carrier maintained excellent redox performance, achieving 61.2% methane conversion, 97.1% CO selectivity, and 1.81 H₂/CO. The material characterization results displayed that the morphology and crystal structure of the oxygen carrier were stable. The results show that La_0.6Sr_0.4Fe_0.8Al_0.2O₃ is an excellent oxygen carrier suitable for chemical looping dry reforming of methane.

A total three-dimensional method for calculating mixing loss is proposed in response to the problem that existing methods for evaluating the mixing loss of film cooling cannot accurately calculate the loss between mainstream and coolant in the area affected by passage secondary flow. Film holes with different compound angles are set in the suction surface based on HS1A turbine guide vane. Under the premise of meeting the requirements from film cooling effectiveness, the loss mechanism is analyzed by controlling the mainstream and coolant parameters. Different models are analyzed by using the entropy creation of mixing as the basis for loss evaluation. The results show that film cooling effectiveness is improved by arranging subregional compound angle holes with entropy creation of mixing increasing. And the loss of film deflection caused by passage secondary flow and film detachment can be reduced by setting compound angles.

Porous media play a key role in the fields of energy conversion, thermal management and energy storage due to their outstanding advantages such as large surface area, light weight and complex channels. In this paper, the homogeneous porosity and gradient porosity W type porous structure (porosity ε=0.3∼0.5, hydraulic diameter dh=1.33∼3.86 mm) were customized based on the triply periodic minimal surface (TPMS) method. A computational model of combustion reaction of porous media was established to investigate the range of flame stabilization within the W type porous structure and the SC-BCC cubic lattice structure, and the results showed that the flame blowout of W type structure limitation is higher (equivalence ratio φ=0.65 and inlet velocity V_in=0.6∼2.2 m·s⁻¹). Furthermore, the regulation of combustion and heat transfer processes by the W type porous structure with continuous gradient porosity was investigated, and the results showed that the continuous gradient porous structure would significantly broaden the range of flame stabilization by adaptively modifying the heat rejection rate. 

Numerical simulations were performed to study the effects of two factors, water content and contact angle of liquid bridges, on the fluidization characteristics, particle temperature, and particle concentration of wet particles in a spouted bed. Results shows that the number of fluidized particles increases with the height of the spouted bed. When the water content is 0, there is no particle agglomeration in the fountain region; as the water content increases, the flow velocity in the Z direction first increases and then decreases, and obvious agglomeration occurs in the fountain region. When the contact angle increases, the flow velocity in the Z direction first decreases and then increases, with the minimum velocity occurring when the contact angle is 30^◦. The larger the contact angle, the closer the core high-temperature region is to the spout region When the contact angle is 30^◦, the particle movement resistance reaches its maximum value, and the least number of particles escape from the bed layer with the airflow.

The straight-through labyrinth seal is an effective and simple sealing technology which is widely used in gas turbines. In this paper, the effect of rotation on the flow heat transfer in a straight-through labyrinth sealing structure is investigated using numerical simulations. All results are obtained at Reynolds numbers of 6000, 10000 and 15000. The results show that the increase in rotational speed reduces the vortex losses and throttling losses in the labyrinth seal channel, which leads to the increase in the discharge coefficient. With the increase of rotational speed, the heat transfer at the tip wall is enhanced while the heat transfer at the tooth cavity is weakened.

As a typical kind of air pollutants, VOCs has hazardous effects on both the gas environment and impair human health. Abatement of VOCs is the major requirement of China’s ecological environment. Simultaneously, it is the significant content of China’s 14th five-year plan. Based on the economic development, emission source of VOCs in China has been illustrated. Furthermore, the national, local and industrial regulations of VOCs emissions have been concluded. Abatement techniques of VOCs have been compared and analyzed, including various recovery and destruction techniques. The multi-technique characteristics has been concluded. Finally, the dominant abatement technique in the future was prospected.

The application, screening, evaluation and funding of National Natural Science Foundation of China programs in Engineering Thermophysics and Energy Utilization Discipline in 2024 are summarized and statistically analyzed. The strategic research, funding proposals in the field of energy and power under the carbon peaking and carbon neutrality goals are introduced. The outstanding achievements funded by the discipline in 2024 and future work in 2025 are introduced as well.

When the distance between solid walls is only a few nanometers, the self-diffusion coefficient of fluid is reduced by 1∼2 orders of magnitude due to the nanoconfinement effect, and the effects of pore size, temperature and pressure are significant and complex. It is of great significance to analyze the diffusion mechanism and law of nanoconfined fluid and establish a concise correlation formula. In this study, a wide range of fluid self-diffusion coefficient data in nanopores were obtained through careful calculation and analysis based on molecular dynamics simulation. The mechanism and law of the adsorption effect on fluid diffusion were analyzed, and a novel dimensionless diffusion coefficient, which can simply describe the fluid diffusion behavior, was proposed, and the corresponding correlation formula was established with Knudsen number, which has strong applicability.

A finite time thermodynamic (FTT) model of HI decomposition membrane reactor under different sweep modes is established. The sweep flow rate, reaction inlet pressure, permeable membrane thickness and reactor length are taken as decision variables, and the multi-objective optimization is carried out to maximize HI conversion rate, H₂ recovery rate and total entropy generation rate. It is found that the HI conversion rate and H₂ recovery rate are consistent to some extent within a given range of decision variables, but they cannot reach the optimum with the total entropy generation rate at the same time. Compared with the co-current sweep mode, the target values of Pareto front have higher HI conversion and H₂ recovery in the countercurrent sweep mode. Different decision methods are used to select the optimal solution. TOPSIS decision point in cocurrent mode and LINMAP decision point in counter-current mode had smaller deviation factors and could be used as the optimal solution for reactor parameter design.

Turbine design is a complex and empirical process. To reduce the difficulty of turbine design, this paper proposed a method with uniform flux along radial direction and verified its advantages by the design case of E3 turbine. The results show that the through-flow design is basically consistent with the spanwise distribution of parameters in 3D CFD simulation. The results show that the through-flow prediction is basically consistent with the 3D CFD prediction. And the new design scheme reduces the secondary flow loss at the root of the first rotor and increases the turbine efficiency by 0.5%. The proposed through-flow design method has some advantages in improving calculating stability, reducing time cost and improving turbine aerodynamic performance.

Aiming at the problems of low power generation efficiency and water high consumption in the chemical recuperated gas turbine cycle, a chemical reinjection gas turbine cycle is proposed in this study. It is proposed that part of the flue gas in the gas turbine was put back into the reactor, the reaction of methane self-reforming reaction, and realize the efficient transformation of methane conversion rate. The technology features fully realize the effect of improving the quality of gas turbine flue gas waste heat, and improving the circulation work. Through the improvement of the fuel energy conversion process and the optimization of the heat transfer process of the system, the power generation efficiency of the new system is increased by 9.12% compared with that of the chemical heat recovery system at the design point, the performance of the power generation of lower pressure ratio and high gas turbine inlet temperature is better. The economic analysis shows that the economic payback period of the system is 2.3 years, which has good economic benefits. 

District heating system is one of the important carriers for coordinating renewable energy and traditional energy and realizing flexible consumption of renewable energy. Considering the impact of the uncertainty of renewable energy output and user cluster heat load on the dynamic transportation process of district heating network, it is necessary to quantitatively analyze the uncertain variables on both sides of the source and load and the dynamic characteristics of the heating network. This paper first established a dynamic transportation model of the heating network to solve its heat loss and transmission delay characteristics. Secondly, the Gram-Chalier A algorithm was applied to calculate the probability distribution semi-analytical expression of the thermal power of the source and load nodes of the system, and Bayesian credible inference method was used to calculate the fluctuation interval of node thermal power. This paper selected a secondary heating network in Beijing for model accuracy validation and case analysis. The system has 90 nodes and 109 pipes. The results show that the proposed model and algorithm can effectively quantify the fluctuation interval of nodes’ thermal power.

Sintering of CaO grains in CaO/CaCO₃ thermochemical energy storage system results in the decrease of energy storage performance. In this study, the selection criterion of dopants for material modification is proposed based on lattice energy. High-valence metallic oxides via the phase combination modification mechanism can promote cyclic stabilities of energy storage materials. Using calcium citrate as the precursor and doping with 10% mole fraction of TiO₂ make the effective conversion of calcium-based material increase to 0.66 after 15 cycles, which is 2.1 times of raw calcium-based materials. Doping TiO₂ results in the partial formation of CaTiO₃ during the material preparation process. The modified energy storage material has lower apparent size distribution and retain complex pore structures, causing its high anti-sintering performance.

The common bulkhead tank is one of the most efficient forms for cryogenic propellants, and the non-loss storage technology in the common bulkhead tank is the key issue. A cold shield model designed for the common bulkhead and zero boil-off storage is built based on the liquid oxygen/liquid methane common bulkhead and zero boil-off storage experiment system. The coefficient of temperature variation that considers the location variation characteristic and time variation characteristic of the temperature during the transient heat transfer and flow process is proposed. Compared with the nonuniformity coefficient, the coefficient of temperature variation is more suitable to evaluate the thermal uniformity of the cold shield. The results show that the cold shield with coiler performs an excellent thermal uniformity, the maximum coefficient of temperature variation of the cold shield is 3.85%.

As derived concepts of exergy, the application of the energy quality and energy grade has achieved great success in the practice of thermal science in the past four decades. However, their concepts and calculations have not been fully clarified, and the summary and induction of relevant applications are very rare in the literature database. Based on literature research, this paper explains the history of the energy quality and energy grade concepts, and expands the calculations and corresponding environmental reference states in detail. Moreover, analysis methods derived from the concepts and their applications are reviewed to better understand the functions of energy quality and energy grade. In addition, the concept of energy potential in the context of thermodynamics is preliminarily sorted out. This paper clearly shows the history and progress of the research in the field of energy quality, energy grade, and energy potential, which is helpful for the definition and thinking of related concepts.

A novel solar energy driven CO₂ and H₂O conversion system based on concentrated spectral splitting technology is proposed for liquid hydrocarbon fuels production. Concentrated solar energy is divided into two parts via spectral splitting technology. The short-wave segment is utilized for photovoltaic power generation to drive a proton exchange membrane cell, by which water can be split into hydrogen and oxygen, providing the reactant necessary for CO₂ hydrogenation into methanol. While, the rest part is converted into heat and stored in heat transfer fluid, and can be used in gaseous reactant preheating, water desalination, and/or direct air carbon dioxide capture, etc. Then, the mathematical theoretical model for this solar energy utilization system is built, and its thermodynamic performance is also investigated. It is found that the system achieves its maximum exergy efficiency of about 23.40%, when the temperature of gallium arsenide photovoltaic panel reaches 105°C, along with a corresponding hydrogen production rate of 2.84 L·min⁻¹. Employing an exhaust gas recirculation reaction mode, the CO₂ hydrogenation reactor achieves a hydrogen conversion rate of 68.19%, methanol product selectivity of 97.90%, and methanol yield of 0.91 g·min⁻¹, at reaction pressure of 5 MPa, reaction temperature of 220°C, and a carbon-hydrogen ratio of 1 : 3. The solar energy-to-methanol fuel efficiency of this system is calculated to be 6.08%. This work offers a viable approach for leveraging solar energy to mitigate carbon emissions and advance carbon neutralization technologies.

The combination of spray flash and mixing evaporation (FME) was one of the most effective way to process wastewater with high aqueous salt concentration. In this paper, aqueous NaCl solution was selected as working fluid, a comprehensive calculation model for flow field of FME, including movement, evaporation and crystallization of droplets was set up on basis of previous experimental results. Numerical simulation was carried out with initial diameter of droplets between 20 and 200 μm, initial temperature between 100 and 120^◦C, initial mass fraction 0.26, initial speed 20 m·s⁻¹, and air speed of 15 m·s⁻¹, air temperature between 100∼300^◦C, and spray angle between 0 and 90◦ Results suggested that, during the increasing of spray angle from 0 to 90^◦ the main location of crystallization moved from the spray axis to the top of spray plume, making crystal easier to be separated. Besides, with the increasing of spray angle, superheat or air temperature, the average mass fraction of crystallization increased, but the average particle size and crystallization distance decreased. In order to measure the effect of crystallization within a given distance, complete crystallization efficiency was defined as the ratio of the mass flow rate of crystal salt of droplets to the mass flow rate of dissolved salt at inlet. Results suggested that this efficiency could be improved by increasing spray angle or air temperature. A semi-empirical formula for complete crystallization efficiency was proposed, and the main error between its calculated value and the simulated value was in ±35%. Above conclusion could provide technical support for design and operation of industrial desalination system.

Lignite has high water content, and direct combustion for power generation has problems of low efficiency and large investment. Therefore, the supercritical CO₂ lignite-fired power plant integrated with heat pump pre-drying lignite system was proposed and the coupled thermal calculation model for the supercritical CO₂ cycle power cycle and the heat pump pre-drying lignite system was established. Then, based on a 660 MW reference case, the thermodynamic performances of the lignite direct-fired power plant and the lignite-fired power plant integrated with heat pump pre-drying lignite system are compared and analyzed. The results show that the heat pump drying lignite can increase the power generation efficiency by 1.44% and reduce the standard coal consumption rate by 8.06 g·(kWh)⁻¹. Exergy analysis further shows that the heat pump drying lignite reduces the combustion exergy destruction of the boiler by 4.47% and reduces the boiler exhaust exergy loss by 0.35%, which increases the exergy efficiency of the lignite-fired power plant by 1.29%. And the reduction of the combustion exergy destruction and exhaust exergy loss of the boiler is the main reason for the energy saving. The higher the dryer efficiency, the lignite drying degree and heat pump COP, the more significant the energy-saving effect of the lignite-fired power plant integrated with heat pump pre-drying lignite system.

This article compares three compression heat recovery methods in cryogenic air separation units, including organic Rankine–assisted air compression heat recovery system (ORC-ACS), organic Rankine-electric driven vapor compression-assisted air compression heat recovery system (ORC-VCR-ACS), and organic Rankine-driven vapor compression-assisted compression heat recovery system (ORVC-ACS). The calculation results show that among the three systems, ORVC-ACS has the highest energy-saving efficiency and economic benefits. The system power saving of the ORVC-ACS reaches 1723 kW, which saves about 15.1 GWh/a, and accounts for about 7.5% of the total power consumption of the air compressors. The ORC-VCR-ACS is next, and its highest system power saving is about 1268 kW, while that of the ORC-ACS is about 783 kW. The highest net present value of the ORVC-ACS is about 99 million CNY with an average discounted payback time less than 3.2 years, indicating great economic potential. The maximum net present value and maximum discounted payback periods of the ORC-VCR-ACS and ORC-ACS are approximately 97 million CNY, 3.4 years and 44 million CNY, 4.8 years, respectively. In addition, the three systems can save up to 10.8, 11.4 and 4.9 million tons of carbon dioxide per year, respectively, showing great environmental potential.

In the context of stringent requirement of energy efficiency and carbon neutrality, there is an increasingly growing trend towards utilizing hydrogen-enriched fuel blend in gas turbine combustors. Introducing hydrogen into the fuel increases the likelihood of flashback during the combustion process, necessitating further investigation. In this research, a premixed bluff-body swirl burner is chosen for numerical simulations of premixed CH₄/H₂-air combustion under various operating conditions. SST k-ω model and Flamelet-Generated Manifold (FGM) are employed, respectively, as turbulence model and combustion model. It is of great importance to accurately calculate turbulent flame speed in order to predict flashback, thus the influence of hydrogen addition on laminar flame speed is considered when concerning the closure of the source term in the transport equation of the progress variable. The numerical models and methods are verified against experimental results. Furthermore, the impact of three key factors on flashback are investigated: equivalence ratio, inlet velocity, and thermal boundary conditions of the bluff-body. Results show that decreasing equivalence ratio, increasing inlet velocity, or enhancing heat transfer at bluff-body wall tend to alleviate the occurrence of flashback. The obtained insight could potentially guide the design of hydrogen-enriched gas turbine combustor.

As a very competitive hydrogen production technology, supercritical water gasification technology has important practical significance in exploring the combination route with hydrogen metallurgy under the carbon peaking and carbon neutrality goals. This paper presents a new technical route: supercritical water gasification for hydrogen production complete the reduction of metal oxide in the same reactor, simplifying the metallurgical process. The feasibility of this technology was verified from thermodynamics and experiment, and Fe₃O₄, Cu and MoO₂ powders could be rapidly reduced to prepare in supercritical water gasification atmosphere with glycerol mass concentrations of 2%, 5% and 10%. Under the mass concentration of formic acid 50%, 60% and 70%, blue tungsten oxide and violet tungsten oxide can be obtained by supercritical water gasification, and under the mass concentration of 50%, the needle structure is obvious and less agglomeration.

To study the effect of the total temperature distortion on the performance and stability of the axial transonic compressor, three-dimensional full annular numerical simulations are performed to compare the frozen rotor steady and transient rotor stator unsteady calculations for the different circumferential ranges of the total temperature distortion region using the Darmstadt transonic compressors as the object of study in this paper. The results show that the larger the circumferential range of the high-temperature distortion region at the design speed conditions, that the difference between the frozen rotor and the transient rotor stator calculations is smaller; According to the analysis of the internal flow field at the highest efficiency point under the 180^◦ circumferential distortion region, it is found that the frozen rotor could not capture accurately the circumferential deflection of the distortion region; The closer the predictions obtained from parallel compressor theory for the near stall operating points to the unsteady calculations when the circumferential range of the high temperature distortion region is larger; The orbit method analysis reveals that the larger the circumferential range of the distortion region, the greater variation in the operating characteristics of the rotor blades at different circumferential positions. 

In this paper, the in-cylinder combustion characteristics and soot generation processes of diesel blended polyoxymethylene dimethyl ethers (PODE)-gasoline reactivity controlled compression ignition mode were collected on an optical engine using high-speed imaging combined with the two-color method for clean fuel replacement of conventional fuels. The test results showed that with the increase of PODE blending at the same premixed ratio, the peak in-cylinder pressure, peak heat release rate and pressure rise rate were reduced, the ignition delay was extended, the combustion duration was increased, the combustion phase was shifted back and the combustion tended to be gentle. At a premixed ratio of 50%, the total heat release volumes from single cycle combustion of P20D80 and P50D50 as direct injection fuel was 97.89% and 95.39% of the total heat release volumes of direct injected P0D100 respectively, the total soot generation from single cycle was 55.22% and 36.55% of that of direct injected P0D100 respectively, the high temperature region of soot was reduced by 52.9% and 73.32% respectively, the stable values of soot average temperature were reduced by 6.65 K and 20.25 K respectively, and the stable values of soot average KL factors were reduced by 10.35% and 16.12% respectively. In contrast, P50D50 as a direct injection fuel ensured high combustion thermal efficiency and effectively suppress the generation of soot.

For the future deep space exploration missions, the cryogenic propulsion engines are expected to be restartable for many times. Prior to the restart, it is required to chill down the cryogenic propellent pipe lines. Besides, for the cryogenic upper stage, there are also demands for the pipe line chilldown. During the chilling down, there is a tradeoff between the propellent consumption and chilldown time. It is necessary to conduct studies on the chilldown characteristics for cryogenic propellent pipe lines and find efficient chilldown methods. However, due to pipe wall temperature variation and cryogenic fluid properties, the chilldown process involves some complicated two-phase flow heat transfer mechanisms such as nucleate boiling, inverted annular film boiling and so on. In addition, some abnormal phenomenons are also observed. For example, transition boiling occurs at around 120 K of wall temperature, which is much higher than steady state experimental results. In this study, a transient numerical model for the liquid hydrogen pipe line chilldown is established. Different heat transfer mechanisms during the chilldown are analyzed. The numerical model simulation results agree well with NASA’s hydrogen experimental data. According to the model, during the chilldow, large amount of liquid hydrogen is entrained by the vapor, which results in the waste of liquid hydrogen. The wall heat flux comes to its local maximum valve when the wall temperature is around 115 K. This is due to the coexistence of nucleate boiling and film boiling. Thermal resistance of the pipe wall cannot be neglected. Using the multistep flow rate chilldown scheme can reduce the liquid hydrogen consumption of 27.5% without evident chilldown time extension.

Abundant compounds with different pore sizes can be synthesized by modifying the metal sites or inorganic pillar of SIFSIX hybrid microporous materials. The mechanism of the influence of TiF₆²⁻ pillar on the adsorption and charge transfer characteristic of CO₂ is still insufficient. In this study, TiFSIX-3-Fe and TiFSIX-3-Zn are constructed. The adsorption and charge transfer characteristic of CO₂ in TiFSIX-3-M (TiF₆²⁻ is inorganic pillar, M is metal atom, M=Fe, Zn) was analyzed by DFT. The adsorption sites, band structure, adsorption energy, differential charge density and Mulliken population were calculated. The results show that: TiFSIX-3-Fe and TiFSIX-3-Zn have smaller volume, which are 1.26×10⁵ nm³ and 1.22×10⁵ nm³, respectively. The addition of Ti atom improves the band structure of TiFSIX-3-Fe but has no obvious effect on Zn. The highest adsorption energy are −0.356 eV and −0.361 eV, respectively. The differential charge density indicates that the charge transfer is from F atom to C atom. Combined with Mulliken’s population, the charge transfer quantity increases to 0.072 e and 0.075 e, respectively. The addition of Ti atom significantly enhances the adsorption energy and charge transfer characteristics of CO₂. This work provides theoretical guidance for the development of hybrid microporous materials with high CO₂ capture characteristic.

Northwest China is rich in wind and solar resources. In recent years, new energy power generation systems based on wind and solar energy have developed rapidly. However, due to factors such as productivity and population, the power load in the northwest region is low. In addition, the random and unstable characteristics of new energy make the phenomenon of abandoning wind power and solar power seriously. Some studies have shown that the problem can be solved by new energy power generation coupled with hydrogen energy storage technology. Based on the concept and structure of the combined wind and solar power generation coupled electrolyzer hydrogen production system, take the hydrogen production rate of system as the research goal. A hybrid architecture model including wind power generation system, photovoltaic power generation system and electrolyzer was established under the MATLAB/Simulink. Through this model, the characteristics and laws of hydrogen production in the natural environment of the northwest region were explored.

In order to accurately understand the thermodynamic characteristics of compressed air in a storage cavern of compressed air energy storage(CAES) plants during the first operation cycle, a numerical model for compressed air thermodynamics in storage cavern was eastablished by considering the feature of conjugate heat transfer between fluid and surrounding rock, and the feature of non-isothermal flow. The model was also validated by Huntorf CAES plant data. Based on the comparison analysis of the thermodynamic process of compressed air in the different shape of caverns with the same volume and of the spatial distribution charateristics of temperature, two ways for temperature control of compressed air in the cavern were proposed, and their effects were also explored by numerical method. The research results show that the average temperatures of compressed air in caverns with different profile shape are significantly different, even if the opreation pressures are the same. The cavern shapes also have a significant influence on the spatial distribution characteristics of compressed air temperature. Local high temperature phenomenon occurs in the tunnel type of cavern. For the tunnel type of cavern, the temperature control effect of the measurement such as extending the initial charging time is poor, while the way of reasonable arrangement of the air inlet and outlet positions has a good effect on temperature control.