Supersonic impinging jets are a versatile configuration that can model the compressible flows of cold-spray manufacturing and vertical take-off-and landing strategy. In this work, rhoCentralFoam, solver of the OpenFOAM framework, and a large-eddy simulation formulation were used to simulate an underexpanded supersonic jet of Mach 1.45 and nozzle pressure ratio of 4, impinging on a flat wall situated at 1.5 nozzle diameters away from the jet outlet. Care was taken in the mesh construction to properly capture the characteristic standoff shock and vortical structures. The grid convergence index was evaluated with three meshes of increasing spatial resolution. All meshes can generally be considered as sufficient in terms of results focused on time-averaged values and mean physical properties such as centerline Mach number profile. However, the highest resolution mesh was found to capture fine shear vortical structures and behaviors that are absent in the coarser cases. Therefore, the notion of adequate grid convergence may differ between analyses of time-averaged and transient information, and so should be determined by the user’s intention for conducting the simulations. To guide the selection of mesh resolution, scaling analyses were performed, for which the current rhoCentralFoam solver displays a good weak scaling performance and maintains a linear strong scaling up to 4096 cores (32 nodes) for an approximately 40 million-cell mesh. Due to the internode communication bottlenecks of OpenFOAM and improvements in central processing units, this work recommends, for future scaling analyses, adopting a “cells-per-node” basis over the conventional “cells-per-core” basis, with particular attention to the interconnect speed and architecture used.