2. Evaporation 7 pts

Thermodynamics · Phase transitions, Diffusion, Heat transfer

Piston–cylinder equilibrium with saturated vapour, and derivation of the wet-bulb temperature of sweating skin in a sauna.

Self-assessment by Claude Opus 4.7. 6.9 / 7.0

Part (i) — 2.0 / 2.0 pts

CriterionPointsResult
Realize that the pressure inside the cylinder equals the saturated vapour pressure of water at temperature T0T_00.8✓ “the gap therefore fills with saturated water vapour at T0=90°CT_0 = 90\,°\text{C}, at the pressure psat(T0)p_{\mathrm{sat}}(T_0)”, with an explicit aside that there is no air in the gap
Read the density ρ\rho from the graph, in the range [400,440]gm3[400,\,440]\,\text{g\,m}^{-3}0.2ρsat(90°C)0.42kgm3=420gm3\rho_{\mathrm{sat}}(90\,°\text{C}) \approx 0.42\,\text{kg\,m}^{-3} = 420\,\text{g\,m}^{-3}
Use the ideal gas law to find an expression for the pressure p1p_1 at temperature T0T_00.4psat=ρsatRT0/μ70kPap_{\mathrm{sat}} = \rho_{\mathrm{sat}} R T_0/\mu \approx 70\,\text{kPa}
Correct expression for the force: S(p0p1)S(p_0 - p_1)0.4Fmin=(p0psat(T0))SF_{\min} = (p_0 - p_{\mathrm{sat}}(T_0))\,S derived from a signed force balance on the piston
Correct numerical answer0.2Fmin300NF_{\min} \approx 300\,\text{N}

Part (ii) — 2.0 / 2.0 pts

CriterionPointsResult
Read the vapour density ρ1\rho_1 from the graph (ρ1[390,420]gm3\rho_1 \in [390,\,420]\,\text{g\,m}^{-3})0.2ρsat(89°C)0.40kgm3=400gm3\rho_{\mathrm{sat}}(89\,°\text{C}) \approx 0.40\,\text{kg\,m}^{-3} = 400\,\text{g\,m}^{-3}
Correct expression for mass of water vapour0.3mv=ρsat(T1)V=ρsat(T1)aSm_v = \rho_{\mathrm{sat}}(T_1)\,V = \rho_{\mathrm{sat}}(T_1)\,a S
Correct expression for the latent heat (mvLm_v L)0.3✓ “latent heat absorbed: mvLm_v L”, with explicit identification of LL as the enthalpy of vaporisation
Correct expression for heat lost by water ((mmv)c(T0T1)(m - m_v)\cdot c\cdot(T_0 - T_1))0.3(mmv)cΔT(m - m_v)\,c\,\Delta T written verbatim before the mvmm_v \ll m simplification
Expression for energy conservation0.4(mmv)cΔT=mvL(m - m_v)\,c\,\Delta T = m_v L
Correct expression for mass of water mm0.3mρsat(T1)aSL/(cΔT)m \approx \rho_{\mathrm{sat}}(T_1)\,a S L / (c\,\Delta T)
Correct numerical answer m[630,680]gm \in [630,\,680]\,\text{g} (with correct dimension)0.2m0.65kg=650gm \approx 0.65\,\text{kg} = 650\,\text{g}

Part (iii) — 2.9 / 3.0 pts

CriterionPointsResult
Heat going away from skin (up) = heat going to skin (down) at equilibrium0.4✓ “In: heat conducted down from the laminar layer, qq per unit area. Out: latent heat carried by evaporating water, LjL\,j per unit area. In steady state these must be equal”
Heat flux down κdT/dx\kappa\,dT/dx0.2q=κ(TTs)/dq = \kappa(T - T_s)/d from the linear T(x)T(x) profile
Heat flux up magnitude DLμRddxPT\frac{DL\mu}{R}\frac{d}{dx}\frac{P}{T} (or equivalent)0.5✓ Solution-2 form Lj=LDdρv/dxL\,j = L D\,d\rho_v/dx, derived via ρv=mH2On\rho_v = m_{\mathrm{H_2O}} n so the mass flux is j=Ddρv/dxj = D\,d\rho_v/dx; the per-NB-note ρ\rho-substitution subsumes the LJmLJm, m=μ/NAm=\mu/N_A and n=P/(TkB)n = P/(Tk_B) subitems
Direction of the heat flow opposite to dn/dxdn/dx (explicit or via the minus sign)0.1✓ “Heat flux (Fourier), magnitude, downward toward the skin” and “Going upward (water leaves the skin and reaches the drier bulk)” — the two fluxes called out as anti-parallel
P=rpP = rp (boundary condition at the bulk and at the skin)0.1✓ "ρv(d)=rρsat(T)\rho_v(d) = r\,\rho_{\mathrm{sat}}(T)" at the bulk and ”ρv(0)=ρsat(Ts)\rho_v(0) = \rho_{\mathrm{sat}}(T_s)… the wetness pins the local humidity to 100%100\%” at the skin
Integrating correctly to κ(TTs)=(DLμ/R)[p(Ts)/Tsrp(T)/T]\kappa(T - T_s) = (DL\mu/R)\,[\,p(T_s)/T_s - r\,p(T)/T\,] (or equivalent in ρ\rho)0.4✓ steady-state plus linear profiles give κ(TTs)=LD(ρsat(Ts)rρsat(T))\kappa(T-T_s) = LD\,(\rho_{\mathrm{sat}}(T_s) - r\,\rho_{\mathrm{sat}}(T)) — the ρ\rho-form is the post-substitution equivalent the NB note explicitly accepts
ρ=pμ/RT\rho = p\mu/RT0.1✓ used in part (i) at 90°C90\,°\text{C} and again in the 110°C110\,°\text{C} cross-check; the part-(iii) derivation works in ρ\rho from the start so the conversion is built into the boundary conditions
Reading ρ\rho correctly (ρ1[800,815]gm3\rho_1 \in [800,\,815]\,\text{g\,m}^{-3} at 110°C110\,°\text{C})0.2~ 0.1 / 0.2 — Claude reports ρsat(110°C)830gm3\rho_{\mathrm{sat}}(110\,°\text{C}) \approx 830\,\text{g\,m}^{-3}, marginally above the keyed band; the in-text ideal-gas cross-check at psat=143kPap_{\mathrm{sat}} = 143\,\text{kPa} would actually yield 808gm3808\,\text{g\,m}^{-3} (in band), so the 830 is a slightly aggressive chart extrapolation past the chart’s right edge at 105°C105\,°\text{C}
Graphical method — noticing ρ(Ts)=rρ(T)+(κ/DL)(TTs)\rho(T_s) = r\rho(T) + (\kappa/DL)(T - T_s) defines a straight line in (T,ρ)(T, \rho)0.8✓ explicit “straight line in TsT_s, of slope κ/(LD)-\kappa/(LD), passing through the point (Ts,ρ)=(T,ρ)(T_s, \rho) = (T, \rho_\infty) on the chart”; the numerical-tabulation root-find is one of the “any other valid numerical method” routes the scheme explicitly accepts
Correct final result Ts[36,47]°CT_s \in [36,\,47]\,°\text{C}0.2Ts43°CT_s \approx 43\,°\text{C}

Discrepancies

ρsat(110°C)\rho_{\mathrm{sat}}(110\,°\text{C}): Claude 830gm3830\,\text{g\,m}^{-3} vs official band [800,815]gm3[800, 815]\,\text{g\,m}^{-3}. The downstream impact is small: Claude’s intercept ρ=rρsat(T)=0.0383025gm3\rho_\infty = r\,\rho_{\mathrm{sat}}(T) = 0.03 \cdot 830 \approx 25\,\text{g\,m}^{-3} versus the official’s 24.3gm324.3\,\text{g\,m}^{-3} — a \sim2.5% difference at the right end of the line that swings by \sim34 g/m³ across the relevant TT-range. Final TsT_s: Claude 43°C43\,°\text{C} vs official 41.5°C41.5\,°\text{C}, both well inside the [36, 47] band; the 1.5°C gap is exactly what the intercept shift predicts.

Overall score: 6.9 / 7.0 pts

Full marks on parts (i) and (ii); a single 0.1 pt dock in part (iii) for an extrapolated chart read of ρsat(110°C)\rho_{\mathrm{sat}}(110\,°\text{C}) that is marginally above the keyed band.

Numerical answers match the official key on every part: Fmin300NF_{\min} \approx 300\,\text{N} (exact match), m650gm \approx 650\,\text{g} (exact match), Ts43°CT_s \approx 43\,°\text{C} vs official 41.5°C41.5\,°\text{C} (within the [36, 47] tolerance band, with the 1.5°C gap traceable to the marginally high ρ(110°C)\rho(110\,°\text{C}) extrapolation).

Commentary

Where this solution goes beyond the grading scheme. The “Overview” front-loads the three structural ideas — (1) liquid–vapour mechanical contact pins the gas pressure to psatp_{\mathrm{sat}}, (2) LL is enthalpy and already includes pΔVp\,\Delta V, (3) conduction and diffusion both scale as 1/d1/d so the layer thickness cancels — and quotes the small parameter 1psat/p00.31 - p_{\mathrm{sat}}/p_0 \approx 0.3 that controls part (i). Part (i) follows the boxed 300N300\,\text{N} with a “physics remark” contrasting the vaporisation route with the alternative of breaking a tensile water column (several MPa for clean water in clean tubing) and three consistency checks: the T0100°CT_0 \to 100\,°\text{C} boiling-point limit where Fmin0F_{\min} \to 0 (matching the everyday observation that boiling water lifts pot lids), the T00°CT_0 \to 0\,°\text{C} cold limit where Fminp0S1000NF_{\min} \to p_0 S \approx 1000\,\text{N} (almost as hard as pulling against a hard vacuum), and an explicit dimensional check. Part (ii) carries an explicit psatV200Jp_{\mathrm{sat}} V \approx 200\,\text{J} versus mvL2700Jm_v L \approx 2700\,\text{J} comparison that quantifies the \sim7% double-counting hazard for a student who treats LL as an internal-energy change, plus an a posteriori check mv/m2×103m_v/m \approx 2 \times 10^{-3} that justifies dropping mvm_v on the right-hand side, an order-of-magnitude estimate (“a litre of water cools by 1K1\,\text{K} when it loses 4200J\sim 4200\,\text{J}”), and a 4%4\% sensitivity calculation showing how using ρ(90°C)\rho(90\,°\text{C}) instead of ρ(89°C)\rho(89\,°\text{C}) would shift the answer to 0.68kg0.68\,\text{kg}. Part (iii) develops the wet-bulb equation pedagogically: first the explicit linear profiles for T(x)T(x) and ρv(x)\rho_v(x) across the laminar layer, then the energy balance with the careful note that the bulk humidity rr does not enter at x=0x = 0 because wetness pins the local humidity to 100%100\%. The “memorable check” — comparing dry sauna (r0r \to 0, line passes through (T,0)(T, 0), TsT_s drops further) with steam sauna (r1r \to 1, line passes through the saturation curve, TsT=110°CT_s \to T = 110\,°\text{C}, fatal) — frames the answer in everyday physical experience. The Lewis number is named (κ/(ρaircpairD)1\kappa/(\rho_{\mathrm{air}} c_p^{\mathrm{air}} D) \approx 1 for water vapour in air) and the answer is contextualised against meteorological wet-bulb data (dry desert at 40°C40\,°\text{C} has Twb20°CT_{wb} \approx 20\,°\text{C}). A plausibility check connects 43°C43\,°\text{C} to protein-denaturation onset at 60°C\sim 60\,°\text{C}, explaining why the experience is painful but survivable.

Where the official solution is sharper. Three places. (1) The official actually overlays the line on the saturation chart and reads the intersection from the printed graph (/figures/nbpho-2025/sol02-fig1.png), which is the most direct realisation of the “graphical method” worth 0.8 pts. Claude works the same equation through a tabulated numerical root-find with an ASCII sketch — equivalent in physical content and accepted by the scheme, but visually less convincing than a real overlay on the supplied chart. (2) The official gives a tighter, less aggressive read of ρsat(110°C)810gm3\rho_{\mathrm{sat}}(110\,°\text{C}) \approx 810\,\text{g\,m}^{-3} (matching exactly the steam-table-based ideal-gas computation psat(110°C)143kPap_{\mathrm{sat}}(110\,°\text{C}) \approx 143\,\text{kPa}, ρ=pμ/(RT)808gm3\rho = p\mu/(RT) \approx 808\,\text{g\,m}^{-3}), placing it inside the keyed [800, 815] band; Claude’s 830 sits 1.85% above the upper band edge, and his own in-text ideal-gas cross-check, if computed exactly, would have produced 808 — the inconsistency cost the 0.1 pt dock. (3) Part (ii): the official’s energy balance includes a one-line justification ”μL4R(T1T0)\mu L \gg 4R(T_1 - T_0) … [so we may] neglect the dependence of LL on temperature and the heat capacity of water vapours”, a clean inequality that shows in two seconds why both omissions are safe. Claude justifies dropping mvm_v from (mmv)(m - m_v) and warns against double-counting pΔVp\,\Delta V, but never explicitly checks whether L(T)L(T)-variation and the vapour heat capacity contribute at the percent level in the same way — the omission is implicit in his ”LL is enthalpy, full stop” framing, where the official spells out the inequality.