thermal conductivity with Muller-Plathe

Dear all,

I want to calculate the thermal conductivity with Muller-Plathe. I use “fix thermal conductivity ” . From this command, I know slab N/2+1 is the “hot” slab, for some reason, I have to change the “hot” slab on slab N. Like this:

Also the thermal conductivity is calculated as

What is different from the command is that I do not need the factor of 2, and I use the free boundary condition in the direction of heat flow. So I wonder if this is okay.

Best wishes,

CL Wang

1AB460DA@...4334....png

FF325672@...4334....png

Not sure what you mean by “is it ok”.
The doc page documents what the code
does. If you want it to do something different
you would have to change the code. If you
want to interpret the results differently (e.g. w/out
1/2 factor), that’s also up to you.

Steve

FF325672@...4334....png

1AB460DA@...4334....png