Difference between revisions of "Newton’s law of Cooling models"

From Department of Mathematics at UTSA
Jump to navigation Jump to search
(fixing link to go straight to pdf)
Line 1: Line 1:
 +
'''Newton's law of cooling''' states that ''the rate of heat loss of a body is directly proportional to the difference in the temperatures between the body and its surroundings.'' The law is frequently qualified to include the condition that the temperature difference is small and the nature of heat transfer mechanism remains the same. As such, it is equivalent to a statement that the heat transfer coefficient, which mediates between heat losses and temperature differences, is a constant. This condition is generally met in heat conduction (where it is guaranteed by Fourier's law) as the thermal conductivity of most materials is only weakly dependent on temperature. In convective heat transfer, Newton's Law is followed for forced air or pumped fluid cooling, where the properties of the fluid do not vary strongly with temperature, but it is only approximately true for buoyancy-driven convection, where the velocity of the flow increases with temperature difference. Finally, in the case of heat transfer by thermal radiation, Newton's law of cooling holds only for very small temperature differences.
 +
 +
When stated in terms of temperature differences, Newton's law (with several further simplifying assumptions, such as a low Biot number and a temperature-independent heat capacity) results in a simple differential equation expressing temperature-difference as a function of time. The solution to that equation describes an exponential decrease of temperature-difference over time. This characteristic decay of the temperature-difference is also associated with Newton's law of cooling.
 +
 +
 +
== Relationship to mechanism of cooling==
 +
Convection cooling is sometimes said to be governed by "Newton's law of cooling."  When the heat transfer coefficient is independent, or relatively independent, of the temperature difference between object and environment, Newton's law is followed. The law holds well for forced air and pumped liquid cooling, where the fluid velocity does not rise with increasing temperature difference. Newton's law is most closely obeyed in purely conduction-type cooling. However, the heat transfer coefficient is a function of the temperature difference in natural convective (buoyancy driven) heat transfer. In that case, Newton's law only approximates the result when the temperature difference is relatively small. Newton himself realized this limitation.
 +
 +
A correction to Newton's law concerning convection for larger temperature differentials by including an exponent, was made in 1817 by Dulong and Petit. (These men are better-known for their formulation of the [[Dulong–Petit law]] concerning the molar specific heat capacity of a crystal.)
 +
 +
Another situation that does not obey Newton's law is radiative heat transfer.  Radiative cooling is better described by the Stefan–Boltzmann law in which the heat transfer rate varies as the difference in the 4th powers of the absolute temperatures of the object and of its environment.
 +
 +
== Mathematical formulation of Newton's law ==
 +
 +
The statement of Newton's law used in the heat transfer literature puts into mathematics the idea that ''the rate of heat loss of a body is proportional to the difference in temperatures between the body and its surroundings''. For a temperature-independent heat transfer coefficient, the statement is:
 +
 +
<math display="block">\dot{Q} = h A \left(T(t) - T_\text{env}\right) = h  A \, \Delta T(t),</math>
 +
where
 +
* <math>\dot{Q}</math> is the rate of heat transfer out of the body (SI unit: watt),
 +
* <math>h</math> is the heat transfer coefficient (assumed independent of ''T'' and averaged over the surface) (SI unit: W/m<sup>2</sup>⋅K),
 +
* <math>A</math> is the heat transfer surface area (SI unit: m<sup>2</sup>),
 +
* <math>T</math> is the temperature of the object's surface (SI unit: K),
 +
* <math>T_\text{env}</math> is the temperature of the environment; i.e., the temperature suitably far from the surface (SI unit: K),
 +
* <math>\Delta T(t) = T(t) - T_\text{env}</math> is the time-dependent temperature difference between environment and object (SI unit: K).
 +
 +
The heat transfer coefficient ''h'' depends upon physical properties of the fluid and the physical situation in which convection occurs. Therefore, a single usable heat transfer coefficient (one that does not vary significantly across the temperature-difference ranges covered during cooling and heating) must be derived or found experimentally for every system that is to be analyzed.
 +
 +
Formulas and correlations are available in many references to calculate heat transfer coefficients for typical configurations and fluids.  For laminar flows, the heat transfer coefficient is usually smaller than in turbulent flows because turbulent flows have strong mixing within the boundary layer on the heat transfer surface. Note the heat transfer coefficient changes in a system when a transition from laminar to turbulent flow occurs.
 +
 +
=== The Biot number ===
 +
 +
The Biot number, a dimensionless quantity, is defined for a body as
 +
 +
<math display="block">\text{Bi} = \frac{h L_{\rm C}}{k_{\rm b}},</math>
 +
where
 +
* ''h'' = film coefficient or heat transfer coefficient or convective heat transfer coefficient,
 +
* ''L''<sub>C</sub> = characteristic length, which is commonly defined as the volume of the body divided by the surface area of the body, such that <math>L_{\rm C} = V_\text{body} / A_\text{surface}</math>,
 +
* ''k''<sub>b</sub> = thermal conductivity of the body.
 +
 +
The physical significance of Biot number can be understood by imagining the heat flow from a hot metal sphere suddenly immersed in a pool to the surrounding fluid. The heat flow experiences two resistances: the first outside the surface of the sphere, and the second within the solid metal (which is influenced by both the size and composition of the sphere). The ratio of these resistances is the dimensionless Biot number.
 +
 +
If the thermal resistance at the fluid/sphere interface exceeds that thermal resistance offered by the interior of the metal sphere, the Biot number will be less than one. For systems where it is much less than one, the interior of the sphere may be presumed always to have the same temperature, although this temperature may be changing, as heat passes into the sphere from the surface. The equation to describe this change in (relatively uniform) temperature inside the object, is the simple exponential one described in Newton's law of cooling expressed in terms of temperature difference (see below).
 +
 +
In contrast, the metal sphere may be large, causing the characteristic length to increase to the point that the Biot number is larger than one. In this case, temperature gradients within the sphere become important, even though the sphere material is a good conductor. Equivalently, if the sphere is made of a thermally insulating (poorly conductive) material, such as wood or styrofoam, the interior resistance to heat flow will exceed that at the fluid/sphere boundary, even with a much smaller sphere. In this case, again, the Biot number will be greater than one.
 +
 +
Values of the Biot number smaller than 0.1 imply that the heat conduction inside the body is much faster than the heat convection away from its surface, and temperature gradients are negligible inside of it.  This can indicate the applicability (or inapplicability) of certain methods of solving transient heat transfer problems.  For example, a Biot number less than 0.1 typically indicates less than 5% error will be present when assuming a lumped-capacitance model of transient heat transfer (also called lumped system analysis). Typically, this type of analysis leads to simple exponential heating or cooling behavior ("Newtonian" cooling or heating) since the internal energy of the body is directly proportional to its temperature, which in turn determines the rate of heat transfer into or out of it. This leads to a simple first-order differential equation which describes heat transfer in these systems.
 +
 +
Having a Biot number smaller than 0.1 labels a substance as "thermally thin," and temperature can be assumed to be constant throughout the material's volume. The opposite is also true: A Biot number greater than 0.1 (a "thermally thick" substance) indicates that one cannot make this assumption, and more complicated heat transfer equations for "transient heat conduction" will be required to describe the time-varying and non-spatially-uniform temperature field within the material body. Analytic methods for handling these problems, which may exist for simple geometric shapes and uniform material thermal conductivity, are described in the article on the heat equation.
 +
 +
== Application of Newton's law transient cooling ==
 +
 +
Simple solutions for transient cooling of an object may be obtained when the internal thermal resistance within the object is small in comparison to the resistance to heat transfer away from the object's surface (by external conduction or convection), which is the condition for which the Biot number is less than about 0.1. This condition allows the presumption of a single, approximately uniform temperature inside the body, which varies in time but not with position. (Otherwise the body would have many different temperatures inside it at any one time.) This single temperature will generally change exponentially as time progresses (see below).
 +
 +
The condition of low Biot number leads to the so-called lumped capacitance model. In this model, the internal energy (the amount of thermal energy in the body) is calculated by assuming a constant heat capacity. In that case, the internal energy of the body is a linear function of the body's single internal temperature.
 +
 +
The lumped capacitance solution that follows assumes a constant heat transfer coefficient, as would be the case in forced convection. For free convection, the lumped capacitance model can be solved with a heat transfer coefficient that varies with temperature difference.
 +
 +
=== First-order transient response of lumped-capacitance objects ===
 +
 +
A body treated as a lumped capacitance object, with a total internal energy of <math>U</math> (in joules), is characterized by a single uniform internal temperature, <math>T(t)</math>. The heat capacitance, <math>C</math>, of the body is <math>C = dU/dT</math> (in J/K), for the case of an incompressible material. The internal energy may be written in terms of the temperature of the body, the heat capacitance (taken to be independent of temperature), and a reference temperature at which the internal energy is zero:  <math>U = C (T - T_\text{ref})</math>.
 +
 +
Differentiating <math>U</math> with respect to time gives:
 +
<math display="block">\frac{dU}{dt}  = C \, \frac{dT}{dt}.</math>
 +
 +
Applying the first law of thermodynamics to the lumped object gives <math display="inline">\frac{dU}{dt} = -\dot{Q}</math>, where the rate of heat transfer out of the body, <math>\dot{Q}</math>, may be expressed by Newton's law of cooling, and where no work transfer occurs for an incompressible material. Thus,
 +
<math display="block">\frac{dT(t)}{dt} = -\frac{hA}{C} (T(t) - T_\text{env}) = -\frac{1}{\tau}\ \Delta T(t),</math>
 +
where the time constant of the system is <math>\tau = C / (hA)</math>. The heat capacitance <math>C</math> may be written in terms of the object's specific heat capacity, <math>c</math> (J/kg-K), and mass, <math>m</math> (kg). The time constant is then <math>\tau = mc / (hA)</math>.
 +
 +
When the environmental temperature is constant in time, we may define <math>\Delta T(t) = T(t) - T_\text{env}</math>. The equation becomes
 +
<math display="block">\frac{dT(t)}{dt} = \frac{d\Delta T(t)}{dt} = -\frac{1}{\tau} \Delta T(t).</math>
 +
 +
The solution of this differential equation, by integration from the initial condition, is
 +
<math display="block">\Delta T(t) = \Delta T(0) \, e^{-t / \tau}.</math>
 +
where <math>\Delta T(0)</math> is the temperature difference at time 0. Reverting to temperature, the solution is
 +
<math display="block">T(t) = T_\text{env} + (T(0) - T_\text{env}) \, e^{-t/\tau}.</math>
 +
 +
The temperature difference between the body and the environment decays exponentially as a function of time.
 +
 +
==References==
 +
# Whewell, William (1866). History of the Inductive Sciences from the Earliest to the Present Times.
 +
# Lienhard, John H., IV; Lienhard, John H., V (2019). "Laminar and turbulent boundary layers". A Heat Transfer Textbook (5th ed.). Mineola, NY: Dover Publications. pp. 271–347. ISBN 9780486837352.
 +
# Frank Incropera; Theodore L. Bergman; David DeWitt; Adrienne S. Lavine (2007). Fundamentals of Heat and Mass Transfer (6th ed.). John Wiley & Sons. pp. 260–261. ISBN 978-0-471-45728-2.
 +
# Lienhard, John H., IV; Lienhard, John H., V (2019). A Heat Transfer Textbook (5th ed.). Mineola, NY: Dover Publications. pp. 419–420. ISBN 9780486837352.
 +
 +
==Resources==
 
* [https://mathresearch.utsa.edu/wikiFiles/MAT1093/Exponential%20growth%20and%20decay%20models/Esparza%201093%20Notes%207.6.pdf Newton’s law of Cooling models]. Written notes created by Professor Esparza, UTSA.
 
* [https://mathresearch.utsa.edu/wikiFiles/MAT1093/Exponential%20growth%20and%20decay%20models/Esparza%201093%20Notes%207.6.pdf Newton’s law of Cooling models]. Written notes created by Professor Esparza, UTSA.
 
* [https://mathresearch.utsa.edu/wikiFiles/MAT1093/Exponential%20growth%20and%20decay%20models/Esparza%201093%20Notes%207.6B.pdf Newton’s law of Cooling models Continued]. Written notes created by Professor Esparza, UTSA.
 
* [https://mathresearch.utsa.edu/wikiFiles/MAT1093/Exponential%20growth%20and%20decay%20models/Esparza%201093%20Notes%207.6B.pdf Newton’s law of Cooling models Continued]. Written notes created by Professor Esparza, UTSA.

Revision as of 00:25, 9 October 2021

Newton's law of cooling states that the rate of heat loss of a body is directly proportional to the difference in the temperatures between the body and its surroundings. The law is frequently qualified to include the condition that the temperature difference is small and the nature of heat transfer mechanism remains the same. As such, it is equivalent to a statement that the heat transfer coefficient, which mediates between heat losses and temperature differences, is a constant. This condition is generally met in heat conduction (where it is guaranteed by Fourier's law) as the thermal conductivity of most materials is only weakly dependent on temperature. In convective heat transfer, Newton's Law is followed for forced air or pumped fluid cooling, where the properties of the fluid do not vary strongly with temperature, but it is only approximately true for buoyancy-driven convection, where the velocity of the flow increases with temperature difference. Finally, in the case of heat transfer by thermal radiation, Newton's law of cooling holds only for very small temperature differences.

When stated in terms of temperature differences, Newton's law (with several further simplifying assumptions, such as a low Biot number and a temperature-independent heat capacity) results in a simple differential equation expressing temperature-difference as a function of time. The solution to that equation describes an exponential decrease of temperature-difference over time. This characteristic decay of the temperature-difference is also associated with Newton's law of cooling.


Relationship to mechanism of cooling

Convection cooling is sometimes said to be governed by "Newton's law of cooling." When the heat transfer coefficient is independent, or relatively independent, of the temperature difference between object and environment, Newton's law is followed. The law holds well for forced air and pumped liquid cooling, where the fluid velocity does not rise with increasing temperature difference. Newton's law is most closely obeyed in purely conduction-type cooling. However, the heat transfer coefficient is a function of the temperature difference in natural convective (buoyancy driven) heat transfer. In that case, Newton's law only approximates the result when the temperature difference is relatively small. Newton himself realized this limitation.

A correction to Newton's law concerning convection for larger temperature differentials by including an exponent, was made in 1817 by Dulong and Petit. (These men are better-known for their formulation of the Dulong–Petit law concerning the molar specific heat capacity of a crystal.)

Another situation that does not obey Newton's law is radiative heat transfer. Radiative cooling is better described by the Stefan–Boltzmann law in which the heat transfer rate varies as the difference in the 4th powers of the absolute temperatures of the object and of its environment.

Mathematical formulation of Newton's law

The statement of Newton's law used in the heat transfer literature puts into mathematics the idea that the rate of heat loss of a body is proportional to the difference in temperatures between the body and its surroundings. For a temperature-independent heat transfer coefficient, the statement is:

where

  • is the rate of heat transfer out of the body (SI unit: watt),
  • is the heat transfer coefficient (assumed independent of T and averaged over the surface) (SI unit: W/m2⋅K),
  • is the heat transfer surface area (SI unit: m2),
  • is the temperature of the object's surface (SI unit: K),
  • is the temperature of the environment; i.e., the temperature suitably far from the surface (SI unit: K),
  • is the time-dependent temperature difference between environment and object (SI unit: K).

The heat transfer coefficient h depends upon physical properties of the fluid and the physical situation in which convection occurs. Therefore, a single usable heat transfer coefficient (one that does not vary significantly across the temperature-difference ranges covered during cooling and heating) must be derived or found experimentally for every system that is to be analyzed.

Formulas and correlations are available in many references to calculate heat transfer coefficients for typical configurations and fluids. For laminar flows, the heat transfer coefficient is usually smaller than in turbulent flows because turbulent flows have strong mixing within the boundary layer on the heat transfer surface. Note the heat transfer coefficient changes in a system when a transition from laminar to turbulent flow occurs.

The Biot number

The Biot number, a dimensionless quantity, is defined for a body as

where

  • h = film coefficient or heat transfer coefficient or convective heat transfer coefficient,
  • LC = characteristic length, which is commonly defined as the volume of the body divided by the surface area of the body, such that ,
  • kb = thermal conductivity of the body.

The physical significance of Biot number can be understood by imagining the heat flow from a hot metal sphere suddenly immersed in a pool to the surrounding fluid. The heat flow experiences two resistances: the first outside the surface of the sphere, and the second within the solid metal (which is influenced by both the size and composition of the sphere). The ratio of these resistances is the dimensionless Biot number.

If the thermal resistance at the fluid/sphere interface exceeds that thermal resistance offered by the interior of the metal sphere, the Biot number will be less than one. For systems where it is much less than one, the interior of the sphere may be presumed always to have the same temperature, although this temperature may be changing, as heat passes into the sphere from the surface. The equation to describe this change in (relatively uniform) temperature inside the object, is the simple exponential one described in Newton's law of cooling expressed in terms of temperature difference (see below).

In contrast, the metal sphere may be large, causing the characteristic length to increase to the point that the Biot number is larger than one. In this case, temperature gradients within the sphere become important, even though the sphere material is a good conductor. Equivalently, if the sphere is made of a thermally insulating (poorly conductive) material, such as wood or styrofoam, the interior resistance to heat flow will exceed that at the fluid/sphere boundary, even with a much smaller sphere. In this case, again, the Biot number will be greater than one.

Values of the Biot number smaller than 0.1 imply that the heat conduction inside the body is much faster than the heat convection away from its surface, and temperature gradients are negligible inside of it. This can indicate the applicability (or inapplicability) of certain methods of solving transient heat transfer problems. For example, a Biot number less than 0.1 typically indicates less than 5% error will be present when assuming a lumped-capacitance model of transient heat transfer (also called lumped system analysis). Typically, this type of analysis leads to simple exponential heating or cooling behavior ("Newtonian" cooling or heating) since the internal energy of the body is directly proportional to its temperature, which in turn determines the rate of heat transfer into or out of it. This leads to a simple first-order differential equation which describes heat transfer in these systems.

Having a Biot number smaller than 0.1 labels a substance as "thermally thin," and temperature can be assumed to be constant throughout the material's volume. The opposite is also true: A Biot number greater than 0.1 (a "thermally thick" substance) indicates that one cannot make this assumption, and more complicated heat transfer equations for "transient heat conduction" will be required to describe the time-varying and non-spatially-uniform temperature field within the material body. Analytic methods for handling these problems, which may exist for simple geometric shapes and uniform material thermal conductivity, are described in the article on the heat equation.

Application of Newton's law transient cooling

Simple solutions for transient cooling of an object may be obtained when the internal thermal resistance within the object is small in comparison to the resistance to heat transfer away from the object's surface (by external conduction or convection), which is the condition for which the Biot number is less than about 0.1. This condition allows the presumption of a single, approximately uniform temperature inside the body, which varies in time but not with position. (Otherwise the body would have many different temperatures inside it at any one time.) This single temperature will generally change exponentially as time progresses (see below).

The condition of low Biot number leads to the so-called lumped capacitance model. In this model, the internal energy (the amount of thermal energy in the body) is calculated by assuming a constant heat capacity. In that case, the internal energy of the body is a linear function of the body's single internal temperature.

The lumped capacitance solution that follows assumes a constant heat transfer coefficient, as would be the case in forced convection. For free convection, the lumped capacitance model can be solved with a heat transfer coefficient that varies with temperature difference.

First-order transient response of lumped-capacitance objects

A body treated as a lumped capacitance object, with a total internal energy of (in joules), is characterized by a single uniform internal temperature, . The heat capacitance, , of the body is (in J/K), for the case of an incompressible material. The internal energy may be written in terms of the temperature of the body, the heat capacitance (taken to be independent of temperature), and a reference temperature at which the internal energy is zero: .

Differentiating with respect to time gives:

Applying the first law of thermodynamics to the lumped object gives , where the rate of heat transfer out of the body, , may be expressed by Newton's law of cooling, and where no work transfer occurs for an incompressible material. Thus,

where the time constant of the system is . The heat capacitance may be written in terms of the object's specific heat capacity, (J/kg-K), and mass, (kg). The time constant is then .

When the environmental temperature is constant in time, we may define . The equation becomes

The solution of this differential equation, by integration from the initial condition, is

where is the temperature difference at time 0. Reverting to temperature, the solution is

The temperature difference between the body and the environment decays exponentially as a function of time.

References

  1. Whewell, William (1866). History of the Inductive Sciences from the Earliest to the Present Times.
  2. Lienhard, John H., IV; Lienhard, John H., V (2019). "Laminar and turbulent boundary layers". A Heat Transfer Textbook (5th ed.). Mineola, NY: Dover Publications. pp. 271–347. ISBN 9780486837352.
  3. Frank Incropera; Theodore L. Bergman; David DeWitt; Adrienne S. Lavine (2007). Fundamentals of Heat and Mass Transfer (6th ed.). John Wiley & Sons. pp. 260–261. ISBN 978-0-471-45728-2.
  4. Lienhard, John H., IV; Lienhard, John H., V (2019). A Heat Transfer Textbook (5th ed.). Mineola, NY: Dover Publications. pp. 419–420. ISBN 9780486837352.

Resources