Newton's Law of Cooling/Heating - Differential Equations

SECTION 1

Historical Context & Motivation

Have you ever noticed how a cup of hot coffee cools quickly at first, then seems to take forever to reach room temperature? In the late 1600s, Sir Isaac Newton began investigating exactly this phenomenon. He realized that the rate at which an object changes temperature depends on the difference between the object's temperature and the temperature of its surroundings. This insight eventually became one of the most widely used applications of first-order ordinary differential equations.

Newton's observation was remarkably practical. Before thermometers were precise, scientists needed a mathematical model to describe heat transfer. Newton's work gave them a clean, elegant equation that applies to everything from forensic science to engineering. Let's trace the key milestones that led to this foundational law.

1701 • Newton Publishes His Cooling Law

Isaac Newton anonymously publishes his observations on cooling in Philosophical Transactions, noting that the rate of heat loss is proportional to the temperature difference between the object and its environment.

1822 • Fourier's Analytical Theory of Heat

Joseph Fourier publishes Théorie analytique de la chaleur, providing a rigorous mathematical framework for heat conduction that places Newton's law within a broader context of partial differential equations.

1900s • Forensic & Engineering Applications

Newton's Law of Cooling becomes a standard tool in forensic science for estimating time of death, and in engineering for modeling the cooling of electronic components, industrial machinery, and food safety protocols.

Today • Modern Classroom & Industry Standard

The law remains a go-to example in differential equations courses worldwide. It illustrates how a simple first-order ODE can model real-world phenomena with surprising accuracy.

The central question Newton addressed is deceptively simple: If you know the current temperature of an object and the temperature of its surroundings, can you predict its temperature at any future time? The answer, as we will see, is a resounding yes — and the mathematics behind it is both elegant and accessible.

SECTION 2

Core Principles & Definitions

Newton's Law of Cooling/Heating rests on a few straightforward ideas. Before we dive into the math, let's make sure you understand the key concepts that underpin the entire model. Each of these principles connects directly to the differential equation you will learn to solve.

Temperature Difference Drives Change

The rate at which an object heats up or cools down is proportional to the difference between its temperature T and the ambient (surrounding) temperature T_s. A bigger gap means faster change.

Ambient Temperature Is Constant

The model assumes that the surrounding environment stays at a constant temperature T_s. For example, the room temperature remains steady while your coffee cools.

Exponential Approach to Equilibrium

The temperature doesn't drop (or rise) linearly. Instead, it follows an exponential curve that approaches T_s asymptotically — meaning it gets closer and closer but never quite reaches it.

The Cooling Constant k

The constant k captures how quickly the object exchanges heat with its environment. Materials, surface area, and insulation all affect k. A larger k means faster cooling or heating.

✦

KEY TAKEAWAY

Think of Newton's Law like a car braking to a stop. When you're going fast (big temperature difference), the brakes are very effective and you slow down quickly. As you approach a standstill (ambient temperature), you're barely decelerating. The temperature exponentially approaches the surrounding temperature, just as a braking car gradually approaches zero speed — fast at first, then slower and slower.

SECTION 3

Visual Explanation

The graph below shows how the temperature of an object changes over time according to Newton's Law of Cooling. Notice how the curve drops steeply at first, then gradually flattens as it approaches the ambient temperature. This shape — a decaying exponential — is the signature of Newton's cooling law.

The pink-to-cyan curve shows how the object's temperature drops from T₀ = 90°C toward the ambient temperature T_s = 20°C (yellow dashed line). The shaded area emphasizes the exponential decay shape. Notice how the curve is steepest near t = 0, where the temperature difference is largest.

This graph captures the essence of Newton's cooling law. At time t = 0, the object starts at its initial temperature T₀. The slope of the curve (how steeply it falls) is determined by the temperature difference between the object and the environment. As that difference shrinks, the curve flattens. The object will never actually reach exactly T_s in finite time — it only approaches it as an asymptote.

SECTION 4

Mathematical Framework

Now let's build the equation from the core principle. We said the rate of temperature change is proportional to the difference between the object's temperature and the ambient temperature. In math, "rate of change" translates directly to a derivative.

NEWTON'S LAW OF COOLING (DIFFERENTIAL FORM)

dT/dt = −k(T − Tₛ)

Where T = temperature of the object at time t, Tₛ = constant ambient (surrounding) temperature, k = positive cooling constant (units: 1/time), and t = time. The negative sign ensures that if T > Tₛ (object is hotter), the temperature decreases.

This is a separable first-order ODE. To solve it, we separate variables and integrate both sides. Let u = T − Tₛ, so du = dT (since Tₛ is constant). The equation becomes du/dt = −ku. Separating variables: du/u = −k dt. Integrating both sides gives ln|u| = −kt + C, which leads to u = Ae^−kt, where A = e^C.

GENERAL SOLUTION

T(t) = Tₛ + (T₀ − Tₛ) × e^(−kt)

Here T₀ is the initial temperature of the object at t = 0. The term (T₀ − Tₛ) represents the initial temperature difference, and $e^{−kt}$ is the exponential decay factor that shrinks that difference over time.

💡

Cooling vs. Heating

The same equation works for both cooling and heating! If the object starts cooler than its surroundings (T₀ < Tₛ), then (T₀ − Tₛ) is negative, so the exponential term adds a negative value that decays toward zero, making T(t) rise toward Tₛ. The math handles both directions automatically.

FINDING THE COOLING CONSTANT k

k = −(1/t) × ln[(T(t) − Tₛ) / (T₀ − Tₛ)]

If you know the temperature at some specific time t, you can solve for k by rearranging the general solution. This formula is essential for real-world problems where k is not given directly.

SECTION 5

Detailed Behavior & Comparison

Let's explore how the cooling constant k affects the shape of the temperature curve. Different materials and conditions produce very different values of k. The diagram below compares three scenarios with different k values, all starting at the same initial temperature and approaching the same ambient temperature.

Three cooling curves with different values of k. The red curve (k = 0.20) cools fastest, reaching near-ambient temperature in about 15 minutes. The violet curve (k = 0.10) takes about 25 minutes. The cyan curve (k = 0.05) is still noticeably above ambient at 25 minutes.

The value of k depends on the physical properties of the system. A thin metal cup of coffee sitting on a windy porch would have a large k because heat escapes quickly through the thin walls and is carried away by the wind. A thick, insulated thermos would have a very small k because it resists heat flow. In every case, the general shape of the curve is the same — only the rate of approach to Tₛ differs.

Scenario	Approximate k (min⁻¹)	Time to Halve ΔT
Hot soup in a bowl (no lid)	0.06 – 0.10	≈ 7 – 12 min
Coffee in a ceramic mug	0.03 – 0.06	≈ 12 – 23 min
Water in a thermos	0.003 – 0.01	≈ 70 – 230 min
Human body after death (forensics)	0.01 – 0.03	≈ 23 – 70 min

Typical k values for common cooling scenarios

SECTION 6

Worked Example

Let's work through a complete problem step by step. Pay attention to how we identify the known values, find the cooling constant, and then use the general solution to predict a future temperature.

Cooling Coffee Problem

1
Step 1 — Read the Problem & Identify Given Values

A fresh cup of coffee is poured at 95°C in a room that is 22°C. After 5 minutes, the coffee has cooled to 70°C. What will its temperature be after 15 minutes total?

T₀ = 95°C, Tₛ = 22°C, T(5) = 70°C, find T(15)
2
Step 2 — Write the General Solution

We start with the solution to Newton's cooling equation: T(t) = Tₛ + (T₀ − Tₛ) × e^−kt. Substituting known values: T(t) = 22 + (95 − 22) × e^−kt = 22 + 73e^−kt.

T(t) = 22 + $73e^{−kt}$
3
Step 3 — Use the Data Point to Find k

We know T(5) = 70°C. Plug in t = 5: 70 = 22 + 73e^−5k. Subtract 22: 48 = 73e^−5k. Divide by 73: e^−5k = 48/73 ≈ 0.6575. Take the natural log: −5k = ln(0.6575) ≈ −0.4198. Divide by −5: k ≈ 0.0840 min⁻¹.

k ≈ 0.0840 min⁻¹
4
Step 4 — Calculate T(15)

Now substitute k and t = 15 into the general solution: T(15) = 22 + 73e^{−0.0840 × 15} = 22 + 73e^−1.260. We calculate e^−1.260 ≈ 0.2837. So T(15) = 22 + 73 × 0.2837 ≈ 22 + 20.71 ≈ 42.7°C.

T(15) ≈ 42.7°C
5
Step 5 — Interpret the Result

After 15 minutes, the coffee has cooled from 95°C to about 42.7°C. That's a comfortable drinking temperature! Notice that the coffee lost 25°C in the first 5 minutes (from 95 to 70), but only lost about 27.3°C more over the next 10 minutes (from 70 to 42.7). This confirms the exponential decay behavior: cooling slows as the temperature gap narrows.

SECTION 7

Strengths & Limitations

Newton's Law of Cooling is powerful and widely applicable, but like every mathematical model, it has boundaries. Understanding when the model works well — and when it breaks down — is an important part of applied mathematics.

Strengths	Limitations
Simple, closed-form solution — no numerical methods needed	Assumes ambient temperature Tₛ stays constant; breaks down if the environment heats up or cools down
Applies to both cooling and heating with the same formula	Ignores radiation; works best for convection-dominated scenarios with moderate temperature differences
Requires only two data points (initial temp + one measurement) to calibrate	Assumes the object's temperature is uniform throughout (no internal temperature gradients)
Widely applicable: forensics, food safety, HVAC, electronics cooling	Becomes inaccurate for very large temperature differences where radiative heat loss dominates

When does Newton's Law of Cooling work well, and when doesn't it?

✦

KEY TAKEAWAY

Newton's Law of Cooling is like a weather forecast for temperature: it's remarkably good for short-term predictions under normal conditions, but it loses accuracy in extreme situations (like predicting the cooling of molten lava). For everyday objects — coffee, food, electronics, body temperature — it's an excellent approximation that balances simplicity with real predictive power.

SECTION 8

Connection to Advanced Theory

Newton's Law of Cooling is actually a special case of more general heat transfer theory. As you move forward in math and physics, you'll encounter these broader frameworks. The table below shows how Newton's model relates to more advanced treatments of heat flow.

Feature	Newton's Law (This Lesson)	Heat Equation (Advanced)
Equation type	First-order ODE	Partial differential equation (PDE)
Temperature distribution	Uniform throughout object	Varies with position inside the object
Heat transfer mode	Convection (surface)	Conduction + convection + radiation
Math required	Algebra and basic calculus	Multivariable calculus and Fourier series
Best for	Small objects, moderate ΔT	Any object, any temperature range

Newton's cooling law vs. the full heat equation

The key assumption that makes Newton's law so simple is called the lumped capacitance model: we treat the object as if it has a single, uniform temperature at every point. This works well for small objects or objects made of highly conductive materials (like metals). For large or insulating objects — say, a thick steak being cooked — the interior heats at a different rate than the surface, and you would need the full heat equation (a PDE studied in more advanced courses).

ℹ

Looking Ahead

The exponential decay pattern you see in Newton's Law of Cooling appears throughout science: radioactive decay, RC circuits in electronics, drug elimination in pharmacology, and population decay. Mastering this pattern here gives you a head start on all of these topics.

SECTION 9

Practice Problems

Test your understanding with these five problems, arranged from conceptual to challenging. Try each one on your own before checking the answer.

Question 1: CONCEPTUAL

A hot pie is placed on a counter in a 25°C kitchen. Without doing any calculation, explain why the pie cools rapidly at first and then more slowly as time goes on. Reference the differential equation dT/dt = −k(T − Tₛ) in your explanation.

Show answer

When the pie is very hot, the term (T − Tₛ) is large, making dT/dt a large negative number — the temperature drops quickly. As the pie cools and T gets closer to Tₛ = 25°C, the difference (T − Tₛ) shrinks, so dT/dt becomes a smaller negative number. The rate of cooling is always proportional to the current temperature gap, which is why the process slows down. This is the mathematical reason behind the characteristic exponential decay curve.

Question 2: BASIC CALCULATION

A cup of tea is brewed at 90°C and placed in a room at 20°C. The cooling constant is k = 0.05 min⁻¹. What is the temperature of the tea after 10 minutes?

Show answer

Use the general solution T(t) = Tₛ + (T₀ − Tₛ) ×

e^{−kt}

. Substituting: T(10) = 20 + (90 − 20) ×

e^{−0.05 × 10}

= 20 + 70 ×

e^{−0.5}

. Since

e^{−0.5}

≈ 0.6065, we get T(10) = 20 + 70 × 0.6065 = 20 + 42.46 ≈ 62.5°C. After 10 minutes, the tea has cooled from 90°C to about 62.5°C — it lost roughly 27.5°C, which is a significant drop because the temperature gap was still fairly large during this interval.

Question 3: INTERMEDIATE

A metal rod heated to 200°C is placed in a room at 25°C. After 8 minutes, the rod's temperature is 120°C. Find the cooling constant k, then determine how long it will take for the rod to reach 50°C.

Show answer

First, find k. Using T(8) = 120: 120 = 25 + (200 − 25) ×

e^{−8k}

, so 95 = 175 ×

e^{−8k}

, giving

e^{−8k}

= 95/175 ≈ 0.5429. Taking ln: −8k = ln(0.5429) ≈ −0.6114, so k ≈ 0.0764 min⁻¹. Next, find t when T(t) = 50: 50 = 25 + 175 ×

e^{−0.0764t}

, so 25 = 175 ×

e^{−0.0764t}

, giving

e^{−0.0764t}

= 25/175 ≈ 0.1429. Taking ln: −0.0764t = ln(0.1429) ≈ −1.9459, so t ≈ 25.5 minutes. It takes about 25.5 minutes for the rod to cool to 50°C.

Question 4: APPLIED

A forensic investigator finds a body at 10:00 PM. The body temperature is 31°C. One hour later at 11:00 PM, the body temperature is 29°C. The room is a steady 21°C. Normal body temperature is 37°C. Estimate the time of death.

Show answer

Use the data at 10:00 PM and 11:00 PM to find k. Let t = 0 be the time of death, with T₀ = 37°C and Tₛ = 21°C. At 10:00 PM (call this time t₁): 31 = 21 + 16 ×

e^{−kt₁}

, so

e^{−kt₁}

= 10/16 = 0.625. At 11:00 PM (time t₁ + 1): 29 = 21 + 16 ×

e^{−k(t₁+1)}

, so

e^{−k(t₁+1)}

= 8/16 = 0.5. Dividing the second equation by the first:

e^{−k}

= 0.5/0.625 = 0.8, so k = −ln(0.8) ≈ 0.2231 hr⁻¹. From

e^{−kt₁}

= 0.625: t₁ = −ln(0.625)/0.2231 ≈ 0.4700/0.2231 ≈ 2.11 hours. So the time of death was about 2.1 hours before 10:00 PM, which is approximately 7:53 PM.

Question 5: CRITICAL THINKING

Newton's Law of Cooling assumes the ambient temperature Tₛ is constant. Suppose instead that the room temperature slowly rises at a constant rate: Tₛ(t) = T_room + rt, where r is a small positive constant. How would you set up the new differential equation? Without solving it, describe qualitatively how the temperature curve would differ from the standard case. Would the object still approach a fixed asymptote?

Show answer

The modified ODE would be dT/dt = −k[T − (T_room + rt)], or equivalently dT/dt = −kT + k(T_room + rt). This is still a first-order linear ODE, but it is no longer separable in the simple way we used — you'd need an integrating factor. Qualitatively, the object would initially cool (or heat) in a similar exponential fashion, but instead of approaching a fixed horizontal asymptote, the temperature curve would eventually begin to track the rising ambient temperature with a slight lag. There is no fixed asymptote; instead, the object's temperature would rise without bound (following the room temperature), but always staying slightly below it. This demonstrates that changing the assumption about Tₛ fundamentally changes the long-term behavior of the solution.

SUMMARY

Lesson Summary

Summary

Newton's Law of Cooling/Heating states that the rate of temperature change of an object is proportional to the difference between its temperature and the ambient temperature. This relationship is expressed by the first-order ODE dT/dt = −k(T − Tₛ), which has the general solution T(t) = Tₛ + (T₀ − Tₛ) × $e^{−kt}$ . The cooling constant k determines how quickly the object approaches the surrounding temperature, and it can be found experimentally using a single temperature measurement at a known time.

The temperature follows an exponential decay curve, dropping (or rising) rapidly at first and then gradually approaching Tₛ as an asymptote. The model works for both cooling and heating and finds real-world applications in forensic science, food safety, and engineering. While it assumes constant ambient temperature and uniform object temperature (the lumped capacitance model), it remains one of the most practical and elegant applications of first-order differential equations.