where the constants q j n,vap and q j n,liq for the n th time step are chosen to reduce the error in steady-state solution associated with operator splitting as discussed in Ref. 30 and in our previous work. Note that the transport equations are decoupled in species, so the transport of each species can be solved separately. Also, the reaction equation is decoupled by trays.