Is it possible to add random height while inserting particles?

Submitted by laudari on Wed, 07/14/2021 - 11:28

Hello everyone, I need some help.

Currently, I am doing simulation in LIGGGHTS. My particles insertion area and surface have height approximately 0.02 in si units.

I would like to add some random height so that my particles drop from different height in the same simulations. Purpose of this is to bring noise on heights.

So I want to drop particles randomly at the height of 0.02 +- 10 % or 0.02 +- 20 %while running simulation.

Is it possible to add this idea in liggghts script?

What will be the command?

Thank you.

laudari | Fri, 07/16/2021 - 09:04

Hello everyone, I need some help.

Currently, I am doing simulation in LIGGGHTS. My particles insertion area and surface have height approximately 0.02 in si units.

I would like to add some random height so that my particles drop from different height in the same simulations. Purpose of this is to bring noise on heights.

So I want to drop particles randomly at the height of 0.02 +- 10 % or 0.02 +- 20 %while running simulation.

Is it possible to add this idea in liggghts script?

What will be the command?

Thank you.

SHUBHAM AGARWAL | Sat, 07/17/2021 - 19:44

Hello
You can try the two ideas:
Idea -1 (Simple): Define region-1 for particle insertion of 0.02 unit height ( the normal one you do), then insert particles. Now define a region -2 from 0.02 - 10 to 0.02 + 10 % or 20%. Insert particles in it randomly. This should give you the distribution you want.

Idea-2 (dirty): Generate particle coordinates separately in MATLAB or python, save in .txt file. Then import the coordinates using 'read data'. I can help you with that, if you need.

Hope it helps
Shubham

laudari | Sun, 07/18/2021 - 05:48

Hi @SHUBHAM really appreciate your reply.

For idea 1:
Currently I am using insertion surface and region as mentioned below(for 0.02 m):
```
fix cad all mesh/surface/stress file meshes/inclined.stl type 1 wear finnie
fix inface all mesh/surface file meshes/insertion_face.stl type 1
group nve_group region domain
region bc cylinder z 0.0 0.0 0.015 0.05 0.12 units box

#particle insertion
fix ins nve_group insert/pack seed 32452867 distributiontemplate pdd1 insert_every 1 overlapcheck yes all_in yes particles_in_region 1 region bc
```
Do you mean by I create another insertion surface with 0.02+10% or - 1-% and define in same script to drop particles?

I will give a try for this one.

For idea 2:

I would be a great help if you I get your help for this one. Its a new idea for me.

Thanks a lot.

SHUBHAM AGARWAL | Sun, 07/18/2021 - 22:26

Hi @Laudari
Yes I meant exactly about the Idea-1 that you do the same for the region-2 as you have done for region-1.
In case you didn't get the expected results, we can try Idea-2. May be we can have a separate zoom call/chat. It would be much easier for me to explain

Shubham

laudari | Mon, 07/19/2021 - 05:57

Thank you @SHUBHAM , I did simulation with 3 regions with different insertion height. Let me analyze it. if I can't solve the problem them I will let you know. Really appreciate your time and help. Thank you so much.