Pub Date : 2026-01-10DOI: 10.1016/j.ijheatmasstransfer.2026.128354
Weipeng Deng , Zhufeng Lin , Huinan Yang , Fengzhong Sun , Huaqiang Chu
Due to the lack of quantitative evaluation methods for plumes, the design of plume control for large wet cooling towers still faces significant challenges. This study uses the outlet surface of cooling towers as the reference and proposes methods for calculating the plume area and the mass flow rate of the plume. Furthermore, a comparative analysis is conducted on the impact of internal (IN_DCZ) and external (EX_DCZ) dry cooling zones on the performance of hybrid cooling towers. Under wind conditions, the performance of IN_DCZ is better than that of EX_DCZ. Compared to EX_DCZ, the outlet water temperature of IN_DCZ reduces by 0.4 °C, the water-saving rate of IN_DCZ increases by 2.2%, and the area and mass flow rate of plume for IN_DCZ are reduced by 45.5 % (663.3 m2) and 50.3 % (6969.7 kg/s), respectively. When the inlet water ratio of the dry cooling zone exceeds 70 %, the plume control capability of IN_DCZ far exceeds that of EX_DCZ, reducing the plume area and mass flow rate by 98.2 % (1054.1 m2) and 98.7 % (9989.2 kg/s) compared to EX_DCZ. This study provides an application exploration for precise plume control in large wet cooling towers.
{"title":"A comparative study on the effects of external and internal dry cooling zones on the performance of hybrid natural draft cooling towers","authors":"Weipeng Deng , Zhufeng Lin , Huinan Yang , Fengzhong Sun , Huaqiang Chu","doi":"10.1016/j.ijheatmasstransfer.2026.128354","DOIUrl":"10.1016/j.ijheatmasstransfer.2026.128354","url":null,"abstract":"<div><div>Due to the lack of quantitative evaluation methods for plumes, the design of plume control for large wet cooling towers still faces significant challenges. This study uses the outlet surface of cooling towers as the reference and proposes methods for calculating the plume area and the mass flow rate of the plume. Furthermore, a comparative analysis is conducted on the impact of internal (IN_DCZ) and external (EX_DCZ) dry cooling zones on the performance of hybrid cooling towers. Under wind conditions, the performance of IN_DCZ is better than that of EX_DCZ. Compared to EX_DCZ, the outlet water temperature of IN_DCZ reduces by 0.4 °C, the water-saving rate of IN_DCZ increases by 2.2%, and the area and mass flow rate of plume for IN_DCZ are reduced by 45.5 % (663.3 m<sup>2</sup>) and 50.3 % (6969.7 kg/s), respectively. When the inlet water ratio of the dry cooling zone exceeds 70 %, the plume control capability of IN_DCZ far exceeds that of EX_DCZ, reducing the plume area and mass flow rate by 98.2 % (1054.1 m<sup>2</sup>) and 98.7 % (9989.2 kg/s) compared to EX_DCZ. This study provides an application exploration for precise plume control in large wet cooling towers.</div></div>","PeriodicalId":336,"journal":{"name":"International Journal of Heat and Mass Transfer","volume":"259 ","pages":"Article 128354"},"PeriodicalIF":5.8,"publicationDate":"2026-01-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145940843","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-10DOI: 10.1016/j.ijheatmasstransfer.2026.128367
Hengyu Hu, Haibo Huang
This study investigates pool boiling enhanced by the well-ordered porous layers using direct numerical simulations based on a coupled lattice Boltzmann method (LBM) and finite difference method (FDM). The simulations capture the dynamic behavior of the vapor–liquid interface during boiling. A physical model is developed to establish a scaling law relating the thickness of porous layers to critical boiling times. Specifically, the inception time , phase change time , and gross boiling time all exhibit a quadratic dependence on . Furthermore, the study shows that under certain conditions, boiling bypasses the nucleate and transition regimes and proceeds directly to film boiling, resulting in the absence of a critical heat flux (CHF). This phenomenon is explained using Rayleigh–Taylor instability theory, from which a criterion is derived to predict the existence of CHF based on the relationship between the critical instability wavelength and the spacing between nucleation sites. The findings offer new insight into boiling mechanisms on porous layers and provide theoretical guidance for designing thermally efficient porous structures.
{"title":"Boiling dynamics on well-ordered porous layers: A lattice Boltzmann simulation study","authors":"Hengyu Hu, Haibo Huang","doi":"10.1016/j.ijheatmasstransfer.2026.128367","DOIUrl":"10.1016/j.ijheatmasstransfer.2026.128367","url":null,"abstract":"<div><div>This study investigates pool boiling enhanced by the well-ordered porous layers using direct numerical simulations based on a coupled lattice Boltzmann method (LBM) and finite difference method (FDM). The simulations capture the dynamic behavior of the vapor–liquid interface during boiling. A physical model is developed to establish a scaling law relating the thickness of porous layers <span><math><mi>H</mi></math></span> to critical boiling times. Specifically, the inception time <span><math><msub><mrow><mi>t</mi></mrow><mrow><mi>i</mi></mrow></msub></math></span>, phase change time <span><math><msub><mrow><mi>t</mi></mrow><mrow><mi>p</mi></mrow></msub></math></span>, and gross boiling time <span><math><msub><mrow><mi>t</mi></mrow><mrow><mi>g</mi></mrow></msub></math></span> all exhibit a quadratic dependence on <span><math><mi>H</mi></math></span>. Furthermore, the study shows that under certain conditions, boiling bypasses the nucleate and transition regimes and proceeds directly to film boiling, resulting in the absence of a critical heat flux (CHF). This phenomenon is explained using Rayleigh–Taylor instability theory, from which a criterion is derived to predict the existence of CHF based on the relationship between the critical instability wavelength and the spacing between nucleation sites. The findings offer new insight into boiling mechanisms on porous layers and provide theoretical guidance for designing thermally efficient porous structures.</div></div>","PeriodicalId":336,"journal":{"name":"International Journal of Heat and Mass Transfer","volume":"259 ","pages":"Article 128367"},"PeriodicalIF":5.8,"publicationDate":"2026-01-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145974941","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-10DOI: 10.1016/j.ijheatmasstransfer.2026.128348
Jun Bi , Huiqi Xia , Sheng Yang , Shengnian Wang , Tingting Wei , Wansheng Pei
Petroleum hydrocarbons (PHs) are compounds containing the elements carbon and hydrogen. The leaks of PHs during the PHs exploration, transportation, refining, and storage process usually contaminate the site, disrupt ecological environment and threaten human health. Therefore, the in-situ thermal desorption technology has been widely used to heat the contaminated sites and remove the pollutants, but the contaminants removal efficiency is controlled by thermal conductivity (TC) in the contaminated site. In this study, a theoretical model was developed by using different preferential invasion pores (i.e., parallel and series pores) and connection patterns between air and liquid mixture of water and PHs (i.e., parallel and series connections). The new model was expressed as a function of porosity, temperature, degree of saturation of PHs, and degree of saturation of water. The theoretical model was evaluated with TC values of diesel-contaminated soils, and good agreements between estimated and measured TC values were obtained. Moreover, all the estimated TC values determined by the new model strictly fall within the Wiener bounds, and most estimated TC results fall within the Hashin-Shtrikman bounds. The theoretical TC model will contribute to the understanding of heat transfer process in the PHs-contaminated site during the remediation process.
{"title":"Modeling on thermal conductivity of petroleum hydrocarbons-contaminated soils considering preferential invasion pores and connection patterns","authors":"Jun Bi , Huiqi Xia , Sheng Yang , Shengnian Wang , Tingting Wei , Wansheng Pei","doi":"10.1016/j.ijheatmasstransfer.2026.128348","DOIUrl":"10.1016/j.ijheatmasstransfer.2026.128348","url":null,"abstract":"<div><div>Petroleum hydrocarbons (PHs) are compounds containing the elements carbon and hydrogen. The leaks of PHs during the PHs exploration, transportation, refining, and storage process usually contaminate the site, disrupt ecological environment and threaten human health. Therefore, the in-situ thermal desorption technology has been widely used to heat the contaminated sites and remove the pollutants, but the contaminants removal efficiency is controlled by thermal conductivity (TC) in the contaminated site. In this study, a theoretical model was developed by using different preferential invasion pores (i.e., parallel and series pores) and connection patterns between air and liquid mixture of water and PHs (i.e., parallel and series connections). The new model was expressed as a function of porosity, temperature, degree of saturation of PHs, and degree of saturation of water. The theoretical model was evaluated with TC values of diesel-contaminated soils, and good agreements between estimated and measured TC values were obtained. Moreover, all the estimated TC values determined by the new model strictly fall within the Wiener bounds, and most estimated TC results fall within the Hashin-Shtrikman bounds. The theoretical TC model will contribute to the understanding of heat transfer process in the PHs-contaminated site during the remediation process.</div></div>","PeriodicalId":336,"journal":{"name":"International Journal of Heat and Mass Transfer","volume":"259 ","pages":"Article 128348"},"PeriodicalIF":5.8,"publicationDate":"2026-01-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145940844","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-09DOI: 10.1016/j.ijheatmasstransfer.2026.128349
M.S. Bidou, J.-G. Bauzin, N. Laraqi
In thermal engineering, the effectiveness of diagnosis, monitoring, and control depends on accurately identifying material properties and locating internal heat sources from limited temperature measurements. This work advances a physics-based hybrid learning framework for 2D thermal conduction, which combines data and physical priors to stabilize inference under limited measurement conditions. We address two complementary objectives: (i) thermophysical identification within an orthotropic medium, estimating the in-plane conductivity tensor and the volumetric heat capacity (with diffusivity reported a posteriori); and (ii) internal volumetric heat-source identification in an isotropic medium, recovering the Gaussian source’s flux amplitude, spatial position, and spatial dispersion, for both fixed and mobile centers. A PINN-based sensitivity analysis guides sensor placement to enhance identifiability with few measurements. Performance is evaluated on noise-free data and under additive Gaussian sensor noise (0.1 K and 0.5 K). Across both settings, the approach delivers high accuracy for material parameters and source localization, and low errors for source flux, supporting robust thermal inverse diagnostics with minimal instrumentation.
{"title":"Physics-informed neural networks for 2D-transient inverse heat conduction problem involving static and moving heat source","authors":"M.S. Bidou, J.-G. Bauzin, N. Laraqi","doi":"10.1016/j.ijheatmasstransfer.2026.128349","DOIUrl":"10.1016/j.ijheatmasstransfer.2026.128349","url":null,"abstract":"<div><div>In thermal engineering, the effectiveness of diagnosis, monitoring, and control depends on accurately identifying material properties and locating internal heat sources from limited temperature measurements. This work advances a physics-based hybrid learning framework for 2D thermal conduction, which combines data and physical priors to stabilize inference under limited measurement conditions. We address two complementary objectives: (i) thermophysical identification within an orthotropic medium, estimating the in-plane conductivity tensor <span><math><mrow><mi>Λ</mi><mo>=</mo><mi>diag</mi><mrow><mo>(</mo><msub><mrow><mi>λ</mi></mrow><mrow><mi>x</mi></mrow></msub><mo>,</mo><msub><mrow><mi>λ</mi></mrow><mrow><mi>y</mi></mrow></msub><mo>)</mo></mrow></mrow></math></span> and the volumetric heat capacity <span><math><mrow><mi>ρ</mi><mi>C</mi></mrow></math></span> (with diffusivity reported a posteriori); and (ii) internal volumetric heat-source identification in an isotropic medium, recovering the Gaussian source’s flux amplitude, spatial position, and spatial dispersion, for both fixed and mobile centers. A PINN-based sensitivity analysis guides sensor placement to enhance identifiability with few measurements. Performance is evaluated on noise-free data and under additive Gaussian sensor noise (0.1 K and 0.5 K). Across both settings, the approach delivers high accuracy for material parameters and source localization, and low errors for source flux, supporting robust thermal inverse diagnostics with minimal instrumentation.</div></div>","PeriodicalId":336,"journal":{"name":"International Journal of Heat and Mass Transfer","volume":"259 ","pages":"Article 128349"},"PeriodicalIF":5.8,"publicationDate":"2026-01-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145940773","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-09DOI: 10.1016/j.ijheatmasstransfer.2026.128342
Songbo Nan , Zeyuan Cheng , Xiang Li , Jingsong Huang , Jianqin Zhu
Flow instability of supercritical hydrocarbon fuels in regenerative cooling channels is a critical issue faced by scramjet engines. Existing research has largely overlooked the influence of pyrolysis on flow instability. This work numerically investigated the instability behaviors and mechanism of supercritical n-decane in a vertical tube, emphasizing the effects of inlet flow rate and pyrolysis. The instability boundaries under different pressures were also established. The results show that the pressure-drop curve exhibits a quintic behavior with two negative-slope regions, caused by sharp density variations in the pseudo-critical and pyrolysis temperature regions. At an inlet flow rate of 0.8 g/s, the fuel exhibits periodic density-wave oscillation at a frequency of 0.23 Hz. Density variation is the fundamental cause of instability, while buoyancy and pyrolysis act as key factors of the periodic oscillations. As the inlet mass flow rate increases from 0.41 g/s to 1.2 g/s, the instability evolves from flow drift coupled with density-wave oscillations to pure density-wave oscillations, and eventually to stable flow, with stability initially decreasing and then increasing. Pyrolysis transforms the pressure drop curve from a cubic to a quintic behavior, which in turn alters the mode of flow instability in the negative-slope regions of the pressure drop curve and shortens the oscillation period from 4.5 s to 4.3 s. The predicted flow instability boundaries under 2.5–4.0 MPa show a maximum deviation of 9.3 % from experimental data. These findings provide theoretical guidance for the safe and reliable design of regenerative cooling channels.
{"title":"Numerical investigation on flow instability of supercritical n-decane in a vertical upward tube with pyrolysis effects","authors":"Songbo Nan , Zeyuan Cheng , Xiang Li , Jingsong Huang , Jianqin Zhu","doi":"10.1016/j.ijheatmasstransfer.2026.128342","DOIUrl":"10.1016/j.ijheatmasstransfer.2026.128342","url":null,"abstract":"<div><div>Flow instability of supercritical hydrocarbon fuels in regenerative cooling channels is a critical issue faced by scramjet engines. Existing research has largely overlooked the influence of pyrolysis on flow instability. This work numerically investigated the instability behaviors and mechanism of supercritical n-decane in a vertical tube, emphasizing the effects of inlet flow rate and pyrolysis. The instability boundaries under different pressures were also established. The results show that the pressure-drop curve exhibits a quintic behavior with two negative-slope regions, caused by sharp density variations in the pseudo-critical and pyrolysis temperature regions. At an inlet flow rate of 0.8 <em>g</em>/s, the fuel exhibits periodic density-wave oscillation at a frequency of 0.23 Hz. Density variation is the fundamental cause of instability, while buoyancy and pyrolysis act as key factors of the periodic oscillations. As the inlet mass flow rate increases from 0.41 g/s to 1.2 g/s, the instability evolves from flow drift coupled with density-wave oscillations to pure density-wave oscillations, and eventually to stable flow, with stability initially decreasing and then increasing. Pyrolysis transforms the pressure drop curve from a cubic to a quintic behavior, which in turn alters the mode of flow instability in the negative-slope regions of the pressure drop curve and shortens the oscillation period from 4.5 s to 4.3 s. The predicted flow instability boundaries under 2.5–4.0 MPa show a maximum deviation of 9.3 % from experimental data. These findings provide theoretical guidance for the safe and reliable design of regenerative cooling channels.</div></div>","PeriodicalId":336,"journal":{"name":"International Journal of Heat and Mass Transfer","volume":"259 ","pages":"Article 128342"},"PeriodicalIF":5.8,"publicationDate":"2026-01-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145940842","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-09DOI: 10.1016/j.ijheatmasstransfer.2026.128327
Anastasia Islamova, Svetlana Kropotova, Andrey Klimenko, Stanislav Shulyaev
Understanding the specific aspects of binary interactions of coal-water slurry droplets is crucial for optimizing the transportation of this fuel and spraying in industrial applications. This research presents experimental findings for binary collisions of coal-water slurry droplets, when varying the ambient air temperature (from 25 to 200°C), concentration of solid coal particles (from 0.5 to 5%), and geometry of collisions. Using high-speed video recording, the following interaction regimes of droplets were recorded: bounce, coalescence, reflexive separation, and stretching separation. Droplet interaction regime maps were plotted to identify the boundaries between regimes depending on the solid phase concentration in the slurry and ambient temperature. The probability of occurrence of a regime was not significantly dependent on the coal concentration in the slurry in the considered range, but was very sensitive to the impact angle and velocity. It was hypothesized that the geometry of the lamella forming in off-center interaction determines the number of secondary fragments (child droplets) after its breakup. The inertial forces of droplets and thermophysical properties of the gas medium rather than the solid phase concentration in the liquid have a decisive effect on the droplet interaction outcome.
{"title":"Binary collisions between slurry droplets at different ambient temperatures","authors":"Anastasia Islamova, Svetlana Kropotova, Andrey Klimenko, Stanislav Shulyaev","doi":"10.1016/j.ijheatmasstransfer.2026.128327","DOIUrl":"10.1016/j.ijheatmasstransfer.2026.128327","url":null,"abstract":"<div><div>Understanding the specific aspects of binary interactions of coal-water slurry droplets is crucial for optimizing the transportation of this fuel and spraying in industrial applications. This research presents experimental findings for binary collisions of coal-water slurry droplets, when varying the ambient air temperature (from 25 to 200°C), concentration of solid coal particles (from 0.5 to 5%), and geometry of collisions. Using high-speed video recording, the following interaction regimes of droplets were recorded: bounce, coalescence, reflexive separation, and stretching separation. Droplet interaction regime maps were plotted to identify the boundaries between regimes depending on the solid phase concentration in the slurry and ambient temperature. The probability of occurrence of a regime was not significantly dependent on the coal concentration in the slurry in the considered range, but was very sensitive to the impact angle and velocity. It was hypothesized that the geometry of the lamella forming in off-center interaction determines the number of secondary fragments (child droplets) after its breakup. The inertial forces of droplets and thermophysical properties of the gas medium rather than the solid phase concentration in the liquid have a decisive effect on the droplet interaction outcome.</div></div>","PeriodicalId":336,"journal":{"name":"International Journal of Heat and Mass Transfer","volume":"259 ","pages":"Article 128327"},"PeriodicalIF":5.8,"publicationDate":"2026-01-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145940848","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-09DOI: 10.1016/j.ijheatmasstransfer.2026.128333
Youping Li , Yue Kuang , Shuang Ma , Xuxia Sun , Zheng Jiang , Han Jiang , Qirong Yang
ZnCl2-NaCl-KCl eutectic salt is considered a promising latent heat phase change material (PCM) that exhibits significant potential for development and application in thermal energy storage systems. However, given its strong corrosiveness and high propensity for leakage, it is typically encapsulated within porous framework materials to enhance stability and safety during practical applications. Different pore structures and framework types can significantly impact the thermal behavior and phase change characteristics in composite phase change materials (CPCMs) through interfacial interactions and structural configurations. However, research on the structure and interface effects of chloride-based salts remains scarce, and their underlying mechanisms on these effects within CPCMs are not yet fully elucidated. In this study, interfacial effects on thermal behavior and phase change characteristics of ZnCl2-NaCl-KCl eutectic salts were systematically explored by integrating the molecular dynamics (MD) simulation with experimental analyses. Additionally, the structural effects on thermal behavior and phase-change characteristics of ternary chloride salts mixture were explored using MD simulations. The results indicate that the interfacial and structural effects are primarily influenced by the ionic structure and interfacial binding energy. Among various interfacial composite models, the salt on Al2O3 composite material exhibits the highest thermal conductivity, specific heat capacity, and latent heat of phase change. Among various structural composite models, the ink-bottle shaped CPCMs achieve higher thermal conductivity than that of the pure chloride salts and other CPCM structures. The specific heat capacity follows a similar trend to the thermal conductivity. The ink-bottle shaped pore structure of the Al2O3 framework is identified as the most appropriate configuration for accommodating the ZnCl2-NaCl-KCl eutectic salt.
{"title":"Interface and structure effects on thermal properties and phase-change characteristics of chloride-based CPCMs toward sustainable thermal storage","authors":"Youping Li , Yue Kuang , Shuang Ma , Xuxia Sun , Zheng Jiang , Han Jiang , Qirong Yang","doi":"10.1016/j.ijheatmasstransfer.2026.128333","DOIUrl":"10.1016/j.ijheatmasstransfer.2026.128333","url":null,"abstract":"<div><div>ZnCl<sub>2</sub>-NaCl-KCl eutectic salt is considered a promising latent heat phase change material (PCM) that exhibits significant potential for development and application in thermal energy storage systems. However, given its strong corrosiveness and high propensity for leakage, it is typically encapsulated within porous framework materials to enhance stability and safety during practical applications. Different pore structures and framework types can significantly impact the thermal behavior and phase change characteristics in composite phase change materials (CPCMs) through interfacial interactions and structural configurations. However, research on the structure and interface effects of chloride-based salts remains scarce, and their underlying mechanisms on these effects within CPCMs are not yet fully elucidated. In this study, interfacial effects on thermal behavior and phase change characteristics of ZnCl<sub>2</sub>-NaCl-KCl eutectic salts were systematically explored by integrating the molecular dynamics (MD) simulation with experimental analyses. Additionally, the structural effects on thermal behavior and phase-change characteristics of ternary chloride salts mixture were explored using MD simulations. The results indicate that the interfacial and structural effects are primarily influenced by the ionic structure and interfacial binding energy. Among various interfacial composite models, the salt on Al<sub>2</sub>O<sub>3</sub> composite material exhibits the highest thermal conductivity, specific heat capacity, and latent heat of phase change. Among various structural composite models, the ink-bottle shaped CPCMs achieve higher thermal conductivity than that of the pure chloride salts and other CPCM structures. The specific heat capacity follows a similar trend to the thermal conductivity. The ink-bottle shaped pore structure of the Al<sub>2</sub>O<sub>3</sub> framework is identified as the most appropriate configuration for accommodating the ZnCl<sub>2</sub>-NaCl-KCl eutectic salt.</div></div>","PeriodicalId":336,"journal":{"name":"International Journal of Heat and Mass Transfer","volume":"259 ","pages":"Article 128333"},"PeriodicalIF":5.8,"publicationDate":"2026-01-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145940850","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-09DOI: 10.1016/j.ijheatmasstransfer.2026.128345
Christian Psenica , Lean Fang , Seth Zoppelt , Mark Leader , Ping He
Conjugate heat transfer (CHT) analysis and optimization is a powerful method for improving thermal management as it simultaneously resolves the temperature distribution in both fluid and solid domains. This paper presents a modular, discrete adjoint-based CHT optimization capability integrated within the OpenMDAO/MPhys framework. A unique feature of the proposed framework is its flexibility to extend to multidisciplinary optimization, including aero-structural-thermal applications. The fluid domain is modeled using a finite-volume Computational Fluid Dynamics (CFD) solver, and the solid domain with a conduction heat transfer solver. A mixed Neumann-Dirichlet boundary condition is developed to enable full submersion of the solid geometry within the fluid domain, while ensuring consistent temperature and heat flux coupling at the CHT interface. Gradient-based optimization is performed; the gradients are efficiently computed using the discrete adjoint solvers implemented in DAFoam. To demonstrate the method, this paper considers two cases related to electric aircraft thermal management: a U-bend heat exchanger and an actively cooled battery pack. The U-bend case aims to minimize pressure loss while maximizing heat flux by changing the pipe geometry. The optimized design reduces pressure loss by 52.7% and increases total heat flux by 2.3%. In the battery pack case, a 3-by-3 cell configuration is cooled by ambient airflow, with constant heat generation prescribed in the cells. The battery casing shape serves as the design variable, and the objective function is a weighted sum of pressure loss and pack weight, subject to a maximum temperature constraint. The optimized design achieves a 44.6% reduction in pressure loss and a 1.5% reduction in weight, while satisfying the thermal constraint. To ensure the reliability of the optimized designs, this study validates coarse-mesh, steady-state predictions against fine-mesh unsteady simulations, demonstrating consistency within acceptable errors. This work demonstrates the potential of the developed framework to enable rapid, high-fidelity design of thermal management systems for electric aircraft.
{"title":"A modular conjugate heat transfer optimization framework for thermal management of electric aircraft","authors":"Christian Psenica , Lean Fang , Seth Zoppelt , Mark Leader , Ping He","doi":"10.1016/j.ijheatmasstransfer.2026.128345","DOIUrl":"10.1016/j.ijheatmasstransfer.2026.128345","url":null,"abstract":"<div><div>Conjugate heat transfer (CHT) analysis and optimization is a powerful method for improving thermal management as it simultaneously resolves the temperature distribution in both fluid and solid domains. This paper presents a modular, discrete adjoint-based CHT optimization capability integrated within the OpenMDAO/MPhys framework. A unique feature of the proposed framework is its flexibility to extend to multidisciplinary optimization, including aero-structural-thermal applications. The fluid domain is modeled using a finite-volume Computational Fluid Dynamics (CFD) solver, and the solid domain with a conduction heat transfer solver. A mixed Neumann-Dirichlet boundary condition is developed to enable full submersion of the solid geometry within the fluid domain, while ensuring consistent temperature and heat flux coupling at the CHT interface. Gradient-based optimization is performed; the gradients are efficiently computed using the discrete adjoint solvers implemented in DAFoam. To demonstrate the method, this paper considers two cases related to electric aircraft thermal management: a U-bend heat exchanger and an actively cooled battery pack. The U-bend case aims to minimize pressure loss while maximizing heat flux by changing the pipe geometry. The optimized design reduces pressure loss by 52.7% and increases total heat flux by 2.3%. In the battery pack case, a 3-by-3 cell configuration is cooled by ambient airflow, with constant heat generation prescribed in the cells. The battery casing shape serves as the design variable, and the objective function is a weighted sum of pressure loss and pack weight, subject to a maximum temperature constraint. The optimized design achieves a 44.6% reduction in pressure loss and a 1.5% reduction in weight, while satisfying the thermal constraint. To ensure the reliability of the optimized designs, this study validates coarse-mesh, steady-state predictions against fine-mesh unsteady simulations, demonstrating consistency within acceptable errors. This work demonstrates the potential of the developed framework to enable rapid, high-fidelity design of thermal management systems for electric aircraft.</div></div>","PeriodicalId":336,"journal":{"name":"International Journal of Heat and Mass Transfer","volume":"259 ","pages":"Article 128345"},"PeriodicalIF":5.8,"publicationDate":"2026-01-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145940846","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-09DOI: 10.1016/j.ijheatmasstransfer.2026.128329
Amin Song , Zhi Yi , Wenchao Ji , Xin Zhang , Xin Guo , Shuangcheng Sun , Linyang Wei , Guojun Li
The Zone Method (ZM) known for high accuracy plays a critical role in solving the radiative heat transfer (RHT) in the high-temperature industrial fields. However, this method still faces challenges in handling high-dimensional integrals for direct exchange areas (DEAs), while having limitations in solving total exchange areas (TEAs) for anisotropic scattering media based on traditional integration method (IM). To address this issue, a reduced integration coupled with Monte Carlo ratios method (RIMCR) is developed to solve the RHT in cylindrical participating media in this study. In this method, the reduced integration scheme (RIS) is used to reduce dimension of DEAs integral formulation for cylindrical geometry, and the proportionality coefficients derived from the DEAs are employed to solve the TEAs iteratively, in which the iterative ratios of gas micro-elements within anisotropic scattering media are determined via tracking the path of energy beams by Monte Carlo method (MCM). A series of cases in two-dimensional axisymmetric configurations with uniform medium properties demonstrates the RIMCR not only has high accuracy, but also significantly improves computational efficiency in solving the two-dimensional cylindrical RHT. Notably, RIMCR achieves a computational efficiency improvement exceeding 75 % compared to the IM. Therefore, RIMCR is a reliable method for effective and accurate solution of the RHT in participating media.
{"title":"Reduced integration coupled with Monte Carlo ratios method for radiative heat transfer in cylindrical participating media","authors":"Amin Song , Zhi Yi , Wenchao Ji , Xin Zhang , Xin Guo , Shuangcheng Sun , Linyang Wei , Guojun Li","doi":"10.1016/j.ijheatmasstransfer.2026.128329","DOIUrl":"10.1016/j.ijheatmasstransfer.2026.128329","url":null,"abstract":"<div><div>The Zone Method (ZM) known for high accuracy plays a critical role in solving the radiative heat transfer (RHT) in the high-temperature industrial fields. However, this method still faces challenges in handling high-dimensional integrals for direct exchange areas (DEAs), while having limitations in solving total exchange areas (TEAs) for anisotropic scattering media based on traditional integration method (IM). To address this issue, a reduced integration coupled with Monte Carlo ratios method (RIMCR) is developed to solve the RHT in cylindrical participating media in this study. In this method, the reduced integration scheme (RIS) is used to reduce dimension of DEAs integral formulation for cylindrical geometry, and the proportionality coefficients derived from the DEAs are employed to solve the TEAs iteratively, in which the iterative ratios of gas micro-elements within anisotropic scattering media are determined via tracking the path of energy beams by Monte Carlo method (MCM). A series of cases in two-dimensional axisymmetric configurations with uniform medium properties demonstrates the RIMCR not only has high accuracy, but also significantly improves computational efficiency in solving the two-dimensional cylindrical RHT. Notably, RIMCR achieves a computational efficiency improvement exceeding 75 % compared to the IM. Therefore, RIMCR is a reliable method for effective and accurate solution of the RHT in participating media.</div></div>","PeriodicalId":336,"journal":{"name":"International Journal of Heat and Mass Transfer","volume":"259 ","pages":"Article 128329"},"PeriodicalIF":5.8,"publicationDate":"2026-01-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145940847","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-09DOI: 10.1016/j.ijheatmasstransfer.2026.128338
Benchi Liu , Jing Wang , Akira Nakayama , Zihao Zheng
In this study, a systematic three-dimensional numerical study and experimental study were conducted for forced convection in a series of isothermally heated bifurcation windward bend structures. The numerical experiment was based on an in-house computer program with the Semi-implicit Method for Pressure-linked Equation (SIMPLE algorithm) and Finite Volume Method (FVM) to solve the continuity equation, Navier-Stokes equation, and energy equation of the fluid phase and solid phase. The heat transfer performance of various bifurcation windward bend structures was investigated by using experimental method and compared the results with numerical simulation. It shows numerical results and experimental results match well, the error between these two is less than 10%. A set of Nusselt numbers under equal Reynolds numbers and pumping powers were calculated to investigate the impact of structural optimization on enhancing heat transfer performance. The results point out the bifurcated windward bend structure with that provides the best heat transfer performance. This structure exhibits excellent heat transfer performance due to its high thermal dispersion and low flow resistance. This numerical study provides evidence that the bifurcation design strategy is an effective method to design and optimize heat exchanger systems.
{"title":"Experimental and numerical studies on the heat transfer enhancement of the bifurcation windward bend structure","authors":"Benchi Liu , Jing Wang , Akira Nakayama , Zihao Zheng","doi":"10.1016/j.ijheatmasstransfer.2026.128338","DOIUrl":"10.1016/j.ijheatmasstransfer.2026.128338","url":null,"abstract":"<div><div>In this study, a systematic three-dimensional numerical study and experimental study were conducted for forced convection in a series of isothermally heated bifurcation windward bend structures. The numerical experiment was based on an in-house computer program with the Semi-implicit Method for Pressure-linked Equation (SIMPLE algorithm) and Finite Volume Method (FVM) to solve the continuity equation, Navier-Stokes equation, and energy equation of the fluid phase and solid phase. The heat transfer performance of various bifurcation windward bend structures was investigated by using experimental method and compared the results with numerical simulation. It shows numerical results and experimental results match well, the error between these two is less than 10%. A set of Nusselt numbers under equal Reynolds numbers and pumping powers were calculated to investigate the impact of structural optimization on enhancing heat transfer performance. The results point out the bifurcated windward bend structure with <span><math><mrow><mi>θ</mi><mo>=</mo><msup><mrow><mn>30</mn></mrow><mo>∘</mo></msup></mrow></math></span> that provides the best heat transfer performance. This structure exhibits excellent heat transfer performance due to its high thermal dispersion and low flow resistance. This numerical study provides evidence that the bifurcation design strategy is an effective method to design and optimize heat exchanger systems.</div></div>","PeriodicalId":336,"journal":{"name":"International Journal of Heat and Mass Transfer","volume":"259 ","pages":"Article 128338"},"PeriodicalIF":5.8,"publicationDate":"2026-01-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145940841","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}