[Physics] How does temperature in a solid sphere change with time when moving through a gas

classical-mechanicsthermal conductivitythermal-radiation

I'm interested in the following problem: There is a solid sphere with radius $r$ and mass $m$ at temperature $T_{s0}$. It is moving at velocity $v_s$ through a gas of temperature $T_g$. How does the temperature of the sphere, $T_s$, change with time?

A sketch:

Sketch of a flying sphere

I think I have to consider the following effects:

  • Radiation losses
  • thermal conduction at the solid/gas interface
  • Transport of hot gas away from the sphere due to $v_s$

My Approach is:

$$ T_s(t) = T_{s0} + m c \Delta E_s(t)$$

Where $c$ is the specific heat capacity and $\Delta E_s(t)$ the change of thermal energy in the sphere. For the latter I've set up an equation which incorporates the Stefan Boltzmann law.

$$ \Delta E_s(t) = – \int_0^t \sigma A (T_s(t)^4-T_g^4) dt + \Delta E_\text{cond}(t) + \Delta E_\text{conv}(t)$$

So, assuming this equation and though process is correct, my question boils down to how to take the thermal conduction (and convection?) represented by $\Delta E_\text{cond}(t)$ and \Delta E_\text{conv}(t) in the upper equation into account.

Edit:
To clarify the magnitude of interest for $v_s$, $T_s$ and $T_g$: $v_s$ is about 0,5 to 5 m/s, the $T_s – T_g$ is in the range of 100 to 300 °C.

Best Answer

So first off, i'm going to have to break to you that this $\Delta E_\text{cond}(t) + \Delta E_\text{conv}(t)$

is by far the most complicated bit. There is a simple equation $\frac{q}{A} = h\Delta T$ where $q$ is the heat flux (W/m^2), $A$ is the surface area of the object, $\Delta T$ is the temperature difference between the bulk of the fluid and the object and $h$ is the heat transfer coefficient. Unfortunately $h$ is quite difficult to determine as it depends on the properties of the fluid and its flow regime as well as the geometry of the problem.

In order to determine the heat transfer coefficient we define a number of dimensionless quantities such as $Pr$; the Prandtl number. This is the ratio of momentum diffusivity to thermal diffusivity for the fluid. Effectively this describes the relative size of the thermal and velocity boundary layers for the fluid. For low $Pr$ conduction through the fluid dominates (Mercury has Pr of around 0.015). In the opposite case of engine oil momentum diffusivity is much larger so convection dominates (Pr of 100 to 40000). The Prandtl number for air and most other gases is around 0.7-0.8.

Obviously the Reynolds number for the fluid, $\Re _\text{L}$ is going to be important and is defined as $\Re = \frac{v L}{\nu}$

This describes the flow regime around the object where $L$ is a characteristic length and can be taken to be the diameter of the sphere in this case. $\nu$ is the kinematic viscosity and $v$ is the velocity of the object relative to the bulk of the fluid. For the appropriate substance and temperature you can find the kinematic viscosity of your gas and then calculate the appropriate Reynolds number.

Lastly we define the Nusselt number. $Nu _\text{L} = \frac{h L}{k}$

where $L$ is again the characteristic length, $h$ is the heat transfer coefficient and $k$ is the thermal conductivity of the fluid which is evaluated at the film temperature. This is defined as the arithmetic mean of the temperature of the bulk of the fluid $T _\text{g}$ also known as $T _\infty$ and the temperature at the solid boundary.

As a small aside I will mention the Biot number. $Bi = \frac{h L}{k _\text{b}}$

where the symbols have their usual meanings except $k _\text{b}$ is the thermal conductivity of the body (your sphere in this case). If this is 0.1 or less then the conduction throughout the solid body is sufficient to assume that the body has no internal temperature gradient (which simplifies your problem).

Now on to solving the problem. We want to express the Nusselt number in terms of the Prandtl and Reynolds numbers in order to find an expression for $h$. Unfortunately this is very difficult but luckily large cohorts of experimentalists have slaved away at this problem for us.

For external flow over a sphere, T. Yuge (I don't have an explicit reference for this see textbook link at the end) found the following correlation $Nu = 2 + 0.43 Re ^\frac{1}{4} $.

This is for Pr approximately equal to 1 and Reynolds number between 1 and 10000.

(I could not find another correlation when Pr is closer to the range expected for air and most other gases. You will have to either accept this approximation or find your own information sorry.)

Decide on appropriate characteristics of your fluid and solid. Calculate the Reynolds and hence Nusselt numbers and then the heat transfer coefficient. Use this to solve for the heat flux for an appropriate temperature difference and then put this in your integral with your radiation term to solve for the total change in energy. I think that about does it. Here is a textbook that goes into a lot of detail:

http://www.unimasr.net/ums/upload/files/2012/Sep/UniMasr.com_919e27ecea47b46d74dd7e268097b653.pdf

I hope this is helpful.

Related Question