Thin-film evaporation in wick structures for cooling high-performance electronic devices is attractive because it harnesses the latent heat of vaporization and does not require external pumping. However, optimizing the wick structures to increase the dry-out heat flux is challenging due to the complexities in modeling the liquid vapor interface and the flow through the wick structures. In this work, we developed a model for thin-film evaporation from micropillar array wick structures and validated the model with experiments. The model numerically simulates liquid velocity, pressure, and meniscus curvature along the wicking direction by conservation of mass, momentum, and energy based on a finite volume approach. Specifically, the three-dimensional meniscus shape, which varies along the wicking direction with the local liquid pressure, is accurately captured by a force balance using the Young Laplace equation. The dry-out condition is determined when the minimum contact angle on the pillar surface reaches the receding contact angle as the applied heat flux increases. With this model, we predict the dry-out heat flux on various micropillar structure geometries (diameter, pitch, and height) in the length scale range of 1-100 mu m and discuss the optimal geometries to maximize the dry-out heat flux. We also performed detailed experiments to validate the model predictions, which show good agreement. This work provides insights into the role of surface structures in thin-film evaporation and offers important design guidelines for enhanced thermal management of high-performance electronic devices.