LIGGGHTS® - User Forum

Hi All,

I want to know id of an atom during computation for that I am using property/atom command like this

compute ID impactor property/atom id
variable aID equal c_ID[1]

here impactor group only has one atom. And I access this with v_aID. But it throws an error :

ERROR: Compute used in variable between runs is not current. Use the update_on_run_end option for computes to avoid this.

After following this error message I change the compute command

compute ID impactor property/atom update_on_run_end yes id

Hi Christoph,

Can you please tell me what this error is?

WARNING: New thermo_style command, previous thermo_modify settings will be lost (../output.cpp:725)
ERROR: Invalid keyword in thermo_style custom command (../thermo.cpp:819)

Hi,
I'm trying to implement a liquid bridge model and wanted to use the available washino model in LIGGGHTS as starting point. So I just renamed the file:
cohesion_model_washino_capillary_viscous.h
cohesion_model_lqb.h

and within the file adjusted accordingly the name. Also for the identifier, gave it a value of 9:
#ifdef COHESION_MODEL
COHESION_MODEL(COHESION_LQB,lqb,9)
#else

When I try to compile the code I get error regarding redefinition of the model parameters. E.g.

Hey all,

I have got a 2D-array of different values for a force field depending on the region from another program (matlab) in a txt.file (other files are also easily possible).
I want to add these values via the addforce command to the particles in my simulation.
Depending on the position of each particle, it should give me an unique forcevector.

Hi all,

I have try to use the 'fix mesh/surface/stress/servo' to apply confining pressure on each sides of a cubic box. (100mm*100mm*100mm)

I want to apply 200kPa confining pressure on each sides with the following command. Take the left side as an example

fix zservo2 all mesh/surface/stress/servo file meshes/rwall4.stl type 1 com 0. 0. 0. ctrlPV force axis 1. 0. 0. target_val 2000 vel_max 100. kp 1e-2. ki 0. kd 0. #left

Hi guys,
I am currently using Red hat .
I tries to compile LIGGGHTS(3.5.0) by using make fedora but error occurs as follows :

/usr/bin/ld: cannot find -l-lsundials_cvode
collect2: ld returned 1 exit status
make[1]: *** [../lmp_fedora] Error 1
make[1]: Leaving directory `/share/apps/LIGGGHTS-PUBLIC-master/src/Obj_fedora'
make: *** [fedora] Error 2

so I tried make spirit instead of make fedora but it gets the same error :

hi to everyone

can any body tell me is there any command that must be used for different size particles ?

Hi,
I am using liggghts data for a macro-scale simulation, where there is aggregation happening leading to increase in the size of the particles. I want to restart the simulations from the previous state but with the new particle sizes. Is there a way to achieve that. Merely changing the sizes in the new infile leads to an insertion error.

Thank you for your help!

-Chai.

aoa

particle are not mixing like discrete . they move like chunk and do not mix properly ... i m uploading my input file plz help me because i have to submit my work this week. i m uploading a pics and code

#Poly dispersed flow

atom_style sphere
atom_modify map array
boundary f f f
newton off

communicate single vel yes

units si

region cyl cylinder z 0. 0. 0.075 0. 0.320 units box
create_box 2 cyl

neighbor 0.003 bin
neigh_modify delay 0

#Material properties required for new pair styles

Hello,

I want to determine average of velocity square < (Vx^2) > and later average it in a band at different heights in the domain.
I tried to do so in the following way using variable and fix ave/spatial command.

variable gvx equal vx*vx
variable gvy equal vy*vy
variable gvz equal vz*vz
fix gvel all ave/spatial 1 1000 1000 z lower 0.01 v_gvx v_gvy v_gvz file meanv.txt