Coupling with interDyMFoam solver.

Submitted by tshan on Mon, 02/28/2011 - 06:48

Hi, Chris!

I have coupled the interDyMFoam solver with LIGGGHTS by using you framework that has previously done in icoFoam solver. However, as you may know, I just modified the UEqn.H and pEqn.H by adding the fp term and the main cpp file interDyMFoam.c by adding the do particle stuff sentences.

However, I am not sure if it is correct because although the impact of fluid on particles is realistic, the impact of particles on fluids is not that good.

Do I need to modify other .H files such as the alphaEqn.H? If yes, how? Just deduct the alpha value by the void fraction?

If I want to consider the boundaries of the particles, is it mean that I should find another solver, such as the IB solver? But the IB solver seems to has not been released...

Thanks in advance!

Best,
Tong

cgoniva's picture

cgoniva | Mon, 02/28/2011 - 10:29

Hi Tong,

Congrats for creating a new solver!

Depending on the volumetric load in your case it might be necessary not to only couple the momentum balances, but also add another voidfraction field representing the granular phase. For single phase flow with particles this is not that difficult and well described in papers on CFD-DEM coupling. It should be possible to do the same approach for a multiphase (VOF) solver - but equations become less pretty.

Immersed Boundary approach might help - but also has its limitations. (less robust!, more computational effort...)

Best regards,
Chris

tshan | Mon, 02/28/2011 - 14:25

Hi, Chris!

Could you please give me some papers that well described the void fraction on CFD-DEM coupling?

Thanks again!

Best,
Tong

cgoniva's picture

cgoniva | Mon, 02/28/2011 - 16:13

Hi Tong,

the group of Prof. Kuipers has published many papers on CFD-DEM, e.g.:

Chemical Engineering Science 62 (2007) 195 – 207
Discrete element study of granulation in a spout-fluidized bed
J.M. Link,W. Godlieb, N.G. Deen∗, J.A.M. Kuipers

Cheers,
Chris

tshan | Tue, 03/01/2011 - 07:55

Hi Chris,

Thanks for your prompt reply!

I am not that clear about how to use the void fraction to generate volumetric load that you said.

In the paper you gave to me, it is an variable used to calculate drag force, and the drag force is used to do momentum transfer stuff. Then the void fraction seems to have not been used in any other fields. Did I misunderstand something?

What I want to see is that when a particle dropped down, it will exclude the water (volumetric load?). But it will be to complex to add the particle boundary into the fluid, so the void fraction is used. However, I do not know how to use the void fraction that represent the particles to generate the volumetric load.

Thanks!

Best,
Tong

cgoniva's picture

cgoniva | Tue, 03/01/2011 - 08:27

Dear Tong,

please see the fluids mass and momentum balance. There the volumefraction field is included to mass and momentum balance. This is how the volume of the particles is accounted for.

best regards,
Chris

maysmech | Mon, 06/17/2013 - 20:32

Dear Chris,
How is the implementation of volume fraction in mass and momentum balance in CFD part of CFDEM?
A solver like pisoFoam doesn't includes volumefraction in its equations i think.
Regards,

cgoniva's picture

cgoniva | Wed, 06/19/2013 - 16:26

Hi Maysmech,

pleas use a terminal and the grep command to find out where the term voidfraction is used - this is the voidfractionField (1-volumefraction).
Alternatively have a look at one of the publications on CFDEMcoupling: see http://cfdem.dcs-computing.com/?q=node/638

Kind Regards,
Chris