[hemph]

Hi,
I am simulating packing of particles in a column using the twoPhaseEulerFoam application.
In the middle of the column, I would like to have some kind of "filter", acting to stop one of the phases, while letting the other phase pass through unaffected. Has anybody tried something similar before?

Essentially, this would need a boundary condition in the middle of the column, but as far as I can tell, this approach is not possible in OF.

Some of the failing tactics used so far. I have tried explicitly setting the flux of phase a to zero at specified faces, but I get a 'seep through' of particles due to the common pressure term between the two phases. Besides, this feels like a somewhat 'messy' approach to the problem. I also tried setting the velocity of phase a to zero, using the setValues() function on the Ua-matrix, but since no equation system is solved for Ua, also this is unsuccesful.

Best Regards
Rasmus H

[mshehata85]

Hi,

Have you found a solution to your problem? I am curious to know how one could implement it in openfoam.
Best Regards

Mahmoud

[ali.m.1]

Hi Rasmus

I've done something similar with Lagrangian particles. I first created an .stl file for the surface (filter). Then created baffles, created regions, and assigned inlets/outlets to the patches. I then used patchToPatchInterpolate to interpolate the phases I wanted to keep, and set an 'outlet' BC for the phases I wanted to disappear.

There should be a tutorial on the chamers university page by the end of the month:
http://www.tfd.chalmers.se/~hani/kurser/OS_CFD_2016

It's the one by me (Alasdair Mackenzie)

Hope that helps