[kuria]

Hey everyone!

I am trying to model the transport of dissolved species due to bubbles motion in the liquid. Essentially it is a scalar transport in the fluid phase.

Since interfoam computes a single velocity field for both the liquid and gas phase (called as phi). So it is important to calculate the flux term for the liquid phase from phi.

I looked into some existing paper (that discusses the treatment of fluxes in alphaEqn.H) like: https://journals.aps.org/pre/abstrac...RevE.79.036306
http://infofich.unl.edu.ar/upload/3b...7523c8ea52.pdf

So my thought was I can calculate the phi associated with the liquid phase (phil) as
phil = phi + (1-r)phir
phir being the compressive velocity.

Did I understand this terminology correctly?


This is how the field was assembled before the MULES loop:
surfaceScalarField phir(phic*mixture.nHatf());

tmp<surfaceScalarField> talphaPhi1Un
(
fvc::flux
(
phiCN(),
cnCoeff*alpha1 + (1.0 - cnCoeff)*alpha1.oldTime(),
alphaScheme
)
+ fvc::flux
(
-fvc::flux(-phir, alpha2, alpharScheme),
alpha1,
alpharScheme
)
);