A numerical investigation of polar lows