* 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