Leaky Boiler and Carnot Engine

Heat escaping through a boiler reduces the power supplied to a Carnot engine. Find two temperatures that give the same output and the temperature that gives the greatest useful power.

Before you begin

You should know conduction through a plane slab, steady energy balance, Carnot efficiency and quadratic equations. The maximum is found by differentiation and checked by algebra.

Conduction through the boiler walls reduces the heat reaching a Carnot engine. Find how this heat loss changes the useful power and the best operating temperature.

A hotter engine is more efficient, but a hotter boiler also loses more heat. We will find two temperatures that give the same output, then locate the temperature that gives the greatest useful power.

The problem

A cubic boiler supplies heat to an ideal Carnot engine, which rejects heat to surroundings at a constant temperature \(T_{0}\). The heater supplies a fixed net power \(P_{0}\) through the base. Heat also escapes directly through the four side walls and the roof. The engine drives a mechanical load, and its operation can be adjusted to obtain different steady boiler temperatures.

A schematic cubic boiler with a heated base, heat-conducting side walls and roof, and cooler surroundings.
Figure 1. The base supplies the stated net heater power. Heat escapes through the four side walls and the roof. The cube is schematic; each of the five loss faces has area \(a^{2}\).
QuantitySymbolValue
Edge length of the cube\(a\)\(20\,\mathrm{cm}\)
Thickness of each loss wall\(d\)\(1.0\,\mathrm{cm}\)
Wall thermal conductivity\(k\)\(1.0\,\mathrm{W\,m^{-1}\,K^{-1}}\)
Net heater power through the base\(P_0\)\(5.0\,\mathrm{kW}\)
Surroundings and cold reservoir\(T_0\)\(27^{\circ}\mathrm{C}\)

The boiler interior has a uniform temperature \(T\). Treat each loss wall as a plane slab of area \(a^{2}\), with its inner face at \(T\) and its outer face at \(T_{0}\). Use steady, one-dimensional conduction; neglect edge effects and additional heat losses. For numerical work, use the school approximation \(0^{\circ }\mathrm{C}=273\ \mathrm{K}\).

Find (a) the engine efficiency and mechanical power when the boiler is at \(100^{\circ }\mathrm{C}\).

(b) All boiler temperatures at which the mechanical power is \(500\ \mathrm{W}\).

(c) The greatest mechanical power obtainable, and the corresponding boiler temperature.

Heat balance and engine efficiency

Count the heat loss faces

Heat escapes through five faces, not six. The base supplies the specified net input. Each of the other faces provides a parallel route from \(T\) to \(T_{0}\), so their areas add:

\[ A=5a^{2}=5(0.20)^{2}=0.20\ \mathrm{m}^{2}. \tag{1}\]

Fourier’s law gives the total wall loss. Define \(G\) as the combined thermal conductance of the five walls:

\[ G=\frac{kA}{d}=\frac{(1.0)(0.20)}{0.010}=20\,\mathrm{W\,K^{-1}}. \tag{2}\]
\[ P_{\ell}=G(T-T_0). \tag{3}\]

The corresponding wall resistance is \(1/G=0.050\ \mathrm{K}\,\mathrm{W}^{-1}\). This is a parallel combination of wall resistances; heat does not pass through all five walls in succession.

Separate the two routes to the surroundings

Heat from the heater divides into direct wall leakage and heat supplied to a Carnot engine; the engine produces work and rejects heat to the same surroundings.
Figure 2. The two routes to the surroundings are distinct. Wall leakage \(P_{\ell}\) bypasses the engine; only \(P_h\) enters it. The engine delivers useful power \(P\) and rejects heat at rate \(P_c\).

Let \(P_{h}\) be the heat input per second to the engine, \(P_{c}\) its heat rejection per second, and \(P\) its useful mechanical power. Steady energy balance for the boiler gives

\[ P_{h}=P_{0}-P_{\ell }=P_{0}-G\left( T-T_{0} \right). \tag{4}\]

For a Carnot engine, use absolute temperatures:

\[ \eta=1-\frac{T_0}{T}. \tag{5}\]
\[ P=\eta P_h=\left(1-\frac{T_0}{T}\right)\left(P_0-G(T-T_0)\right). \tag{6}\]

The Carnot factor multiplies the engine’s heat input \(P_{h}\), not the heater power \(P_{0}\).

Part a at 100 degrees Celsius

Here \(T=373\ \mathrm{K}\) and \(T_{0}=300\ \mathrm{K}\). The wall temperature difference is \(73\ \mathrm{K}\):

\[ P_{\ell }=20(73)=1460\ \mathrm{W},\quad\quad P_{h}=5000-1460=3540\ \mathrm{W}. \tag{7}\]
\[ \eta=1-\frac{300}{373}=\frac{73}{373}=0.19571,\quad\quad P=3540\left( \frac{73}{373} \right)=692.82\ \mathrm{W}. \tag{8}\]

Thus \(\eta\simeq 19.6\%\) and \(P\simeq 693\ \mathrm{W}\), or \(0.69\ \mathrm{kW}\) to two significant figures. Only \(P/P_{0}=13.9\%\) of the heater input becomes work.

Two temperatures for the same power

Use the temperature rise as the unknown

Put \(x=T-T_{0}\). A difference in kelvin has the same numerical value as a difference in degrees Celsius. Since \(T=T_{0}+x\), the power relation becomes

\[ P(x)=\frac{x}{T_{0}+x}\left( P_{0}-Gx \right). \tag{9}\]

For the engine to deliver positive power, it must receive heat and be hotter than its cold reservoir:

\[ 0<x<\frac{P_{0}}{G}=250\ \mathrm{K}. \tag{10}\]

The two endpoints give zero output for different reasons: at \(x=0\), the Carnot efficiency is zero; at \(x=250\ \mathrm{K}\), all heater power leaks through the walls.

Solve part b without discarding a root

For a specified output \(P\), rearrangement gives

\[ P\left( T_{0}+x \right)=P_{0}x-Gx^{2}\quad\Rightarrow \quad Gx^{2}-\left( P_{0}-P \right)x+PT_{0}=0. \tag{11}\]

With \(P=500\ \mathrm{W}\), the quadratic formula gives

\[ x=\frac{4500\pm \sqrt{4500^{2}-4(20)(500)(300)}}{2(20)}\ \mathrm{K}=\frac{225\pm 25\sqrt{33}}{2}\ \mathrm{K}. \tag{12}\]

Thus \(x_{1}=40.69297\ \mathrm{K}\) and \(x_{2}=184.30703\ \mathrm{K}\). Add the outside Celsius temperature, \(27^{\circ }\mathrm{C}\), to these temperature rises:

\[ \theta_{1}={67.69297}^{\circ }\mathrm{C},\quad\quad\theta_{2}={211.30703}^{\circ }\mathrm{C}. \tag{13}\]

Both answers are admissible: both rises lie between \(0\) and \(250\ \mathrm{K}\), so each state has positive engine heat input and positive efficiency.

Why the hotter state gives the same output

QuantityCooler stateHotter state
Boiler temperature\(67.69^{\circ}\mathrm{C}\)\(211.31^{\circ}\mathrm{C}\)
Direct wall loss \(P_{\ell}\)\(813.86\,\mathrm{W}\)\(3686.14\,\mathrm{W}\)
Engine heat input \(P_h\)\(4186.14\,\mathrm{W}\)\(1313.86\,\mathrm{W}\)
Carnot efficiency \(\eta=P/P_h\)\(11.94\%\)\(38.06\%\)
Useful mechanical power \(P\)\(500\,\mathrm{W}\)\(500\,\mathrm{W}\)

The cooler engine converts a smaller fraction of a larger heat input into work. The hotter engine converts a larger fraction of a smaller heat input. These effects give the same output at two different temperatures. The question concerns adjustable steady operating states; their response to a disturbance would require a separate control model.

The greatest mechanical power

Find the peak of the power function

Define the limiting boiler temperature \(T_{\mathrm{lim}}\) as the temperature at which wall leakage alone consumes the entire heater input:

\[ T_{\mathrm{lim}}=T_{0}+\frac{P_{0}}{G}=300+250=550\ \mathrm{K}. \tag{14}\]

Then \(P_{h}=G\left( T_{\mathrm{lim}}-T \right)\), so the useful power can be written as

\[ P=G\left( T_{\mathrm{lim}}+T_{0}-T-\frac{T_{0}T_{\mathrm{lim}}}{T} \right). \tag{15}\]

This form separates two competing temperature terms. Differentiating with respect to \(T\) gives

\[ \frac{dP}{dT}=G\left( \frac{T_{0}T_{\mathrm{lim}}}{T^{2}}-1 \right). \tag{16}\]

The derivative is positive below \(\sqrt{T_{0}T_{\mathrm{lim}}}\) and negative above it. Therefore there is one maximum within the allowed interval \(T_{0}<T<T_{\mathrm{lim}}\):

\[ T_{*}=\sqrt{T_{0}T_{\mathrm{lim}}}=\sqrt{(300)(550)}=406.20192\ \mathrm{K}. \tag{17}\]
\[ \theta_{*}=406.20192-273={133.20192}^{\circ }\mathrm{C}. \tag{18}\]

Substituting this temperature gives a useful general result:

\[ P_{\mathrm{max}}=G\left( \sqrt{T_{\mathrm{lim}}}-\sqrt{T_{0}} \right)^{2}=751.92319\ \mathrm{W}. \tag{19}\]

Part c: the greatest output is about \(752\ \mathrm{W}\), at a boiler temperature of about \({133.2}^{\circ }\mathrm{C}\). This is well below the temperature \(277^{\circ }\mathrm{C}\) at which the wall loss would consume all \(5.0\ \mathrm{kW}\).

An algebraic method using the two roots

The official solution uses the quadratic from part b. As the requested output rises, its two physical temperature roots approach one another. At the maximum they coincide, so the discriminant is zero:

\[ \left( P_{0}-P \right)^{2}-4GPT_{0}=0. \tag{20}\]

Solving this quadratic in \(P\) gives

\[ P=P_{0}+2GT_{0}\pm \sqrt{\left( P_{0}+2GT_{0} \right)^{2}-P_{0}^{2}}. \tag{21}\]

The plus sign gives a power larger than the heater input and a negative value of the double temperature-rise root \(x=\left( P_{0}-P \right)/(2G)\). It is physically inadmissible. The minus sign gives \(751.92319\ \mathrm{W}\), agreeing with the calculus result.

The maximum can also be checked without solving a quadratic: for positive \(T\), the arithmetic-geometric mean inequality gives \(T+T_{0}T_{\mathrm{lim}}/T\geq 2\sqrt{T_{0}T_{\mathrm{lim}}}\), with equality precisely at \(T=T_{*}\).

Physical checks and conclusions

Mechanical power rises from zero at 27 degrees Celsius to a peak near 752 watts at 133.2 degrees Celsius, then falls to zero at 277 degrees Celsius; 500 watts occurs at two temperatures.
Figure 3. Both intersections with the horizontal line give an output of 500 watts. The highest point gives about 752 watts at 133.2 degrees Celsius. The curve is plotted only over the physical temperature interval.

Check where the full heater power goes

At the maximum, the temperature rise is \(106.20192\ \mathrm{K}\). Hence

\[ P_{\ell }=2124.0384\ \mathrm{W},\quad\quad P_{h}=2875.9616\ \mathrm{W}. \tag{22}\]

The engine rejects \(P_{c}=P_{h}-P_{\mathrm{max}}=2124.0384\ \mathrm{W}\). Consequently the complete energy balance is

\[ P_0=P_{\ell}+P_{\mathrm{max}}+P_c. \tag{23}\]
\[ (2124.0384+751.9232+2124.0384)\,\mathrm{W}=5000.0000\,\mathrm{W}. \tag{24}\]

Equal wall leakage and engine heat rejection emerge at the optimum; this was not assumed. The Carnot efficiency is \(26.1\%\), while \(15.0\%\) of the heater input becomes work.

Answers at a glance

PartResult
a\(\eta\simeq19.6\%\) and \(P\simeq693\,\mathrm{W}\)
b\(\theta\simeq67.7^{\circ}\mathrm{C}\) or \(211.3^{\circ}\mathrm{C}\)
c\(P_{\mathrm{max}}\simeq752\,\mathrm{W}\) at \(\theta_*\simeq133.2^{\circ}\mathrm{C}\)

Points to remember

The five loss faces provide parallel conduction paths.

Use kelvin in temperature ratios; kelvin and Celsius differences are numerically equal.

Apply Carnot efficiency to \(P_{h}=P_{0}-P_{\ell }\), and check each root’s physical range.

Increasing efficiency can accompany decreasing useful power.

Sources and further reading

The problem is adapted from the 59th Slovenian High-School Physics Competition, national round, 24 April 2021, Group III, Problem 4. The numerical data are retained, and the explanation and all three diagrams are independently prepared. DMFA question and official solution, PDF pages 7 and 21.

Source check: The question gives \(d=1.0\ \mathrm{cm}\). The official solution header prints \(1.0\ \mathrm{dm}\), but its calculations use centimetres. This handout follows the question and rounds only at the end. Extra digits support checking; the final results agree with the official rounded answers.