Fouling of heat-transfer surfaces by fireside deposits can be of significant concern in industrial boilers burning poor-quality fuel. It is commonly controlled by sootblowers that blast deposits with high-pressure supersonic steam or air jets. However, sootblowing is expensive, which motivates efforts to fundamentally understand how sootblower jets behave and how they interact with heat-exchanger geometries and deposits, to guide efforts to improve and optimize sootblower use. Here, we report on the development of a computational fluid dynamics (CFD) model to predict the flow behavior of sootblower jets, work that began with the customization of a research code but has more recently led to an implementation using the commercial CFD software ANSYS Fluent, which makes the model more accessible to the wider engineering community. CFD model results are compared to experimental data that we obtained for jet flow within model tube bank geometries, which are representative of superheaters and generating banks in industrial boilers. The results quantify the deposit removal effectiveness of sootblower jets in the different geometries: the centerline rate of decay of the jet peak impact pressure as a function of the relative position of the sootblower nozzle and tube geometry, the strength of the secondary jets that form when a sootblower jet deflects off of a tube, and the force imposed on various tube positions in the different configurations.