Pub Date : 2026-09-01DOI: 10.1016/j.ijheatfluidflow.2026.110648
Yangyang Jin, Chong Zhai, Jiahao Gong, Jiayi Wang, Guoxing Wang, Yao Nai, Dong Li
The transient flow physics and conjugate heat transfer mechanisms of highly compressible superheated steam in long-distance pipelines remain incompletely understood. Most existing dynamic models simplify the non-ideal thermophysical properties of steam and represent pipeline thermal inertia using lumped parameters, thereby failing to distinguish the respective roles of the inner metal wall and outer insulation layers. To address this limitation, the present study develops a transient thermo-hydraulic coupled model that resolves multi-layer fluid–solid heat storage. The model accounts for the nonlinear thermophysical variations of superheated steam and employs a high-order TVD-MUSCL scheme to minimize numerical diffusion at the thermal front. Validation against field measurements yields maximum relative errors of − 0.98% for temperature and − 0.38% for pressure. Results indicate that increasing inlet pressure enhances the volumetric heat capacity of steam, which accelerates thermal wave propagation and reduces dynamic pressure drops. Compared with a step change, a heating ramp of 1 °C/min extends the thermal response time by approximately 50% while decreasing the thermal power ramp-up rate by a factor of six. Structural analysis further shows that the metal wall primarily governs transient heat storage and thermal response delay, whereas the insulation layer mainly controls steady-state heat dissipation. These findings clarify the distinct contributions of multi-layer solid boundaries to transient conjugate heat transfer and provide a theoretical basis for the structural optimization and dynamic regulation of large-scale steam transport systems.
{"title":"Unsteady thermo-hydraulic modeling and conjugate heat transfer mechanisms in long-distance superheated steam pipelines","authors":"Yangyang Jin, Chong Zhai, Jiahao Gong, Jiayi Wang, Guoxing Wang, Yao Nai, Dong Li","doi":"10.1016/j.ijheatfluidflow.2026.110648","DOIUrl":"10.1016/j.ijheatfluidflow.2026.110648","url":null,"abstract":"<div><div>The transient flow physics and conjugate heat transfer mechanisms of highly compressible superheated steam in long-distance pipelines remain incompletely understood. Most existing dynamic models simplify the non-ideal thermophysical properties of steam and represent pipeline thermal inertia using lumped parameters, thereby failing to distinguish the respective roles of the inner metal wall and outer insulation layers. To address this limitation, the present study develops a transient thermo-hydraulic coupled model that resolves multi-layer fluid–solid heat storage. The model accounts for the nonlinear thermophysical variations of superheated steam and employs a high-order TVD-MUSCL scheme to minimize numerical diffusion at the thermal front. Validation against field measurements yields maximum relative errors of − 0.98% for temperature and − 0.38% for pressure. Results indicate that increasing inlet pressure enhances the volumetric heat capacity of steam, which accelerates thermal wave propagation and reduces dynamic pressure drops. Compared with a step change, a heating ramp of 1 °C/min extends the thermal response time by approximately 50% while decreasing the thermal power ramp-up rate by a factor of six. Structural analysis further shows that the metal wall primarily governs transient heat storage and thermal response delay, whereas the insulation layer mainly controls steady-state heat dissipation. These findings clarify the distinct contributions of multi-layer solid boundaries to transient conjugate heat transfer and provide a theoretical basis for the structural optimization and dynamic regulation of large-scale steam transport systems.</div></div>","PeriodicalId":335,"journal":{"name":"International Journal of Heat and Fluid Flow","volume":"121 ","pages":"Article 110648"},"PeriodicalIF":4.8,"publicationDate":"2026-09-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148934408","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-09-01Epub Date: 2026-09-02DOI: 10.1016/j.ijheatfluidflow.2026.110629
Wojciech Sadowski, Hakan Demir, Francesca di Mare
The work considers a fully turbulent flow with heat transfer in a channel half-filled with an array of cubes based on the work of Breugem and Boersma (2005) and Chandesris et al. (2013), at and three different Prandtl numbers, . The temperature is modelled as a passive scalar and two different boundary condition configurations are simulated. The influence of the Prandtl number on the mean temperature, its variance and the terms of the temperature budget is highlighted, including the analysis of the distribution and relative importance of the turbulent heat transfer, molecular diffusion, tortuosity and Brinkman terms near the porous-fluid interface. The latter two have been found to be insignificant for the highest . A set of terms, typically neglected during the upscaling procedure (related to the Taylor expansion of the filtered variables), is analysed for the first time for turbulent heat transfer at the porous-fluid interface, and these terms are found to be significant at low . The upscaled fields are evaluated with three different cellular-average kernels: linear (i.e., tent), quadratic and cubic, and the influence of the chosen filter is additionally studied.
{"title":"Influence of Prandtl number on heat transfer over a permeable wall","authors":"Wojciech Sadowski, Hakan Demir, Francesca di Mare","doi":"10.1016/j.ijheatfluidflow.2026.110629","DOIUrl":"10.1016/j.ijheatfluidflow.2026.110629","url":null,"abstract":"<div><div>The work considers a fully turbulent flow with heat transfer in a channel half-filled with an array of cubes based on the work of Breugem and Boersma (2005) and Chandesris et al. (2013), at <span><math><mrow><msub><mrow><mi>Re</mi></mrow><mrow><mi>bulk</mi></mrow></msub><mo>=</mo><mn>5485</mn></mrow></math></span> and three different Prandtl numbers, <span><math><mrow><mi>Pr</mi><mo>=</mo><mn>0.71,0.1,0.05</mn></mrow></math></span>. The temperature is modelled as a passive scalar and two different boundary condition configurations are simulated. The influence of the Prandtl number on the mean temperature, its variance and the terms of the temperature budget is highlighted, including the analysis of the distribution and relative importance of the turbulent heat transfer, molecular diffusion, tortuosity and Brinkman terms near the porous-fluid interface. The latter two have been found to be insignificant for the highest <span><math><mi>Pr</mi></math></span>. A set of terms, typically neglected during the upscaling procedure (related to the Taylor expansion of the filtered variables), is analysed for the first time for turbulent heat transfer at the porous-fluid interface, and these terms are found to be significant at low <span><math><mi>Pr</mi></math></span>. The upscaled fields are evaluated with three different cellular-average kernels: linear (i.e., tent), quadratic and cubic, and the influence of the chosen filter is additionally studied.</div></div>","PeriodicalId":335,"journal":{"name":"International Journal of Heat and Fluid Flow","volume":"121 ","pages":"Article 110629"},"PeriodicalIF":4.8,"publicationDate":"2026-09-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148934367","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-09-01Epub Date: 2026-08-24DOI: 10.1016/j.ijheatfluidflow.2026.110639
Kazem Mashayekh, Kevin Pope, Yuri S. Muzychka
In this paper, the thermo-hydraulic performance of spherical dimpled tubes is numerically investigated, and their geometric parameters are optimized using a statistical design approach. A validated three-dimensional numerical model is developed for (i) smooth and (ii) dimpled tubes to accurately predict turbulent flow and heat transfer characteristics. The heat transfer enhancement mechanisms are first examined through detailed flow and thermal field analyses using the one-factor-at-a-time (OFAT) method over a Reynolds number (Re) range of 5000–15,000. The investigated geometric parameters include dimple diameter, dimple pitch, and number of dimple stars. Subsequently, Response Surface Methodology (RSM) combined with the Box-Behnken Design (BBD) is employed at Re = 5000 to evaluate the effects of these parameters on the Nusselt number (Nu), friction factor (f), and performance evaluation criterion (PEC). Quadratic and modified quadratic regression models are developed and evaluated using analysis of variance (ANOVA). The ANOVA results reveal that dimple diameter has a significant positive effect on Nu but also increases f. Increasing the dimple pitch reduces Nu, while the interaction between dimple diameter and number of dimple stars significantly affects f and PEC. Multi-objective optimization based on the desirability function identifies the optimal geometric configuration, with prediction errors of less than 3% for all response variables. A comparison with the OFAT approach demonstrates that the BBD-based optimization provides a more efficient and systematic framework for identifying the optimal dimple geometry while substantially reducing the computational effort required for design optimization.
{"title":"Multi-objective optimization of spherically dimpled tubes for enhanced thermal-hydraulic performance","authors":"Kazem Mashayekh, Kevin Pope, Yuri S. Muzychka","doi":"10.1016/j.ijheatfluidflow.2026.110639","DOIUrl":"10.1016/j.ijheatfluidflow.2026.110639","url":null,"abstract":"<div><div>In this paper, the thermo-hydraulic performance of spherical dimpled tubes is numerically investigated, and their geometric parameters are optimized using a statistical design approach. A validated three-dimensional numerical model is developed for (i) smooth and (ii) dimpled tubes to accurately predict turbulent flow and heat transfer characteristics. The heat transfer enhancement mechanisms are first examined through detailed flow and thermal field analyses using the one-factor-at-a-time (OFAT) method over a Reynolds number (<em>Re</em>) range of 5000–15,000. The investigated geometric parameters include dimple diameter, dimple pitch, and number of dimple stars. Subsequently, Response Surface Methodology (RSM) combined with the Box-Behnken Design (BBD) is employed at <em>Re</em> = 5000 to evaluate the effects of these parameters on the Nusselt number (<em>Nu</em>), friction factor (<em>f</em>), and performance evaluation criterion (<em>PEC</em>). Quadratic and modified quadratic regression models are developed and evaluated using analysis of variance (ANOVA). The ANOVA results reveal that dimple diameter has a significant positive effect on <em>Nu</em> but also increases <em>f</em>. Increasing the dimple pitch reduces <em>Nu</em>, while the interaction between dimple diameter and number of dimple stars significantly affects <em>f</em> and <em>PEC</em>. Multi-objective optimization based on the desirability function identifies the optimal geometric configuration, with prediction errors of less than 3% for all response variables. A comparison with the OFAT approach demonstrates that the BBD-based optimization provides a more efficient and systematic framework for identifying the optimal dimple geometry while substantially reducing the computational effort required for design optimization.</div></div>","PeriodicalId":335,"journal":{"name":"International Journal of Heat and Fluid Flow","volume":"121 ","pages":"Article 110639"},"PeriodicalIF":4.8,"publicationDate":"2026-09-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148853931","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
This numerical study investigates the transition process in a Falkner–Skan–Cooke boundary layer on a homogeneously-roughened swept flat plate, with varying roughness wavelengths and heights. Using direct numerical simulations, we determine the onset and development of stationary crossflow vortices (CFVs) and their eventual breakdown via secondary instabilities (SIs). The results highlight the wavelength-selective nature of roughness-induced instability: roughness matching the most unstable mode for the stationary CFV triggers the earliest CFV development and collapse. In contrast, both longer and shorter wavelengths alter the growing mode, ultimately delaying the collapse of CFVs and the transition to turbulence. Through spectral and statistical analyses, we identify clear milestones—occurrences of CFVs and SIs, and collapse of CFVs—which correspond well with flow visualizations. While canonical SI types (Type-1, 2, 3) are observed on short-wavelength roughness, mid- and long-wavelength cases reveal atypical unsteady behavior lacking clear spectral peaks. The observed delay in transition due to wavelength mismatches would have practical relevance to turbulent transition control.
{"title":"DNS study on the roughness wavelength selectivity of crossflow vortex transition due to distributed surface roughness","authors":"Kotaro Minamio, Kosuke Nakagawa, Ryo Araki, Takahiro Ishida, Takahiro Tsukahara","doi":"10.1016/j.ijheatfluidflow.2026.110632","DOIUrl":"10.1016/j.ijheatfluidflow.2026.110632","url":null,"abstract":"<div><div>This numerical study investigates the transition process in a Falkner–Skan–Cooke boundary layer on a homogeneously-roughened swept flat plate, with varying roughness wavelengths and heights. Using direct numerical simulations, we determine the onset and development of stationary crossflow vortices (CFVs) and their eventual breakdown via secondary instabilities (SIs). The results highlight the wavelength-selective nature of roughness-induced instability: roughness matching the most unstable mode for the stationary CFV triggers the earliest CFV development and collapse. In contrast, both longer and shorter wavelengths alter the growing mode, ultimately delaying the collapse of CFVs and the transition to turbulence. Through spectral and statistical analyses, we identify clear milestones—occurrences of CFVs and SIs, and collapse of CFVs—which correspond well with flow visualizations. While canonical SI types (Type-1, 2, 3) are observed on short-wavelength roughness, mid- and long-wavelength cases reveal atypical unsteady behavior lacking clear spectral peaks. The observed delay in transition due to wavelength mismatches would have practical relevance to turbulent transition control.</div></div>","PeriodicalId":335,"journal":{"name":"International Journal of Heat and Fluid Flow","volume":"121 ","pages":"Article 110632"},"PeriodicalIF":4.8,"publicationDate":"2026-09-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148934366","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-09-01DOI: 10.1016/j.ijheatfluidflow.2026.110644
M.M. Molinari, D. Cavalieri, L. Lucchese, J. Liberatori, M. Blandino, M. Valorani, P.P. Ciottoli
Drop-in sustainable aviation fuels coupled with high overall pressure ratio gas turbine engines constitute a promising pathway toward near-term carbon neutrality in the aviation sector. Under typical high-power operating conditions, however, the pressure in the combustor exceeds the fuel’s critical pressure, so injection occurs in a transcritical regime. In this regime, the thermophysical and transport properties of the working fluid exhibit pronounced non-linear variations, which significantly influence the mixing processes. While the effects of pressure and temperature on transcritical injection are well established for cryogenic rocket propellants, substantial uncertainty persists regarding how surrogate fuel composition affects numerical predictions for sustainable aviation fuels. In this work, we first validate a numerical framework for representing transcritical, multicomponent aviation fuels by coupling OpenFOAM with Cantera. We then investigate the effect of the number of chemical species in SAF surrogates on turbulent transcritical mixing under aero-engine relevant conditions through the diffuse-interface large-eddy simulation approach. The study focuses on the alcohol-to-jet fuel POSF-11498, classified as a category C fuel within the National Jet Fuels Combustion Program. The results demonstrate that a two-component surrogate, designed to capture the dominant and highly branched iso-paraffins of POSF-11498, reproduces the behavior of a more detailed seven-component surrogate with no appreciable loss of accuracy, at less than half its computational cost. In contrast, a single-component representation departs appreciably from the multicomponent surrogates in the prediction of density, heat capacity and thermal conductivity. In all cases the mixing field itself is nearly surrogate-independent, so the composition sensitivity stems from the real-fluid thermophysical properties rather than from the mixing dynamics. The proposed numerical framework provides a versatile foundation for high-fidelity LES investigations of transcritical injection of sustainable aviation fuels.
{"title":"Composition effects in transcritical turbulent mixing of sustainable aviation fuel surrogates","authors":"M.M. Molinari, D. Cavalieri, L. Lucchese, J. Liberatori, M. Blandino, M. Valorani, P.P. Ciottoli","doi":"10.1016/j.ijheatfluidflow.2026.110644","DOIUrl":"10.1016/j.ijheatfluidflow.2026.110644","url":null,"abstract":"<div><div>Drop-in sustainable aviation fuels coupled with high overall pressure ratio gas turbine engines constitute a promising pathway toward near-term carbon neutrality in the aviation sector. Under typical high-power operating conditions, however, the pressure in the combustor exceeds the fuel’s critical pressure, so injection occurs in a transcritical regime. In this regime, the thermophysical and transport properties of the working fluid exhibit pronounced non-linear variations, which significantly influence the mixing processes. While the effects of pressure and temperature on transcritical injection are well established for cryogenic rocket propellants, substantial uncertainty persists regarding how surrogate fuel composition affects numerical predictions for sustainable aviation fuels. In this work, we first validate a numerical framework for representing transcritical, multicomponent aviation fuels by coupling OpenFOAM with Cantera. We then investigate the effect of the number of chemical species in SAF surrogates on turbulent transcritical mixing under aero-engine relevant conditions through the diffuse-interface large-eddy simulation approach. The study focuses on the alcohol-to-jet fuel POSF-11498, classified as a category C fuel within the National Jet Fuels Combustion Program. The results demonstrate that a two-component surrogate, designed to capture the dominant <span><math><msub><mrow><mi>C</mi></mrow><mrow><mn>12</mn></mrow></msub></math></span> and <span><math><msub><mrow><mi>C</mi></mrow><mrow><mn>16</mn></mrow></msub></math></span> highly branched iso-paraffins of POSF-11498, reproduces the behavior of a more detailed seven-component surrogate with no appreciable loss of accuracy, at less than half its computational cost. In contrast, a single-component representation departs appreciably from the multicomponent surrogates in the prediction of density, heat capacity and thermal conductivity. In all cases the mixing field itself is nearly surrogate-independent, so the composition sensitivity stems from the real-fluid thermophysical properties rather than from the mixing dynamics. The proposed numerical framework provides a versatile foundation for high-fidelity LES investigations of transcritical injection of sustainable aviation fuels.</div></div>","PeriodicalId":335,"journal":{"name":"International Journal of Heat and Fluid Flow","volume":"121 ","pages":"Article 110644"},"PeriodicalIF":4.8,"publicationDate":"2026-09-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148934368","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Deep learning, as a primary driver of the current technological revolution, is driving AI chips and data center servers toward ultra-high levels of integration. The resulting extreme thermal management challenges have rendered conventional cooling technologies inadequate. To address the critical heat dissipation bottlenecks of high-power AI server chips, this study presents an experimental investigation into the enhanced flow boiling heat transfer characteristics of R134a in microchannels. Microchannel heat sinks with three different channel spacings (0.4, 0.6 and 1.0 mm) were fabricated by a skiving process and then integrated with a novel large area controllable electrodeposited copper porous structure. A heating source with an area as large as 22.5 cm2 was employed in the experiment. Under conditions of 5 K inlet subcooling and a range of low mass fluxes, the dynamic flow pattern transitions, heat transfer coefficients (HTC), critical heat flux (CHF), and pressure drop characteristics were evaluated. Flow visualization confirmed that the electrodeposited micro-nano porous morphology provides high nucleation site density, accelerating phase change, and inducing an earlier transition to the high efficiency annular flow pattern characterized by thin film evaporation. The 0.6 mm porous channel (P-0.6) effectively optimized the balance between bubble expansion and liquid film spreading in the confined space. Attributed to the capillary wicking of the porous network, wall dryout was significantly delayed, enabling superior heat transfer performance across the entire tested mass flux range. At a mass flux of 83.1 kg·m−2·s−1, the P-1 heat sink successfully dissipated an extreme heat load of 1823 W, maintaining the wall temperature below 60℃ and the heater surface temperature below 62°C. While incurring a minor penalty in peak Coefficient of Performance (COP), the P-0.6 configuration significantly extended the high-efficiency operating range and cooling capacity of the heat sink. These findings provide critical data and theoretical guidance for the development of next-generation two-phase cooling architectures for high performance servers.
{"title":"Two phase heat transfer performance of porous microchannels for kilowatt level high performance computing chips","authors":"Yali Feng, Chengyu Hu, Xiaoping Yang, Yonghai Zhang, Zhiwei Zhang, Fangjiao Hui, Jinjia Wei","doi":"10.1016/j.ijheatfluidflow.2026.110651","DOIUrl":"10.1016/j.ijheatfluidflow.2026.110651","url":null,"abstract":"<div><div>Deep learning, as a primary driver of the current technological revolution, is driving AI chips and data center servers toward ultra-high levels of integration. The resulting extreme thermal management challenges have rendered conventional cooling technologies inadequate. To address the critical heat dissipation bottlenecks of high-power AI server chips, this study presents an experimental investigation into the enhanced flow boiling heat transfer characteristics of R134a in microchannels. Microchannel heat sinks with three different channel spacings (0.4, 0.6 and 1.0 mm) were fabricated by a skiving process and then integrated with a novel large area controllable electrodeposited copper porous structure. A heating source with an area as large as 22.5 cm<sup>2</sup> was employed in the experiment. Under conditions of 5 K inlet subcooling and a range of low mass fluxes, the dynamic flow pattern transitions, heat transfer coefficients (HTC), critical heat flux (CHF), and pressure drop characteristics were evaluated. Flow visualization confirmed that the electrodeposited micro-nano porous morphology provides high nucleation site density, accelerating phase change, and inducing an earlier transition to the high efficiency annular flow pattern characterized by thin film evaporation. The 0.6 mm porous channel (P-0.6) effectively optimized the balance between bubble expansion and liquid film spreading in the confined space. Attributed to the capillary wicking of the porous network, wall dryout was significantly delayed, enabling superior heat transfer performance across the entire tested mass flux range. At a mass flux of 83.1 kg·m<sup>−2</sup>·s<sup>−1</sup>, the P-1 heat sink successfully dissipated an extreme heat load of 1823 W, maintaining the wall temperature below 60℃ and the heater surface temperature below 62°C. While incurring a minor penalty in peak Coefficient of Performance (COP), the P-0.6 configuration significantly extended the high-efficiency operating range and cooling capacity of the heat sink. These findings provide critical data and theoretical guidance for the development of next-generation two-phase cooling architectures for high performance servers.</div></div>","PeriodicalId":335,"journal":{"name":"International Journal of Heat and Fluid Flow","volume":"121 ","pages":"Article 110651"},"PeriodicalIF":4.8,"publicationDate":"2026-09-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148934415","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-09-01Epub Date: 2026-08-30DOI: 10.1016/j.ijheatfluidflow.2026.110646
Yaping Sun, Qilong Liu, Shijie Li, Shizhen Qi, Shaohua Han, Runsheng Zhang, Leping Zhou, Li Li, Hui Zhang, Xiaoze Du
This study investigates the fundamental flow physics and thermodynamic irreversibilities governing the performance of bio-inspired, shark-skin rib structures for internal cooling applications. A multi-objective optimization framework integrating Response Surface Methodology (RSM) and a Genetic Algorithm (GA) was employed to optimize four geometric parameters. For turbine blade cooling, a balance was sought between maximizing heat transfer efficiency (Nu/Nu0) and minimizing flow resistance (f/f0) under a Re = 10,000. This methodology enables not only a substantial reduction in computation time but also efficient identification of Pareto-optimal solutions by the RSM-GA framework within the investigated design space. The Pareto-optimal designs show two performance patterns (RA and RB), which are named primarily according to the initials of the ribs and their order. Design RA has a better balance, with an 81% higher thermo-hydraulic performance factor. Design RB maximizes heat transfer (77% gain) but causes excessive pressure loss. Detailed flow diagnostics show that the optimized geometries disrupt coherent streamwise vortices, generating complex secondary flows and multi-scale turbulence that enhance fluid mixing. In contrast to existing rib designs, this study adopts an integrated approach combining multi-objective optimization of bionic ribs with entropy generation analysis. The entropy analysis reveals that the peak entropy production of RB was 67% and 44% higher than that of RO and RA. This work links geometric design to underlying vortex dynamics and second-law efficiency, offering fundamental insights for advanced thermal management.
{"title":"Vortex dynamics and entropy generation in the multi-objective optimization of bio-inspired ribs for internal cooling","authors":"Yaping Sun, Qilong Liu, Shijie Li, Shizhen Qi, Shaohua Han, Runsheng Zhang, Leping Zhou, Li Li, Hui Zhang, Xiaoze Du","doi":"10.1016/j.ijheatfluidflow.2026.110646","DOIUrl":"10.1016/j.ijheatfluidflow.2026.110646","url":null,"abstract":"<div><div>This study investigates the fundamental flow physics and thermodynamic irreversibilities governing the performance of bio-inspired, shark-skin rib structures for internal cooling applications. A multi-objective optimization framework integrating Response Surface Methodology (RSM) and a Genetic Algorithm (GA) was employed to optimize four geometric parameters. For turbine blade cooling, a balance was sought between maximizing heat transfer efficiency (<em>Nu</em>/<em>Nu</em><sub>0</sub>) and minimizing flow resistance (<em>f</em>/<em>f</em><sub>0</sub>) under a <em>Re</em> = 10,000. This methodology enables not only a substantial reduction in computation time but also efficient identification of Pareto-optimal solutions by the RSM-GA framework within the investigated design space. The Pareto-optimal designs show two performance patterns (RA and RB), which are named primarily according to the initials of the ribs and their order. Design RA has a better balance, with an 81% higher thermo-hydraulic performance factor. Design RB maximizes heat transfer (77% gain) but causes excessive pressure loss. Detailed flow diagnostics show that the optimized geometries disrupt coherent streamwise vortices, generating complex secondary flows and multi-scale turbulence that enhance fluid mixing. In contrast to existing rib designs, this study adopts an integrated approach combining multi-objective optimization of bionic ribs with entropy generation analysis. The entropy analysis reveals that the peak entropy production of RB was 67% and 44% higher than that of RO and RA. This work links geometric design to underlying vortex dynamics and second-law efficiency, offering fundamental insights for advanced thermal management.</div></div>","PeriodicalId":335,"journal":{"name":"International Journal of Heat and Fluid Flow","volume":"121 ","pages":"Article 110646"},"PeriodicalIF":4.8,"publicationDate":"2026-09-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148934401","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-09-01Epub Date: 2026-08-31DOI: 10.1016/j.ijheatfluidflow.2026.110643
Wenchang Wu, Xiaoliang Deng, Xingsi Han
The study seeks an effective high-fidelity numerical simulation method for compressible turbulent flows. To this end, an efficient unified seamless hybrid Reynolds-Averaged Navier-Stokes/Large eddy simulation (RANS/LES) method coupled with high-order finite differencing scheme is implemented and evaluated. The high-order PHengLEI flow solver is used to perform the numerical simulations. To obtain the high resolution of multiscale turbulence, a newly developed six-point multi-order candidate (MOC) nonlinear weighting with an adaptive Riemann solver is applied within the framework of high-order weighted compact nonlinear scheme (WCNS). To be applicable to compressible separated flows with affordable computation grids, a two-layer self-adaptive model is developed to improve the original Self-Adaptive Turbulence Eddy Simulation (SATES) method, resulting in the SATES-II model, in which the filter length scale has a dynamic coefficient coupling local turbulent flow features. The proposed methods are validated through supersonic base flow and a multi-element airfoil at high Reynolds number. Results demonstrate that the SATES-II model with adaptive filter length scale outperforms both the previous SATES-SM model and the IDDES model. Furthermore, the adaptive filter length scale in the present SATES-II model is extended to formulate a new-version of DES method, i.e. AIDDES model, which shows notable improvements.
{"title":"Hybrid RANS/LES simulations of compressible turbulent flows using high-order finite differencing method","authors":"Wenchang Wu, Xiaoliang Deng, Xingsi Han","doi":"10.1016/j.ijheatfluidflow.2026.110643","DOIUrl":"10.1016/j.ijheatfluidflow.2026.110643","url":null,"abstract":"<div><div>The study seeks an effective high-fidelity numerical simulation method for compressible turbulent flows. To this end, an efficient unified seamless hybrid Reynolds-Averaged Navier-Stokes/Large eddy simulation (RANS/LES) method coupled with high-order finite differencing scheme is implemented and evaluated. The high-order PHengLEI flow solver is used to perform the numerical simulations. To obtain the high resolution of multiscale turbulence, a newly developed<!--> <!-->six-point multi-order candidate (MOC) nonlinear weighting with an adaptive Riemann solver is applied within the framework of high-order weighted compact nonlinear scheme (WCNS). To be applicable to compressible separated flows with affordable computation grids, a two-layer self-adaptive model is developed to improve the original Self-Adaptive Turbulence Eddy Simulation (SATES) method, resulting in the SATES-II model, in which the filter length scale has a dynamic coefficient coupling local turbulent flow features. The proposed methods are validated through supersonic base flow and a multi-element airfoil at high Reynolds number. Results demonstrate<!--> <!-->that the SATES-II model with adaptive filter length scale outperforms<!--> <!-->both the previous SATES-SM model and the IDDES model. Furthermore, the adaptive filter length scale in the present SATES-II model is extended to formulate a new-version of DES method, i.e. AIDDES model, which shows notable improvements.</div></div>","PeriodicalId":335,"journal":{"name":"International Journal of Heat and Fluid Flow","volume":"121 ","pages":"Article 110643"},"PeriodicalIF":4.8,"publicationDate":"2026-09-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148934494","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-09-01Epub Date: 2026-05-05DOI: 10.1016/j.ijheatfluidflow.2026.110449
Tatsuya Yonemura, Takashi Ohta
The flow of a cryogenic fluid such as liquid hydrogen can become turbulent, and phase changes can occur owing to variations in wall temperatures, resulting in significant alterations in the flow field. Direct numerical simulations were conducted to investigate turbulent liquid hydrogen with near-wall boiling in a parallel-plate channel. The simulations showed that heating one wall induced boiling within fully developed wall turbulence. The turbulence statistics and the relationship between near-wall turbulence structures and phase change were analyzed in detail. The results demonstrate that boiling phenomena accelerate wall-normal and spanwise velocities that promote the sweep and ejection near the wall, temporarily increasing the turbulent kinetic energy by approximately 1.5 times. Additionally, turbulence influences boiling after an initial transient period, thereby enhancing phase-change. The turbulence statistics indicate that boiling modulates wall turbulence toward isotropy. The local phase change also increases the streamwise velocity near the wall, which in turn reduces the velocity gradient, thereby decreasing the turbulent kinetic energy. Furthermore, the combined effects of these mechanisms reveal a turbulence modulation due to the interaction between the near-wall turbulence structure and the phase change, characterized by turbulence enhancement, followed by its suppression over time. Accurately predicting these phenomena using statistical modeling techniques that treat the turbulence structure as the Reynolds stress is fundamentally challenging, and an appropriate turbulence model is required to reproduce the phenomena observed in this study.
{"title":"Direct numerical simulation of turbulent flow of liquid hydrogen with wall boiling","authors":"Tatsuya Yonemura, Takashi Ohta","doi":"10.1016/j.ijheatfluidflow.2026.110449","DOIUrl":"10.1016/j.ijheatfluidflow.2026.110449","url":null,"abstract":"<div><div>The flow of a cryogenic fluid such as liquid hydrogen can become turbulent, and phase changes can occur owing to variations in wall temperatures, resulting in significant alterations in the flow field. Direct numerical simulations were conducted to investigate turbulent liquid hydrogen with near-wall boiling in a parallel-plate channel. The simulations showed that heating one wall induced boiling within fully developed wall turbulence. The turbulence statistics and the relationship between near-wall turbulence structures and phase change were analyzed in detail. The results demonstrate that boiling phenomena accelerate wall-normal and spanwise velocities that promote the sweep and ejection near the wall, temporarily increasing the turbulent kinetic energy by approximately 1.5 times. Additionally, turbulence influences boiling after an initial transient period, thereby enhancing phase-change. The turbulence statistics indicate that boiling modulates wall turbulence toward isotropy. The local phase change also increases the streamwise velocity near the wall, which in turn reduces the velocity gradient, thereby decreasing the turbulent kinetic energy. Furthermore, the combined effects of these mechanisms reveal a turbulence modulation due to the interaction between the near-wall turbulence structure and the phase change, characterized by turbulence enhancement, followed by its suppression over time. Accurately predicting these phenomena using statistical modeling techniques that treat the turbulence structure as the Reynolds stress is fundamentally challenging, and an appropriate turbulence model is required to reproduce the phenomena observed in this study.</div></div>","PeriodicalId":335,"journal":{"name":"International Journal of Heat and Fluid Flow","volume":"121 ","pages":"Article 110449"},"PeriodicalIF":4.8,"publicationDate":"2026-09-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148958107","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-09-01Epub Date: 2026-06-12DOI: 10.1016/j.ijheatfluidflow.2026.110505
Lin Chen, Atsuki Komiya, Bengt Sunden, Xiaohu Yang
{"title":"Editorial for the special issue: Thermos-fluidics in low-carbon/green power conversion systems","authors":"Lin Chen, Atsuki Komiya, Bengt Sunden, Xiaohu Yang","doi":"10.1016/j.ijheatfluidflow.2026.110505","DOIUrl":"10.1016/j.ijheatfluidflow.2026.110505","url":null,"abstract":"","PeriodicalId":335,"journal":{"name":"International Journal of Heat and Fluid Flow","volume":"121 ","pages":"Article 110505"},"PeriodicalIF":4.8,"publicationDate":"2026-09-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148958108","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}