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]

• 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]

• Some ambiguities in experimental test results. • Some ambiguities in the simulation results. • Lack of definition of some symbols. This paper presents some comments on "Chen, Bo, et al. "DPM-LES investigation on flow field dynamic and acoustic characteristics of a twin-fluid nozzle by multi-field coupling method." International Journal of Heat and Mass Transfer 192 (2022): 122927. "We thank the authors for conducting this valuable research. However, some points should be noted to improve their paper. In general, there are some serious concerns about the accuracy and validity of experimental and simulation results. In addition, some less important things such as the lack of definition of some symbols in their paper should be corrected. Details are mentioned below in comments 1-8. [ABSTRACT FROM AUTHOR]

• Flame length of moving fires is 15.2%-45.7% more than that of stationary fires. • Temperature upstream moving fire sources is significantly elevated. • Entrainment and heat transfer mechanism of moving and stationary fires compared. • Backwards heating behavior of moving fires can bring threatening thermal hazard. Moving fires have been paid great attention by researchers in recent years, while the behavior of this special combustion phenomenon has not been fully investigated. In this paper, a 1:10 scale burning car was innovatively designed based on Froude similarity criterion, and a series of moving model experiments were conducted. The flame morphological characteristics and heating behaviors between the stationary and moving fires were thoroughly compared. The results manifest that the distorted flame of moving fires is 15.2%-45.7% more than that of stationary fires, while this behavior is not accurately predicted by current models for wind-blown and wall-attached flames, which suggests that the wind tunnel test has limitations in studying moving fires. Besides, compared with stationary fires, the temperature upstream of the fire source is considerably elevated, and the heat transfer from the flame to the top surface of the burning car is dramatically enhanced. A physical model is proposed to describe the flame distortion and backwards heating behaviors in moving fires, which reflects the important role of the air entrainment and the cool fuel-rich cone in the combustion process. This paper represents a first step toward understanding the behavior of moving fires, which can provide a framework to deeply investigate the combustion mechanism of moving fires in the future. [ABSTRACT FROM AUTHOR]

• In this paper is provided an calculation expression for heat transfer during condensation inside tubes that provides an adequate fit with experimental data available from a total of 22 fluids, including refrigerants, water, and a variety of substances organic. This paper presents a mathematical deduction of a new improved model for heat transfer during condensation inside tubes. This new model has been developed with the aid of the Gaussian Equation for an infinite straight line, considering this relation with the differential equations that govern the heat transfer process. The proposed model was verified by comparison with available experimental data of 22 different fluids, including various refrigerants, water, and organic substances, which condense inside vertical, inclined and horizontal tubes. The proposed model is valid for an reduced pressure values ranging from 0.0008 to 0.91, values of Reynolds number for single-phase between 68 and 84,827 and for Reynolds number for two-phase between 900 and 59,4373, a range of internal diameters ranging from 2 to 50 mm, vapor quality from 0.01 to 0.99, P r values for single-phase from 1 to 18 and mass flux rates in the ranges of 3 to 850 kg/(m2s). The mean deviation found for the analyzed data for horizontal tubes was 11.8%, while for the vertical and inclined tubes data the mean deviation was 13.0%. In all cases, the agreement of the proposed model is good enough to be considered satisfactory for practical design. [ABSTRACT FROM AUTHOR]

• An analytical model for heat transfer inside deep borehole heat exchangers (DBHE). • Convolution theorem to address vertically increasing borehole wall temperature. • Relationships between DBHE thermal performance and various parameters. • An effective tool fast and easy enough for the design and optimization of DBHE. Deep borehole heat exchangers (DBHE) provide an effective solution for ground coupled heat pump (GCHP) systems in cold climate region where heating is dominant. Concerning the analytical heat transfer models, the simplification that borehole wall temperature is constant along the depth of ground heat exchangers (GHE) which the existing quasi-three-dimensional models have assumed in the application of shallow borehole GHE, can no longer be accepted in the application of DBHE due to the geothermal gradient in deep ground. Making this simplification cannot give the real temperature distribution of circulating fluid along the depth of DBHE. Therefore, this paper developed a new analytical model for DBHE with coaxial pipes by successfully addressing the increasing borehole wall temperature using the convolution theorem, so that the widely employed quasi-three-dimensional models for shallow borehole GHE is extended for DBHE with coaxial pipes for the first time. The new analytical model is validated by comparing with an existing numerical model. Using the newly developed analytical model, the trends revealing the relationships between thermal performance of DBHE and various parameters are firstly plotted. Because of the high accuracy and quick calculation, this new analytical model can be used as a benchmark for numerical models. More importantly, the proposed analytical model can be an effective tool for the design and optimization of DBHE, since current numerical models are always calculation-demanding, time-consuming and difficult for engineers and designers to use. Also, the method this paper proposed to address the varying borehole wall temperature can certainly be employed to improve the existing quasi-three-dimensional models for shallow borehole GHE so that they can be applicable in some other cases, for example, GHE installed in layered soils, or affected by underground seepage flow in partial depth. [ABSTRACT FROM AUTHOR]

• A theoretical framework is proposed to consider non-equilibrium configurations at the macroscopic scale. • The new formulation is validated using a full dual-scale model of coupled heat and mass transfer. • With the new model, a very good accuracy is obtained with a huge reduction of CPU time. • Dual-scale diffusion and molecular relaxation can be considered additively in the modified macroscopic formulation. • This framework can be applied to other configurations, such as coupled transfer in packedbeds. This paper proposes a macroscopic formulation of coupled heat and mass transfer that can consider non-local equilibrium often encountered in biosourced building materials (wood- and plant-fiber based materials). Transferring dual-scale effects and molecular relaxation at the macroscopic level involves a kernel function acting in a convolution product. To ease the computational solution of the set of equations, the memory function is decomposed as a series of exponential functions. Each function yields an internal variable that obeys a simple ordinary differential equation (ODE). This paper first describes the macroscopic and dual-scale formulations used as reference solutions. Subsequently, the modified macroscopic formulation of coupled transfer and its computational solution are presented in detail. The major outcomes of the present study, validated against reference solutions obtained with a comprehensive dual-scale model, are as follows: • Dual-scale diffusion can be approached accurately by two exponential functions, • Even though the dual-scale phenomenon and molecular relaxation do not occur at the same scale, both can be considered in the modified macroscopic formulation of coupled transfer additively, • The new macroscopic formulation, together with the computational procedure proposed in this study, can be applied to various configurations, namely coupled heat and mass transfer in packed beds. [ABSTRACT FROM AUTHOR]

• This paper presents a thermodynamic assessment of different low GWP refrigerants. • The thermodynamic assessment involves an experimental and a theoretical analysis. • The theoretical analysis is based on the Total Temperature Penalisation (TTP) criteria. • R32 is an effective substitute for R410A, R290 and R1270 are valuable substitutes for R404A. • R152a, R1234yf and R1234ze(E) are effective substitutes for R134a and R1234ze(Z) is the most effective substitute for R236fa. This paper presents the heat-transfer assessment of different low GWP substitutes for traditional HydroFluoroCarbon (HFC) refrigerants obtained by applying both an experimental analysis based on the direct measurement of thermal (heat transfer coefficients) and hydraulic (pressure drops) performances and a theoretical analysis based on a specific Performance Evaluation Criteria (PEC), the Total Temperature Penalisation (TTP), in the specific case of boiling and condensation in a commercial Brazed Plate Heat Exchanger (BPHE). The results of both the experimental and the theoretical assessment confirm that the refrigerant R32 has comparable heat transfer characteristics to R410A. Similarly the refrigerants R290 and R1270 exhibit superior heat transfer performance to R404A; the refrigerants R152a, R1234yf and R1234ze(E) show heat transfer characteristics comparable to R134a and R1234ze(Z) exhibits superior heat transfer performance to R236fa. The coupling of experimental and theoretical assessments provides a sound procedure for the heat-transfer assessment of the low GWP alternatives for traditional HFC refrigerants. [ABSTRACT FROM AUTHOR]

• Hypersingular heat flux interface integral equation is derived. • Interface fluxes are calculated for multi-medium heat transfer analysis using IIBEM. • IIBEM is a promising BEM for solving complicated practical engineering problems. In this paper, based on a direct method for evaluating hypersingular integrals, a heat flux interface integral equation is presented to calculate the hypersingular interface heat fluxes in multi-medium heat transfer analysis by using the interface integral boundary element method (IIBEM). The recently proposed IIBEM has a significant advantage of a single integral equation method for solving multi-medium problems, which makes up for the lack of boundary integral equations to solve multi-medium problems. Compared with conventional multi-domain boundary element method (MDBEM), the IIBEM is more efficient in input data preparation, program coding and computational time. However, a big issue encountered is the determination of interface heat fluxes in multi-medium heat transfer analysis through using IIBEM, since the derived flux interface integral equation involves hypersingular integrals on the interface, and no additional information can be utilized to recover the heat fluxes by using the in-direct "traction-recovery" methods. In this paper, based on a new efficient direct method for numerical evaluation of high order singular curved boundary integrals, the hypersingularities in flux interface integral equation are precisely evaluated, and based on the successful implementation of flux hypersingular integral equation, the interface heat fluxes are firstly calculated in the multi-medium heat transfer analysis by using IIBEM. Several two- and three-dimensional numerical examples are given to validate the correctness of the presented method. [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]

Highlights • Constructal design of shell-and-tube evaporator is conducted. • Evaporator is used in Kalina cycle with ammonia-water working fluid. • A complex function minimization is taken as optimization objective. • It composed of heat transfer rate and total pumping power. • External diameter of the heat transfer tube and mass flow rates are optimized. Abstract On the basis of constructal theory, structure design of a shell-and-tube (ST) evaporator with ammonia-water working fluid is conducted in this paper. A complex function composed of heat transfer rate (HTR) and total pumping power (TPP) is taken as optimization objective, and the heat transfer area of the tubes is taken as the constraint. The complex function of the ST evaporator is minimized, and the corresponding optimal external diameter (ED) of the heat transfer tube (HTT) and optimal mass flow rates (MFRs) of the ammonia-water and hot water are obtained. Compared with the initial design, the complex function, HTR and TPP are decreased by 17.24%, 22.88% and 64.14% after constructal optimization, respectively. It illustrates that the complex function is a tradeoff between the heat transfer and fluid flow performances of the evaporator. As the ammonia concentration increases, the overall performance of the evaporator is improved. The inlet temperature of hot water has certain influences on the optimal MFRs of the ammonia-water and hot water, but has no influence on the optimal ED of the HTT. The optimal results obtained in this paper can provide some guidelines for the optimal designs of the ST evaporators, and can be further extended to the optimal designs of the ST condensers and even the whole systems of Kalina cycles. [ABSTRACT FROM AUTHOR]

Highlights • A new periodical solution with reverse-central-symmetrical-structure has been discovered. • The flow is divided into seven stages by ten critical buoyancy ratios. • Dynamical behavior, flow patterns and heat and mass transfer efficiencies in each stage are investigated. • The governing equations are discrete by using high-order high-resolution schemes. Abstract The paper investigates numerical results of double-diffusive convection in a slim rectangular enclosure with horizontal temperature and concentration gradients. The case with the thermal Rayleigh number, Ra T = 10 5 , the Prandtl number, Pr = 1 , the Lewis number, Le = 2 and height-to-width aspect ratio, A = 4 is considered. Effects of the buoyancy ratio of composition to temperature, N , on dynamical behavior, flow structure and heat and mass transfer efficiency are investigated carefully. The present results suggest that with the buoyancy ratio N increasing from 1 to 2, the flow of transition is a complex series of changes from the periodic to quasi-periodic, chaotic, periodic, and finally to laminar with concentration dominated. There are three types of periodic flows, the first one is simple periodic flow with single fundamental frequency (FF), the second one is periodic flow with FF and its integer multiple of frequencies (INTMF), the last one is the more complex periodic flow with FF, INTMF and its sub-harmonic frequency multiple (SFM). A new type periodical solution with reverse central symmetrical structure is captured in the present paper. Multi cells structures traversing between two vertical walls will produce salt figure phenomenon and benefit to enhance heat and mass transfer efficiency. [ABSTRACT FROM AUTHOR]

Highlights • A 30–50 K deterioration of thermal insulation performance of functionally graded thermal barrier coatings (FGTBCs) compared to traditional TBCs is reported. • Thermophysical properties of Functional Graded Materials(FGMs) are experimentally determined. • A comprehensive heat transfer model was built to predict the stationary temperature distribution and heat flux inside FGTBCs. Conduction heat flux and thermal radiation flux were analyzed separately. Abstract To further study the heat transfer process inside the functionally graded thermal barrier coatings (FGTBCs) beyond the companion part of the present paper concerning their radiative properties, this paper takes both conductive and radiative heat transfer into account and establishes a comprehensive heat transfer model considering both conductive and radiative heat transfer to estimate the thermal insulation performance of FGTBCs. Thermal conductivities of more than thirty yttria-stabilized zirconia (YSZ)/NiCoCrAlY FGTBC samples were determined based on their thermal diffusivities determined by the laser flash analysis. To predict the overall thermal insulation performance of 8YSZ/NiCoCrAlY FGTBCs, a comprehensive heat transfer model based on the two-flux method and multi-band model was built to solve the stationary temperature distribution and heat flux inside thermal barrier coatings (TBCs) and FGTBCs. Results indicate that FGTBCs show worse thermal insulation performance compared to conventional YSZ TBCs; we also investigated the influence of the number of layers. The temperature of engine blade is increased by 30 – 50 K if coated with an FGTBC instead of a conventional TBC. Conduction heat flux is significantly greater than radiative flux for both TBCs and FGTBCs. This paper gives a detailed analysis by simultaneously considering the heat conduction and thermal radiation in FGTBCs, which is helpful for achieving better thermal design of YSZ/NiCoCrAlY FGTBCs. Along with the Part I, this two-part article presents a deep and thorough study on the thermal insulation performance of YSZ/NiCoCrAlY FGTBCs. [ABSTRACT FROM AUTHOR]

Highlights • The theory of Karman vortex street was adopted to explain the reason of weak region formation. • Two tube configurations of which the section was semicircular and semielliptical were discussed in this paper. • The velocity vector, the flow streamline and the temperature contour of the modified configurations were displayed. Abstract The main objectives were to develop a method for designing an effective finned heat exchanger through changing the configuration to ensure the tubes be streamlined shape. Kinds of inserts including, vortex generator, louver fin, and serrated fin in the heat transfer was used to increase the effect of thermal performance, but also was harmful to the flow characteristics. Consequently, the comprehensive performance was obviously lower than the heat exchangers with no inserts. In order to solve this problem, streamlined tube configuration was adopted in this paper to search a new method to enhance the comprehensive performance by reducing the Karman vortex street effect. Two modified configurations of tube were investigated in this paper. The modified configurations were both made of semicircular tube and semielliptical tube. The simulation was validated by Wang's experiment and Arafat's simulation. The numerical data were presented at Reynolds number in the range between 1300 and 2000. Based on the numerical results, the effect of elliptical axis ratio on the thermal hydraulic performance was discussed. The results showed that the comprehensive performance of the modified configurations was greater than that of the circular tube. In addition, the velocity vector, flow streamline and temperature contour of the configurations were displayed. The Karman vortex street effect was adopted to explain the reason of wake region formation. It can be concluded that long and thin section of heat transfer tube, which was coincident with streamlined design, was beneficial to improve the comprehension performance of the heat exchanger. [ABSTRACT FROM AUTHOR]

Highlights • Existing EHD convective boiling heat transfer correlations are for specific geometries, fluids and test conditions and are fully empirical having no physical basis. • The performance of EHD convective boiling devices, in particular pressure drop, remains largely unpredictable. • Phenomenological correlations for EHD convective boiling heat transfer coefficient and pressure drop, based on two phase flow pattern, are presented. • Considerable improvement over previous EHD convective boiling performance correlations was found both in error and the physical significance of coefficients. Abstract EHD is an active heat transfer enhancement technique for convective boiling or condensing dielectric fluids whose performance varies widely in the literature for even the same geometry and fluids. As such, the performance of EHD convective boiling devices remains largely unpredictable. Two empirical EHD convective boiling heat transfer coefficient correlations exist in the literature with good performance when compared to their exact test dataset and geometry, however they have been shown to have poor correlations in predicting the performance for external test datasets. In this paper, a phenomenological approach is taken in the development of a new EHD convective boiling heat transfer coefficient correlation and a novel EHD convective boiling pressure drop correlation. The performance correlations are based on widely-used free-field convective boiling correlations for heat transfer coefficient and pressure drop with phenomenological enhancement factors based on the two phase flow pattern as predicted using the EHD two phase flow pattern map which accounts for the effect of electric field strength on flow pattern redistribution, as described in a previously published paper (Nangle-Smith and Cotton, 2018). A review of the EHD convective boiling experimental literature was conducted to determine confounding factors that can impact the large variance in the reported performance. Flow pattern and applied heat flux were identified as common parameters not maintained constant in the experimental data. Two datasets, including data from the present study, in which flow pattern and applied heat flux are maintained constant were used in the development of the performance correlations. Considerable improvement over previous EHD convective boiling performance correlations was found both in error and the physical significance of coefficients. The correlations were developed for thermodynamic qualities in the range of 20–60% below the onset of dryout or mist flow patterns and for applied heat fluxes < 30 kW/m2. Furthermore, test conditions in this study were chosen to focus on the dielectrophoretic effect and therefore studies with significant electrophoresis were not included. Thus, the correlation is limited to positive applied voltages, saturated flow boiling conditions, and field strengths below the onset of charge injection. The authors recommend this approach for EHD two phase performance correlation development as it is more mechanistic, is analogous to the state-of-the-art approaches for free-field two phase performance, and is likely to yield more accurate models as more experimental data becomes available. [ABSTRACT FROM AUTHOR]

Highlights • Method of images is successfully adopted to predict the LVG vortex dynamics. • A mixed LVG configuration is suggested to enhance the heat transfer performance. • LVG aspect ratio and channel height are very critical to the LVG effectiveness. Abstract This paper reports a numerical investigation of fluid flow and heat transfer in a rectangular channel with delta-shaped winglet longitudinal vortex generators (LVGs) under different configurations. In addition to the conventional common-flow-down configuration and the common-flow-up configuration, a unique mixed configuration is suggested. The "method of images" is successfully adopted to analyze the dynamics of the longitudinal vortices due to wall interference. The Nusselt number, friction factor and overall performance coefficient for the three configurations are compared at various Reynolds numbers (all less than 2200), LVG row numbers, channel heights, and LVG aspect ratios. It is found that the channel height and LVG aspect ratio are the two most critical factors influencing the effectiveness of the different LVG configurations, which may explain why inconsistent conclusions have been presented by previous studies comparing the performance of the common-flow-down and common-flow-up configurations. When the channel height is relatively small, the unique mixed configuration is found to be the most effective at enhancing fluid mixing (and therefore improving heat transfer) at streamwise cross sections. When the LVG aspect ratio becomes large, the common-flow-down configuration outperforms the common-flow-up configuration. This paper sheds some insight on how to design the optimal LVG configuration for enhancing heat transfer performance. [ABSTRACT FROM AUTHOR]

Highlights • Composite fin with high thermal conductivity coatings deliver large heat transfer enhancement. • At high Re typical of Engine Section Stator vanes, the coating thickness should be around 100 μm. • Composite vanes would be more effective at low Re, as predicted also for composite fins. • Using the internal web structure of vanes can improve heat transfer in long vanes. Abstract This paper main goal is aimed at a paradigm shift in the enhancement of heat transfer rate between finned surfaces and surrounding fluid by presenting a novel approach in composite fins. This approach consists in using high thermally conductive coatings on top of the finned substrate in order to increase the local temperature along the fin washed surface. In the present paper a high thermal conductivity coating has been applied to an aeroengine vane of different shape and dimensions subject to icing conditions at high Reynolds numbers and where the main aim was to keep the vane warm as to avoid icing in aerospace applications. Numerical simulations have been carried out to ascertain the range of thickness of the coatings to be used to maximise the wanted effect. Both short Engine Section Stator (ESS) and longer fan Outlet Guide Vane (OGV) have been modelled, having different final goals. In the case of longer OGV, an additional novel design modification has been suggested to enhance further the heat transfer along the vane length by the use of the internal webs. The experimental validation also carried out at much higher Reynolds numbers than that reported in Buffone et al. (2005), demonstrate that the novel concept of heat transfer enhancement in composite fins is a simple and yet powerful strategy in a wide range of Reynolds numbers. A fin analysis has been performed of both the present ESS vane the fins tested in Buffone et al. (2005) at much lower Reynolds numbers and shows that the improvement obtained with the coated fins tested in Buffone et al. (2005) is much larger than the coated ESS vanes investigated in the present study. This said, the present study demonstrates that the use of high conductive coatings in composite fins can keep the ESS vanes ice free, something that was not possible with uncoated vanes. It is important to note that the actual optimal thickness of the thermally conductive coatings is a function of Biot number, fin shape, dimensions and thermal conductivities of fin substrate and coating; depending on the application, a proper design of the fin substrate and coating should be carried out. [ABSTRACT FROM AUTHOR]

Highlights • Advection-diffusion problems with multi-term time-fractional derivatives and of variable coefficients are considered in higher dimensions. • These problems have applications mainly in heat conduction or heat transfer, mass and chemical transports into porous or nonporous media. • A meshless spectral radial point interpolation method (MSRPIM) using radial basis functions (RBFs) and point interpolation, is proposed and implemented for obtaining accurate solutions of these problems. • Stability analysis of the proposed method is carried out which is justified computationally. • Computed solutions obtained via different RBFs are matched with each other as well matched with the exact and earlier works for various fractional order derivatives. Comparisons reveal accuracy and efficiency of the proposed method in obtaining accurate solutions for such problems. Abstract In this paper, a weighted meshless spectral radial point interpolation method (MSRPIM) is implemented for the numerical solutions of a class of multi-term time fractional advection-diffusion problems with spatial and temporal dispersion. In the proposed methodology, meshless shape functions are used for discretization in space. These functions are generated via radial basis functions (RBFs) and point interpolation method. Weighted θ -scheme and first order finite differences together with a simple quadrature rule are employed for the discretization of time fractional derivatives. Efficiency and accuracy of the formulated method is validated in terms of number of nodal points N , time-step size τ and weight parameter θ ∈ [ 0 , 1 ]. Quality of the computed solutions are checked using E ∞ , E 2 and E rms error norms. Stability of the proposed method is analyzed theocratically and validated computationally, which is an important task of the current paper. Simulations reveal very good accuracy for both 1D and 2D test problems having constants and variable coefficients. [ABSTRACT FROM AUTHOR]

Graphical abstract Highlights • Flow boiling of R1233zd(E) exploring different operating conditions. • HTC increase with mass flux. Low effect of saturation temperature and heat flux. • HTC and dry-out well predicted by Wojtan et al. [1,2] methods. • Asymmetry on annular flow heat transfer dependent on vapor Froude number. • Modified Cioncolini and Thome [3] model to catch also asymmetric conditions. Abstract Flow boiling heat transfer and dry-out vapor quality data for the low-GWP refrigerant R1233zd(E) are presented in this paper. Tests have been performed in a single horizontal stainless steel tube with an internal diameter of 6.0 mm for the whole range of vapor qualities (from 0.008 to 0.977) by varying the mass flux from 147 to 300 kg/m2 s, the heat flux from 2.4 to 40.9 kW/m2 and the saturation temperature from 24.2 to 65.2 °C, corresponding to a reduced pressure range of 0.036–0.125. The assessment of some available predictive methods has been performed: the best ones for the heat transfer and the dry-out predictions are both by Wojtan et al. (2005). 68.3% of the data points are falling into an error band of ±30% for the heat transfer coefficients, while the MAE for the dry-out predictions is 14.4%. The range of operating conditions explored allowed to segregate the almost pure convective boiling results from those affected by nucleation. The effect of the non-uniform liquid film distribution along the perimeter has been related to the Froude vapor number: for Fr V > 3 the pure convective boiling data are symmetric around the tube and the correlation of Cioncolini and Thome (2011) returns the best agreement (MAE = 25.4%). By using the local expression of the authors' original model, a new equation has been developed in order to capture the experimental trend also for asymmetric flows. The corrected Cioncolini and Thome (2011) method provides an overall MAE of 7.40% when compared to the entire annular flow database of this paper. [ABSTRACT FROM AUTHOR]

Highlights • A quantitative study of the frost distribution and growth on FTHX was reported. • Frost pattern between windward and leeward sides of FTHX was compared. • Frost thickness, density and mass on different parts of FTHX were analyzed. • The effects of frosting conditions on frost distribution and growth were studied. Abstract Recently, the authors have published a paper which provided the frost distribution and growth along the airflow direction of finned tube heat exchangers (FTHXs), and demonstrated that some frost was accumulated on the edge of windward fins at a given frosting condition. As a follow-up to the work reported in the previous paper, in the current paper, the effects of frosting conditions on frost distribution and growth characteristics of FTHXs with different fin pitches were experimentally investigated. The results showed that the protruded thickness of frost on the edge of windward fins and the ratio of the frost mass on the edge of windward fins to that on the entire FTHX were both increased with the increases in air temperature (T a) and relative humidity (RH a), but decreased with an increase in initial air face velocity (v a). In addition, the effects of RH a on the mass ratio were more significant than those of T a and v a. Furthermore, the averaged frost density on the edge of windward fins was increased with an increase in RH a ; and that on the surfaces of fins and tubes was increased with an increase in v a , but decreased with an increase in RH a. [ABSTRACT FROM AUTHOR]

Highlights • Microchannel heat sinks with triangular ribs on sidewalls. • Three-dimensional numerical model of laminar flow and heat transfer. • The triangular ribs effectively prevent the decline of local heat transfer. • Different arrangements result in different local thermohydraulic characteristics. • Geometry of rib shows significant influence on local thermohydraulic performance. Abstract Triangular ribs are mounted on the parallel sidewalls of microchannels in order to reinitialize the thermal boundary layer and improve the mixing of cold and hot fluids. This paper presents a detailed numerical study on local laminar fluid flow and heat transfer characteristics in microchannel heat sinks with tandem triangular ribs for Reynolds number of 443. Three-dimensional conjugate heat transfer models considering entrance effect, viscous heating, as well as temperature-dependent thermophysical properties are employed. Water and silicon are respectively used as fluid and solid for the computational domain. Triangular ribs are attached in microchannels with either aligned or offset arrangement. Four non-dimensional geometry parameters relative to the width, height, converging-diverging ratio and spacing of triangular ribs are proposed to investigate the influence on local fluid flow and heat transfer characteristics. Velocity contour, pressure and temperature distributions are examined to demonstrate the basic fluid flow and heat transfer mechanism. Local pressure and temperature profiles are studied to show the influence of the triangular ribs on fluid flow and heat transfer process. Local friction factor and Nusselt number for different non-dimensional geometry variables are further investigated to comprehensively indicate the impact of triangular ribs. Results shows that the triangular ribs can significantly reduce the temperature rise of the heat sink base and efficiently prevent the drop of local heat transfer coefficient along the flow direction, but also result in higher local friction factor than the straight microchannel. For the studied operation conditions and geometry parameters of flow passage, the heat sink base temperature varies in the range of from 301.90 to 324.31 K, the computed pressure drop and heat transfer coefficient fluctuate from one triangular rib to the next and their amplitude and wavelength significantly depend on the geometry and arrangement of triangular ribs. Compared to the reference straight microchannel heat sink, a superior configuration considered in this paper can yield an improvement of up to 2.15 times higher of average Nusselt number. [ABSTRACT FROM AUTHOR]

Highlights • Turbulent convection heat transfer over several artificial and realistic rough surfaces is investigated by means of full-geometry-resolving DNS. • The effects of morphological properties (e.g. density and slope) of the roughness as well as its height (in both inner and outer scales) on Reynolds analogy factor are studied independently and systematically. • A correlation for Reynolds analogy factor is extracted using the present DNS results and the data already published in the literature. Abstract Direct Numerical Simulations are used to solve turbulent flow and heat transfer over a variety of rough walls in a channel. The wall geometries are exactly resolved in the simulations. The aim is to understand the effect of roughness morphology and its scaling on the augmentation of heat transfer relative to that of skin friction. A number of realistic rough surface maps obtained from the scanning of gas turbine blades and internal combustion engines as well as several artificially generated rough surfaces are examined. In the first part of the paper, effects of statistical surface properties, namely surface slope and roughness density, at constant roughness height are systematically investigated, and it is shown that Reynolds analogy factor (two times Stanton number divided by skin friction coefficient) varies meaningfully but moderately with the surface parameters except for the case with extremely low slope or density where the Reynolds analogy factor grows significantly and tends to that of a smooth wall. In the second part of the paper, the roughness height is varied (independently in both inner and outer units) while the geometrical similarity is maintained. Considering all the simulated cases, it is concluded that Reynolds analogy factor correlates fairly well with the equivalent sand roughness scaled in inner units and asymptotically tends to a plateau. [ABSTRACT FROM AUTHOR]

Graphical abstract Ammonia-sodium thiocyanate (NH 3 -NaSCN) solution is an efficient absorption refrigeration working pair. This paper presents the experimental results of the flow boiling heat transfer coefficients (HTCs) of ammonia-sodium thiocyanate solution in horizontal tubes with three different inner diameters (4, 6 and 8 mm). The experimental measurements were conducted with different heat fluxes, mass fluxes, boiling temperatures and vapor qualities to determine their impact on the flow boiling heat transfer coefficients of the tested working pair. The results of the present paper will help to extend the basic researching results of NH 3 -NaSCN working pair, which will be of great significance to the design and optimization of generators in the absorption refrigeration units, and will ultimately be beneficial to energy saving for refrigeration area. Highlights • The experimental results of Flow boiling HTCs of NH 3 -NaSCN solution are presented. • Influences of the heat flux, mass flux, exit vapor quality, inner diameter are analyzed. • Comparisons between the NH 3 -NaSCN and NH 3 -LiNO 3 are conducted. • Prediction validation of an updated correlation are provided. • The present results are significant for solar and exhaust heat driven generators. Abstract Ammonia-sodium thiocyanate (NH 3 -NaSCN) solution is an efficient absorption refrigeration working pair. This paper presents the experimental results of the flow boiling heat transfer coefficients (HTCs) of ammonia-sodium thiocyanate solution in horizontal tubes with three different inner diameters (4, 6 and 8 mm). The experimental measurements were conducted with different heat fluxes, mass fluxes, boiling temperatures and vapor qualities to determine their impact on the flow boiling heat transfer coefficients of the tested working pair. Comparison of NH 3 -NaSCN and NH 3 -LiNO 3 solution has been performed and higher flow boiling HTCs of NH 3 -NaSCN solution than that of NH 3 -LiNO 3 solution are obtained under similar working conditions. The results also indicate that the heat flux has the most significant influence on the flow boiling HTCs. The updated flow boiling HTCs prediction correlation for ammonia-salt working pair proposed by the authors’ previous work has been validated by the present experimental results. The present study help to extend the basic researching results of NH 3 -NaSCN working pair, which is of great significance to the design and optimization of generators in the absorption refrigeration units, and will ultimately be beneficial to energy saving for refrigeration area. [ABSTRACT FROM AUTHOR]

Highlights • Gas injection increases a maximal velocity within the boundary layer. • Gas suction decreases an overshoot phenomenon till his absolute vanishing. • Separation of thermal boundary layer at gas injection. • Conventional point of thermal boundary layer separation is proposed. Abstract The paper reports on numerical simulation of an accelerated laminar boundary layer over the heated permeable wall. Xenon was used as a working gas. The finite difference method was applied to solve the system of laminar boundary layer equations with variable gas properties. It has been shown that the gas injection in the mentioned conditions increases maximal velocity within the boundary layer and decreases skin-friction coefficient and thermal Stanton number. On the contrary, gas suction decreases an overshoot phenomenon down to its vanishing and increases skin-friction coefficient and thermal Stanton number up to the value at asymptotic suction. The paper considers a separation of thermal boundary layer at gas injection. The proposed conventional point of this separation is situated farther from the leading edge than the separation point of the dynamic boundary layer. [ABSTRACT FROM AUTHOR]

Highlights • CFD analyses of heat transfer to supercritical carbon dioxide are performed. • AHFM is adopted as a valid tool for modelling the turbulent heat flux. • Phenomena reported in the experimental data were reproduced by the adopted model. Abstract The present paper reports the results of analyses concerning heat transfer to supercritical pressure fluids, performed by adopting a k − ε turbulence model modified in association with the AHFM model, here used for calculating both buoyancy effects and the turbulent heat flux. The promising capabilities of this approach were already highlighted in past studies and the present paper represents a further step in this line of research. Experimental data concerning supercritical carbon dioxide flowing in tubes are here considered, with operating conditions involving both high and low mass flux values and spanning from relatively low inlet temperatures to values higher than the pseudocritical threshold. Some of the interesting features appearing in the experimental data are correctly reproduced by the model, which manages to predict reliable wall temperature trends, both qualitatively and quantitatively. The performed analyses, though reporting successes in a sufficiently wide range of operating conditions, suggest that some parameters of the proposed model should be varied in accordance with the boundary conditions, e.g. considering the mass flux, in order to improve predictions in the most challenging situations. [ABSTRACT FROM AUTHOR]

Highlights • A steam injection physical simulation experimental device is designed independently. • Wellbore stock liquid is taken into consideration to study variable mass flow in horizontal wells. • A novel model is proposed to predict the thermophysical properties of steam along the horizontal wellbore. • The influence of parameters on the thermophysical properties distribution of steam is analyzed in detail. Abstract In this paper, a novel model is proposed to study the variable mass flow process in a long horizontal well and to predict the distribution of thermophysical properties along the horizontal wellbore. First, a physical simulation device is designed to carry out a steam-injection experiment. Second, considering both the uneven distribution of steam and the effect of wellbore stock liquid, a steam-absorption model and a pressure drop model are proposed to predict the distribution of steam and pressure along the horizontal wellbore. Third, the effects of different parameters on steam distribution are analyzed in detail. The results indicate that: (1) an uneven employment phenomenon exists along the wellbore; (2) the pressure distribution along the wellbore is higher when wellbore stock liquid is taken into consideration; (3) the length of the unemployed section along the wellbore increases with an increase in the wellbore stock liquid viscosity; (4) the length of the unemployed section along the wellbore decreases with an increase in the steam injection rate; and (5) the reasonable wellbore length can prevent an unemployed section along the wellbore. This paper presents a basic reference for engineering for parameter optimization as well as for the prediction of steam distribution along the horizontal wellbore. [ABSTRACT FROM AUTHOR]

Graphical abstract Highlights • A numerical modeling framework for thermo-fluids is applied to investigate liquid hydrogen. • Dynamic evolution of liquid hydrogen cavitating flow is numerically investigated. • Thermal transition process and transition temperature in liquid hydrogen are presented. • A modified C -factor is proposed to evaluate and predict thermal transition in thermos-fluids. Abstract The objectives of this paper are to (1) validate the present numerical modeling framework for liquid hydrogen cavitating flow, (2) investigate the dynamic evolution of liquid hydrogen cavitating flows in a wide range of free-stream conditions and (3) propose a thermal parameter to evaluate and predict the transition process of two typical cavitation dynamics in thermo-fluids. The dynamic evolutions of liquid hydrogen cavitating flows in a wide range of free-stream conditions (T ∞ = 14–33 K, U ∞ = 63.9 m/s, 90 m/s, σ ∞ = 0.38, 0.8, 1.2) are numerically investigated. General agreements are obtained between the numerical results and the experimental measurements, including the pressure distribution, temperature drop as well as the cavity structures. The results show that the cavitation behaviors near the triple point are similar to cavitating flows without thermodynamic effects. Two typical cavitation dynamics named the quasi-isothermal mode and the thermo-sensitive mode in liquid hydrogen are observed under the same cavitation number and flow velocity. When the free-stream T ∞ is below 30 K (T N ≤ 0.85), as the temperature increases, the cavity area increases to the maximum around the thermal transition temperature and then decreases when the cavitation dynamics transits from the quasi-isothermal mode to the thermo-sensitive mode. Further analysis indicates that the thermal transition temperature for liquid hydrogen is approximately at T ∞ = 17 K and the transition temperature slightly increases with the increasing free-stream velocity and cavitation number. When the free-stream T ∞ is above 30 K (T N > 0.85), the thermodynamic effects significantly decrease due to the suddenly changed physical properties. And then, the cavitation dynamics turns from the thermo-sensitive mode to the quasi-isothermal mode with the increasing temperature. The modified C -factor including both the thermodynamic effects and the effects of turbulent pressure fluctuations proposed in this paper could quantitatively evaluate and predict dynamics transition from the quasi-isothermal mode to the thermo-sensitive mode in thermos-fluids cavitating flows. It could be utilized as a parameter to design the operating conditions for thermos-fluids apparatus or system to avoid the maximum cavitation aggressiveness around the thermal transition temperature. [ABSTRACT FROM AUTHOR]

Highlights • Effect of external magnetic field on convection heat transfer of ferrofluid flow is investigated. • Magnetic field changes the streamline patterns and enhances the heat transfer. • Increasing magnetic field strength would not necessarily improve the heat transfer. • Magnetic field strength should be large enough to create large circulation zones inside the duct. • Adding line dipoles on high-temperature regions closer to the duct inlet is beneficial. Abstract In this paper, a numerical method is developed to simulate the effect of an external magnetic field on convection heat transfer of a ferrofluid flow inside a rectangular duct. This model is two-dimensional, steady-state, laminar and incompressible. Simple finite volume method is used to couple and solve continuity, momentum and energy equations. In the first part of this paper, it is observed that ferrofluid particles movement along the magnetic field of a line dipole changes the streamline patterns and enhances heat transfer. In the second part of this paper, a sensitivity analysis is performed to investigate the effects of Reynolds number, magnetic field strength, location, and number of line dipoles on the convection heat transfer. It is observed that cooling rate is higher at larger Reynolds numbers while increasing magnetic field strength would not necessarily result in noticeable improvement in heat transfer. It is also noticed that magnetic field strength should be large enough to merge small vortexes and to create large circulation zones inside the duct. Our simulation also revealed that heat transfer could be augmented by adding line dipoles on high-temperature regions closer to the duct inlet. [ABSTRACT FROM AUTHOR]