[HectorLee]

Hello, everyone. Currently, I am simulating transport phenomena in multiphase flow and the solute only exists in one phase, which means that the solute cannot cross the interface between two phases (e.g., salt in the water-air system). With the help from this link:
InterFoam: Add an equation to only one phase
I added the transport equation to the solver as below:

Info<< "solve CEqn" << endl;
{
    surfaceScalarField Cf = fvc::interpolate(C);
    fvScalarMatrix CEqn
    (
        fvm::ddt(C)
      + fvm::div(phi, C)
      - fvm::laplacian(DT, C)
      == -fvc::laplacian(DT*Cf, alpha1)
    );
    
    CEqn.relax();
    fvOptions.constrain(CEqn);
    CEqn.solve();
    fvOptions.correct(C);
    
}

However, the solute still propagates through the interface somehow during the simulation as the attached figure shows.
To solve this problem, how should I modify the equation of the transport? Is there anyone who can enlighten me on this problem?
Thanks a lot!

[Z.Q. Niu]

Hi Li,
Have you solved this problem? Thanks!

[Blumenkind]

This problem has been solved.


It is implemented in OpenFOAM v1612
https://www.openfoam.com/releases/op...nd-physics.php
(see at the bottom of the page)



and in OpenFOAM-dev (for three weeks now)
https://twitter.com/CFDdirect/status...17062073167872
https://github.com/OpenFOAM/OpenFOAM...561205476f8b42