Radiative cooling time derivation (for ambient temperature = 0 °K)

[javascripter]

Homework Statement

As part of a space experiment, a small instrumentation package is released from a space vehicle. It can be approximated as a solid aluminum sphere, 4 cm in diameter. The sphere is initially at 30°C and it contains a pressurized hydrogen component that will condense and malfunction at 30 K. If we take the surrounding space to e 0 K, how long may we expect the implementation package to function properly? Is it legitimate to use the lumped-capacity method in solving the problem? (Hint: see the directions for problem 1.17.) [Time = 5.8 weeks.]

I am specifically only interested in the time part of the question. For completeness, here is the relevant parts of 1.17:

A blackened copper sphere 2 cm in diameter and uniformy at 200°C is introduced into an evacuated black chamber that is maintained at 20°C.
- Write a differential equation that expresses T(t) for the sphere, assuming lumped thermal capacity.
- Integrate your differential equation and plot the temperature response for the sphere

Relevant Equations

$$Q = FεAσ(T_∞^4-T^4)$$
Q is heat (watts)
F is view factor (for this problem it is 1 but I leave it in because it doesn't make the math any more difficult)
ε is emmisivity (I'm not sure what the problem uses, but I think it is close to 1)
A is surface area
σ is a boltzman constant being 5.69*10^-8
T(∞) is ambient temperature (°K)
T is the temperature (°K)

$$Q = \frac{dU}{dt} = mc \frac{dT}{dt}$$
m is mass (kg)
c is specific heat capacity (J/kg°K)

as $T_∞ = 0$, I use $Q=−F\epsilon A\sigma T^4$ for this problem as the integration is much easier, so we start with
$mcdTdt=−F\epsilon A\sigma T^4$
rearrange so T is on one side and t is on the other:
$−mcF\epsilon A\sigma T^4dT=dt$

on the left side I integrate from the initial temperature (Ti) to the final temperature (Tf), on the write I integrate from 0 to t with t being the time for the sphere to cool from initial temperature to final temperature

$$\frac{mc}{3F\epsilon A\sigma} (\frac{1}{T_f^3}−\frac{1}{T_i^3}) = t$$

now I calculate all the constants and plug them in:

density of aluminum = 2707 kg/m^3
volume of sphere = 2.68*10^-4 m^3 (radius is .04 m)
m = .725476 kg

F = 1
A = .02 m^2
c = 905 J/kg°K

ε = 1 (is not correct, but should be close to the correct answer and t should be slightly lower than 5.8 weeks with this)

making the substitutions and calculating, I get:
t = 7.116*10^6 seconds = 11.77 weeks

I've been trying to find what the equation should look like and have found this: http://hyperphysics.phy-astr.gsu.edu/hbase/thermo/cootime.html
it looks pretty close to what I got, but they use 1.5Nk instead of mc; plugging my numbers into that calculator gives about 6 weeks, while I would expect it to be less than 5.8 weeks, but it is still much closer than my result

another thing I'm not sure about is that the specific heat is given at a constant pressure instead of constant volume, and is given for 20°C, so I'm not sure if I should or how to factor that in somehow. I know specific heat can be specified at constant volume, but I can't find any value of it for aluminum.

thanks for any help