Pub Date : 2025-12-15DOI: 10.1016/j.ijheatfluidflow.2025.110204
Pavel Panek , Davy Brouzet , Mohsen Talei
The ability of large-eddy simulation (LES) to capture premixed flame acoustics is studied for a turbulent, premixed, round jet flame. Central finite differencing schemes are used with artificial damping. The flame surface density (FSD) approach is used as the combustion model for the reaction progress variable. The effects of combustion modelling and the spatial discretisation scheme order on the flame shape and the sound pressure level (SPL) spectra are examined. It is shown that the FSD approach does not correctly capture flame annihilation. This is where the progress variable gradient decreases as two flame surfaces approach each other. A significant impact of FSD modelling on the SPL is also observed. Comparing the SPL spectra obtained with FSD modelling with those obtained with single-step, Arrhenius-rate chemistry on a fine grid shows that the SPL is underestimated below and overestimated above , where is a non-dimensional frequency based on the jet diameter and mean centreline jet velocity known as the Strouhal number. Coarse grid simulations approximately match the fine grid results at frequencies up to but overestimate the SPL even more strongly above this level. It is concluded that FSD modelling needs further adjustments to accurately capture combustion noise.
{"title":"A posteriori analysis of flame surface density modelling for combustion noise prediction in premixed flames","authors":"Pavel Panek , Davy Brouzet , Mohsen Talei","doi":"10.1016/j.ijheatfluidflow.2025.110204","DOIUrl":"10.1016/j.ijheatfluidflow.2025.110204","url":null,"abstract":"<div><div>The ability of large-eddy simulation (LES) to capture premixed flame acoustics is studied for a turbulent, premixed, round jet flame. Central finite differencing schemes are used with artificial damping. The flame surface density (FSD) approach is used as the combustion model for the reaction progress variable. The effects of combustion modelling and the spatial discretisation scheme order on the flame shape and the sound pressure level (SPL) spectra are examined. It is shown that the FSD approach does not correctly capture flame annihilation. This is where the progress variable gradient decreases as two flame surfaces approach each other. A significant impact of FSD modelling on the SPL is also observed. Comparing the SPL spectra obtained with FSD modelling with those obtained with single-step, Arrhenius-rate chemistry on a fine grid shows that the SPL is underestimated below <span><math><mrow><mi>S</mi><mi>t</mi><mo>≈</mo><mn>4</mn></mrow></math></span> and overestimated above <span><math><mrow><mi>S</mi><mi>t</mi><mo>≈</mo><mn>4</mn></mrow></math></span>, where <span><math><mrow><mi>S</mi><mi>t</mi></mrow></math></span> is a non-dimensional frequency based on the jet diameter and mean centreline jet velocity known as the Strouhal number. Coarse grid simulations approximately match the fine grid results at frequencies up to <span><math><mrow><mi>S</mi><mi>t</mi><mo>≈</mo><mn>3</mn></mrow></math></span> but overestimate the SPL even more strongly above this level. It is concluded that FSD modelling needs further adjustments to accurately capture combustion noise.</div></div>","PeriodicalId":335,"journal":{"name":"International Journal of Heat and Fluid Flow","volume":"118 ","pages":"Article 110204"},"PeriodicalIF":2.6,"publicationDate":"2025-12-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145797334","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 : 2025-12-14DOI: 10.1016/j.ijheatfluidflow.2025.110202
Hyowon Bang , Huichan Lee , Giyoung Park , Sungjae Kim , Seangwock Lee
As aerospace propulsion steadily transitions toward electrification, effectively dissipating heat in compact electric tail rotor motors has become a pressing design challenge. In this study, a high-fidelity computational fluid dynamics (CFD) framework is employed to investigate an air-cooled propulsion system specifically configured for helicopter tail rotor applications. Three passive geometric modifications are examined in detail: widening of lateral duct openings, introducing a 15° inclination to the duct sidewalls, and incorporating rotating blades along the central shaft. The optimized configuration demonstrated substantial thermal gains, reducing peak temperatures by 10.38% in the winding, 8.23% in the magnet, and 19.0% in the shaft. This enhancement was primarily attributed to enhanced airflow penetration, boundary-layer thinning, and the development of favorable secondary recirculation zones. Convective heat transfer performance was also elevated, with Nusselt numbers increasing by up to 11.6% across key components. Importantly, these thermal benefits were achieved while maintaining aerodynamic integrity, as the simulated thrust deviated by only 1.34% from experimental rig measurements. The results highlight how targeted, flow-physics-based passive design strategies can simultaneously improve cooling performance and preserve propulsion efficiency. This approach offers a scalable and integration-friendly pathway for next-generation electric rotorcraft requiring lightweight, thermally resilient tail rotor systems.
{"title":"CFD investigation of thermal management and flow behavior in an optimized air-cooled tail rotor propulsion system","authors":"Hyowon Bang , Huichan Lee , Giyoung Park , Sungjae Kim , Seangwock Lee","doi":"10.1016/j.ijheatfluidflow.2025.110202","DOIUrl":"10.1016/j.ijheatfluidflow.2025.110202","url":null,"abstract":"<div><div>As aerospace propulsion steadily transitions toward electrification, effectively dissipating heat in compact electric tail rotor motors has become a pressing design challenge. In this study, a high-fidelity computational fluid dynamics (CFD) framework is employed to investigate an air-cooled propulsion system specifically configured for helicopter tail rotor applications. Three passive geometric modifications are examined in detail: widening of lateral duct openings, introducing a 15° inclination to the duct sidewalls, and incorporating rotating blades along the central shaft. The optimized configuration demonstrated substantial thermal gains, reducing peak temperatures by 10.38% in the winding, 8.23% in the magnet, and 19.0% in the shaft. This enhancement was primarily attributed to enhanced airflow penetration, boundary-layer thinning, and the development of favorable secondary recirculation zones. Convective heat transfer performance was also elevated, with Nusselt numbers increasing by up to 11.6% across key components. Importantly, these thermal benefits were achieved while maintaining aerodynamic integrity, as the simulated thrust deviated by only 1.34% from experimental rig measurements. The results highlight how targeted, flow-physics-based passive design strategies can simultaneously improve cooling performance and preserve propulsion efficiency. This approach offers a scalable and integration-friendly pathway for next-generation electric rotorcraft requiring lightweight, thermally resilient tail rotor systems.</div></div>","PeriodicalId":335,"journal":{"name":"International Journal of Heat and Fluid Flow","volume":"118 ","pages":"Article 110202"},"PeriodicalIF":2.6,"publicationDate":"2025-12-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145797236","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}
Wind energy has become increasingly crucial among renewable energy resources, and as the prerequisite for effective wind energy utilization, accurate wind resource assessments are hence prominent. Topography significantly affects near-surface wind fields and is thus an essential factor to be considered for wind resource assessment. However, few studies have addressed the wind resource distribution downwind of spanwise-heterogeneous terrains, which are frequently encountered in both natural and urban environments. To fill this gap, large eddy simulations were performed to investigate wind distributions over a uniform terrain downwind of terrains featuring spanwise-alternating strips of low and high surface roughness, and the high roughness surface coverage ratio is varied from 12.5% to 100%. The analysis focuses on the wind speed within the turbine blades region, and two dominant mechanisms are identified: the rough-to-smooth transition (RST) and spanwise heterogeneity of momentum (SHM). The dominance of the two mechanisms depends strongly on , with RST dominating over SHM as increases. Moreover, RST initially increases the wind speed due to the growth of the internal boundary layer but leads to a decrease of wind speed further downwind as the internal boundary layer rearranges. SHM results from secondary flows induced by the heterogeneous terrain upwind, which decay downwind. As increases, secondary flow structures transition from asymmetric to symmetric patterns, and lower leads to asymmetric flows, sustaining momentum heterogeneity farther downwind. Hence, upwind spanwise heterogeneous terrains can significantly impact the wind speed distribution, and suggestions are made for wind farm planning based on the results.
{"title":"Effects of the upwind spanwise heterogeneous terrain on the wind resource distribution","authors":"Haosen H.A. Xu , Tianxiang Yu , Raúl Bayoán Cal , Xiaowei Zhu","doi":"10.1016/j.ijheatfluidflow.2025.110182","DOIUrl":"10.1016/j.ijheatfluidflow.2025.110182","url":null,"abstract":"<div><div>Wind energy has become increasingly crucial among renewable energy resources, and as the prerequisite for effective wind energy utilization, accurate wind resource assessments are hence prominent. Topography significantly affects near-surface wind fields and is thus an essential factor to be considered for wind resource assessment. However, few studies have addressed the wind resource distribution downwind of spanwise-heterogeneous terrains, which are frequently encountered in both natural and urban environments. To fill this gap, large eddy simulations were performed to investigate wind distributions over a uniform terrain downwind of terrains featuring spanwise-alternating strips of low and high surface roughness, and the high roughness surface coverage ratio <span><math><mi>λ</mi></math></span> is varied from 12.5% to 100%. The analysis focuses on the wind speed within the turbine blades region, and two dominant mechanisms are identified: the rough-to-smooth transition (RST) and spanwise heterogeneity of momentum (SHM). The dominance of the two mechanisms depends strongly on <span><math><mi>λ</mi></math></span>, with RST dominating over SHM as <span><math><mi>λ</mi></math></span> increases. Moreover, RST initially increases the wind speed due to the growth of the internal boundary layer but leads to a decrease of wind speed further downwind as the internal boundary layer rearranges. SHM results from secondary flows induced by the heterogeneous terrain upwind, which decay downwind. As <span><math><mi>λ</mi></math></span> increases, secondary flow structures transition from asymmetric to symmetric patterns, and lower <span><math><mi>λ</mi></math></span> leads to asymmetric flows, sustaining momentum heterogeneity farther downwind. Hence, upwind spanwise heterogeneous terrains can significantly impact the wind speed distribution, and suggestions are made for wind farm planning based on the results.</div></div>","PeriodicalId":335,"journal":{"name":"International Journal of Heat and Fluid Flow","volume":"118 ","pages":"Article 110182"},"PeriodicalIF":2.6,"publicationDate":"2025-12-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145797335","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 : 2025-12-14DOI: 10.1016/j.ijheatfluidflow.2025.110199
Joseph Farrell , Evan T. Hurlburt , Arganthaël Berson , Allison Mahvi
Nucleate boiling in thin film annular flow is a complex phenomenon influenced by a variety of heat transfer mechanisms. While the fundamentals are rooted in the comparatively well-understood fields of pool and flow boiling, thin film flow boiling introduces additional challenges that make the heat transfer behavior hard to predict. Specifically, the relative importance of the heat transfer mechanisms controlling nucleate boiling in thin annular films and their roles in flow regime transitions remains inadequately understood. In light of this, the foundational aspects of nucleate boiling and the relationship between flow parameters – such as liquid film thickness, system pressure, surface tension, channel geometry, mass flow rate, and applied heat flux – and the intensity of nucleate boiling in thin annular films are analyzed. The importance of these parameters is demonstrated through the presentation of multiple correlations, with a critical examination of the limitations inherent in correlation-based modeling. Additionally, recent research has identified disturbance waves as a significant factor in enhancing bubble nucleation, yet the underlying mechanisms driving this phenomenon remain undefined. Therefore, the insights and deficiencies of three theories seeking to explain wave-based nucleation are extensively considered. This research aims to introduce and categorize the wide breadth of literature regarding thin film annular flow nucleate boiling to identify gaps in understanding and facilitate future physics-based modeling efforts.
{"title":"A comprehensive review of nucleate boiling fundamentals applied to thin film annular flow","authors":"Joseph Farrell , Evan T. Hurlburt , Arganthaël Berson , Allison Mahvi","doi":"10.1016/j.ijheatfluidflow.2025.110199","DOIUrl":"10.1016/j.ijheatfluidflow.2025.110199","url":null,"abstract":"<div><div>Nucleate boiling in thin film annular flow is a complex phenomenon influenced by a variety of heat transfer mechanisms. While the fundamentals are rooted in the comparatively well-understood fields of pool and flow boiling, thin film flow boiling introduces additional challenges that make the heat transfer behavior hard to predict. Specifically, the relative importance of the heat transfer mechanisms controlling nucleate boiling in thin annular films and their roles in flow regime transitions remains inadequately understood. In light of this, the foundational aspects of nucleate boiling and the relationship between flow parameters – such as liquid film thickness, system pressure, surface tension, channel geometry, mass flow rate, and applied heat flux – and the intensity of nucleate boiling in thin annular films are analyzed. The importance of these parameters is demonstrated through the presentation of multiple correlations, with a critical examination of the limitations inherent in correlation-based modeling. Additionally, recent research has identified disturbance waves as a significant factor in enhancing bubble nucleation, yet the underlying mechanisms driving this phenomenon remain undefined. Therefore, the insights and deficiencies of three theories seeking to explain wave-based nucleation are extensively considered. This research aims to introduce and categorize the wide breadth of literature regarding thin film annular flow nucleate boiling to identify gaps in understanding and facilitate future physics-based modeling efforts.</div></div>","PeriodicalId":335,"journal":{"name":"International Journal of Heat and Fluid Flow","volume":"118 ","pages":"Article 110199"},"PeriodicalIF":2.6,"publicationDate":"2025-12-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145797421","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 : 2025-12-12DOI: 10.1016/j.ijheatfluidflow.2025.110192
W. Lu , W.K. Lam , T. Zahtila
This study presents a simple autoencoder-based bi-fidelity method, where the latent space is constrained to be invariant between both low- and high-fidelity datasets. The method is applied to the flow past a cylinder confined in a duct, where model fidelity is governed by computational model resolution. Significant cost savings are observed without significant compromise to approximation accuracy. The influence of several hyperparameters are considered including: latent space dimension, sample number and distribution. A method of generating synthetic solutions is also explored. This is done by performing regression on the common latent space of low- and high-fidelity models, which is then passed through the high-fidelity decoder. As a large number of samples are available from the low-fidelity model, significant improvements are observed compared with the traditional method of interpolating only high-fidelity data.
{"title":"An autoencoder-based bi-fidelity method utilising frozen latent spaces. Application to flow past a confined cylinder","authors":"W. Lu , W.K. Lam , T. Zahtila","doi":"10.1016/j.ijheatfluidflow.2025.110192","DOIUrl":"10.1016/j.ijheatfluidflow.2025.110192","url":null,"abstract":"<div><div>This study presents a simple autoencoder-based bi-fidelity method, where the latent space is constrained to be invariant between both low- and high-fidelity datasets. The method is applied to the flow past a cylinder confined in a duct, where model fidelity is governed by computational model resolution. Significant cost savings are observed without significant compromise to approximation accuracy. The influence of several hyperparameters are considered including: latent space dimension, sample number and distribution. A method of generating synthetic solutions is also explored. This is done by performing regression on the common latent space of low- and high-fidelity models, which is then passed through the high-fidelity decoder. As a large number of samples are available from the low-fidelity model, significant improvements are observed compared with the traditional method of interpolating only high-fidelity data.</div></div>","PeriodicalId":335,"journal":{"name":"International Journal of Heat and Fluid Flow","volume":"118 ","pages":"Article 110192"},"PeriodicalIF":2.6,"publicationDate":"2025-12-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145748815","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 : 2025-12-11DOI: 10.1016/j.ijheatfluidflow.2025.110189
Huan Liu , Pavel Skripov , Lin Chen
The supercritical fluid region is generally considered a complex phase region, which has a lot on the thermodynamic transitions and fluid transport behaviors of a system. The present study introduces a phase-field method using density as the order parameter, which is combined with a free energy minimization technique to predict the dynamic responses of supercritical CO2 under parameter equilibrium procedures. During this phase, field simulation analysis, cases from the close-to-critical region, liquid region, and low-density fluid region have been examined, which procedures then show a system equilibrium stage from ‘quasi-uniform’ state to an ‘overall equilibrium’ state. The simulation process follows a continuous and gradual change in macroscopic properties rather than the abrupt transition observed in the conventional phase change process. By statistically tracking the temporal evolution of ‘pseudo-liquid regions’, the local non-uniformities emergence process during the formation of supercritical phases has been revealed, which in turn gives the dynamic development of spatial heterogeneity. The extrema of density fluctuations form ridge lines near the pseudo-critical line, with a maximum spatial fluctuation error of 0.06% compared to the pseudo-critical line. The boundary identified with intense density fluctuations resembles the dynamic transition across the pseudo-critical triangle.
{"title":"Preliminary study on the dynamic equilibrium process of supercritical region CO2 fluid by phase field method","authors":"Huan Liu , Pavel Skripov , Lin Chen","doi":"10.1016/j.ijheatfluidflow.2025.110189","DOIUrl":"10.1016/j.ijheatfluidflow.2025.110189","url":null,"abstract":"<div><div>The supercritical fluid region is generally considered a complex phase region, which has a lot on the thermodynamic transitions and fluid transport behaviors of a system. The present study introduces a phase-field method using density as the order parameter, which is combined with a free energy minimization technique to predict the dynamic responses of supercritical CO<sub>2</sub> under parameter equilibrium procedures. During this phase, field simulation analysis, cases from the close-to-critical region, liquid region, and low-density fluid region have been examined, which procedures then show a system equilibrium stage from ‘quasi-uniform’ state to an ‘overall equilibrium’ state. The simulation process follows a continuous and gradual change in macroscopic properties rather than the abrupt transition observed in the conventional phase change process. By statistically tracking the temporal evolution of ‘pseudo-liquid regions’, the local non-uniformities emergence process during the formation of supercritical phases has been revealed, which in turn gives the dynamic development of spatial heterogeneity. The extrema of density fluctuations form ridge lines near the pseudo-critical line, with a maximum spatial fluctuation error of 0.06% compared to the pseudo-critical line. The boundary identified with intense density fluctuations resembles the dynamic transition across the pseudo-critical triangle.</div></div>","PeriodicalId":335,"journal":{"name":"International Journal of Heat and Fluid Flow","volume":"118 ","pages":"Article 110189"},"PeriodicalIF":2.6,"publicationDate":"2025-12-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145748811","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 : 2025-12-10DOI: 10.1016/j.ijheatfluidflow.2025.110190
Anthuvan Ezhilarasi P., Dhivya Mohanavel
Arterial stenosis causes the arteries to narrow, a change that can significantly impact blood flow. Identifying the conditions related to arterial blockage is crucial for prompt diagnosis and treatment. Keeping this fact in mind, the purpose of this study is to analyse the heat and mass transfer properties of a Casson hybrid blood flow in a constricted artery inhabited by oxytactic microorganisms, and also by employing the Arrhenius energy. The mathematical assessment of copper and titanium oxide in the stenosis region, where blood circulates, considers temperature-dependent viscosity and thermal conductivity. The hybrid (Cu/TiO) nanoparticles used in the bloodstream synergise to enhance heat and mass transfer, offers antithrombotic and antimicrobial benefits which promote better vascular health and prevent infection. The governing coupled partial differential equations are converted into a system of nonlinear ordinary differential equations via the similarity transformation technique, and the resulting ODE is then numerically solved using the MATLAB solver. The influences of various factors, including variable viscosity, variable thermal conductivity, activation energy, chemical reaction, bioconvection Schmidt number, and Peclet number, have been examined in equations for velocity, temperature, concentration, and microorganisms. The primary findings indicate that varying the viscosity parameter slows the flow rate, whereas raising the values of thermal conductivity parameter makes heat transmission more reliable. Additionally, the reduction in motile density profiles was attributed to both Lewis number and bioconvection Schmidt number. Furthermore, this study was compared to previous results, showing a high degree of correspondence. The sensitivity analysis used in this model examines the impact of different factors on the effectiveness of convective and diffusive transport processes in mass transfer. This research can improve arterial disease diagnostic techniques by simulating intricate blood flow scenarios.
{"title":"Mathematical modelling and heat-mass transport analysis of a bioconvective Casson hybrid nanofluid in a stenosed artery","authors":"Anthuvan Ezhilarasi P., Dhivya Mohanavel","doi":"10.1016/j.ijheatfluidflow.2025.110190","DOIUrl":"10.1016/j.ijheatfluidflow.2025.110190","url":null,"abstract":"<div><div>Arterial stenosis causes the arteries to narrow, a change that can significantly impact blood flow. Identifying the conditions related to arterial blockage is crucial for prompt diagnosis and treatment. Keeping this fact in mind, the purpose of this study is to analyse the heat and mass transfer properties of a Casson hybrid blood flow in a constricted artery inhabited by oxytactic microorganisms, and also by employing the Arrhenius energy. The mathematical assessment of copper and titanium oxide in the stenosis region, where blood circulates, considers temperature-dependent viscosity and thermal conductivity. The hybrid (Cu/TiO<span><math><msub><mrow></mrow><mrow><mn>2</mn></mrow></msub></math></span>) nanoparticles used in the bloodstream synergise to enhance heat and mass transfer, offers antithrombotic and antimicrobial benefits which promote better vascular health and prevent infection. The governing coupled partial differential equations are converted into a system of nonlinear ordinary differential equations via the similarity transformation technique, and the resulting ODE is then numerically solved using the MATLAB solver. The influences of various factors, including variable viscosity, variable thermal conductivity, activation energy, chemical reaction, bioconvection Schmidt number, and Peclet number, have been examined in equations for velocity, temperature, concentration, and microorganisms. The primary findings indicate that varying the viscosity parameter slows the flow rate, whereas raising the values of thermal conductivity parameter makes heat transmission more reliable. Additionally, the reduction in motile density profiles was attributed to both Lewis number and bioconvection Schmidt number. Furthermore, this study was compared to previous results, showing a high degree of correspondence. The sensitivity analysis used in this model examines the impact of different factors on the effectiveness of convective and diffusive transport processes in mass transfer. This research can improve arterial disease diagnostic techniques by simulating intricate blood flow scenarios.</div></div>","PeriodicalId":335,"journal":{"name":"International Journal of Heat and Fluid Flow","volume":"118 ","pages":"Article 110190"},"PeriodicalIF":2.6,"publicationDate":"2025-12-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145748814","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 : 2025-12-10DOI: 10.1016/j.ijheatfluidflow.2025.110178
Md. Mahbub Alam , Zia Ullah , Hongjun Zhu , Chunning Ji , Md. Islam , Mostafa Zeinoddini
This study numerically investigates the fluid flow, aerodynamic forces, and heat transfer characteristics around a circular cylinder (diameter D) positioned in the wake of another cylinder (diameter d), with varying diameter ratios d/D (= 0.4, 0.6, 0.8, and 1.0) and inter-cylinder spacing ratios L* (= L/D = 1.2 –5.0). Significant effects of d/D and L* are identified on forces, heat transfer, and fluid dynamics at ReD = 200. For a given L*, time-mean drag, fluctuating lift and drag, and Nusselt number all are generally smaller for larger d/D, with the time-mean drag force even becoming negative for d/D = 1.0. In the reattachment flow regime, the time-mean drag and Nusselt number decrease with increasing L* for smaller d/D = 0.4 and 0.6 but increase for larger d/D = 0.8 and 1.0. In the coshedding flow regime, the Nusselt number and forces increase with L*, regardless of d/D. The critical L* separating the reattachment and coshedding flows increases with increasing d/D from 0.4 to 0.8 before declining from d/D = 0.8 to 1.0. Streamwise velocity fields and local Nusselt number distributions on the cylinder surface are analyzed to explain the observed heat transfer phenomena. Local heat transfer is maximal on the front surface of the cylinder, corresponding to the inflection point on the near-surface velocity profile.
本文通过数值模拟研究了不同直径比D /D(= 0.4、0.6、0.8和1.0)和缸间间距比L* (= L/D = 1.2 -5.0)下,一个圆柱体(直径D)与另一个圆柱体(直径D)的尾迹之间的流体流动、气动力和传热特性。在ReD = 200时,确定了d/ d和L*对力、传热和流体动力学的显著影响。对于给定的L*,当d/ d较大时,时间平均阻力、波动升力和阻力以及努塞尔数一般都较小,当d/ d = 1.0时,时间平均阻力甚至变为负值。当d/ d = 0.4和0.6时,时间平均阻力和Nusselt数随L*的增大而减小,当d/ d = 0.8和1.0时,时间平均阻力和Nusselt数随L*的增大而增大。在共脱落流态中,努塞尔数和力随L*的增加而增加,与d/ d无关。分离重附流和共脱落流的临界L*随着d/ d的增加从0.4增大到0.8,然后从d/ d = 0.8减小到1.0。分析了圆柱表面的流向速度场和局部努塞尔数分布来解释观察到的换热现象。局部换热在柱体前表面最大,与近表面速度剖面的拐点相对应。
{"title":"Impact of cylinder diameter and spacing on fluid flow, forces, and heat transfer in tandem cylinder configuration","authors":"Md. Mahbub Alam , Zia Ullah , Hongjun Zhu , Chunning Ji , Md. Islam , Mostafa Zeinoddini","doi":"10.1016/j.ijheatfluidflow.2025.110178","DOIUrl":"10.1016/j.ijheatfluidflow.2025.110178","url":null,"abstract":"<div><div>This study numerically investigates the fluid flow, aerodynamic forces, and heat transfer characteristics around a circular cylinder (diameter <em>D</em>) positioned in the wake of another cylinder (diameter <em>d</em>), with varying diameter ratios <em>d/D</em> (= 0.4, 0.6, 0.8, and 1.0) and inter-cylinder spacing ratios <em>L*</em> (= <em>L/D</em> = 1.2 –5.0). Significant effects of <em>d/D</em> and <em>L*</em> are identified on forces, heat transfer, and fluid dynamics at <em>Re<sub>D</sub></em> = 200. For a given <em>L*</em>, time-mean drag, fluctuating lift and drag, and Nusselt number all are generally smaller for larger <em>d/D</em>, with the time-mean drag force even becoming negative for <em>d/D</em> = 1.0. In the reattachment flow regime, the time-mean drag and Nusselt number decrease with increasing <em>L*</em> for smaller <em>d/D</em> = 0.4 and 0.6 but increase for larger <em>d/D</em> = 0.8 and 1.0. In the coshedding flow regime, the Nusselt number and forces increase with <em>L*</em>, regardless of <em>d/D</em>. The critical <em>L*</em> separating the reattachment and coshedding flows increases with increasing <em>d/D</em> from 0.4 to 0.8 before declining from <em>d/D</em> = 0.8 to 1.0. Streamwise velocity fields and local Nusselt number distributions on the cylinder surface are analyzed to explain the observed heat transfer phenomena. Local heat transfer is maximal on the front surface of the cylinder, corresponding to the inflection point on the near-surface velocity profile.</div></div>","PeriodicalId":335,"journal":{"name":"International Journal of Heat and Fluid Flow","volume":"118 ","pages":"Article 110178"},"PeriodicalIF":2.6,"publicationDate":"2025-12-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145748810","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}
Optimal thermal regulation in battery modules is critical for maintaining the efficient and reliable operation of battery packs. This study designed an active immersion cooling system for prismatic lithium-ion battery modules, demonstrating maximum temperature reductions of 34 % and 47.7 % relative to static submersion cooling and free convection methods, respectively. The impacts of coolant flow rate, flow direction, and module arrangement patterns on forced-flow immersion cooling performance were numerically investigated. Results indicate that parallel battery arrangements show superior cooling performance compared to staggered configurations. Among the five flow patterns (top to bottom, bottom to top, top to top, bottom to bottom, and center to center), the top to bottom layout exhibits the optimal cooling efficiency. As the inlet flow rate increases, both the maximum battery temperature and the temperature difference across the battery pack first drop sharply and then slowly after the flow rate reaches 0.023 kg/s. Power consumption demonstrates a positive relationship with inlet velocity, whereas the cooling index shows an inverse relationship. As the horizontal spacing between batteries varies from 1 to 6 mm, the maximum temperature of the battery and temperature difference show a “U” shaped trend, the power consumption decreases monotonically, while the cooling index shows a “И” shaped variation trend. Similarly, with the increase of longitudinal distance between batteries, Both the peak battery temperature and temperature differential present a “U” shaped variation trend, and the power dissipation decreases monotonically. Conversely, the cooling index increases monotonously. Finally, a 4-mm horizontal and 5-mm longitudinal spacing were identified as the optimal configuration. This study addressed both operational safety and thermal management efficiency for prismatic lithium-ion batteries, and established design guidelines for high-performance immersion cooling systems.
{"title":"Design and performance optimization of liquid immersion cooling system for prismatic lithium-ion battery modules","authors":"Luyao Zhao, Jiafeng Wang, Minxue Zheng, Mingyi Chen","doi":"10.1016/j.ijheatfluidflow.2025.110181","DOIUrl":"10.1016/j.ijheatfluidflow.2025.110181","url":null,"abstract":"<div><div>Optimal thermal regulation in battery modules<!--> <!-->is critical for maintaining<!--> <!-->the efficient and reliable operation of battery packs. This study designed an active immersion<!--> <!-->cooling system for prismatic lithium-ion battery modules, demonstrating maximum temperature reductions of 34 % and 47.7 % relative to static submersion cooling and free convection methods, respectively. The impacts of coolant flow rate, flow direction, and module arrangement patterns on forced-flow immersion cooling performance were numerically investigated. Results indicate that parallel battery arrangements show superior cooling performance compared to staggered configurations. Among the five flow patterns (top to bottom, bottom to top, top to top, bottom to bottom, and center to center), the top to bottom layout exhibits the optimal cooling efficiency. As the inlet flow rate increases, both the maximum battery temperature and the temperature difference across the battery pack first drop sharply and then slowly after the flow rate reaches 0.023 kg/s. Power consumption demonstrates a positive relationship with inlet velocity, whereas the cooling index shows an inverse relationship. As the horizontal spacing between batteries varies from<!--> <!-->1 to 6 mm, the maximum temperature of the battery and temperature difference show a “U” shaped trend, the power consumption decreases monotonically, while the cooling index shows a “И” shaped variation trend. Similarly, with the increase of longitudinal distance between batteries, Both the peak battery temperature and temperature differential present a “U” shaped variation trend, and the power dissipation decreases monotonically. Conversely, the cooling index increases monotonously. Finally, a 4-mm horizontal and 5-mm longitudinal spacing were identified as the optimal configuration. This study addressed both operational safety and thermal management efficiency for prismatic lithium-ion batteries, and established design guidelines for high-performance immersion cooling systems.</div></div>","PeriodicalId":335,"journal":{"name":"International Journal of Heat and Fluid Flow","volume":"118 ","pages":"Article 110181"},"PeriodicalIF":2.6,"publicationDate":"2025-12-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145748812","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 : 2025-12-06DOI: 10.1016/j.ijheatfluidflow.2025.110170
Pengtao Wang , Rui Liang , Shimin Wang , Junyu Feng , Jie Zhang , Shengshan Bi
This study utilizes a geothermal heating project located in Xi’an as a case study, developing full-scale numerical models for medium-deep coaxial borehole heat exchangers (DCBHE) and medium-deep U-type borehole heat exchangers (DUBHE) at equivalent depths. Through extensive heat transfer simulations conducted over five cycles, the research provides a comprehensive comparison of the heat extraction efficiency and economic viability of the two types of ground heat exchangers. This comparison employs a ratio method that accounts for pump energy losses, the levelized cost of heat (LCOH), and net present value (NPV) metrics. The findings indicate that the heat exchange per unit length for DCBHE and DUBHE ranges from 123.1 to 149.5 W/m and 93.9 to 170.71 W/m, respectively, with a critical flow rate of 7.21 kg/s. When incorporating pump power consumption, the critical flow rate is determined to be 6.62 kg/s, which corresponds to a comprehensive heat exchange per unit length of 127.7 W/m. Beyond this flow rate, the comprehensive heat transfer performance of DUBHE surpasses that of DCBHE. Therefore, this study divides the heat extraction zones into low-flow zones and high-flow zones based on the critical flow rate of 6.62 kg/s, providing a qualitative identification basis for selecting appropriate flow rate conditions for different types of heat exchangers. The economic analysis reveals that at the critical flow rate, the cost per unit of heat extraction for DCBHE exceeds that of DUBHE, and as the heat extraction duration increases, the NPV of DUBHE becomes more favorable compared to DCBHE. Consequently, DUBHE is determined to be more economically viable than DCBHE.
{"title":"Comparative study on heat extraction performance and economic feasibility between medium-deep coaxial and U-type ground heat exchangers","authors":"Pengtao Wang , Rui Liang , Shimin Wang , Junyu Feng , Jie Zhang , Shengshan Bi","doi":"10.1016/j.ijheatfluidflow.2025.110170","DOIUrl":"10.1016/j.ijheatfluidflow.2025.110170","url":null,"abstract":"<div><div>This study utilizes a geothermal heating project located in Xi’an as a case study, developing full-scale numerical models for medium-deep coaxial borehole heat exchangers (DCBHE) and medium-deep U-type borehole heat exchangers (DUBHE) at equivalent depths. Through extensive heat transfer simulations conducted over five cycles, the research provides a comprehensive comparison of the heat extraction efficiency and economic viability of the two types of ground heat exchangers. This comparison employs a ratio method that accounts for pump energy losses, the levelized cost of heat (LCOH), and net present value (NPV) metrics. The findings indicate that the heat exchange per unit length for DCBHE and DUBHE ranges from 123.1 to 149.5 W/m and 93.9 to 170.71 W/m, respectively, with a critical flow rate of 7.21 kg/s. When incorporating pump power consumption, the critical flow rate is determined to be 6.62 kg/s, which corresponds to a comprehensive heat exchange per unit length of 127.7 W/m. Beyond this flow rate, the comprehensive heat transfer performance of DUBHE surpasses that of DCBHE. Therefore, this study divides the heat extraction zones into low-flow zones and high-flow zones based on the critical flow rate of 6.62 kg/s, providing a qualitative identification basis for selecting appropriate flow rate conditions for different types of heat exchangers. The economic analysis reveals that at the critical flow rate, the cost per unit of heat extraction for DCBHE exceeds that of DUBHE, and as the heat extraction duration increases, the NPV of DUBHE becomes more favorable compared to DCBHE. Consequently, DUBHE is determined to be more economically viable than DCBHE.</div></div>","PeriodicalId":335,"journal":{"name":"International Journal of Heat and Fluid Flow","volume":"118 ","pages":"Article 110170"},"PeriodicalIF":2.6,"publicationDate":"2025-12-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145748813","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}