ERROR volume expectancy too small

Submitted by chow2003 on Fri, 06/07/2019 - 13:06

Hello,
As a part of my project i'm dealing with sand of less particle size. I got an error saying that volume expectancy too low and i tried to change units but I'm not succeeded. Here i attach my part of code.
Thank you.

### Initialization

# Preliminaries
units si
atom_style sphere
boundary f f f
newton off
communicate single vel yes

# Domain
region reg block -0.320 0.320 -1.960 1.460 -0.290 0.335 units box
create_box 3 reg

### Setup

# Material and interaction properties required
fix m1 all property/global youngsModulus peratomtype 2.5e7 2.5e7 2.5e7
fix m2 all property/global poissonsRatio peratomtype 0.25 0.25 0.25
fix m3 all property/global coefficientRestitution peratomtypepair 3 1 0.5 1 0.5 1 1 1 1 1
fix m4 all property/global coefficientFriction peratomtypepair 3 1 0.5 1 0.5 1 1 1 1 1

###Particle Insertion

region bc1 cylinder z 2 3 5 -5.0 1 units box
fix pts1 all particletemplate/sphere 12347 atom_type 1 density constant 1800.0 radius constant 4.7e-5
fix pts2 all particletemplate/sphere 12373 atom_type 1 density constant 1800.0 radius constant 3.45e-5
fix pts3 all particletemplate/sphere 12377 atom_type 1 density constant 1800.0 radius constant 2.55e-5
fix pts4 all particletemplate/sphere 12379 atom_type 1 density constant 1800.0 radius constant 1.9e-5
fix pts5 all particletemplate/sphere 12391 atom_type 1 density constant 1800.0 radius constant 1.4e-5
fix pts6 all particletemplate/sphere 12401 atom_type 1 density constant 1800.0 radius constant 1e-5
fix pts7 all particletemplate/sphere 12409 atom_type 1 density constant 1800.0 radius constant 7.5e-6
fix pts8 all particletemplate/sphere 12413 atom_type 1 density constant 1800.0 radius constant 5.5e-6
fix pts9 all particletemplate/sphere 12421 atom_type 1 density constant 1800.0 radius constant 4e-6
fix pts10 all particletemplate/sphere 12433 atom_type 1 density constant 1800.0 radius constant 3e-6
fix pts11 all particletemplate/sphere 12437 atom_type 1 density constant 1800.0 radius constant 2e-6
fix pts12 all particletemplate/sphere 12451 atom_type 1 density constant 1800.0 radius constant 1.5e-6
fix pts13 all particletemplate/sphere 12457 atom_type 1 density constant 1800.0 radius constant 1e-6
fix pts14 all particletemplate/sphere 12569 atom_type 1 density constant 1800.0 radius constant 5e-7
fix pdd1 all particledistribution/discrete/numberbased 15485959 14 pts1 0.007 pts2 0.045 pts3 0.1004 pts4 0.1322 pts5 0.1244 pts6 0.1007 pts7 0.0794 pts8 0.0655 pts9 0.0608 pts10 0.0601 pts11 0.0556 pts12 0.0449 pts13 0.0397 pts14 0.0843

achuth1992 | Tue, 06/11/2019 - 07:24

Hello,

It seems like you have the insertion region (region bc1 cylinder z 2 3 5 -5.0 1 units box) larger than the simulation domain itself (region reg block -0.320 0.320 -1.960 1.460 -0.290 0.335 units box).
Your insertion region should be within the simulation domain. Adjusting these parameters will fix the problem.

With regards,
Achuth