A colloidal dispersion under steady (DC) zigzag shear flow. †

• Simulation parameters
  • CFD Lattice: 128 x 128 x 128 mesh
  • Number of Particles: N = 500
  • Diameter of Particles: D = 8 mesh
  • Volume Fraction: Φ = 0.064
  • Temperature: kBT = 1
  • Shear Rate: 0.01
  • Number of Iteration Steps: 20,000
    
    

• Data plot

   > gnuplot
   ...
   gnuplot> set yrange [0:10]
   gnuplot> plot "out1" u 1:5 t "viscosity" w line

A colloidal dispersion under oscillatory (AC) zigzag shear flow. †

• Simulation parameters
  • CFD Lattice: 128 x 128 x 128 mesh
  • Number of Particles: N = 500
  • Diameter of Particles: D = 8 mesh
  • Volume Fraction: Φ = 0.064
  • Temperature: 0.01
  • Shear Rate: 0.01
  • Frequency: 0.05
  • Number of Iteration Steps: 5,000

• Data plot

   > gnuplot
   ...
   gnuplot> set yrange [-0.05:0.05]
   gnuplot> plot "out2" u 1:6 t "Shear stress" w line, "out2" u 1:2 t "Shear rate" w line

• Files needed for reproducing this simulation
  • Input UDF file: shear_dc.udf shear_ac.udf
  • Define UDF file: define.udf [v3.02]
  • Command (Remove "./" if you use Windows command prompt. Replace ">&" with "2>" if you use sh, bash or Windows command prompt.)

   > mkdir ./avs_dc
   > mkdir ./avs_dc/avs
   > ./kapsel -Ishear_dc.udf -Ooutput_dc.udf -Ddefine.udf -Rrestart_dc.udf >& out1
   > ./kapsel -Ishear_ac.udf -Ooutput_ac.udf -Ddefine.udf -Rrestart_ac.udf >& out2

• Visualization
  • GOURMET: particleshow_velocity.py
  • AVS/Express:　shear_avs.v

• Related Papers
  • Takuya Iwashita, Ryoichi Yamamoto, Direct numerical simulations for non-Newtonian rheology of concentrated particle dispersions, cond-mat arXiv:0905.0130