Planer shear flow using Lees-Edwards PBC

  • Simulation parameters for Input UDF file: filecolloid_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: 20000
    • Computation Time: 5 hour
  • Simulation parameters for Input UDF file: filecolloid_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: 20000
    • Computation Time: 5 hour
  • 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
ss.jpg

  • Related Papers