Particle growth rate edit

Submitted by ullassssss on Wed, 06/01/2016 - 09:40

A part of my syntax is given above.

#parameters for gradually growing particle diameter
variable alphastart equal 0.0001
variable alphatarget equal 0.000125
variable growts equal 50000
variable growevery equal 40
variable relaxts equal 20000

# Set the dumps
dump dmp all custom 10000 Ullas/shearcell/dump.1 id type type x y z ix iy iz vx vy vz fx fy fz &
omegax omegay omegaz radius
#dump dumpmesh all mesh/vtk 10000 Ullas/shearcell/dump*.vtk stress
#dump mov all movie 10000 Ullas/shearcell/movie.m4v type type size 1280 720

#calculate grow rate
variable Rgrowrate equal (${alphatarget}/${alphastart})^(${growevery}/(3.*${growts}))
print "The radius grow rate is ${Rgrowrate}"

#do the diameter grow
compute rad all property/atom radius

variable dgrown atom ${Rgrowrate}*2.*c_rad
fix grow all adapt ${growevery} atom diameter v_dgrown

#run
run ${growts}

#let the packing relax
unfix grow
run ${relaxts}

How can i edit the particle growth rate? Since the given syntax doesnt give me complete growth in the stipulated time step.