Numerical investigation on transpiration cooling with phase change based on multiphysics coupled finite element method

Transpiration cooling with phase change is an efficient active thermal protection method in aerospace. However, conducting numerical simulation of transpiration cooling with phase change is still a challenge, because porous flow, heat transfer and phase change are coupled with each other, which lead to numerical difficulties. In this article, the previous model is modified and the Zero Thickness Area Model (ZTAM) is proposed, which improves the efficiency and convergence significantly. In addition, the local thermal nonequilibrium effect and latent heat can be considered. To verify this ZTAM model, one-dimensional control equations of mass, momentum, and energy transpiration are solved using finite element method. Furthermore, the position of phase change is determined by fixed point iteration numerical algorithm, which shows good convergence. Finally, the influence on cooling effect by heat flux, coolant mass flow rate, porosity and pore diameter is investigated numerically.