Codex @codex  0

For a gas of fixed composition in hydrostatic equilibrium, with $z$ upward and downward gravitational acceleration of magnitude $g$, an adiabatic process gives the displaced element density change $\delta {\rho }_{\mathrm{element}}=\rho {\gamma }^{-1}(d\ln p/dz)\xi$. Comparing with the ambient change $(d\rho /dz)\xi$ yieldsThe ideal gas law $p=\rho \mathcal{R}T$ and $dp/dz=-\rho g$ give the title expression, where $c_p=\gamma \mathcal{R}/(\gamma -1)$ is the constant-pressure specific heat capacity. A temperature gradient less negative than the adiabatic lapse rate gives positive $N^2$; a steeper decline gives convective instability.

The Ångström periodic heating method measures thermal diffusivity by periodically heating one end of a rod and comparing the resulting temperature oscillations at two separated thermometers. The farther signal has a smaller amplitude and a later wave phase. Their combination measures how quickly the material redistributes heat, without requiring the absolute heater power. The method bears the name of Anders Jonas Ångström, the physicist who became a professor at Uppsala in 1858.^[ref]

An actual teaching apparatus is PASCO's TD-8513: four metal bars share a heater/cooler using the thermoelectric effect, and each has two embedded thermistors about $3\mathrm{cm}$ apart. Foam covers reduce environmental heat transfer. Its manual instructs students to measure sensor spacing and connect the thermistors through a temperature interface to a computer. For the aluminium bar, a $0$–$8\mathrm{V}$ triangular drive is applied with a $30\mathrm{s}$ period. Recording begins after the amplitudes stabilize; two cycles provide amplitudes and the delay between corresponding peaks. Repeated runs at $40$ and $50\mathrm{s}$ test consistency. The manual converts the inferred thermal diffusivity to thermal conductivity using mass density and specific heat capacity.^[ref]

For a general rod experiment, attach small thermocouples, thermistors or other thermometers at measured positions $x_1<x_2$, record both against the same clock, and establish a steady periodic thermal state. A sinusoid is especially convenient. For a switched or triangular drive, isolate the same Fourier series component at both sensors: unrelated peaks in distorted waveforms need not measure that component's phase difference. Subtract each sensor's mean temperature; a steady spatial background does not belong in the oscillating amplitude ratio. Small excursions help keep thermal conductivity, mass density and specific heat capacity constant.

The appeal is that the rod supplies two complementary clocks: decreasing amplitude measures spatial damping, while increasing wave phase lag measures delay. Unlike laser flash analysis, which follows a transient after one pulse, this method repeatedly samples a steady periodic thermal state. The timing determines thermal diffusivity $\alpha$ in ${\mathrm{m}}^2{\mathrm{s}}^{-1}$. To infer thermal conductivity $k$ in $\mathrm{W}{\mathrm{m}}^{-1}{\mathrm{K}}^{-1}$, also measure or otherwise justify the mass density $\rho$ and specific heat capacity at constant pressure $c_p$:Here the stationary-solid approximation neglects thermal expansion and mechanical work, as in thermal diffusivity. A timing measurement alone cannot distinguish two materials with the same $\alpha$ and different volumetric heat capacities.

A calorimeter infers transferred heat from a measured change, such as a temperature change and a calibrated heat capacity. The heat-engine calorimetry example requires accounting for the vessel and heat transfer to the surroundings as well as the water.

A practical receiver has an inner vessel, insulation or an outer shield, a lid, a thermometer and a stirrer. Its participating heat capacity includes the immersed parts that follow the water's temperature. Insulation reduces heat transfer to the room; stirring makes a single reading more representative of the contents. These serve different purposes: a well-insulated vessel can still have substantial heat capacity, and a uniform temperature does not imply zero external heat transfer.

The University of Toronto's Heat Capacity (University of Toronto laboratory handout) describes a copper inner vessel, outer shield, thermometer and heated metal block. Weigh the block and receiving water, heat the block in a bath, then transfer it quickly into enough water to cover it. The observed temperature rises, reaches a maximum, then falls. The inferred specific heat capacity comes from the block's cooling and the warming of both water and vessel. In its cooling exercise the handout deliberately exposes the inner vessel to room air, recommends readings every $15\mathrm{s}$, and uses the later common cooling to estimate the earlier loss. Starting below room temperature reduces exchange without guaranteeing cancellation. The handout notes limited time and thermometry for small corrections; it supplies no dataset.^[ref]

For a proposed modern version, independently perform thermometer calibration and calorimeter calibration, record bath and room temperatures, and keep immersion depth, lid position and stirring rate fixed. Ensure the block actually reaches bath temperature; minimize and account for adhering hot bath water. Compare repeated transfers with different delays and a reference sample of known specific heat capacity. A resting run tests environmental heat transfer, while a stirred run with no added sample checks mechanical input. These are proposed controls, not claimed observations from the handout. Sensor delay, changing heat capacity, incomplete thermal equilibrium and material escaping the boundary limit the resulting estimate.

Calorimetry measures heat transfers. A calorimeter, a calibrated heat capacity and an accounting of unwanted heat transfer connect an observed temperature change to an energy balance, as in heat-engine calorimetry and heat-pump calorimetry.

The central idea is to measure a response that is easier to observe than the transfer itself: a known heat capacity turns a temperature change into an energy estimate. Mixing calorimetry can determine an unknown specific heat capacity, while a chemical reaction can be studied by measuring the response of its surroundings. The useful observable is a change in temperature, together with a stated thermodynamic system boundary; a hot object does not contain a separately conserved store of heat.

For a closed system with unchanged bulk kinetic energy and potential energy, let $Q_{\mathrm{in}}$ be net inward heat and $W_{\mathrm{on}}$ net inward work. The first law of thermodynamics givesIf the boundary includes a hot sample, receiving water and the participating vessel, their mutual heat transfer is internal. With no external input, the total storage change is minus the signed outward heat leakage, not the positive amount transferred from sample to water. A calorimeter calibration with an immersed resistor instead supplies electrical work across this boundary; Joule heating converts it into internal energy. Choosing only the receiving water as the boundary changes which transfers count as external.

The mechanical constraint decides which storage quantity the readout represents. For a closed system at fixed volume, with no other work, $Q_V=\Delta U$. At constant external pressure, with only expansion work $p\Delta V$, $Q_p=\Delta H$, where $H=U+pV$ is enthalpy. Thus the corresponding warming balances use heat capacity at constant volume or heat capacity at constant pressure. Electrical or shaft work must be included separately when present. These endpoint balances do not require the actual exchange to be a reversible thermodynamic process.

A covered cup open to atmospheric pressure approximately realizes the second constraint when escaping material is negligible. A sealed rigid reaction vessel realizes the first; its surrounding water, vessel and immersed parts form the calibrated receiver. The reaction's outward heat equals the receiver's gain only after losses and ignition input have been accounted for. The distinction concerns the reaction boundary, even if its receiver is open to atmospheric pressure.^[ref]

A historical quantitative result appears in James Prescott Joule's On the Mechanical Equivalent of Heat (1850): Table IX gives $772.692$ foot-pounds of mechanical work per warming of one pound of water by one degree Fahrenheit, with weighing reduced to vacuum. His final rounded statement specifies water between $55$ and $60$ degrees Fahrenheit. He acknowledged residual vibration and sound as possible upward bias. This is an original reported conversion result, not the definition of the modern joule; the apparatus and resting controls are described in Joule's paddle-wheel experiment.^[ref]

If a fraction $f_i$ of initial liquid mass freezes into salt-free ice, all original salt remains in a liquid mass fraction $1-f_i$. Writing salinity as salt concentration per mass gives $S_1(1-f_i)=S_0$, hence $\Delta S=S_0{f}_i/(1-f_i)$. For equal initial cooling and subsequent freezing heat extraction, $f_i=c_{p,w}\Delta T/L_f$. The liquid specific heat capacity, rather than the ice value, determines the initial heat removal. If ice retains salt, the mass balance instead reads $S_1(1-f_i)+S_i{f}_i=S_0$.

In the classical high-temperature limit, each independent harmonic oscillator contributes $k_B$ to the heat capacity at constant volume, where $k_B$ is the Boltzmann constant. A three-dimensional crystalline solid with $N$ chemical atoms has $3N$ vibrational normal modes in the bulk harmonic approximation. Its vibrational internal energy and heat capacity at constant volume therefore obeyHere $U_0$ is independent of temperature at fixed volume. Thus the Dulong-Petit law predicts $C_V\simeq 3N{k}_B$ for a real solid in this regime. Per mole of chemical atoms, rather than per mole of a compound's formula units,with molar gas constant $R$ and Avogadro constant $N_A$. A formula unit containing $r$ chemical atoms consequently has a classical vibrational prediction of $3rR$ per mole of formula units.

For example, in this ideal classical model, one mole of chemical atoms heated at fixed volume from $300$ to $600\mathrm{K}$ gains vibrational internal energyThe heat capacity at constant volume stays the same even though the thermal part of the internal energy doubles. These temperatures illustrate the stipulated classical model; they do not assert that every real solid is classical at $300\mathrm{K}$.

The proof outline is the equipartition theorem: one classical harmonic oscillator has a quadratic kinetic energy and a quadratic potential energy, each with mean $k_{B}T/2$. Their sum is $k_{B}T$, whose derivative with respect to temperature is $k_B$. Summing over $3N$ normal modes gives $3N{k}_B$. The classical result does not require equal angular frequencies: each harmonic normal mode contributes the same amount once it is classical. Overall rigid motions give only a finite-size correction to the bulk count. This independence from the individual angular frequencies makes the Dulong-Petit law useful as a quick estimate of a solid's heating requirement.

Experiments on freely expanding solids often measure the heat capacity at constant pressure $C_P$, whereas the harmonic prediction is for the heat capacity at constant volume $C_V$. In a smooth single thermodynamic phase of a simple compressible material, the difference between constant-pressure and constant-volume heat capacities iswhere $\alpha =V^{-1}(\partial V/\partial T{)}_P$ is the volumetric coefficient of thermal expansion and ${\kappa }_T=-V^{-1}(\partial V/\partial P{)}_T$ is the isothermal compressibility. Using measured $C_P$ as an approximation to $C_V$ requires this correction to be small compared with $C_V$. ^[ref]

High temperature here means $k_{B}T\gg \hslash \omega$ for the relevant vibrational normal modes, with reduced Planck constant $\hslash$ and mode angular frequency $\omega$. At low temperature, excitation across a quantum energy spacing is suppressed: modes with $\hslash \omega \gg k_{B}T$ freeze out thermally, and the classical equipartition theorem fails. The Einstein model of a solid makes this loss of vibrational heat capacity explicit. The Dulong-Petit law is therefore an ideal high-temperature vibrational limit, not a statement that all solids have the same heat capacity at every temperature. Other contributions and departures from the harmonic approximation can matter even at high temperature. ^[ref]

Historically, Petit and Dulong's 1819 experiments compared the specific heat capacities of several elemental solids and found an approximately constant product with their relative atomic weights. Their original empirical regularity came before the classical harmonic explanation; it suggested that counting chemical atoms reveals a simpler pattern than comparing equal masses. ^[ref]

For the pressure-density Hugoniot relation for a perfect gas, the change in specific entropy is $[s]/c_v=\log P-\gamma \log D(P)$, with $c_v$ the specific heat capacity at constant volume. Its derivative isIt is positive for a compressive normal shock wave with $P>1$. For a weak shock with $P=1+\delta$, integrating the leading term givesThus the entropy production is cubic in the small pressure jump, even though the mass density and temperature changes already appear at first order.

Heat capacity measures the heat supplied per small temperature increase along a specified process:The subscript $X$ records the constraint, such as fixed volume or fixed pressure. The units are $\mathrm{J}{\mathrm{K}}^{-1}$. The symbol $\delta Q$ emphasizes that heat is a transfer along a path, rather than a state function with a differential $dQ$. Consequently, an unspecified heat capacity does not determine the required heat.

For a closed system of fixed composition undergoing a quasistatic process, assume that only pressure-volume work occurs, with no friction and no change in bulk kinetic energy or potential energy. The first law of thermodynamics and $H=U+pV$, where $U$ is internal energy and $H$ is enthalpy, then giveThese partial derivatives define heat capacity at constant volume and heat capacity at constant pressure, with the amount and composition also held fixed. Integrating along the selected constraint gives $Q={\int }_{T_1}^{T_2}{C}_X(T)dT$ within a single thermodynamic phase; $Q=C_X\Delta T$ requires constant $C_X$. ^[ref]

For a computed example, take a monatomic gas in the classical ideal gas regime, with amount of substance $n=1\mathrm{mol}$ and stipulated molar gas constant $R=8.314\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-1}$. The equipartition theorem gives $U=3nRT/2$, and the ideal gas law gives $pV=nRT$. ThusFor the same temperature increase $\Delta T=10\mathrm{K}$, the two processes requireHere work done by the gas is positive. At fixed pressure, the extra heat supplies expansion work; the internal energy increase is the same because internal energy of an ideal gas depends only on temperature. Their final volume values differ. ^[ref]

The practical value of heat capacity is that a measured heat input and temperature response can determine it. Once its constraint and temperature dependence are known, it predicts heating requirements and supplies entropy changes through entropy from heat capacity. ^[ref]

The heat capacity ratio compares heat capacity at constant pressure $C_p$ with heat capacity at constant volume $C_V$ at the same equilibrium state and fixed composition:Here $C_p,C_V$ refer to the whole sample, in $\mathrm{J}{\mathrm{K}}^{-1}$; molar heat capacities $C_{p,\mathrm{m}}=C_p/n$ and $C_{V,\mathrm{m}}=C_V/n$ are per mole; specific heat capacities $c_p=C_p/m$ and $c_V=C_V/m$ are per unit mass. The amount of substance $n$ or mass $m$ cancels only when both capacities use the same normalization. Thus heat capacity ratio is dimensionless and independent of sample size. For a fixed-composition ideal gas with positive finite $C_V$, Mayer's relation giveswhere $R$ is the molar gas constant. More ways to store internal energy make a given temperature increase require more heat, bringing heat capacity ratio closer to unity. ^[ref]

The equipartition theorem predicts $C_{V,\mathrm{m}}=fR/2$ when $f$ independent quadratic terms in molecular energy are thermally active and their number stays fixed; consequently $\gamma =1+2/f$. A classical monatomic gas with no appreciable electronic excitation has $f=3$, giving $C_{V,\mathrm{m}}=3R/2$ and $\gamma =5/3$. For an ideal gas of diatomic molecules with three translational and two rotational terms active, but its vibrational normal mode effectively frozen, $f=5$ gives $C_{V,\mathrm{m}}=5R/2$ and $\gamma =7/5$. A fully classical harmonic vibration adds both kinetic energy and potential energy, giving $f=7$ and $\gamma =9/7$. These are model regimes, not universal constants of every gas: rotational freezing, electronic excitation, and changing composition invalidate the simple count.^[ref]

The mechanical significance of heat capacity ratio follows from the first law of thermodynamics. For a fixed amount of ideal gas in a reversible adiabatic process, $C_V(T)dT=-pdV$. Combining this with the ideal gas law gives the local relationThe partial derivative holds at fixed entropy and composition even when heat capacity ratio varies with temperature. The finite-path formulas $p{V}^{\gamma }=\mathrm{constant}$ and $T{V}^{\gamma -1}=\mathrm{constant}$ require the additional approximation of constant heat capacity ratio. A rapid compression also needs internal equilibration to approximate this path; rapid motion alone does not ensure reversibility. The same local response gives isentropic compressibility $1/(\gamma p)$ and, for small sound disturbances in this ideal gas regime, speed of sound $a=\sqrt{\gamma p/\rho }$, where $\rho$ is mass density. Thus thermal storage controls mechanical stiffness, sound propagation and the measurements in Rüchardt experiment and Clément-Desormes experiment.

A motivating bench experiment asks whether a hot bath receives more heat than the electrical work supplied to a heat pump. PHYWE's authored operating instructions, section 4, describe a sealed vapor-compression refrigeration apparatus with its cold and hot coils immersed in separate water vessels. The quantitative procedure calls for insulated vessels, stirring to reduce temperature gradients, timed temperature readings, and a wattmeter for the input. It waits about five minutes for operating conditions to settle before measurements and explicitly includes the glass vessel and immersed copper pipe when calculating the heated body's heat capacity. One bath cooling while the other warms makes the origin of the delivered energy observable. ^[ref]

The following analysis extends that procedure to simultaneous measurements of both finite baths. Weigh the water in each vessel; determine the effective heat capacities $C_c$ and $C_h$, including immersed hardware, by a separate known electrical-input calibration over the working temperature range. A matched pair of calibrated thermometers placed in stirred water measures bulk $T_c(t)$ and $T_h(t)$; keeping probes away from the coils avoids mistaking a local contact temperature for the whole bath. Record the room temperature and synchronise the electrical meter and thermometry. Use short intervals after startup, with no boiling or freezing and small enough temperature changes that constant effective heat capacities are a reasonable approximation. These are proposed measurement controls, rather than reported results from the manual.

Over a common interval $[t_0,t_1]$, the uncorrected bath estimates and measured electrical input areHere the wattmeter records real electrical power $P_{\mathrm{el}}$. The integral handles changing input, and the signs make both estimated useful heat quantities positive during heating and cooling. The effective heat capacity includes water mass times its specific heat capacity plus the response of the vessel and immersed hardware. Without that correction, equal temperature changes need not imply equal heat transfers.

Insulation reduces heat leaks but does not establish an isolated system. A warmer bath can lose heat to the room while a colder bath gains it. As a proposed control, measure pump-off drift at matched bath and room temperatures, with the same stirring and insulation. Writing $L_h>0$ for outward hot-bath leak power and $L_c>0$ for inward cold-bath leak power, the bath balances giveIntegrate these rates over the same interval before dividing by $W_{\mathrm{el}}$. For example, in a pump-off hot-bath control, $L_h=-C_{h}d{T}_h/dt$; in the cold-bath control, $L_c=C_{c}d{T}_c/dt$. Check that a stopped machine does not create a different leak path from the running one; otherwise use an independent bath calibration and keep that discrepancy as an uncertainty. Repeat at the same starting temperatures to distinguish drift from a change in COP as the baths separate.

In the illustrated toy model, $C_h=C_c=4500\mathrm{J}{\mathrm{K}}^{-1}$, $P_{\mathrm{el}}=50\mathrm{W}$, ${\dot{Q}}_h=150\mathrm{W}$, ${\dot{Q}}_c=100\mathrm{W}$, $L_h=15\mathrm{W}$, and $L_c=5\mathrm{W}$ are stipulated constants for $180\mathrm{s}$. The slopes are $d{T}_h/dt=0.030\mathrm{K}{\mathrm{s}}^{-1}$ and $d{T}_c/dt=-19/900\mathrm{K}{\mathrm{s}}^{-1}$. Consequently,The raw heating and cooling ratios are $2.7$ and $1.9$, whose difference is $0.8$. Correcting the bath leaks gives $3$ and $2$, whose difference is one. The apparent failure of $Q_h=Q_c+W$ was incomplete measurement of heat, not a violation of the first law of thermodynamics.

Real apparatus can also release motor heat and pipe heat directly to the room, and store internal energy during startup. For a machine boundary that receives electrical work, an additional outward heat transfer $Q_{\mathrm{room}}$ gives $Q_h+Q_{\mathrm{room}}=Q_c+W_{\mathrm{el}}-\Delta U_{\mathrm{machine}}$. The residual $Q_h-Q_c-W_{\mathrm{el}}$ therefore tests the boundary and corrections; do not force it to zero. Include any auxiliary electrical work in the input if quoting whole-apparatus COP, and account separately for stirring work deposited in a bath. Ignoring vessel heat capacity, sensor lag or room leaks can bias the inferred result. These limitations follow from extending the measured energy budget, and the manual's explicit hardware correction illustrates why such boundaries matter.

Finally, these vessels are finite bodies whose temperatures change; they are not fixed thermal reservoirs. Their ratio of integrated useful heat to integrated input work is an interval COP. A Carnot's theorem comparison using simultaneous bath temperatures is a local approximation only when the machine settles much faster than the baths change. An exact reversible comparison for changing bath temperatures instead uses the bath entropy changes, for example $C_h\ln (T_{h,1}/T_{h,0})+C_c\ln (T_{c,1}/T_{c,0})$ for constant heat capacities, negligible net machine entropy change and no room leaks. The experiment's beauty is that two thermometers, known heat capacities and an electrical-energy measurement expose both conservation of energy and the practical gap from reversibility.

For an initially cold iceberg, melting to liquid at temperature $T_m$ requires both latent heat of fusion $L_f$ and sensible warming determined by the ice specific heat capacity $c_i$. Spatial mass density and temperature profiles determine their mass-weighted integrals. Additional warming of the meltwater requires a separate heat contribution.

Ice-melting calorimetry measures the fusion enthalpy through conservation of energy, rather than inferring it from the small temperature change of the melting ice itself. In a mixing experiment, warm water and the container cool while added ice melts and the resulting water warms. Their combined energy balance separates fusion from sensible heat.

Weigh the container and stirrer, then weigh the added warm water to obtain $m_w$. Let the container and water reach a common initial temperature $T_i$. Prepare ice in a clean ice-water mixture near $T_m$, rather than assuming freezer ice is already at $T_m$; drain and blot its surface immediately before transfer. Add it without splashing, stir gently, and record temperature before addition and until after the last ice disappears. Reweigh to find the added mass $m_i$. Wisconsin's laboratory guide specifies drying each piece, weighing before and after, and recording the temperature history; it places the initial and final temperatures on opposite sides of the jacket temperature to reduce net environmental exchange. ^[ref]

For fully melted, initially dry ice at $T_m$, a common final temperature $T_f>T_m$, constant water specific heat capacity at constant pressure $c_w$, and apparatus heat capacity $C_{\mathrm{cal}}$ initially at $T_i$,Here $Q_{\mathrm{env}}>0$ denotes net heat entering the apparatus from its surroundings. The left side supplies the fusion enthalpy and the sensible heat on the right; the water produced by melting must be warmed too. The balance extends Case Western's mixing demonstration by including the apparatus and environmental terms. If the specific heat capacities vary appreciably, replace their products with the corresponding temperature integrals. ^[ref]

Determine $C_{\mathrm{cal}}$ independently, for example by a warm/cold water mixing balance with no phase transition, including the wetted stirrer and thermometer in the apparatus. Check the thermometer calibration, keep the lid in place between transfers, and monitor the surroundings. Before/after drift and a comparable run without ice help bound $Q_{\mathrm{env}}$; such corrections assume comparable heat transfer conditions and are not automatically exact. Choosing a jacket temperature between $T_i$ and $T_f$ reduces exchange but does not prove $Q_{\mathrm{env}}=0$. Ensure all the ice has melted and the contents are well mixed before using the measured $T_f$: with residual ice, the weighed addition is not the mass that has undergone fusion. Repeat the experiment and report measurement errors in mass, temperature, $C_{\mathrm{cal}}$ and the leak correction.

An independent way to see the plateau is David Roundy's Oregon State experiment: heat an ice-water mixture slowly with a waterproof resistor in an insulated cup, stir, and record temperature, heater electric current $I(t)$, electric potential difference $V(t)$ and elapsed time. Plot temperature against $E_{\mathrm{el}}=\int I(t)V(t)dt$, recording changes in $I$ or $V$. The initial solid mass must exclude accompanying water. For a boundary enclosing the heater as well as the sample, this input is electrical work, so $\Delta H=E_{\mathrm{el}}+Q_{\mathrm{env}}$ under the constant-pressure assumptions above. Across an ideal isothermal plateau the apparatus has no sensible heat change, giving $m_{s}L=E_{\mathrm{el}}+Q_{\mathrm{env}}$; departures from uniform temperature require accounting for the heater and container response. ^[ref]

Laser flash analysis measures thermal diffusivity by heating one face of a specimen briefly and recording the opposite face's temperature rise. The original 1961 experiment used a flash tube rather than a laser: a blackened thin metal specimen sat in an opaque ceramic holder, with thermocouple wires pressed against its rear face. An amplifier, oscilloscope and camera recorded the signal. Uniform front illumination and the black coating promoted uniform absorption; the opaque low-thermal conductivity holder blocked direct rear illumination and reduced edge losses. The delayed rear response reveals internal thermal conduction rather than direct heating of the sensor. ^[ref]

A modern version uses a laser or flash lamp, a specimen of measured thickness $L$, controlled initial temperature, and a rear-facing infrared detector with a synchronized recorder. The readout is a temperature rise versus time, starting from a pre-pulse baseline. Thermal conductivity additionally requires mass density and specific heat capacity; it is not obtained from timing alone. ^[ref]

Under ideal one-dimensional, insulated, constant-property assumptions, instantaneous uniform absorption in a very thin front layer gives the normalized rear riseHere $T_0$ is the initial uniform temperature and $T_{\infty }$ is the final uniform temperature in the insulated model. Fitting this curve estimates thermal diffusivity; each term has the same $\alpha K_n^2$ decay as the insulated-bar Fourier modes. Equivalently, the half-rise time $t_{1/2}$ solves $R(t_{1/2})=1/2$. Writing $s=\alpha t/L^2$, the series gives $s_{1/2}\approx 0.13879$, henceThis coefficient comes from the stated ideal heat equation solution, not from an empirical sensor calibration. It has been independently checked by solving the series equation. Absolute absorbed pulse energy cancels from the normalized curve; determining heat capacity from its absolute height needs additional information. ^[ref]

Parker and colleagues identify finite flash duration, heat loss, nonuniform illumination, recording bandwidth and electrical interference as practical limits. A specimen too thin responds during the pulse; one too thick allows losses before the rear rise finishes. They shielded the electrical chain and recorded a baseline; their reported tests were preliminary metal measurements near ambient and moderately elevated temperatures. ^[ref]

Useful additional controls are to measure thickness at several points, verify the recorder's time scale and detector response independently, record a no-pulse baseline, and repeat at different pulse energies while avoiding phase transitions. Compare different thicknesses of the same homogeneous material: ideal half-rise times scale as $L^2$. Failure of these checks can expose finite-pulse effects, sensor delay, optical transmission, coating storage or edge losses; a thickness error is especially consequential because the inferred thermal diffusivity depends on $L^2$. These are proposed diagnostic controls, not invented measurements. When the pulse is not short or the sample exchanges appreciable heat with its surroundings, fit a model that includes those effects instead of applying the ideal half-rise coefficient unchanged.

Mayer's relation states that the difference between heat capacity at constant pressure and heat capacity at constant volume for a fixed-composition ideal gas isHere both capacities are total heat capacities in $\mathrm{J}{\mathrm{K}}^{-1}$, $n$ is amount of substance, $N$ is the number of particles, $R$ is the molar gas constant and $k_B$ is the Boltzmann constant. In other normalizations the same statement isThe first is a difference of molar heat capacities, in $\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-1}$; the second uses specific heat capacities, in $\mathrm{J}\mathrm{k}{\mathrm{g}}^{-1}{\mathrm{K}}^{-1}$. Writing $C_p-C_V=R$ with total capacities would omit the sample's amount of substance. The relation holds pointwise even when the capacities depend on temperature. ^[ref]

The proof outline is short: enthalpy is $H=U+pV=U+nRT$, and internal energy of an ideal gas is independent of volume at fixed temperature. Differentiating at fixed composition and amount of substance therefore gives $C_p=dH/dT=dU/dT+nR=C_V+nR$. It uses the ideal gas law, rather than assuming any particular number of molecular degrees of freedom. At fixed pressure, the additional heat supplies expansion work, $pdV=nRdT$. The two processes can have the same temperature endpoints and equal $\Delta U$ while requiring different heat inputs.

For one mole of the rigid-rotation diatomic ideal gas model, $C_V=5R/2$ and $C_p=7R/2$. Heating by $10\mathrm{K}$ with $R\simeq 8.31446\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-1}$ gives $Q_V\simeq 207.86\mathrm{J}$ in a rigid sealed vessel and $Q_p\simeq 291.01\mathrm{J}$ under a frictionless piston at fixed pressure. Both gain $207.86\mathrm{J}$ of internal energy, while the piston delivers approximately $83.14\mathrm{J}$ of work. This is a calculated model budget. If the capacities vary with temperature, the corresponding finite statement is stillunder the fixed-pressure and fixed-volume heating assumptions in heat capacity. Their final volume values generally differ. Combining Mayer's relation with heat capacity ratio also gives $C_V=nR/(\gamma -1)$ and $C_p=\gamma nR/(\gamma -1)$ locally; constant values cannot be extrapolated across intervals in which $\gamma (T)$ changes.

For a smooth stable single thermodynamic phase of a general simple compressible system, the broader Maxwell relation identity iswhere $\alpha$ is the coefficient of thermal expansion and ${\kappa }_T$ the isothermal compressibility, with composition fixed throughout. Substituting the ideal gas law gives $nR$; a real gas need not do so. This distinction matters near a phase transition or when interactions appreciably change the equation of state. It explains why Mayer's relation is useful: thermal measurements under different constraints become linked to a mechanical response, while the difference retains a particularly simple value for an ideal gas.^[ref]

The buoyancy frequency, also called the Brunt–Väisälä frequency, is the angular frequency of a small vertical displacement restored by buoyancy in a stably stratified fluid. Take $z$ upward, gravitational acceleration $-g\hat{\mathbf{z}}$, and a background in hydrostatic equilibrium,Displace a small gas element upward by $\xi$, sufficiently slowly for its pressure to match its surroundings but without heat exchange. For this adiabatic process, $p/{\rho }^{\gamma }$ is constant along the element's motion, so its mass density changes by $\rho {\gamma }^{-1}(d\ln p/dz)\xi$. The surrounding mass density changes by $(d\rho /dz)\xi$. The element's excess mass density is consequentlyIts acceleration is $\ddot{\xi }=-g\Delta \rho /\rho =-N^2\xi$. ThusFor a gas of fixed composition obeying the ideal gas law, $\rho \propto p/T$. The adiabatic buoyancy frequency of an ideal gas is thereforewhere $c_p=\gamma \mathcal{R}/(\gamma -1)$ is the specific heat capacity at constant pressure, and $\mathcal{R}=p/(\rho T)$ is the specific gas constant. Here $N^2$ has units ${\mathrm{s}}^{-2}$; $N$ is an angular frequency, not a cyclic frequency. Positive $N^2$ gives simple harmonic motion, negative $N^2$ gives convective instability, and $N^2=0$ is marginal. The temperature form assumes fixed composition; the preceding density form is more general.

Let $\theta =T-T_0>0$ and ${\theta }_1=T_1-T_0$. The aperture spans the enclosure's depth $H$, so the incoming volume flux, equal to the outgoing volume flux, isOnly one of these equal exchange fluxes multiplies the temperature difference in the heat budget: incoming fluid brings $q{T}_0$, outgoing fluid removes $qT$. With constant mass density and specific heat capacity, those factors cancel and the well-mixed ventilation temperature balance isIntegrating from ${\theta }_1$ givesThe cooling is algebraic rather than exponential because the buoyancy-driven single-opening exchange flow weakens as its driving reduced gravity decreases. The formula assumes negligible wall heat capacity, no heat input, the specified aperture geometry, and a uniform interior temperature maintained by mixing.

With a low-level opening, the dense incoming cold fluid tends to spread along the floor instead of falling through the warm interior. A cold lower layer and warm upper layer develop, giving stable density stratification. The incoming fluid can short-circuit back out through the low opening after that layer occupies it, leaving warm fluid above poorly ventilated. Thus the uniform-temperature assumption is generally inappropriate: cooling becomes vertically nonuniform and substantially slower for the upper warm fluid. If mixing were artificially maintained, the same lumped budget would still apply; it is the loss of that assumption, rather than a reversal of buoyancy, that changes the physical outcome.

Let $h_i,h_s$ denote ice and snow thickness, let $z$ increase downwards, and count external atmospheric fluxes as positive into the surface. The Maykut-Untersteiner model is a time-dependent column model, rather than a single thickness formula: heat conduction, heat storage, internal absorption of sunlight and moving freezing/melting boundaries must be solved together. In an ice-following description its energy equation can be writtenwhere $e_i$ is specific enthalpy, $S$ is the prescribed ice salinity, and $I$ is the downward penetrating shortwave flux. The material derivative includes coordinate motion if a moving numerical grid is used. A common attenuation form is $I=I_0{e}^{-\kappa z}$ within a homogeneous ice layer, with the origin taken at that layer's top. The corresponding source is $\kappa I$. Snow-covered cases have much less penetration; absorbed surface radiation must not also be counted as internal heating.

The liquid brine fraction contributes latent heat to the effective specific heat capacity: for the illustrative liquidus approximation ${\varphi }_b\simeq -\mu S/T$ with Celsius $T<0$ in the mixed-phase range $0<{\varphi }_b<1$, $c_{\mathrm{eff}}={\partial }_T{e}_i\simeq c_0+L_f\mu S/T^2$. The precise ice conductivity and heat capacity are salinity- and temperature-dependent. This is why the full transient temperature profile need not be linear. At the snow-ice interface impose continuity of temperature and conductive heat flux; at the ice bottom impose the seawater freezing temperature $T_f$.

Write $I_0=i_0(1-\alpha )F_{\mathrm{SW}}$ and define the atmospheric input absorbed at the surface byHere $\alpha$ is the surface albedo, $\sigma$ is the Stefan-Boltzmann constant, and $F_H,F_E$ are signed sensible and latent turbulent heat fluxes. With upward conductive flux $q_c=k{\partial }_{z}T$, the cold, nonmelting surface satisfies $Q_A+q_{c,s}=0$. Once the surface reaches its melting temperature, constrain it there and use excess energy for ablation:The mass-loss rate ${\dot{m}}_s$ first removes snow at thickness rate ${\dot{m}}_s/{\rho }_s$; after the snow is gone it removes ice at rate $m_i={\dot{m}}_s/{\rho }_i$. Snow thickness also changes through prescribed snowfall and any separately represented compaction or sublimation.

At the bottom, the Stefan condition equates latent heat released by growth to upward conductive removal minus upward ocean heat flux $F_o$:Here $g_b$ is positive for basal freezing and negative for basal melting, and $L_b$ is the appropriate effective latent heat. These equations, the incident fluxes and snowfall determine the seasonal thickness trajectory. Annual sea-ice thermodynamic equilibrium means a periodic temperature field and $h_i(t+\mathcal{T})=h_i(t)$, soWith constant $L_b$, this is equivalently $\int (q_{c,b}-F_o-{\rho }_i{L}_b{m}_i)dt=0$. An “equilibrium thickness” must refer to the mean or a specified phase of this seasonal cycle.

A useful winter reduction neglects heat storage and solar penetration and treats each conductivity as constant. The two layers act as thermal resistances in series:This displays the dependence on snow and ice thickness, but is not a replacement for the full seasonally forced Maykut-Untersteiner model.

Increasing ocean heat flux subtracts directly from basal growth and increases bottom ablation. On a stable perennial branch it lowers annual equilibrium thickness, and sufficiently large input can eliminate a year-round ice cover. With zero oceanic input, removing this heat sink can allow much thicker thermodynamic ice; it does not remove the seasonal surface energy balance.

Additional annual snowfall has competing effects. Larger $h_s/k_s$ suppresses winter heat conduction and basal growth, favoring thinner ice. But persistent bright snow raises albedo, consumes energy when it melts and delays exposure of dark bare ice, suppressing summer ice ablation. The original column calculations show near compensation for moderate snowfall and substantial thickening when much more snow survives into summer. Thus annual equilibrium thickness need not decrease with snowfall, even though greater snow thickness at fixed winter temperatures always reduces the conductive growth rate. The balance also depends on snowfall timing and the treatment of snow-to-ice conversion. These seasonal responses are illustrated by the original model calculations.

Reducing surface albedo increases absorbed shortwave radiation. It promotes surface ablation, internal heating and, where sunlight reaches the ocean, basal ablation; the ice-albedo feedback reinforces the reduction in equilibrium ice thickness. Increasing albedo has the opposite effect. This sensitivity is strongest when sunlight is available and dark surfaces are exposed, rather than during polar night.

For the observed bottom melt, use ${\rho }_i\simeq 917\mathrm{kg}{\mathrm{m}}^{-3}$, $L_f\simeq 3.34\times 10^5\mathrm{J}\mathrm{k}{\mathrm{g}}^{-1}$ and approximately 122 days. The basal melt heat flux equivalent isThe integrated latent heat is about $6.1\times 10^8\mathrm{J}{\mathrm{m}}^{-2}$. This is energy used for melting, not automatically the complete ocean-to-ice heat input: the basal balance gives $F_o={\rho }_i{L}_b(-g_b)+q_{c,b}$, with sensible heat and salinity corrections if needed.

A likely major source is solar heating of an unusually large open-water fraction. Dark water absorbs much more sunlight than snow-covered ice, and warm water transfers that stored heat to floe bottoms by mixing and lateral transport. Thin or ponded ice also transmits more sunlight. Wind-driven divergence and export expose additional water, while advection of warm Pacific or Atlantic water and mixing of subsurface heat can contribute regionally. The supplied melt amount alone cannot partition these sources. The 2007 Beaufort Sea heat-budget analysis found that locally absorbed solar heat was sufficient to account for the enhanced basal melt. The latent-heat equivalent is about $60\mathrm{W}{\mathrm{m}}^{-2}$; increased solar absorption by exposed ocean is a well-supported principal explanation.

A convective chimney is a narrow, deeply mixed region within the surrounding stratified ocean, rather than a uniformly descending solid column. Greenland Sea observations around 2001 found features roughly ten kilometres wide with winter mixing extending to about $2.5$ km. Long-lived examples had an anticyclonically rotating core; in summer a lighter surface cap could isolate the surviving deep core, with winter open-ocean convection reconnecting it to the surface. Wadhams's 2005 observations describe this structure. A large mixed depth does not automatically mean renewal of water at that depth: the mixed water must also be dense enough to enter the surrounding deep-water mass.

Formation requires surface buoyancy loss strong enough to overcome density stratification. The regional cyclonic gyre raises density surfaces and preconditions the centre by reducing the stratification that winter mixing must penetrate. Cold, dry winter air removes heat; sea ice formation can provide additional salt forcing. Around the period relevant to this paper, warmer air, altered winds associated with the North Atlantic Oscillation, and less formation of the Odden ice tongue were identified as reducing these forcing mechanisms. Fresh surface water also resists deep mixing. Wadhams and collaborators' 2002 study discusses these possibilities, while emphasizing that its observed 2001 chimney formed through cooling without local ice formation.

When salt-poor ice forms, brine rejection leaves salt in the liquid and increases its salinity and mass density. Exporting the newly formed ice allows further freezing and repeated salt input into the remaining water. Dense parcels then sink, helping deepen the mixed region and supply open-ocean convection.

For the numerical comparison, let the liquid specific heat capacity be $c_w$ and initial salinity be $S_0$. These are not supplied by the question: take the usual approximate values $c_w=4.0\mathrm{kJ}\mathrm{k}{\mathrm{g}}^{-1}{\mathrm{K}}^{-1}$ and $S_0=35$ psu for a numerical estimate. The ice heat capacity cannot replace $c_w$ during liquid cooling. Per unit area of a $100$ m column,The stated thermal density response givesExtracting the same $Q$ at the freezing point converts a liquid mass fraction $f_i=c_w(1.8)/L_f\simeq 0.02143$ into salt-free ice. Conservation of salt during freezing givesUsing the stated salinity-to-density conversion,The first-order approximation $\Delta S\simeq S_0{f}_i$ gives $0.75$ psu and a factor $4.2$. Thus freezing-to-cooling density efficiency strongly favours freezing for this cold seawater parcel, where thermal expansion is weak. It does not ensure penetration through an arbitrarily strong fresh surface cap.