A multiscale approach is used to model and analyze the ablation of porous materials. Models are developed for the oxidation of a carbon preform and of the char layer of two phenolic impregnated carbon ablators with the same chemical composition but with different structures. Oxygen diffusion through the pores of the materials and in depth oxidation and mass loss are first modeled at the microscopic scale. The microscopic model is then averaged to yield a set of partial differential equations describing the macroscopic behavior of the material. Microscopic and macroscopic approaches are applied with progressive degrees of complexity to gain a comprehensive understanding of the ablation process. Porous medium ablation is found to occur in a zone of the char layer that we call the ablation zone. The thickness of the ablation zone is a decreasing function of the Thiele number. The studied materials are shown to display different ablation behaviors, a fact not captured by current models that are based on chemical composition only. Applied to Stardust's phenolic impregnated carbon ablator, the models explain and reproduce the unexpected drop in density measured in the char layer during Stardust postflight analyses [Stackpoole, M., Sepka, S., Cozmuta, I., and Kontinos, D., "Post-Flight Evaluation of Stardust Sample Return Capsule Forebody Heat-Shield Material," AIAA Paper 2008-1202, Jan. 2008]. Nomenclature C = oxygen concentration, mol=m 3 D = diffusion coefficient, m 2 =s D dis = dispersion coefficient, m 2 =s d p = mean pore diameter, m f = rarefaction function J = molar ablation rate, mol=m 2 =s Kn = Knudsen number k = reactivity, m=s L D = diffusion length, m L R = reaction length, m L 0 = sample initial length, or thickness, m M = molar mass, kg=mol n = vector normal to the surface P = pressure, Pa Pe = Péclet number P f = local fiber perimeter, m q = local molar flux density, mol=m 2 =s R = ideal gas constant, 8:314 J=mol=K Sx; y; z; t = surface function S p = local horizontal section of a pore, m 2 s = volume surface, m 2 =m 3 T = temperature, K V = averaging volume, m 3 v = velocity, m=s v = molecular agitation mean velocity, m=s x, y, z, t = space (m) and time (s) coordinates = exponential coefficient for tortuosity law of the pore-filling matrix ERSA = effective reactive surface-area coefficient " = porosity = tortuosity = mean free path, m h 2 j i = mean square displacement in direction j, m = density, kg=m 3 = Thiele number = solid molar volume, m 3 =mol Subscripts B = bulk e, eff = effective ERSA = effective reactive surface area f = fiber g = pyrolysis gases i = indicator or index-number K = Knudsen m = matrix ref = reference w = wall