A 2-D Lattice Boltzmann Method, designed to ensure stability at high Reynolds numbers, is combined with an Actuator Line Model to compute the loads on a two-bladed vertical axis wind turbine. Tests on the kernel size at a high mesh resolution reveal that a size equal to half of the full chord length yields the most accurate results. The aerodynamic load solution is validated against a fully resolved Scale-Adaptive Simulation (SAS) output, demonstrating high correlation, and enabling an assessment of near wake and downstream effects. The model’s adaptability to various rotor operating conditions is confirmed through tests at high and low tip-speed ratios. Additionally, a Biot–Savart-based Vortex Model (VM) is employed for further comparison, showing good agreement with the Lattice Boltzmann output. The results indicate that the Highly Stable Lattice Boltzmann Method integrated with the Actuator Line Model enhances the accuracy of flow field resolution and effectively captures complex aerodynamic phenomena, making it a valuable tool for simulating vertical axis wind turbines.