Source term to VOF equation
Dear Foamers,
Having been working on solid particles - free surface flow interaction for a while, so far I got some good results in terms of coupled particle-fluid interaction. Nevertheless, an unnoticeable issue was found on the volume conservation of fluid phase. The modified solver is based on interFoam where an open-source DEM (discrete element method) code is called as a shared library to handle dynamics of solid particles. The interaction is realized via a momentum exchange term (Sp_f) to the R.H.S of N-S equation, and empirical drag models to calculate fluid forces acting on solid particles. The VOF and continuity equations used in the solver are given by: https://36.media.tumblr.com/efea4a8e...ho2_r3_540.png hence the Source term which stands for the displaced volume by solid particles reads as: https://40.media.tumblr.com/872018eb...ho3_r2_250.png here "alpha_f" is the volume fraction of fluid phase in a cell (= 1- V_particles/V_cell), and "U_f" the fluid velocity. The continuity equation is already implemented into the pressure equation instead of the original one. The alplaEqn.H is implemented as: Code:
{ https://36.media.tumblr.com/ad0e53d4...saho1_1280.png However. what i got from the coupled particle-fluid simulation is odd: the water level remains at 0.1 m which is not true after several hundred of particles poured into the water column. https://40.media.tumblr.com/712f4266...o4_r1_1280.png Where is the problem in the VOF implementation, especially the source term and the MULES solver. I read the MULES.c and MULESTemplates.C but still could not fully understand. Can we directly implement the alpha equation just like in UEqn.H without using MULES (e.g. bubbleFoam)? I would appreciate any helpful advice:) Cheers, David |
Update
Update
The direct implementation of the modified VOF method seems working, but still suffer some mass loss. Code:
{ thanks! |
Dear David Long
this model is also used for the lake modelin ?.For example in the lake water and solid particle are coming together and get settled at the bottom of the lake/river. any advice? thanks! |
Quote:
|
Hello,
I know this is an old thread but I'm dealing with exactly the same problem with water/air/particles solver (i.e. no volume displacement due to addition of particles). I'm able to maintain the volume fraction of the liquid phase and keep alpha1 between 0 and 1, but the water level doesn't rise accordingly. I've tried MULES and direct implementation. Have you managed back then to implement the alpha equation correctly? Could you share your results? Your help would be much appreciated. Thanks! |
Hi, did you try adding your alpha source term in alphaSuSp.H? That will make it get used by MULES correctly. There are a few equations where consistent source terms need to be added in for MULES, which is why they seem to have added that alphaSuSp file.
|
All times are GMT -4. The time now is 14:00. |