sandbox/ecipriano/src/fixedflux.h

Fixed Flux Phase Change Model

This phase change model does not compute the vaporization rate from a mass or energy balance, instead, it uses a constant value mEvapVal defined by the user.

extern double mEvapVal;

Fields

We define the field mEvap with the vaporization rate per unit of surface in every interfacial cell, and we decleare the mEvapList used by evaporation.h.

scalar mEvap[];
scalar * mEvapList = {mEvap};

Phase Change Event

We loop over the interfacial cells and we assign the value of mEvapVal in the field mEvap.

event phasechange (i++)
{
  foreach() {
    mEvap[] = 0.;
    if (f[] > F_ERR && f[] < 1.-F_ERR)
      mEvap[] = mEvapVal;
  }
}