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: 20000
Computation Time: 5 hour
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: 20000
Computation Time: 5 hour
Files needed for reproducing this simulation
> 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