LIGGGHTS® - User Forum

Hi All,

I meant to ask if it is possible to include a command in the job input script to perform post-processing.

For instance, I have included in my script;

variable a universe 1 2 3 4 5 6 7 8 9

and defined results post directories as;

post_0.$a

My dump files have the format dump*.data so;

shell cd post_0.$a (I am sure gets me into the directories) and now I want to execute;

lpp dump*.data (Which I normally execute from the terminal - How do I include this in my job input script).

Thank you.

Hello all,

I am new in LIGGGHTS so would need some help:

As part of my degree final project (related to modeling of electro-filters) I am trying to simulate the interaction between charged particles with/without an external electric field.
I have been researching the manuals, the forum, etc, but I could not find a solution to my problem:

I do not know how to add coulomb interaction between particles. I found on this forum that -in previous versions- this could work:

Hi All,

Is there a way of conducting parametric study with one input script and saving dump files in different post directories. I mean say vary coefficient of rolling friction such as 0.1, 0.2,...,0.8, 0.9 and save the dump files in different folders say post_0.1, post_0.2,...,post_0.8, post_0.9.

Thank you.

Hi

I am trying a mulisphere break case.

In the LIGGGHTS Tutorials-public "multisphere_stone_restitution" example in.multisphere file, I add

"fix m6 all property/atom fz scalar yes no no 0.02"
before integrator, and change integrator to be

"fix integr all multisphere/break trigger_threshold 0.1 trigger_timeStep 0 trigger_fixName m6"

After type ./lmp_fedora -echo both < in.multisphere2, I get the ERROR of

I keep getting an error during bond creation:

creating bond btw atoms 1501 and 1503 (i has now 6 bonds) at step 1
creating bond btw atoms 1502 and 1501 (i has now 1 bonds) at step 1
creating bond btw atoms 1502 and 1503 (i has now 2 bonds) at step 1
creating bond btw atoms 1502 and 1504 (i has now 3 bonds) at step 1
ERROR: Internal fix bond/create error (../fix_bond_create_gran.cpp:465)

Hi All,

Attached is the insertion face stl and give embedded is a full script where I want to have coordination number of particles in a multiple particle size distribution system;

# Testing Script

atom_style sphere
atom_modify map array
boundary f f f
newton off

communicate multi vel yes

units si

region reg block -0.1 0.1 0. 0.40 -0.1 0.1 units box
create_box 3 reg

neighbor 0.03 multi
neigh_modify delay 0 every 1 check yes

# Material properties required for new pair styles

Hi ,
I am trying to set up a simulation with multiple spheres clumped together to form rod shape particle.
As I am still very new in LIGGGHTS, I am having dificulties in doing this simulation.

Can someone share example solved exactly as seen in LIGGGHTS Features > Non-spherical particles > Example - pouring a continuous stream of non-spherical particles, where can be seen particles made of 3 spheres stick together.

Thanks in advance.
Kashif

Hello,

Can I have an example of a simulation where both spheres and multispheres are used?
In particular I would like to specify more types of spheres (different diameters) and multispheres (different clusters).

I feel a bit lost when I have to specify the material properties of all the types.
Do I have to specify a set of types for spheres starting from 1 and another set of types starting again from 1?

Hi all.

I'm running a script where I have a mesh/surface shutter that is rotating counterclockwise from a starting angle. I want to unfix this mesh once it reaches a certain angle and I've been trying to think about the best way to do this in LIGGGHTS.

The best idea I had was to track the co-ordinate of a point on the mesh (i.e. a bottom corner, which is one of the co-ordinates I inputted to make the mesh) and use an if statement when this co-ordinate exceeds a certain value.

Hello,

I am having various problems while implementing SPH cases in LIGGGHTS. Firstly, I am not able to include .stl geometries into the fix wall/region/sph command as it does not accept a fix mesh/surface. I then converted the .stl into an ASCII VTK format and tried to introduce it as a region, which was successful. However, implementing this region in fix wall/region/sph is causing problems when it is an open geometry, such as a cylinder with one face open.