Adding radiation to the surface of a reacting particle
 

I want to add a radiation heat as a source to the boundary of a reacting particle.
Q rad = ε*sigma*F (Tw ^4 − Ts^ 4 )

how can i write the programme to identify the surface area(varies with time) to do this?