It is hypothesized that spallation could affect the aerodynamic heating rates of reentry vehicles using ablative heat shields. To investigate spallation effects, a code is developed to compute the dynamics of spalled particles. The code uses a finite-rate chemistry model to study the chemical interactions of the particles with the flow field. The spallation code is one-way coupled to a computational fluid dynamics solver that models the hypersonic flow field around an ablative sample. Spalled particle behavior is numerically studied in both argon and air flow fields. Nomenclature A p = area of cross section of the particle, m 2 A s = surface area of the particle, m 2 C D = drag coefficient c i = local concentration of species i at particle surface, mol∕m 3 c v = specific heat at constant volume, J∕mol · K d p = diameter of the particle, m E p = internal energy of the particle, J∕kg E R = energy barrier for reaction, J∕mol F D = drag force, N J C i = vapor flux of C i species, kg∕m 2 · s k f = forward reaction rate, m∕s M = Mach number M w = molar weight, kg∕mol m p = mass of the particle, kg _ m C∕i = mass rate of the particle due to reactions, kg∕s Nu = Nusselt number P = vapor pressure, N∕m 2 p drag = drag power, J∕s _ q = heat rate, J∕s R = universal gas constant Re = Reynolds number based on relative velocity S = molecular speed ratio T = temperature, K U = state vector U ∞ = freestream velocity, m∕s u; v; w = velocity components in axial, radial, and z directions, m∕s V r = relative velocity between particle and fluid, m∕s W = source vector x; y; z = position components in axial, radial, and z directions, m Y i = mass fraction, kg∕kg α v = vaporization coefficient γ 0 = reaction efficiency ΔG f = Gibbs free energy of formation, J∕mol Δh = enthalpy of formation, J∕mol ϵ = emissivity of the particle κ = thermal conductivity, W∕m · K μ = dynamic viscosity, kg∕m · s ν = mean thermal speed, m∕s ρ = density, kg∕m 3 σ = Stefan-Boltzmann constantSubscripts conv = convection f = flow field g = gas phase p = particle rad = radiation rxn = chemical reactions s = solid phase (graphite) sub = sublimation tr = translational-rotational energy mode ve = vibrational-electronic energy mode