Radius uniform mass

Submitted by billsouls on Thu, 12/04/2014 - 17:56

Hi all,

I am simulating a sand foundation with a linear distribution in mass. When I use "radius uniform mass 0.00015 0.0003", it reminded me "mass distribution not implemented for uniform". What does it mean? How could I solve this problem?

Many thanks,

Su

Forums:

JoshuaP | Thu, 12/04/2014 - 19:08

try "uniform number", that should work. But I would recommend to create an own discrete distribution. You can do it like this:
fix pts1 all particletemplate/sphere 1 atom_type 1 density constant 2720 radius constant ${radius1}
fix pts2 all particletemplate/sphere 1 atom_type 1 density constant 2720 radius constant ${radius2}
fix pts3 all particletemplate/sphere 1 atom_type 1 density constant 2720 radius constant ${radius3}
fix pts4 all particletemplate/sphere 1 atom_type 1 density constant 2720 radius constant ${radius4}
fix pts5 all particletemplate/sphere 1 atom_type 1 density constant 2720 radius constant ${radius5}
fix pts6 all particletemplate/sphere 1 atom_type 1 density constant 2720 radius constant ${radius6}
fix pts7 all particletemplate/sphere 1 atom_type 1 density constant 2720 radius constant ${radius7}
fix pdd1 all particledistribution/discrete 1 7 pts1 ${phi1} pts2 ${phi2} pts3 ${phi3} pts4 ${phi4} pts5 ${phi5} pts6 ${phi6} pts7 ${phi7}

where $phi is mass %

regards
Joshua

billsouls | Fri, 12/05/2014 - 03:33

Yes. This should work. We can also use "uniform number" and discrete distribution together to create a continuum distribution. Small modification on your codes:
fix pts1 all particletemplate/sphere 1 atom_type 1 density constant 2720 radius uniform number ${radius1} ${radius2}
fix pts2 all particletemplate/sphere 1 atom_type 1 density constant 2720 radius uniform number ${radius2} ${radius3}
fix pts3 all particletemplate/sphere 1 atom_type 1 density constant 2720 radius uniform number ${radius3} ${radius4}
fix pts4 all particletemplate/sphere 1 atom_type 1 density constant 2720 radius uniform number ${radius4} ${radius5}
fix pts5 all particletemplate/sphere 1 atom_type 1 density constant 2720 radius uniform number ${radius5} ${radius6}
fix pts6 all particletemplate/sphere 1 atom_type 1 density constant 2720 radius uniform number ${radius6} ${radius7}
fix pdd1 all particledistribution/discrete 1 7 pts1 ${phi1} pts2 ${phi2} pts3 ${phi3} pts4 ${phi4} pts5 ${phi5} pts6 ${phi6}

JoshuaP | Fri, 12/05/2014 - 11:09

you just have to consider that uniform number is distribution by number % and the other is by mass %. In your code now its mixed.