Understanding the properties of grain boundaries in polycrystalline superconductors is essential for optimizing their critical current density. Here, we provide computational simulations of 2D Josephson junctions (JJs) in low magnetic fields using time--dependent Ginzburg--Landau theory, since they can be considered a proxy for a grain boundary between two grains.