The performance of catalytic membrane reactors (CMRs) depends on the specific details of interactions at different levels between catalytic and separation parts. A clear understanding of decisive factors affecting their operational parameters can be provided via mathematical simulations. In the present paper, main results of numerical studies of ethanol steam reforming, followed by downstream hydrogen permeation through an asymmetric supported membrane, are reported. The membrane module consists of a thin selective layer supported on a substrate with graded porous structure. One-dimensional isothermal reaction–transport model for the CMR has been developed, and its validation has been carried out by using performance data from a lab-scale reactor with a disk-shaped membrane. Simulations demonstrate the model’s capabilities to analyze local concentrations gradients, as required to provide accurate estimates of the relationship between structure–property–performance. It was shown that transport properties of multilayer asymmetric membranes are highly related to the structural properties of each single layer.

Membrane reactors (MR) with an appropriate catalyst are considered to be an innovative and intensified technology for converting a fuel into the hydrogen-rich gas with the simultaneous recovery of high-quality hydrogen. Characteristics of an asymmetric membrane disk module consisting of a gas-tight nanocomposite functional coating (Ni + Cu/Nd_5.5WO_11.25-δ mixed proton-electron conducting nanocomposite) deposited on a gas-permeable functionally graded substrate has previously been extensively studied at lab-scale using MRs, containing the catalyst in a packed bed and in the form of a monolith. The catalytic monolith consisted of a FeCrAl substrate with a washcoat and an Ni + Ru/Pr_0.35Ce_0.35Zr_0.35O₂ active component. It has been shown that the driving potential for hydrogen permeation across the same membrane in a monolithic catalyst –assisted MR is greater compared to the packed bed catalyst. This paper presents results of the study where a one-dimensional isothermal model was used to interrelate catalytic and permeation phenomena in a MR with ethanol steam reforming over the monolith, operating at atmospheric pressure and in the temperature range of 700–900 °C. The developed mathematical reaction–transport model for the constituent layers of the catalyst-asymmetric membrane assembly together with a Sieverts’ equation for the functional dense layer, taking also into account the effect of boundary layers, was implemented in a COMSOL Multiphysics environment. Good agreement with the experimental data of the lab-scale MR with reasonable parameters values is provided. In numerical experiments, concentration profiles along the reactor axis were obtained, showing the effect of the emerging concentration gradient in the boundary layer adjacent to the membrane. Studies have shown that a MR with a catalytic monolith along with appropriate organization of a stagnant feed flow between the monolith and the membrane surface may enhance production and flux of hydrogen, as well as the efficiency characteristics of the reactor compared to a reactor with packed beds.