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: 100,000
- Computation Time: 5 hours on single core
- 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: 100,000
- Computation Time: 5 hour on single core
- Files needed for reproducing this simulation
- Define UDF file:
define.udf [v3.02]
- GOURMET:
particleshow.py
particleshow_version2_cone.py
- AVS/Express:
- Define UDF file:
- 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
- Related Papers
Attach file:






