Water-gas shift membrane reactors (WGS-MRs) offer a pathway to affordable blue H2 generation/purification from gasified feedstock or reformed fuels. To exploit their cost benefits for blue hydrogen production, WGS-MRs’ performance needs to be optimized, which includes navigating the multidimensional design space (e.g., temperature, feed pressures, space velocity, membrane permeance and selectivity, catalytic performance). This work describes an equation-oriented modeling framework for WGS-MRs in the Pyomo ecosystem, with an emphasis on model scaling and multi-start initialization strategies to facilitate reliable convergence with nonlinear optimization solvers. We demonstrate, through sensitivity analysis, that our model converges rapidly (< 1 CPU second on a laptop computer) under a wide range of operating parameters (e.g., feed pressures of 1-3 MPa, reactor temperatures of 624-824 K, sweep-to-feed ratios of 0-0.5, and steam/carbon ratios of 1-5). Ongoing work includes (1) validation and calibration of the WGS-MR model using benchtop laboratory data and (2) design, intensification, and optimization of blue H2 processes using the WGS-MR model.