Aerodynamics of flying snakes using an immersed boundary method