Three-dimensional pore-scale simulations of dynamic wicking processes on micro-structured wicks

Capillary wicking characteristics play an important role in two-phase thermal management devices including heat pipes and vapor chambers, yet three-dimensional (3D) pore-scale simulations of the dynamic capillary wicking process on various micro-structured surfaces have been rare. In this paper, we conduct 3D pore-scale simulations of capillary wicking on three commonly used micro-structured wicks including micro-pillar array, micro-channel, and sintered particles. The micro-scale liquid propagation dynamics and the “stick-slip” behavior of the propagating liquid front are captured using a 3D pseudo-potential multiple-relaxation-time lattice Boltzmann method. Based on the Lucus–Washburn approach and a work-energy approach, we theoretically analyze wickabilities of different micro-structured wicks. Effects of wick geometry and structural parameters on the capillary wicking characteristics are discussed. We demonstrate that an optimal pillar pitch distance exists, which maximizes the wickability of the micro-pillar array. We show that when the porosity is relatively low, the wickability of the micro-channel is higher than that of the micro-pillar array and the sintered particles. When the porosity is large, however, the sintered particles exhibit higher wickability than the micro-pillar array and the micro-channel. We also demonstrate that the capillary pressure of the sintered particles is always higher than that of the micro-pillar array and the micro-channel throughout the porosity range investigated. The numerical simulation results are compared with theoretical predictions. Findings in this work provide guidelines for the designs of porous wick in various two-phase thermal management systems for high heat flux devices.