This study presents the results obtained by performing a set of direct numerical simulations (DNS) of periodic channel flows over an isothermal surface subjected to thermo-chemical ablation. A specific boundary condition to handle the surface ablation is developed and allows to represent any kind of ablation scheme. By nature, the periodic channel flow configuration is statically unsteady when ablation is present, which requires a particular attention when investigating flow statistics. An ensemble average over several flow realizations is then employed to generate converged statistics. Hence, this procedure allows to investigate different features of the flow when focussing on species/momentum/energy/atomic conservation. The analysis reveals that the flow can be considered at chemical equilibrium under the conditions investigated. Moreover, the convective effect introduced by the Stephan velocity at ablative surface appears to have no influence on momentum conservation whereas it strongly participates to the surface cooling effect. Finally, an atomic mass conservation equation of Shvab-Zel'dovich type is formulated, and highlights passive scalar like conservation mechanism of atomic elements.