* Planer shear flow using Lees-Edwards PBC [#q90146c5]
#youtube(AAsNKH68RIw)
#youtube(P7T8rabkvU8)

- Simulation parameters for Input UDF file: &ref(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


#youtube(lEqJRGogDnw)

- Simulation parameters for Input UDF file: &ref(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: &ref(define.udf); [v3.02]
-- GOURMET: &ref(particleshow.py); &ref(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
#ref(ss.jpg,,75%);

----

- Related Papers