30 May 2025, Volume 46 Issue 6
    

  • Select all
    |
  • WUSIMAN Kuerbanjiang, SHI Lin, DAI Xiaoye
    Journal of Engineering Thermophysics. 2025, 46(6): 1729-1737.
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    Sugar alcohols are gaining attention as highly promising medium-temperature phase change materials (PCMs) due to their superior overall performance. Currently, the evaluation of the comprehensive performance of sugar alcohol composite PCMs optimized for a single objective remains incomplete. In this study, erythritol (Ery) was used as a PCM, combined with expanded graphite (EG) and Al2O3 as additives. A ternary composite PCM was prepared using the melt blending method, and its comprehensive performance was assessed. The results showed that the components of the ternary composite PCM exhibited excellent physical compatibility. Among them, the ternary composite PCM with 0.5%(wt) Al2O3 / 1.5%(wt) EG/Ery demonstrated the best overall performance. This ternary composite PCM maintained a high melting enthalpy of 312.9 kJ/kg and, compared to Ery, achieved a thermal conductivity of 1.083 W/(m·K), an increase of 54%. In the isothermal cooling test, its supercooling degree was reduced by 9.3°C, in the non-isothermal cooling process, its supercooling degree was reduced by 16.7°C. The pyrolysis temperature increased by 20°C, and the maximum pyrolysis rate temperature increased by 8.3°C, significantly enhancing thermal stability. Additionally, after 40 isothermal cycling tests, its supercooling degree and melting enthalpy remained essentially unchanged, demonstrating excellent cycling stability.
  • SUN Deping, QIAO Guangchao, SHI Feixiong, LI Yiran, ZHANG Dazhi, FENG Xing, LIU Changxin
    Journal of Engineering Thermophysics. 2025, 46(6): 1738-1446.
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    The International Maritime Organization (IMO) and port States have increasingly stringent requirements for ship energy efficiency, and the improvement of ship energy efficiency based on waste heat recovery is one of the most effective ways to meet this challenge. TEG-ORC combined cycle is a new method to realize the utilization of multiple waste heat steps in ships, but the effect of bottom cycle ratio on the performance of combined cycle system has not been studied in detail. The theoretical model is optimized and the influence of variable bottom cycle ratio on the performance of the main parameters of the combined cycle system is studied experimentally. The experimental results show that under the conditions of ORC bottom cycle working medium R245fa, working medium mass flow rate m = 0.079 kg/s and evaporation pressure P = 0.7 MPa, with the gradual increase of bottom cycle ratio, the system output power and flue gas waste heat utilization rate of the main engine increase, and the cost of combined power generation decreases. When the bottom cycle ratio of TEG/ORC is 0.885, the flue gas waste heat utilization rate of the main engine is 85.07%, the total output power of the system is 688.4 W, the thermal efficiency of the system is 7.07%, and the power generation cost of the combined cycle system is 3.338 CNY/kWh.
  • ZHENG Zhimin, SHEN Haoqi, FENG Shuaishuai, LU Lehao, WU Haibo
    Journal of Engineering Thermophysics. 2025, 46(6): 1747-1759.
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    The efficiency of the supercritical carbon dioxide (SCO2) Brayton cycle is closely related to the parameters of the components in the cycle system. Based on the MATLAB/Simulink platform, a modularized calculation procedure for the SCO2 Brayton cycle was established. The design and analysis of a large capacity coal-fired power generation system were carried out using split flow recompression and secondary reheating schemes, and the influence of relevant parameters of components was investigated, such as turbines, compressors, and flue gas coolers. The results showed that the optimal split ratio of the recompressor decreased with the increase of turbine inlet pressure; The power distributions of the heater exchangers in boiler were significantly affected by the split ratio of the flue gas cooler; The circulation efficiency of the system reached its highest level at 35 MPa and 600°C , when the primary reheat pressure and secondary reheat pressure were 25 MPa and 16 MPa, respectively, and the split ratio of the recompressor and the flue gas cooler is 0.3 and 0.05, respectively; As the inlet temperature and pressure of the compressor increased, the circulation efficiency gradually decreased. The research results of this article have important reference value for the establishment and optimization of the large-scale coal-fired SCO2 power generation systems.
  • XU Xin, WU Chuang, LIU Chao
    Journal of Engineering Thermophysics. 2025, 46(6): 1760-1771.
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    Owing to their favorable thermodynamic performances, transcritical power cycles utilizing carbon dioxide (CO2)-based mixtures have emerged as an optimal choice for the power cycle part in medium to low-temperature trough-type concentrating solar power (CSP) systems. This study introduces a two-layer decision-making framework grounded in multi-objective optimization to facilitate the selection of appropriate CO2-based mixtures. The analysis assesses the effects of diverse key parameters on the thermodynamic performance, economic feasibility, and environmental dimensions of the system. The research findings reveal that the CO2/R32 mixed working fluid maximizes the net power output, achieving 253.77 kW. Simultaneously, the CO2/R161 mixed working fluid exhibits superior thermal efficiency (9.80%), economic viability (0.6444 USD/kWh), and reduction in carbon emissions (0.0329 kg/kWh). Moreover, TOPSIS decision-making results suggest that, across various thermal storage tank temperatures, the selection of the CO2/R161 mixed working fluid is of paramount importance.
  • ZHAN Junnan, LIU Taixiu, LI Zhulian, GAO Shuo, LIU Qibin
    Journal of Engineering Thermophysics. 2025, 46(6): 1772-1779.
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    To address the challenges of carbon dioxide capture in hydrogen production and electricity generation from fossil fuels, as well as the supply-demand mismatch in the utilization of intermittent and unstable solar energy, this paper proposes a solar-driven hydrogen and electricity cogeneration system with source decarbonization. The system employs an iron-nickel mixed oxygen carrier, which reduces the reduction reaction temperature while achieving efficient carbon dioxide capture. The results indicate that under a solar irradiance of 750 W/m2, the system achieves a solar energy input ratio of 30.80%, with energy utilization rates of 67.23% and 59.47% for hydrogen production and electricity generation under design conditions, respectively. In addition, by regulating the sensible heat of the inlet oxygen carrier and its oxidation degree during the reaction process, the system can operate stably under varying irradiation conditions and flexibly adjust hydrogen and electricity outputs. When the solar irradiation intensity is 750 W/m2, the adjustable hydrogen-toelectricity ratio ranges from 0.17 to 0.59 m3/kWh, demonstrating excellent adaptability to varying operating conditions. This study offers a novel technological pathway for highly efficient hydrogen and electricity co-production and source decarbonization through the complementary use of solar energy and natural gas.
  • QIAN Qian, LUO Zhongyang, SHI Jingkang, WANG Yuanlin, WEI Qi, SHI Yunmei
    Journal of Engineering Thermophysics. 2025, 46(6): 1780-1790.
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    Efficient use of lignocellulose is essential for reducing carbon emissions and dependence on fossil fuels. While advances have been made in (hemi)cellulose conversion, improving its economic feasibility remains a pressing need. Valorization of lignin is essential for achieving profitability in biorefineries and remains a major technical challenge. This study proposes innovative cross-coupling routes that integrate hydrothermal-pyrolysis and hydrothermal-biochemical conversion to efficiently separate the primary components of lignocellulose, converting lignin and its derivatives into a relatively narrow product slate. Following the lignin-first strategy, two poplar graded utilization schemes were developed, modeled in Aspen Plus, and evaluated. While both systems demonstrated comparable carbon efficiency, there is a notable difference in their energy recovery. Techno-economic analysis indicates that the minimum fuel selling prices for the two schemes are 7456 CNY/t and 6990 CNY/t, respectively, which are comparable to domestic gasoline market prices. Future technological breakthroughs to effectively enhance ideal product yield are key to improving the economic viability of the system. Additionally, decreasing fixed capital investment by reducing reactor size and/or pressure, as well as controlling major consumable costs is expected to further enhance the product’s market competitiveness.
  • LI Xinlu, DUAN Liqiang, WANG Xiaomeng, YANG Chaoyun, WANG Qiushi
    Journal of Engineering Thermophysics. 2025, 46(6): 1791-1800.
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    A combined cooling, heating and power(CCHP) system with solid oxide fuel cell (SOFC) coupled with chemical chain combustion (CLC) carbon capture is proposed. After methanol reforming, H2 is separated by membrane separation technology to drive SOFC to generate electricity, and the remaining syngas is burned into the chemical chain to achieve carbon capture, and supercritical carbon dioxide (SCO2) Brayton cycle, absorption refrigeration system (ARS) and heating system are integrated to realize cascade utilization of energy. The performance analysis models of key components and the whole CCHP system are constructed, and the effects of key parameters such as reforming reactor temperature, water-carbon ratio, SOFC operating pressure and temperature on the system performance are studied. A combination of energy analysis and analysis was used to comprehensively evaluate the performance of the system. The energy efficiency and exergy efficiency of the new system reach 66.29% and 57.70%. Compared with the reference system of traditional CCHP system with SOFC without integrated with chemical chain combustion, after achieving a CO2 capture rate of more than 99%, the power generation efficiency of new system is about 1.95% higher than that of the reference system and the exergy efficiency of new system has an increase of 0.75%.
  • LUO Jing, SUN Yanlei, CHEN Yanyan, ZHANG Limin, HU Jianying, LUO Ercang
    Journal of Engineering Thermophysics. 2025, 46(6): 1801-1811.
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    As the key component of free-piston Stirling generator to achieve the acoustic-to-electrical conversion, it is of great significance to carry out independent and detailed testing on the linear alternator. In this paper, a system model for dynamic testing of a high-temperature and high-power linear alternator is established based on SAGE software, which adopts an opposed compressor as the driving source and introduces acoustic impedance into the model to efficiently match the compressor with the alternator to be tested, so as to investigate its dynamic characteristics under different operating conditions. The results show that under the rated operating conditions of mean pressure of 15 MPa and operating frequency of 70 Hz, the total power consumption of the compressor has to reach 26.8 kW in order to achieve 18 kW of electric power output from the alternator to be tested, with 90.5% of the compressor’s electric-to-acoustic efficiency, 86.4% of the acoustic power transfer efficiency in the thermal buffer tube, and 85.9% of the acoustic-to-electric efficiency of the alternator to be tested. By comparing with the experimental results, the accuracy of the calculation model is verified. This paper has a certain guiding significance for the design work of the dynamic testing system for linear alternators.
  • CAI Jiahao, SU Bosheng, LI Liang, YUAN Shuo, GONG Xiaofeng
    Journal of Engineering Thermophysics. 2025, 46(6): 1812-1824.
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    Due to the non-uniform concentration characteristics, the thermochemical reactor driven by the dish solar collector is prone to the problem of material life attenuation and even structural damage. The ‘two-dimensional’ concentrating design idea of the heating surface of the traditional reactor is extended to the ‘three-dimensional’ concentrating design idea, which can dynamically adjust the concentration ratio of the heating surface to promote the thermal matching of the thermochemical reaction, thus having great potential to reduce the irreversible loss of solar thermochemical conversion. In this paper, a spiral flow channel solar hydrogen production reactor based on inverse cone structure is proposed. By designing the cone angle to regulate the concentration ratio, the reasonable matching between the heating temperature of the heating surface of the flow channel and the heating temperature required for the thermochemical reaction process is realized. Compared with the ‘two-dimensional’ heating surface concentration design, the solar thermochemical efficiency is significantly improved by 38.2%.
  • TANG Jingchun, DAI Zikun, LI Jing, LI Pengcheng, MA Zhenyu, NIE Sifan, JIE Desuan, LIU Xunfen
    Journal of Engineering Thermophysics. 2025, 46(6): 1825-1834.
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    Supercritical carbon dioxide (sCO2) cycles have a low expansion ratio and high turbine outlet temperature and need recuperators. The thermodynamic irreversibility in the recuperators is large due to the mismatch between the high-pressure and low-pressure sCO2. The power generation efficiency of the sCO2 cycle can be improved by coupling an organic Rankine cycle (ORC) at its bottom and constituting a cascade system. However, current research focuses on conventional organic working fluids with critical temperatures <200°C, and the unidirectional heat transfer from sCO2 to ORC results in a low overall efficiency. This paper proposes a novel cascade system coupling sCOwith an ORC using a high-temperature working fluid biphenyl-diphenyl oxide (its evaporation temperature is 400°C). The high-temperature sCO2 from the turbine outlet is used to drive the ORC. Meanwhile, the heat released from the ORC is used to preheat sCO2, thereby increasing the sCO2 preheating temperature and the thermal efficiency. This paper proposes four coupling configurations. Thermodynamic analysis shows that the first type has the highest efficiency. The maximum efficiency of 45.45% is achieved at the sCO2 turbine inlet temperature of 650◦C and inlet pressure of 20 MPa. It is higher than the efficiency of existing sCO2-ORC cascade system at the same operating conditions. In addition, the optimized exergy efficiency can reach a maximum of 72.76%.
  • WANG Xinyu, MA Teng, LI Zhuoling, LIU Yang, YU Bo, CHEN Yumin, ZHAI Ming, ZHOU Huaichun
    Journal of Engineering Thermophysics. 2025, 46(6): 1835-1843.
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    Biomass, as a renewable carbon-neutral energy source with abundant reserves, is a good choice for future energy supply. Its coupled combustion with NH3 can effectively reduce CO2 emissions, which is the main way to realize emission reduction. In this paper, the combustion characteristics of the co-combustion of corn straw char/NH3 were investigated and the conversion mechanism of elemental N was elucidated by reactive force field molecular dynamics (ReaxFF-MD) simulations. The effects of temperature, oxygen equivalence ratio (Ω) and NH3 mixing ratio (CR) on the products were investigated. The results showed that HNOx· played a decisive role in the production of H2O and NOx. Most of the free radicals such as HNO·, HNO2·and CH2O· react with OH· to produce H2O. the highest frequency of the reaction is HNO· + OH· →H2O + NO with 35 times. This reaction is also the main reaction to generate NO2. Among the CO2 generation reactions CNO3· →CO2 + NO had the highest frequency of 23 reactions. This study provides a reference for subsequent pollutant emissions from biomass power generation and its utilization.
  • YU Jiazheng, QIU Jiahui, YANG Chen, ZHANG Min, DU Juan, ZHANG Lei
    Journal of Engineering Thermophysics. 2025, 46(6): 1844-1857.
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    In order to quickly evaluate the characteristics of axial flow compressors, a one-dimensional prediction approach based on the mean-line method was developed, and an adaptive calibration method was developed for the empirical model to improve the prediction accuracy. Based on these, the off-design characteristics of three compressors were analyzed and verified. Meanwhile, characteristics of a 1.5 stage transonic compressor with various angles of the inlet guide vane was analyzed. Results are compared with 3D numerical simulations and experiments, and they further validates the prediction accuracy of the developed prediction method under different conditions of rotating speeds and inlet guide vane angles. Finally, based on the above one-dimensional analysis method, the NSGA II algorithm is adopted to carry out an optimization 0f inlet guide vane regulation law for the 1.5-stage compressor at different rotating speeds. Also, three-dimensional numerical simulations are performed to validate the optimized results. The results show that the developed coupling method is expected to provide an efficient and accurate prediction method for rapid screening of variable guide vane/stator angle schemes of multistage axial compressor.
  • QIN Kai, HUANG Diangui
    Journal of Engineering Thermophysics. 2025, 46(6): 1858-1864.
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    In order to eliminate the influence of the discontinuous curvature on the performance of the compressor airfoil and improve the optimal design capability of compressor cascade, a parameterization method with a continuous-curvature surface was proposed. The curvature continuity condition of the two Bezier curves at the connection point can be achieved, by using the connection of the endpoint guide vector and the control point of the Bezier curve to make the corresponding control point satisfy a specific relationship. Based on this property of the Bezier curve, the leading edges and trailing edges of the compressor blade are reconfigured. In order to remove discontinuous curvature at the connection point between the leading edge and the airfoil,the least squares method was used to fit the suction section and pressure section of compressor airfoil in the parameterization process. The numerical simulation method is used to calculate the flow field, and the aerodynamic performance of the parameterization airfoil and the original airfoil are compared. The results show that the parameterization method can obtain a continuous-curvature surface, which can effectively improve the flow condition of the leading edge and improve the aerodynamic performance of the compressor airfoil.
  • LI Liying, WU Hong, QUE Xiaobin, QIU Ying, HE Liu
    Journal of Engineering Thermophysics. 2025, 46(6): 1865-1874.
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    In order to reduce the total pressure loss of inlet manifold and inlet casing, and decrease compressor inlet distortion in heavy-duty gas turbine, the effect of circumferentially non-uniform design of compressor inlet casing struts was researched, and the design concept of non-consistent strut profile has been proposed and studied. Different design of the strut profiles are utilized at different circumferential positions, which can reduce the impact of strut wake and improve the compressor inlet non-uniformity. The analysis indicates that the circumferentially non-uniform and non-consistent profile of struts can reduce the total pressure loss of compressor inlet casing and decrease the total pressure and swirl distortion of compressor inlet.
  • SHI Xuyang, WU Yanhui, LI Ziliang, LI Haohua
    Journal of Engineering Thermophysics. 2025, 46(6): 1875-1882.
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    A representative transonic turbine stage TTM was selected as the research object to examine the effects of Circumferential Multiple Grooves (CMG) design parameters on the efficacy of tip leakage flow control. A validated high-fidelity numerical methodology was employed to analyze the impact of key design parameters-including the width of the shaping area (WSA), the number of circumferential grooves (N), and CMG distribution ratio (P)-on the control efficacy of tip leakage flow in the TTM turbine stage. The findings demonstrate that the CMG’s flow control effectiveness exhibits a characteristic trend of initial decline followed by an increase, culminating in gradual stabilization with WSA increases when P maintained constant. For CMG with fixed groove number N, the control effectiveness under varying P shows an initial decrease, followed by an increase, then a slight decrease until stabilization with increasing shaping width WSA. Flow field visualization reveals that CMG implementation induces substantial complexity in turbine blade tip clearance flow structures. Since the leakage flow through the CMG (CMG-TLF) is in the opposite direction to the original leakage flow (TLF), the introduction of a CMG with appropriate parameters can make the net leakage flow rate across the rotor blade tip region significantly reduced. This configuration effectively minimizes tip leakage flow magnitude, resulting in enhanced turbine stage efficiency. When the design parameters of the CMG are set to P0.8WSA10N2 the flow rate of the rotor tip leakage flow is reduced by 89.76%, and the turbine stage efficiency is increased by 0.53%.
  • ZHANG Lidong, TIE Hao, LIU Huiwen, TIAN Wenxin, ZHAO Xiuyong, CHANG Zihan, LI Qinwei
    Journal of Engineering Thermophysics. 2025, 46(6): 1883-1890.
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    The staggered height layout has a great influence on the overall wake structure of wind farms. In order to study the effect of small wind turbine operation on the wake in a staggered-height wind farm, a wind tunnel was used to study the evolution of the wake in the smooth flow by small wind turbines operating at different location in downstream of the large wind turbine downstream. The results show that the operation of a small wind turbine exacerbates the cleavage of the vortex structure inside the wake and reduces the contribution of the vortex turbulence kinetic energy to the whole field energy to 1/4 of that in the case of no small wind turbine operation, and the extent of the effect in the vertical direction diminishes with the increase of the height, with the maximum extent of the effect being no higher than that of the upper tip position of the large wind turbine. In the horizontal direction, the effect of small wind turbine operation only occurs within two times the turbine diameter downstream of the turbine, but can still be observed farther upstream. 
  • FENG Yuhao, WANG Xiaodong, WANG Liwen, ZENG Tao, XU Na
    Journal of Engineering Thermophysics. 2025, 46(6): 1891-1897.
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    As wind turbine technology advances, research on pre-bent blades is increasing globally. However, the complex geometry of pre-bent blades also increases the error in aerodynamic performance analysis, posing challenges to theoretical research. This paper conducts research on the aerodynamic characteristics of pre-bent blades using BEM calculations. Based on the existing 3D geometric model of pre-bent blades, a self-developed program determines the pre-bending curve of the blade. This curve is then incorporated into the traditional BEM method, and by modifying the existing BEM theory, it alters the inflow angle and the force situation on each blade element to obtain more accurate aerodynamic performance of pre-bent blades. Moreover, for different wind speed ranges, the tip loss correction factor is divided into two parts, adopting the Prandtl tip loss model and the Shen tip loss model, to better account for the tip effect. The results show that compared to traditional methods, the error of the modified BEM method is significantly reduced, and its accuracy is improved, providing a theoretical basis for the rapid calculation of the aerodynamic performance of pre-bent blades.
  • YUAN Jiaqi, YE Kaijie, WANG Yijin, HE Denghui, YANG Lina, CAO Ming, BAI Bofeng
    Journal of Engineering Thermophysics. 2025, 46(6): 1898-1905.
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    To investigate the effect of inlet gas content of the return pump in the aviation engine lubricating oil system on its working characteristics, the vibration signals of the rotary vane oil pump under different inlet gas fraction conditions were measured to obtain the influence of changes in inlet gas fraction on the vibration characteristics of the rotary vane pump. The vibration signals were processed and analyzed based on the Complete Ensemble Empirical Mode Decomposition (CEEMD) method. The results indicate that an increase in inlet gas fraction will weaken the vibration intensity of the rotary vane pump, and the horizontal and vertical vibration intensities are significantly higher than that in the axial direction. The horizontal vibration of the front bearing seat at low gas fraction is mainly induced by pressure pulsation during the operation of the rotary vane pump, and the dominant factor of vibration at high gas fraction is mechanical factors. The sample entropy value of the vibration acceleration signal decreases with the increase of the inlet gas fraction, showing a good regularity, and can be used for predicting the inlet gas volume fraction of the rotary vane pump and fault diagnosis.
  • JIANG Lijie, ZHANG Renhui, CHEN Xuebing, GUO Guangqiang
    Journal of Engineering Thermophysics. 2025, 46(6): 1906-1911.
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    In order to reveal the complex unsteady flow characteristics of liquid ring pump ejector, flow visualization and pressure testing techniques were used to reveal the spatial evolution and pressure pulsation characteristics of the high-speed jet in the liquid ring pump ejector under different vacuum. The results show that the distribution characteristics and dynamic evolution of the highspeed jet structure in the liquid ring pump ejector can be described by glycerol atomization tracer particles. The shear layer between high-speed jet and low speed fluid has obvious fold characteristics because of enrolling low speed secondary fluid, and the high-speed jet has obvious fluctuation phenomenon under the influence of low speed fluid. The time domain variation of the pressure signal of each measuring point in the ejector has unstable fluctuation under different vacuum. Due to the influence of rotor-stator interaction of the liquid ring pump, the flow in the liquid ring pump ejector presents the characteristics of multi-frequency coupling, and the dominant characteristic frequencies are axis frequency fn, blade frequency fBPF and high harmonic frequency.
  • TAN Jie, LOU Jia, XIA Lei, JIANG Dongyue
    Journal of Engineering Thermophysics. 2025, 46(6): 1912-1914.
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    The achievement of distant, directed, and high-speed movement of liquid droplets on solid surfaces holds significant importance in applications such as water collection and energy recovery. Currently, researchers have made substantial progress in magnetic, thermal, and electric fielddriven methods. Among these, the dielectric wetting driving mechanism has garnered widespread attention due to its advantages of fast response, low energy consumption, and simple device setup. However, current electric wetting driving is predominantly conducted on hydrophobic surfaces, imposing certain limitations on speed. In recent years, the excellent performance of superhydrophobic surfaces has been extensively studied. Yet, there is limited literature on the electric wetting-driven motion of liquid droplets on superhydrophobic surfaces. This paper proposes a method for preparing a superhydrophobic surface and further achieves ultra-high-speed electric wetting-driven droplet motion on this surface. Experimental results indicate stable control over saline water, acids, bases, and certain organic reagents. This is of significant importance for the further expansion of electric wetting-based digital microfluidic technologies.
  • MA Teng, LI Mingjia, XUE Hanwen, NI Jingwei
    Journal of Engineering Thermophysics. 2025, 46(6): 1915-1921.
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    In this paper, the S-CO2 cycle dynamic model considering the thermal inertia of the recuperator is constructed. The flow inertia is characterized by the total channel volume and the fluid volume flow. The heat transfer inertia is characterized by the heat transfer quantity and the thermal conductivity. The dynamic performances of the S-CO2 cycle under mass flow control and temperature control are compared. The changes of S-CO2 cycle dynamic performance parameters such as real-time efficiency and load change rate with the key design parameters of the recuperator are analyzed. The mathematical expression of S-CO2 cycle load change rate, fluid channel volume and thermal conductivity are fitted by multiple regression analysis method. The results show that the mass flow control is more effective. The efficiency loss is 0.72% and the load change rate is 10 %/min. When V=1.5V0 and kA=0.5kA0, the load change rate of S-CO2 cycle is the highest, which is 11.32 %/min.
  • QIAO Lanqing, LAI Qingzhi, TAN Jianyu
    Journal of Engineering Thermophysics. 2025, 46(6): 1922-1930.
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    The integration of chips reduces the heat transfer space, while the pressure drop and pump power of small heat transfer structures are generally large, which is not conducive to energy saving and emission reduction. In order to balance the pressure drop and heat transfer, aiming at the heat dissipation problem of high power single chip in confined space, this paper carried out multiple sets of water cooling experiments on the air channel and the fully filled copper foam jet impinged heat transfer device, and carried out optimization design and experiments on the device filled with copper foam with needle rib structure based on evaluation factors. The results show that when the porosity of copper foam is 96%, the device filled with 20 PPI copper foam at large flow rate has the best comprehensive effect. When the inlet Reynolds number is less than 1492, the comprehensive performance of the device filled with 20 PPI is optimal.
  • WANG Xin, LIU Hongwei, LIANG Bing, YANG Xinle, WANG Fang, SUN Weiji, LI Weizhong, SONG Yongchen
    Journal of Engineering Thermophysics. 2025, 46(6): 1931-1941.
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    Methane hydrate, composed of methane molecules and water molecules, represent a clean energy source. The maintenance of hydrate stability under specific temperature and pressure conditions is referred to as the phase equilibrium of hydrate. Capillary pressure significantly influences the phase equilibrium conditions of hydrate, particularly the equilibrium pressure. This study improves existing phase equilibrium pressure models and proposes a novel hydrate decomposition model, incorporating existing hydrate decomposition models. The proposed model is validated through numerical simulations of experimental studies, demonstrating good consistency between numerical simulations and experimental results, thus verifying the reliability of the proposed model. It distinctly highlights the necessity of considering capillary pressure between hydrate and water when investigating methane hydrate decomposition within micropores. The proposed model further elucidates the evolution of methane generation rate, phase distributions, temperature distributions, and equilibrium pressure distributions. The research findings underscore the importance of incorporating capillary pressure between hydrate and water into hydrate decomposition models, as it significantly affects their accuracy. Specifically, neglecting capillary pressure between hydrate and water leads to underestimated cumulative gas production and methane generation rates. The coupling effect of capillary pressure, hydrate decomposition endothermicity, and inter-phase heat transfer influences the distribution of equilibrium pressure, methane generation rate, phase distributions, temperature distributions, and ultimately, the distribution of equilibrium pressure. Capillary pressure promotes an increase in the phase equilibrium pressure of hydrate, enhances the driving force for hydrate decomposition, and facilitates hydrate decomposition. The proposed model contributes to the theoretical foundation for enhancing hydrate decomposition models.
  • CAI Qingfeng, GENG Wenguang, LIU Fang
    Journal of Engineering Thermophysics. 2025, 46(6): 1942-1946.
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    In order to alleviate the increasingly serious environmental problems caused by livestock and poultry waste, this paper adopts the microwave method to prepare the composite activated carbon of cow dung/steel slag from the perspective of resource utilization of “Treating waste of waste”, and at the same time conducts degradation research on farm wastewater to investigate the degradation rate of chemical oxygen demand (COD) in farm wastewater under different reaction conditions. The results showed that: the appropriate amount of steel slag mixing ratio can improve the COD degradation rate to offset the negative impact of the reduction of the specific surface area of composite carbon; the degradation rate of COD with the increase of microwave power and microwave time will show a gradual increase in the trend; in the optimal reaction conditions, the COD degradation of aquaculture wastewater was up to 55.4%, and the operating cost of treating the unit volume of aquaculture wastewater was about 384.98 CNY.
  • LI Yi, LI Yanxia, LIU Zhongliang, WANG Zhenqiang, MA Weiqi, YU Fengjiao
    Journal of Engineering Thermophysics. 2025, 46(6): 1947-1955.
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    Using a liquid nitrogen refrigeration low-temperature frosting experimental setup, the growth process of initial frost crystals on cryogenic surfaces under natural convection conditions was experimentally studied. The influences of the cold surface temperature, the relative humidity and temperature air on the growth of initial frost crystals was investigated. The experimental results revealed that when the cold surface temperature ranges from −50°C to −100°C, the percentage of frost coverage area increases with the decrease of cold surface temperature, with the increase of air relative humidity and temperature. However, when the cold surface temperature is −100∼−170°C, the percentage of frost coverage area decreases as the cold surface temperature further decreases. In contrast, the relative humidity and temperature of air hardly affect the initial frost crystal growth. It can be seen that the temperature of the cold surface plays a dominant role in the growth of initial frost crystals on the cryogenic cold surface.
  • YU Jiatong, CHEN Zhihao, UTAKA Yoshio, JIN Xuhao
    Journal of Engineering Thermophysics. 2025, 46(6): 1956-1963.
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    With the development of electric vehicles and electronic devices, the two-phase immersion cooling technology using dielectric coolants has gradually attracted attention due to its significant advantages in temperature control performance and integration. The demand for the application of this technology in low-temperature conditions is also increasing. However, there is currently limited research on the boiling heat transfer of dielectric coolants, especially at low temperatures. This study investigates the effect of liquid static pressure on the boiling heat transfer characteristics of Novec 7100 on a smooth copper surface under low saturation pressure conditions using an airtight experimental system with controllable saturation pressure (temperature). The experiment reveals a unique trend of critical heat flux decreasing first and then increasing with decreasing saturation pressure. After the saturation pressure drops below 8.34 kPa, the pool boiling process exhibits “intermittent boiling” and “self-induced subcooled boiling” as two distinct heat transfer modes, and the critical heat flux increases with the liquid height. When the saturation pressure is 0.741 kPa, the critical heat flux values at liquid levels of 23 mm and 35 mm show the greatest difference, approximately 39 kW/m2.
  • LIU Xiaoliang, ZHANG Xuan, ZHANG Long, SONG Mengjie, GAO Yubo, LI Kailiang, WU Longping
    Journal of Engineering Thermophysics. 2025, 46(6): 1964-1968.
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    To explore the effect of Bond number on the freezing characteristics of water droplets, experimental and modeling studies are carried out. The effects of contact angle and cold plate temperature are analyzed. The results show that: 1) When the Bond number is larger than 1, the freezing characteristics of sessile and pendant water droplets are significantly different, and the freezing profile of pendant water droplets is high and the freezing time is long. As the Bond number increases, the difference between them increases. 2) With the increase of contact angle, the freezing profile of water droplets becomes high and the freezing time increases, and the difference of freezing characteristics between sessile and pendant water droplets becomes more obvious. The freezing profile is independent of the cold plate temperature. 3) The correlation formula of freezing time calculation is proposed, and the prediction error is less than 20%. The research has guiding significance for the optimization of anti-icing and de-icing technology.
  • GUO Yangyu
    Journal of Engineering Thermophysics. 2025, 46(6): 1969-1974.
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    It remains an open question in thermophysics to understand and model inelastic phonon scattering in heat transport at solid/solid interface. In this work, we develop a theoretical model to compute and spectrally decompose the energy exchange at the interface based on anharmonic phonon non-equilibrium Green’s function method, to study the interfacial heat transport at Si/Ge interface. For high-frequency phonons, local anharmonic scattering dominates the interfacial energy exchange, while both local and non-local scattering are important for moderate- and low-frequency phonons. The anharmonic decay of interfacial phonon modes plays a crucial role in bridging the bulk phonon modes at both sides of the interface. This work will promote the understanding of interfacial heat transport mechanism and future development of theoretical models.
  • RUAN Rongsheng, LIU Zhengang, CHENG Long, ZHANG Yixuan, LIU Zhenxia, WU Dingyi
    Journal of Engineering Thermophysics. 2025, 46(6): 1975-1990.
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    In order to investigate the deposition characteristics of particles on the surface of plate, a simulating experiment was conducted at room temperature, in which the atomized wax was utilized as the particles and the test model was a flat plate. Six different types of film cooling hole specimens were selected for the experiment, including cylindrical, conical, CONSOLE, RTSH, Laidback fan-shaped, and fan-shaped film cooling holes. The deposition distribution and deposition amount variation on the surface of the different hole specimens were studied. The experimental results showed that the deposition distribution on the pressure surface was concentrated downstream of the film cooling holes. Compared with the cylindrical holes, the five types of air film cooling holes in this study could effectively reduce the deposition of particles on the surface of the flat plate. The deposition rate of the CONSOLE hole decreased by about 69%, and the deposition rate of the other hole types decreased by about 80%. But the cooling efficiency remained relatively constant before and after deposition. Further experiments with cylindrical specimens revealed that the deposition rate tends to increase as the cooling air temperature drops. Moreover, the lower the cooling air temperature, the more pronounced the change in cooling efficiency becomes between the pre- and post-deposition stages.
  • CHEN Yuluo, MA Han, CHEN Qiang, FEI Qingguo
    Journal of Engineering Thermophysics. 2025, 46(6): 1991-2002.
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    Integral thermal protection structures are key components of new reusable high-speed vehicles. Timely and accurate acquisition of three-dimensional temperature field characteristics of thermal protection structures under time-varying aerodynamic thermal loads is an important basis for monitoring the health status of the vehicle. The paper proposes a method for online reconstruction of three-dimensional temperature field of integral thermal protection structure based on spatial-temporal surrogate model. Firstly, the time-varying full-field temperature field dataset of thermal protection structures is doubly downgraded in time and space by means of an proper orthogonal decomposition to reduce the number of online reconstructed model parameters; then, based on the reduced-order temperature field data and the deep neural network method, an efficient surrogate model is established for the inversion of the external thermal loads of the structure and the reconstruction of the full-field temperature field; finally, the time-varying three-dimensional temperature field simulation analysis of the structure under typical working conditions is carried out, and the load inversion and temperature field reconstruction models are linked to achieve the online reconstruction  of the three-dimensional temperature field of the thermal protection structure, which takes the local temperature response of the structure as the input. The results show that the online reconstruction method proposed in this paper is capable of accurately predicting the three-dimensional temperature field of thermal protection structures under time-varying non-uniform loading, based on the temperature response of a small number of local measurement points. 
  • WANG Lu, ZANG Yuanqin, LU Xin, DAI Ziwen, WANG Yibo, WANG Liangbi
    Journal of Engineering Thermophysics. 2025, 46(6): 2003-2011.
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    Train braking process is an important part of train running process and it is an important guarantee of train safety. The rail vehicle converts kinetic energy into thermal energy through the sliding friction between the brake disc and the brake pad during braking process, and it is necessary to adequately cool the brake disc to prevent various braking problems caused by excessive temperature. Therefore, the research on the heat dissipation capacity of brake disc is very important for train safety. This paper introduce a new method to simulate the heat transfer characteristics of brake disc, simplifies the test specimen and device in the heat transfer test of the internal passage of the brake disc, and establishes the motion model of the brake disc rolling at uniform speed on the rolling wall. We simulate the rolling process of the brake disc by the moving grid method, and obtain the flow field distribution around the pin fin brake disc test specimen and the convection heat transfer characteristics in the internal passage of the brake disc test specimen.
  • WANG Jingyu, WANG Lei, SHEN Leilei
    Journal of Engineering Thermophysics. 2025, 46(6): 2012-2020.
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    The operating parameters have significant effect on the performance of the methane reforming membrane reactor driven by the trough solar collector, and the present paper optimized the operating parameters by the response surface methodology (RSM). The multi-physics coupled numerical model of the reactor was established, the species transfer and reaction process inside the reactor was analyzed, and the influence of key parameters on the reactor performance was explored. The central composite face-centered design (CCF) method was used to design numerical experiment cases, and the quadratic regression equations of the methane conversion rate XCH4 , hydrogen yield YH2 , hydrogen recovery ratio RH2 and the energy storage efficiency η were fitted with the variations of the inlet temperature Tin, the inlet mass flow rate qm, the pressure difference Δp between the reaction side and the permeate side, and the pressure in the permeate side pper. All R2 values were higher than 0.97 in RSM models. The contribution rate and interaction influence of each factor on the response values were then obtained through the RSM models. Finally, multi-objective optimization was carried out by taking XCH4 , YH2 , RH2 and η as optimization objectives, and with the maximum temperature inside the reactor not exceeding 873 K as a constraint, the optimal parameters obtained were Tin=673.1 K, qm=0.005 kg·s−1, Δp=0.42 MPa, pper=0.02 MPa.
  • CHANG Zhijuan, LANG Xufeng, LIU Mengyao, WU Xuehong, LÜ Cai, LIU Yong
    Journal of Engineering Thermophysics. 2025, 46(6): 2021-2028.
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    In order to study the effect of the nanoconfinement space formed by graphene on the phase transformation behavior of paraffin, the molecular dynamics (MD) simulation method was used to investigate the confinement effect of 5.1 nm thickness of graphene interlayer on the phase transformation behavior of n-octadecane. The simulation results show that the graphene interlayer can inhibit the thermal motion of n-octadecane molecules, making it more difficult to melt and easier to solidify. Compared with the pure n-octadecane system, the graphene interlayer of this scale increases the solidification temperature and melting temperature of n-octadecane by 51.39 K and 19.33 K, respectively. The solidification latent heat and melting latent heat were decreased by 7.41% and 8.10% respectively compared with the theoretical value. Under the interaction between graphene and n-octadecane, graphene as a two-dimensional plane can provide a template for n-octadecane crystallization, so that it preferentially forms a crystalline layer at the graphene interface. With the decrease of temperature, n-octadecane molecules gradually oriented crystallize outwards from the graphene interface and the graphene interlayer can force the n-octadecane to form a crystalline layer parallel to it.
  • ZHANG Heng, WANG Hu, SUI Bangjie, ZHAN Zhigang
    Journal of Engineering Thermophysics. 2025, 46(6): 2029-2036.
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    To deeply investigate the distribution of liquid water in the proton exchange membrane fuel cell, a two-dimensional model of non-isothermal two-phase flow including gas-water-heatelectricity-force multi-physical fields was developed in this paper to simulate the cell performance and liquid water distribution of two different types of gas diffusion layer (GDL) materials at the operating conditions of 80°C, 200 kPa abs, and 80% RH. Performance tests and neutron radiograph experiments were also conducted and compared with the simulation results. The results show that the cell with Freudenberg GDL have a better performance compared to the cell with Toray GDL, especially at high current density. The peak liquid water saturation using Toray GDL occurs in catalyst layer under the ribs, and it can reach to 0.4. The peak liquid water saturation using Freudenberg GDL occurs in the gas diffusion layer under the ribs, and it can reach to 0.25.
  • HUANG Jingyu, WU Xinzhuang, XIA Shuan
    Journal of Engineering Thermophysics. 2025, 46(6): 2037-2044.
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    In this paper, a typical solver is developed for passive containment water film cooling performance analysis of third-generation nuclear power plants based on open-source computational fluid dynamics software OpenFOAM. The numerical results of water evaporation rate are compared with testing data and calculation error is within ±10%, which verifies the solver is suitable for engineering calculations of passive containment water film cooling performances and complicated thermal-hydraulic conditions. Based on the solver, effects of water film mass flow rate, air flow velocity in the annulus and containment wall temperature on water film cooling performances are simulated separately. Results indicate that simulation error of water evaporation rate is within ±10%; it has little effect on water evaporation rate when changes of water mass flow rate and air flow velocity in the annulus do not exceed ±25%, and changes of containment wall temperature do not exceed 10°C.
  • TANG Qinglong, YAO Mingfa
    Journal of Engineering Thermophysics. 2025, 46(6): 2045-2054.
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    Pre-chamber jet combustion technology is a crucial pathway for achieving breakthroughs in the thermal efficiency of spark-ignition engines. However, there are still many research gaps regarding the combustion mechanisms of pre-chamber engines. This paper reviews the development of pre-chamber engine jet combustion technology and reports recent advancements in laser diagnostics and numerical simulations of pre-chamber engine jet combustion. Research indicates that the ignition modes of pre-chamber engines can be categorized into “flame ignition” and “jet ignition”. In the flame ignition mode, flame quenching does not occur at the pre-chamber nozzle, resulting in good combustion stability and improved engine thermal efficiency. Due to the large flame quenching distance of ammonia fuel, pre-chamber ignition tends to exhibit the jet ignition mode, which can lead to reduced combustion stability and thermal efficiency of the engines.
  • ZHANG Wei, LIU Weijie, ZHOU Yi, SU He, XU Chuanlong
    Journal of Engineering Thermophysics. 2025, 46(6): 2055-2063.
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    Experiments of the free radical chemiluminescence properties of ethanol-fueled spray swirling flame are carried out through an UV camera and an enhanced spectrometer. The effect of free radical chemiluminescence on the spray swirling flame structural characteristics, equivalence ratio and heat release rate are investigated. A model for characterizing the parameters of the spray swirling flame based on the chemiluminescence intensities of the radicals has been established by OH, CH, and C2 at 310 nm, 431 nm, and 516 nm bands. The relationships between the chemiluminescence intensity ratios of various radical combinations (OH/CH, OH/C2, CH/C2), and the equivalence ratios and heat release rates are analyzed and characterized. Results show that the chemiluminescence intensities of OH*, CH*, and C2* radicals exhibit a significant linear correlation with the equivalence ratios and heat release rates under diverse gas flow velocities The intensity ratio of OH/C2 chemiluminescence in radical combinations exist the strongest correlation with flame parameters and is capable of characterizing the equivalence ratio and heat release rate of spray swirling flames.
  • SHAO Jianing, YANG Li, CAO Yunqi, SONG Chen, LIU Fang
    Journal of Engineering Thermophysics. 2025, 46(6): 2064-2072.
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    This study, based on chemical looping technology, investigates the recycling and reduction characteristics of spent LiCoO2 batteries. Hydrogen reduction was employed to reduce spent LiCoO2 in a chemical looping reactor, and the reduction rate, recovery rate, and reaction mechanism of LiCoO2 were analyzed under different temperature and time conditions. The experimental results indicate that temperature is the most important factor influencing the reduction reaction of LiCoO2, with higher temperatures significantly improving the reduction rate. The optimal reduction conditions were found to be 800°C for 120 min. Thermodynamic analysis revealed the specific process of LiCoO2 reduction. Simulation results showed that an appropriate hydrogen concentration can reduce energy consumption during the reduction process. This study provides new insights for the efficient recycling of spent LiCoO2 batteries and lays the experimental foundation for the application of chemical looping hydrogen reduction.
  • ZHANG Haiyang, TU Yaojie, LIU Zirui, LIU Hao
    Journal of Engineering Thermophysics. 2025, 46(6): 2073-2080.
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    In this study, the experimental test of CH4/NH3 swirling premixed flames with axial air-staging was conducted based in a 5 kW model combustor, and the effects of air-staging ratio, staged-air injection height and the number of staged-air injector (N) on the combustion and emission characteristics of CH4/NH3 swirling premixed flame were investigated. The results show that the flame length becomes longer as XNH3 rises, and NO emission firstly increases and then drops, which exhibits an inflection at XNH3 of 50%. NO emission can be effectively reduced by using axial airstaging to make the primary combustion chamber in a rich combustion state under 50% ammonia doping condition, in which the optimal proportion of staged-air is 30%∼50% (equivalence ratio 1.21∼1.70), and it tends to decrease as the proportion of staged-air increases. In addition, under the best air-staging ratio and height conditions (= 80 mm, SAR = 30%), increasing the number of staged-air injector can mitigate the influence of staged-air impinging on the fuel-rich condition of the upstream main combustion zone, and further reduces NO emission.
  • DONG Wenlong, YANG Yuhang, YAO Jinfang, YAN Liang, HONG Run, CHU Huaqiang
    Journal of Engineering Thermophysics. 2025, 46(6): 2081-2089.
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    In this paper, n-heptane, iso-octane and 2,5-dimethylfuran were selected as the objects of study, and the effects of different system densities (0.01∼0.2 g·cm−3) and simulation temperatures (2000∼3500 K) on the generation of products, the generation and breakage of bonds, and the characteristics of the largest molecules during the pyrolysis of gasoline alternative fuel doped with 2,5-dimethylfuran were analyzed on the atomistic level by means of reactive force field molecular dynamics simulations. It was shown that the increase of system density inhibited the generation of C1, C2, and C3 species quantities, slowed down the C-C breakage and inhibited the C-C bond breakage. The lower simulation temperature was not favorable to study the process of soot particle formation, while the higher simulation temperature led to the decomposition of soot particles. The H/C ratios of the largest molecules were 0.314, 0.181, 0.204, 0.184 with the increase of system density. The H/C ratios of the largest molecules were 1.467, 0.902, 0.126, 0.0006 with the increase of simulation temperature. The effect of simulation temperature on the morphology and H/C ratio of soot was more obvious than the system density.