Understanding fix insert/rate/region - particlerate & insert_every -.

mdasilva's picture
Submitted by mdasilva on Thu, 03/23/2017 - 20:48

Dear All, I've running binflow simulation available at LIGGGHTS Tutorial and couldn't understand the meaning of insert_every option of insert/rate/region fix.

I have to run a few simulations quite similar to this one implemented in the binflow example, and then went through the input file reviewing the documentation of each command used (by the way, this and other examples are very illustrative, congratulations).

This is the command I have in the input file:

fix ins all insert/rate/region seed 32452843 distributiontemplate pdd &
nparticles 150000 particlerate 150000 insert_every 1000 &
overlapcheck yes vel constant 0. 0. -1.0 region factory ntry_mc 10000

Two options are yet a bit confusing to me, particlerate and insert_every. What do they mean in this example?

According to the manual 150000 will be inserted by time unit (I am using SI units system), so it'll be 150000 particles/second. On the other hand, at every 1000 timesteps particles are inserted, and at each insertion 150000 / 1000 = 150 particles are inserted.

The number of particles inserted at each insertion is calculeted according to what I read in the insert/rate/region documentation.

Now, the last question (my apologies if this one is a bit silly), how should particlerate, insert_every and timestep size be set such that they keep a correct correlation? I mean, let timestep size = 1 x 10^(-5), particlerate = 150000 and insert_every = 1000. It'll be necessary 1000000 timesteps to insert all of the particles, but after such amount of timesteps the time elapsed is 10 seconds.

Is that allright?

Greetings,

ckloss's picture

ckloss | Sun, 04/23/2017 - 12:21

Hi,

>>Two options are yet a bit confusing to me, particlerate and insert_every. What do they mean in this example?

insert_every is how often particles are inserted, and the particle rate determines how many particles will be inserted with each event!

Best wishes
Christoph