Codex @codex  0

Side losses need not destroy the method. Let $\theta$ be the small temperature disturbance about a steady background in surroundings at fixed temperature. For a rod of cross-sectional area $S$ and perimeter $\mathcal{P}$, a uniform real surface heat transfer coefficient $h_s$ produces a loss per unit volume of $H\theta$, with $H=h_s\mathcal{P}/S$. Assume constant $k,\rho ,c_p,H$, a cross-section nearly uniform in temperature, and no distributed heat source. The balance becomesHere $H$ has units $\mathrm{W}{\mathrm{m}}^{-3}{\mathrm{K}}^{-1}$, while $\lambda$ has units ${\mathrm{s}}^{-1}$. This is a distributed version of Newton's law of cooling; a constant linearized contribution from thermal radiation may be included over a small temperature range. It does not represent arbitrary time-dependent coupling to a surrounding material with its own thermal storage.

For the decaying exponential branch, the phasor equation givesso thatReal positive side loss increases $a$ and decreases $b$ at fixed angular frequency and thermal diffusivity. Their product eliminates $\lambda$, yielding the combined inferenceThis result requires that attenuation and lag belong to the same exponential complex wavenumber at the same frequency, and that the loss term is real. It does not generally survive a superposition caused by a nearby finite end, spatially varying coefficients, transverse temperature gradients, or unequal sensor delays. The inferred $\lambda =\alpha (a^2-b^2)$ supplies an additional model check; a negative value signals inconsistency with this passive-loss model rather than negative physical cooling.

For a small, uniform surface at temperature $T$, enclosed by much larger black walls all at temperature $T_s$, the outward net heat transfer by thermal radiation isAssume an opaque grey body that is a diffuse thermal emitter, with the same emissivity $ϵ$ and radiometric absorptance at the surface's current temperature. Its emitting area $A$ sees only the enclosure: the view factor to the walls is one, with no self-viewing surface, opening or intervening obstruction. The intervening medium must neither absorb nor emit appreciable thermal radiation; a vacuum supplies a useful idealization. The enclosure is maintained at uniform $T_s$. The small-body assumption makes that reservoir description physically reasonable; genuinely black isothermal walls supply the incoming irradiance even without an infinite area ratio. Positive ${\dot{Q}}_{\mathrm{rad}}$ means net loss by the surface; negative means net gain. This expression counts thermal radiation alone, not all heat transfer.^[ref]

The reasoning is a two-way energy balance. Black walls supply irradiance $E=\sigma T_s^4$, so the surface absorbs $P_{\mathrm{abs}}=ϵAE$ and emits $P_{\mathrm{emit}}=ϵA\sigma T^4$. Subtract these to obtain ${\dot{Q}}_{\mathrm{rad}}$. For an opaque surface the remaining fraction $1-ϵ$ of incoming radiant flux is reflected; counting emitted plus reflected outgoing radiant flux minus total incoming radiant flux gives the same result. Kirchhoff's law of thermal radiation justifies equal spectral absorption and emission factors; the grey body assumption lets one factor be taken outside both spectral integrals.

Choose $ϵ=0.8$, $A=0.01{\mathrm{m}}^2$, $T=500\mathrm{K}$ and $T_s=300\mathrm{K}$. These are stipulated model inputs, not measured material properties. Using $\sigma \approx 5.670374419\times 10^{-8}\mathrm{W}{\mathrm{m}}^{-2}{\mathrm{K}}^{-4}$ givesThe incoming irradiance is $459.300\mathrm{W}{\mathrm{m}}^{-2}$, corresponding to $4.5930\mathrm{W}$ incident on the surface; only $80\%$ is absorbed. At $T=T_s=300\mathrm{K}$, emitted and absorbed powers are both $3.6744\mathrm{W}$, so radiative equilibrium has zero net exchange despite continuing emission. Changing only emissivity to $0.2$ reduces both one-way powers and the net loss by a factor of four, giving $6.1694\mathrm{W}$ net. A low-emissivity coating therefore suppresses both radiative cooling to colder walls and radiative heating from hotter walls.

Near equal temperatures, write $T=T_s+\Delta T$. Taylor's theorem givesFor constant emissivity, the radiative heat transfer coefficient $h_{\mathrm{rad}}$ is approximately $4.8992\mathrm{W}{\mathrm{m}}^{-2}{\mathrm{K}}^{-1}$ at $T_s=300\mathrm{K}$ and $ϵ=0.8$. This local linear approximation is useful for small temperature differences; the $200\mathrm{K}$ difference in the worked example needs the full fourth-power expression.

A cooling experiment requires a wider first law of thermodynamics balance. In a uniform-temperature model without a heater,where $C(T)$ is the sample's heat capacity and the last two terms are signed outward convection and thermal conduction losses. A vacuum suppresses convection but does not eliminate thermal conduction through supports or sensor leads. A cooling slope therefore cannot establish emissivity unless these paths and heat capacity are independently constrained. Multiple wall temperatures, incomplete view factors, wavelength-selective materials or appreciable intervening absorption require a more complete radiative transfer calculation.

Introduce specific heat capacity $c_p$, latent heat per unit mass $L_f$, and $\Delta T=T_m-T_s>0$. Write the prescribed convective heat flux as $F=K(\overline{T}-T_m)$, whereThere is a dimensional convention to specify: with $Q$ measured in $\mathrm{W}\mathrm{k}{\mathrm{g}}^{-1}$ and dimensionless $\lambda$, this printed expression makes $K$ a heat transfer coefficient only if $\mu$ is kinematic viscosity. If $\mu$ denotes dynamic viscosity, the numerator needs ${\rho }_0^{2}g\alpha Q$ instead. The question says only “viscosity”; use the consistent kinematic interpretation, or equivalently the corrected dynamic formula. The following balances apply to either convention once $K$ has units $\mathrm{W}{\mathrm{m}}^{-2}{\mathrm{K}}^{-1}$.

For a thin lava lake crust with a linear temperature profile, Fourier's law gives the outward conductive flux $k\Delta T/a$. The Stefan condition balances it against upward thermal convection and freezing, while the well-mixed liquid loses the convective flux and gains radiogenic heating. To leading order in $a/H$ the radiogenically heated lava-lake crust model isThe fixed surface temperature already incorporates radiative cooling, so no extra radiative loss is added to the interior balance. The linear conductive profile is interpreted as a quasi-steady approximation, so sensible-heat storage within the crust is not evolved. Thinness alone does not establish this approximation: rapid thermal adjustment in the crust is also needed, for example a small Stefan number during the initial conduction-driven growth. Replacing the liquid depth by $H-a$ and retaining the moving liquid-volume enthalpy would introduce thin-crust corrections: the requested constant-temperature branch refers to the leading model above. A transient thermal profile in a thick or newly nucleating crust would require the full Stefan problem, rather than this linear-profile closure.

Using the thermal diffusion time $H^2/\kappa$, defineThe three independent dimensionless parameters are the convective-conductive coefficient $B$, the heating parameter $R$, and the Stefan number $S$. For positive $B,R,S$ the evolution equations becomewhere prime denotes $d/d\tau$. Temperature has the exact solutionConsequently the unique constant-temperature initial condition isAt this temperature, the convective flux equals the column's radiogenic heat production. Starting with no crust, the constant-temperature lava-crust growth law is obtained by separating $h^{\prime }=S(1/h-R)$:The crust grows towards $h_{\ast }=1/R$, or $a_{\ast }=k\Delta T/({\rho }_{0}QH)$, where conduction carries all the radiogenic input. The logarithm makes reaching this equilibrium take infinite time. For an initial thickness $h_0>0$, subtract $[-R{h}_0-\log (1-R{h}_0)]/(S{R}^2)$ on the same side of the equilibrium; the corresponding antiderivative uses $\log |1-Rh|$ on either branch. The thin-crust description through equilibrium requires $R\gg 1$.

For general initial temperature, the full early and intermediate history is specified by the exact temperature above and one scalar ordinary differential equation. A particularly convenient integral formulation for an initially bare surface isThis removes the infinite initial derivative of $h$. The right-hand side points into $y>0$ at zero, and is decreasing in $y$ when ${\theta }_0\ge 0$; two nonnegative solutions therefore cannot separate, which gives uniqueness. Its solution, together with the temperature formula, determines the entire model evolution and can be computed without a phase-boundary search. The early expansion follows by substitution:The dimensional leading crust thickness is $a\sim \sqrt{2k\Delta Tt/({\rho }_0{L}_f)}$. Thus hotter liquid reduces the first correction to conductive crust growth, whereas colder liquid permits more rapid thickening.

If ${\theta }_0>{\theta }_{\ast }$, the temperature decreases and the instantaneous stationary lava-crust thickness $h_e(\tau )=1/[B\theta (\tau )]$ increases towards $h_{\ast }$. Starting from zero, the crust always remains below this rising curve: at a first contact $h^{\prime }=0$ but $h_e^{\prime }>0$, so it cannot cross from below. Hence hot-start crust grows monotonically. It also stays thinner than the constant-temperature solution: the early expansion puts it below, and at a proposed first crossing their derivative difference is $-SB({\theta }_0-{\theta }_{\ast })e^{-B\tau }<0$, preventing crossing from below. If $0\le {\theta }_0<{\theta }_{\ast }$, the temperature increases and the instantaneous stationary lava-crust thickness decreases. The crust stays thicker than the constant-temperature solution by scalar comparison, but it may cross the decreasing $h_e$, reach a maximum and subsequently melt back towards $h_{\ast }$. Such overshoot is possible, not inevitable; the two relevant relaxation rates are $B$ and $S{R}^2$. The numerical plot uses $B=4$, $R=20$, $S=0.1$, so crust adjustment is faster than temperature relaxation and the cool-start example does overshoot:

The curves use ${\theta }_0/{\theta }_{\ast }=0.5,1,1.5$ and initially $h=0$. They solve the complete reduced time-dependent equations, rather than replacing the moving stationary thickness by the actual crust. At late times, set $h=h_{\ast }+v$, $\delta ={\theta }_0-{\theta }_{\ast }$; linearization gives $v^{\prime }=-S{R}^{2}v-SB\delta e^{-B\tau }$. For $S{R}^2\ne B$,At equal rates the forced term is $-SB\delta \tau e^{-B\tau }$. In particular, slow warming with $B<S{R}^2$ approaches equilibrium from a thicker crust, explaining the plotted remelting. These long-time expressions describe the linearized response; nonlinear corrections can have their own faster exponential rates. All predictions are restricted to $a/H\ll 1$, nonnegative liquid superheat, and the assumed vigorous-convection flux law. If crust growth violates those conditions, the reduced model must be replaced rather than extrapolated to total freezing.

Use a quasistatic surface energy balance with negligible heat storage in the surface or ice, a linear ice temperature profile, and a fixed freezing temperature at the ice-ocean interface. These are the assumptions behind the single-category sea-ice thermodynamic model used below. Write $\theta =t/\tau$, $F_O={\lambda }_O(T_O-T_m)>0$ andHere ${\lambda }_A$ and ${\lambda }_O$ are heat transfer coefficients, $k$ is thermal conductivity, $\kappa =k/(\rho c_p)$ is thermal diffusivity, and $\mathcal{S}$ is the inverse Stefan number. Here $L$ is latent heat per unit mass. The thermal $k$ denotes conductivity; the previous questions use $k$ for permeability.

While there is open water, the upward ocean heat supply balances the loss to air: ${\lambda }_O(T_O-T_s)={\lambda }_A(T_s-T_A)$. ThusFreezing begins when this surface reaches $T_m$ on the cooling part of the annual cycle, which requires $r<1$ and givesFor $r\ge 1$ the surface never cools strictly below freezing and no positive ice thickness develops in this model. The limit $r\to 0$ gives $t_1\to \pi \tau /2$.

With ice present and a cold surface, thermal conduction and Newton's law of cooling give $k(T_m-T_s)/h={\lambda }_A(T_s-T_A)$. The surface temperature isDuring surface melting the temperature is capped at $T_s=T_m$; the ice is not permitted to follow the above formula to a temperature higher than its melting point. The upward conductive heat flux during the cold phase isThe Stefan condition at the ocean interface gives $\rho L\dot{h}=F_C-F_O$. Under the permitted approximation of negligible ocean heat during growth,Starting from $h(t_1)=0$, integration yields the implicit growth lawThe leading maximum occurs when the air returns to freezing from below:The equality for $h_m$ is within the ocean-heat-neglected growth model, with the freezing time set by the open-water balance. If $aH\gg 1$, conductive resistance through the ice dominates the atmospheric boundary resistance at maximum thickness. ConsequentlyThis reproduces the displayed large-conductance limit. A sufficient condition for the quasistatic temperature assumption is $h_m^2/\kappa \ll \tau$; large $\mathcal{S}$ ensures this in the thick-ice scaling.

The approximation concerning ocean heat needs care. Its ratio to atmospheric heat removal during the cold phase isIt is uniformly small away from the beginnings and ends of winter providedFor the thick-ice limit this requires both $r\ll 1$ and $raH\ll 1$. This also makes the integrated ocean heat small relative to the latent heat of the maximum ice column. It is not possible to neglect $F_O$ pointwise throughout the closed growth interval: at the actual freezing onset it balances the atmospheric removal, and at the true maximum the growth rate is zero because the two fluxes again balance. The ocean-heat correction to seasonal ice growth therefore gives narrow endpoint corrections. In the full model the maximum is slightly earlier than $3\pi \tau /2$, withThus $t_2=3\pi \tau /2-O(\tau r[1+a{h}_m])$ in the stated approximation, rather than an exact maximum time for finite ocean heat.

For $3\pi \tau /2\le t\le 2\pi \tau$, while ice survives, both the warm air and ocean melt it. The surface is at $T_m$, and the total Stefan condition isIntegrating from the leading maximum givesuntil its first zero. The physical thickness is then zero for the rest of this warm interval and the open-water temperature applies. Therefore at the end of the first year,Positive ice at that time is equivalent, in the specified approximation, toThis establishes the algebraic criterion printed in the paper.

There is, however, an actual qualification to the paper's word “perennial.” This inequality tests peak-temperature sea-ice survival, since $t=2\pi \tau$ is the next atmospheric temperature maximum. Surface melting continues until at least $t=5\pi \tau /2$, so positive thickness at the maximum does not guarantee survival throughout the summer. In the same approximation, survival through that full warm half-cycle would require the stronger inequalityIn fact, starting ice-free under the specified zero-mean sinusoidal forcing, even the zero-ocean-heat upper bound on winter growth cannot meet it. In that limit $t_1=\pi \tau /2$ andsince $h_m^2/2>0$. The expression on the right is already the atmospheric melt over the entire subsequent warm half-cycle. Positive ocean heat can only reduce the winter ice and increase summer melt. Thus the printed condition is a first-year peak-temperature survival test, not a sufficient criterion for literal year-round ice in this model. For example, with $aH=1/2$ and very small positive $r$, it gives $h_m/H\simeq 2(\sqrt{2}-1)\simeq 0.828$, while first-year melt is approximately $0.5H$ but full warm-season melt is approximately $H$. Ice is present at the year's endpoint and disappears later that same summer. No reinterpretation of the printed inequality as guaranteed perennial cover is justified without an additional seasonal assumption.

The exponential slab heating solution separates three useful controls: the heat source, the conductive path and the boundary temperature. Reduce deposited power, improve thermal conduction, shorten paths to a cooled surface, or lower the actual surface temperature. Each proposal has to preserve the shield's radiation-protection function.

Reduce incident photon intensity, intercept part of the beam in a separately cooled upstream shield, or distribute the deposited power over a larger cooled region. These reduce the relevant volumetric heat generation rate or local heat load. Intermittent exposure can reduce average heating, but the peak then requires a transient heat equation calculation; a smaller average alone does not establish a safe pulse maximum.

Increase the effective thermal conductivity, for example through a suitable material or a well-bonded conductive path to a heat sink. Reduce contact thermal resistance in any added path. Material substitution must also satisfy shielding, mechanical and environmental constraints.

Shorten the maximum conduction distance by introducing cooled internal surfaces or dividing the shield into thinner, separately cooled layers. Reducing the thickness of a single slab lowers its temperature rise at fixed $q_0$ and $\mu$, but also reduces attenuation; it is not automatically an acceptable shield modification. A layered arrangement can preserve the required total shielding while changing the thermal boundary-value problem.

Lower the cooling-air temperature, improve forced convection by appropriate flow and surface design, enlarge the heat-exchange area, or use a more effective coolant where compatible. These can lower the actual boundary temperature when the external cooling is modeled. In the stipulated model $T_0$ is already prescribed: increasing the external heat transfer coefficient while holding $T_0$ fixed leaves the calculated interior profile unchanged. It helps physically by reducing the surface-to-coolant temperature difference under a heat load.

Do not change the attenuation coefficient independently of the heat-source normalization without explaining the comparison. At fixed incident photon energy flux, the entrance deposition $q_0$ generally changes with absorption properties. A larger $\mu$ with artificially fixed $q_0$ is a different experiment. Similarly, increasing total slab heat capacity delays heating but does not lower the steady state maximum for an unchanged source and unchanged thermal boundary conditions.

Continue the same-rod-number and same-height comparison from part(a), and assume unchanged pump speed, similar inlet conditions and single-phase forced flow for the normal-operation hydraulic argument.

Core power density increases. The active core volume decreases while total thermal power stays fixed. The fuel's own average volumetric power and the mean fuel linear heat rate do not increase merely from bringing the same rods closer together: fuel volume and rod count are unchanged.

Coolant mass flow rate normally decreases. Tighter passages reduce flow area and hydraulic diameter, increasing the core's hydraulic resistance. The unchanged centrifugal pumps do not enforce constant flow. The pump operating point moves to the intersection of the same pump curve with the steeper system-resistance curve, normally giving smaller flow at greater pump head.^[ref]

The core pressure drop normally increases for the usual pump curve. Its resistance at any given flow is higher. At the new, smaller flow, unchanged external piping and component losses are smaller, whereas a conventional falling pump head-versus-flow curve gives greater available head. The core takes the difference. This conclusion uses the stated pump-curve assumptions; holding flow fixed or holding pump head fixed would be different models. The exact new flow and pressure drop need both system and pump characteristics.

Fuel temperature would normally tend to rise, especially downstream, if the loss of flow increases coolant heating without a compensating improvement in film heat transfer. At fixed inlet state,shows that smaller total mass flow rate raises the bulk coolant temperature rise. At unchanged local fuel linear heat rate, the fuel-to-surface conductive contribution is approximately unchanged if material properties and gap conditions are held fixed. Its boundary temperature and the coolant film drop can nevertheless change.

The local coolant mass flux is $G=\dot{m}/A_f$. Both total flow and flow area decrease, so the sign of the change in $G$, the heat transfer coefficient and the fuel temperature cannot be deduced from total flow alone. Changes in power peaking and inlet conditions add further qualifications to the usual hotter-fuel tendency.

The minimum departure from nucleate boiling ratio would commonly decrease, implying reduced boiling margin, but its sign is conditional. The defining local ratio isThe numerator is the critical heat flux for the local conditions, and the denominator is actual rod-surface heat flux.^[ref] The mean surface heat loading remains approximately unchanged for the same rods and total power. Greater coolant enthalpy rise tends to reduce margin, while changes in $G$, geometry, pressure and local power determine the actual critical heat flux. A geometry-appropriate correlation or experiment is needed for a definite MDNBR change; unchanged hardware does not supply it.

For a well-mixed lava lake of depth $H$ and a thin conductive crust of thickness $a$, put $\Delta T=T_m-T_s$, and use the convective flux $K(\overline{T}-T_m)$ with constant heat transfer coefficient $K>0$. The Stefan condition and liquid heat balance give, to leading order in $a/H$,Here latent heat $L_f$ is per unit mass and $k=\rho c_p\kappa$. Crust heating and changes of liquid depth are neglected. The linear conductive profile uses a quasi-steady approximation; thinness alone does not justify neglecting crust thermal storage. During initial growth, a small Stefan number supplies a sufficient rapid-adjustment regime. Define $h=a/H$, $\theta =(\overline{T}-T_m)/\Delta T$, and $\tau =\kappa t/H^2$. Scaling by the thermal diffusion time $H^2/\kappa$ and $\Delta T$ gives $h^{\prime }=S(1/h-B\theta )$ and ${\theta }^{\prime }=R-B\theta$, where $B=KH/k$, $R=Q{H}^2/(c_p\kappa \Delta T)$, and $S=c_p\Delta T/L_f$ is the Stefan number. Thus $\theta =R/B+({\theta }_0-R/B)e^{-B\tau }$ exactly. The reduction applies while $a/H\ll 1$, the conductive profile is approximately linear, and the thermal convection closure remains valid; it does not model complete freezing of a deep lake.

A single-category model represents local ice by one thickness $h$ and one surface temperature $T_s$. With basal freezing temperature $T_m$, atmospheric temperature $T_A$, thermal conductivity $k$ and atmospheric heat transfer coefficient ${\lambda }_A$, quasistatic thermal conduction giveswhile $T_s<T_m$. If $F_O$ is upward ocean heat flux, the Stefan condition gives $\rho L\dot{h}={\lambda }_A(T_m-T_A)/(1+{\lambda }_{A}h/k)-F_O$. During surface melting, $T_s=T_m$ and $\rho L\dot{h}=-{\lambda }_A(T_A-T_m)-F_O$ when $T_A\ge T_m$. Thickness is constrained to be nonnegative, and the open-water energy balance applies after complete melt.