LIGGGHTS® - User Forum

LIGGGHTS® related topics can be discussed here: discussion about models, installation, feature requests and general discussion

[SOLVED] Explosions in particles interactions

Submitted by xpt on Fri, 06/24/2016 - 06:37

Dears,
I'm simulating a vibrating box with particles inside (see box.png) that it moves in the y-axis.
In my simulation particles explodes at certain instant time (see exp.gif). I saw some forum topics about this, but it isn't help me.

I try changing some parameters and when I change Young's modulus(from 23 MPa to 5 MPa) it isn't exploit.
I attach input script and geometries.

How i can solve this without change the Young's modulus?

Best regards,

LIGGGHTS 3.3.1 - Public make errors with openmpi

Submitted by Stephanie.Poole on Fri, 06/24/2016 - 04:47

I am attempting to make LIGGGHTS with openmpi and I am getting the following error near the end:

atom.o: In function `LAMMPS_NS::Atom::new_avec(char const*, char*, int&)":
atom.cpp: (.text+0x12a8): undefined reference to `LAMMPS_NS::AtomVecTri::AtomVecTri(LAMMPS_NS::LAMMPS*)'
atom.cpp: (.text+0x13b7): undefined reference to `LAMMPS_NS::AtomVecTri::AtomVecTri(LAMMPS_NS::LAMMPS*)'
colect2: error ld returned 1 exit status

Could someone provide some guidance on how to resolve this issue? What step/steps am I missing?

Thank you!

Access to Atom information and change the drag force with velocity

Submitted by Kazuhisa Otsubo on Thu, 06/23/2016 - 16:13

Hi all

I'm a new user of LIGGGHTS 2.2. I'm trying to build a simulation to calculate slury flow with large particles.
Unfortunately I have one problem. I need to change the drag force on the particle with the velocity. (Ex. drag = 0.0 when v = 0.0, drag = 1.0 when v > 0.0 ) So I try to use the if-statement in the script but I can not access to the velocity information of the particle in the if statement. I don' why. Maybe I guess if-statement part is no good.

Please give me advise and helpful information !

# Slary Flow in Inclined Pipe

LIGGGHTS compute contact/atom command

Submitted by domenico on Thu, 06/23/2016 - 12:39

Hi everyone, I have a question.
The command compute contact/atom define a computation that calculates the number of contacts for each atom in a group.

The contact number is defined for finite-size spherical particles as the number of neighbor atoms which overlap the central particle, meaning that their distance of separation is less than or equal to the sum of the radii of the two particles.

For spherical particles, all with the same diameter, the maximum contact number should be 6, but I've got also 7,8 or 9 contacts. Is it possible?

filling of cylinder with raw particles

Submitted by Robert Blyth on Thu, 06/23/2016 - 11:52

Hello,

My final year part time master thesis is based on the shrinkage problem during the HIP (Hot Isostatic Pressing) process, which is based on my actual working experiences using the HIP process for an offshore oil and gas project.

Consequently I have a constructed the same dimensional cylinder (large volume cylinder) in autoCad Inventor replicating the same overall dimensions as previously used for this project. My thesis aim is to determine the various packing density values in the cylinder when raw particles are poured into the cylinder.

How does the processor distrib react to a change in change_box

Submitted by vkoppejan on Mon, 06/20/2016 - 22:06

Hi Everyone,

For my simulations I'm slowly growing the simulation domain using a change_box command with scaling in the z direction.

The docs tell me that "the mapping of processors to the simulation box is not changed from its initial 3d configuration".

Does this mean that the processors on the z boundaries of the domain get the extra simulation domain? If so, is there a clever way to perform a "fools" load balancing step?

Thanks in advance,

Victor

Geometry conversion from Solid Works to Binary/ASCII STL

Submitted by Mreza on Mon, 06/20/2016 - 20:06

Hello, I’m trying to incorporate LIGGGHTS for simulating a specific type of powder mixer. I’m using SolidWorks 2016 for creating the 3D geometries involved. In addition, I’m utilizing Paraview for post-processing purposes. Having read the user manual for LIGGGHTS, from what I understood the geometry files used in the execution script must be saved as ASCII STL format and that LIGGGHTS does not support Binary STL formats. When I attempted to open the ASCII files in Paraview the software stopped responding. Therefore, I can only view Binary files in Paraview.

Center of mass in multisphere (SOLVED)

Submitted by Westbrink on Mon, 06/20/2016 - 13:38

Hi all,

I am using the particletemplate/multisphere command to create clumps. The command allows the use of self-defined inertia tensor and mass.
I need to adjust also the center of mass. Currently LIGGGHTS uses the monte carlo method to calculate the center of mass.
Does anyone knows how to set the center of mass individually?

Thanks

Pages

Subscribe to RSS - LIGGGHTS® - User Forum