We discuss the design and implementation of a preliminary software package written in Python 3 that is intended to represent complex multiphase reactors as networks of ideal continuous stirred tank reactors. This software also implements statistical design of experiments, uncertainty quantification, and global sensitivity analysis. These advanced features can provide important qualitative and quantitative insights into the effect of operating conditions and model parameters on predicted reactor performance. We demonstrate the utility of the program by modeling the vapor phase catalytic upgrading of bio-oil in a bubbling fluidized bed reactor.