Evaporation–Condensation Model

• Evaporation–condensation model is a class of formulations that treat phase change as a coupled evaporation and condensation process, with applications varying by domain.
• The model incorporates methods from kinetic jump conditions to diffuse-interface hydrodynamics to capture asymmetries and nonlocal effects in phase transitions.
• Its interdisciplinary applications span solar filament formation, capillary phase transitions, and finite-size statistical mechanics, offering actionable insights into material transport and stability.

Searching arXiv for recent/relevant papers on evaporation-condensation model usage across domains. arxiv_search.query: "all:evaporation condensation model OR all:evaporative-condensation model" “Evaporation-condensation model” is used in several technically distinct literatures, including solar filament formation, interfacial kinetic theory, capillary phase transitions, microcanonical first-order transitions, and thin-film or diffuse-interface dynamics (Yang et al., 2021, Malijevsky et al., 2014, Nogawa et al., 2011). A plausible unifying description is that these formulations all treat phase change or phase selection through a coupled evaporation–condensation process rather than as a single local constitutive step: material may be evaporated, transported, and later condensed, or evaporation and condensation may appear as asymmetric branches with different critical singularities, different kinetic boundary conditions, or different finite-size thresholds (Chen, 2022, Teshigawara et al., 2010).

1. Scope of the term

In the literature considered here, the term does not denote one universal equation. It denotes a class of models in which evaporation and condensation are treated as the central dynamical or thermodynamic mechanism, but the underlying state variables differ strongly by field. In solar physics, the model links magnetic reconnection, chromospheric evaporation, coronal cooling, and in situ condensation into filament material (Yang et al., 2021). In interfacial kinetic theory, it denotes a set of jump conditions relating mass flux, heat flux, and density or pressure discontinuities at a liquid–vapor boundary (Chen, 2022). In confined-fluid theory, it denotes condensation and evaporation transitions of a capillary meniscus controlled by groove geometry and long-ranged wall-fluid forces (Malijevsky et al., 2014). In finite-size statistical mechanics, it denotes the threshold between an evaporated supersaturated state and a condensed droplet state in the microcanonical ensemble (Nogawa et al., 2011).

Domain	Principal variable or object	Representative formulation
Solar filaments	Footpoint heating, evaporation, coronal condensation	Reconnection-driven evaporative-condensation cycle (Yang et al., 2021)
Liquid–vapor interfaces	$m$, $q$, $T_l$, $T_v$, $P_v$	Kinetic jump conditions from Boltzmann/BGK theory (Chen, 2022)
Capillary grooves	Meniscus position $\ell_C,\ell_E$	Asymmetric condensation and evaporation singularities (Malijevsky et al., 2014)
Microcanonical Potts model	Condensation ratio $\lambda$	Droplet threshold between supersaturation and phase separation (Nogawa et al., 2011)
Thin films and droplets	Interface profile $h$ or density field $n$	Diffuse-interface or variational evolution with phase change (Teshigawara et al., 2010)

This diversity matters because the same phrase can refer to rather different mechanisms. In some papers, evaporation and condensation are literal mass exchange between liquid and vapor; in others they are effective descriptions of meniscus unbinding, droplet formation under an extensive constraint, or immediate clipping of a transported scalar to a saturation profile. The commonality is therefore structural rather than purely phenomenological. This suggests that “evaporation-condensation model” functions as a cross-disciplinary label for coupled phase-loss/phase-gain dynamics rather than a single canonical formalism (Sukhatme et al., 2011).

2. Core physical mechanisms and mathematical structures

A first major class of formulations is microscopic or kinetic. In the SPC/E water study, evaporation is described as ballistic escape from a deep potential well with no additional barrier beyond the cohesive strength of the liquid, using the molecule–interface distance $a$ and the normal velocity component as reduced coordinates (Varilly et al., 2012). The transition state ensemble is well captured by

$q$0

and the outgoing normal-speed distribution is the barrierless flux form

$q$1

The reverse interpretation is a near-unity probability for a water molecule impinging upon a liquid droplet to condense, so the condensation or mass accommodation coefficient is near unity (Varilly et al., 2012).

A different microscopic picture is proposed in the atomistic barrier model, where evaporation is treated as escape of condensed-phase particles out of a potential well of depth $q$2, while condensation is treated as trapping of incoming vapor particles whose kinetic energy is lower than $q$3 (Semak, 2014). In that formulation, the evaporation flux is restricted by the threshold

$q$4

and the saturated vapor pressure is derived from the same interface-barrier picture: $q$5 This replaces the standard free-emission picture by an escape-and-trapping picture at the condensed-gas boundary (Semak, 2014).

A second major class is kinetic-theory jump modeling of nonequilibrium interfaces. Starting from the Boltzmann transport equation, one paper derives three interfacial boundary conditions—one for the mass flux, one for the heat flux, and a third for the density discontinuities—so that the interfacial temperature jump becomes quantitatively solvable rather than merely postulated (Chen, 2022). The mass-flux condition recovers the Schrage equation,

$q$6

but the model adds a heat-flux condition and a density or pressure jump condition, thereby coupling $q$7, $q$8, $q$9, $T_l$0, and $T_l$1 in a closed interface formulation (Chen, 2022).

A third class is continuum or diffuse-interface modeling, where evaporation and condensation are not imposed by an external kinetic law. In the dynamic van der Waals model for a one-component fluid, the density varies continuously across a finite-thickness interface, and phase change emerges self-consistently from coupled transport of mass, momentum, and entropy or heat (Teshigawara et al., 2010). In that framework the relevant equations are the compressible hydrodynamic equations with a Korteweg-type gradient stress, and local phase conversion is diagnosed a posteriori by the normal mass flux

$T_l$2

This shifts the modeling burden from an imposed interfacial flux law to a full nonisothermal diffuse-interface evolution (Teshigawara et al., 2010).

3. Solar evaporative-condensation model of filament formation

In solar physics, the evaporation-condensation model addresses the longstanding prominence or filament mass-supply problem by linking magnetic-structure formation and thermal mass loading in a single causal chain. The 2014 February 2 event in NOAA AR 11967 is interpreted as the first clear observation of a filament formed by magnetic reconnection and associated chromospheric evaporation and subsequent coronal condensation (Yang et al., 2021). In that event, sustained photospheric shear during flux emergence drove sequential tether-cutting reconnection, producing an M1.3 confined flare, an inverse-S sigmoid, and conjugate compact footpoint brightenings rooted in sunspots $T_l$3 and $T_l$4 (Yang et al., 2021).

The crucial spectroscopic diagnostic is explosive chromospheric evaporation at those conjugate footpoints. During 21:31:12–21:39:58 UT, Hinode/EIS observations showed redshifted He II 256.32 Å together with blueshifted Fe XV 284.16 Å and Fe XVI 262.98 Å at both footpoints. The measured He II redshifts were 14.1 and 14.0 km s$T_l$5, the Fe XV blueshifts were $T_l$6 and $T_l$7 km s$T_l$8, and the Fe XVI blueshifts were $T_l$9 and $T_v$0 km s$T_v$1 (Yang et al., 2021). Under the standard flare spectroscopic criterion, high-temperature blueshifts combined with chromospheric or transition-region redshifts diagnose explosive rather than gentle evaporation. The authors therefore infer that reconnection heated the chromospheric footpoints of the newly formed sigmoid and drove upward injection of hot plasma into it (Yang et al., 2021).

The subsequent condensation phase is spatially and temporally distinct. Initial condensation appeared at 22:50:38 UT, about 34 minutes after flare end, near the middle portion of the postflare loops in AIA 335 Å, then moved in opposite directions along the structure (Yang et al., 2021). Along the lower bright loops, condensations moved toward the positive end at about 43.3 km s$T_v$2 and toward the negative end at about 39.4 km s$T_v$3; along the inverse-S bright loops the corresponding mean speeds were about 36.6 and 37.1 km s$T_v$4 (Yang et al., 2021). Because the cool material appeared first in the middle of the coronal structure rather than being directly injected as already cool plasma from below, the event favors the evaporation-condensation picture over direct injection or pure levitation (Yang et al., 2021). The authors are explicit, however, that the plasma “might suffer from thermal instability or thermal non-equilibrium,” and they do not present a formal instability criterion (Yang et al., 2021).

This solar usage also includes reduced flare-hydrodynamic formulations. In the conduction-driven flare model, the post-heating transition region is treated as an isothermal Riemann problem that generates a downward chromospheric condensation shock, an upward evaporation shock, and a rarefaction wave (Longcope, 2014). The most important fitted relations are

$T_v$5

for the maximum evaporation velocity and

$T_v$6

for the condensation velocity, while the flare temperature obeys

$T_v$7

These relations were shown to fit the paper’s simulations and a variety of previously published conduction-driven flare simulations (Longcope, 2014).

The model has also been extended to time-dependent multidimensional filament dynamics. In a quadrupolar arcade, periodic localized heating before, during, and after filament formation drives cyclic evaporation-condensation and magnetic stretching, producing synthetic H$T_v$8 and 304 Å signatures of a “winking filament” (Zhou et al., 2023). In that case the imposed heating period is 20 min, but the oscillation or winking period is about $T_v$9 min, which the authors interpret with a forced-oscillator model rather than a direct one-to-one tracking of the heating period (Zhou et al., 2023).

4. Interfacial transport, confinement, and asymmetry between evaporation and condensation

Between two parallel plates, the evaporation-condensation model becomes a coupled liquid–vapor transport problem with kinetic interfacial jump conditions and continuum bulk transport. In one formulation, the bulk vapor satisfies

$P_v$0

while the liquid films are pure conduction layers,

$P_v$1

These bulk relations are closed with Hertz–Knudsen–Schrage mass-transfer laws and kinetic heat-flux laws at the two interfaces (Chen, 2023). The analytical criterion for vapor temperature inversion is rederived as

$P_v$2

and the paper proves that the temperature jump at each interface is in the same direction as the externally applied temperature difference: $P_v$3 at the evaporating interface and $P_v$4 at the condensing interface (Chen, 2023). The same model predicts that the vapor phase temperature can be much lower than the lowest wall temperatures and much higher than the highest wall temperature imposed, and that for sufficiently small latent heat evaporation can occur at the low-temperature side while condensation occurs at the high-temperature side, opposing the temperature gradient (Chen, 2023).

A more kinetic version of the parallel-plate problem solves the steady BGK equation in the vapor gap and heat conduction in the two liquid films, using accommodation and saturated half-range Maxwellian emission as interface boundary conditions (Chen et al., 14 Apr 2025). That model predicts an evaporative refrigeration effect: for $P_v$5 and $P_v$6, the vapor temperature can drop to about $P_v$7, below the cold wall temperature (Chen et al., 14 Apr 2025). The paper identifies two mechanisms: the dominant one is asymmetry between outgoing and incoming molecular distributions at the evaporating interface, and the second is additional cooling inside the Knudsen layer due to sudden expansion similar to the Joule-Thomson effect, although not isenthalpic (Chen et al., 14 Apr 2025).

In confined simple fluids, the asymmetry appears as a difference in the order and critical singularity of the transitions themselves. For a deep capillary groove with long-ranged wall-fluid forces and short-ranged fluid-fluid attraction, condensation and evaporation are controlled by meniscus unbinding from opposite ends of the groove (Malijevsky et al., 2014). Below the wetting temperature $P_v$8, condensation is first-order while evaporation is continuous; above $P_v$9, both are continuous, but with different exponents,

$\ell_C,\ell_E$0

The evaporation branch above $\ell_C,\ell_E$1 is covariant with planar complete wetting: $\ell_C,\ell_E$2 The asymmetry arises from long-ranged forces combined with the different substrate geometry at the groove bottom and at the opening (Malijevsky et al., 2014).

At moving contact lines, evaporation-condensation has also been proposed as a regularization mechanism for the hydrodynamic contact line paradox. In partial wetting with a volatile liquid wedge, Kelvin-effect-driven condensation near one part of the meniscus and compensating evaporation elsewhere can remove the classical no-slip singularity mathematically, while conserving total liquid mass in the diffusion-controlled problem (Janeček et al., 2016). The resulting matched asymptotics recovers a Cox–Voinov intermediate region, but numerical estimates for ethanol, water, and glycerol at ambient conditions give characteristic regularization lengths that are tiny. The paper therefore concludes that the mechanism is mathematically singularity-removing, yet likely not a practically sufficient continuum-scale resolution under ambient conditions (Janeček et al., 2016).

5. Finite-size statistical mechanics and stochastic moisture models

In finite-size statistical mechanics, the evaporation-condensation model describes a threshold between a homogeneous supersaturated state and a phase-separated state containing a macroscopic droplet. In the two-dimensional ferromagnetic Potts model with constrained internal energy, Wang–Landau sampling reveals a discrete finite-size transition between a supersaturation state and a phase-separation state in the microcanonical ensemble (Nogawa et al., 2011). The microcanonical inverse temperature

$\ell_C,\ell_E$3

shows the expected backbending and a negative-specific-heat coexistence region. The excess energy is decomposed as

$\ell_C,\ell_E$4

and the phenomenological theory predicts a condensation threshold at

$\ell_C,\ell_E$5

with a discontinuous jump of the condensation ratio to

$\ell_C,\ell_E$6

in the asymptotic finite-size scaling limit (Nogawa et al., 2011).

A later phenomenological theory modifies the Potts finite-size scaling by introducing an exponent $\ell_C,\ell_E$7 relating droplet interface cost to droplet volume (Ibáñez-Berganza, 2015). In that formulation,

$\ell_C,\ell_E$8

so the Ising-like compact-droplet value $\ell_C,\ell_E$9 gives $\lambda$0, whereas the 2D Potts proposal $\lambda$1 gives $\lambda$2 (Ibáñez-Berganza, 2015). The paper argues that $\lambda$3 is compatible with simulations up to $\lambda$4 and $\lambda$5 sites, and that the perimeter-versus-area behavior of Fortuin–Kasteleyn clusters is more compatible with $\lambda$6 than with $\lambda$7 (Ibáñez-Berganza, 2015). This is a direct controversy within the evaporation-condensation literature: the paper explicitly places its $\lambda$8 scaling in disagreement with previous $\lambda$9 claims (Ibáñez-Berganza, 2015).

A related stochastic usage appears in the advection-condensation model of atmospheric humidity, where condensation is represented as immediate removal of supersaturation,

$h$0

while moistening enters only through a reset source at a boundary (Sukhatme et al., 2011). Parcel advection is modeled as Brownian motion on an isentropic surface, and an exact Fokker–Planck solution gives broad, non-Gaussian, spatially inhomogeneous PDFs of specific humidity and relative humidity. Near the source the humidity PDF is bimodal; farther away it is unimodal; and the domain-averaged relative humidity PDF is bimodal with distinct moist and dry peaks (Sukhatme et al., 2011). Here the “evaporation-condensation” aspect is indirect: evaporation is not represented as an interior thermodynamic process, but the combination of localized moistening and rapid condensation creates the same structural competition between vapor retention and condensational clipping (Sukhatme et al., 2011).

6. Diffuse-interface, dewetting, and multiscale droplet formulations

In diffuse-interface hydrodynamics for a one-component fluid, evaporation and condensation are produced by the coupled evolution of density, velocity, and temperature fields rather than by an imposed interfacial jump law. The dynamic van der Waals model uses the entropy density

$h$1

and the corresponding Korteweg-type reversible stress

$h$2

to regularize a finite-thickness liquid–gas interface (Teshigawara et al., 2010). In complete wetting on a cooled substrate, a precursor film spreads and condensation localizes near the advancing film edge; the released latent heat produces a hot spot in the gas near that edge; and in late stages the film growth can be dominated by direct condensation onto the edge rather than by lateral flow from the droplet body (Teshigawara et al., 2010). Under weak heating, a steady-state circular thin film can form; under stronger heating, evaporation dominates and the liquid eventually disappears (Teshigawara et al., 2010).

In solid-state dewetting, evaporation-condensation can be formulated as a sharp-interface variational evolution. For an elastic thin film with moving contact lines, the evolution law is the $h$3-gradient flow of surface energies in the presence of epitaxial strain, and the intrinsic normal-velocity law is

$h$4

The graph PDE is fourth order rather than the sixth-order equation associated with surface diffusion, and the endpoint dynamics yields contact-line ODEs coupled to the bulk elastic equilibrium (Indulekha, 11 Aug 2025). The main result is short-time existence of weak solutions with moving contact lines in the dewetting regime $h$5 (Indulekha, 11 Aug 2025). This shows that, in materials science, “evaporation-condensation model” can denote attachment-detachment-controlled interface motion rather than mass-conserving surface diffusion (Indulekha, 11 Aug 2025).

A hybrid microscopic–mesoscopic droplet model combines lattice-gas kinetic Monte Carlo with a thin-film equation. In that framework, single-particle KMC moves represent diffusion, evaporation, and condensation as exchange of particles between droplets and the surrounding vapour, while collective thin-film moves represent advective hydrodynamic fluid motion (Areshi et al., 2019). The microscopic Hamiltonian

$h$6

is coupled to the continuum equation

$h$7

Evaporation and condensation are explicit only in the KMC sector; the thin-film sector is conserved and purely hydrodynamic (Areshi et al., 2019). The model distinguishes coarsening through vapour from coarsening over the surface and shows a temperature-controlled crossover: at higher temperature, vapour-mediated transfer dominates; at lower temperature, surface transport and droplet translation dominate (Areshi et al., 2019).

Across these formulations, the most persistent open issues are not merely technical. In solar applications, the decisive thermodynamic trigger is left as “thermal instability or thermal non-equilibrium” rather than a formally evaluated criterion (Yang et al., 2021). In kinetic interfacial theory, the magnitude and microscopic interpretation of accommodation coefficients remain debated, and classical simulations of water continue to support a near-unity condensation probability despite some experimental reports of lower values (Varilly et al., 2012). In contact-line theory, evaporation-condensation is mathematically elegant yet physically confined to nanometric or subnanometric scales under ambient conditions (Janeček et al., 2016). The literature therefore supports a broad conclusion: evaporation-condensation models are most successful when the relevant transport mechanism, interfacial closure, and asymptotic regime are specified explicitly, because the phrase itself spans fundamentally different mathematical objects.