Planer shear flow using Lees-Edwards PBC †

• Simulation parameters for Input UDF file: colloid_LE_v31_t01.udf
  • CFD Lattice: 64 x 64 x 64 mesh
  • Number of Beads: N = 300
  • Diameter of Particles: D = 8 mesh
  • Volume Fraction: Φ = 0.307
  • Initial configuration: Uniform random
  • Temperature: kBT = 0.1
  • Shear Rate: s = 0.001
  • Number of Iteration Steps: 100,000
  • Computation Time: 5 hours on single core

• Simulation parameters for Input UDF file: colloid_LE_v55_t50.udf
  • CFD Lattice: 64 x 64 x 64 mesh
  • Number of Beads: N = 540
  • Diameter of Particles: D = 8 mesh
  • Volume Fraction: Φ = 0.552
  • Initial configuration: BCC with hollow space
  • Temperature: kBT = 5.0
  • Shear Rate: s = 0.001
  • Number of Iteration Steps: 100,000
  • Computation Time: 5 hour on single core

• Files needed for reproducing this simulation
  • Define UDF file: define.udf [v3.02]
  • GOURMET: particleshow.py particleshow_version2_cone.py
  • AVS/Express:

• Run

 > mkdir ./avs_ss
 > mkdir ./avs_ss/avs
 > ../../kapsel -Icolloid_LE_v31_t01.udf -Ooutput_dc.udf -Ddefine.udf -Rrestart.udf >& out1_v31_t01
 > ../../kapsel -Icolloid_LE_v55_t50.udf -Ooutput_dc.udf -Ddefine.udf -Rrestart.udf >& out1_v55_t50

• Data plot (temporal viscosity vs. strain)

 > gnuplot
 ...
 gnuplot> set logscale y
 gnuplot> plot "out1_v31_t01" u 3:7 t "viscosity at VF=0.31 and T=0.1" w line,"out1_v55_t50" u 3:7 t "viscosity at VF=0.55 and T=5.0" w line

• Related Papers