LAMMPS WWW Site - LAMMPS Documentation - LAMMPS Mailing List Archives
Re: [lammps-users] velocity plot from lammps output file
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [lammps-users] velocity plot from lammps output file


From: Sisir Das <sisirdas@...24...>
Date: Wed, 6 Sep 2017 11:22:50 +0530

It is solved . Need to change the fix ID.

Thanks
Sisir

On Wed, 6 Sep 2017 at 9:43 AM, Sisir Das <sisirdas@...24...> wrote:
HI Vanessa / Lammps User,

I run the code adding the velocity average lines . But I'm getting below errors :

LAMMPS (26 Jan 2017-ICMS)
OMP_NUM_THREADS environment is not set. Defaulting to 1 thread. (../comm.cpp:90)
  using 1 OpenMP thread(s) per MPI task
Lattice spacing in x,y,z = 1.28436 2.22457 1.28436
Created orthogonal box = (0 0 -0.321089) to (25.6871 22.2457 0.321089)
  1 by 1 by 1 MPI processor grid
Created 420 atoms
60 atoms in group lower
60 atoms in group upper
120 atoms in group boundary
300 atoms in group flow
Setting atom values ...
  60 settings made for type
Setting atom values ...
  60 settings made for type
ERROR: Replacing a fix, but new style != old style (../modify.cpp:798)
Last command: fix     4 all ave/chunk 20 250 5000 layers vx file profile.nemd.2d

I modified the code for averaging with below lines:

compute         layers all chunk/atom bin/1d y center 0.05 units reduced
fix             4 all ave/chunk 20 250 5000 layers vx file profile.nemd.2d

Sisir
On 5 September 2017 at 23:01, Vanessa Oklejas <voklejas@...24...> wrote:
See 'in.nemd.2d' file in the /examples/VISCOSITY in your lammps distribution. 

These 2 commands produce a velocity profile. See the fix ave/chunk documentation http://lammps.sandia.gov/doc/fix_ave_chunk.html to see how to specify how the velocity profile is averaged.

compute         layers all chunk/atom bin/1d y center 0.05 units reduced
fix             4 all ave/chunk 20 250 5000 layers vx file profile.nemd.2d



On Tue, Sep 5, 2017 at 10:01 AM Sisir Das <sisirdas@...24...> wrote:
Hi Team,

I was trying to get velocity profile from 2D Poiseulle flow.
I added

dump   5 all custom 500 velocities1.txt type x y vx vy

and created velocity of atoms as per attached dump.
How to get velocity profile from the dump ?
I understand create_box within the region , is it possible to get velocities from it with some modification ? or do we need to use some other tool with lammps output file ?

------------------------------------------------------------------------------
Check out the vibrant tech community on one of the world's most
engaging tech sites, Slashdot.org! http://sdm.link/slashdot_______________________________________________
lammps-users mailing list
lammps-users@...396...sourceforge.net
https://lists.sourceforge.net/lists/listinfo/lammps-users