In this paper, numerical results of thermal blooming with 4-D code has been presented. Numerical simulations of phase compensation and field compensation have been carried out by means of using the 4-D code, which includes the beacon. The calculation results show that simple phase compensation is unstable, but field compensation is completely effective.