Discontinuous shear thickening on dense non-brownian suspensions via lattice Boltzmann method
|
|
- Todd Phelps
- 5 years ago
- Views:
Transcription
1 Discontinuous shear thickening on dense non-brownian suspensions via lattice Boltzmann method Pradipto and Hisao Hayakawa Yukawa Insitute for Theoretical Physics Kyoto University Rheology of disordered particle part 1, YITP, June 22nd /38
2 Outline Introduction Outline of our model Hydrodynamics: lattice Boltzmann method Particle contacts Electrostatic repulsive forces Sheared suspensions simulations Results Discussions Future prospect 2/38
3 Suspensions 3/38
4 Suspensions Solid particles suspended in solvent fluid 3/38
5 Suspensions Solid particles suspended in solvent fluid Fluids are described by Stokes equation (Re 0): u = 0 p = η 2 u 3/38
6 Suspensions Solid particles suspended in solvent fluid Fluids are described by Stokes equation (Re 0): u = 0 p = η 2 u Suspensions rheology have unique behavior! 3/38
7 Discontinuous shear thickening 4/38
8 Discontinuous shear thickening (experimental observations) H. A. Barnes, J. Rheol. 33, 329 (1989). 5/38
9 Discontinuous shear thickening (experimental observations) R. G. Egres and N. J. Wagner J. Rheol. 49 3, (2015) H. A. Barnes, J. Rheol. 33, 329 (1989). 5/38
10 Discontinuous shear thickening (Stokesian Dynamics simulations) Seto. et al, Phys. Rev. Lett 111, (2013) Mari. et al, J. Rheol.58(6): (2014) 6/38
11 Discontinuous shear thickening (Stokesian Dynamics simulations) Seto. et al, Phys. Rev. Lett 111, (2013) Mari. et al, J. Rheol.58(6): (2014) 6/38
12 Discontinuous shear thickening (Stokesian Dynamics simulations) Seto. et al, Phys. Rev. Lett 111, (2013) Mari. et al, J. Rheol.58(6): (2014) 6/38
13 Discontinuous shear thickening 7/38
14 Discontinuous shear thickening No Theory 7/38
15 Discontinuous shear thickening No Theory Possible explanations: Lubrication (Brady s group) 7/38
16 Discontinuous shear thickening No Theory Possible explanations: Lubrication (Brady s group). Friction (R. Seto, et al, 2013) 7/38
17 Discontinuous shear thickening No Theory Possible explanations: Lubrication (Brady s group). Friction (R. Seto, et al, 2013) Granular systems (Otsuki and Hayakawa, 2011) 7/38
18 Discontinuous shear thickening No Theory Possible explanations: Lubrication (Brady s group). Friction (R. Seto, et al, 2013) Granular systems (Otsuki and Hayakawa, 2011) Brownian suspensions (R. Mari, et al, 2015) 7/38
19 Discontinuous shear thickening No Theory Possible explanations: Lubrication (Brady s group). Friction (R. Seto, et al, 2013) Granular systems (Otsuki and Hayakawa, 2011) Brownian suspensions (R. Mari, et al, 2015) Boundary effect (Brown and Jaeger, 2012) and (Allen, et al, 2017) (Brown s group) the contributions beneath the DST is still unclear Motivations: Recover DST using LBM 7/38
20 Discontinuous shear thickening No Theory Possible explanations: Lubrication (Brady s group). Friction (R. Seto, et al, 2013) Granular systems (Otsuki and Hayakawa, 2011) Brownian suspensions (R. Mari, et al, 2015) Boundary effect (Brown and Jaeger, 2012) and (Allen, et al, 2017) (Brown s group) the contributions beneath the DST is still unclear Motivations: Recover DST using LBM Decompose all contributions to the shear stress 7/38
21 Discontinuous shear thickening No Theory Possible explanations: Lubrication (Brady s group). Friction (R. Seto, et al, 2013) Granular systems (Otsuki and Hayakawa, 2011) Brownian suspensions (R. Mari, et al, 2015) Boundary effect (Brown and Jaeger, 2012) and (Allen, et al, 2017) (Brown s group) the contributions beneath the DST is still unclear Motivations: Recover DST using LBM Decompose all contributions to the shear stress Helpful for the theoretical construction 7/38
22 Lattice Boltzmann vs Stokesian Dynamics Lattice Boltzmann Stokesian Dynamics Particle motion obey the Newton equations. Hydrodynamic fields are calculated by solving the resistance matrix Separate lubrication calculation. 8/38
23 Lattice Boltzmann vs Stokesian Dynamics Lattice Boltzmann (susp3d) Stokesian Dynamics Particle motion obey the Newton equations. Hydrodynamic fields are calculated by solving the resistance matrix Separate lubrication calculation. Hydrodynamic fields are calculated locally at each lattice point. mesoscopic, possible to have simple local rules between fluid and solid Particle motion obey the Newton equations Separate lubrication calculation. 8/38
24 Outline of our model Equation of motions: m d ( ) U = ( ) Fα dt Ω T α α F α = F h + F c + F R α T α = T h + T c + T R α U(t + t) = U(t) + t m F(t) Ω(t + t) = Ω(t) + t m T(t) 9/38
25 Outline of our model Equation of motions: m d ( ) U = ( ) Fα dt Ω T α α F α = F h + F c + F R α Hydrodynamics (lattice Boltzmann) T α = T h + T c + T R α U(t + t) = U(t) + t m F(t) Ω(t + t) = Ω(t) + t m T(t) 9/38
26 Outline of our model Equation of motions: m d ( ) U = ( ) Fα dt Ω T α α F α = F h + F c + F R α Hydrodynamics (lattice Boltzmann) Contact between particles T α = T h + T c + T R α U(t + t) = U(t) + t m F(t) Ω(t + t) = Ω(t) + t m T(t) 9/38
27 Outline of our model Equation of motions: m d ( ) U = ( ) Fα dt Ω T α α F α = F h + F c + F R α Hydrodynamics (lattice Boltzmann) Contact between particles T α = T h + T c + T R α U(t + t) = U(t) + t m F(t) Ω(t + t) = Ω(t) + t m T(t) Electrostatic repulsive forces 9/38
28 Outline of our model Equation of motions: m d ( ) U = ( ) Fα dt Ω T α α F α = F h + F c + F R α Hydrodynamics (lattice Boltzmann) Contact between particles T α = T h + T c + T R α U(t + t) = U(t) + t m F(t) Ω(t + t) = Ω(t) + t m T(t) Electrostatic repulsive forces 9/38
29 Hydrodynamics: lattice Boltzmann method Hydrodynamics (lattice Boltzmann) 10/38
30 Hydrodynamics: lattice Boltzmann method Hydrodynamics (lattice Boltzmann) By using LBM, we need to discretize the unit length a (particle radius) to lattice unit x! 10/38
31 Hydrodynamics: lattice Boltzmann method Hydrodynamics (lattice Boltzmann) By using LBM, we need to discretize the unit length a (particle radius) to lattice unit x! Our simulation used x = 0.1a 10/38
32 Hydrodynamics: lattice Boltzmann method Hydrodynamics (lattice Boltzmann) By using LBM, we need to discretize the unit length a (particle radius) to lattice unit x! Our simulation used x = 0.1a All quantities below are written in lattice units! 10/38
33 History of LBM Initially devised as an extension of the Lattice Gas Automata by McNamara and Zanetti (McNamara, Zanetti, 1988) 11/38
34 History of LBM Initially devised as an extension of the Lattice Gas Automata by McNamara and Zanetti (McNamara, Zanetti, 1988) Now, LBM is widely used for various computational fluid dynamics simulation. 11/38
35 The discrete distribution function (He and Luo, 1997) n particle velocity (v) distribution function on space time point (r, t) 12/38
36 The discrete distribution function (He and Luo, 1997) n particle velocity (v) distribution function on space time point (r, t) time evolution of n in continuous space: ( ) dn t n + v n = dt coll 12/38
37 The discrete distribution function (He and Luo, 1997) n particle velocity (v) distribution function on space time point (r, t) time evolution of n in continuous space: ( ) dn t n + v n = dt coll we can project to the Hermite bases: n(r, v, t) = ω(v) l=0 1 l! A(l) H (l) (v) 12/38
38 The discrete distribution function (He and Luo, 1997) n particle velocity (v) distribution function on space time point (r, t) time evolution of n in continuous space: ( ) dn t n + v n = dt coll we can project to the Hermite bases: n(r, v, t) = ω(v) l=0 1 l! A(l) H (l) (v) the expansion coefficients: A (l) (r, t) = dvn(r, v, t)h (l) (v) 12/38
39 The discrete distribution function (He and Luo, 1997) n particle velocity (v) distribution function on space time point (r, t) time evolution of n in continuous space: ( ) dn t n + v n = dt coll we can project to the Hermite bases: n(r, v, t) = ω(v) l=0 1 l! A(l) H (l) (v) the expansion coefficients: A (l) (r, t) = dvn(r, v, t)h (l) (v) which is linear combination of the moments of n The first few of the Hermite polynomials (H (l) (v)): H (0) (v) = 1 H (1) (v) = v α H (2) (v) = v α v β c 2 s δ αβ the Gaussian weight function: ( ) 1 ω(v) = exp v2 (2π) D/2 2 12/38
40 The discrete distribution function (He and Luo, 1997) n particle velocity (v) distribution function on space time point (r, t) time evolution of n in continuous space: ( ) dn t n + v n = dt coll we can project to the Hermite bases: n(r, v, t) = ω(v) l=0 1 l! A(l) H (l) (v) the expansion coefficients: A (l) (r, t) = dvn(r, v, t)h (l) (v) which is linear combination of the moments of n The first few of the Hermite polynomials (H (l) (v)): H (0) (v) = 1 H (1) (v) = v α H (2) (v) = v α v β c 2 s δ αβ the Gaussian weight function: ( ) 1 ω(v) = exp v2 (2π) D/2 2 12/38
41 From Boltzmann equation to LBE Truncate: (to several orders of K) n(r, v, t) = ω(v) K l=0 1 l! A(l) H (l) (v) 13/38
42 From Boltzmann equation to LBE Truncate: (to several orders of K) n(r, v, t) = ω(v) K l=0 only finite sets of velocities are needed! 1 l! A(l) H (l) (v) 13/38
43 From Boltzmann equation to LBE Truncate: (to several orders of K) n(r, v, t) = ω(v) K l=0 1 l! A(l) H (l) (v) only finite sets of velocities are needed! we can compute the integral for A with Gauss-Hermite quadrature. (p(v) arbitrary polynomial of v, c i discrete lattice velocity, w i quadrature weight ) dvω(v)p(v) = b w i p(c i ) i=1 13/38
44 From Boltzmann equation to LBE Truncate: (to several orders of K) n(r, v, t) = ω(v) K l=0 1 l! A(l) H (l) (v) only finite sets of velocities are needed! we can compute the integral for A with Gauss-Hermite quadrature. (p(v) arbitrary polynomial of v, c i discrete lattice velocity, w i quadrature weight ) dvω(v)p(v) = b w i p(c i ) i=1 dvn(r, v, t)v v = dv ω(v) n(r, v, t)v v ω(v) b n(r, v, t) = w i c i c i ω(v) = i=1 b n i (r, c i, t)c i c i i=1 13/38
45 From Boltzmann equation to LBE Truncate: (to several orders of K) n(r, v, t) = ω(v) K l=0 1 l! A(l) H (l) (v) only finite sets of velocities are needed! we can compute the integral for A with Gauss-Hermite quadrature. (p(v) arbitrary polynomial of v, c i discrete lattice velocity, w i quadrature weight ) dvω(v)p(v) = b w i p(c i ) i=1 dvn(r, v, t)v v = dv ω(v) n(r, v, t)v v ω(v) b n(r, v, t) = w i c i c i ω(v) = i=1 b n i (r, c i, t)c i c i discrete distribution function n i (r, c i, t) = w i n(r, v, t)/ω(v) i=1 13/38
46 the lattice Boltzmann method Evolving equation: n i (r + c i, t + 1) = n i (r, t) + i (r, t) 14/38
47 the lattice Boltzmann method Evolving equation: n i (r + c i, t + 1) = n i (r, t) + i (r, t) streaming 14/38
48 the lattice Boltzmann method Evolving equation: n i (r + c i, t + 1) = n i (r, t) + i (r, t) streaming collision 14/38
49 the lattice Boltzmann method Evolving equation: n i (r + c i, t + 1) = n i (r, t) + i (r, t) streaming collision Hydrodynamic fields: mass density ρ = i n i momentum density j = i n i c i momentum flux Π = i n i c i c i 14/38
50 LBM: Collision operator Linearized collision operator i (n) = i (n eq i ) + j L ij n neq j 15/38
51 LBM: Collision operator Linearized collision operator i (n) = i (n eq i ) + j L ij n neq j First, we find the discrete equilibrium distribution function n eq i, 15/38
52 LBM: Collision operator Linearized collision operator i (n) = i (n eq i ) + j L ij n neq j First, we find the discrete equilibrium distribution function n eq i, Maxwell-Boltzmann distribution: ( ) 3 ( ) n eq m 2 m(v u) (v u) = ρ exp 2πk b T 2k b T 15/38
53 LBM: Collision operator Linearized collision operator i (n) = i (n eq i ) + j L ij n neq j First, we find the discrete equilibrium distribution function n eq i, Maxwell-Boltzmann distribution: ( ) 3 ( ) n eq m 2 m(v u) (v u) = ρ exp 2πk b T 2k b T Expanding with Hermite polynomials up to 2nd order, with coefficients: A 0 eq = ρ A 1 eq = ρu A 2 eq = ρ(u c 2 s ) 15/38
54 LBM: Collision operator Linearized collision operator i (n) = i (n eq i ) + j L ij n neq j First, we find the discrete equilibrium distribution function n eq i, Maxwell-Boltzmann distribution: ( ) 3 ( ) n eq m 2 m(v u) (v u) = ρ exp 2πk b T 2k b T Expanding with Hermite polynomials up to 2nd order, with coefficients: A 0 eq = ρ A 1 eq = ρu A 2 eq = ρ(u c 2 s ) discrete equilibrium distribution function: ( ) n eq i = a c i ρ+ j c i + (ρuu) : (c ic i cs 2 1) cs 2 2cs 4 15/38
55 LBM: Collision operator Linearized collision operator i (n) = i (n eq i ) + j L ij n neq j First, we find the discrete equilibrium distribution function n eq i, Maxwell-Boltzmann distribution: ( ) 3 ( ) n eq m 2 m(v u) (v u) = ρ exp 2πk b T 2k b T Expanding with Hermite polynomials up to 2nd order, with coefficients: A 0 eq = ρ A 1 eq = ρu A 2 eq = ρ(u c 2 s ) discrete equilibrium distribution function: ( ) n eq i = a c i ρ+ j c i + (ρuu) : (c ic i cs 2 1) cs 2 2cs 4 Weight coefficients: a 0 = 12 a 1 = 2 a 2 = 1 c s = k b T lattice sound speed A : B = Tr(AB) 15/38
56 LBM: Collision operator Linearized collision operator i (n) = i (n eq ) + j L ij n neq j 16/38
57 LBM: Collision operator Linearized collision operator i (n) = i (n eq ) + j L ij n neq j not necessary to construct and calculate L ij, 16/38
58 LBM: Collision operator Linearized collision operator i (n) = i (n eq ) + j L ij n neq j not necessary to construct and calculate L ij, use its eigen equation instead! L ij = 0 c i L ij = 0 c i c i L ij = λc j c j ci 2 L ij = λ ν cj 2 i i i i c j c j traceless 16/38
59 LBM: Collision operator Linearized collision operator i (n) = i (n eq ) + j L ij n neq j not necessary to construct and calculate L ij, use its eigen equation instead! L ij = 0 c i L ij = 0 c i c i L ij = λc j c j ci 2 L ij = λ ν cj 2 i c j c j traceless i discrete equilibrium distribution function: ( n eq i = a c i ρ + j c i c 2 s i + (ρuu) : (c ic i c 2 s 1) 2c 4 s ) i 16/38
60 LBM: Collision operator Linearized collision operator i (n) = i (n eq ) + j L ij n neq j not necessary to construct and calculate L ij, use its eigen equation instead! L ij = 0 c i L ij = 0 c i c i L ij = λc j c j ci 2 L ij = λ ν cj 2 i c j c j traceless i discrete equilibrium distribution function: ( n eq i = a c i Only the second moment of n eq i ρ + j c i c 2 s i + (ρuu) : (c ic i c 2 s 1) 2c 4 s that is affected by collision ) i 16/38
61 LBM: Collision operator Linearized collision operator i (n) = i (n eq ) + j L ij n neq j not necessary to construct and calculate L ij, use its eigen equation instead! L ij = 0 c i L ij = 0 c i c i L ij = λc j c j ci 2 L ij = λ ν cj 2 i c j c j traceless i discrete equilibrium distribution function: ( n eq i = a c i ρ + j c i c 2 s i + (ρuu) : (c ic i c 2 s 1) 2c 4 s Only the second moment of n eq i that is affected by collision Post collision distribution function: n i = a c i ( ρ + j c i c 2 s ) + (ρuu + Πneq, ) : (c i c i c 2 s 1) 2c 4 s i ) 16/38
62 LBM: Collision operator Linearized collision operator i (n) = i (n eq ) + j L ij n neq j not necessary to construct and calculate L ij, use its eigen equation instead! L ij = 0 c i L ij = 0 c i c i L ij = λc j c j ci 2 L ij = λ ν cj 2 i i i i c j c j traceless discrete equilibrium distribution function: ( n eq i = a c i ρ + j c i c 2 s + (ρuu) : (c ic i c 2 s 1) 2c 4 s Only the second moment of n eq i that is affected by collision Post collision distribution function: n i = a c i ( ρ + j c i c 2 s ) + (ρuu + Πneq, ) : (c i c i c 2 s 1) 2c 4 s Calculate Π neq, ) 16/38
63 LBM: Collision operator Π neq, Π neq = Π Π eq Π = i n ic i c i 17/38
64 LBM: Collision operator Π neq, Π neq = Π Π eq Π = i n ic i c i nonequilibrium second moments obtained from the eigen equation of L ij : L ij = 0 c i L ij = 0 c i c i L ij = λc j c j ci 2 L ij = λ ν cj 2 i i i i 17/38
65 LBM: Collision operator Π neq, Π neq = Π Π eq Π = i n ic i c i nonequilibrium second moments obtained from the eigen equation of L ij : L ij = 0 c i L ij = 0 c i c i L ij = λc j c j ci 2 L ij = λ ν cj 2 i i i i Π neq, = (1 + λ) Π neq (1 + λ ν)(π neq : 1)1 17/38
66 LBM: Collision operator Π neq, Π neq = Π Π eq Π = i n ic i c i nonequilibrium second moments obtained from the eigen equation of L ij : L ij = 0 i c i L ij = 0 i c i c i L ij = λc j c j i Π neq, = (1 + λ) Π neq (1 + λ ν)(π neq : 1)1 λ and λ ν are related to shear η and bulk η ν viscosities from the multiscale analysis: ( 1 η = ρcs 2 t λ + 1 ) ( 2 η ν = ρcs 2 t + 1 ) 2 3λ ν 3 ci 2 L ij = λ ν cj 2 i 17/38
67 LBM: Collision operator Π neq, Π neq = Π Π eq Π = i n ic i c i nonequilibrium second moments obtained from the eigen equation of L ij : L ij = 0 i c i L ij = 0 i c i c i L ij = λc j c j i Π neq, = (1 + λ) Π neq (1 + λ ν)(π neq : 1)1 λ and λ ν are related to shear η and bulk η ν viscosities from the multiscale analysis: ( 1 η = ρcs 2 t λ + 1 ) ( 2 η ν = ρcs 2 t + 1 ) 2 3λ ν 3 What will happen on the solid boundary conditions (surface of the particles)? ci 2 L ij = λ ν cj 2 i 17/38
68 Solid-fluid boundary condition Anthony Ladd s bounce-back rule: n b (r, t + t) = n b (r, t) 2ac bρ 0 u b c b c 2 s Velocity of the boundary nodes: u b = U + Ω (r b R) R is the center of mass of the particle 18/38
69 Solid-fluid boundary condition Forces exerted at the boundary nodes: [ ] f(r b, t+ 1 3 x t) = 2nb (r, t) 2ac b ρ 0u b c b 2 t cs 2 c b 19/38
70 Solid-fluid boundary condition Sum over all boundary nodes within a particle: Forces exerted at the boundary nodes: [ ] f(r b, t+ 1 3 x t) = 2nb (r, t) 2ac b ρ 0u b c b 2 t cs 2 c b 19/38
71 Solid-fluid boundary condition Sum over all boundary nodes within a particle: Forces exerted at the boundary nodes: [ ] f(r b, t+ 1 3 x t) = 2nb (r, t) 2ac b ρ 0u b c b 2 t cs 2 c b F h = b T h = b σ h = b f(r b ) r b f(r b ) r b f(r b ) 19/38
72 Shared nodes 20/38
73 Shared nodes need separated lubrication forces calculation! 20/38
74 Lubrications If gap between particle is less than 1 lattice unit.. 21/38
75 Lubrications If gap between particle is less than 1 lattice unit.. Grand-resistance formulation (Nguyen and Ladd, 2002) 21/38
76 Lubrications If gap between particle is less than 1 lattice unit.. Grand-resistance formulation (Nguyen and Ladd, 2002) F 1 A 11 B 11 B 22 T 1 B 11 C 11 C 12 T 2 S 1 S 2 = B 22 C 12 C 22 G 11 H 11 H 12 G 22 H 21 H 22 (Kim and Karilla, 1991) U 12 Ω 1 Ω 2 21/38
77 Lubrications If gap between particle is less than 1 lattice unit.. Grand-resistance formulation (Nguyen and Ladd, 2002) F 1 A 11 B 11 B 22 T 1 B 11 C 11 C 12 T 2 S 1 S 2 U 12 = U 1 U 2 relative velocity = B 22 C 12 C 22 G 11 H 11 H 12 G 22 H 21 H 22 (Kim and Karilla, 1991) U 12 Ω 1 Ω 2 21/38
78 Lubrications If gap between particle is less than 1 lattice unit.. Grand-resistance formulation (Nguyen and Ladd, 2002) F 1 A 11 B 11 B 22 T 1 B 11 C 11 C 12 T 2 S 1 S 2 U 12 = U 1 U 2 relative velocity = B 22 C 12 C 22 G 11 H 11 H 12 G 22 H 21 H 22 (Kim and Karilla, 1991) U 12 Ω 1 Ω 2 Each coefficients can be expressed in terms of scalar function i.e 21/38
79 Lubrications If gap between particle is less than 1 lattice unit.. Grand-resistance formulation (Nguyen and Ladd, 2002) F 1 A 11 B 11 B 22 T 1 B 11 C 11 C 12 T 2 S 1 S 2 U 12 = U 1 U 2 relative velocity = B 22 C 12 C 22 G 11 H 11 H 12 G 22 H 21 H 22 (Kim and Karilla, 1991) U 12 Ω 1 Ω 2 Each coefficients can be expressed in terms of scalar function i.e H 12 = Y H 12(ɛ αγδ d δ d β + ɛ βγδ d δ d α ) d displacement unit vector along axis. ɛ Levi-Civita symbol 21/38
80 Lubrications If gap between particle is less than 1 lattice unit.. Grand-resistance formulation (Nguyen and Ladd, 2002) F 1 A 11 B 11 B 22 T 1 B 11 C 11 C 12 T 2 S 1 S 2 U 12 = U 1 U 2 relative velocity = B 22 C 12 C 22 G 11 H 11 H 12 G 22 H 21 H 22 (Kim and Karilla, 1991) U 12 Ω 1 Ω 2 Each coefficients can be expressed in terms of scalar function i.e H 12 = Y H 12(ɛ αγδ d δ d β + ɛ βγδ d δ d α ) d displacement unit vector along axis. ɛ Levi-Civita symbol Each scalar function is a function of gap h and β = a i a j ( ) 1 Y12 H 2β 2 (1 + 7β) = 8 log πηa i h 5(1 + β) 5 i.e 21/38
81 Lubrications Cutoff length δ = acontact a hydro a contact to allow contact 22/38
82 Lubrications Cutoff length δ = acontact a hydro a contact to allow contact We used δ = /38
83 Particle contacts 23/38
84 Contact model Linear spring dashpot model (Luding, 2008) 24/38
85 Contact model Linear spring dashpot model (Luding, 2008) F c ij = Fnor ij + F tan ij (Fleischmann, 2015) 24/38
86 Contact model Coulomb friction rules: F tan ij µ( F nor ij ) 25/38
87 Contact model Coulomb friction rules: F tan ij µ( F nor ij ) slip 25/38
88 Contact model Coulomb friction rules: F tan ij µ( F nor ij ) slip F tan ij µ( F nor ij ) 25/38
89 Contact model Coulomb friction rules: F tan ij µ( F nor ij ) slip F tan ij µ( F nor ij ) stick 25/38
90 Contact model Coulomb friction rules: F tan ij µ( F nor ij ) slip F tan ij µ( F nor ij ) stick Stress contribution from contact: Normal: σ nor αβ = 1 2V i j i (r ij,α F nor ij,β + r ij,βf nor ij,α) 25/38
91 Contact model Coulomb friction rules: F tan ij µ( F nor ij ) slip F tan ij µ( F nor ij ) stick Stress contribution from contact: Normal: σ nor αβ = 1 2V i j i (r ij,α F nor ij,β + r ij,βf nor ij,α) Tangential: σ tan αβ = 1 V i j i r ij,α F tan ij,β 25/38
92 Contact model Coulomb friction rules: F tan ij µ( F nor ij ) slip F tan ij µ( F nor ij ) stick Stress contribution from contact: Normal: σ nor αβ = 1 2V i j i (r ij,α F nor ij,β + r ij,βf nor ij,α) Tangential: σ tan αβ = 1 V i j i r ij,α F tan ij,β 25/38
93 Electrostatic repulsive forces (Israelachvili, 2001) 26/38
94 Electrostatic repulsive forces (Israelachvili, 2001) ( ) F R ij = 1 a i a j λ a i + a j F exp( h/λ)ˆn ij λ Debye length 26/38
95 Electrostatic repulsive forces (Israelachvili, 2001) ( ) F R ij = 1 a i a j λ a i + a j F exp( h/λ)ˆn ij λ Debye length we used λ = 0.2a 26/38
96 Electrostatic repulsive forces (Israelachvili, 2001) ( ) F R ij = 1 a i a j λ a i + a j F exp( h/λ)ˆn ij λ Debye length we used λ = 0.2a Stress contribution: σ R αβ = 1 V R ij,α F R ij,β i j i 26/38
97 Sheared suspensions 27/38
98 Sheared suspensions wall moves with velocity u wall to x and x directions All simulation use N=512 particles. 27/38
99 Rheology Shear stress: σ αβ = σ h αβ + σtan αβ + σnor αβ + σr αβ 28/38
100 Rheology Shear stress: σ αβ = σ h αβ + σtan αβ + σnor αβ + σr αβ Apparent viscosity: η = σ αβ / γ 28/38
101 Rheology Shear stress: σ αβ = σ h αβ + σtan αβ + σnor αβ + σr αβ Apparent viscosity: η = σ αβ / γ Dimensionless shear rate: γ = 6πη 0a 2 γ F 28/38
102 Rheology Shear stress: σ αβ = σ h αβ + σtan αβ + σnor αβ + σr αβ Apparent viscosity: η = σ αβ / γ Dimensionless shear rate: γ = 6πη 0a 2 γ F 28/38
103 Sheared Suspensions: (Mari, et al, 2014) η η 0 vs γ, compared with results from 29/38
104 Sheared Suspensions: (Mari, et al, 2014) η η 0 vs γ, compared with results from 29/38
105 η η 0 Sheared Suspensions: vs σ/η 0 γ, compared with results from (Mari, et al, 2014) 30/38
106 η η 0 Sheared Suspensions: vs σ/η 0 γ, compared with results from (Mari, et al, 2014) 30/38
107 Contributions to shear viscosity for φ = /38
108 Contributions to shear viscosity for φ = /38
109 Contributions to shear viscosity for φ = 0.54 and /38
110 Contributions to shear viscosity for φ = 0.54 and /38
111 Time evolution of shear viscosity φ = 0.57 and φ = /38
112 Time evolution of shear viscosity φ = 0.57 and φ = /38
113 Time evolution of contact number for φ = /38
114 Time evolution of contact number for φ = /38
115 Evolution of the contact network (φ = 0.57) high shear rate low shear rate 35/38
116 Discussions Slight deviation from Seto s 36/38
117 Discussions Slight deviation from Seto s Need to allow separate switching length between boundary nodes and lubrication calculation for normal, tangential, and rotation modes 36/38
118 Discussions Slight deviation from Seto s Need to allow separate switching length between boundary nodes and lubrication calculation for normal, tangential, and rotation modes We need to find a suitable contact spring constant for every φ 36/38
119 Discussions Slight deviation from Seto s Need to allow separate switching length between boundary nodes and lubrication calculation for normal, tangential, and rotation modes We need to find a suitable contact spring constant for every φ Contact force dominates the shear thickening regime 36/38
120 Discussions Slight deviation from Seto s Need to allow separate switching length between boundary nodes and lubrication calculation for normal, tangential, and rotation modes We need to find a suitable contact spring constant for every φ Contact force dominates the shear thickening regime Number of contact increased on high shear rate, creates force chains from one side to another. 36/38
121 Discussions Slight deviation from Seto s Need to allow separate switching length between boundary nodes and lubrication calculation for normal, tangential, and rotation modes We need to find a suitable contact spring constant for every φ Contact force dominates the shear thickening regime Number of contact increased on high shear rate, creates force chains from one side to another. Electsrotatic repulsive force shear thinning at low shear rate 36/38
122 Discussions Slight deviation from Seto s Need to allow separate switching length between boundary nodes and lubrication calculation for normal, tangential, and rotation modes We need to find a suitable contact spring constant for every φ Contact force dominates the shear thickening regime Number of contact increased on high shear rate, creates force chains from one side to another. Electsrotatic repulsive force shear thinning at low shear rate The computational cost highly depends on the size. 36/38
123 Discussions Slight deviation from Seto s Need to allow separate switching length between boundary nodes and lubrication calculation for normal, tangential, and rotation modes We need to find a suitable contact spring constant for every φ Contact force dominates the shear thickening regime Number of contact increased on high shear rate, creates force chains from one side to another. Electsrotatic repulsive force shear thinning at low shear rate The computational cost highly depends on the size. Our simulation used x = 0.1a, to let Re < 1, one need typical particle speed U < x t. 36/38
124 Discussions Slight deviation from Seto s Need to allow separate switching length between boundary nodes and lubrication calculation for normal, tangential, and rotation modes We need to find a suitable contact spring constant for every φ Contact force dominates the shear thickening regime Number of contact increased on high shear rate, creates force chains from one side to another. Electsrotatic repulsive force shear thinning at low shear rate The computational cost highly depends on the size. Our simulation used x = 0.1a, to let Re < 1, one need typical particle speed U < x t. To travel a distance of its radius, one particles needs > 6200 time steps 36/38
125 Future prospects Good starting point for the theoretical works 37/38
126 Future prospects Good starting point for the theoretical works friction scenario. 37/38
127 Future prospects Good starting point for the theoretical works friction scenario. Analyze the percolation of the contact network 37/38
128 Future prospects Good starting point for the theoretical works friction scenario. Analyze the percolation of the contact network Implement contact with rolling friction. 37/38
129 Future prospects Good starting point for the theoretical works friction scenario. Analyze the percolation of the contact network Implement contact with rolling friction. Use more realistic boundary conditions 37/38
130 References Ladd, A.J.C Numerical simulations of particulate suspensions via a discretized Boltzmann equation. Part 1. Theoretical foundation, Journal of fluid mechanics, vol 271, , Ladd, A.J.C Numerical simulations of particulate suspensions via a discretized Boltzmann equation. Part 2. Numerical Results, Journal of fluid mechanics, vol 271, , N.Q. Nguyen and Ladd, A.J.C, Lubrication corrections for lattice-boltzmann simulations of particle suspensions. Physical Review E 66, , Ladd, A.J.C, Susp3D, an open source code for simulating suspensions based on Lattice Boltzmann Method, Romain Mari, Ryohei Seto, Jeffrey F. Morris and Morton M. Denn, Shear thickening, frictionless and frictional rheologies in non-brownian suspensions, Journal of Rheology 58(6), , Stefan Luding, Cohesive, frictional powders: contact models for tension, Granular Matter 10:235246, Ryohei Seto,Romain Mari,Jeffrey F. Morris, and Morton M. Denn, Discontinuous Shear Thickening of Frictional Hard-Sphere Suspensions, Physical Review Letters 111, , S. Kim and S. J. Karrila, Microhydrodynamics: Principles and Selected Applications, Butterworth-Heinemann, Boston, MA, /38
Discontinuous Shear Thickening
Discontinuous Shear Thickening dynamic jamming transition Ryohei Seto, Romain Mari, Jeffrey F. Morris, Morton M. Denn Levich Institute, City College of New York First experimental data Williamson and Hecker
More informationA simulation study on shear thickening in wide-gap Couette geometry. Ryohei Seto, Romain Mari Jeffery Morris, Morton Denn, Eliot Fried
A simulation study on shear thickening in wide-gap Couette geometry Ryohei Seto, Romain Mari Jeffery Morris, Morton Denn, Eliot Fried Stokes flow: Zero-Reynolds number fluid mechanics Repulsive Attractive
More informationShear-induced organization of forces in dense suspensions: signatures of discontinuous shear thickening
Shear-induced organization of forces in dense suspensions: signatures of discontinuous shear thickening The essential idea underlying the jamming-based rheology model [3, 6] for DST is that it requires
More informationShear-induced organization of forces in dense suspensions: signatures of discontinuous shear thickening
Shear-induced organization of forces in dense suspensions: signatures of discontinuous shear thickening The MIT Faculty has made this article openly available. Please share how this access benefits you.
More informationEmergence of collective dynamics in active biological systems -- Swimming micro-organisms --
12/08/2015, YITP, Kyoto Emergence of collective dynamics in active biological systems -- Swimming micro-organisms -- Norihiro Oyama John J. Molina Ryoichi Yamamoto* Department of Chemical Engineering,
More informationMicromechanical Modeling of Discontinuous Shear Thickening in Granular Media-Fluid
1 2 Micromechanical Modeling of Discontinuous Shear Thickening in Granular Media-Fluid Suspension 3 4 Daniel H. Johnson a, Farshid Vahedifard b, Bohumir Jelinek c, John F. Peters d 5 6 7 8 9 10 11 12 13
More informationarxiv: v2 [cond-mat.soft] 3 Dec 2015
Discontinuous Shear Thickening in Brownian Suspensions by Dynamic Simulation Romain Mari Benjamin Levich Institute, City College of New York, New York, NY 10031, USA Ryohei Seto Mathematical Soft Matter
More informationFluid Equations for Rarefied Gases
1 Fluid Equations for Rarefied Gases Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc 23 March 2001 with E. A. Spiegel
More informationSEMICLASSICAL LATTICE BOLTZMANN EQUATION HYDRODYNAMICS
The Seventh Nobeyama Workshop on CFD: To the Memory of Prof. Kuwahara Sanjo Kaikan, the University of Tokyo, Japan, September. 23-24, 2009 SEMICLASSICAL LATTICE BOLTZMANN EQUATION HYDRODYNAMICS Jaw-Yen
More informationOn the bulk viscosity of suspensions
J. Fluid Mech. (2006), vol. 554, pp. 109 12. c 2006 Cambridge University Press doi:10.1017/s002211200600948 Printed in the United Kingdom 109 On the bulk viscosity of suspensions By JOHN F. BRADY, ADITYA
More informationResearch of Micro-Rectangular-Channel Flow Based on Lattice Boltzmann Method
Research Journal of Applied Sciences, Engineering and Technology 6(14): 50-55, 013 ISSN: 040-7459; e-issn: 040-7467 Maxwell Scientific Organization, 013 Submitted: November 08, 01 Accepted: December 8,
More informationLattice Boltzmann Method
3 Lattice Boltzmann Method 3.1 Introduction The lattice Boltzmann method is a discrete computational method based upon the lattice gas automata - a simplified, fictitious molecular model. It consists of
More informationFluid equations, magnetohydrodynamics
Fluid equations, magnetohydrodynamics Multi-fluid theory Equation of state Single-fluid theory Generalised Ohm s law Magnetic tension and plasma beta Stationarity and equilibria Validity of magnetohydrodynamics
More informationGeneralized Local Equilibrium in the Cascaded Lattice Boltzmann Method. Abstract
Accepted for publication on Physical Review E (R), code ERR1034 Generalized Local Equilibrium in the Cascaded Lattice Boltzmann Method Pietro Asinari Department of Energetics, Politecnico di Torino, Corso
More informationOrigins of Mechanical and Rheological Properties of Polymer Nanocomposites. Venkat Ganesan
Department of Chemical Engineering University of Texas@Austin Origins of Mechanical and Rheological Properties of Polymer Nanocomposites Venkat Ganesan $$$: NSF DMR, Welch Foundation Megha Surve, Victor
More informationLATTICE BOLTZMANN MODELLING OF PULSATILE FLOW USING MOMENT BOUNDARY CONDITIONS
6th European Conference on Computational Mechanics (ECCM 6) 7th European Conference on Computational Fluid Dynamics (ECFD 7) 5 June 28, Glasgow, UK LATTICE BOLTZMANN MODELLING OF PULSATILE FLOW USING MOMENT
More informationSmoothed Dissipative Particle Dynamics Model for Predicting Self-Assembled Nano-Cellulose Fibre Structures
Smoothed Dissipative Particle Dynamics Model for Predicting Self-Assembled Nano-Cellulose Fibre Structures David Vidal and Tetsu Uesaka FPInnovations, Pointe-Claire, Québec, CANADA Nano-cellulose fibres
More informationFluid Equations for Rarefied Gases
1 Fluid Equations for Rarefied Gases Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc 21 May 2001 with E. A. Spiegel
More information(Crystal) Nucleation: The language
Why crystallization requires supercooling (Crystal) Nucleation: The language 2r 1. Transferring N particles from liquid to crystal yields energy. Crystal nucleus Δµ: thermodynamic driving force N is proportional
More informationBROWNIAN DYNAMICS SIMULATIONS WITH HYDRODYNAMICS. Abstract
BROWNIAN DYNAMICS SIMULATIONS WITH HYDRODYNAMICS Juan J. Cerdà 1 1 Institut für Computerphysik, Pfaffenwaldring 27, Universität Stuttgart, 70569 Stuttgart, Germany. (Dated: July 21, 2009) Abstract - 1
More informationLattice Boltzmann approach to liquid - vapour separation
Lattice Boltzmann approach to liquid - vapour separation T.Biciușcă 1,2, A.Cristea 1, A.Lamura 3, G.Gonnella 4, V.Sofonea 1 1 Centre for Fundamental and Advanced Technical Research, Romanian Academy Bd.
More informationLecture 8: Tissue Mechanics
Computational Biology Group (CoBi), D-BSSE, ETHZ Lecture 8: Tissue Mechanics Prof Dagmar Iber, PhD DPhil MSc Computational Biology 2015/16 7. Mai 2016 2 / 57 Contents 1 Introduction to Elastic Materials
More informationSimulation of floating bodies with lattice Boltzmann
Simulation of floating bodies with lattice Boltzmann by Simon Bogner, 17.11.2011, Lehrstuhl für Systemsimulation, Friedrich-Alexander Universität Erlangen 1 Simulation of floating bodies with lattice Boltzmann
More informationPhysical Modeling of Multiphase flow. Boltzmann method
with lattice Boltzmann method Exa Corp., Burlington, MA, USA Feburary, 2011 Scope Scope Re-examine the non-ideal gas model in [Shan & Chen, Phys. Rev. E, (1993)] from the perspective of kinetic theory
More informationLattice Boltzmann Method for Moving Boundaries
Lattice Boltzmann Method for Moving Boundaries Hans Groot March 18, 2009 Outline 1 Introduction 2 Moving Boundary Conditions 3 Cylinder in Transient Couette Flow 4 Collision-Advection Process for Moving
More informationTheory of Force Transmission and Rigidity in Granular Aggregates. A Different kind of NonEquilibrium System
Theory of Force Transmission and Rigidity in Granular Aggregates A Different kind of NonEquilibrium System Dapeng Bi Sumantra Sarkar Kabir Ramola Jetin Thomas Bob Behringer, Jie Ren Dong Wang Jeff Morris
More informationViscoelasticity. Basic Notions & Examples. Formalism for Linear Viscoelasticity. Simple Models & Mechanical Analogies. Non-linear behavior
Viscoelasticity Basic Notions & Examples Formalism for Linear Viscoelasticity Simple Models & Mechanical Analogies Non-linear behavior Viscoelastic Behavior Generic Viscoelasticity: exhibition of both
More informationMultiscale modeling of active fluids: selfpropellers and molecular motors. I. Pagonabarraga University of Barcelona
Multiscale modeling of active fluids: selfpropellers and molecular motors I. Pagonabarraga University of Barcelona Introduction Soft materials weak interactions Self-assembly Emergence large scale structures
More informationarxiv: v1 [nucl-th] 9 Jun 2008
Dissipative effects from transport and viscous hydrodynamics arxiv:0806.1367v1 [nucl-th] 9 Jun 2008 1. Introduction Denes Molnar 1,2 and Pasi Huovinen 1 1 Purdue University, Physics Department, 525 Northwestern
More informationSummary PHY101 ( 2 ) T / Hanadi Al Harbi
الكمية Physical Quantity القانون Low التعريف Definition الوحدة SI Unit Linear Momentum P = mθ be equal to the mass of an object times its velocity. Kg. m/s vector quantity Stress F \ A the external force
More informationCHAPTER 9. Microscopic Approach: from Boltzmann to Navier-Stokes. In the previous chapter we derived the closed Boltzmann equation:
CHAPTER 9 Microscopic Approach: from Boltzmann to Navier-Stokes In the previous chapter we derived the closed Boltzmann equation: df dt = f +{f,h} = I[f] where I[f] is the collision integral, and we have
More informationEFFECT OF A CLUSTER ON GAS-SOLID DRAG FROM LATTICE BOLTZMANN SIMULATIONS
Ninth International Conference on CFD in the Minerals and Process Industries CSIRO, Melbourne, Australia 10-12 December 2012 EFFECT OF A CLUSTER ON GAS-SOLID DRAG FROM LATTICE BOLTZMANN SIMULATIONS Milinkumar
More informationarxiv:comp-gas/ v1 28 Apr 1993
Lattice Boltzmann Thermohydrodynamics arxiv:comp-gas/9304006v1 28 Apr 1993 F. J. Alexander, S. Chen and J. D. Sterling Center for Nonlinear Studies and Theoretical Division Los Alamos National Laboratory
More informationFrictional shear thickening in suspensions: The effect of rigid asperities.
Frictional shear thickening in suspensions: The effect of rigid asperities Adam K Townsend and Helen J Wilson Department of Mathematics, University College London, Gower Street, London WCE 6BT (Dated:
More informationThe lattice Boltzmann method for contact line dynamics
The lattice Boltzmann method for contact line dynamics Sudhir Srivastava, J.H.M. ten Thije Boonkkamp, Federico Toschi April 13, 2011 Overview 1 Problem description 2 Huh and Scriven model 3 Lattice Boltzmann
More informationShear thickening in highly viscous granular suspensions
OFFPRINT Shear thickening in highly viscous granular suspensions Qin Xu, Sayantan Majumdar, Eric Brown and Heinrich M. Jaeger EPL, 07 (204) 68004 Please visit the website www.epljournal.org Note that the
More information4. The Green Kubo Relations
4. The Green Kubo Relations 4.1 The Langevin Equation In 1828 the botanist Robert Brown observed the motion of pollen grains suspended in a fluid. Although the system was allowed to come to equilibrium,
More informationModeling of Suspension Flow in Pipes and Rheometers
Modeling of Suspension Flow in Pipes and Rheometers Nicos S. Martys, Chiara F. Ferraris, William L. George National Institute of Standards and Technology Abstract: Measurement and prediction of the flow
More informationFrom a Mesoscopic to a Macroscopic Description of Fluid-Particle Interaction
From a Mesoscopic to a Macroscopic Description of Fluid-Particle Interaction Carnegie Mellon University Center for Nonlinear Analysis Working Group, October 2016 Outline 1 Physical Framework 2 3 Free Energy
More informationLecture 5: Kinetic theory of fluids
Lecture 5: Kinetic theory of fluids September 21, 2015 1 Goal 2 From atoms to probabilities Fluid dynamics descrines fluids as continnum media (fields); however under conditions of strong inhomogeneities
More informationSimulation of 2D non-isothermal flows in slits using lattice Boltzmann method
Simulation of 2D non-isothermal flows in slits using lattice Boltzmann method João Chambel Leitão Department of Mechanical Engineering, Instituto Superior Técnico, Lisboa, Portugal Abstract: The present
More information12. MHD Approximation.
Phys780: Plasma Physics Lecture 12. MHD approximation. 1 12. MHD Approximation. ([3], p. 169-183) The kinetic equation for the distribution function f( v, r, t) provides the most complete and universal
More information1. The Properties of Fluids
1. The Properties of Fluids [This material relates predominantly to modules ELP034, ELP035] 1.1 Fluids 1.1 Fluids 1.2 Newton s Law of Viscosity 1.3 Fluids Vs Solids 1.4 Liquids Vs Gases 1.5 Causes of viscosity
More informationAnalysis of the trajectory of a sphere moving through a geometric constriction
Analysis of the trajectory of a sphere moving through a geometric constriction Sumedh R. Risbud, Mingxiang Luo, Joëlle Fréchette, and German Drazer Citation: Phys. Fluids 25, 062001 (2013); doi: 1063/1.4809729
More informationMicro-Macro transition (from particles to continuum theory) Granular Materials. Approach philosophy. Model Granular Materials
Micro-Macro transition (from particles to continuum theory) Stefan Luding MSM, TS, CTW, UTwente, NL Stefan Luding, s.luding@utwente.nl MSM, TS, CTW, UTwente, NL Granular Materials Real: sand, soil, rock,
More informationTurbulent Flows. g u
.4 g u.3.2.1 t. 6 4 2 2 4 6 Figure 12.1: Effect of diffusion on PDF shape: solution to Eq. (12.29) for Dt =,.2,.2, 1. The dashed line is the Gaussian with the same mean () and variance (3) as the PDF at
More informationIntroduction to Marine Hydrodynamics
1896 1920 1987 2006 Introduction to Marine Hydrodynamics (NA235) Department of Naval Architecture and Ocean Engineering School of Naval Architecture, Ocean & Civil Engineering First Assignment The first
More informationSimulation of Particulate Solids Processing Using Discrete Element Method Oleh Baran
Simulation of Particulate Solids Processing Using Discrete Element Method Oleh Baran Outline DEM overview DEM capabilities in STAR-CCM+ Particle types and injectors Contact physics Coupling to fluid flow
More informationStrategy in modelling irregular shaped particle behaviour in confined turbulent flows
Title Strategy in modelling irregular shaped particle behaviour in confined turbulent flows M. Sommerfeld F L Mechanische Verfahrenstechnik Zentrum Ingenieurwissenschaften 699 Halle (Saale), Germany www-mvt.iw.uni-halle.de
More informationExternal and Internal Incompressible Viscous Flows Computation using Taylor Series Expansion and Least Square based Lattice Boltzmann Method
Available online at http://ijim.srbiau.ac.ir/ Int. J. Industrial Mathematics (ISSN 2008-5621) Vol. 10, No. 2, 2018 Article ID IJIM-00726, 8 pages Research Article External and Internal Incompressible Viscous
More informationAnomalous effect of turning off long-range mobility interactions in Stokesian Dynamics
Anomalous effect of turning off long-range mobility interactions in Stokesian Dynamics Adam K. Townsend, a), b) and Helen J. Wilson Department of Mathematics, University College London, Gower Street, London
More informationAn electrokinetic LB based model for ion transport and macromolecular electrophoresis
An electrokinetic LB based model for ion transport and macromolecular electrophoresis Raffael Pecoroni Supervisor: Michael Kuron July 8, 2016 1 Introduction & Motivation So far an mesoscopic coarse-grained
More information7 The Navier-Stokes Equations
18.354/12.27 Spring 214 7 The Navier-Stokes Equations In the previous section, we have seen how one can deduce the general structure of hydrodynamic equations from purely macroscopic considerations and
More informationConvection. forced convection when the flow is caused by external means, such as by a fan, a pump, or atmospheric winds.
Convection The convection heat transfer mode is comprised of two mechanisms. In addition to energy transfer due to random molecular motion (diffusion), energy is also transferred by the bulk, or macroscopic,
More informationParticle-Simulation Methods for Fluid Dynamics
Particle-Simulation Methods for Fluid Dynamics X. Y. Hu and Marco Ellero E-mail: Xiangyu.Hu and Marco.Ellero at mw.tum.de, WS 2012/2013: Lectures for Mechanical Engineering Institute of Aerodynamics Technical
More informationPREDICTION OF INTRINSIC PERMEABILITIES WITH LATTICE BOLTZMANN METHOD
PREDICTION OF INTRINSIC PERMEABILITIES WITH LATTICE BOLTZMANN METHOD Luís Orlando Emerich dos Santos emerich@lmpt.ufsc.br Carlos Enrique Pico Ortiz capico@lmpt.ufsc.br Henrique Cesar de Gaspari henrique@lmpt.ufsc.br
More informationConnection Between the Lattice Boltzmann Equation and the Beam Scheme
NASA/CR-1999-209001 ICASE Report No. 99-10 Connection Between the Lattice Boltzmann Equation and the Beam Scheme Kun Xu Hong Kong University of Science and Technology, Kowloon, Hong Kong Li-Shi Luo ICASE,
More informationFree energy concept Free energy approach LBM implementation Parameters
BINARY LIQUID MODEL A. Kuzmin J. Derksen Department of Chemical and Materials Engineering University of Alberta Canada August 22,2011 / LBM Workshop OUTLINE 1 FREE ENERGY CONCEPT 2 FREE ENERGY APPROACH
More informationHydrodynamics, Thermodynamics, and Mathematics
Hydrodynamics, Thermodynamics, and Mathematics Hans Christian Öttinger Department of Mat..., ETH Zürich, Switzerland Thermodynamic admissibility and mathematical well-posedness 1. structure of equations
More informationAndrés Santos Universidad de Extremadura, Badajoz (Spain)
Andrés Santos Universidad de Extremadura, Badajoz (Spain) Outline Moment equations molecules for Maxwell Some solvable states: Planar Fourier flow Planar Fourier flow with gravity Planar Couette flow Force-driven
More information2. Molecules in Motion
2. Molecules in Motion Kinetic Theory of Gases (microscopic viewpoint) assumptions (1) particles of mass m and diameter d; ceaseless random motion (2) dilute gas: d λ, λ = mean free path = average distance
More informationLattice Boltzmann Method for Fluid Simulations
1 / 16 Lattice Boltzmann Method for Fluid Simulations Yuanxun Bill Bao & Justin Meskas Simon Fraser University April 7, 2011 2 / 16 Ludwig Boltzmann and His Kinetic Theory of Gases The Boltzmann Transport
More informationMULTIPHASE DEBRIS FLOW SIMULATIONS WITH THE DISCRETE ELEMENT METHOD COUPLED WITH A LATTICE-BOLTZMANN FLUID
III International Conference on Particle-based Methods Fundamentals and Applications PARTICLES 2013 M. Bischoff, E. Oñate, D.R.J. Owen, E. Ramm & P. Wriggers (Eds) MULTIPHASE DEBRIS FLOW SIMULATIONS WITH
More informationPhysics 221. Exam III Spring f S While the cylinder is rolling up, the frictional force is and the cylinder is rotating
Physics 1. Exam III Spring 003 The situation below refers to the next three questions: A solid cylinder of radius R and mass M with initial velocity v 0 rolls without slipping up the inclined plane. N
More informationMicrofluidic crystals: Impossible order
Microfluidic crystals: Impossible order Tsevi Beatus, Roy Bar-Ziv, T. T. Weizmann Institute International Symposium on Non-Equilibrium Soft Matter Kyoto 2008 1 Outline Micro-fluidic droplets: micron sized
More informationGranular materials (Assemblies of particles with dissipation )
Granular materials (Assemblies of particles with dissipation ) Saturn ring Sand mustard seed Ginkaku-ji temple Sheared granular materials packing fraction : Φ Inhomogeneous flow Gas (Φ = 012) Homogeneous
More informationMechanics Answers to Examples B (Momentum) - 1 David Apsley
TOPIC B: MOMENTUM ANSWERS SPRING 2019 (Full worked answers follow on later pages) Q1. (a) 2.26 m s 2 (b) 5.89 m s 2 Q2. 8.41 m s 2 and 4.20 m s 2 ; 841 N Q3. (a) 1.70 m s 1 (b) 1.86 s Q4. (a) 1 s (b) 1.5
More informationFlow simulation of fiber reinforced self compacting concrete using Lattice Boltzmann method
Flow simulation of fiber reinforced self compacting concrete using Lattice Boltzmann method 1 Oldřich Švec 1 * 1 Technical University of Denmark, Department of Civil Engineering, Lyngby, Denmark 2 Jan
More informationPolymer Dynamics and Rheology
Polymer Dynamics and Rheology 1 Polymer Dynamics and Rheology Brownian motion Harmonic Oscillator Damped harmonic oscillator Elastic dumbbell model Boltzmann superposition principle Rubber elasticity and
More informationAnalysis and boundary condition of the lattice Boltzmann BGK model with two velocity components
Analysis and boundary condition of the lattice Boltzmann BGK model with two velocity components Xiaoyi He and Qisu Zou arxiv:comp-gas/9507002v1 22 Jul 1995 Abstract In this paper, we study the two dimensional
More informationComputer Simulation of Particle Suspensions
Computer Simulation of Particle Suspensions Jens Harting, Martin Hecht, Hans J. Herrmann, and Sean McNamara Institute for Computational Physics, University of Stuttgart, Pfaffenwaldring 27, D-70569 Stuttgart
More informationSimulation of T-junction using LBM and VOF ENERGY 224 Final Project Yifan Wang,
Simulation of T-junction using LBM and VOF ENERGY 224 Final Project Yifan Wang, yfwang09@stanford.edu 1. Problem setting In this project, we present a benchmark simulation for segmented flows, which contain
More informationTable of Contents. Preface... xiii
Preface... xiii PART I. ELEMENTS IN FLUID MECHANICS... 1 Chapter 1. Local Equations of Fluid Mechanics... 3 1.1. Forces, stress tensor, and pressure... 4 1.2. Navier Stokes equations in Cartesian coordinates...
More informationDrag Law for Bidisperse Gas-Solid Suspensions Containing Equally Sized Spheres
Ind. Eng. Chem. Res. 009, 48, 7 41 7 Drag Law for Bidisperse Gas-Solid Suspensions Containing Equally Sized Spheres Xiaolong Yin and Sankaran Sundaresan* Chemical Engineering Department, Princeton UniVersity,
More informationPerfect entropy functions of the Lattice Boltzmann method
EUROPHYSICS LETTERS 15 July 1999 Europhys. Lett., 47 (2), pp. 182-188 (1999) Perfect entropy functions of the Lattice Boltzmann method I. V. Karlin( ), A. Ferrante and H. C. Öttinger ETH Zürich, Department
More informationarxiv:cond-mat/ v1 8 Nov 2005
Test of the Stokes-Einstein relation in a two-dimensional Yukawa liquid Bin Liu and J. Goree Department of Physics and Astronomy, The University of Iowa, Iowa City, Iowa 52242 arxiv:cond-mat/0511209 v1
More informationEntanglements. M < M e. M > M e. Rouse. Zero-shear viscosity vs. M (note change of slope) Edwards degennes Doi. Berry + Fox, slope 3.4.
Entanglements Zero-shear viscosity vs. M (note change of slope) M < M e Rouse slope 3.4 M > M e Edwards degennes Doi slope 1 Berry + Fox, 1968 Question: Which factors affect the Me: T, P, M, flexibility,
More informationPrinciples of Convection
Principles of Convection Point Conduction & convection are similar both require the presence of a material medium. But convection requires the presence of fluid motion. Heat transfer through the: Solid
More informationSteady Flow and its Instability of Gravitational Granular Flow
Steady Flow and its Instability of Gravitational Granular Flow Namiko Mitarai Department of Chemistry and Physics of Condensed Matter, Graduate School of Science, Kyushu University, Japan. A thesis submitted
More informationLectures on basic plasma physics: Introduction
Lectures on basic plasma physics: Introduction Department of applied physics, Aalto University Compiled: January 13, 2016 Definition of a plasma Layout 1 Definition of a plasma 2 Basic plasma parameters
More informationWay to Success Model Question Paper. Answer Key
A Way to Success Model Question Paper,aw;gpay; / PHYSICS Answer Key (Based on new Question pattern 2019) gphpt I / SECTION I 1 (c) 9 (b) 2 jpirntfk; velocity 10 (d) 3 AB cosθ; 11 (c) 4 v = f λ 12 (b) 5
More informationJ07M.1 - Ball on a Turntable
Part I - Mechanics J07M.1 - Ball on a Turntable J07M.1 - Ball on a Turntable ẑ Ω A spherically symmetric ball of mass m, moment of inertia I about any axis through its center, and radius a, rolls without
More informationChristel Hohenegger A simple model for ketchup-like liquid, its numerical challenges and limitations April 7, 2011
Notes by: Andy Thaler Christel Hohenegger A simple model for ketchup-like liquid, its numerical challenges and limitations April 7, 2011 Many complex fluids are shear-thinning. Such a fluid has a shear
More informationReview of fluid dynamics
Chapter 2 Review of fluid dynamics 2.1 Preliminaries ome basic concepts: A fluid is a substance that deforms continuously under stress. A Material olume is a tagged region that moves with the fluid. Hence
More informationAccelerated Stokesian Dynamics simulations
J. Fluid Mech. (2001), vol. 448, pp. 115 146. c 2001 Cambridge University Press DOI: 10.1017/S0022112001005912 Printed in the United Kingdom 115 Accelerated Stokesian Dynamics simulations By ASIMINA SIEROU
More informationMicrostructural theory and the rheology of. concentrated colloidal suspensions
Under consideration for publication in J. Fluid Mech. 1 Microstructural theory and the rheology of concentrated colloidal suspensions E H S S A N N A Z O C K D A S T 1 and J E F F R E Y F. M O R R I S
More informationCoarse-graining limits in open and wall-bounded dissipative particle dynamics systems
THE JOURNAL OF CHEMICAL PHYSICS 124, 184101 2006 Coarse-graining limits in open and wall-bounded dissipative particle dynamics systems Igor V. Pivkin and George E. Karniadakis a Division of Applied Mathematics,
More informationModeling the combined effect of surface roughness and shear rate on slip flow of simple fluids
Modeling the combined effect of surface roughness and shear rate on slip flow of simple fluids Anoosheh Niavarani and Nikolai Priezjev www.egr.msu.edu/~niavaran November 2009 A. Niavarani and N.V. Priezjev,
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationDSMC-Based Shear-Stress/Velocity-Slip Boundary Condition for Navier-Stokes Couette-Flow Simulations
DSMC-Based Shear-Stress/Velocity-Slip Boundary Condition for Navier-Stokes Couette-Flow Simulations J. R. Torczynski and M. A. Gallis Engineering Sciences Center, Sandia National Laboratories, P. O. Box
More informationParticles in Fluids. Sedimentation Fluidized beds Size segregation under shear Pneumatic transport Filtering Saltation Rheology of suspensions
Particles in Fluids Sedimentation Fluidized beds Size segregation under shear Pneumatic transport Filtering Saltation Rheology of suspensions Sandstorm Fluidized Bed Equation of motion v of v fluid v vx
More informationDynamic Simulation of Shear-induced Particle Migration in a Two-dimensional Circular Couette Device *
Chin. J. Chem. Eng., 15(3) 333 338 (2007) Dynamic Simulation of Shear-induced Particle Migration in a Two-dimensional Circular Couette Device * YU Zhaosheng( 余钊圣 ) a,b, **, SHAO Xueming( 邵雪明 ) a and Roger
More informationConservation of Mass. Computational Fluid Dynamics. The Equations Governing Fluid Motion
http://www.nd.edu/~gtryggva/cfd-course/ http://www.nd.edu/~gtryggva/cfd-course/ Computational Fluid Dynamics Lecture 4 January 30, 2017 The Equations Governing Fluid Motion Grétar Tryggvason Outline Derivation
More informationAtoms, electrons and Solids
Atoms, electrons and Solids Shell model of an atom negative electron orbiting a positive nucleus QM tells that to minimize total energy the electrons fill up shells. Each orbit in a shell has a specific
More informationBAE 820 Physical Principles of Environmental Systems
BAE 820 Physical Principles of Environmental Systems Stokes' law and Reynold number Dr. Zifei Liu The motion of a particle in a fluid environment, such as air or water m dv =F(t) - F dt d - 1 4 2 3 πr3
More informationParticle-particle interactions and models (Discrete Element Method)
Particle-particle interactions and models (Discrete Element Method) Stefan Luding MSM, TS, CTW, UTwente, NL Granular Materials Real: sand, soil, rock, grain, rice, lentils, powder, pills, granulate, micro-
More informationIMPLEMENTING THE LATTICE-BOLTZMANN
IMPLEMENTING THE LATTICE-BOLTZMANN METHOD A RESEARCH ON BOUNDARY CONDITION TECHNIQUES by S.C. Wetstein in partial fulfillment of the requirements for the degree of Bachelor of Science in Applied Physics
More informationCandidacy Exam Department of Physics February 6, 2010 Part I
Candidacy Exam Department of Physics February 6, 2010 Part I Instructions: ˆ The following problems are intended to probe your understanding of basic physical principles. When answering each question,
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationREGULARIZATION AND BOUNDARY CONDITIONS FOR THE 13 MOMENT EQUATIONS
1 REGULARIZATION AND BOUNDARY CONDITIONS FOR THE 13 MOMENT EQUATIONS HENNING STRUCHTRUP ETH Zürich, Department of Materials, Polymer Physics, CH-8093 Zürich, Switzerland (on leave from University of Victoria,
More information