The performance of an axisymmetric multi-fluid Hall thruster code that incorporates a self-consistent, data-driven closure model for the anomalous electron transport is investigated. Five different operating conditions of the H9 magnetically shielded Hall thruster are simulated with the Jet Propulsion Laboratory’s Hall2De. In order to capture the inherent uncertainty associated with the closure model, O(100) simulations are run for each condition, each of using a coefficient set sampled randomly from a probability distribution. The results of these simulations provide probabilistic predictions of thruster performance quantities including thrust, and discharge current, as well as several component efficiencies and centerline plasma properties. The model is found to yield converged solutions at all conditions, with large 10 kHzrange oscillations and performance trends with voltage and flow rate similar to experiment. The model under-predicts the thrust by 15-25% and over-predicts the discharge current by 20% on average compared to experiments at the same discharge voltage and mass flow rate. This performance discrepancy is due to lower beam utilization, mass utilization, and divergence efficiency than experiment, resulting from high Hall parameters in the acceleration region, which lead to a protracted ion acceleration region. The physical processes underlying this result are discussed in the context of future data-driven modeling efforts.