FourField Solver

FourField Solver is cooler.

class Solvers.FourField.Wave(film)

Bases: Solvers.AbstractFilm

WAVE Summary of this class goes here Detailed explanation goes here

Constructor Summary

Wave(film): WAVE Creates a Wave(), wave Detailed explanation goes here

Property Summary

NZ = 0: % Solver properties

NTIME = 0: [-] Number of time steps

TIME = 0: [s] Time series

DT = 0: [s] Time step size

TIDX = 1: [-] Time step index

Z = 1.: [m] Elevation

HFLUX = 1.: [W/m^2] Film heat flux

W = 1.: Flow properties

U = 1.: [m/s] Velocity

H = 1E6: [J/kg] Enthalpy

FREQUENCY = 1: [Hz] Wave frequency

ITR: % Iteration properties

Method Summary

WL(wave, zIdx): WL() wave mass flow rate per unit perimeter

THICK(wave, zIdx): THICK() wave thickness

BETA(wave, zIdx): BETA() Wave film interfacial fraction

EPSILON(wave, zIdx): BETA Wave mass flow fraction

BETAP(wave, zIdx): BETAP() Wave film heat flux fraction

ETA(wave, zIdx): ETA() Wave film deposition fraction

MENT(wave, zIdx): MENT() Wave entrainment mass flux

MEVAP(wave, zIdx): MEVAP() Wave evaporation mass flux

MTURB(wave, zIdx): MTURB() Turbulent mass exchange

MBASE(wave, drop, zIdx): MWAVE Mass flux interaction with base

MTOT(wave, drop, zIdx): MTOT() Total mass flux

FWALL(wave, zIdx): FWALL() Film wall shear stress

FBASE(wave, zIdx): FBASE() Base interfacial force

FBASEMASS(wave, drop, zIdx): FBASEMASS() Wave mass exchange force

FVAPOR(wave, zIdx): FVAPOR() Film vapor shear stress

FDRAG(wave, zIdx): FDRAG() Film vapor shear stress due to drag

FSHEAR(wave, zIdx): FSHEAR() Film vapor shear stress

FTOT(wave, drop, zIdx): FTOT() Total

SHAPEFACTOR(wave, zIdx)

EQSTROUHAL(wave, zIdx): STROUHAL Correlation of equilibrium Strouhal number

EQFREQUENCY(wave, zIdx)

EQPERIOD(wave, zIdx): EQPERIOD() Inverse of EQFREQ

SPACING(wave, zIdx): SPACING() Wave spacing

PERIOD(wave, zIdx): PERIOD() Wave time period Inverse of wave frequency

N(wave, zIdx): N() Wave number density Inverse of wave spacing

WIDTH(wave, zIdx): WIDTH() Wave width Function of amplitude, Shape factor

AMP(wave, zIdx): AMP() Wave amplitude Function of WL, rho, Shape factor, frequency

DELTAFREQ(wave, zIdx): DELTAFREQ() Wave frequency source/sink

DRAGCOEF(wave, zIdx): DRAGCOEF() Wave drag coefficient

class Solvers.FourField.Drop

Bases: Solvers.ThreeField.Drop

DROP Summary of this class goes here Detailed explanation goes here

class Solvers.FourField.Base(film)

Bases: Solvers.AbstractFilm

BASE Summary of this class goes here Detailed explanation goes here

Constructor Summary

Base(film): BASE Creates a Base(), base Detailed explanation goes here

Property Summary

NZ = 0: % Solver properties

NTIME = 0: [-] Number of time steps

TIME = 0: [s] Time series

DT = 0: [s] Time step size

TIDX = 1: [-] Time step index

Z = 1.: [m] Elevation

HFLUX = 1.: [W/m^2] Film heat flux

W = 1.: Flow properties

U = 1.: [m/s] Velocity

H = 1E6: [J/kg] Enthalpy

ITR: % Iteration properties

Method Summary

EPSILON(base, zIdx): BETA Base film mass flow fraction

BETA(base, zIdx): BETA() Base film interfacial fraction

BETAP(base, zIdx): BETAP() Base film heat flux fraction

ETA(base, zIdx): ETA() Base film deposition fraction

MENT(base, zIdx): MENT() Base entrainment mass flux TODO: use coefficient later

MEVAP(base, zIdx): MEVAP() Base evaporation mass flux

MTURB(base, zIdx): MTURB() Turbulent mass exchange

MWAVE(base, drop, zIdx): MWAVE() Net Mass flux interaction with wave Defined as a source term, Mwave is the

MTOT(base, drop, zIdx): MTOT() Total mass flux of the base film

FWAVE(base, zIdx): FWAVE() Wave interfacial force

FWAVEMASS(base, drop, zIdx): FWAVEMASS() Base mass exchange force

FBASEVAPOR(base, zIdx): FBASEVAPOR() Film base vapor shear stress TODO: consider different way of doing this…

FDEP(base, drop, zIdx): FWAVE Droplet

FTOT(base, drop, zIdx): FTOT() Total

EQTHICK(base, zIdx): EQTHICK() Base equilibrium thickness

TBASE(base, zIdx): TBASE() Base film exposed time The time in which the base film is exposed to the vapor (eq.24)

TDRY(base, drop, zIdx): TDRY() Time duration in which the base is expected to dry out. Not included in the reference paper.

FDRY(base, drop, zIdx): FDRY() Dry fraction

WMIN(base, drop, zIdx): WMIN() Minimum base film mass flow rate (eq. 27) Base film mass flow rate at the end of the intermittent exposed time (TBASE).

WMINL(base, drop, zIdx): WMINL() Minimum base film mass flow rate per perimeter

THICKMIN(base, drop, zIdx): THICKMIN() Minimum base film thickness (eq. 28)

copyFlowProperties(srcObj, targetObj, opts): COPYFLOWPROPERTIES

class Solvers.FourField.Film(inputSet, fluid)

Bases: Solvers.AbstractFilm

FILM Summary of this class goes here Detailed explanation goes here

Constructor Summary

Film(inputSet, fluid): FILM Creates a Film(), film Detailed explanation goes here

Property Summary

NZ = 0: Solver properties

NTIME = 0: [-] Number of time steps

TIME = 0: [s] Time series

DT = 0: [s] Time step size

TIDX = 1: [-] Time step index

Z = 1.: [m] Elevation

HFLUX = 1.: [W/m^2] Film heat flux

MEVAP = -1.: [kg/s/m^2] Evaporation mass flux

ITR: Iteration properties

W: Flow properties

U

) double {mustBeNumeric} = 1. % [m/s] Velocity

Type:: (
Type:: ,

H

) double {mustBeNumeric} = 1E6 % [J/kg] Enthalpy

Type:: (
Type:: ,

Method Summary

distributeOAFW(film, W, zIdx): DISTRIBUTEOAFW Distributes film mass flow rate at onset of annular flow to base and wave, based on option: OAFFILMSPLIT

EQUIL(film, Wf, zIdx): EQUIL() calculate base/wave split ratio in equilibrium

initializeBaseAndWave(film, WIN, ITR): INITIALIZEBASEANDWAVE Initialize base and wave arrays given WTOT

copyFlowProperties(srcObj, targetObj, opts): COPYFLOWPROPERTIES