# Variables used in multiple scripts # domain variables variable H0 equal 11.5 # [cm] variable EF equal 3.0 # [-] variable hb equal v_H0*v_EF # [cm] variable height equal 1.1*v_H0*v_EF # [cm], add 10% extra height to domain variable width equal 1.0 # [cm] variable depth equal 1.0 # [cm] variable eps equal 1.-$((1.-0.4)/v_EF) # average vol fraction, needed for particle insertion eqn variable rhop equal 2.8 # [g/cm3] # particle insertion variables variable dslice equal 0.2 # thicknes of a single slice variable slices equal floor(${hb}/${dslice}) # nr of insertion slices variable vslice equal $(v_width)*$(v_width)*v_dslice # volume of a single slice variable dred equal 0.5 # reduction factor for insertion size variable CG equal 2. # particle coarse graining factor # Growth variables variable dred equal 0.1 # reduction factor for insertion size variable growts equal 25000 variable growevery equal 1 variable relaxts1 equal 10000 variable relaxts2 equal 100000 # granular variables (based on eqns from Landry 2003 ) variable dpav equal 0.025 variable mpav equal $(v_rhop*PI/6.*v_dpav*v_dpav*v_dpav) variable g equal 981. variable kn equal $(2.e5*v_mpav*v_g/v_dpav) variable kt equal $(2./7.*v_kn) variable gamman equal 50*$(sqrt(v_g/v_dpav)) variable gammat equal ${gamman} variable mu equal 0.15 print "included variables"