Per Particle Power Density

Submitted by mardussa on Wed, 07/18/2018 - 05:00

Hi everyone,

I have been attempting to look at the functionality of CFDEM in terms of being able to define the temperature (specifically the power density) of a particle that is generating its own head internally. I have a case where there is a large number of particles (1000s to 10,000s) that all have varying power densities, and these powers will change over time. I have been looking at past forum posts, and through the example problems but it appears that this is not something that is possible in CFDEM. I have a couple questions:

1) Is it possible to define a power density instead of a particle temperature?
2) Can a simulation contain particles of different power densities (ie temperatures), or is it possible to only have a uniform value over every particle in a system?
3) Can a CFDEM simulation be paused so that information can be written to a file to be passed to a separate solver (a neutronics solver), and then have that solver pass back updated power densities to the CFDEM simulation so that it can continue to a steady state simulation?

I posted this here and not the -dev forum, but if it is more appropriate to be placed there, please let me know.

Any help or suggestions on this would be greatly appreciated.

Kind regards,
Rob

mardussa | Wed, 08/15/2018 - 03:12

Thank you for your reply paul. Ok I guess not a large number of particles.
I have not had much experience looking through the CFDEM/LIGGGHTS src so I will have to look in to it and see if I can figure out how to write a fix that includes pebblewise power.
Thank you also for linking me the information on running LIGGGHTS as a library that can be coupled with an external code. I know that the external code I want to use has an example minimum coupling python script, I will look in to how I can get LIGGGHTS and CFDEM to fit in to that existing interface.
Regards,
Rob