 Fourth-power thermal boundary model

Choose a positive constant $\Sigma$ and the illustrative interface lawHere $T_h$ and $T_c$ are the hot and cold interface-side temperatures, not distant thermometer readings. This is a low-temperature phonon model under restrictive assumptions, not a universal fit to Kapitza's measured exponent or a material-property table. Its coefficient $\Sigma$ is not the Stefan-Boltzmann constant: the carriers in this model are phonons, rather than photons crossing a vacuum.

The fourth power has a simple counting explanation. In the three-dimensional low-temperature Debye model, a linear phonon dispersion relation gives a density of states proportional to ${\omega }^2$. A phonon has energy $\hslash \omega$ and Bose-Einstein distribution occupation $[\exp (\hslash \omega /k_{B}T)-1{]}^{-1}$. If propagation speeds and the angle-averaged transmission introduce no additional frequency scale over the occupied modes, the outgoing energy flux has the formExtending the upper integration limit assumes occupied modes lie well below the cutoff. Balancing the two outgoing fluxes, with coefficients satisfying zero net heat flux at equal temperatures, gives the chosen difference of fourth powers. Frequency-dependent transmission, different excitations or additional transport channels can change this conclusion.

The finite-drop area-specific thermal resistance follows by factoring a difference of powers:Near thermal equilibrium, write $T_h=T+\delta T$ and $T_c=T$. Expansion gives $q=4\Sigma T^3\delta T+O((\delta T{)}^2)$, so the differential area-specific Kapitza resistance isAt a finite bias, varying only the hot side gives $(\partial q/\partial T_h{)}^{-1}=1/(4\Sigma T_h^3)$; this is generally different from $r_{K,\sec }$. Stating what is held fixed makes a differential resistance reproducible.

For an original normalized example, choose reference scales $T_0>0$ and $A_0>0$, and setFor contact area $A=2A_0$, the total heat transfer rate is $\dot{Q}=2.1472A_0\Sigma T_0^4$. The finite-drop resistance satisfies $\Sigma T_0^3{r}_{K,\sec }=0.2/1.0736=125/671\approx 0.18629$. Using the equilibrium tangent at $T_0$ instead would predict $j=4(0.2)=0.8$, illustrating why a substantial temperature rise requires the nonlinear law. Doubling the common equilibrium temperature reduces $R_K$ by eight in this chosen model.

A series check can be done without guessing any material properties. Stipulate a hot-side slab with dimensionless areal bulk resistance $b_s=\Sigma T_0^3{L}_s/k_s=0.10$, and neglect the cold-side bulk drop. At the same $j=1.0736$, the hot reservoir must have $T_H/T_0=1.2+b_{s}j=1.30736$. The imposed reservoir drop splits into $0.10736$ in the bulk and $0.20$ at the interface. Applying the fourth-power law directly to $T_H$ would assign the bulk drop to the interface and overestimate the heat flux.

Finally, transfer heat at rate $Aq$ between maintained interface-side temperatures $T_h>T_c>0$. The hot side loses entropy at rate $Aq/T_h$, and the cold side gains it at rate $Aq/T_c$. With no interface entropy accumulation, the entropy production per unit area isFor the normalized example it is ${\dot{S}}_{\mathrm{prod}}/(A\Sigma T_0^3)=1.0736(1-1/1.2)=0.178933\dots$. The same energy crosses both sides in a steady state, while total entropy increases, consistently with the Second law of thermodynamics. A finite temperature jump during heat transfer therefore expresses dissipation, not a failure of energy conservation.