Methodology and software for pore-scale modeling of multi-phase heterogeneous catalytic reactions is described. Its main features include the modeling of internal wetting with simultaneous phase transition and chemical reaction, and the calculation of liquid- and vapor-phase distribution inside the pore space of a computer-reconstructed catalyst pellet or washcoat layer, under periodically varying external conditions.