Study on heat transfer performance of cold plate with grid channel | Scientific Reports

Scientific Reports volume 14, Article number: 4913 (2024) Cite this article

1630 Accesses

Metrics details

The utilization of cold plate radiators as a prevalent method for indirect liquid cooling has been extensively investigated and implemented in server cooling systems. However, there is a lack of comprehensive study on the application of this technology at the chip size, indicating a need for more development and exploration in this area. A proposal was made for a grid-channel chip cold plate heat sink to facilitate the dissipation of heat from a chip. tests were conducted to investigate the impact of the flow rate of the cold plate and the layout of the inlet and outlet on various thermal parameters, including the average temperature, maximum temperature, thermal resistance, and uniformity coefficient of the cold plate. The tests were specifically conducted under a chip power of 150W, and the accuracy of the simulation was confirmed through the use of FLUENT. The findings indicate that the cold plate effectively regulates the temperature of the chip, ensuring it remains below 85 °C throughout all experimental groups. In contrast to the single in single out configuration, the single-in multiple-out layout exhibits a higher degree of temperature uniformity within the cold plate. Nevertheless, it is important to note that augmenting the quantity of exits does not guarantee an improvement in heat transfer efficiency. This outcome is contingent upon the presence of a longitudinal flow channel shared by the outlet and intake, as well as the dispersion characteristics of the outlet. Enhancing the dispersion of the exit can significantly enhance the thermal transfer efficiency of the cold plate. Furthermore, a strategy for adjusting the aperture of the orifice is proposed as a solution to address the challenges related to flow uniformity and the issue of high pressure drop in the cold plate.

The introduction of domestically developed chips from China has resulted in a notable enhancement in the performance of these chips. Moreover, the heat flux of these chips has experienced an increase as a result of improvements in technology. The issue of overheating presents a significant challenge in contemporary electronic equipment, namely in the domains of high-performance computing, artificial intelligence, and mobile electronics. The issue of chip overheating has significant implications for the stability and integrity of electronic components and materials, as well as the potential for component failure or damage1,2. As a result, the occurrence of server failure and subsequent loss of data might have detrimental effects on system security3,4. Consequently, the use of efficient temperature management and heat dissipation techniques are imperative in order to mitigate overheating and improve the operational stability and dependability of the microchip5.

The primary approaches for heat dissipation in servers are air cooling and liquid cooling6. Air cooling, while having a simple structure, has limited heat dissipation capacity7. In contrast, liquid cooling employs liquid media and radiators as a means to transport heat to the external environment, resulting in enhanced heat dissipation performance and improved energy efficiency8. Indirect liquid cooling technology, namely cold plate technology, has garnered considerable interest within the realm of liquid cooling technology owing to its straightforward implementation and dependable system performance. The configuration of channels within the cold plate is a crucial factor in determining the effectiveness of heat transfer. The analysis conducted by Huo et al.9 investigated the effects of channel number, flow direction, flow rate, and ambient temperature on a mini-channel cold plate with a parallel flow pattern. In a comparable manner, an analogous investigation was conducted in Reference10 to optimize the effectiveness of heat transfer. The primary differentiation lies in the convergence of the inlet and outlet inside the parallel channels.

The cold plate is a vital component in the field of indirect liquid cooling heat transfer technology, and has attracted considerable attention11,12,13. Previous research has explored the heat dissipation of the cold plate, with a focus on thermal economic analysis14, energy saving potential15, and especially structural optimization12. Porter Harris et al.16 has developed a series of spoiler structures for cold plates used in lithium battery thermal management. These structures effectively regulate the battery's temperature to ensure it remains within the optimal range throughout operation. Alfaryjat et al.17 conducted a comparative experiment on flow channels with varying cross sections as part of their investigation into the traditional study of flow channel topology. It was discovered that at a Reynolds number of 2200, the convective heat transfer coefficient of the circular channel may reach 1500 \({{\text{W}}/{\text{cm}}}^{2}\), while the maximum temperature of the heat source is below 75 °C. However, as compared to other experimental groups, the circular channel exhibits the most pressure drop. Song et al.18 examined the flow channels consisting of six trapezoidal sections. By altering the aspect ratio and flow direction of the channel section, it is determined that the heat transfer efficiency is higher in the countercurrent flow direction compared to the downstream flow direction due to the presence of greater disturbance. Raj et al.19 modified the topology of the flow channel to incorporate more intricate stepped and parabolic shapes. The two new types of flow channels, with their increased number of cross section angles, facilitate the formation of a gasification core compared to the traditional rectangular flow channel. This poses a challenge for the rectangular flow channel to promptly release the boiling gasification fluid under specific superheat conditions. The heat transfer rates of parabolic and stepped flow channels exhibit an increase of 88% and 169% respectively, compared to standard rectangular flow channels. Pulugundla et al.20 employed the time-accurate computational fluid dynamics (CFD) technique to conduct a numerical simulation of the thermal management system (TMS) for lithium batteries. A comparative analysis is conducted on the heat transfer properties of 8 distinct flow channels. Ultimately, the thermal conductivity of different cold plates is comparable, but the pressure loss caused by fluid flow varies significantly. Hence, it is imperative to take into account the decrease in pressure during the fluid flow process while constructing the structure of the flow channel. To address the issue of pressure drop in the cold plate. Wang et al.21 introduced many transverse gaps in the fins. These gaps serve to connect neighboring channels and improve the fluid mixing between them. These investigations have demonstrated that the cold plate has promising heat exchange capabilities and offers advantages compared to air cooled heat exchange systems.

Currently, the variety of flow channel types in cold plates is limited, with a lack of innovative structures in cold plate channels. There is a very limited number of studies examining the impact of inlet/outlet configurations on the heat transfer performance of cold plates. Research is also lacking in the characteristics and spatial arrangement of the number of outlets in cold plates. Additionally, there are few studies on methods to optimize the uniformity of fluid flow within cold plates.

This study investigates the heat dissipation performance of a grid-type-channel cold plate for a chip-scale heat source. It evaluates the influence of different inlet/outlet configurations on the cold plate's heat transfer efficiency and suggests an flow channel optimal design.

Figure 1 illustrates the physical diagram of the cold plate, along with the structural and physical diagrams of the analogue chip. The cold plate under consideration is a cube made of copper, with dimensions of 150 mm × 150 mm × 15 mm. The cold plate is encompassed by nine equidistant through holes.

Structure of cold plate and analog chip.

Based on the trial group’s findings, the through holes of the cold plate have been effectively sealed, and the screw holes have been appropriately positioned according to the specified criteria. Additionally, the inlet and outlet of the cold plate have been connected to the circulating pump using rubber pipes. The analogue chip situated at the center of the cold plate will undergo repair using thermal grease (HY500). The thermocouple (GG-K-30) should be positioned according to the configuration depicted in Fig. 1c. Proceed by activating the power supply, adjusting the regulator to manipulate the pump and heating power of the analogue chip, and initiating the experimental procedure.

The ambient temperature of the cold plate is set to 25 °C. The pump is activated to facilitate water circulation within the cold plate. The flow rate is regulated by adjusting the voltage on the pump using a regulator. The heating chip is employed for the purpose of emulating heat generation, while the regulator is calibrated to maintain a constant power reading of 150W on the heating chip power. The inlet and outlet arrangement of the cold plate should be modified to Y, U, type I, as illustrated in Fig. 2. The experiment must be repeated by adjusting the distance between the two outlets of the Y-arrangement, ranging from 30 to 120 mm.

One inlet and one outlet Settings (a) Type T, (b) Type Y, (c) Type U, (d) Type I.

This section employs SolidWorks for modeling research and provides a description of the cold plate model, governing equations, and boundary conditions.

The schematic representation of the cold plate radiator simulation model is illustrates in Fig. 3. The entire cold plate model is a cubic structure measuring 150 mm × 150 mm × 15 mm. The model consists of four sides, each of which is perforated with a total of nine holes. These holes are evenly spaced at a distance of 15 mm from each other and have a diameter of 6 mm. The chip is located at the central region of the cold plate.

Schematic diagram of the physical model.

For the numerical computations the below assumptions are given here: (1) the fluid(water) is considered to be uncompressible; (2) the flow process is steady-state, constant; (3) disregard the transfer of heat through convection between the system and its surroundings. (4) In the proposed model, the heat source distributes thermal energy homogeneously across the system, and the rate of heat supply remains invariant throughout the process.

The governing differential equation of the model is shown as follows:

Continuity Equation:

momentum equation:

energy equation:

where \(V\) is the volume; \(\rho \) is the density; ▽ Divergence; μ is viscosity; P is pressure; \({C}_{P}\) is the specific heat capacity at constant pressure; T is the temperature; k is thermal conductivity; S is the source entry.

In the present numerical study, the ANSYS CFX software (The copyrights for all commercial software applications employed in this manuscript have been duly acquired.) is used to model the turbulent flow and conjugate heat transfer. The Reynolds number (Re) from the experimental study is used to determine the nature of the flow (laminar or turbulent) based on Eq. (4).

where \(\rho \) denotes the fluid density, \(v\) denotes the mean fluid velocity, \(D\) denotes the hydraulic diameter and \(\mu \) denotes the dynamic fluid viscosity. The Re of inlet18 is \({4.18\times 10}^{4}\). Therefore, turbulent flow is used in this study. The \(k-\varepsilon \) turbulence model is used.

The turbulent kinetic energy and eddy viscosity equations are provided as follows:

where \({C}_{1}{,C}_{2},{C}_{3}\) denote the model constants, \({S}_{M}, {S}_{\varepsilon }\) denote the UDF source terms, \({Y}_{M}\) denotes the contribution of the fluctuating dilatation in the compressible turbulence to the overall dissipation rate, \({G}_{b}\) denotes the generation of turbulence kinetic energy due to buoyancy, \({G}_{k}\) denotes the generation of turbulence kinetic energy due to the mean velocity gradients, and \({\sigma }_{k}\) and \({\varepsilon }_{k}\) denote the turbulent Prandtl numbers for \(k\) and \(\varepsilon \) respectively.

The turbulent viscosity is calculated as follows:

Material and boundary conditions are set as shown in Table 1.

The convergence criteria of temperature, velocity, pressure and turbulent valuables are set as \({10}^{-3}\) in the numerical simulation.

To ensure both calculation accuracy and efficiency, the appropriate number of mesh elements is determined, and a local mesh encryption is applied to the fluid region. For the verification of grid independence, the cold plate model in Fig. 3 was selected, and the fluid region was individually coded with a minimum mesh size of 2.97 × \({10}^{-12}{{\text{m}}}^{3}\).The present study aimed to evaluate the impact of varying mesh numbers on the computed results of the volume-weighted average temperature of the fluid within the cold plate. This investigation was conducted while maintaining a constant inlet flow rate of 43 ml/s. It is observed that thermal equilibrium is achieved when the number of grid nodes reaches 2,145,400 as shown in Fig. 4.

Grid independence verification result.

The pressure–velocity coupling in numerical solution processes is realized by the pressure correction method based on the Coupled algorithm. In the spatial discretization of computational domains, the pressure term is discretized using the PRESTO algorithm. Figure 5 shows the trend of thermal resistance variation with flow rate in experimental and simulation test results. Thermal resistance can be defined as the resistance of the heat source in the system to transfer heat to the circulating liquid. The expression is shown as follows:

where \({T}_{source}\) is the average temperature of the heat source, °C, \({T}_{f,in}\) is the liquid temperature at the inlet of the cold plate, °C, Q is heat source power, W.

Comparison of experimental and simulated thermal resistance with flow rate.

The trend observed in the results indicates a consistent relationship between the simulation and experimental outcomes. Specifically, the simulation calculations neglect the convection heat transfer between the cold plate and the external environment, leading to a higher surface temperature of the chip compared to the experimental group. Consequently, this discrepancy results in an overall increase in thermal resistance. The maximum difference observed between the two approaches is only 7.65%. Therefore, the simulation is capable of providing a realistic depiction of the flow and heat transfer properties in this particular scenario.

In this paper, the heat transfer performance of the cold plate with different flow rates is compared by simulation.

To ensure optimal operating performance, the maximum temperature of the chip is maintained below 85 °C. It is because that the higher temperatures have been observed to negatively impact its functionality4. Figure 6 illustrates the relationship between the temperature of the chip and the fluid flow within the cold plate. The results indicate that as the flow rate increases, the temperature of the chip decreases. Specifically, when the flow rate increased from 1 to 2 m/s, the chip temperature exhibited a decrease by 3.89%. However, when the flow rate further increased to 7 m/s, the chip temperature only decreased by 0.61%. Based on these findings, a flow rate of 7 m/s is selected for experimental purposes in this study.

Variation of temperature of heat source with flow rate.

The findings from the experiment align closely with the simulated results reported in previous scholarly works22,23, which investigate the impact of chip temperature on the flow rate of cold plates. The reason for this phenomenon is attributed to the significant flow resistance within the grid channel, which results in an increase in flow rate. Consequently, the pressure drop within the cold plate becomes excessively high, hindering the circulation of the fluid.

A comparison between SIMO (T and Y) and SISO (U and I) was made for the proposed grid-type-channel cold plate, and the liquid flow in the cold plate was studied by FLUENT software.

Various configurations of the inlet/outlet are simulated, resulting in the generation of a temperature and velocity cloud image (Figs. 7, 8). The uniformity of water temperature in the cold plate is enhanced while employing T-shaped and Y-shaped arrangements. From the velocity cloud image, it can be seen that the T-shaped and Y-shaped inlet and arrangements can better transport fluid to the cold plate.

Temperature cloud image of cold plate section.

Velocity cloud image of cold plate section.

The mean of uniformity coefficient is one of the indexes to evaluate the uniformity of the surface temperature of an object24. A higher value indicates a poorer mean temperature on the covered side. The calculation method is shown as follows:

where N is the number of observation point; \({T}_{i}\) is the temperature of each point; \({T}_{ave}\) is the average temperature of each point.

This study assesses the four essential performance indicators (KPIs) listed in Table 2, namely chip temperature, cold plate temperature, cold plate uniformity coefficient, and thermal resistance. The objective is to analyze the impact of inlet/outlet configuration on the heat.

As can be seen in Table 2, SIMO system has obvious advantages in all parameters compared with SISO system, and the average temperature coefficient of cold plate has a maximum difference of 73.83%. The maximum thermal resistance difference is 48.15%.

This section will change the outlet quantity of cold plate and the setting mode of outlet. The specific setting mode is shown in Table 3 (this study only focuses on the axisymmetric outlet arrangement mode). The objective is to investigate the influence of these modifications on the heat transfer performance of cold plate. Standard deviations are introduced to assess the dispersion of outlets in various Settings. Specific methods are shown as follows.

Standard deviation:

where d is the distance between the outlet and the symmetry axis of the cold plate which is shown in Fig. 9.

Schematic diagram of each outlet number of cold plate.

The influence of varying outlet quantity on heat transfer was assessed by computing the average temperature coefficient of the cold plate and the average temperature of the chip for each outlet quantity group in Table 3, using all setting techniques.

As illustrated in Figs. 10 and 11, there is a discernible trend wherein the average temperature coefficient of the cold plate and the temperature of the chip initially increase and thereafter decrease as the number of outlets is augmented.

Influence of the number of outlets on chip temperature.

Influence of outlet quantity on uniformity coefficient of cold plate.

Figures 10 and 11 demonstrate that the heat transfer efficiency of the cold plate with an odd number of outlets is much lower compared to the cold plate with an even number of outlets.

According to the continuity equation, if the entrance flow is determined, the outlet flow is assured. However, in this investigation, different quantities of exits will have varied configurations. These diverse layouts have a considerable influence on the flow of the fluid in the cold plate. As indicated above, since this work only covers symmetrical exit Settings, when the number of cold plate exits is odd, the cold plate exit No. 0 (Fig. 9) will certainly be opened. In Table 4, it can be shown that in all odd exit Settings, outlet 0 in all cases accounts for more than 66% of the entire Mass flow rate. The outlet and the entry share a longitudinal flow channel, resulting in a large volume of flow cannot flow through the cold plate and out of other outlets, which deteriorates the temperature equalization and heat exchange capacity of the cold plate.

The influence of output dispersion on chip temperature is presented in Fig. 12. The graph represents the single setting method for 8 and 9 outlets through the use of horizontal lines. It can be observed that the chip temperature remains relatively close to 73 °C for each outlet quantity. Furthermore, it is observed that the cooling effect of the cold plate positioned at 2–7 outlets on the chip diminishes to varying extents as the outlet dispersion degree increases. One notable optimization among the options is the implementation of the cold plate set designed for the 2 outlets. It is evident that this configuration, known as the I configuration, has the poorest relative heat transfer performance. However, it is capable of reducing the temperature by a maximum of 4.89%. This observation demonstrates a positive correlation between the degree of dispersion at the outlet of the cold plate and its heat transfer capability.

Influence of outlet dispersion degree on chip temperature.

In Fig. 7, where the T and Y configurations of the 0th flow channel (Fig. 9) exhibit minimal fluid flow on both sides. This is attributed to the uniform aperture size across the grid-style cold plate, preventing effective fluid guidance for uniform flow within the cold plate. We elected to block specific channels within the cold plate to direct the fluid in an S-shape (as shown in Fig. 13), but this resulted in a significant increase in pressure drop. According to the research by Wang et al.21, adding gaps to the fins helps alter the overall pressure drop of the heat exchanger. In subsequent studies, we reduced the diameter of the pipes between the S-shaped flow channels, treating them as gaps, which achieved the desired uniformity in cold plate flow. The temperature variation in each chip remained relatively minor (Table 5). It was observed that when the gap exceeded 3 mm, the effectiveness of flow guidance greatly diminished, with the lowest pressure drop reaching 9632 Pa (Table 5).

Velocity contour map for flow channel optimization.

Based on the analysis conducted, it can be observed that the distribution of water temperature in the cold plate is more homogeneous when the arrangement is either T-shaped or Y-shaped. The preference for the SIMO configuration is attributed to its ability to facilitate a consistent liquid flow within the cold plate. The chip has a comparatively low temperature. The limited proximity of the inlet/outlet to the chip component results in a minimal direct flow of fluid through the chip. Consequently, the circulation of liquid in the U-shaped configuration of the cold plate is more comprehensive compared to the I-shaped configuration. This enhanced circulation contributes to the lower temperature of the liquid in the U-shaped cold plate. In the SISO type design, the majority of the liquid exhibits challenges in effectively engaging in the circulation process to facilitate heat dissipation. Consequently, this leads to the accumulation of heat in the bottom region of the chip, exacerbating the overall heat dissipation issue.

The disparity observed in Table 2 can be attributed to the distinct arrangement of multiple outlets in the SIMO system. This arrangement allows for the fluid to be directed to various sections of the cold plate, enabling multiple outlets to share the burden of discharging heat. In contrast, the SISO system relies on a single outlet to handle the entire load of discharged heat. Therefore, the SISO system faces challenges in effectively dissipating heat in the cold plate, leading to increased temperatures and a deterioration in heat transfer efficiency.

Based on the analysis of Figs. 10 and 11, it is observed that when the outlet reaches a height of 4–6 units, the primary flow of the mainstream fluid is concentrated near the outlet and in close proximity to the symmetry axis. This concentration of flow restricts the diffusion of heat throughout the entire cold plate. As the number of outlets increases, the nearest outlet to the symmetry axis moves further away. This means that the flow of the mainstream fluid becomes less concentrated near the symmetry axis, allowing for better distribution of heat throughout the cold plate. The increased dispersion of outlets helps to enhance the overall heat transfer capability of the cold plate. It is important to note that the inflow rate remains constant throughout this analysis, indicating that the changes in heat distribution and heat transfer efficiency are primarily influenced by the arrangement and dispersion of the outlets on the cold plate.

As previously mentioned, the diffusion of heat within the entire cold plate is impeded. As the number of outlets continues to increase, the proximity of the nearest outlet to the symmetry axis decreases. Consequently, the pressure drop experienced during outflow decreases, causing the flow near the symmetry axis to become saturated. This saturation forces the fluid to exit through outlets located at a greater distance from the symmetry axis. Furthermore, in the case of a cold plate with an odd number of outlets, the outlet that has a longitudinal flow route with the entry will be opened, resulting in a number of water was discharged from outlet 0. As a result, a significant quantity of fluid (more than 66%) was observed to outlet through outlet 0, thereby exacerbating the heat transfer efficiency of the cold plate. Additionally, it can be inferred from this observation that the presence of an odd number of outlets in a cold plate has a comparatively lesser impact on temperature uniformity and chip temperature when compared to cold plates with an even number of outlets.

To improve the uniformity of fluid flow within the grid cold plate, we reduced the diameter of certain flow channels as shown in Fig. 13, effectively utilizing these narrowed sections as gaps to guide the direction of fluid flow within the cold plate. This approach, under conditions where the temperature differential of the heat source was not significant, effectively resolved the issue of excessive pressure drop in the S-shaped flow channels. It reduced the pressure drop in the cold plate to as low as 9632 Pa and significantly enhanced the uniformity of flow and temperature distribution within the cold plate.

Through the evaluation of the heat transfer performance of a grid-type-channel cold plate chip heat exchanger, the following results have been derived:

The grid-type-channel cold plate developed in this study underwent experimental testing and fluid simulation. The simulation results validated the accuracy of the simulation by demonstrating a maximum residual error of only 7.65%.

As the flow rate increases, the rate at which the chip temperature decreases is reduced. Specifically, when the flow rate is increased from 7 to 8 m/s, the chip temperature experiences a decrease of only 0.61%.

A comparison was made between the heat exchange performance of SIMO (Single Inlet Multiple Outlet) and SISO (Single Inlet Single Outlet) cold plates. Due to the inability of the SISO configuration to uniformly distribute high temperature fluid across the whole cold plate, the heat exchange efficiency was significantly lower compared to that of the SIMO system. Consequently, the average temperature coefficient of the cold plate exhibited a maximum deviation of 73.83%. The disparity in heat resistance reaches a maximum of 48.15%.

The rise in the quantity of outlets has the potential to negatively impact the heat exchange efficiency of the cold plate within the range of 4–6 outlets. Conversely, as the appropriate number of outlets continues to increase within the range of 6–9, the heat exchange performance of the cold plate is expected to improve. The heat transfer performance of the cold plate is compromised when the inlet and outlet are aligned along the same longitudinal flow path. This limitation arises from the predominant flow of fluid down this single path, preventing effective distribution across the entire cold plate.

The concept of outlet dispersion was introduced, and a comparison was made about the variations in chip temperature under different conditions of outlet quantities and designs. The study revealed that as the dispersion of the outlet increased, the heat emitted by the chip heat source dropped to varying extents in all cases. Among the various configurations, it was seen that the cold plate set designed for outlet 2 exhibited the most significant optimization. Furthermore, the type I arrangement, which displayed the poorest relative heat transfer performance, achieved a temperature reduction of up to 4.89%. This observation demonstrates a positive correlation between the degree of dispersion at the outlet of the cold plate and its heat transfer capability.

By altering the diameters of flow channels in different parts of the cold plate to guide fluid flow, and employing small-diameter gaps between these channels, we can effectively reduce the flow pressure drop in the cold plate. This approach also significantly enhances the uniformity of fluid flow and temperature distribution within the cold plate.

In summary, this study assessed the heat transfer performance of the grid cold plate and determined the impact of flow rate, arrangement mode, inlet/outlet configurations on chip heat dissipation and optimization of the flow channel. Additionally, optimal arrangement conditions were proposed to achieve optimal heat dissipation efficiency. The cold plate under consideration demonstrates efficient heat transfer capabilities to chip-scale heat sources, hence offering a promising approach for addressing semiconductor heat dissipation challenges.

The datasets used and/or analysed during the current study available from the corresponding author on reasonable request.

Sharma, C. S., Zimmermann, S., Tiwari, M. K., Michel, B. & Poulikakos, D. Optimal thermal operation of liquid-cooled electronic chips. Int. J. Heat Mass Transf. 55, 1957–1969 (2012).

Article CAS Google Scholar

Yeh, L. T. et al. Thermal Management of Microelectronic Equipment: Heat Transfer Theory Analysis Methods and Design Practices 144 (ASME International, 2002).

Book Google Scholar

Amon, C. H. et al. MEMS-enabled thermal management of high-heat-flux devices EDIFICE: Embedded droplet impingement for integrated cooling of electronics. Exp. Therm. Fluid Sci. 25, 231–242 (2001).

Article Google Scholar

Chethana, G. D. & Gowda, S. Influence of server density in racks on data centre thermal management: A numerical study. Int. J. Sci. Res. Sci. Eng. Technol. 1, 6 (2022).

Google Scholar

Zhang, Z. Y., Wang, F., Qu, Y., Wang, H. J. & Huang, Y. Y. Analysis of heat dissipation effect and energy consumption on a new type of cooling device used in server CPU. J. Shanghai Univ. Eng. Sci. 36(4), 317–321 (2014).

Google Scholar

Shrigondekar, H., Lin, Y. C. & Wang, C. C. Investigations on performance of single-phase immersion cooling. Int. J. Heat Mass Transf. 206, 123961 (2023).

Article CAS Google Scholar

Cheng, R. Study on heat dissipation techniques of computer CPU chips. J. Liaoning Teach. Coll. 16(4), 90–93 (2014).

Google Scholar

Sharma, R. D., Sharma, D., Awasthi, K. & Shamsi, N. A. Power management and energy optimization in hybrid electric vehicle—A review. Adv. Fluid Therm. Eng. 13, 978–981 (2019).

Google Scholar

Huo, Y. et al. Investigation of power battery thermal management by using mini-channel cold plate. Energy Convers. Manag. 89, 387–395 (2015).

Article Google Scholar

Qian, Z., Li, Y. & Rao, Z. Thermal performance of lithium-ion battery thermal management system by using mini-channel cooling. Energy Convers. Manag. 126, 622–631 (2016).

Article CAS Google Scholar

Qin, Y., Cai, C., Biao, C. & Kun, Z. Experimental research on liquid cooling performance of a domestic dual-socket server. Comput. Eng. Sci. 45(4), 607–612 (2023).

Google Scholar

Akbarzadeh, M. et al. A novel liquid cooling plate concept for thermal management of lithium-ion batteries in electric vehicles. Energy Convers. Manag. 231, 113862 (2021).

Article Google Scholar

Mo, X., Zhi, H., Xiao, Y., Hua, H. & He, L. Topology optimization of cooling plates for battery thermal management. Int. J. Heat Mass Transf. 178, 121612 (2021).

Article Google Scholar

Bao, Y. & Shao, S. Numerical study on ultrathin wide straight flow channel cold plate for Li-ion battery thermal management. J. Energy Storage 64, 107263 (2023).

Article Google Scholar

Dubey, P., Pulugundla, G. & Srouji, A. K. Direct comparison of immersion and cold-plate based cooling for automotive Li-ion battery modules. Energies 14, 1259 (2021).

Article CAS Google Scholar

Porter Harris, W., et al. Cold plate blade for battery modules. U.S. Patent No. 11557800 B2. U.S. Patent and Trademark Office (2023).

Alfaryjat, A. A. et al. Influence of geometrical parameters of hexagonal, circular, and rhombus microchannel heat sinks on the thermohydraulic characteristics. Int. Commun. Heat Mass Transf. 52, 121–131 (2014).

Article Google Scholar

Song, J. et al. Numerical studies on the hydraulic and thermal performances of trapezoidal microchannel heat sink. Int. J. Therm. Sci. 161, 106755 (2021).

Article Google Scholar

Raj, S., Pathak, M. & Khan, M. K. Flow boiling characteristics in different configurations of stepped microchannels. Exp. Therm. Fluid Sci. 119, 110217 (2020).

Article Google Scholar

Pulugundla, G., Dubey, P. & Srouji, A. Time-Accurate CFD Analysis of Liquid Cold Plates for Efficient Thermal Performance of Electric Vehicle Li-Ion Battery Modules. SAE Technical Paper. 2019-01-0500, 8 (2019).

Wang, S. L. et al. Heat transfer enhancement of symmetric and parallel wavy microchannel heat sinks with secondary branch design. Int. J. Therm. Sci. 171, 107229 (2022).

Article Google Scholar

Yang, S. et al. Investigation of Z-type manifold microchannel cooling for ultra-high heat flux dissipation in power electronic devices. Int. J. Heat Mass Transf. 218, 124792 (2024).

Article CAS Google Scholar

Nangle-Smith, S. & Cotton, J. S. Phenomenological correlations for two phase heat transfer and pressure drop performance of electrohydrodynamic horizontal convective boiling of dielectric fluids. Int. J. Heat Mass Transf. 132, 1018–1028 (2019).

Article CAS Google Scholar

Huang, Y., Ge, J., Chen, Y. & Zhang, C. Natural and forced convection heat transfer characteristics of single-phase immersion cooling systems for data centers. Int. J. Heat Mass Transf. 207, 124023 (2023).

Article Google Scholar

Download references

This work was supported by a grant from the Educational Commission of Guangdong Province of China with No. 2021ZDZX1106.

School of Mechanical and Electrical Engineering, Shenzhen Polytechnic University, Shenzhen, Guangdong, China

Kaihua Zhang & Zhijian Hou

School of Optical Information and Energy Engineering, Wuhan Institute of Technology, Wuhan, Hubei, China

Kaihua Zhang, Xiaojiang Ye, Yahui Ren & Qiuyi Shi

School of Energy and Environment, City University of Hong Kong, Hong Kong, China

Qiuyi Shi

You can also search for this author in PubMed Google Scholar

K.Z.: Conduct experiments and data processing, schematic drawing, simulation model testing. X.Y.: Thesis guidance teacher. Z.H.: Provision of test sites and funding. Y.R.: processing of reference1,2,3,4,5,6,7,8,9. Q.S.: processing of reference9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24.

Correspondence to Zhijian Hou.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

Zhang, K., Ye, X., Hou, Z. et al. Study on heat transfer performance of cold plate with grid channel. Sci Rep 14, 4913 (2024). https://doi.org/10.1038/s41598-024-52238-6

Download citation

Received: 30 October 2023

Accepted: 16 January 2024

Published: 28 February 2024

DOI: https://doi.org/10.1038/s41598-024-52238-6

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative