[gpumd] Re: Calculating thermal conductivity contribution of some monolayer on substrates.

  • From: Alex Gabourie <agabourie47@xxxxxxxxx>
  • To: Bruce Fan <brucenju@xxxxxxxxx>
  • Date: Thu, 11 Feb 2021 18:46:20 -0800

Hello All,

While Zheyong's suggestion is the correct/safe approach and should be used,
I still want to correct my original response in case someone decides to use
it anyways. My original suggestion to change N at line 85 of measure/
hnemd_kappa.cu was incorrect. That line should remain unchanged. In that
file, line 52 should be changed, replacing N with the number of atoms you
wish to consider for the thermal conductivity. For example, if you wish to
consider the first 1000 atoms of your structure, line 52 would look like: *if
(n < N) {*  --> *if (n < 1000) {*. The rest of my answer is still correct.
Again, it is better to use Zheyong's suggestion so mistakes like these
don't happen.

Regards,
  Alex

On Wed, Nov 25, 2020 at 6:08 AM Alex Gabourie <agabourie47@xxxxxxxxx> wrote:

Soonsung,

Zheyong is correct. His way is safer and will give you the correct result.
I believe that means you will also need to calculate the thermal
conductivity from the heat current, but that should be relatively
straightforward.

Good luck,
  Alex

On Wed, Nov 25, 2020 at 6:14 AM Bruce Fan <brucenju@xxxxxxxxx> wrote:

Using the HNEMD method, and output group heat currents using the compute
heyword, you can get what you want without changing source code.

Zheyong


"소순성" <soonsung2001@xxxxxxxxxx> 于 2020年11月25日周三 11:42写道:

Thank you Alex~!
I'm sure it will be very helpful advice for me!

And after I got your answer, I have additional question.

Following the analogous logic,
can I calculate the thermal conductivity contribution of each graphene
layer from bilayer graphene?


Other related posts: