A reduced kinetic method (RKM) with a first-principles collision operator is introduced in a 1D2V planar geometry and implemented in a computationally inexpensive code to investigate non-local ion heat transport in multi-species plasmas. The RKM successfully reproduces local results for multi-species ion systems and the important features expected to arise due to non-local effects on the heat flux are captured. In addition to this, novel features associated with multi-species, as opposed to single species, cases are found. Effects of non-locality on the heat flux are investigated in mass and charge symmetric and asymmetric ion mixtures with temperature, pressure, and concentration gradients. In particular, the enthalpy flux associated with diffusion is found to be insensitive to sharp pressure and concentration gradients, increasing its significance in comparison to the conductive heat flux driven by temperature gradients in non-local scenarios. The RKM code can be used for investigating other kinetic and non-local effects in a broader plasma physics context. Due to its relatively low computational cost it can also serve as a practical non-local ion heat flux closure in hydrodynamic simulations or as a training tool for machine learning surrogates.