Calcium oxide/Calcium hydroxide can be utilized as a reaction system for thermochemical heat storage. It features a high storage capacity, is cheap, and does not involve major environmental concerns. Operationally, different fixed-bed reactor concepts can be distinguished; direct reactor are characterized by gas flow through the reactive bulk material, while in indirect reactors, the heat-carrying gas flow is separated from the bulk material. This study puts a focus on the indirectly operated fixed-bed reactor setup. The fluxes of the reaction fluid and the heat-carrying flow are decoupled in order to overcome limitations due to heat conduction in the reactive bulk material. The fixed bed represents a porous medium where Darcy-type flow conditions can be assumed. Here, a numerical model for such a reactor concept is presented, which has been implemented in the software DuMux. An attempt to calibrate and validate it with experimental results from the literature is discussed in detail. This allows for the identification of a deficient insulation of the experimental setup. Accordingly, heat-loss mechanisms are included in the model. However, it can be shown that heat losses alone are not sufficient to explain the experimental results. It is evident that another effect plays a role here. Using Bayesian inference, this effect is identified as the reaction rate decreasing with progressing conversion of reactive material. The calibrated model reveals that more heat is lost over the reactor surface than transported in the heat transfer channel, which causes a considerable speed-up of the discharge reaction. An observed deceleration of the reaction rate at progressed conversion is attributed to the presence of agglomerates of the bulk material in the fixed bed. This retardation is represented phenomenologically by mofifying the reaction kinetics. After the calibration, the model is validated with a second set of experimental results. To speed up the calculations for the calibration, the numerical model is replaced by a surrogate model based on Polynomial Chaos Expansion and Principal Component Analysis.