wall heat flux boundary condition

Submitted by msandli on Wed, 09/23/2015 - 20:40

I'm trying to model some particles flowing over a heated surface. In my lab, I know how much power my heater consumes, and therefore know the wall heat flux.

Is it possible to apply this boundary condition to a wall in LIGGGHTS? I know how to give a wall a certain temperature, but as stated, I don't technically know that. If this is possible, does the wall have to be an imported mesh, or can it be a primitive? I see that I can compute the heat flux of a particle in contact with a wall, so I thought that making this a fixed value would be easy to do.

ckloss's picture

ckloss | Thu, 09/24/2015 - 10:02

Hi,

currently a heat flux BC is not implemented. That would be possible but a little tricky since you'd have to count the #of particles in contact and distribute the heat flux to them. Also, there would have to be some time-averaging as the # particle-wall contacts typically fluctuates considerably. You could implement some kind of P controller or PI controller to do this.

Best wishes
Christoph

msandli | Fri, 09/25/2015 - 00:05

Drat. If I knew the heat flux per area for a surface, it almost looks like it would just take a couple of lines of code to add the feature in - or at least the calculation. I don't know near enough about c++ code to add in all the support features (such as another keyword for fix wall/gran) myself.

If anyone ever reads this and decides to tackle this project, let me know!