Hi
They updated the LAMMPS version on the cluster in my University from v30Jan14 to stable_17Nov2016 . However, when I run a test file (see below) I get the following error
LAMMPS (17 Nov 2016)
dimension 2
boundary p p p
units si
atom_style hybrid meso bond
ERROR: Invalid atom_style command (…/atom_vec.cpp:71)
The SPH package is installed and the test file works perfectly on my PC (I installed the latest Windows version). Thus, I suppose there is something wrong with the cluster installation. The problem seems to arise with the hybrid atom_style; test cases with only meso work fine.
I know I should probably attach the log files of the installation, but I don’t have those because the update is done centrally. Despite this, I’m wondering if anybody has encounter something similar and has any suggestion/tip?
Best Regards
Alessio
The test file is the following
dimension 2
boundary p p p
units si
atom_style hybrid meso bond
newton off
comm_modify vel yes
read_data pipe.initial
reset_timestep 0
variable R equal 0.01
… a bunch of variables here…
variable Nout equal 10/10/${Dt} #ogni 1. sec
group fluid type 1
group wall type 2
set group all meso/rho ${rho}
set group fluid meso/e ${e0}
set group wall meso/e ${ew}
set group all meso/cv ${Cp}
velocity all set 0.0 0.0 0.0
pair_style hybrid/overlay sph/rhosum 1 sph/taitwater mie/cut ${cut} sph/heatconduction
pair_coeff * * sph/taitwater {rho} {c} {alpha} {h}
pair_coeff * * sph/rhosum ${hm}
pair_coeff * * sph/heatconduction {DL} {h}
bond_style none
fix 10 wall meso/stationary
fix 20 fluid meso
fix 30 all enforce2d
fix 40 fluid addforce 0.0 v_nuova_forza 0.0
compute rho_peratom all meso/rho/atom
compute ie_peratom all meso/e/atom
compute de_wall wall property/atom de
compute ie_ave fluid reduce ave c_ie_peratom
compute de_tot wall reduce sum c_de_wall
variable T_fluid equal c_ie_ave/${Cp}
variable nuova_forza atom {m}*{g}${beta}(v_T_fluid-v_temperatura)
thermo ${Nout}
thermo_style custom step time v_T_fluid c_de_tot
thermo_modify norm no
neighbor ${skin} bin
dump dump_id all custom ${Nout} dump.lammpstrj id type x y z vx vy c_rho_peratom c_ie_peratom v_temperatura c_de_wall
timestep ${Dt}
run ${Nt}
Dr. Alessio Alexiadis
Lecturer in Chemical Engineering
School of Chemical Engineering
University of Birmingham
Edgbaston
Birmingham, B15 2TT
United Kingdom
Telephone +44 (0) 121 414 5305