LIGGGHTS® - Developer Forum

Topics related to developing with LIGGGHTS® can be discussed here: discussion about implementation details, C++, MPI and debugging tools

MODIFY FIX EFIELD SOURCE CODE

Submitted by pain08 on Fri, 03/05/2021 - 16:44

Hello everyone , i want to add the region command to the fix efield i did modify the .h and .cpp files in src folder but i still can not use an electrostatic field in a specified region with at an angle needed for my simulations i hope you can help me and thanks

this is the .h code

#ifdef FIX_CLASS

FixStyle(efield,FixEfield)

#else

#ifndef LMP_FIX_EFIELD_H
#define LMP_FIX_EFIELD_H

#include "fix.h"

namespace LAMMPS_NS {

voronoi computation called in pair style

Submitted by Jeff_DONG on Mon, 02/08/2021 - 07:36

Hello everyone,

I am calling the Voronoi compute in pair_gran_base.h to use for some pressure computation.

// setup voronoi info
Compute *voronoi_info;
char id_voro[] = "voroinfo";
int icompute = modify->find_compute(id_voro);
if (icompute < 0)
error->all(FLERR,"Could not find voro compute ID");
voronoi_info = modify->compute[icompute];

Then I just add "compute voroinfo all voronoi/atom" into the .in file

It is all good at this moment.

Triaxial Test

Submitted by hosseinhaeri on Sun, 02/07/2021 - 09:17

I would be really happy that any of you share a triaxial code with me as I am beginner in this field.
I have already wrote this code :

#####################################################################Triaxial Test######################################################################

#################################################################### Granular model ####################################################################

atom_style granular

atom_modify map array

boundary m m m

newton off

communicate single vel yes

how to release powder from a height

Submitted by nada on Thu, 01/07/2021 - 14:36

Hello,
I am trying to run a simulation where particles released freely from a height of 700 µm, due to gravitational and contact forces, these released powders freely deposited into the box and formed a powder bed. A roller is then applied to spread a 50 µm high layer of metal particles.
Can anyone share a code for that ? With which tutorials I can start or modify it or which command I should use ?
I am actually a total beginner. I would be glad if anyone can help through with the steps involved in the simulation.
thnx.

write_dump of gzipped file hangs simulation when running on several nodes

Submitted by iluvatar on Sun, 12/20/2020 - 04:24

Hello,
(Edit : Moved to liggghts users forum)
My simulations have been hanging when running them on several nodes. They work fine in a single node with several cpus. I was able to get a Minimal Working Example based on the packing example. Just by adding a single line
write_dump all custom out/dump_initial.gz id

Reset array from fix

Submitted by jagan1mohan on Fri, 12/11/2020 - 13:51

Hello friends, I created an array in hooke model to dump vdw force components. The values are accumulating over time. Let us say, I have two particles that are coming closer, colliding, and moving away in 10 steps. Here are a few details.

step fvdw fc
1 0 0
2 0 0
3 non-zero 0
4 non-zero non-zero // Collision
5 non-zero non-zero // Collision
6 non-zero non-zero // Collision
7 non-zero 0
8 non-zero 0
9 non-zero 0
10 0 0

Different initial condition for particle set

Submitted by mofazli on Fri, 12/11/2020 - 00:08

Hi,
I am working with liggghts to simulate the motion of a number of particles in a defined block. I have obtained some results for that, but I wanna see the initial condition effect on my result. To this end, I have changed the seed number in my insertion command

fix ins all insert/pack seed 32452843 distributiontemplate pdd1 vel constant 0. 0. 0. insert_every once overlapcheck yes all_in yes particles_in_region 417 region reg

Pages

Subscribe to RSS - LIGGGHTS® - Developer Forum