[Nicholas_Langbein]

Hello,

I need help to add gravity flux to alphaEqn.H for multiphase solver (interFoam, interMixingFoam)

I have volScalarField fsl (velcoity computed using strain rate)

i need flux in gravity direction with fsl*alpha1*alpha2 .. please help

volScalarField Us = 0.26 * shear_rate * d_s ; //unit will be m/s
volScalarField term = Us * alpha1 * (1.0 - alpha2);
volScalarField gravity_flux = fvc::grad(term).component(vector::Z);

solve
(
fvm::ddt(alpha1) // Time derivative term
- fvm::laplacian(DT, alpha1) // Diffusion term
+ fvm::div(phi, Cs) // Convection term
==
- gravity_flux // flux
);

But this dosent work meaning it dosent yield expected results... any idea