An implicit ablation and thermal response program is presented for simulation of one-dimensional transient thermal energy transport in a multilayer stack of isotropic materials and structure that can ablate from a front surface and decompose in depth. The governing equations and numerical procedures for solution are summarized. Solutions are compared with those of an existing code, CMA, and also with arcjet data. Numerical experiments show that the new code is numerically more stable and solves a much wider range of problems compared with the older code. To demonstrate its capability, applications for thermal analysis and sizing of aeroshell heatshields for planetary missions of Stardust, Mars Microprobe (Deep Space II), Saturn Entry Probe, and Mars 2001, using advanced lightweight ceramic ablators developed at NASA Ames Research Center, are presented and discussed. Nomenclature a = absorption coef cient, m ¡1 B 0 = dimensionless mass blowing rate, P m=½ e u e C M B a = pre-exponential constant in Eq. (8), s ¡1 C H = Stanton number for heat transfer C M = Stanton number for mass transfer c p = speci c heat, J/kg-K E a = activation temperature in Eq. (8), K F = view factor g = outward pyrolysis mass ux, kg/m 2 -s h = enthalpy, J/kg N h = partial heat of charring, de ned in Eq. (6), J/kg I 0 = radiation source function in Eq. (2), W/m 2 -sr i C = radiant intensity in Cx direction, W/m 2 -sr i ¡ = radiant intensity in ¡x direction, W/m 2 -sr K = extinction coef cient, a C ¾ s , m ¡1 k = thermal conductivity, W/m-K P m = mass ux, kg/m 2 -s P = pressure, N/m 2 q C = conductive heat ux, W/m 2 q R = radiative heat ux, W/m 2 R = universal gas constant, J/kmol-K s = surface recession, m P s = surface recession rate, m/s T = temperature, K u = velocity, m/s x = moving coordinate, y ¡ s, m y = stationary coordinate, m Z ¤ = coef cient in Eq. (9), de ned in Ref. 10 ® = surface absorptance 0 = volume fraction of resin " = surface emissivity µ = time, s · = optical thickness · D = optical thickness for path of length Ḑ = blowing reduction parameter ½ = density, kg/m 3 ¾ = Stefan-Boltzmann constant, W/m 2 -K 4 ¾ s = scattering coef cient, m ¡1 ¿ = mass fraction of virgin material, de ned in Eq. (5) 9 = decomposition reaction order in Eq. (8) Subscripts c = char e = boundary-layeredge g = pyrolysis gas i = density component (A, B, and C ) j = surface species v = virgin w = wall