Search

• A comprehensive review on numerical simulations of boiling heat transfer from atomistic to mesoscopic and macroscopic scales is presented. • Critical issues related to nanoscale bubble nucleation mechanisms, pool boiling, and flow boiling are highlighted. • Potential solutions and future research in the field of boiling heat transfer using the MDS, LB and CFD methods are proposed. Boiling is an efficient heat transfer mode with significant potential for thermal management in high-power electronic equipment. However, a comprehensive understanding of the boiling process, which encompasses bubble nucleation, growth, coalescence, slipping, and detachment across various scales, remains challenging. Molecular dynamics simulation, lattice Boltzmann, and computational fluid dynamics methods are popular and powerful tools for investigating boiling heat transfer phenomena at microscopic, mesoscopic, and macroscopic scales. These methods enable researchers to uncover the underlying boiling mechanisms and propose heat transfer enhancement techniques. Therefore, this paper provides a comprehensive review of boiling heat transfer, spanning from atomistic to mesoscopic and macroscopic scales, utilizing these three numerical methods. It addresses critical issues related to nanoscale bubble nucleation mechanisms, pool boiling, and flow boiling, and proposes potential solutions and future researches, supplementing our previous review [Some advances in numerical simulations of multiscale heat transfer problems and particularly for boiling heat transfer, Annu. Rev. Heat Transf., 6 (2022) 217–269]. Besides, by shedding light on the characteristics of these numerical methods in studying boiling heat transfer, this paper aims to foster their development and advance enhanced heat transfer technologies. [ABSTRACT FROM AUTHOR]

• Local condensation heat transfer in smooth channels and pin fin array is presented. • The effect of operating parameters on the heat transfer characteristics was analyzed. • Sensitivity analysis is performed of heat transfer to identify the critical parameters. • Condensation heat transfer increases with heat flux in the range of vapor quality above 0.4. In this paper, an experimental study of flow condensation was carried out of R134a in a multiport mini channel and a mini channel with pin fin array. The former comprised 20 parallel rectangular channels with an equivalent diameter of 0.64 mm. The latter is a narrow rectangular mini channel containing 10 rows of diamond pin fins with a staggered configuration, and height and longitudinal/transverse pitch of 0.5 mm and 2 mm respectively. To acquire the local value of heat flux and heat transfer coefficient, air jet impingement cooling devices were adopted to condense the vapor of refrigerants. The effects of vapor quality, mass flux, heat flux, and saturation pressure on flow condensation heat transfer coefficient were investigated with the operating conditions: vapor quality from 1 to 0, mass flux from 160 to 450 kg/(m2s), heat flux from 10.0 to 39.8 kW/m2, and saturation pressure from 600 to 1500 kPa. The experimental results indicated that the pin fin array significantly improves the condensation heat transfer coefficient up to nearly 186 % compared to the multiport mini channel in the high vapor quality region. For both channels, the local heat transfer coefficient increases with an increase in vapor quality, mass flux, and heat flux whereas decreases with increase in saturation pressure. The influence of heat flux and mass flux on heat transfer coefficient was more pronounced in the high vapor quality region than in the low vapor quality region. A sensitivity analysis was performed at different vapor quality levels. The most influential parameters in the smooth channel and the pin fin array are saturation pressure and mass flux, respectively. The relative contribution of heat flux is more obvious in the high vapor quality region. Some of the existing correlations have good prediction performance for the smooth channel experimental data in this paper, and the corresponding MAD is within 20 %. However, their deviations are much greater for the applications of pin fin array. [ABSTRACT FROM AUTHOR]

• The liquid metal cross flow tube bundle under pulsation condition is simulated. • The effect of different pulsation frequency and amplitude are analyzed. • Circumferential flow field and temperature field of the single tube are analyzed. • The vortex characteristics and turbulent structure are analyzed with PSD and CWT. Liquid metal has a broad application prospect in tube bundle heat exchangers because of its excellent thermal conductivity. However, conventional in-line tube bundle arrangements often exhibit limited secondary flow, so new methods are urgently needed to improve heat transfer. The topic of this paper is to combine two methods, liquid metal cross flow tube bundle and pulsating flow, to improve the heat transfer performance in-line tube bundles. Firstly, the applicability of the model was validated through experimental data from existing literature, encompassing overall flow dynamics, heat transfer performance, and circumferential pressure distribution around the tubes. Subsequently, numerical simulations are conducted involving liquid metal cross flow in-line tube bundles at varying pulsating frequencies, amplitudes, and Reynolds numbers. Circumferential pressure distributions and temperature profiles for individual tubes were analyzed and compared. The amplitude of reduction in thermal resistance compared to the no-oscillation condition increases gradually from a minimum of 8 % in the first row to a maximum of 40 % in the third row, while the reduction in the rear row of tubes is similar to that of the third row. Remarkably, a reversal in the trend of thermal resistance for the first three rows of tubes under pulsating flow were observed, contrasting conventional flow except for the case with an amplitude of 0.1. By analyzing the local transient flow field distribution near single tube and the corresponding circumferential heat transfer performance, the influence of the vortex evolution characteristics on the heat transfer performance at different phases of the pulsating velocity is revealed, which transient Nu is up to three times that of flow with no pulsating. Finally, a series of coherence analyses reveal that as frequency rises, turbulence exhibits a richer spectrum of small-scale vortex structures, intricately linked to enhanced energy cascade efficiency. Amplitude augmentation intensifies velocity fluctuations, particularly at high amplitudes, resulting in substantial energy peak amplification. The integrated heat transfer performance PEC factor varies in the range of 1.25 to 1.51 for all the conditions simulated in this paper, which indicates that the liquid metal coupled pulsating flow approach can significantly increase the comprehensive performance of the tube bundle heat exchangers. These findings hold significant promise for advancing heat exchanger design and enhancing thermal efficiency, with implications for various engineering applications. [ABSTRACT FROM AUTHOR]

• HFE-7100 boiling curves from the onset of boiling to dry-out, in a vertical and rectangular mini-channel. • Local heat transfer coefficient (HTC) calculated with an inverse heat conduction method. • Influence of subcooling and pressure on heat transfer, pressure drop and flow patterns. • Increasing pressure improves the critical heat flux (CHF) and HTC, delays the dry-out, and reduces pressure drops. • Increasing subcooling improves CHF and HTC, delays the dry-out, and modifies the flow pattern. The background of this study is the problem of thermal runaway in the battery cells of electric vehicles. Increasingly powerful and compact battery packs require highly efficient thermal management. One potential solution is the direct liquid cooling of lithium-ion battery cells using a dielectric liquid which boils over a safety limit temperature. This approach intends to prevent in a first step the cell thermal runaway and if it fails in a second step its propagation to its neighbors: the liquid boiling helps to mitigate the cell temperature rise and to transport and spread the heat. The arrangement of the cells in the pack and their close spacing induce a scientific problem of convective boiling in a mini-channel. To enhance designs using this solution, a comprehensive understanding of flow boiling and a precise local heat transfer quantification is essential. The HFE-7100 has been selected as a working fluid for its attractive thermophysical properties and a suitable saturation temperature. The present paper introduces new experimental results of flow boiling dielectric fluid, the HFE-7100, in a vertical rectangular mini-channel. Measurements are carried out from start of boiling to dry-out. Boiling curves are obtained for a mass flux of 391 kg/(m².s), at different pressures (0.7, 1, 1.5 bar) and at different subcoolings (Δ T s u b = 0 ; 15 ; 30 ; 45 ∘ C). The test section is an aluminum block instrumented with three lines of five K-type thermocouples located all along the flow. The visualization of the flow is achieved through a glass pane with a high-speed camera to identify different flow patterns. The local heat transfer coefficient is calculated by a resolution of a 2D inverse heat conduction method. The experiment mainly aims at determining the heat transfer coefficient, the critical heat flux (CHF) and the pressure drop and at characterizing both the flow regimes and the dry-out phenomenon. The influence of subcooling and pressure on flow boiling is analyzed. An increase in subcooling enhances the CHF and the heat transfer coefficient, and delays the dry-out. Subcooling has an impact on the flow pattern and modifies the type of CHF. In the operating conditions tested, a decrease in pressure enhances CHF and delays dry-out. On the other hand, the heat transfer coefficient increases with rising pressure. Similarly, pressure drops are lower with higher pressure. [ABSTRACT FROM AUTHOR]

• The mathematical model can accurately predict the melting heat transfer process of PCMs. • The Rayleigh number increased with the heat flux, and the higher the heat flux, the larger the amplitude. • Increasing the heat flux improved the integrated heat transfer coefficient and thermal storage performance of the CPCMs. • Conduction was the primary heat transfer mechanism; however, the proportion of natural convection increased with the heat flux. To investigate the effect of heat flux on the enhanced heat transfer in phase change materials (PCMs) composed of copper metal foam (CMF) and paraffin wax (PX), a semi-cylindrical visualized heat storage device with a varied heat flux of 7.3 kW/m2, 8.5 kW/m2, 9.7 kW/m2, and 10.9 kW/m2 was established in this paper. Moreover, a 3-D mathematical model of composite phase change materials (CPCMs) was established using hybrid grid to study the temperature distribution, imbalance effect, flow, heat transfer mechanisms, and heat storage performance during the melting process. The results showed that the melting time of the CPCMs decreased as the heat flux increased. Compared with a heat flux of 7.3 kW/m2, the melting time was shortened by 10.05 %, 20.71 %, and 27.21 %, respectively. In addition, the Rayleigh number (Ra) increased with the increase in heat flux, and the amplitudes of Ra were 1.45×108,1.54×108,1.61×108, and 1.74×108, respectively, which indicated that the higher the heat flux, the larger the amplitudes. Moreover, the heat transfer mechanism during the melting process was dominated by conduction; however, with the increase in heat flux, the maximum flow velocity of the liquid paraffin increased, which caused the proportion of natural convection to increase from 27.80 % to 31.30 %. Hence, the integrated heat transfer coefficient of CPCMs increased from 5.69 W/(m·K) to 5.81 W/(m·K). However, the temperature imbalance effect of the CPCMs was exacerbated, as evidenced by the increased maximum temperature difference in the vertical direction of the center of the CPCMs from 24.7 K to 35.6 K. Besides, the heat storage performance was enhanced. Compared with a heat flux of 7.3 kW/m2, the heat storage capacities of the CPCMs increased by 4.21%, 9.46%, and 12.66 % with increasing heat flux, and the heat storage rates of the CPCMs increased by 15.3 %, 34.9 %, and 52.3 %, respectively. Furthermore, the relative error of the overall melting time was 3.4 %, indicating that the model provided reasonable predictions. [ABSTRACT FROM AUTHOR]

• 1 Flow and boiling heat transfer in a 5 mm tube were studied experimentally. • 2 Heat transfer enhancement is obviously founded in the region where x e < 0.05. • 3 New empirical correlations are developed based on the experimental data. • 4 New correlations have higher prediction accuracy than the existing correlations. Experimental investigation on flow boiling heat transfer in a small tube under high system pressure has been conducted in this study. The test-section is constructed with a 316L stainless steel tube. The inner diameter and heated length of the test-section are 5 mm and 930 mm, respectively. The working fluid is deionized water. The operating conditions of the experiments are as follows, the mass flux from 548 to 978 kg/(m2·s), the wall heat flux from 136 to 238 kW/m2, and the system pressure from 1.0 to 2.9 MPa. The experimental results showed that the inlet subcooling mainly affects the onset of nucleate boiling, but its impact on the heat transfer coefficient at the same vapor quality is almost negligible. The heat transfer coefficient increases with the increase of mass flux and wall heat flux, especially in the region where x e < 0.02. The average heat transfer coefficient is higher under high system pressure owing to the variation of the thermophysical properties. Additionally, the existing flow and heat transfer correlations cannot predict the flow boiling heat transfer coefficient and frictional pressure drop in the experiments accurately enough. Therefore, the new flow and heat transfer correlations are developed based on the experimental data in this paper, and the mean absolute deviation of the flow boiling heat transfer correlation and frictional pressure drop correlation are 13.3% and 12.1%, respectively. [ABSTRACT FROM AUTHOR]

• Bionic method developed to solve the issue of large melting dead zones of TES units. • The lotus root-inspired structure significantly enhanced the heat transfer. • Shorter melting time and higher power density achieved without reducing PCM volume. • A design complying with industrial standards optimized for low-cost mass production. Thermal energy storage (TES) is crucial in the efficient utilization and stable supply of renewable energy. This study aims to enhance the performance of shell-and-tube latent heat thermal energy storage (LHTES) units, particularly addressing the issue of the significant melting dead zones at the bottom, which are responsible for the long charging time. This paper proposed a new heat transfer enhancement technique inspired by the air channel distribution inside the root of the lotus. Numerical simulations are used to explore its melting behavior and heat storage performance, and a comparison is made with conventional shell-and-tube TES units. The results indicate that compared to the single-tube type, the bionic-lotus root type reduced the total melting time by 89.1 %, increased the average temperature by 13.2 °C, and enhanced the average effective power density by 7.6 times. Subsequently, multi-objective optimization was conducted based on response surface methodology, and an industrial standard type was constructed according to existing pipeline products to meet the manufacturing standards for mass production. The melting time was further decreased by 5.8 %, and the average effective power density increased by 6.6 %. The results of this study provide a promising and prospective solution for enhancing shell-and-tube LHTES units, with the potential to increase the efficiency, reduce the footprint and manufacturing costs of energy systems equipped with TES. [Display omitted] [ABSTRACT FROM AUTHOR]

• Algorithm for PLIC-based interfacial area concentration on 3D non-orthogonal grids. • Thermal phase change model with geometrical basis and free of empirical input. • Oscillating behavior of GOX jet plume in LOX crossflow is reproduced in simulation. • Imbalance between the mass transfer rate and the supply causes the oscillation. • Methodology developed is verified by the Stefan problem and the bubble experiment. In order to enhance the physical basis of the mass transfer model in numerical simulation of gas-liquid direct contact condensation and improve the fidelity, a set of methods and models are developed. Based on the concept of the piecewise-linear interface calculation (PLIC), an interface reconstruction algorithm is proposed to calculate the interfacial area concentration in three-dimensional non-orthogonal structured grids. The actual interface in the cell is replaced by a plane, whose position is determined iteratively to satisfy the local volume fraction and orientation of the actual interface. The area of the polygonal cut by the plane is taken as the interfacial area in the cell. The proposed algorithm is validated by calculating the surface area of a series of basic geometries. Obvious superiority of the proposed algorithm over the traditional models (gradient model and algebraic model) is seen. Then, the reconstructed area is used to determine the interfacial area concentration in the thermal phase change model. The process is implemented by ANSYS Fluent user-defined functions. The Stefan problem and steam bubble condensation experiment are used to test the performance of the developed methods and models. The simulation results agree well with the theoretical and experimental results with the quantitative relative errors no more than 4.3 % and 15 %, respectively. The oxygen jets condensation experiments are also used as test cases. The primary dominant frequencies predicted by the simulations are 163.01 Hz and 131.65 Hz, which agree well with the 139.72 Hz and 125.75 Hz in the experiments. The relative errors are only 17 % and 4.7 %. The unstable flow pattern of oxygen jet in liquid oxygen crossflow is explained based on the dynamic characteristics of turbulence, heat and mass transfer. The periodic variation of mass transfer rate is considered as the cause of condensation oscillation. The method for calculating the interfacial area proposed in this paper can serve as an optional and practical improvement to the numerical simulation of direct contact condensation. [ABSTRACT FROM AUTHOR]

• The dynamics of the thermal runaway front in large-scale Li-ion batteries are deduced theoretically. • The velocity of the front yields a 1/2 power relationship between the thermal conductivity and reaction rate. • Freezing thermal runaway front is possible in theory and validated by experiments. • Heat dissipation enlarged by heat transfer time delay is critical for the freeze. Large-scale lithium-ion batteries are favored in electric vehicles and energy storage stations; for instance, BYD blade batteries and CATL Kirin batteries are popular. A tiny defect will trigger thermal runaway from a local point towards the other parts of the battery, and the boundary between the failure region and the intact region is called the "thermal runaway front," of which the dynamics are still unclear. Herein, this paper deduces two formulas for the moving speed (v TR) of the thermal runaway front by applying the principles of energy conservation (Formula I) and flame front (Formula II); formulas for the thermal runaway front velocity (v TR) are derived. The velocity v TR conforms to the 1/2 power of the heat conductivity and reaction rate in both formulas. Data regression for the experimental data demonstrated that the overall estimation error was 10 % using Formula I and 25 % using Formula II. As the zero velocity value indicates the termination of thermal runaway propagation, further derivations reveal that forced cooling can suppress thermal runaway as long as the heat dissipation power is higher than the heat generation power. Sufficient heat dissipation time is achieved by the time delay of heat transfer in the long run of thermal runaway propagation. The formula for the dynamics of thermal runaway propagation universally applies to various large-scale lithium-ion batteries, assisting the safe design of high-energy batteries and providing new insights into suppressing battery thermal runaway by the heat transfer delay effect. [ABSTRACT FROM AUTHOR]

• A Bayesian optimization is conducted utilizing an automated framework to enhance the trade-off between heat transfer and pressure drop using spline-based two-dimensional pin-fin arrays. • The optimization objectives include minimization of pressure drop, maximization of heat transfer, and maximization of efficiency. • The maximization of efficiency optimization yields designs which balanced aerodynamic characteristics with total surface area, leading to enhanced thermal performance. • Pin-fin arrays using four splines outperform all other geometries due to the symmetric nature of the shape generation. Pin-fin arrays are crucial in the cooling of aerospace components. Specifically, pin-fins find extensive usage in the trailing edge of turbine blades, primarily due to the harsh operating conditions these blades experience. With the development of additive manufacturing, complex internal geometries that were previously too expensive to manufacture are now feasible. This advancement has created a new pin-fin design space to be explored. In this paper, Bayesian optimization (BO) is conducted to find the pin-fin array that optimizes the trade-off between heat transfer and pressure drop in a two-dimensional channel. The complexity of the pin-fin geometry varies, and is controlled by the number of connecting piece-wise cubic splines used to close the shape. Different designs are obtained using 3-, 4-, 5-, and 6-splines in a staggered pin-fin array orientation. The optimization objectives included the minimization of pressure drop, the maximization of heat transfer, and the maximization of efficiency. It is found that the pin-fins which utilize 4-splines outperform all others due to their symmetric characteristics. The performances of the pin-fin array plateaus after increasing the complexity beyond 5-splines. It is concluded that the optimization objective significantly influences the final shape characteristics and behavior. As expected, the pressure drop minimization finds aerodynamic shapes while heat transfer maximization increases the total surface area in the channel. The maximization of efficiency results in a balance of aerodynamic shapes and total surface area. [ABSTRACT FROM AUTHOR]

• A Reference-Standard-Volume optimization method is proposed to improve CESPB. • The best aspect ratio of CESPB is 2–4 considering efficiencies and pressure drop. • The volume multiple is given at 1.1 with the highest exergy efficiency of 63.72 %. • The thermal and exergy efficiencies are increased by 2.03 % and 3.65 %, respectively. Cryogenic liquids (e.g., liquid air, liquid hydrogen, liquid carbon dioxide) have gained popularity in electricity storage due to their high energy density, no geographical constraints, and environmental friendliness. The cryogenic energy storage packed bed (CESPB) is widely employed as a cold recovery device to enhance the round-trip efficiency of cryogenic energy storage systems. Nonetheless, the cycle efficiencies of CESPB remain relatively low, with limited research investigating efficient methods for determining the design parameters. To address this, a novel optimization method named as the Reference-Standard-Volume (RSV) approach is introduced to enhance the thermodynamic performance of CESPB. In this paper, the complete optimization process of CESPB is illustrated, presenting the optimization outcomes in the view of heat and mass transfer. The findings reveal that a larger aspect ratio leads to the higher efficiencies but increased pressure drops. An optimal aspect ratio of 3 and a volume multiple of 1.1 is identified, achieving the highest exergy efficiency of 63.72 %. Meanwhile, increasing charging/discharging time enhances heat transfer in CESPB, however, it leads to an increase in external cold dissipation, which gives an optimal charging/discharging time of 3/3 h. By implementing the RSV optimization method, the percentage increase of discharging efficiency, thermal efficiency, and exergy efficiency of CESPB is 2.08 %, 2.03 %, and 3.65 %, respectively. These research outcomes offer valuable reference and guidance for future CESPB design and research. [ABSTRACT FROM AUTHOR]

• The new BTMS has two heat transfer paths, including OHP and liquid cooling plate. • Pulsating flow improves OHP stability through periodic enhanced thermal loads. • Optimal OHP stability lies in 38% or 62% duty cycles, 0.03–0.05 Hz frequencies. • Stably operating OHP enhances the new system's heat transfer performance. After a comprehensive review of oscillating heat pipe (OHP) based on battery thermal management system (BTMS), a novel battery liquid-cooling system with a ⊥-shaped OHP is presented to increase the volume utilization efficiency of battery module and the total amount of heat dissipation. The new system retains the heat transfer path of the OHP while incorporating an additional heat dissipation route where the liquid cooling plate is in direct contact with the cells. In this paper, pulsating flow is used to tackle the problem of the OHP failing to operate stably due to the insufficient heat load achieved in the new system. The parametric effects of pulsating flow frequencies, cooling initiation temperatures, and duty cycles on the operational performance of the OHP and the average surface temperature of the battery module ( T b ) within the new system were investigated. The results indicate that pulsating flow periodically elevates the heat load on the battery surface, enhancing the OHP working properly. The smaller the pulsating flow frequency, the greater the impact on the temperature fluctuations of the new system. The T b is lowest when the pulsating flow frequencies range from 0.03 Hz to 0.05 Hz. Beyond a frequency of 0.06 Hz, the OHP cannot operate stably. Additionally, moderately increasing the cooling initiation temperature is beneficial for enhancing the average temperature oscillation amplitude of the OHP. The cooling initiation temperature between 45 °C and 47 °C results in a better T b . Both pulsating flow duty cycles of 38% and 62% contribute to OHP running effectively. The new system with OHP under optimum pulsating flow condition has the optimal T b , eventually rising by 17.5 °C, which is a 9.8% decrease in temperature rise compared to the system without OHP. And the battery module's surface temperature uniformity has been improved through a reduction of 1.3 °C in the maximum temperature difference ( ΔT max ), amounting to a 17.8% decrease. It has been proved that the improved thermal performance of the new system under pulsating flow is attributed to the stable operation of the OHP. [ABSTRACT FROM AUTHOR]

• A reservoir-assisted loop thermosyphon with a micro-grooved evaporator (MGTPLT) was fabricated. • Micro-grooved structures allowed faster startup and more stable operation of the MGTPLT. • The convection heat transfer ratio was proposed to evaluate the heat transfer enhancement. • The MGTPLT reduced evaporator temperatures and achieved thermal performance enhancement. The two-phase loop thermosyphon (TPLT) offers significant advantages in high-efficiency and long-distance heat transport for various thermal and energy-related applications without external power consumption. In this paper, the startup and heat transfer characteristics of a reservoir-assisted loop thermosyphon with a micro-grooved evaporator (MGTPLT) were experimentally investigated and compared with a loop thermosyphon having a smooth-tube evaporator (STPLT). R245fa was used as a working fluid at volumetric filling ratios of 50–80 %. The results demonstrated that the additional microgrooves allowed a shorter startup time characterized by a lower maximum temperature with respect to the temperature overshoot. Compared to the STPLT, the MGTPLT exhibited lower evaporator temperatures and better temperature uniformities; the micro-grooved structures could alleviate local intermittent dry-out at moderate and high levels of power input for filling ratios of 50 and 60 %. The addition of micro-grooved structures enhanced the heat transfer performance of the TPLT characterized by lower thermal resistances. The convection heat transfer ratio (CHTR) was proposed to evaluate the enhancement of the convection heat transfer within the evaporator region after using microgrooves. The CHTR showed a decreasing tend with increasing the power input. At moderate filling ratios, the value of the CHRT was greater than 1.1, and an approximately 18 % increase in the convection heat transfer was achieved under the combined condition of 60 % filling ratio and 250 W power input. The fundamental mechanism of heat transfer performance enhancement with the aid of micro-grooved structures was analyzed and elucidated. [ABSTRACT FROM AUTHOR]

• Analysis of single and double pass solar air heater-phase change material. • Evaluation based of PCM thickness to length ratio. • Improved outlet temperature during non-sunshine hours. • Improvement in useful energy during discharging period. • System's maximum efficiency of 63.34% achieved during PCM discharging. This paper presents a two-dimensional transient model for a solar air heater with phase change material (SAH-PCM), focusing on the thickness-to-length ratio (t/L) of the PCM container. Verified through experiments, the model considers single (SP) and double pass (DP) flow configurations, assessing liquid fraction, dead length, outlet temperature, energy, exergy, and system efficiency during PCM charging and discharging. Increasing the t/L ratio improves peak liquid fraction, reduces dead length, and delays PCM discharge in both configurations, enhancing system energy and exergy efficiencies. In a SP flow system, optimizing t/L at 0.042 yields a peak liquid fraction of 0.928, a 92.28% dead length reduction, and a 2 h and 22 min delay for PCM discharge (base-case t/L ratio: 0.0182), achieving a highest discharging efficiency of 39.1%. For a DP flow system, the optimized t/L ratio of 0.058 results in a peak liquid fraction of 0.976 and a 1 hour and 34 min delay for PCM discharge, attaining the highest discharging efficiency of 63.34%. Moreover, an improvement in the energy payback time (EPBT) is observed for both SP and DP flow as the t/L is increased. For the DP flow, a minimum EPBT of 0.68 years (energy basis) and 6.55 years (exergy basis) is observed at a t/L ratio of 0.058. In terms of energy cost, SP flow reaches a minimum of Rs. 20.824 Rs/kWh at a t/L ratio of 0.042. For DP flow, the minimum energy cost is 13.033 Rs/kWh at a t/L ratio of 0.058, representing a noTable 37.42% improvement compared to SP flow. These findings underscore the critical role of optimizing system parameters, specifically the t/L ratio, for enhanced energy and exergy performance, and economic viability of SAH-PCM. [ABSTRACT FROM AUTHOR]

• Thermal behavior of rectangular fins located within a square cavity were analyzed. • Inverse CFD method and experimental data were adopted for heat transfer prediction. • Suitable turbulence model based on the RMSE of the flow model were selected. • Effect of fin geometry on the flow and heat transfer within a cavity was evaluated. • Effect of PCM heat sink on the flow and heat transfer within a cavity was studied. This study examines the fluid flow and heat transfer characteristics of rectangular fins positioned within a square cavity through the use of the inverse three-dimensional computational fluid dynamics (CFD) method, experimental analysis, and the constant heat transfer rate assumption. This paper proposes a heat sink that uses a phase change material (PCM) made of paraffin within a small cavity. The inverse CFD method, combined with the least squares method, root mean square error, and excessive experimental data, is adopted to predict the unknown heat transfer rate Q and absorption heat Q ab. One of the objectives of this study is to investigate the effect of the PCM heat sink on fluid flow and heat transfer characteristics within a cavity. A sequence of tests on various flow models indicates that employing the RNG k-ε turbulence model with the standard wall function is the most suitable choice for all scenarios in the three-fin model. Conversely, the zero-equation model proved to be a more adequate fit for the PCM heat sink. Another objective of this study is to study the effect of the height of the PCM heat sink on Q and Q ab. The incorporation of a PCM heat sink results in a further improvement in the efficiency of heat dissipation. Part of the heat created is absorbed by the PCM heat sink, which absorbs 22.5 % of the thermal energy. Under the same volume, the PCM heat sink with a lower height (H p = 0.012 m) absorbs 14 % more heat compared to H p = 0.02 m. Thus, a lower-height PCM heat sink has a better heat dissipation effect. [ABSTRACT FROM AUTHOR]

• Verified rivulet behavior during dropwise condensation on vertical surface • Identified effects of surface length and rivulet flow on heat transfer • Related rivulet-covered area and dropwise condensation heat transfer Due to its superior heat transfer performance, dropwise condensation has been studied for industrial applications. This paper experimentally reports the effects of vertical length on dropwise condnsation in the presence of noncondensable gas. The formation of rivulet flows, which reduce heat transfer due to droplet coalescence, was observed. Compared to the upper side, rivulet flows on the lower side of the condensing surface were faster and more frequent. Lower heat flux was also measured on the lower side, indicating a correlation between heat transfer and rivulet flow parameters during dropwise condensation on the vertical surface. These findings represent the first study of the effects of rivulet flows on dropwise condensation heat transfer, and they can be applied to evaluate condensation heat transfer on long vertical surfaces, such as those found in small modular reactor containment systems as well as large containment building with passive containment cooling systems. [ABSTRACT FROM AUTHOR]

· Different change rate of cross-sections along the flow field is investigated. · Microchannel condenser is used in loop heat pipe. · Effects of operating height situation and microchannel structure are studied. · Smooth operation and anti-gravity performance are effectively improved with the help of tapered channel. The microchannel condenser has the advantages of high condensation efficiency and low flow resistance. The vapor can liquefy and return to the evaporator rapidly, and the heat transfer performance and the smooth operation of the loop heat pipe can be improved. Suitable change rate of the cross-section along the microchannel flow field can adjust the pressure drop and flow rate of working fluid, and improve the anti-gravity performance of the loop heat pipe. In this paper, microchannels with different change rates of cross-section along the flow field, including 0.5 mm straight channel, 1 mm straight channel, divergent channel, and tapered channel, are designed and applied as condensers in loop heat pipes. The entrance size of the gradient microchannel is 0.5 mm or 1 mm, which uniformly expands to 1 mm or shrinks to 0.5 mm along the flow field. The vapor and cooling water form a countercurrent direction of motion in the microchannel, and the heat exchange is carried out through the bottom of the channel. The influence of the microchannel condenser structure and operating height situation between the evaporator and condenser on the heat transfer performance of the loop heat pipe are experimentally investigated. The results show that the loop heat pipe achieves an extremely low thermal resistance of 0.004℃/W with gravity assist. Under the gravity-free operation situation, the maximum evaporator temperature of the loop heat pipe with 0.5 mm straight channel and 1000 W heat load is 128.1 °C, and the vapor pressure drop fluctuated is in the range of -1000 ∼ 2500 Pa. Under the anti-gravity operation, the loop heat pipe with the tapered channel has effective transfer performance at 500 W, and the influence factor is only 28.1 %. [ABSTRACT FROM AUTHOR]

• A new wall boiling heat transfer model is developed to offer a fresh approach for CFD simulations within the Eulerian two-fluid framework. • Each parameter within the new model is calculated based on local physical properties and macroscopic parameters at cell level. • The new model reasonably predicts the wall temperature in flow boiling cases. • The sensitivity of the coefficients introduced in the proposed model are analyzed. Accurate prediction of flow boiling heat transfer is prominently dependent on the modeling of wall heat flux partitioning. In this paper, a new wall boiling heat transfer model was developed for three-dimensional Computational Fluid Dynamics (CFD) code to predict the wall heat flux and wall temperature. The proposed model partitioned the wall heat flux into convective heat flux and nucleate boiling heat flux, which were further modified by two correction factors. The key feature is that the new wall boiling heat transfer model was derived from bubble growth mechanism, incorporating reasonable assumptions, and each parameter within the model was calculated based on local physical properties and macroscopic parameters at the cell level. On this basis, the new wall boiling heat transfer model was coupled into ANSYS-Fluent and validated against various public experiments as well as the KIMOF experiments conducted under different conditions. Simulation results indicated that the proposed model could predict reasonable results for wall temperature and cross-section average void fraction. Finally, a comprehensive investigation was carried out to assess the sensitivity of the computational grids and the coefficients introduced in the new model. [ABSTRACT FROM AUTHOR]

• Presents an analytical solution for thermal spreading in a cylindrical body with spatially varying side wall cooling. • A series solution is derived, with coefficients shown to be the solution of a set of linear algebraic equations. • Impact of key non-dimensional parameters on temperature field and thermal resistance is investigated. • Results are shown to be in excellent agreement with past work for special cases. • Results are of relevance for edge cooling of microprocessor chips as well as batteries. Thermal spreading and constriction resistance play an important role in thermal management of a variety of engineering systems, including semiconductor devices. Moreover, thermal contact resistance at interfaces between materials is also governed by spreading and constriction processes. Most of the past literature on theoretical modeling of thermal spreading from a source into a larger body assumes an adiabatic side wall, whereas several recent thermal management problems involve active cooling at the side wall. This work addresses such scenarios by solving the thermal spreading problem in a finite thickness cylinder with spatially varying convective heat transfer on the side wall. An eigenfunction-based series solution for the temperature field is written, the coefficients of which are shown to be governed by a set of coupled algebraic equations. Results are shown to agree well with past work for special cases of the general analysis presented here. Thermal spreading resistance is computed for a number of practical scenarios, including step function, periodic and sinusoidal convective cooling on the side wall. In each case, it is shown that the temperature distribution is consistent with the nature of the spatially varying boundary condition. A key theoretical insight from this work is that under certain conditions, the classical definition of thermal spreading resistance may become negative due to the presence of two heat sinks in the problem. Therefore, analysis based on total resistance may be more appropriate. This work offers a significant generalization of past papers, making it possible to analyze realistic thermal management problems, in which, side wall cooling plays an important role. [ABSTRACT FROM AUTHOR]

• Heat transfer coefficient of highly subcooled water during flow boiling. • Challenges in minichannels during subcooled flow boiling. • Wall temperature distribution of mini and conventional channels during flow boiling. • Temporal temperature variations in minichannel at near to zero vapor quality. • Wall temperature fluctuations in minichannel for bo > 1.76 × 10−4. Channels of different sizes are employed in diverse engineering applications to enhance the efficiency of heat transfer system. Minichannels are preferred for their compactness, enhanced heat transfer, larger surface area to volume ratio, and high heat transfer coefficient values. There is a dearth of literature addressing flow boiling in minichannels with highly subcooled water at atmospheric system pressure. This research paper analyses the flow boiling phenomenon in a 2 mm diameter minichannel, comparing it with conventional 4 mm and 11.7 mm diameter tubes while maintaining a constant length to diameter ratio of 150. The analysis investigates the wall temperature distribution, local heat transfer coefficient, and temporal heat transfer fluctuations at a mass flux of 150 – 2400 kg/m2s and 109 – 1080 kW/m2 heat flux. The study explores the impact of boiling number on wall temperature fluctuations. Results reveal that mass flux doesn't impact heat transfer in subcooled flow boiling but does in saturated flow boiling in conventional channel. There are no significant temporal fluctuations of wall temperature in a 4 mm tube. However, vapor bubble size affects the circumferential area of minichannels, leading to abrupt flow disturbance. Wall temperature fluctuations occur in the minichannel for the boiling number exceeding 1.76 × 10−4. Temperature fluctuations in minichannel increase with boiling number, but decrease as mass flux is increased. Delving into the depths of minichannel heat transfer, this intriguing study unveils the pivotal role of the boiling number in ensuring process stability. [ABSTRACT FROM AUTHOR]

In this paper, we perform a detailed numerical study of the linear dynamics of Rayleigh–Bénard convection in a fluid layer confined in a Hele–Shaw cell and subjected to a temperature gradient with a stationary and a periodic time-dependent component. Using Floquet theory for time and the spectral method for space, we show that the system exhibits three classes of time-dependent solutions, namely harmonic, sub-harmonic and quasi-periodic instabilities. These compete into the primary bifurcation driven by two different mechanisms depending on the modulation amplitude. Indeed, harmonic convective rolls, driven by the stationary thermal gradient, arise below critical modulation amplitudes beyond which the well-known parametric resonance mechanism is responsible for the development of convective rolls. In addition, it is revealed that the convective flow reversal could identify this exchange of instability mechanisms occurring via either bicritical states with codimension-two bifurcation points between harmonic and subharmonic modes or via displacement of Floquet multipliers in the Floquet-complex plane. This latter scenario sets in between subharmonic and quasiperiodic instability modes in the high frequency limit by varying the phase-lag between the temperatures of the lower and upper boundaries. Moreover, it is shown that the nature of the instability mechanism alters considerably the heat transfer dynamics of the system. • A Floquet-based analysis on time-modulated Rayleigh-Bénard convection in a Hele-Shaw cell is presented. • Three classes of time-dependent instability modes are found namely harmonic, subharmonic and quasiperiodic instabilities. • Flow reversal is used to identify the primary bifurcation mechanism. • The effect of the instability mechanism on the heat transfer of the system is also presented. [ABSTRACT FROM AUTHOR]

• The process of melting-reconsolidation of the frost layer is considered. • Improvement of the existing Eulerian multiphase flow model of frost formation. • The whole process of frost layer growth under prolonged periods of time was studied. • Changes of ice-phase volume fraction and melt-reconsolidation region with time are discussed. Frost growth is a complex heat and mass transfer process, which involves the vapor-liquid-solid three-phase interconversion problem, and the investigation of the frost growth process helps to understand the mechanism of frost growth. In this paper, frost model is considered, taking into account the role of the frost layer melting-reconsolidation. The aim is to improve the existing Eulerian multiphase flow model. Using the improved frost model, the whole process of the frost layer growth is simulated. In addition, an investigation is presented dealing with, the frost layer growth in various stages of the ice phase volume fraction and melting-reconsolidation area, considering the change and the proportion of the situation. The results showed that the volume fraction of the ice phase increased with time. Additionally, the melting-reconsolidation process starts from the front area of the frost layer surface and then gradually spreads to the whole surface. Moreover, it is shown that the proportion of the melting area increases and then decreases with time before attaining a constant value. On the other hand, the proportion of the solidifying area increases with time and finally tends to 0. [ABSTRACT FROM AUTHOR]

This paper analyzes the influence of the anomalous diffusive effects caused by micro-scale heterogeneity and kinetic and inertial thermal relaxations on the optically induced thermoelastic bending component of the photoacoustic response. We calculated the temperature distribution for a one-dimensional heat transfer problem with planar and periodic excitation, neglecting the influence of thermoelastic strains on the temperature profile. Thermoelastic bending was evaluated using a theoretical approximation of a thin plate, while pressure fluctuations in the photoacoustic cell were obtained by assuming adiabatic changes in the closed air. The model analysis shows that the relaxation processes could significantly affect the mechanical piston component of the photoacoustic response at frequencies higher than the minima of the inverse of two thermal relaxation times, while the influence of micro-scale heterogeneity is observable in the whole frequency range. • The hyperbolic DPL model enhances accuracy of thermoelastic effect at higher frequencies. • Sensitivity to fractional order of DPL model emphasizes precision in TE control. • Anomalous DPL heat diffusion shifts and attenuates the peak of the hyperbolic result. • Maximum damping occurs to identical kinetic and thermal relaxation times. [ABSTRACT FROM AUTHOR]

• An equivalent interface thin layer temperature field model (ETTM) is proposed considering the effect of TCR. • ETTM can describe the temperature gradient caused by rough contact interfaces. • The validity of the TCC model and ETTM is verified by comparing the model proposed with the previous research findings and ANSYS analysis results. • The effects of multiple factors on thermal contact conductivity, equivalent thin layer parameters and temperature field are analyzed. Improving the accuracy and efficiency of predicting the temperature field is a crucial requirement for iteratively calculating and online monitoring the temperature field, particularly for combined structures. The heat flow movement is impeded by incomplete interface contact, resulting in a temperature jump. To address this, the equivalent thin-layer temperature field model (ETTM) is proposed to replace the interface thermal contact conductance (TCC) model. Firstly, the fractal TCC model of a bolt connection is derived based on fractal geometry. Subsequently, the contact surface of the circumferential bolt connection is considered to be the parallel connection of multiple bolt thermal resistances. The TCC of the equivalent thin-layer is obtained. Building on this, the ETTM of the built-up rotor structure is developed. Finally, a temperature field model of the bolted rotor is derived. The validity of models is verified by comparing the model proposed in this paper with the previous research findings and ANSYS analysis results. Numerical simulation results reveal that TCC is influenced by several factors: contact pressure, the nominal contact area within the fractal domain, fractal dimension, fractal roughness, and the root mean square (RMS) height. Furthermore, the effects of the thickness of the thin layer model and the thermal conductivity of the material on heat transfer are examined. Ultimately, the findings indicate that heat transfer at the interface is significantly impeded by thermal contact resistance across various temperature boundary conditions. Notably, lower contact pressure intensifies this impedance. [ABSTRACT FROM AUTHOR]

• Optimize the thermosyphon loop thermal performance by controlling loop wettability. • Long-term effective hydrophilic and hydrophobic modification method was determined. • The optimum wall wettability under different filling ratios was clarified. • The optimized R loop is reduced by 6.9% to 53.5% under different filling ratios. The behavior of vapor-liquid two-phase transport within the riser of the two-phase thermosyphon loop significantly impacts the heat transfer process from the evaporation to the condensation, as well as the circulation velocity, thus playing a crucial role in the heat transfer performance. However, existing research has predominantly focused on optimizing the performance of the evaporator, overlooking the optimization of vapor-liquid two-phase flow patterns and transport characteristics within the loop. This paper aims to optimize the heat transfer performance of the two-phase thermosyphon loop by controlling the wall wettability of the loop. The strategy employed involves modifying the wettability of the riser to regulate the two-phase flow within the loop, while visualization techniques are utilized to investigate the influence of riser wettability on the vapor-liquid two-phase distribution and heat transfer performance of the two-phase thermosyphon loop. The findings reveal that the impact of wall wettability on flow behavior and heat transfer performance within the loop varies significantly with changes in the filling ratio. This paper sheds some light on the riser wettability of two-phase thermosyphon loops to optimize the thermal performance. [ABSTRACT FROM AUTHOR]

• A compact jet impingement boiling heat transfer device achieves a high-powder cooling power of 1532 w. • The variation trends of critical heat flux, heat transfer coefficient and pressure drop with mass flux, saturated pressure and liquid subcooling are obtained. • Compared with smooth surface, a micro-pin-finned surface can increase CHF and HTC by up to 60 % and 32 %, while only increasing a pressure drop of 14 %. In this paper, jet impingement boiling heat transfer performance was studied using a novel new environmentally friendly refrigerant medium (HP-1) based on a high-power chip heat sink. A series of experimental investigations of jet impingement boiling were conducted on smooth and micro-pin-finned surfaces (PF0.5–0.5–1, PF0.5–0.5–2 and PF0.3–0.3–1). Heat transfer performances were evaluated over the saturated pressures from 450 to 590 kPa, liquid subcoolings from 2 to 15 K, nozzle sizes from 1 to 2 mm, mass fluxes from 71 to 858 kg/(m2s). The results show that critical heat flux (CHF) and heat transfer coefficient (HTC) increases with the increase of mass flux and saturated pressure. At ΔT sub =15 K, the CHF of the micro-pin-finned surfaces can be increased by 63 % compared with smooth surface, in which the CHF of PF0.5–0.5–2 surface is up to 1532 kW/m2. Boiling heat transfer coefficient first increases and then decreases with the increase of heat flux. Due to the increase of the inlet liquid subcooling, leading to a large temperature difference, so the heat transfer coefficient is smaller at a high liquid subcooling. Several factors, such as mass flux, saturated pressure and liquid subcooling, have a great impact on pressure drop. With the increase of vapor quality, the shear action between the vapor-liquid phase is enhanced, and the pressure drop will gradually increase. Bubble behaviors have also been captured by a high-speed camera to help analyze the mechanism of jet impingement boiling heat transfer in a confined space for high-power chips. [ABSTRACT FROM AUTHOR]

• Physics-Informed neural networks (PINNs) are used for single and two-phase flows. • The study investigates the rising bubble problem with and without heat transfer. • The pertinent physics of wake interaction and thermal boundary layer are discussed. • Physics-Informed neural networks are agnostic to geometry and fluid properties. • The maximum error in the position of the centre of mass by PINN is less than 2.8 %. This paper presents data-driven simulations of two-phase fluid processes with heat transfer. A Physics-Informed Neural Network (PINN) was applied to capture the behaviour of phase interfaces in two-phase flows and model the hydrodynamics and heat transfer of flow configurations representative of established numerical test cases. The developed PINN approach was trained on simulation data derived from physically based Computational Fluid Dynamics (CFD) simulations with interface capturing. The present study considers fundamental problems, including tracking the rise of a single gas bubble in a denser fluid and exploring the heat transfer in the wake of a bubble rising close to a heated wall. Tracking of a rising bubble phase interface of fluids with disparate properties was performed, revealing a maximum error of only 5.2% at the interface edge and a maximum error of 2.8% at the position of the centre of mass. Inferred (hidden variable) flows are studied in addition to a purely extrapolative inverse isothermal bubble case. When no velocity data was supplied, velocity field predictions remained accurate. Rise of an inferred isothermal bubble with unseen fluid properties was found to produce a maximum mean-squared error of 0.28 and centre of mass error of 1.25%. For the case of the rising bubble with a hot wall, the maximum error in the temperature domain using specified boundary conditions was 6.8%, while the bubble position analysis reveals a maximum positional error of 3.6%. These results demonstrate that PINN is agnostic to geometry and fluid properties when studying the combined effects of convection and buoyancy on two-phase flows for the first time. This work serves as a starting point for PINN in multiphase cases involving heat transfer over a range of geometries. Eventually, PINN will be used in such cases to provide solutions for forward, inverse, and extrapolative cases. Each of which represent a dramatic saving in computational cost compared to traditional CFD. Contours of volume fraction (α) and temperature (T*) predicted by the CFD (Top) and PINN (Bottom) for the evolution of flow showing how the rising bubble influences the thermal boundary layer on the wall. [Display omitted] [ABSTRACT FROM AUTHOR]

• Heat and mass transfer model of pyrolysis gases in a confined space was developed. • Growth rates of pyrolysis product concentrations were determined. • Synergistic effects were identified for different fire hazard sources. • Role of ventilation in the growth of pyrolysis gases concentrations was defined. • Correlation of room sizes with concentrations of pyrolysis gases was established. Numerical research findings are presented on the patterns of spread of thermal decomposition products of wood-based materials in confined spaces isolated from the external environment. Based on the experimental data, a physical problem statement was formulated and a mathematical model of heat and mass transfer was developed for rooms with the most typical sources of fire hazard. The experimental and numerical research findings were obtained for the conditions corresponding to compartments isolated from the external environment (e.g., rooms at dedicated facilities, like ships, small and medium-sized rooms) and to underventilated rooms. The performed experiments and calculations revealed satisfactory agreement in the findings and demonstrated the validity of the proposed model of heat and mass transfer. Wood was used as the main material in fires. A comparative analysis was also carried out for paper of different quality. The modeling was performed for the most frequent causes of fires: careless handling of fire, unsafe operation of heating equipment and electrical short circuits. The rates of growth of thermal decomposition product concentrations were determined for different distances from the material. The impact of a number of factors on the characteristics under study was determined. These are the fire and room sizes, fire cause, supply and exhaust ventilation. The data were extrapolated to the conditions corresponding to business, warehouse and residential premises. [Display omitted] [ABSTRACT FROM AUTHOR]

• Performance evaluation of three different cooling methods for data centers is conducted. • Single-phase immersion cooling has advantages in heat transfer and energy consumption. • Correlation between oil properties and immersion cooling performance is established. • Dynamic viscosity has the greatest impact on the flow and thermal performance of oil coolants. This paper conducts a comprehensive performance evaluation and optimization of single-phase immersion cooling (SPIC) systems to address challenges in data center server cooling. A three-dimensional model of blade servers with different cooling methods is constructed and numerically solved to compare thermal transport performance, flow characteristics, and energy consumption performance. Furthermore, the influence mechanisms of different oil coolants and operating conditions on SPIC system performance are studied. The results indicate that the SPIC scheme exhibits significant advantages in terms of heat transfer capability, heat transfer uniformity, and flow resistance loss. Moreover, the power usage effectiveness (PUE) of SPIC systems decreases by 20.8 % and 17.6 % compared to forced air cooling and water cooling plate solutions, respectively. The dynamic viscosity of oil coolants has the greatest impact on the SPIC performance. Moreover, high-viscosity silicone oil increases the flow resistance of SPIC systems and weakens the convective heat transfer capability between the coolant and components. By ensuring safe operation, increasing the inlet temperature of oil coolants reduces flow resistance losses and enhances heat transfer capacity. Additionally, initially increasing the mass flow rate improves convective heat transfer capacity. However, the heat transfer improvement gradually diminished beyond a certain threshold while flow resistance losses increased significantly. [ABSTRACT FROM AUTHOR]

• A review of experimental and CFD methods used to analyze the thermal-hydraulic characteristics of sub-channels in rod bundle geometries is presented. • Both adiabatic and boiling two-phase flow studies are examined. • The effects of the design and flow parameters on the two-phase flow hydrodynamics and heat transfer performance are summarized. • Based on previous numerical studies, a basic CFD model for two-phase flows in rod bundle geometries is proposed. Thermal-hydraulic analysis of two-phase flow plays a crucial role in the optimal design of nuclear fuel assemblies and in nuclear safety. This paper presents the state of the art in two-phase flow hydrodynamics and heat transfer in nuclear fuel assemblies. A comprehensive review of experimental and CFD methods used in the last 30 years to analyze the thermal-hydraulic characteristics of sub-channels in rod bundle geometries under different flow conditions (adiabatic and boiling two-phase flow) is presented. The effects of design parameters, e.g. rod-pitch length, mixing vane angle, and flow parameters, e.g. pressure, mass flux and heat flux, on flow hydrodynamics and heat transfer performance are summarized. In addition, the capabilities of CFD modelling of two-phase in rod bundles are discussed. Moreover, based on the existing studies, some recommendations for future studies are given. [ABSTRACT FROM AUTHOR]

• Numerical study on falling film flow and heat transfer on porous tubes. • Compared falling film flow and heat transfer on plain and porous tubes. • Studied porosity, layer thickness, and material effects on flow and heat transfer. • Porous layer leads lower flow, thinner film, fewer rupture, better heat transfer. • Smaller porosity, thicker layer, aluminum and copper facilitate heat transfer. Porous surfaces often exhibit exceptional performance in liquid spreading and heat transfer, making them promising candidates for enhancing falling film heat transfer on solid surfaces. This paper delves into the flow and heat transfer characteristics occurring within falling films on horizontal tubes coated with a porous layer. To achieve this, a three-dimensional numerical model that couples the Volume of Fluid (VOF) method with porous media considerations was employed. The investigation comprehensively covers film spreading behavior, flow patterns, film velocity, film thickness, temperature distribution, and heat transfer coefficients, for both plain and porous tubes. Furthermore, we delve into the effects of key parameters, such as porosity, layer thickness, and material properties, on the overall heat transfer process. The results demonstrate the accuracy of our model in reproducing film thickness and heat transfer coefficients within falling films on horizontal porous tubes. When comparing plain tubes to porous ones, some key differences emerge. Specifically, the liquid film on porous tubes spreads more slowly, leading to increased film thickness and fewer occurrences of dry patches. Additionally, porous tubes result in lower film velocity and a more uniform temperature field compared to plain tubes. Interestingly, the findings reveal that smaller porosities yield thicker films and higher heat transfer coefficients, particularly when the porosity is below 0.4. Moreover, thicker porous layers result in thinner films, higher wetting ratios, and increased heat transfer coefficients. Among the materials we tested, aluminum and copper porous layers exhibit the highest heat transfer performance, while steel demonstrates the lowest. [ABSTRACT FROM AUTHOR]

• A novel laboratory-scale experimental setup to measure changes in moisture content of textiles during drum drying. • A simple mathematical heat transfer model to assess moisture changes in drum drying. • Quantification of the impact of interfacial thermal contact resistance on drying time. Drum dryers are widely employed in various industrial fields for their remarkable capability to efficiently dry thin porous materials such as textiles and paper. The heat transfer process in drum drying involves an interaction between the drum surface and the materials being dried. However, measuring interfacial thermal contact conductance during drum drying is very challenging, and little data has been reported in the literature. It is therefore unclear how changes in interfacial thermal contact conductance quantitatively affect drum drying in terms of drying time. We created an innovative laboratory-scale experimental setup designed to measure changes in moisture content of textiles at various contact pressures during drum drying. By combining our experiments with a simple yet effective heat transfer model, we were able to quantify the impact of interfacial thermal contact resistance on drying time. Specifically, in drying tests using thin cotton textiles, the experimental results revealed that when the contact pressure was increased from 120 Pa to 480 Pa, the thermal contact conductance increased by approximately 30% and the drying time decreased by approximately 20%. The outcomes of our study significantly advance the current understanding of drum drying processes and reveal vital insights into techniques for optimizing energy consumption. [ABSTRACT FROM AUTHOR]

• The 3D hybrid surface with a tall HB pillar array surrounded by a short HP pillar array is proposed to enhance boiling heat transfer. • The 3D hybrid surface could enhance CHF and HTC simultaneously. • Compared with the plain SiO 2 surface, the 3D hybrid surface gave a 107% higher CHF and a 360% higher HTC. • The height difference between the HB and HP pillar arrays could mitigate liquid–vapor interference and increase the CHF. • The hydrodynamic instability limit corresponded closely with the experimental CHF values of the surfaces. Boiling heat transfer is widely used in numerous industrial systems because of its excellent heat dissipation performance. The critical heat flux (CHF) represents the upper limit of efficient heat removal. The most widely applied CHF model is Zuber's hydrodynamic instability model, which posits that the interference between vapor and liquid counterflows results in the CHF. In this paper, a three-dimensional (3D) hybrid micropillar array surface for hindering the liquid–vapor interface is proposed to enhance the CHF and heat transfer coefficient (HTC). The 3D hybrid surface comprises a tall hydrophobic (HB) pillar array surrounded by a short hydrophilic (HP) pillar array. The tall HB pillar array acts as a potential nucleation site for increasing the HTC, and the short HP pillar array acts as a wicking structure for facilitating liquid rewetting. In addition, the height difference between the HB and HP pillar arrays can potentially reduce the interference between liquid and vapor counterflows and increase the CHF. Relative to a plain silica surface, the 3D hybrid surface increases the CHF and HTC by 107% and 360%, respectively. Theoretical modeling suggests that the CHF on the 3D hybrid surface was still limited by pool hydrodynamics. [Display omitted] [ABSTRACT FROM AUTHOR]

The influence of thermal interface resistance is an important design problem for the optimisation of heat transfer systems. For thermoelectric generators (TEG), it can limit the ability to convert thermal energy directly to electrical energy. This paper describes an experimental and numerical study of the influence of thermal interface resistance on TEG module performance. A test apparatus was constructed to evaluate TEGs under varying temperature and heat source/sink conditions. The effect of heat source/sink interface finish grades, corresponding to medium quality (N8-N6) and very fine (N5-N4) surface finishes, at varying clamp forces, is evaluated. Then the application of cold spray additive manufactured copper, as a novel method of TEG module thermal interfacing is explored. Experimental data is accompanied with numerical analysis using TEG ideal equations and the effective material properties approach with the incorporation of temperature dependency and fabrication quality factors to improve modelling accuracy. A system of simultaneous heat flow equations is then established to numerically solve for the presence of thermal interface resistance. The results show that high clamp forces (3600 N) and very fine finish grades were sufficient to minimise thermal interface resistance to the range 0–5.98 × 10 − 5 m 2 K W − 1 , representing a loss of 0–9.9% in maximal achievable power. For the cold spray thermally interfaced TEG module, resistance estimates were lower and in the range 0–5.4 × 10 − 6 m 2 K W − 1 , corresponding to a 0–1.0% loss in power. • Heat source/sink surfaces, with varying finishes (N8-N4) and clamp forces, are evaluated for influence on thermoelectric generator (TEGs) performance. • Cold spray additive manufacturing is proposed and experimentally evaluated as a novel method of TEG heat source/sink thermal interfacing. • A system of heat flow equations are used to solve for the presence of thermal interface resistance (R T I R) at TEG heat source/sink interfaces. • High clamp forces (3600 N) and very fine finishes (N4) minimised R T I R to the range 0–5.98 × 10 − 5 m 2 K W − 1 , representing a loss of 0–9.9% in maximal achievable TEG power. • R T I R estimates for the cold spray thermally interfaced TEGs were lower and in the range 0–5.4 × 10 − 6 m 2 K W − 1 , corresponding to a 0–1% loss in power. [ABSTRACT FROM AUTHOR]

Impinging droplets on a hot surface is a significant phenomenon in many industrial applications, such as spray cooling. Previous studies reported various morphological behaviours affecting the droplet's wetting, heat transfer, and evaporation characteristics. However, there is an ambiguity in understanding the influence of the surrounding medium on the resulting boiling regimes. The present work aims to study this aspect through an experimental investigation using high-speed photography and infrared thermography. During the experiments, the degassed FC-72 liquid droplets have impinged onto a heated chromium-plated calcium fluoride (C a F 2) surface with a constant impact velocity. Three ambient conditions, such as saturated vapour (heated vapour), air at liquid saturation temperature (heated air), and air at room temperature (ambient air), are used. The surface temperature is varied in each ambience until the complete droplet dewetting, i.e., the Leidenfrost phenomenon has appeared. In this study, the stages of boiling are perceived as film evaporation, excessive and nominal bubbly boiling, transition, and Leidenfrost. With different surrounding media, the onset of a regime and its boiling morphology are distinct. During the film evaporation regime, the droplet heat transfer is high for the impact in ambient air, and the contact line evaporation is dominant in the vapour medium. In contrast, heat transfer is improved in a vapour medium during the bubbly boiling process, reaching up to 30% enhancement compared to the droplets in air surroundings. It is attributed to increased nucleation sites and excessive bubbling compared to air media. Correlations from the literature are also revisited to understand the onset of bubbly boiling and the droplet breakup at the Leidenfrost point. An overall comparison of droplet heat transfer during the boiling regimes in surrounding media is discussed in this paper. • The effect of surrounding media on boiling regimes is experimentally investigated. • Boiling regimes are categorised into film evaporation, bubbly boiling, transition, and Leidenfrost point. • In each regime, the droplet's spreading and heat transfer are analysed using shadow photography and infrared thermometry. • The significance of contact line evaporation in various surroundings during the film evaporation regime is presented. • Correlations related to the onset of bubbly boiling and droplet breakup at Leidenfrost point are revisited. [ABSTRACT FROM AUTHOR]

• Application of MCrAlY coating to increase heat transfer during water spray cooling. • The use of the coating intensifies heat transfer compared to the uncoated surface. • The thickness of the coating affects the intensity of the heat removal. The main purpose of the paper was to determine the influence of NiCoCrAlY coating application and its thickness on the heat transfer during water spray cooling under surface boiling conditions. The scope of the research included: sensor preparation, measurements of the sensor temperature during cooling, and numerical calculations to determine the intensity of sensor cooling. The intensity of heat removal was determined by the average values of the heat flux on the cooled surfaces. The inverse method for the heat conduction equation was used to determine the heat transfer boundary conditions. Temperature measurements during cooling were carried out for an initial sensor temperature of 500°C, 700°C and 900°C and a liquid pressure 0.1 MPa, 0.2 MPa. Due to this, it was possible to select a coating thickness that contributed to the greatest extent to the increase in the intensity of heat dissipation from the hot surface. The highest values of the heat flux were obtained during the cooling of the sensor with a 134.5 μm coating. [ABSTRACT FROM AUTHOR]

• Obtained local surface temperature (T w), boiling curves and liquid loss rate (m loss). • Heat transfer rate (Q w) and m loss increase markedly at onset of nucleate boiling. • Jet outlet velocity (v) affects Q w and m loss more obviously than impact height (H). • Effects of v and H on T w are more significant at non-phase transition stage. • Growing of Q w and M loss slows down at anaphase of nucleate boiling. Combining the benefits of jet impact and phase-transition heat transfer, jet impact boiling has been widely adopted in industrial scenarios. Notably, owing to limitations of high-velocity jets, such as high impact stress and pumping pressure in liquid circuits, low-velocity jet impacting deserves more attention and explorations. In this paper, under an atmospheric environment with wall heating temperatures ranging from 36 °C to 128 °C, the two-phase flow and heat transfer characteristics (local surface temperature, boiling curves and liquid loss rate) of low-velocity jet impacting on a cylindrical surface are investigated by experiment, followed by the analysis on the effects of jet outlet velocity and impact height. The findings reveal that, with the increase of heat transfer rate, the heat transfer regimes in sequence are non-phase transition (single-phase convection and air-liquid two-phase convection) and nucleate boiling. The influences of jet impact height and outlet velocity on local surface temperatures are pronounced at non-phase transition stage. In impact region, along with average superheat, the growth rates of heat transfer rate and liquid loss rate increase significantly from non-phase transition to nucleate boiling stage, and slow down at anaphase of nucleate boiling stage. Interestingly, the increase in outlet velocity delays the onset of slowdown in growth of liquid loss rate to some extent. Overall, the heat transfer rate and liquid loss rate are significantly enhanced at higher jet outlet velocity of 0.60 m/s, but not affected significantly by impact height. [ABSTRACT FROM AUTHOR]

• A revised self-induced jet impingement device is introduced in pool boiling using R1336mzz(Z). • Liquid-vapor separation promoted by the revised device is confirmed from the visualization results. • The revised device does nearly no effort on the nucleate boiling heat transfer coefficient until a relatively high heat flux. • Parametric investigations on the characteristics of the revised device are conducted. • The revised device gains up to 69.5 % and 36.7 % in the critical heat flux and the maximum nucleate boiling heat transfer coefficient, respectively. Pool boiling performance enhancement has attracted considerable interest on high power electronic cooling. In this paper, we introduce a self-induced jet impingement device for practical boiling enhancement, achieved by integrating the guidance tube with an orifice plate. With R1336mzz(Z) as the working fluid, we conduct visualization and parametric investigations on the pool boiling performance, considering the varying characteristics of this device. Visualization snapshots confirm the liquid-vapor separation and indicate that the two-phase flow pattern in the guidance tube eventually stabilizes as a churn-annular flow. Our findings suggest that self-induced liquid jet impingement has a minimal effect on the nucleate boiling heat transfer coefficient (h NB) under a certain heat flux, around 75 % CHF of the standard pool boiling, as the single-phase heat transfer is weakened due to the low thermal conductivity of R1336mzz(Z). As the heat flux intensifies, we observe substantial enhancements in both the Critical Heat Flux (CHF) and the Maximum Nucleate Boiling Heat transfer coefficient (h MNB) resulted from the delayed onset of boiling crisis, owing to the additional liquid supply and the promotion on vapor exhausting. The parametric studies reveal that the flow rate of the liquid jet, amplified by either extending the length or the inner diameter of the guidance tube, or by increasing the number of jet holes, positively impacts the boiling performance. However, these improvements show diminishing returns. Our study shows enhancements in CHF and h MNB of up to 69.5 % and 36.7 %, respectively, compared to the standard pool boiling. [ABSTRACT FROM AUTHOR]

• The scale inhibition mechanism of HAS was explored from three perspectives. • The main responsible for scale inhibition is carboxyl functional group of HAS. • There is a strong mutual attraction between HAS and the surface of calcite crystal. • The effect of HAS concentration on scale inhibition performance increases first and then decreases. The scale inhibition mechanism of sodium humate (HAS) on Calcium Carbinate (CaCO 3) was investigated by experiments, quantum chemical calculations and molecular dynamics simulations. The changes of electrochemical impedance spectroscopy were analyzed by electrochemical methods at different concentrations (0, 10, 30, 50, 70, 100 mg/L) of HAS. The molecular structure of HAS was optimized by quantum chemical calculation, and the charge distribution, electrostatic potential distribution and frontier orbital energy of HAS were obtained. The interaction simulation system between HAS and the crystal planes of calcite was constructed based on molecular dynamics simulation. The interaction mechanism between specific functional groups in HAS and calcite crystal plane was analyzed. The results show that the capacitance arc and charge transfer resistance increase with the increase of HAS concentration in the range of 0–50 mg. However, In the concentration range of 50–100 mg, the capacitance arc and charge transfer resistance decrease with the increase of concentration. In the optimized structure of HAS, the oxygen atom on the carboxyl group is more negatively charged. There is a strong mutual attraction between HAS and calcite crystal plane. The primary functional group responsible for achieving scale inhibition is formed through the bonding between the carboxyl oxygen atom within its molecular structure and the calcium atom on the calcite crystal plane. The research in this paper provides a strong guidance for the treatment of scale and the development of scale inhibitors. [ABSTRACT FROM AUTHOR]

• The usefulness of heat transfer ability of an object for heating/cooling is studied. • Useful entransy originates from entransy and describes useful heat transfer ability. • Heat transfer ability described by useful entransy fits heat/cold analysis better. • Useful entransy efficiency is more accurate and practical than entransy efficiency. • Useful entransy efficiency analysis and entransy dissipation analysis are complementary. The existing entransy theory does not consider the usefulness of heat transfer ability. This leads to the fact that the heat transfer ability described by entransy is not completely consistent with the heat/cold utilization, and the entransy efficiency analysis also has certain limitations. Hence, this paper first studies the usefulness of heat transfer ability. Under the ambient air with infinite heat capacity, only the exothermic ability of the object above ambient temperature T 0 is useful for heating, and only the endothermic ability of the object below T 0 is useful for cooling. On this basis, with T 0 as the reference point, the useful entransy is proposed to describe the useful heat transfer ability of an object for heating/cooling. It is demonstrated that heating/cooling mainly utilizes the heat transfer ability of an object relative to T 0 (useful entransy), rather than its heat transfer ability relative to 0 K (entransy). Moreover, taking latent heat storage as an example, it is proved that the useful entransy efficiency is more accurate than the entransy efficiency to evaluate the heat transfer irreversibility. And the guidance based on useful entransy efficiency is more practical. Furthermore, useful entransy efficiency analysis and entransy dissipation analysis are complementary. The combination of the two can evaluate and analyze the heat transfer performance more comprehensively. In summary, the useful entransy is an important improvement of entransy theory. [ABSTRACT FROM AUTHOR]

· The sheet-network and solid-network TPMS structures considering porosity and other parameters are constructed. · The solid-network TPMS structures with larger through-hole and skeleton concentration ratio have higher overall performance. · Under the same pressure drop, the heat transfer performance of TPMS with D-type sheet structure is better. Triply periodic minimal surface (TPMS) structures are widely used in the field of flow and heat transfer. In this paper, four TPMS structures, Schwartz Diamond-sheet (D-sheet), Schwartz Diamond-solid (D-solid), Schoen IWP-sheet (IWP-sheet), and Schoen IWP-solid (IWP-solid), are designed based on two structure generation strategies. The above structures all are manufactured through additive manufacturing, and the material is AlSi10Mg. The samples were characterized by X-ray scanning. The flow and heat transfer performance of different TPMS structures was studied by computational fluid dynamics (CFD) and experimental methods. The results show that the pressure drop of sheet-network is larger than that of solid-network for the same TPMS type. With the same structure generation strategy, the pressure drop of Schwartz Diamond (D-type) is greater than that of Schoen IWP (IWP-type). The heat transfer performance of the TPMS structure shows the same results as the pressure drop. The performance enhancement criterion (PEC) of the solid-network TPMS is better than that of the sheet-network structure, and the D-type TPMS structure is better than the IWP-type structure. The average PEC of D-solid and IWP-solid is 31 % and 44 % larger than that of D-sheet and IWP-sheet, respectively. The average PEC of D-sheet and D-solid is 11.8 % and 1.7 % larger than that of IWP-sheet and IWP-solid, respectively. The mechanism of the higher PEC of the solid-network structure is the low pressure drop caused by the high through-hole ratio and relatively concentrated skeleton. Under the same pressure drop, the heat transfer performance of TPMS with D-type sheet structure is better. The research content of the article is mainly applied to fin-enhancement, which can replace fins in traditional heat exchangers. In actual engineering applications, the TPMS structures type is reasonably selected according to the requirements. [ABSTRACT FROM AUTHOR]

• This paper clarifies the analogy principle of the conjugate heat transfer analysis and thermal stress analysis for film cooling. • Increasing the compound angle of the film hole increases the stress concentration factor while enhancing the cooling performance. • Changing the cylinder film hole into an elliptical film hole can both enhance the cooling performance and reduce the stress concentration factor. Film cooling is a commonly used cooling technology in gas turbines, which could reduce heat load by spreading cooling air. However, it faces the problem of stress concentration in structures with holes. The present work focuses on conjugate heat transfer and thermal stress characteristics for a flat film cooling plate through an experimentally validated numerical approach. It discusses the effect of three geometry factors: inclination angle, compound angle, and ellipticity. An engine case simulation and a laboratory case simulation are compared to validate the analogy principle of conjugate heat transfer and thermal stress analysis derived from theoretical analysis. According to the analogy principle, the laboratory flat plate film cooling cases are set to conduct numerical studies. The maximum Mises stress of the film hole is equal to the product of the reference stress and the stress concentration factor. The reference Mises stress is proportional to the overall cooling effectiveness. The stress concentration factor is mainly related to the geometry of the film hole. Reducing the inclination angle of the film hole from 90° to 30° increases the overall cooling effectiveness and drastically increases the stress concentration factor from 2.7 to 5.5. Based on a cylinder hole with a 30° inclination angle, increasing the compound angle enhances the overall cooling effectiveness and slightly increases the stress concentration factor. Based on a cylinder hole with a 30° inclination angle, increasing the ellipticity increases the overall cooling effectiveness and reduces the stress concentration factor. With the increase of ellipticity from 0 to 4, the stress concentration factor decreases from 5.5 to 2. Elliptical film holes are ideal for future low-stress, high-cooling-efficiency turbine blades to extend the life of cooling vanes. [ABSTRACT FROM AUTHOR]

• A correct formulation of average Nusselt number for the case of linearly heated wall is derived. • Thermal balance analysis provides the verification of the correct formulation. • Comparison of the existing and the correct formulation is included. The comment below addresses an inaccurate formulation in the calculation of the Nusselt number for linear heating as presented in the paper titled "Steady natural convection flows in a square cavity with linearly heated side wall(s)" by M. Sathiyamoorthy, T. Basak, S. Roy, and I. Pop, Int. J. Heat Mass Transfer 50 (2007) 766–775. Consequently, this analysis provides a mathematical reformulation for the local and mean Nusselt numbers of a linearly heated side wall. It regenerates a few quantitative results for a specific case solely intended for comparative purposes. [ABSTRACT FROM AUTHOR]

• Deep learning (DL) was used to predict supercritical flow field and heat transfer with drastically changing thermo-physical properties. • The deep neural network with a Multilayer Perceptron (MLP) model was established. • The relative errors between the DL prediction and the CFD calculation were within 4 %. • The prediction by the DL model was of several orders (106) faster than the CFD simulation. • DL method provided an efficient tool for design and rapid control of the regenerative cooling system. Flow dynamics and heat transfer of supercritical hydrocarbon fuels have significant influences on the regeneratively cooled advanced aero-engines. In this paper, based on the principle of deep learning (DL), a deep neural network with a Multilayer Perceptron (MLP) model was developed to predict the flow field and heat transfer characteristics of supercritical hydrocarbon fuels. The network was trained by the calculation of Computation Fluid Dynamics (CFD). The analysis shows that the predicted values of bulk fuel temperature, fuel velocity and heat transfer coefficient by the trained DL model match well with the CFD simulations. The linear correlation coefficients (R) between the DL model predicted values and the CFD calculated results are all bigger than 0. 99. Most of the relative errors are lower than 2 %. The consumed time of the DL model is of several orders (106) less than the CFD simulation. The effects of inlet Reynolds numbers, pressure and heating rate on the flow field and heat transfer of supercritical hydrocarbon fuels inside the regenerative cooling channel have been successfully learned by the DL model. Therefore, the proposed deep learning model provides an efficient tool for the sensitivity analysis and design optimization of the regenerative cooling system. [ABSTRACT FROM AUTHOR]

• RBF-DG was a novel prediction for absorption process in falling film and membrane based absorbers. • RBF-DG method provided excellent solutions. • Comparison with analytical solutions and experiments provided good agreements. • It becomes possible to predict accurately the absorption rate of falling film and membrane. The present investigation is intended to develop an understanding of how the vapor pressure and inlet temperature of the solution affect the absorption rate, which is helpful for the design of absorption heat pump systems. Many conjugated heat and mass transfer processes in falling film and membrane during lithium bromide (LiBr) absorption have not been well addressed. This paper proposes combining the meshless Radial Basis Function (RBF) and Discontinuous Galerkin (DG) methods to solve the conjugated heat and mass transfer in falling film and membrane-based absorption problems. The derivatives in the transverse direction are discretized using the RBF method. Owing to the parabolic nature of the governing equations, the DG methods can be used to discretize the derivatives in the streamwise direction. The efficacy of the proposed method is shown by comparing its results with analytical solutions and experimental data available in the literature. Furthermore, the proposed method is used to evaluate the analytical solutions of the heat and mass transfer problem in a membrane-based absorber. The previous published works simplified the problem to obtain the analytical solutions by providing the mean temperature and concentration values as the boundary conditions at the membrane-solution interface. The comparison of the RBF-DG method and the analytical solutions reveals that there are minor changes in temperature and local concentration values along the streamwise direction, but the profiles and trends are relatively similar across the film. A comparison of mass absorption rate calculations shows that RBF-DG and analytical solutions are compatible. The advantage of the numerical method is that various types of boundary conditions can be handled easily without any simplifications. Moreover, the numerical results of the proposed method are compared with experimental data, and they present good agreement. As a results, the absorption rate of falling film and membrane becomes capable of being accurately predicted. [ABSTRACT FROM AUTHOR]

• Performance of advanced turbulence models was compared. • Flow acceleration numbers are used to determine model applicability. • The structure effect of buoyancy in mixed convection cannot be ignored. • Provide reference for the selection and improvement of turbulence model. The extreme thermophysical property changes of fluids near the pseudocritical point lead to strong buoyancy and acceleration effects, posing challenges for turbulence modelling. This paper aims to assess the performance and applicability of various turbulence models, with the intention of providing a useful reference for quickly selecting appropriate models and for models' further improvement. Additional Favre-averaged turbulence models were implemented through user-defined scalars in ANSYS Fluent and compared to published experimental and high-fidelity simulation data. This study found that the tested turbulence models have different sensitivities to molecular Prandtl numbers, buoyancy effects, and flow acceleration. For forced convection, all turbulence models can be applied in the case of the flow acceleration number Kv < 1.37 × 10 − 8 , and the ZXY model has the widest application range. None of the tested turbulence models showed generalizability under mixed convective conditions. The Lien-PSA model qualitatively reproduced the wall temperature trends for all tested cases. [ABSTRACT FROM AUTHOR]

• 3D micro-droplet printing of electronic devices based on heterogeneous wettable surfaces. • Manipulation of wetting differences to obtain four droplet equilibrium contact line morphologies. • While ensuring the 3D printing accuracy, the technique effectively curbs its rebound and oscillatory actions. • This strategy effectively reduces the droplet axial height and improves the droplet heat transfer capacity. Microdroplet printing technology can effectively realize the precision preparation of complex electronic devices at the micro-scale, relying on the mutual interaction between the tiny micro-droplets, which can improve the organizational microstructure of the fabricated parts and effectively enhancement the mechanical properties of the fabricated parts. In this paper, the interaction process between micro-droplets and patterned wetted modified surfaces during the microdroplet printed process is studied based on the coupled level collective integral number method (CLSVOF) and equivalent heat capacity method. The influence of patterned schemes for non-uniform surface wetting modification, regional wetting differences and impingement process parameters on the impingement spreading flow of metal droplets and their heat transfer cooling process is investigated. The evolution mechanisms of the droplet geometry and the distribution characteristic of the droplet internal temperature and heat flow density during the impingement spreading process are discussed. The optimization mechanism of heterogeneous wetting modified surfaces for molten droplet forming is revealed. The results show that the four-way homogeneous pattern modification scheme can effectively improve the microdroplet stabilization rate and enhance the preparation efficiency by ensuring the shape size and positioning accuracy. Four three-phase contact lines morphology patterns for the equilibrium state of droplets, namely, circular, rectangular, concave shape and rhombic, were obtained by adjusting the wetting differences between the areas. A scheme to regulate the surface behavior of printed microdroplets with lyophilic high-adhesion region to anchor the droplets and lyophobic low-adhesion environment to provide a physical barrier to confine the droplet morphology was developed. In addition, the pattern modification scheme has an obvious optimization effect for the droplet forming and its heat exchange process in the process of metal droplet printing. Basis on the ensuring the radial dimensional accuracy and sphericity of the molten droplet molding, the rebound and oscillation behavior of the molten droplet is effectively suppressed. The stabilizing efficiency of droplet printing is greatly enhanced and further improves the heat transfer performance of the droplet. The results further enrich the theory of precise regulation of the dynamic behavior of droplet surfaces. Not only can it be used as a new strategy for universal high-precision printing and droplet control, but it will also promote the application and development of technologies in numerous fields, such as liquid transport and microfluidics. [ABSTRACT FROM AUTHOR]

• The chemical source term is the main control factor of thermal runaway. • The higher the temperature intensity, the shorter the explosion delay time. • The critical temperature of thermal runaway for AP/HTPB propellant is 458 K. • If the propellant is thick enough, the temperature distribution remains constant. For the safe use of solid rocket motors, it is crucial to investigate the thermal effect of AP/HTPB propellants. Transient high temperature may cause AP/HTPB propellants to burn or explode, resulting in serious damage. The fundamental equations that describe the thermal effect of AP/HTPB propellants are the heat conduction differential equation and Arrhenius law. In this paper, the thermal safety of AP/HTPB at transient high temperature is obtained by solving the heat conduction differential equation. The results show that the chemical source term S (T) is the main control factor of thermal runaway. When the action duration of the external high-temperature is short (ms-s magnitude), the first explosion occurs closest to the high-temperature source. When the action duration of the external high-temperature action is long (min-h magnitude), the first explosion occurs inside the propellant. The higher the high-temperature intensity of external action, the shorter the explosion delay time, but the explosion temperature does not change accordingly. The critical temperature of thermal runaway is 458 K under the condition of no shell and a high-temperature action duration of 2 h. At high-temperature intensities of 458 K, 459 K, and 500 K, the propellant explosion temperatures are 533 K, 539 K, and 540 K, and the explosion delay time is 1592s, 1102s, and 160 s, respectively. The critical duration of thermal runaway is 430 ms, 650 ms, 1350 ms, and 9820 ms, and the stable explosion delay time is 450 ms, 680 ms, 1400 ms, and 9920 ms, respectively, under the conditions of 2000 K, 1500 K, 1000 K, and 600 K without the shell. The propellant explosion temperature remains constant as the shell thickness grows while the explosion delay time gradually increases. When the propellant is thick enough, the temperature distribution of the system remains constant. The results of this study may be useful as a guide for assessing the thermal safety of AP/HTPB propellant. [ABSTRACT FROM AUTHOR]

This paper focuses on designing surrogate models that have uncertainty quantification capabilities to effectively improve the thermal performance of rib-turbulated internal cooling channels. To construct the surrogate, we use the deep operator network (DeepONet) framework, a novel class of neural networks designed to approximate mappings between infinite-dimensional spaces using relatively small datasets. The proposed DeepONet takes an arbitrary rib geometry as input and outputs continuous detailed pressure and heat transfer distributions around the profiled ribs. The datasets needed to train and test the proposed DeepONet framework were obtained by simulating a 2D rib-roughened internal cooling channel. To accomplish this, we continuously modified the input rib geometry by adjusting the control points according to a simple random distribution with constraints, rather than following a predefined path or sampling method. The studied channel has a hydraulic diameter, D h , of 66.7 mm, and a length-to-hydraulic diameter ratio, L / D h , of 10. The ratio of rib center height to hydraulic diameter (e / D h), which was not changed during the rib profile update, was maintained at a constant value of 0.048. The ribs were placed in the channel with a pitch-to-height ratio (P / e) of 10. In addition, we provide the proposed surrogates with effective uncertainty quantification capabilities. This is achieved by converting the DeepONet framework into a Bayesian DeepONet (B-DeepONet). B-DeepONet samples from the posterior distribution of DeepONet parameters using the novel framework of stochastic gradient replica-exchange MCMC. Finally, we demonstrate the performance of the proposed DeepONet-based surrogate models with uncertainty quantification by incorporating them into a constrained, gradient-free optimization problem that enhances the thermal performance of the rib-turbulated internal cooling channel. • Develop a deep operator learning-based surrogate model. • Develop a Bayesian DeepONet model with uncertainty quantification capability. • Optimize internal cooling channel rib profiles using the surrogate model. [ABSTRACT FROM AUTHOR]

• Hydrodynamic instability will induce natural frequencies in steady flow and Kelvin-Helmholtz instability appears only at low Reynolds numbers. • A periodic process of vortex generation-shedding-fusion is found behind each block in pulsating flow. • As the forcing frequency St F of the pulsating flow equals the natural frequency St N in steady flow, pulsating flow improves heat transfer, especially at low Reynolds numbers. This paper presents a numerical investigation on flow and heat transfer in a grooved channel formed by 9 rectangular blocks with steady and pulsating flow. Large eddy simulation is used to simulate the hydrodynamic instability in steady flow and heat transfer enhancement for pulsating flow. The simulations for steady flow at Re =700, 1500, 2500, are firstly conducted, showing that there are 23Hz, 70Hz and 130Hz natural frequencies f N induced by hydrodynamic instability (corresponding to Strouhal number St N =0.876, 1.243, 1.385) behind the first block, which disappear quickly with the flow development. A natural frequency f N of 34Hz (St N =1.295) is found only at Re =700 and it can be maintained and amplified by groove structures. Subsequently, a sinusoidal pulsating flow is applied, with the forcing frequency St F equal to the natural frequency St N in steady flow and the dimensionless amplitude A fixed at 0.2. Flow and thermal performance of the grooved channel is studied under the pulsating condition. A periodic process of vortex generation-shedding-fusion is found behind each block and the fluid in the channel has higher turbulent kinetic energy than that under the steady condition. Ultimately, compared with the thermal and pressure-drop performance between steady and pulsating flow, the enhancement of time-average Nusselt number ranges from 48.08%, 15.64%, 8.61% and time-average friction factor increases by about 40% at Re =700, 1500, 2500. The thermal performance factor is increased by 30.96% and 4.20% at Re =700, 1500 but is reduced slightly by 3.48% at Re =2500. It indicates that the augmentation of heat transfer through pulsating flow is significant at low Reynolds numbers. [ABSTRACT FROM AUTHOR]