Dear all:
In the manual, there is an example LAMMPS input script for solid Ar-----utilizing compute heat/flux command to calculate average thermal conductivity. I make some modification to simulate iron’s thermal conductivity. But the result puzzled me. Besides some necessary modification, I only amend the N of “run N” command, nevertheless the results were vastly different.
So I write to ask for some help. Is it correct what I have modified or the example input script couldn’t be use like this? Or why in the same temperature, the results is so different?
The result of the first attachment is 124.5W/Km, the second is 2715W/Km, the third is 4545W/Km
the input script:
atom_style atomic
units metal
variable kB equal 1.3806504e-23 # [J/K] Boltzmann
variable kCal2J equal 4186.0/6.02214e23
variable T equal 300
variable V equal vol
variable dt equal 0.0001
variable p equal 200 # correlation length
variable s equal 10 # sample interval
variable d equal $p$s # dump interval

Dear all:
In the manual, there is an example LAMMPS input script for solid Ar-----utilizing compute heat/flux command to calculate average thermal conductivity. I make some modification to simulate iron’s thermal conductivity. But the result puzzled me. Besides some necessary modification, I only amend the N of “run N” command, nevertheless the results were vastly different.
So I write to ask for some help. Is it correct what I have modified or the example input script couldn’t be use like this? Or why in the same temperature, the results is so different?
The result of the first attachment is 124.5W/Km, the second is 2715W/Km, the third is 4545W/Km.
Maosheng
the input script:
atom_style atomic
units metal
variable kB equal 1.3806504e-23 # [J/K] Boltzmann
variable kCal2J equal 4186.0/6.02214e23
variable T equal 300
variable V equal vol
variable dt equal 0.0001
variable p equal 200 # correlation length
variable s equal 10 # sample interval
variable d equal $p$s # dump interval

Mail list readers (like me) generally don't have
time to read thru complicated scripts to figure
out what you are doing or what might be going wrong.
You need to post a more specific question. Someone
familiar with thermal conductiity may want to comment
on whether you can use this technique for Fe, since
the thermal conductivity of metals presumably depends
on electronic contributions.

Dear all:
In the manual, there is an example LAMMPS input script for solid
Ar-----utilizing compute heat/flux command to calculate average thermal
conductivity. I make some modification to simulate iron's thermal
conductivity. But the result puzzled me. Besides some necessary
modification, I only amend the N of "run N" command, nevertheless the
results were vastly different.
So I write to ask for some help. Is it correct what I have modified or
the example input script couldn't be use like this? Or why in the same
temperature, the results is so different?

if you run for different length and get different results, then
this is often an indication that your results are not converged
within the time you simulate and/or the size of the system you
are looking at. increasing the values of N in a systematic
way, should have them converge at some point. you can also
break down longer trajectories into chunks and analyse them
independently to get a feeling for how large fluctuations you
get on your results.