Passive autocatalytic recombiners (PAR) is the only used method for hydrogen removal from the containment buildings in modern nuclear reactors. Numerical models of such devices, based on the CFD approach, are the subject of this paper. The models may be coupled with two types of computer codes: the lumped parameter codes, and the computational fluid dynamics codes. This work deals with 2D numerical model of PAR and its validation. Gaseous hydrogen may be generated in water nuclear reactor systems in a course of a severe accident with core overheating. Therefore, a risk of its uncontrolled combustion appears which may be destructive to the containment structure.