ABSTRACT
Numerical modeling of phase change requires accurate estimations of heat and mass transport at the interface. The present work develops a model in ANSYS-Fluent with user-defined functions to address phase change with a planar interface. An interface boundary method determines the heat fluxes with the exact location of the interface without interpolation functions. Five cases are analyzed based on the classic Stefan problem for validating the model. The numerical model is validated against closed form theoretical solutions with agreement within 0.55%. This work can be extended to include curvature effects and the interaction between the interface and heterogeneous surfaces.
Nomenclature
Aint | = | interface area in the cell (m2) |
cp | = | specific heat (J/kg K) |
dT | = | superheat/subcooled level (K) |
dy | = | cell length (m) |
F1 | = | phase 1 volume of fraction |
= | source term in the enthalpy equation (W/m3) | |
k | = | thermal conductivity (W/m-K) |
L | = | latent heat of evaporation (J/kg) |
l | = | domain length (m) |
p | = | pressure (N/m2) |
T | = | temperature (K) |
t | = | time (s) |
u | = | velocity (m/s) |
Vcell | = | cell volume (m3) |
Y | = | analytic interface displacement (m) |
y | = | y-coordinate (m) |
α | = | thermal diffusivity (m2/s) |
λ | = | growth rate constant |
μ | = | dynamic viscosity (Pa s) |
ν | = | ratio of thermal diffusivities between phases |
ρ | = | density (kg/m3) |
= | mass source term (kg/s m3) | |
ϕ | = | ratio of densities between phases |
Nomenclature
Aint | = | interface area in the cell (m2) |
cp | = | specific heat (J/kg K) |
dT | = | superheat/subcooled level (K) |
dy | = | cell length (m) |
F1 | = | phase 1 volume of fraction |
= | source term in the enthalpy equation (W/m3) | |
k | = | thermal conductivity (W/m-K) |
L | = | latent heat of evaporation (J/kg) |
l | = | domain length (m) |
p | = | pressure (N/m2) |
T | = | temperature (K) |
t | = | time (s) |
u | = | velocity (m/s) |
Vcell | = | cell volume (m3) |
Y | = | analytic interface displacement (m) |
y | = | y-coordinate (m) |
α | = | thermal diffusivity (m2/s) |
λ | = | growth rate constant |
μ | = | dynamic viscosity (Pa s) |
ν | = | ratio of thermal diffusivities between phases |
ρ | = | density (kg/m3) |
= | mass source term (kg/s m3) | |
ϕ | = | ratio of densities between phases |