LAMMPS WWW Site - LAMMPS Documentation - LAMMPS Mailing List Archives
Re: [lammps-users] Voronoi/atom issue
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [lammps-users] Voronoi/atom issue

From: Eric Murphy <murphyericjames@...24...>
Date: Wed, 8 Nov 2017 22:33:41 -0600

Answering because nobody else has.

I cannot give an idea as to why it's >2x as big (more worryingly its not constant) - you should get the volume of the box if everything's done correctly.  The voronoi tesselation doesn't know anything about where the cluster ends and the void begins.

In the mean time, try the scientific method, test it on some simpler systems where things should behave better (non-clustered) and see what you get.

Eric Murphy
Graduate Research Assistant, Mechanical Engineering
Iowa State University

On Wed, Nov 8, 2017 at 12:01 PM, Rajdeep Behera <rajdeepbehera87@...24...> wrote:
Hi all,
I want to calculate the volume of a spherical nano-particle resides inside of a cubical simulation box (see the attached figure) utilizing Voronoi tessellation. My input commands are given below,

units      metal
dimension  3
boundary   p p p
atom_style  atomic

lattice      fcc 4.247
region       whole  block 0 50 0 50 0 50 units box
create_box   2 whole
region       Domain sphere 25 25 25 10 units box 
lattice      fcc 4.247
create_atoms 2 region Domain

region         Core sphere 25 25 25 6 units box 
group          Core region Core 

pair_style       eam/alloy
pair_coeff      * * AlGa.eam Al Ga
neighbor        2.0 bin
neigh_modify    delay 10 check yes

minimize        1e-15 1e-15 1000 100000

compute         v Core voronoi/atom only_group
compute         tv Core  reduce sum c_v[1]

dump            1 all custom 100 dat.test.*.cfg id type xs ys zs vx vy vz

thermo_style    custom step time temp vol c_tv

thermo          10
timestep        0.001

fix             relax all nvt temp 2.5 2.5 0.001

run 1000
write_data      Eql-NVT-1 

But, interestingly the volume of the nanoparticle is computed almost double of the simulation box. Here is my output,

Step Time Temp Volume c_tv 
      64            0            0                       125000    265676.44 
      70        0.006    50.253752             125000    265675.37 
      80        0.016 1.4260517e-010       125000    265682.03 
      90        0.026 3.2938985e-006       125000    265682.01 
     100        0.036  0.040481345          125000    265679.68 

Am I missing anything? 

Please make your comment/suggestion in this regard.



Check out the vibrant tech community on one of the world's most
engaging tech sites,!
lammps-users mailing list