---
title: Thermal Non-Equilibrium (TNE) Dynamics
url: https://www.emergentmind.com/topics/thermal-non-equilibrium-tne
type: topic
---

# Thermal Non-Equilibrium (TNE) Dynamics

Thermal non-equilibrium (TNE) denotes a class of nonequilibrium thermal states for which steady thermal closure fails under the imposed heating, cooling, transport, or collisional conditions. In solar-coronal physics, TNE refers to the absence of any steady fluid equilibrium in a coronal loop or flux tube despite quasi-steady heating, so the plasma undergoes repeated evaporation–condensation cycles rather than settling into a static or steady-flow state [1905.09767]. In kinetic and discrete Boltzmann formulations, TNE denotes the departure of the distribution function from local equilibrium, quantified through non-conserved kinetic moments such as non-organized momentum fluxes and non-organized energy fluxes [2205.13809]. A broader thermodynamic literature extends the term to discrete systems, nucleation, multiphase flows, and open coronal outflows, where temperature variables, internal variables, or mesoscale gradients cannot be reduced to a single equilibrium description [2010.14241][2302.08358][2411.10215].

## 1. Terminological scope and operative definitions

The modern literature uses TNE in several technically distinct but structurally related ways. In all of them, the common feature is that local or global thermal balance cannot be represented by a steady equilibrium state or by low-order constitutive closure.

| Domain | Operative definition of TNE | Representative manifestations |
|---|---|---|
| Solar coronal loops | No steady fluid equilibrium exists under quasi-steady, highly stratified heating | Evaporation–condensation cycles, coronal rain, prominences, long-period EUV pulsations |
| Discrete Boltzmann and compressible-flow modeling | Deviation of \(f\) from \(f^{eq}\) measured by non-conserved moments of \(f-f^{eq}\) | Viscous stress, heat flux, higher-order fluxes, shock-layer nonequilibrium |
| Discrete non-equilibrium thermodynamics | Enlarged state space with contact temperature \(\Theta\) and internal variables \(\xi\) | Accompanying reversible processes, embedding theorem, generalized cyclic efficiency |
| Kinetic nucleation | Different temperatures for gas, cluster translational motion, and cluster internal modes | Modified forward and backward nucleation rates, order-of-magnitude changes in cluster abundances |

In the solar literature, TNE is explicitly distinguished from thermal instability in the strict sense. Instability presupposes an equilibrium that can be perturbed, whereas TNE means that no equilibrium exists for the chosen heating and geometry; the local cooling episode is therefore better described as thermal runaway within a global limit cycle [1905.09767]. In discrete Boltzmann work, by contrast, TNE is the kinetic departure from local equilibrium itself, often contrasted with hydrodynamic non-equilibrium and organized through central moments of \(f-f^{eq}\) [1503.08450][2205.13809].

This suggests a broad unifying interpretation: across fields, TNE marks the failure of a reduced steady or near-equilibrium thermal description, but the mathematically relevant object differs by discipline—field-aligned fluid balance in coronal loops, kinetic moments in compressible flows, or non-equilibrium state variables in discrete thermodynamics.

## 2. Solar-coronal TNE as a limit-cycle process

In coronal loops, TNE occurs when heating is both highly stratified and quasi-constant or quasi-steady, usually concentrated toward the loop footpoints [1911.09710][1701.01309]. Under those conditions, footpoint heating drives chromospheric evaporation, the corona fills with denser plasma, radiative losses strengthen, the upper loop cools, and condensations form. The cycle can repeat as long as the heating remains sufficiently stable compared with the cooling time [1802.04010].

A standard heating prescription used in 1D loop studies is
\[
H(s) = H_0 + H_1\left(e^{-g(s)/\lambda_1}+e^{-g(L-s)/\lambda_2}\right),
\]
with
\[
g(s)=\max(s-\Delta,0), \qquad \Delta=5~\mathrm{Mm},
\]
where \(H_0\) is a weak background term, \(H_1\) sets footpoint heating strength, and \(\lambda_1,\lambda_2\) are the heating scale lengths in the two legs [1802.04010]. In this formulation, small \(\lambda\) concentrates heating low in the atmosphere, driving evaporation without adequately heating the apex. The loop then fills until radiative losses exceed available coronal heating and runaway cooling begins.

The global TNE cycle may terminate in either incomplete or complete condensations. In the parameter survey of Froment and collaborators, complete condensations were identified when temperature somewhere in the coronal part of the loop dropped below \(0.5\) MK, whereas incomplete condensations remained at coronal temperatures [1802.04010]. This distinction is observationally consequential: complete condensations are associated with coronal rain, while incomplete condensations can still generate long-period EUV pulsations without obvious chromospheric condensates [1701.01309].

Klimchuk and Luna derived two practical criteria that summarize the coronal-loop regime. As a rule of thumb, the ratio of apex to footpoint heating rates must be less than about \(0.1\), and asymmetries must be less than about a factor of \(3\), although the precise values are case dependent [1905.09767]. Their formulation emphasizes that strongly bottom-heavy heating is necessary, but large enough asymmetry in heating and/or cross-sectional area can replace TNE with a steady end-to-end flow whose enthalpy transport restores energy balance.

## 3. Observational manifestations in the corona

Long-period EUV pulsations and coronal rain are treated in the recent solar literature as two observational expressions of the same TNE mechanism [1802.04010]. Froment and collaborators argued that the observed pulsations are not wave signatures: for periods of \(3\)–\(16\) hr, standard MHD coronal loop modes are implausible, and the power spectra are better matched by periodic trains of pulses with random amplitudes superimposed on a power-law background [1608.03789]. Their model,
\[
f(t)=\sum_{m=0}^{M-1} a_m\,p(t-mT),
\]
predicts a power spectral density with both a continuum term from amplitude randomness and a harmonic comb from periodic repetition, consistent with observed Fourier and wavelet signatures [1608.03789].

A multi-scale observational demonstration was given with SDO/AIA and the Swedish 1-m Solar Telescope. Within the same loop bundle, a dominant coronal periodicity of about \(6\) hr was detected, with close peaks at \(42.2~\mu\mathrm{Hz}\) and \(46.1~\mu\mathrm{Hz}\), corresponding to \(6.6\) and \(6.0\) hr [1911.09710]. Time-lag analysis over the pulsating region showed positive lags for the coronal sequence \(335\rightarrow171\), \(211\rightarrow171\), \(193\rightarrow171\), and \(94\rightarrow171\), with median lags of \(88\), \(52\), \(31\), and \(18\) minutes, respectively, consistent with cooling from roughly \(2.5\) MK through \(2.0\), \(1.5\), and \(1.0\) MK [1911.09710]. DEM analysis found the DEM-weighted temperature varying between about \(2.0\) and \(2.5\) MK and the total emission measure varying from \(6.3\times10^{26}\) to \(15.7\times10^{26}\,\mathrm{cm^{-5}}\), with a relative amplitude of about \(60\%\) [1911.09710].

The same event linked the coronal cycle to chromospheric-scale rain. SST/CRISP and SST/CHROMIS resolved rain into narrow strands with typical widths of about \(0.2''\), or roughly \(140\) km, while overall rain blob widths ranged from \(100\) to \(600\) km with an average around \(330\) km [1911.09710]. Doppler velocities of \(20\)–\(60~\mathrm{km\,s^{-1}}\) and plane-of-sky velocities up to about \(150~\mathrm{km\,s^{-1}}\) were measured; average total velocities along pulsating-loop paths were around \(88\)–\(94~\mathrm{km\,s^{-1}}\) [1911.09710]. Combined H\(\alpha\) and Ca II K diagnostics yielded inferred average temperatures of roughly \(10.0\)–\(10.5\) kK for pulsating-loop paths, with inferred densities typically \(10^9\)–\(10^{12}\,\mathrm{cm^{-3}}\) and average values around \(7.8\)–\(8.5\times10^{10}\,\mathrm{cm^{-3}}\) [1911.09710]. Using blob densities around \(1.1\)–\(1.4\times10^{-13}\,\mathrm{g\,cm^{-3}}\), the observed blob number, and widths, the derived mass fluxes were roughly \(1.3\)–\(3.5\times10^9\,\mathrm{g\,s^{-1}}\), which the authors interpreted as significant mass circulation through the corona–chromosphere system [1911.09710].

Spectroscopic confirmation of the predicted coronal-temperature flows is more difficult. In a Hinode/EIS search across \(11\) raster datasets cross-matched with known AIA pulsation events, only \(2\) datasets showed periodic Doppler variations compatible with TNE downflows [1912.02538]. The clearest cases showed loop-leg velocity fluctuations of \(3\pm0.4~\mathrm{km\,s^{-1}}\) and \(2\pm0.6~\mathrm{km\,s^{-1}}\), much smaller than the \(\sim30~\mathrm{km\,s^{-1}}\) projected velocities expected from simulation [1912.02538]. The paper attributed the low detection rate to line-of-sight dilution, projection effects, low signal-to-noise ratio, limited cadence, and residual Hinode/EIS orbital drift. It explicitly concluded that the absence of a Doppler signature does not rule out TNE [1912.02538].

## 4. Geometry, asymmetry, and numerical control parameters

Large simulation surveys show that TNE occupies restricted regions of parameter space rather than arising generically in all loops. In a study of \(1{,}020\) 1D hydrodynamic simulations over three loop geometries, TNE occurred only for specific combinations of loop geometry, heating strength, and heating stratification [1802.04010]. In the semi-circular loop, TNE appeared mainly near symmetric heating; in one skewed observed loop, it shifted toward asymmetric heating with \(\lambda_1>\lambda_2\); and in another observed loop the preferred TNE region changed with heating strength [1802.04010]. The simulated periods ranged from about \(2.4\) to \(15.5\) hr, with longer loops generally showing longer periods, and the lag between temperature and density was typically about \(20\)–\(30\%\) of the cycle period [1802.04010].

A realistic loop study based on linear force-free field extrapolations showed that, for a loop of length \(367\) Mm and apex altitude about \(87\) Mm, TNE only appeared when the heating was asymmetric, with the eastern leg more strongly heated than the western one [1701.01309]. The successful simulation used
\[
H_0 = 10^{-7}\,\mathrm{W\,m^{-3}},\qquad
H_1 = 1.28\times10^{-4}\,\mathrm{W\,m^{-3}},
\]
\[
\lambda_1 = 50~\mathrm{Mm},\qquad
\lambda_2 = 20~\mathrm{Mm},
\]
and reproduced a \(9\)-hour periodic TNE regime with incomplete condensations [1701.01309]. The cross-correlation delay between temperature and density was about \(111\) minutes in the simulation, compared with \(119\) minutes from the observations [1701.01309]. The same paper emphasized that EUV cooling signatures and time lags are not uniquely diagnostic of impulsive heating, because a constant heating model can produce the same ordering across AIA channels [1701.01309].

The role of asymmetry in rain formation was addressed with a much larger survey of \(9{,}000\) simulations. Using the reconstructed geometry of the “rain bow” event and comparison cases from Froment and collaborators, the study found that both prominences and TNE cycles, with or without coronal rain, can form within the same magnetic structure [2110.09975]. Asymmetric loops were overall less likely to produce coronal rain regardless of the heating. In symmetric loops, coronal rain formed when the heating was also symmetric; in asymmetric loops, rain formed only when the heating compensated the asymmetry [2110.09975]. For the six stereoscopically based geometries denoted S10–S98, TNE occurred in \(78.5\%\) of the explored parameter space; among TNE cases, \(69.0\%\) produced coronal rain, \(20.0\%\) produced prominence-like structures, and \(11.0\%\) produced TNE without rain [2110.09975].

Numerical treatment of the transition region is also decisive. For a \(180\) Mm loop, inadequate transition-region resolution can suppress TNE entirely or distort its periodicity [1904.07287]. In HYDRAD, runs with refinement levels RL \(=1,2,4\) relaxed to static equilibrium rather than entering TNE, RL \(=5\)–\(8\) produced TNE but with substantially different cycle periods, and convergence required RL \(\ge 9\), corresponding to spatial resolution better than \(1.95\) km [1904.07287]. The converged reference period was about \(2.75\) hr. The paper further showed that TNE is suppressed unless the background heating is sufficiently small: for \(Q_{bg}\gtrsim9\times10^{-5}\,\mathrm{J\,m^{-3}\,s^{-1}}\), the loop relaxed to static equilibrium [1904.07287]. An unresolved-transition-region jump-condition method implemented in LareJ reproduced the converged behavior at far lower cost, with average coronal density and temperature errors of about \(3\%\) and \(4\%\), respectively [1904.07287].

## 5. Kinetic-moment formulations in compressible and multiphase flows

In discrete Boltzmann modeling, TNE is defined by non-conserved kinetic moments of the deviation from local equilibrium. A representative set of central-moment measures is
\[
\bm{\Delta}^{\sigma*}_{2},\qquad
\bm{\Delta}^{\sigma*}_{3,1},\qquad
\bm{\Delta}^{\sigma*}_{3},\qquad
\bm{\Delta}^{\sigma*}_{4,2},
\]
where \(\bm{\Delta}^{\sigma*}_{2}\) is the non-organized momentum flux (generalized viscous stress), \(\bm{\Delta}^{\sigma*}_{3,1}\) is the non-organized energy flux (generalized heat flux), \(\bm{\Delta}^{\sigma*}_{3}\) is the flux of \(\bm{\Delta}^{\sigma*}_{2}\), and \(\bm{\Delta}^{\sigma*}_{4,2}\) is the flux of \(\bm{\Delta}^{\sigma*}_{3,1}\) [2303.16692]. At first order, these quantities connect directly to macroscopic gradients through transport coefficients:
\[
\Delta^{\sigma*}_{2,\alpha\beta}
=
-\mu^\sigma \left(
\frac{\partial u_\alpha^\sigma}{\partial r_\beta}
+
\frac{\partial u_\beta^\sigma}{\partial r_\alpha}
-
\frac{2}{D}\frac{\partial u_\gamma^\sigma}{\partial r_\gamma}\delta_{\alpha\beta}
\right),
\]
\[
\Delta^{\sigma*}_{3,1,\alpha}
=
-\kappa^\sigma \frac{\partial T^\sigma}{\partial r_\alpha},
\]
so non-organized momentum flux depends on velocity gradients and viscosity, while non-organized energy flux depends on temperature gradients and thermal conductivity [2303.16692].

Several papers use scalar or vector summaries of TNE intensity. One recurring choice is
\[
D^*=\sqrt{
\mathbf{\Delta }_{2}^{\ast 2}
+
\mathbf{\Delta }_{3}^{\ast 2}
+
\mathbf{\Delta }_{3,1}^{\ast 2}
+
\mathbf{\Delta }_{4,2}^{\ast 2}
},
\]
or its density-weighted spatial average \(\bar D^*\), which acts as a compact measure of departure from equilibrium [1503.08450][2205.04027]. Another formulation introduces
\[
\mathbf{S}_{TNE}=(\tau,\Delta \mathbf{u},\Delta T,\bm{\Delta}_2^*,\bm{\Delta}_{3,1}^*),
\]
to emphasize that no single scalar fully characterizes nonequilibrium depth [2205.13809]. This literature explicitly distinguishes hydrodynamic non-equilibrium from thermodynamic non-equilibrium: non-central moments contain both organized flow and thermal deviations, whereas central moments isolate the thermal contribution [1503.08450][2502.01446].

As the degree of nonequilibrium deepens, higher-order kinetic moments become necessary. A family of ES-BGK-based discrete Boltzmann models was constructed to retain first- through sixth-order TNE effects, with D2V16, D2V25, D2V36, D2V49, D2V64, and D2V81 discrete velocity sets [2205.13809]. The study concluded that to obtain satisfying hydrodynamic quantities in a two-dimensional free jet, the model should include at least up to the third-order TNE effects [2205.13809]. In three-dimensional high-speed compressible flows, a second-order Burnett-level D3V55 model preserved the equilibrium moment set
\[
\bm{\Phi}=(\mathbf{M}_0,\mathbf{M}_1,\mathbf{M}_{20},\mathbf{M}_2,\mathbf{M}_{31},\mathbf{M}_3,\mathbf{M}_{42},\mathbf{M}_4,\mathbf{M}_{53}),
\]
allowing recovery of nonlinear constitutive relations for viscous stress and heat flux beyond the Navier–Stokes level [2502.01446]. A later super-Burnett D3V91 construction extended this to third-order TNE in 3D and showed that odd-order central moments, including heat flux and viscous-stress flux, are primarily governed by temperature gradients, whereas even-order central moments, including viscous stress and heat-flux-related flux, are dominated by velocity gradients; these leading dependences are substantially modified by density gradients, secondary gradients, higher derivatives, and cross-couplings [2606.05938].

In multiphase systems, TNE diagnostics are used not only to describe local kinetics but also to segment process stages. During phase separation, the TNE strength \(D\) attains its maximum at the end of the spinodal decomposition stage and then decreases during domain growth [1503.08450]. In shock–bubble interaction, viscosity has only limited influence on macroscopic parameters during the rapid shock-compression stage but significantly affects the TNE features, strengthening both non-organized momentum flux and non-organized energy flux and increasing entropy production from both channels [2303.16692].

## 6. Extensions to coalescence, open plasmas, nucleation, and discrete thermodynamics

The same mesoscopic language has been extended to droplet and bubble coalescence. In bubble coalescence, the averaged \(xx\)-component of the non-organized momentum flux, \(\bar\Delta_{2xx}^*\), identifies two characteristic instants: the time when the mean coalescing speed reaches its maximum, at which the ratio of minor to major axes is about \(1/2\), and the time when the ratio \(L_d=L_y/L_x\) reaches \(1\) for the first time [2205.04027]. The study reported that TNE intensity, coalescence acceleration, and the negative slope of boundary length show a high degree of correlation and attain their maxima simultaneously [2205.04027]. A related droplet-coalescence analysis found that the temporal evolutions of \(\bar D^*(t)\) and the mean entropy production rate \(\bar{\dot S}_{\rm pr}(t)\) provide concise and consistent stage criteria: maxima correspond to the fastest decrease of liquid–vapor interface length, and valleys correspond to the most elongated elliptical state of the merged droplet [2311.06546].

Under non-isothermal droplet coalescence, latent heat introduces a competition between temperature-gradient-driven and velocity-gradient-driven nonequilibrium effects. Before contact, \(\bar D_3^*\) and \(\bar D_{3,1}^*\) dominate; after contact, \(\bar D_2^*\) and \(\bar D_{4,2}^*\) become dominant [2502.17111]. The paper concluded that latent heat-induced temperature rise significantly refrains the TNE intensity in the thermal case, so the non-isothermal system has a lower peak TNE intensity, longer duration, and more complex spatiotemporal evolution than the isothermal case [2502.17111]. The same work reported a much longer coalescence initiation time in the non-isothermal case, with \(t_c=1.8\) for the isothermal case and \(t_c=5.6\) for the non-isothermal case [2502.17111].

Solar-plasma applications are no longer limited to closed loops. A one-dimensional HYDRAD study showed that TNE can also occur on open field lines in a transonic solar wind when sufficiently concentrated foot-point heating creates a radiation-dominated thermal sink in the low corona [2411.10215]. The imposed heating profile was
\[
\mathcal{H}(s)=A(s)^{-1}\left(\mathcal{H}_b e^{-(s-s_0)/\ell_b}+\mathcal{H}_g e^{-(s-s_0)/\ell_g}\right),
\]
with \(\ell_b=0.01R_\odot\), \(\ell_g=R_\odot\), and \(s_0=10\) Mm [2411.10215]. In the simulation sequence, steady solutions existed up to about \(\alpha\approx3.45\), while a tiny increase to \(\alpha=3.5\) triggered TNE cycles [2411.10215]. The marginal case had a cycle period of about \(24\) hr; as \(\alpha\) increased, the period dropped to about \(5\) hr by \(\alpha\approx4\) and approached about \(3\) hr for \(\alpha\gtrsim10\) [2411.10215]. This result directly generalizes the loop-centered view of TNE by showing that appreciable mass and energy fluxes do not preclude nonequilibrium cycling on open structures.

In kinetic nucleation, TNE means that the gas temperature, cluster kinetic temperature, and cluster internal temperature are not equal:
\[
T_{\mathrm{gas}},\qquad T_N^{\mathrm{kin}},\qquad T_N^{\mathrm{int}}.
\]
The forward rate coefficient is modified through the temperature-weighted reduced mass \(\mu_T\), and the backward rate is derived from the law of mass action under thermal non-equilibrium [2302.08358]. For homogeneous TiO\(_2\) nucleation at a gas temperature of \(1250\) K, differences in kinetic cluster temperatures as small as \(20\) K increased the formation of larger TiO\(_2\) clusters by over an order of magnitude; at \(1000\) K, an increase in cluster temperature of around \(20\) K reduced the formation of a larger TiO\(_2\) cluster by over an order of magnitude [2302.08358]. These results show that small thermal offsets can materially alter the stability of large clusters.

A more axiomatic formulation appears in the thermodynamics of discrete Schottky systems. There the non-equilibrium state space is
\[
Z=(a,n,U,\Theta,\xi),
\]
where \(\Theta\) is the contact temperature and \(\xi\) denotes internal variables [2010.14241]. Projection onto the equilibrium subspace generates a reversible accompanying process,
\[
PZ(t)=Z^{eq}(t),
\]
and the embedding theorem ensures compatibility between non-equilibrium entropy and equilibrium entropy along this projection [2010.14241]. The contact temperature is defined operationally by the inequality
\[
\left(\frac{1}{\Theta}-\frac{1}{T^\Box}\right)\dot Q \ge 0,
\]
and the paper shows that if entropy production is independent of internal energy, the contact temperature reduces to the thermostatic equilibrium temperature even though non-equilibrium is still under consideration [2010.14241]. This is a conceptually different use of TNE from the loop or DBM literatures, but it shares the same central theme: equilibrium temperature alone is not always an adequate state descriptor.

Across these literatures, TNE is therefore not a single model but a family of precise nonequilibrium constructs. In solar loops it constrains where and how heating is deposited; in discrete Boltzmann modeling it organizes the hierarchy of non-conserved kinetic moments; in multiphase and coalescence problems it tracks stage transitions and entropy production; in nucleation it controls forward and backward rates through small temperature offsets; and in discrete-system thermodynamics it motivates an enlarged state space with measurable non-equilibrium temperature. This suggests that the strongest common denominator of TNE is not a particular equation set, but the requirement to treat thermal evolution with variables or closures beyond steady equilibrium.

Source: https://www.emergentmind.com/topics/thermal-non-equilibrium-tne