Computational fluid dynamics (CFD) models were employed to investigate flow conditions inside a model reactor in which yield stress non-Newtonian liquid is mobilized using submerged recirculating jets. The simulation results agree well with the experimental results of active volume in the reactor obtained using flow visualization by the authors in a previous study. The models developed are capable of predicting a critical jet velocity (v(c)) that determines the extent of active volume obtained due to jet mixing. The v(c) values are influenced both by the rheological properties of the liquid and the nozzle orientation. The liquid with higher effective viscosity leads to higher v(c) for a downward facing injection nozzle. However, an upward facing injection nozzle along with a downward facing suction nozzle generates enhanced complementary flow fields which overcome the rheological constraints of the liquid and lead to lower v(c).