The objective of this research is to estimate two dimensional ablating and charring of heat shield materials in severe aero-thermal heat transfer. This estimation requires an accurate and rapid technique for its serious heat transfer with a moving boundary. Aerodynamic heating is obtained by an explicit relation which is a function of Mach number and air condition, while a fully implicit method is used for heat transfer calculations. Moving boundary is captured by FLIAR method which is a subgroup of VOF. Thickness of ablating and charring of heat shield, temperature of the moving surface and rate of radiation heat are calculated and compared with references. The results are in good agreement with other calculations.