Tag Archives: drop evaporation model

Drop evaporation and phase change

the equations we solve are two: Mass conservation for each chemical species i and phase k :

\frac{\partial\rho_{k}\Upsilon_{ki}}{\partial t} + \nabla\cdot \left(\rho_{k}\Upsilon_{ki}\vec{u_{k}} - \rho_{k}D_{ki}\nabla\Upsilon_{ki}\right)=0

and energy conservation for each phase:

\frac{\partial\rho_{k} E_{k}}{\partial t} + \nabla\cdot\left( \rho_{k} E_{k}\vec{u_{k}}+p_{k}\vec{u_{k}}+\vec{q_{k}}\right)= \bar{\tau}:\nabla\vec{u_{k}}

 

where \Upsilon_{ki} is the mass fraction of phase k of species i and the others are the standard notation for density, pressure,velocity,time,energy,heat flux and diffusion coefficient. Moreover the flow is considered isobaric and the viscous dissipation is neglected.

The algorithm that solves the problem is included in the document.

phase-change-submodel

Leave a comment