The design of the water-cooled ceramic breeder (WCCB) blanket includes beryllium multiplier layers. At elevated temperature, beryllium reacts with steam in an exothermic reaction producing beryllium oxide and hydrogen. Such situation may occur in WCCB in the case of the rupture of one of the cooling pipes in the blanket module. This process occurs locally in a complex 3-D geometry of the blanket containing several different granular levels and a network of cooling pipes and structural supports. The process is also inherently unsteady since reaction rate depends on concentration of steam and pure beryllium which changes in time. In order to perform the detailed analysis of the process, the model of the reacting flowthrough porous media was developed and introduced into 3-D computational fluid dynamics code. In this model, granular beds are introduced as porous solids simplifying the model geometry and reducing typical mesh size to manageable amount of tens of millions of elements. A reaction rate between solid beryllium and steam is obtained from experimental results, and depends on temperature and concentration of the reactants. The differential equation for beryllium oxide fraction is introduced, allowing obtaining distributions of beryllium oxide in space and time. Multicomponent flow consisting of a homogenous mixture of steam and hydrogen is considered flowing through the porous solid with variable properties. Sink and source terms for steam and hydrogen fractions are determined by local beryllium oxide mass fraction source according to the molar ratios of beryllium steam reaction. The conjugated heat transfer approach is applied to calculate heat transfer in support structures as well as coolant flow, simultaneously with the porous medium steam flow in a blanket’s granular beds. The model is validated using experimental data on beryllium steam reaction for granular bed samples.