4
$\begingroup$

I'm trying to find the temperature of a CubeSat in orbit at 140km. I'm trying to use the Stefan-Boltzmann law but don't know how to incorporate the energy from sunlight. I'm using this as a check for a code that models temperature of the spacecraft through time and am considering the case where the plate is directly between the Earth and the Sun.

If the plate is at some angle, would I have to consider the dot product of the surface vector and the sunlight to consider the "strength" of the heat transferred? (i.e. 90 degrees = max heat transfer, 0 degrees = no heat transfer)

Assume orbit at 140km, the sheet is thin and has dimensions 10x10cm, emissivity of 0.09 (aluminum), reflectivity of 0.9, and the power from the sunlight is 1360 W/m^2.

enter image description here

$\endgroup$
8
  • $\begingroup$ Reflectivity? Emissivity? $\endgroup$ Commented Apr 17, 2019 at 21:12
  • $\begingroup$ You may find some helpful bits here and here but this is a different problem. By the way heating from the Earth is a consideration as well, as it takes up a considerable amount of solid angle. $\endgroup$
    – uhoh
    Commented Apr 18, 2019 at 0:58
  • 4
    $\begingroup$ At 140km, re-entry is nigh! This means temperature is pretty hot due to convection heat exchange, I wouldn't even bother with radiation/Stefan-Boltzmann. $\endgroup$
    – Mefitico
    Commented Apr 18, 2019 at 12:52
  • $\begingroup$ I adjusted the solar constant from 1000 to 1360 W/m^2 but I didn't notice the orbit altitude. That's a strangely low altitude to choose for a calculation about sunlight. $\endgroup$
    – uhoh
    Commented Apr 19, 2019 at 0:30
  • 1
    $\begingroup$ @uhoh is that adjusted "solar" constant meant to include the radiation from Earth? The increment is basically the same as I roughly calculate for Earth emission, including an emissivity fairly close to 1.0 form Earth's surfaces. $\endgroup$ Commented Apr 20, 2019 at 17:39

1 Answer 1

3
$\begingroup$

As you want to apply Stefan-Boltzmann, I will assume that you are interested in vacuum heat and radiation dynamics. As some comments have pointed out, at 140 km, interactions with the atmosphere would be very significant or dominant. That doesn't quite look like what you are after though, so I will completely ignore the atmosphere.

For a starting point, let's first also ignore the Earth. In a vacuum, exposed to the Sun, any object will reach an equilibrium where $\Delta W = 0$, that is

$$W_{in} = -W_{out}$$

For $W_{in}$, that's the solar constant of 1360 W/m², times the cross section area and what portion of the light gets absorbed (for real satellites, you may also want to add a term for the heat produced onboard).

$$W_{in} = G_{SC}A_1(1 - R)$$

$A_1$, the cross section area, is for the case of planar object, $\cos(\alpha)A$, where $\alpha$ is complementary to the angle in your diagram.

For the radiation, we can apply Stefan-Boltzmann directly

$$W_{out} = -A_2\epsilon \sigma T^4$$

$A_2$, the area that radiates heat, does not have to be equal to $A_1$. In this case it's both the front and the back of the aluminium sheet, so $2A$.

This gives us

$$ G_{SC}A_1(1 - R) = A_2\epsilon \sigma T^4$$

Solving for $T$ yields the equilibrium temperature

$$ T = \sqrt[\huge{4}]{\frac{G_{SC}A_1(1 - R)}{A_2\epsilon \sigma}}$$

But this is still just an object floating in space. The next thing to factor in is the night side of the Earth.

For an object with large thermal mass, the temperature will stay approximately the same over the time it takes to complete an orbit. For the simplest case of an object orbiting in the ecliptic plane, the shadow angle of the earth can be expressed as:

$$\beta = 2\sin^{-1}\left(\frac{r_e}{r}\right)$$

Where $r$ is the orbital radius and $r_e$ the radius of the Earth. For large orbits, you may want to take into account some shadow cone effects due to the physical size of the Sun, but at that point the Earth's shadow isn't really relevant any more.

You could make use of this immediately, and calculate an "average" adjusted solar constant, but there's still one important factor we have not yet taken into account:

Earth shine

The Earth both reflects and emits radiation, in fact, exactly as much as it receives, because it too is in a thermal equilibrium. As a first order approximation, the Earth receives radiation in a circular cross section, and scatters it in a sphere, the area ratio being 1 : 4 (in practice, since the Earth does not have perfect heat transport, the number will be somewhere between $\pi$ and 4 for the equator) For a low Earth orbit, this means around 300 W/m². You can scale this by multiplying by how much the radiation spreads, a factor of $\frac{r^2}{r_e^2}$

When adjusting for this, you have a usable model, at least for the case of large thermal mass. A piece of aluminium foil doesn't, so we have to go one step further.

For a small thermal mass, it will quickly reach a local thermal equilibrium, so it will have a day-side temperature and a night side temperature. Again, there would be some twilight effects and delay in practice, but at that point we would require numerical simulations. The temperature variation would in any case be between the two extreme temperatures and only dampen them somewhat.

The Earth does not directly absorb all the radiation that hits it. We have to subtract what's reflected throught eh bond albedo, which leaves about 2/3 absorbed. It's also a little colder at night most places, so a surface flux of 200 W/m² on the night side should be okay for calculating that equilibrium.

On the day side, we will still have this have this flux on top of the Sun light, and we will also have reflected light from the Earth, by area ratio of diffuse reflection about half of the bond albedo portion, or another 200 W/m²

$\endgroup$
2
  • $\begingroup$ What is the value of (1-R) in the 2nd equation? If the satellite is directly facing the sun, does that mean the value of (1-R) be equal to 1? $\endgroup$
    – Auberron
    Commented Jan 17, 2022 at 4:57
  • 1
    $\begingroup$ @Auberron R is the reflectivity, so this could have any value from 0 to 1. $\endgroup$ Commented Jan 17, 2022 at 9:04

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.