[kuria]

Hey everyone!

I am using interFoam solver for my simulation of species transport from the walls due to bubble rising.

I took the example case of bubble rising in a column and added a scalar transport equation :
solve
(
fvm::ddt(Y_K)
- v_K*fvm::laplacian(D_K, Y_K)
+ v_K*fvm::div(phi,Y_K)
);

As far as I understand interFoam solver, volume fraction of the phases is computed by a simple advection equation and the fluid properties like density and flow field in the whole domain is computed based on weighing of the fluid and gas fractions.

Due to this, when I run my equation for scalar transport, I end up solving it for the whole domain (both liquid and gas).

I also tried
solve
(
fvm::ddt(Y_K)
- v_K*fvm::laplacian(alpha1*D_K, Y_K)
+ v_K*fvm::div(phi,Y_K)
);
But this also gives diffusion into the bubble.

The only way I could fix it was by forcing the solution to be zero in the bubble like:
solve
(
fvm::ddt(Y_K)
- v_K*fvm::laplacian(D_K, Y_K)
+ v_K*fvm::div(phi,Y_K)
);
Y_K=Y_K*alpha1;


Has someone found a better way to address this problem?

[Z.Q. Niu]

Hi Kurian, have you found the solution? Please see this post:interFoam with transport: add a transport equation to only one phase . I hope it may help you!