Immersed boundary method for high-order flux reconstruction based on volume penalization