Convection coupling + nuclear heating of grains

Submitted by jtvanlew on Wed, 06/26/2013 - 21:46

I read what I could through the src file of cfd_coupling_convection and the convection coupling is being applied through an addendum of the DEM heat equation with a source term. In my earlier DEM simulations, I've been applying nuclear heating to pebbles and now want to consider an interstitial coolant gas. The nuclear heating is of course also with a source term added to the pebbles.

My question is just whether or not LIGGGHTS and/or the CFDEM coupling of it can handle a second source heating term. I'm concerned that if I apply nuclear heating it will overwrite the convection source or vice versa. I'm working on making a small sample case to experiment with adding a source term before or after the coupling fix to see if it either works at all, and if so if it works only when given in the right order.

But I thought I'd pose the question to the community overall to see if someone can answer concretely just from a better understanding of the coupling algorithm.


jtvanlew | Fri, 07/19/2013 - 23:33


Picking through the code and witnessing my own results, I can say confidently that the nuclear heating, as done with a
set group all property/atom heatSource 0.25
kind of command in LIGGGHTS will continue to work with CFD coupling. The coupling energy term is brought into the DEM with the heat flux, not as a body source term (as I had assumed it would in analogy to the momentum).

I'm happy that that's the case and I wasn't forced to do any internal coding as that would have been a huge mess :)