Institute for Condensed Matter Theory and Beckman Institute, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
|
|
- Irene Waters
- 5 years ago
- Views:
Transcription
1 Shear-thinning of molecular fluids in Couette flow Bharath V. Raghavan 1 and Martin Ostoja-Starzewski 2 1) Department of Mechanical Science and Engineering,University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA 2) Department of Mechanical Science and Engineering and Institute for Condensed Matter Theory and Beckman Institute, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA (Dated: 27 January 2017) We use Non-equilibrium Molecular Dynamics (NEMD) simulations, the Boltzmann equation, and continuum thermomechanics to investigate and characterize the shearthinning behavior of molecular fluids undergoing Couette flow, interacting via a Lennard-Jones (LJ) potential. In particular, we study the shear-stress under steadystate conditions and its dependency on fluid density and applied shear-strain rate. Motivated by kinetic theory, we propose a rheological equation of state that fits observed system responses exceptionally well and captures the extreme shear-thinning effect. We notice that beyond a particular strain-rate threshold, the fluid exhibits shear-thinning, the degree of which is dependent on the density and temperature of the system. In addition, we obtain a shear-rate dependent model for the viscosity which matches the well established Cross viscosity model. We demonstrate how this model arises naturally from the Boltzmann equation, and possesses an inherent scaling parameter that unifies the rheological properties of the LJ fluid. We compare our model with those in the literature. Finally, we formulate a dissipation function modeling the LJ fluid as a quasilinear fluid. 1
2 I. INTRODUCTION Current developments in the field of micro and nanoscale systems have led to observations of interesting properties in fluids. At the molecular level, classical (continuum) fluid dynamic models (i.e. the Navier-Stokes equations) often breakdown due to their inability to account for the discrete nature of the fluid 1, and thus its observable properties. In such cases, one often has to resort to Molecular Dynamics (MD) simulations in order to account for the discrete nature of the fluid. Non-equilibrium Molecular Dynamics (NEMD) and in particular the SLLOD algorithm are useful techniques for studying the rheological properties of molecular fluids. Shearthinning has been widely observed in molecular fluids under SLLOD dynamics. Heyes 2 used macro-rheological equations to model the viscosity of a LJ fluid undergoing SLLOD dynamics. Additionally, Heyes 3 also studied the structural features of the flow that result in shear-thinning such as the shear-rigidity modulus and the structural anisotropy of the flow. He investigated the transition from an amorphous liquid to a string phase that induces shear-thinning. Furthermore, it was found that the viscosity was proportional to γ 0.5 in the nonlinear regime of the flow 4,5, and via mode-coupling analysis 6. Recently, experimental verification has been presented for shear-thinning in NEMD simulations where the authors fitted experimental and simulation results for the viscosity of squalane to a single Carreau equation 7. NEMD has seen widespread applicability in determining the rheological properties of alkanes 8 11, polymer melts 12, and colloidal suspensions 13. These studies have extensively investigated the nonlinear shear-thinning regime of the flow, and the structural features at the onset of shear-thinning. A framework of using Monte-Carlo simulations and NEMD simulations using the SLLOD equations to study the non-newtonian shear viscosities of simple liquids under shear has been proposed 14. The results from Monte-Carlo simulations and MD simulations compare excellently. This theory was extended by removing the limitation of requiring isothermal conditions 15, and was subsequently used to extensively study the rheological properties of the LJ fluid. Applying the Generalized Boltzmann equation, a nonlinear expression for the viscosity of the LJ fluid similar to the Ree-Eyring formula, but absent empirical parameters, was derived. It was discovered that the scaled shear-viscosities at various state-points collapsed onto a single curve given by this expression, alluding to the 2
3 existence of a universal scaling function for all temperatures and densities 15. The classical Navier-Stokes (NS) equations for Couette flow consider a linear dependence between shear stress and deformation rate. As a consequence, in studies covering the transition to turbulence, there are instances where transition occurs at a lower Reynolds number compared to that predicted by linear stability analyses 16. This is believed to be due to a drop in internal friction between fluid layers that results in shear-thinning. Conversely, it has been found that the onset of turbulence can be delayed in Couette flow experiments by dissolving a small amount of polymer molecules in water without altering its Newtonian nature 17. Motivated by this need for a more general constitutive stress-deformation model, we resort to kinetic theory and MD to study the steady-state shear stress and its dependency on the density and temperature of the system, and the applied shear strain rate. We propose a rheological equation of state to capture the shear-thinning feature observed in the stress-deformation plot. Finally, we compare our model to other equations of state established in the literature. From this equation, we obtain a model for the viscosity of the LJ fluid that depends entirely on state parameters (i.e. the density and temperature) and the applied strain-rate. II. SIMULATION METHODOLOGY A. Theoretical background The problem of simulating homogeneous flows driven by boundaries in real physical systems (e.g:-couette or elongational flows) is that a microscopic simulation explicitly including the walls invariably induces density inhomogeneities into the fluid 18. This profound shortcoming gave rise to NEMD, which currently is well established as an accurate method for simulating homogeneous shear flows. Two well-known algorithms in this category are the DOLLS and SLLOD algorithms. These produce a microscopic description consistent with macroscopic observations. The first of the homogeneous NEMD algorithms was based on the DOLLS Hamiltonian 19 defined as 3
4 (a) (b) FIG. 1. (a) Sliding brick and (b) deforming cube representations of the Lees-Edwards periodic BCs used in NEMD simulations of molecular fluids undergoing Couette flow. H dolls (r N, p N, t) = ϕ(r N )+ i 2m i p 2 i + i r i. v.p i Θ(t) (1) where ϕ(r N ) is the system potential energy due to the interaction between N atoms, r i and p i are the laboratory position and momentum of atom i, and v is the gradient of the streaming velocity v. Furthermore, we enforce that the flow begins at t = 0 as indicated by Θ(t) which is the Heaviside step function. This Hamiltonian generates the DOLLS equations of motion ṙ i = ṗi m + r i. v 4 (2a)
5 ṗ i = F ϕ i + v.p i (2b) where F ϕ i is the interatomic force on atom i due to all the other atoms. However, to account for linear and nonlinear system responses, a modification was proposed as ṙ i = ṗi m + r i. v ṗ i = F ϕ i + p i. v (3a) (3b) which are the so-called SLLOD equations of motion to indicate the transpose v.p i p i. v. In addition, it should be noted that the SLLOD equations cannot be derived from a system Hamiltonian. To obtain an exact representation of planar Couette flow, we apply the Lees-Edwards periodic boundary conditions 6,23 and inducing a streaming velocity field of the form v(r, t) = i γy v(r, t) = (v x, v y, v z ) = ( γy, 0, 0) (4) where i is the unit vector, y is the channel height, and γ = v x / y is the strain rate such that v(r, t) = γy 0 0 (5) Then, the SLLOD equations are expressed as ṙ i = ṗi m + i γy i ṗ i = F ϕ i i γp yi (6a) (6b) A schematic of the MD simulation methodology is shown in Fig. 1. The two representations are equivalent 18. B. Molecular Dynamics model The molecular dynamics simulations are carried out in two-dimensions in LAMMPS 24, with particles interacting via the Lennard-Jones (LJ) potential, the interaction potential for which is given by 4ϵ[( σ ϕ(r N r ) = )12 ( σ r )6 ], r r c (7) 0, r > r c 5
6 where ϵ is the depth of the potential well, and σ is the finite distance at which the interparticle potential is zero. The dynamics of the simulation is governed by the SLLOD equations [Eq. (6)] with Lees-Edwards periodic boundary conditions. To induce isothermal conditions, a Nose-Hoover thermostat in the canonical ensemble (NVT) is implemented, to yield the modified set of equations as 5 ṙ i = ṗi m + i γy i ṗ i = F ϕ i i γp yi αp i α = 1 [ ] T (t) 1 Q T 0 (8a) (8b) (8c) where α is the thermostatting multiplier, T (t) is the actual kinetic temperature at time t, T 0 is the preset kinetic temperature, and Q is a constant. The equilibrium distribution function for the Nose-Hoover thermostat in the canonical ensemble is given by 6 f c (Γ, ζ) = exp[ β(h Qζ2 )] dγdζ exp[ β(h Qζ2 )] (9) where Q and ζ are thermostat parameters. III. RESULTS AND DISCUSSION A. Shear-thinning model based on kinetic theory We now employ the Boltzmann equation for the Couette flow under investigation. The general expression for the Boltzmann equation is f t + v. f + F f m v = J[f, f] (10) where f is the distribution function, v is the velocity, F is the external force, m is the mass, and J[f, f] is the collision operator 25. For simplicity, we consider the BGK model 26 of the form f t + v. f + F f m v = ν(f f L) (11) where f L is the equilibrium distribution function, and ν is an effective collision parameter. For isothermal, uniform shear flow conditions, we have u i (r) = d ij r j (12a) 6
7 d = d ij = j u i = γδ ix δ jy T = 0 (12b) (12c) In order to analyze SLLOD dynamics under Lees-Edwards periodic BCs using the Boltzmann equation, we use a Galilean transformation to transition to the local rest frame (Lagrangian frame) of the fluid. Consequently, we define the peculiar velocities of the fluid particles as 14,15,25 V = v u(r) and r = r ut. In the local rest frame of the fluid, f(r, v, t) = f(v, t). Using Eqs. (12) in the BGK approximation in Eq. (11), and assuming the external thermostatting force is of the form F = mαv,where α is the thermostatting mulitplier related to the heating rate, we obtain which upon simplification yields f t [ f f t + v. r r r + f V αv f v + [I dt] V f r ] V r V v = ν(f f L) (13) [dv + αv ] f V = ν(f f L) (14) Under the conditions of statistical spatial homogeneity, we can neglect spatial derivatives in Eq. (14) to obtain f t [d ijv j + αv i ] f V i = ν(f f L ) (15) At this point, we note that the pressure tensor is the second moment of the peculiar velocity and is expressed as P := mv V f(v, t)dv (16) where denotes the vector outer product. Multiplying Eq. (15) by V i V j and integrating yields P ij t + d ik P jk + d jk P ik + 2αP ij = ν(p ij pδ ij ) (17) Expanding Eq.(17) and bearing in mind the following; P 13 = P 23 = P 33 = 0 under steady-state conditions p = 1 3 [P 11 + P 22 + P 33 ] = 1 3 trp 7
8 all time derivatives vanish under steady-state conditions we replace P 11 = p to obtain the set of coupled differential equations governing the timeevolution of the pressure tensor as P 22 Neglecting all the time-derivatives yields, p t γp αp = 0 (18a) + 2αP 22 + νp 22 = νp (18b) t P 12 + γp αP 12 = νp 12 (18c) t 2 3 γp αp = 0 (19a) 2αP 22 + νp 22 = νp (19b) γp αP 12 = νp 12 (19c) Eq. (19c) yields In addition, Eq. (19a) yields the thermostatting multiplier as P 12 = γ P 22 ν 1 + 2α (20) ν α = γp 12 3p (21) Subsituting Eq. (21) in Eq. (20), noting the BGK approximation that ν = p η 0, where η 0 is the Newtonian viscosity of the fluid, and that under steady-state conditions P 22 = p we have η 0 γ P 12 = [ ν γp 12 3p Substituting P 12 = η 0 γ in Eq. (22), and noting that P 12 = τ 12 yields τ 12 = ] (22) η 0 γ ν 2 γ 2 (23) where τ 12 is the shear stress (which is the negative of the shear component of the pressure tensor). To find a suitable expression for the effective collision term ν in terms of reduced LJ parameters (i.e. the temperature T, density ρ, molecular mass m, and viscosity η 0, we know that 27 ν = 1 ( ) 1 8kB T 2 λ πm (24) 8
9 where λ is the molecular mean free path. The mean free path for extended potentials is defined as 28,29 λ = η 0 ρ ( ) 1 πm 2 2k B T Substituting Eq. (25) in Eq. (24) and converting to reduced LJ units 15, we have (25) ν = νσ (m/ϵ) 1 2 = 4ρ T πη 0m (26) Thus, in reduced LJ units, the shear-thinning equation of state is τ 12 = τ 12 ( σ 3 /ϵ ) = η 0 γ ( γ /ν ) 2 (27) B. Comparison of Numerical and Analytical results The data of primary relevance in this study is the mean value of the steady-state shear stress (i.e. the negative of the shear component of the pressure tensor, the latter being the default output in LAMMPS). After the initial equilibration period, we extract mean values for time-steps at t = increments. We divided the total simulation into blocks of M = 1, 000 time-steps which yields a total number of N = 24, 500 blocks. For each block, we calculate the section average which uses the 1, 000 data points associated with that block given by P 12 ( γ) B = M i=0 P (i) 12 (28) and the cumulative average which accounts for all the previous section averages expressed as P 12 ( γ) C = N i=0 P 12 ( γ) (i) B (29) The steady-state is achieved when fluctuations in the cumulative average are no longer statistically significant. We run simulations for various state-points (Table I) and compare the proposed shear-thinning model, and the associated viscosity model, with those available in the literature. The simulations are carried out for a system size of N = 200 particles, and strain rates 0.0 γ = γσ m/ϵ 2.5 with γ = 0.5 (in LJ reduced units), and a cutoff radius of r c = 2 1/6 σ. The system is run for a total of time steps with t = ; it is sufficiently equilibrated for time-steps and then run for time-steps for data extraction. 9
10 TABLE I. LJ state points and parameter estimates in Eq. (27) obtained via curve fitting. ρ T η0 2/(3ν 2 ) Consistent with well-established observations made in the literature, the LJ fluid exhibits shear-thinning, the degree of which depends on the fluid state point as illustrated in Fig. 2. At low densities, the fluid only experiences moderate shear-thinning (Fig. 2a) for the range of strain-rates considered. With increasing density, the fluid exhibits extreme shear-thinning as shown by its stress-deformation response in Fig. 2b. The level of shear-thinning tends to set in quicker as the fluid density increases. For instance, extreme shear-thinning occurs at γ 1.0 at the triple point (i.e ρ = and T = 0.722), but only at γ 1.5 for ρ = 0.8 and T = 1.0. Furthermore, an increase in temperature tends to delay the onset of extreme shear-thinning. Comparing two state points at ρ = 0.8 and T = 1.0 and ρ = 0.8 and T = 1.1, the latter requires a strain rate of γ 2.0 compared to γ 1.5 for the former. In either case, we see that Eq. (27) is able to capture he stress-deformation response extremely well. In order to evaluate the effectiveness of our proposed model, we compare it with widely used empirical formulae available in the literature and its goodness-of-fit measured by the r 2 value. The models of interest are; The Sisko equation τ 12 = a γ b + c γ (30) The Ree-Eyring equation τ 12 = k BT V τ sinh 1 ( τ γ) (31) where V τ is the shear stress activation volume and τ is the relaxation time. 10
11 3 2.5 (0.4,1.5) (0.6,1.0) (0.6,1.5) 2 τ γ (a) τ (0.8442,0.722) (0.8,1.0) (0.8,1.1) (0.9,1.1) γ (b) FIG. 2. (a) Moderate and (b) extreme shear-thinning in the stress-deformation curve for the LJ fluid at various state points. The curve-fit is performed using Eq. (27) and is represented by the dotted lines. The Bair-Winer equation τ 12 = τ L 12 ( 1 e η 0 γ τ L 12 ) (32) where τ L 12 is the limiting shear stress; a concept frequently used in high-pressure lubri- 11
12 cation. In all the cases, the rheological equations of state from the literature studied in this article provide excellent models for cases of moderate shear-thinning, as illustrated in Fig. 3, as observed at state points of (ρ, T ) = (0.4, 1.5), (0.6, 1.0) and (0.6, 1.5). However, these models fail to accurately capture the dramatic shear-thinning feature of the LJ fluid at higher densities of (ρ, T ) = (0.8, 1.0), (0.8, 1.1), (0.8442, 0.722) and (0.9, 1.1). We observe that for moderate shear-thinning, the Sisko equation of state, Eq. (30), is comparable in effectiveness to Eq. (27) and fits the data with r However, the former fails to fit the data at state points where dramatic shear-thinning occurs. On the other hand, Eq. (27) provides an exceptionally good fit for all the LJ state points considered, and in particular is able to capture the extreme shear-thinning effect. The Ree-Eyring [Eq. (31)] and Bair-Winer [Eq. (32)] equations of state, despite being widely considered more versatile than the former equations, still suffer from the same drawbacks as being unable to capture dramatic shearthinning. In general, these equations are excellent models at low strain-rates, i.e. below the threshold where extreme shear-thinning sets in. The models are able to successfully capture the stress-deformation curve with r for moderate shear-thinning, which falls to around 0.7 r for extreme shear-thinning stress-deformation behavior. The poor performance of these models at higher strain rates can be attributed to their empirical nature,where parameters are obtained by fitting experimental data, or their derivation from macro-viscoelastic models. In order to obtain a strain-dependent viscosity model, we note that τ 12 = η( γ) γ (33) Thus, using Eq. (27) in Eq. (33) we find that η 1 = η0 ( γ ) 2 (34) 3 ν We compare our model with the well-established Cross model given in the literature as η = η 0 η 1 + β γ n + η (35) where β and n are fitting constants, and η is the infinite shear-rate limiting viscosity. Interestingly, Eq. (34) yields the Cross viscosity model, with n = 2 and η = 0, absent 12
13 3.5 3 Eq. (30) Eq. (31) Eq. (32) 2.5 τ γ FIG. 3. Demonstration of the ability of existing rheological equations [Eq. (30)-(32)] to capture moderate shear-thinning in the LJ fluid and their inability to model extreme shear-thinning. The data markers are the same as used to identify the state points in Fig. 2 empirical parameters since the parameters in Eq. (34) are completely determined by the fluid s density (ρ ), temperature (T ), and Newtonian viscosity (η0). To serve as a validation for our simulations, Eq. (34) yields a viscosity of η 0 = 3.76 at the LJ triple point of (ρ, T ) = (0.8442, 0.722), which is comparable to the value of η 0 = 3.73 calculated in the literature 30. The Cross model has been previously used to model the strain-rate dependent viscosities of simple fluids, but in an empirical sense 3,31. It should be noted that our model has been completely derived from hydrodynamic and kinetic theories. A similar result has been obtained for the Ree-Eyring viscosity formula, which arose naturally from the Generalized Boltzmann equation. In our model, we find that the effective collision parameter, ν, behaves as a universal scaling parameter that facilitates collapsing all the data onto a universal curve, suggesting the existence of a universal scaling law for the LJ fluid, as shown in Fig. 4. The slight scatter is attributed to the inaccuracy of NEMD in accessing small strain-rates for viscosity calculations. Ideally, one would resort to the use of Transient Time Correlation functions 32. Additionally, we attribute the discrepancy between our calculations and those avaliable in the literature to the truncation distance. In most studies, a truncation of r c = 2.5σ was used which would account for the attractive component in the LJ potential as well, compared to r c = 2 1/6 σ in this work. 13
14 10 1 (0.6,1.0) (0.8442,0.722) (0.4,1.5) (0.6,1.5) (0.8,1.0) (0.8,1.1) (0.9,1.1) η /η γ /ν FIG. 4. Scaled viscosities for the LJ fluid at various state points for a system size of N = 200 particles. The fit is given by Eq. (34). C. Quasilinearity based on continuum thermomechanics Note that Eq. (27) may allow one to determine a dissipation function which should equal the irreversible rate of entropy production 33. Within the framework of thermomechanics, we know that the dissipation function is defined as A (d) kl (d ij)d ij = Φ(d ij ) (36) Where the dissipative forces A (d) kl define a symmetric second-order tensor and d ij is the deformation rate. For uniform shear flow, the dissipative forces are the shear stress σ (d) ij that σ (d) ij (d ij)d ij = Φ(d ij ) = so η 0 γ 2 ( γ ) 2 (37) 3 ν Adopting thermodynamic orthogonality (or principle of maximal rate of entropy production ), we have σ (d) ij = µ Φ d ij (38) where µ is a constant of proportionality to be determined by the condition of thermodynamic orthogonality. Under uniform shear, the invariants of the strain-rate tensor are d (1) = 0, 14
15 d (2) = γ 2, and d (3) = 0. Therefore, on account of these side conditions, Eq. (38) becomes σ (d) ij = µ [ ] Φ λ1 d (1) λ 2 d (3) (39) d ij where λ 1 and λ 2 are lagrange multipliers. Since d ij is a deviatoric tensor, then so is σ (d) ij. Thus, we have tr[σ (d) ij (39), we obtain which simplifies to On account of the constraints, we have ] = 0. Furthermore, we also have det [σ(d)] = 0. Upon expanding Eq. [ σ (d) Φ d (1) ij =µ + Φ d (2) d (1) d ij d (2) d ij + Φ d (3) d (3) d ij λ 1 d (1) d ij λ 2 d (3) d ij [ σ (d) Φ ij =µ (d ij d (1) δ ij ) d (12) ] λ 1 δ ij λ 2 (d ik d kj d 2 δ ij ) ij ] (40) (41) 3λ 1 λ 2 (d ik d kj 3d 2 ) = 0 3λ 1 + λ 2 d 2 = 0 λ 1 + λ 2 d (2) = 0 (42a) (42b) Solving the system of equations, Eq. (42a) and (42b), we obtain λ 1 = λ 2 = 0, so that σ (d) ij Multiplying both sides of Eq. (43) by d ij yields For uniform shear, we obtain = µ Φ d (2) d ij (43) Φ = µ Φ d (2) d ij d ij (44) Φ = 2µ Φ ( γ 2 ) γ2 (45) Substituting the explicit form for Φ in Eq. (45) and taking the derivative with respect to γ 2, we obtain µ = 1 + (2/3) ( γ2 /ν 2 ) 2 (46) Thus we see that a dissipation function of the form in Eq. (37) along with Eq. (43) suggest that a molecular fluid behaves at least like a quasi-linear fluid 33 with a constitutive relationship of the form σ d ij = 2η(d (2) )d ij σ d 12 = 2η( γ 2 ) γ (47) where the fluid viscosity η needs to atleast depend on d (2) = γ 2. 15
16 IV. CONCLUSION In this article, we simulate LJ fluids at different state points and employ the Boltzmann equation to study the behavior of the steady-state shear stress with applied strain rate. We observe that, in all cases, the fluid exhibits shear-thinning. At low densities, the shearthinning effect is moderate for the range of strain-rates considered. At higher densities, dramatic shear-thinning sets in at much lower strain-rates. In light of these observations, we are motivated to use kinetic theory to find an equation of state that is able to capture this wide range of behaviors. We compare our proposed equation with other models available in the literature to validate its effectiveness. In general, we find that all the models available in the literature performed well when the degree of shear- thinning is moderate. However, for cases of extreme shear-thinning, the Ostwalde-de Waele and Sisko equations of state completely breakdown, while the Ree- Eyring and Bair-Winer equations experience a significant drop in modeling capability. In contrast, the proposed model performs consistently well across all state-points. Additionally, the formalism involving the Boltzmann equation naturally gives rise to a form of the Cross viscosity model as a steady-state solution. Moreover, the effective collision frequency used in the BGK approximation acts as a scaling parameter that seems to unify the behavior of the LJ fluid across all the state-points considered. Finally, a thermomechanical framework involving the dissipation function suggests that the LJ fluid behaves at least as a quasilinear fluid. This offers an insight into the minimum required approximations to be used in the process of constructing a suitable constitutive model to represent a molecular fluid. The derivation shows that the proposed equation provides a suitable shear-thinning model while preserving thermodynamic orthogonality in molecular fluids, which is an important feature of continuum thermomechanics. The framework of thermodynamic orthogonality provides a means for transitioning from molecular to continuum scales and for a link to the fluctuation theorem which grasps spontaneous violations of the second law of thermodynamics on very small length and/or time scales 34. Thus, for the purposes of modeling such fluids as continua, we need to assume that the viscosity is dependent on the strain-rate in a manner consistent with the postulates of continuum thermomechanics. 16
17 ACKNOWLEDGMENTS This work has been supported by the NSF under grant CMMI The authors would like to thank the University of Illinois for computational resources necessary for undertaking this project. Additionally, the authors would like to thank Dr. R. E. Khayat of the University of Western Ontario and three anonymous reviewers for constructive comments on this manuscript. REFERENCES 1 J. Delhommelle and D. J. Evans, Poiseuille flow of a micropolar fluid, Molecular Physics 100, (2002). 2 D. M. Heyes, Shear thinning of the Lennard-Jones fluid by molecular dynamics, Physica A: Statistical Mechanics and its Applications 133, (1985). 3 D. M. Heyes, Non-Newtonian behaviour of simple liquids, Journal of non-newtonian fluid mechanics 21, (1986). 4 D. J. Evans, Molecular dynamics simulations of the rheological properties of simple fluids, Physica A: Statistical Mechanics and its Applications 118, (1983). 5 A. Baranyai and P. T. Cummings, Fluctuations close to equilibrium, Physical Review E 52, 2198 (1995). 6 G. P. Morriss and D. J. Evans, Statistical Mechanics of Nonequilbrium Liquids (ANU Press, 2007). 7 S. Bair, C. McCabe, and P. T. Cummings, Comparison of nonequilibrium molecular dynamics with experimental measurements in the nonlinear shear-thinning regime, Physical review letters 88, (2002). 8 S. Cui, S. Gupta, P. Cummings, and H. Cochran, Molecular dynamics simulations of the rheology of normal decane, hexadecane, and tetracosane, The Journal of chemical physics 105, (1996). 9 S. Cui, P. Cummings, H. Cochran, J. Moore, and S. Gupta, Nonequilibrium molecular dynamics simulation of the rheology of linear and branched alkanes, International journal of thermophysics 19, (1998). 10 A. Berker, S. Chynoweth, U. C. Klomp, and Y. Michopoulos, Non-equilibrium molecular 17
18 dynamics (NEMD) simulations and the rheological properties of liquid n-hexadecane, Journal of the Chemical Society, Faraday Transactions 88, (1992). 11 H.-C. Tseng, J.-S. Wu, and R.-Y. Chang, Shear thinning and shear dilatancy of liquid n-hexadecane via equilibrium and nonequilibrium molecular dynamics simulations: Temperature, pressure, and density effects, The Journal of chemical physics 129, (2008). 12 Y.-R. Jeng, C.-C. Chen, and S.-H. Shyu, A molecular dynamics study of lubrication rheology of polymer fluids, Tribology Letters 15, (2003). 13 X. Xu, S. A. Rice, and A. R. Dinner, Relation between ordering and shear thinning in colloidal suspensions, Proceedings of the National Academy of Sciences 110, (2013). 14 H. Farhat and B. C. Eu, Monte Carlo method and the nonequilibrium structure and non-newtonian viscosity of a sheared simple liquid, The Journal of chemical physics 110, (1999). 15 R. Laghaei, A. E. Nasrabad, and B. C. Eu, Statistical-mechanical theory of rheology: Lennard-Jones fluids, The Journal of chemical physics 123, (2005). 16 K. Volokh, An investigation into the stability of a shear thinning fluid, International Journal of Engineering Science 47, (2009). 17 K. Volokh, Navier-Stokes model with viscous strength, Computer Modeling in Engineering & Sciences 92, (2013). 18 B. Todd and P. J. Daivis, Homogeneous non-equilibrium molecular dynamics simulations of viscous flow: techniques and applications, Molecular Simulation 33, (2007). 19 W. G. Hoover, D. J. Evans, R. B. Hickman, A. J. Ladd, W. T. Ashurst, and B. Moran, Lennard-Jones triple-point bulk and shear viscosities. Green-Kubo theory, Hamiltonian mechanics, and nonequilibrium molecular dynamics, Physical Review A 22, 1690 (1980). 20 D. J. Evans and G. P. Morriss, Non-Newtonian molecular dynamics, Computer Physics Reports 1, (1984). 21 D. J. Evans and G. Morriss, Nonlinear-response theory for steady planar Couette flow, Physical Review A 30, 1528 (1984). 22 A. J. Ladd, Equations of motion for non-equilibrium molecular dynamics simulations of viscous flow in molecular fluids, Molecular Physics 53, (1984). 23 A. Lees and S. Edwards, The computer study of transport processes under extreme con- 18
19 ditions, Journal of Physics C: Solid State Physics 5, 1921 (1972). 24 S. Plimpton, Fast parallel algorithms for short-range molecular dynamics, Journal of computational physics 117, 1 19 (1995). 25 J. W. Dufty, A. Santos, J. J. Brey, and R. F. Rodriguez, Model for nonequilibrium computer simulation methods, Physical Review A 33, 459 (1986). 26 R. Zwanzig, Nonlinear shear viscosity of a gas, The Journal of Chemical Physics 71, (1979). 27 W. G. Vincenti and C. H. Kruger, Introduction to physical gas dynamics, (Wiley, New York, 1965). 28 D. L. Morris, L. Hannon, and A. L. Garcia, Slip length in a dilute gas, Physical review A 46, 5279 (1992). 29 S. K. Prabha, P. Sreehari, S. P. Sathian, et al., The effect of system boundaries on the mean free path for confined gases, AIP Advances 3, (2013). 30 W. Ashurst and W. Hoover, Argon shear viscosity via a Lennard-Jones potential with equilibrium and nonequilibrium molecular dynamics, Physical Review Letters 31, 206 (1973). 31 I. Borzsak, P. T. Cummings, and D. J. Evans, Shear viscosity of a simple fluid over a wide range of strain rates, Molecular Physics 100, (2002). 32 G. P. Morriss and D. J. Evans, Application of transient correlation functions to shear flow far from equilibrium, Physical Review A 35, 792 (1987). 33 H. Ziegler and C. Wehrli, The derivation of constitutive relations from the free energy and the dissipation function, Advances in applied mechanics 25, (1987). 34 D. J. Evans, E. Cohen, and G. Morriss, Probability of second law violations in shearing steady states, Physical Review Letters 71, 2401 (1993). 19
20
21
22 3 2.5 (0.4,1.5) (0.6,1.0) (0.6,1.5) 2 τ γ
23 (0.8442,0.722) (0.8,1.0) (0.8,1.1) (0.9,1.1) τ γ
24 3.5 3 Eq. (30) Eq. (31) Eq. (32) 2.5 τ γ
25 10 1 (0.6,1.0) (0.8442,0.722) (0.4,1.5) (0.6,1.5) (0.8,1.0) (0.8,1.1) (0.9,1.1) η /η γ /ν
Statistical Mechanics of Nonequilibrium Liquids
1. Introduction Mechanics provides a complete microscopic description of the state of a system. When the equations of motion are combined with initial conditions and boundary conditions, the subsequent
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 informationBehaviour of microscale gas flows based on a power-law free path distribution function
Behaviour of microscale gas flows based on a power-law free path distribution function Nishanth Dongari, Yonghao Zhang and Jason M Reese Department of Mechanical Engineering, University of Strathclyde,
More informationNonequilibrium Molecular Dynamics Simulations of Shear Viscosity: Isoamyl Alcohol, n-butyl Acetate and Their Mixtures
Draft Date: October 5, 1999 Nonequilibrium Molecular Dynamics Simulations of Shear Viscosity: Isoamyl Alcohol, n-butyl Acetate and Their Mixtures Y. Yang 1, T. A. Pakkanen, 2 and R. L. Rowley 1,3 Draft
More informationNon-equilibrium phenomena and fluctuation relations
Non-equilibrium phenomena and fluctuation relations Lamberto Rondoni Politecnico di Torino Beijing 16 March 2012 http://www.rarenoise.lnl.infn.it/ Outline 1 Background: Local Thermodyamic Equilibrium 2
More informationA Nonequilibrium Molecular Dynamics Study of. the Rheology of Alkanes. S.A. Gupta, S. T. Cui, P. T. Cummings and H. D. Cochran
A Nonequilibrium Molecular Dynamics Study of the Rheology of Alkanes S.A. Gupta, S. T. Cui, P. T. Cummings and H. D. Cochran Department of Chemical Engineering University of Tennessee Knoxville, TN 37996-2200
More informationOleg A. Mazyar, Guoai Pan and Clare McCabe*
Molecular Physics Vol. 107, No. 14, 20 July 2009, 1423 1429 RESEARCH ARTICLE Transient time correlation function calculation of the viscosity of a molecular fluid at low shear rates: a comparison of stress
More informationHydrodynamics. Stefan Flörchinger (Heidelberg) Heidelberg, 3 May 2010
Hydrodynamics Stefan Flörchinger (Heidelberg) Heidelberg, 3 May 2010 What is Hydrodynamics? Describes the evolution of physical systems (classical or quantum particles, fluids or fields) close to thermal
More informationTwo recent works on molecular systems out of equilibrium
Two recent works on molecular systems out of equilibrium Frédéric Legoll ENPC and INRIA joint work with M. Dobson, T. Lelièvre, G. Stoltz (ENPC and INRIA), A. Iacobucci and S. Olla (Dauphine). CECAM workshop:
More informationRheology and Constitutive Equations. Rheology = Greek verb to flow. Rheology is the study of the flow and deformation of materials.
Rheology and Constitutive Equations Rheology = Greek verb to flow Rheology is the study of the flow and deformation of materials. The focus of rheology is primarily on the study of fundamental, or constitutive,
More informationIntroduction to Nonequilibrium Molecular Dynamics
Introduction to Nonequilibrium Molecular Dynamics Dr Billy D. Todd Centre for Molecular Simulation Swinburne University of Technology PO Box 218, Hawthorn Victoria 3122, Australia Email: btodd@swin.edu.au
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 informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationLifshitz Hydrodynamics
Lifshitz Hydrodynamics Yaron Oz (Tel-Aviv University) With Carlos Hoyos and Bom Soo Kim, arxiv:1304.7481 Outline Introduction and Summary Lifshitz Hydrodynamics Strange Metals Open Problems Strange Metals
More informationMath 575-Lecture Viscous Newtonian fluid and the Navier-Stokes equations
Math 575-Lecture 13 In 1845, tokes extended Newton s original idea to find a constitutive law which relates the Cauchy stress tensor to the velocity gradient, and then derived a system of equations. The
More informationRevisit to Grad s Closure and Development of Physically Motivated Closure for Phenomenological High-Order Moment Model
Revisit to Grad s Closure and Development of Physically Motivated Closure for Phenomenological High-Order Moment Model R. S. Myong a and S. P. Nagdewe a a Dept. of Mechanical and Aerospace Engineering
More informationChapter 1. Continuum mechanics review. 1.1 Definitions and nomenclature
Chapter 1 Continuum mechanics review We will assume some familiarity with continuum mechanics as discussed in the context of an introductory geodynamics course; a good reference for such problems is Turcotte
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 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 informationShear viscosity of molten sodium chloride
Shear viscosity of molten sodium chloride Jerome Delhommelle and Janka Petravic Citation: The Journal of Chemical Physics 118, 2783 (2003); doi: 10.1063/1.1535213 View online: http://dx.doi.org/10.1063/1.1535213
More informationEfficient viscosity estimation from molecular dynamics simulation via momentum impulse relaxation
JOURNAL OF CHEMICAL PHYSICS VOLUME 113, NUMBER 6 8 AUGUST 2000 Efficient viscosity estimation from molecular dynamics simulation via momentum impulse relaxation Gaurav Arya, Edward J. Maginn, a) and Hsueh-Chia
More informationFluid Mechanics Theory I
Fluid Mechanics Theory I Last Class: 1. Introduction 2. MicroTAS or Lab on a Chip 3. Microfluidics Length Scale 4. Fundamentals 5. Different Aspects of Microfluidcs Today s Contents: 1. Introduction to
More informationMolecular Dynamics Simulation Study of Transport Properties of Diatomic Gases
MD Simulation of Diatomic Gases Bull. Korean Chem. Soc. 14, Vol. 35, No. 1 357 http://dx.doi.org/1.51/bkcs.14.35.1.357 Molecular Dynamics Simulation Study of Transport Properties of Diatomic Gases Song
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 informationCE 530 Molecular Simulation
1 CE 530 Molecular Simulation Lecture 1 David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu 2 Time/s Multi-Scale Modeling Based on SDSC Blue Horizon (SP3) 1.728 Tflops
More informationNosé-Hoover chain method for nonequilibrium molecular dynamics simulation
PHYSICAL REVIEW E VOLUME 61, NUMBER 5 MAY 000 Nosé-Hoover chain method for nonequilibrium molecular dynamics simulation A. C. Brańka Institute of Molecular Physics, Polish Academy of Sciences, Smoluchowskiego
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 informationLes Houches School of Foam: Rheology of Complex Fluids
Les Houches School of Foam: Rheology of Complex Fluids Andrew Belmonte The W. G. Pritchard Laboratories Department of Mathematics, Penn State University 1 Fluid Dynamics (tossing a coin) Les Houches Winter
More informationChapter 2: Fluid Dynamics Review
7 Chapter 2: Fluid Dynamics Review This chapter serves as a short review of basic fluid mechanics. We derive the relevant transport equations (or conservation equations), state Newton s viscosity law leading
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 informationContinuum Mechanics Fundamentals
Continuum Mechanics Fundamentals James R. Rice, notes for ES 220, 12 November 2009; corrections 9 December 2009 Conserved Quantities Let a conseved quantity have amount F per unit volume. Examples are
More informationViscosity of confined inhomogeneous nonequilibrium fluids
JOURNAL OF CHEMICAL PHYSICS VOLUME 121, NUMBER 21 1 DECEMBER 2004 Viscosity of confined inhomogeneous nonequilibrium fluids Junfang Zhang and B. D. Todd a) Centre for Molecular Simulation, Swinburne University
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 informationEntropy generation and transport
Chapter 7 Entropy generation and transport 7.1 Convective form of the Gibbs equation In this chapter we will address two questions. 1) How is Gibbs equation related to the energy conservation equation?
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 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 informationPredicting the viscosity of alkanes using nonequilibrium molecular dynamics: Evaluation of intermolecular potential models
Brigham Young University BYU ScholarsArchive All Faculty Publications 1997-06-22 Predicting the viscosity of alkanes using nonequilibrium molecular dynamics: Evaluation of intermolecular potential models
More informationMD Thermodynamics. Lecture 12 3/26/18. Harvard SEAS AP 275 Atomistic Modeling of Materials Boris Kozinsky
MD Thermodynamics Lecture 1 3/6/18 1 Molecular dynamics The force depends on positions only (not velocities) Total energy is conserved (micro canonical evolution) Newton s equations of motion (second order
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 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 informationNON-EQUILIBRIUM MOLECULAR DYNAMICS CALCULATION OF THE
FluidPhaseEquilibria, 53 (1989) 191-198 Elsevier Science Publishers B.V., Amsterdam - Printed in The Netherlands 191 NON-EQUILIBRIUM MOLECULAR DYNAMICS CALCULATION OF THE TRANSPORT PROPERTIES OF CARBON
More informationPoiseuille flow of Lennard-Jones fluids in narrow slit pores
JOURNAL OF CHEMICAL PHYSICS VOLUME 112, NUMBER 4 22 JANUARY 2000 Poiseuille flow of Lennard-Jones fluids in narrow slit pores Karl P. Travis a) and Keith E. Gubbins Department of Chemical Engineering,
More informationAssessment of phenomenological models for viscosity of liquids based on nonequilibrium atomistic simulations of copper
THE JOURNAL OF CHEMICAL PHYSICS 123, 104506 2005 Assessment of phenomenological models for viscosity of liquids based on nonequilibrium atomistic simulations of copper Peng Xu, Tahir Cagin, a and William
More informationMadrid, 8-9 julio 2013
VI CURSO DE INTRODUCCION A LA REOLOGÍA Madrid, 8-9 julio 2013 NON-LINEAR VISCOELASTICITY Prof. Dr. Críspulo Gallegos Dpto. Ingeniería Química. Universidad de Huelva & Institute of Non-Newtonian Fluid Mechanics
More informationIn the name of Allah the most beneficent the most merciful
In the name of Allah the most beneficent the most merciful Transient flows of Maxwell fluid with slip conditions Authors Dr. T. Hayat & Sahrish Zaib Introduction Non-Newtonian Newtonian fluid Newtonian
More informationONSAGER S VARIATIONAL PRINCIPLE AND ITS APPLICATIONS. Abstract
ONSAGER S VARIAIONAL PRINCIPLE AND IS APPLICAIONS iezheng Qian Department of Mathematics, Hong Kong University of Science and echnology, Clear Water Bay, Kowloon, Hong Kong (Dated: April 30, 2016 Abstract
More informationModelling the Rheology of Semi-Concentrated Polymeric Composites
THALES Project No 1188 Modelling the Rheology of Semi-Concentrated Polymeric Composites Research Team Evan Mitsoulis (PI), Professor, NTUA, Greece Costas Papoulias (Research Student), NTUA, Greece Souzanna
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 informationAnisotropic fluid dynamics. Thomas Schaefer, North Carolina State University
Anisotropic fluid dynamics Thomas Schaefer, North Carolina State University Outline We wish to extract the properties of nearly perfect (low viscosity) fluids from experiments with trapped gases, colliding
More informationColloidal Suspension Rheology Chapter 1 Study Questions
Colloidal Suspension Rheology Chapter 1 Study Questions 1. What forces act on a single colloidal particle suspended in a flowing fluid? Discuss the dependence of these forces on particle radius. 2. What
More informationModeling the Rheology and Orientation Distribution of Short Glass Fibers Suspended in Polymeric Fluids: Simple Shear Flow
Modeling the Rheology and Orientation Distribution of Short Glass Fibers Suspended in Polymeric Fluids: Simple Shear Flow Aaron P.R. berle, Donald G. Baird, and Peter Wapperom* Departments of Chemical
More informationFluid Dynamics Exercises and questions for the course
Fluid Dynamics Exercises and questions for the course January 15, 2014 A two dimensional flow field characterised by the following velocity components in polar coordinates is called a free vortex: u r
More informationScientific Computing II
Scientific Computing II Molecular Dynamics Simulation Michael Bader SCCS Summer Term 2015 Molecular Dynamics Simulation, Summer Term 2015 1 Continuum Mechanics for Fluid Mechanics? Molecular Dynamics the
More informationAdvances in nonequilibrium molecular dynamics simulations of lubricants and additives
Friction ISSN 2223-7690 https://doi.org/10.1007/s40544-018-0207-9 CN 10-1237/TH REVIEW ARTICLE Advances in nonequilibrium molecular dynamics simulations of lubricants and additives J. P. EWEN, D. M. HEYES,
More informationMicrostructural Randomness and Scaling in Mechanics of Materials. Martin Ostoja-Starzewski. University of Illinois at Urbana-Champaign
Microstructural Randomness and Scaling in Mechanics of Materials Martin Ostoja-Starzewski University of Illinois at Urbana-Champaign Contents Preface ix 1. Randomness versus determinism ix 2. Randomness
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 information(2.1) Is often expressed using a dimensionless drag coefficient:
1. Introduction Multiphase materials occur in many fields of natural and engineering science, industry, and daily life. Biological materials such as blood or cell suspensions, pharmaceutical or food products,
More informationRHEOLOGY Principles, Measurements, and Applications. Christopher W. Macosko
RHEOLOGY Principles, Measurements, and Applications I -56081-5'79~5 1994 VCH Publishers. Inc. New York Part I. CONSTITUTIVE RELATIONS 1 1 l Elastic Solid 5 1.1 Introduction 5 1.2 The Stress Tensor 8 1.2.1
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 informationWall Effects in Convective Heat Transfer from a Sphere to Power Law Fluids in Tubes
Excerpt from the Proceedings of the COMSOL Conference 9 Boston Wall Effects in Convective Heat Transfer from a Sphere to Power Law Fluids in Tubes Daoyun Song *1, Rakesh K. Gupta 1 and Rajendra P. Chhabra
More informationMechanical properties of polymers: an overview. Suryasarathi Bose Dept. of Materials Engineering, IISc, Bangalore
Mechanical properties of polymers: an overview Suryasarathi Bose Dept. of Materials Engineering, IISc, Bangalore UGC-NRCM Summer School on Mechanical Property Characterization- June 2012 Overview of polymer
More informationAGITATION AND AERATION
AGITATION AND AERATION Although in many aerobic cultures, gas sparging provides the method for both mixing and aeration - it is important that these two aspects of fermenter design be considered separately.
More informationKinetic Theory for Binary Granular Mixtures at Low Density
10 Kinetic Theory for Binary Granular Mixtures at Low Density V. Garzó Departamento de Física, Universidad de Extremadura, E-06071 Badajoz, Spain vicenteg@unex.es Many features of granular media can be
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 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 informationElectrical Transport in Nanoscale Systems
Electrical Transport in Nanoscale Systems Description This book provides an in-depth description of transport phenomena relevant to systems of nanoscale dimensions. The different viewpoints and theoretical
More informationFluid Mechanics II Viscosity and shear stresses
Fluid Mechanics II Viscosity and shear stresses Shear stresses in a Newtonian fluid A fluid at rest can not resist shearing forces. Under the action of such forces it deforms continuously, however small
More informationComputation of Shear Viscosity: A Systems Approach
Computation of Shear Viscosity: A Systems Approach Joshua L. Hurst and John T. Wen Dept. of Mechanical, Aerospace, & Nuclear Eng. Dept. of Electrical, Computer, & Systems Eng. Rensselaer Polytechnic Institute
More information1. Introduction, tensors, kinematics
1. Introduction, tensors, kinematics Content: Introduction to fluids, Cartesian tensors, vector algebra using tensor notation, operators in tensor form, Eulerian and Lagrangian description of scalar and
More informationSummary We develop an unconditionally stable explicit particle CFD scheme: Boltzmann Particle Hydrodynamics (BPH)
Summary Steps in BPH Space is divided into Cartesian cells Finite number of particles ~10 6-10 7 Particles fly freely between t n and t n+1 Mass, momentum and total energy are conserved Relax into a (LTE)
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 informationEuler equation and Navier-Stokes equation
Euler equation and Navier-Stokes equation WeiHan Hsiao a a Department of Physics, The University of Chicago E-mail: weihanhsiao@uchicago.edu ABSTRACT: This is the note prepared for the Kadanoff center
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationAdvantages of a Finite Extensible Nonlinear Elastic Potential in Lattice Boltzmann Simulations
The Hilltop Review Volume 7 Issue 1 Winter 2014 Article 10 December 2014 Advantages of a Finite Extensible Nonlinear Elastic Potential in Lattice Boltzmann Simulations Tai-Hsien Wu Western Michigan University
More informationMolecular Dynamics. Timothy O Sullivan Supervised by Naida Lacevic & John Sader The University of Melbourne
Molecular Dynamics Timothy O Sullivan Supervised by Naida Lacevic & John Sader The University of Melbourne 1 1 Introduction Molecular Dynamics (MD) is a method of simulating systems of nanoscale volumes,
More informationG : Statistical Mechanics
G25.2651: Statistical Mechanics Notes for Lecture 15 Consider Hamilton s equations in the form I. CLASSICAL LINEAR RESPONSE THEORY q i = H p i ṗ i = H q i We noted early in the course that an ensemble
More informationShear Thinning Near the Rough Boundary in a Viscoelastic Flow
Advanced Studies in Theoretical Physics Vol. 10, 2016, no. 8, 351-359 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/astp.2016.6624 Shear Thinning Near the Rough Boundary in a Viscoelastic Flow
More informationStructural changes and viscoplastic behavior of a generic embedded-atom model metal in steady shear flow
Structural changes and viscoplastic behavior of a generic embedded-atom model metal in steady shear flow Igor Stankovic* and Siegfried Hess Institut für Theoretische Physik, PN 7-1, Fakultät II, Technische
More informationThe Truth about diffusion (in liquids)
The Truth about diffusion (in liquids) Aleksandar Donev Courant Institute, New York University & Eric Vanden-Eijnden, Courant In honor of Berni Julian Alder LLNL, August 20th 2015 A. Donev (CIMS) Diffusion
More informationRheology of Soft Materials. Rheology
Τ Thomas G. Mason Department of Chemistry and Biochemistry Department of Physics and Astronomy California NanoSystems Institute Τ γ 26 by Thomas G. Mason All rights reserved. γ (t) τ (t) γ τ Δt 2π t γ
More informationTurbulentlike Quantitative Analysis on Energy Dissipation in Vibrated Granular Media
Copyright 011 Tech Science Press CMES, vol.71, no., pp.149-155, 011 Turbulentlike Quantitative Analysis on Energy Dissipation in Vibrated Granular Media Zhi Yuan Cui 1, Jiu Hui Wu 1 and Di Chen Li 1 Abstract:
More informationThe Effect of Model Internal Flexibility Upon NEMD Simulations of Viscosity
Draft: September 29, 1999 The Effect of Model Internal Flexibility Upon NEMD Simulations of Viscosity N. G. Fuller 1 and R. L. Rowley 1,2 Abstract The influence of model flexibility upon simulated viscosity
More informationContinuum Model of Avalanches in Granular Media
Continuum Model of Avalanches in Granular Media David Chen May 13, 2010 Abstract A continuum description of avalanches in granular systems is presented. The model is based on hydrodynamic equations coupled
More informationAuthors: Correspondence: ABSTRACT: Keywords:
Implementation of a material model with shear rate and temperature dependent viscosity Authors: Mathias Vingaard, Benny Endelt, Jesper declaville Christiansen Department of Production Aalborg University
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 informationExplicit algebraic Reynolds stress models for boundary layer flows
1. Explicit algebraic models Two explicit algebraic models are here compared in order to assess their predictive capabilities in the simulation of boundary layer flow cases. The studied models are both
More informationStochastic Particle Methods for Rarefied Gases
CCES Seminar WS 2/3 Stochastic Particle Methods for Rarefied Gases Julian Köllermeier RWTH Aachen University Supervisor: Prof. Dr. Manuel Torrilhon Center for Computational Engineering Science Mathematics
More informationScale invariant fluid dynamics for the dilute Fermi gas at unitarity
Scale invariant fluid dynamics for the dilute Fermi gas at unitarity Thomas Schaefer North Carolina State University Fluids: Gases, Liquids, Plasmas,... Hydrodynamics: Long-wavelength, low-frequency dynamics
More informationThe Polymers Tug Back
Tugging at Polymers in Turbulent Flow The Polymers Tug Back Jean-Luc Thiffeault http://plasma.ap.columbia.edu/ jeanluc Department of Applied Physics and Applied Mathematics Columbia University Tugging
More informationDivergence of the gradient expansion and the applicability of fluid dynamics Gabriel S. Denicol (IF-UFF)
Divergence of the gradient expansion and the applicability of fluid dynamics Gabriel S. Denicol (IF-UFF) arxiv:1608.07869, arxiv:1711.01657, arxiv:1709.06644 Frankfurt University 1.February.2018 Preview
More informationPart III. Polymer Dynamics molecular models
Part III. Polymer Dynamics molecular models I. Unentangled polymer dynamics I.1 Diffusion of a small colloidal particle I.2 Diffusion of an unentangled polymer chain II. Entangled polymer dynamics II.1.
More informationNon-Newtonian Fluids and Finite Elements
Non-Newtonian Fluids and Finite Elements Janice Giudice Oxford University Computing Laboratory Keble College Talk Outline Motivating Industrial Process Multiple Extrusion of Pastes Governing Equations
More informationViscous Fluids. Amanda Meier. December 14th, 2011
Viscous Fluids Amanda Meier December 14th, 2011 Abstract Fluids are represented by continuous media described by mass density, velocity and pressure. An Eulerian description of uids focuses on the transport
More informationRheology of two- and three-dimensional granular mixtures under uniform shear flow: Enskog kinetic theory versus molecular dynamics simulations
Rheology of two- and three-dimensional granular mixtures under uniform shear flow: Enskog kinetic theory versus molecular dynamics simulations José María Montanero Departamento de Electrónica e Ingeniería
More informationMicromechanics of Colloidal Suspensions: Dynamics of shear-induced aggregation
: Dynamics of shear-induced aggregation G. Frungieri, J. Debona, M. Vanni Politecnico di Torino Dept. of Applied Science and Technology Lagrangian transport: from complex flows to complex fluids Lecce,
More informationShear viscosity of liquid rubidium at the triple point
J. Phys. F: Met. Phys. 18 (1988) 1439-1447. Printed in the UK Shear viscosity of liquid rubidium at the triple point P T Cummingsi- and G P Morriss-l: i Department of Chemical Engineering, Thornton Hall,
More informationHydrodynamics. l. mahadevan harvard university
Thermodynamics / statistical mechanics : Hydrodynamics N>>1; large number of interacting particles! Hydrodynamics : L >> l; T >> t; long wavelength, slow time - average over (some) microscopic length and
More informationAvailable online at ScienceDirect. Procedia Engineering 103 (2015 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 103 (2015 ) 237 245 The 13 th Hypervelocity Impact Symposium Strength of Granular Materials in Transient and Steady State Rapid
More informationCH.9. CONSTITUTIVE EQUATIONS IN FLUIDS. Multimedia Course on Continuum Mechanics
CH.9. CONSTITUTIVE EQUATIONS IN FLUIDS Multimedia Course on Continuum Mechanics Overview Introduction Fluid Mechanics What is a Fluid? Pressure and Pascal s Law Constitutive Equations in Fluids Fluid Models
More informationOn the Asymptotic Convergence. of the Transient and Steady State Fluctuation Theorems. Gary Ayton and Denis J. Evans. Research School Of Chemistry
1 On the Asymptotic Convergence of the Transient and Steady State Fluctuation Theorems. Gary Ayton and Denis J. Evans Research School Of Chemistry Australian National University Canberra, ACT 0200 Australia
More information