CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   twoPhaseEulerFoam: dispersed phase not conservative? (https://www.cfd-online.com/Forums/openfoam-solving/199582-twophaseeulerfoam-dispersed-phase-not-conservative.html)

RL-S March 9, 2018 10:11

twoPhaseEulerFoam: dispersed phase not conservative?
 
1 Attachment(s)
Dear Foamers,

when simulating a liquid-solid system with twoPhaseEulerFoam, it appears that the mass of the solid phase within the system increases - without actual influx of sediment.

To visualise this phenomenon, a simple rectangular, two-dimensional mesh was created (height: 0.5m, width: 0.6m): The lateral BCs are defined as cyclic, and the top and bottom as walls (no slip at the bottom, slip at the top, for velocities of both phases). Thus, no influx of sediment into the system should be possible.
The lower 0.2m of the domain are filled with a solid phase fraction of 0.6. This defines the initial volume of the solid phase as 0.6*0.2m*0.6m = 0.072 m/m.

All other BCs and initial values can be seen in the case, but are nothing fancy. The turbulence model for the solid phase is kinetic theory. Runs with OpenFOAM 3.0.

The expectation now would be that only a sedimentation process occurs, resulting in a stable phase fraction profile after some time. However, instead, the solid bed surface continually rises.

This can be seen in this video, and the attached phase fraction profile plot (alphaProf-time.png). The latter shows that the total solid phase volume in the domain increases, as the area below the curve clearly grows.

Does anyone know what might cause this?

Kind regards,
Lennart

P.S: I also tried setting all boundaries to walls, i.e. velocities to fixedValue (0 0 0). No significant difference in the described phenomenon.

RL-S March 16, 2018 12:56

Hmm. Feedback is not exactly overwhelming. Should I instead write this as a bug report or is that useless if I'm not using the current OF version?

axjjienn May 3, 2018 03:33

Post your case files so we can get a better idea of what is going on.

RL-S May 3, 2018 05:06

I wrote a link in the OP where you can download the case. Here's that link again:
https://www.dropbox.com/s/4ehmjvo3td...ation.zip?dl=0

BlnPhoenix May 8, 2018 04:49

Have you tested it with 3D? Is gravity direction correctly applied?

RL-S May 8, 2018 09:32

I have not tested it in 3D. The case is 2D.
I am quite sure that gravity is correctly applied, but even if it weren't, that would not explain the observable phenomenon: Apparently, judging by the concentration profiles, a phase change between liquid and solid phase occurs.

BlnPhoenix May 8, 2018 09:46

Quote:

Originally Posted by RL-S (Post 691612)
I have not tested it in 3D. The case is 2D.
I am quite sure that gravity is correctly applied, but even if it weren't, that would not explain the observable phenomenon: Apparently, judging by the concentration profiles, a phase change between liquid and solid phase occurs.

If it's 2D then you are correct. If the case was 3D a concentration profile would not be enough information.

What is your mesh resolution? Have you checked refinement? What about convergence, time step size? What about boundedness of your alpha field during calculation?

I suspect something in the setup, but i can not download the case and check everything.


All times are GMT -4. The time now is 20:47.