Shear viscosity of liquid rubidium at the triple point
|
|
- Margery Greer
- 6 years ago
- Views:
Transcription
1 J. Phys. F: Met. Phys. 18 (1988) 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, University of Virginia, Charlottesville. VA 22901, USA $ Research School of Chemistry, Australian National University, GPO Box 4, Canberra, ACT Australia Received 11 August in final form 5 January 1988 Abstract. In a recent paper we reported non-equilibrium molecular dynamics simulations of liquid rubidium using the isokinetic sllod algorithm for the dynamics and the Price potential for the interactions. New and more extensive results are presented which confirm that the results exhibit normal square-root behaviour of the shear viscosity as a function of strain rate at the triple point. I. Introduction The non-equilibrium molecular dynamics (NEMD) technique has developed into a powerful tool for the study of transport properties of both simple and molecular fluids (Hoover and Ashurst 1975, Evans and Morriss 1984, Evans 1986, Evans and Hoover 1986). Recent developments, such as the sllod algorithm for shear viscosity and the derivation of synthetic techniques for calculating such properties as the thermal conductivity and self-diffusion, the use of Gauss s principle of least constraint as a tool for constraining simulations in physically meaningful ways (yielding, for example, isokinetic and/or isobaric ensemble simulations) (Hoover er af 1982, Evans and Morriss 1983) and the extension of linear response theory to test the validity of NEMD techniques (Morriss and Evans 1985), are reviewed by Evans (1986). In a recent paper (Cummings and Morriss 1987) we reported non-equilibrium molecular dynamics simulations of liquid rubidium using the isokinetic sllod algorithm for the dynamics and the Price potential for the interactions. In that paper, some concern was expressed about the results obtained at the triple point since the simulation results for shear viscosity available at the time appeared to deviate a little from the expected square-root dependence on the strain rate (see equation (6) below). This raised the possibility that at the triple point density and temperature in the absence of shear the equilibrium state of the simulated system (rubidium pseudoatoms interacting via the Price potential) is a solid rather than a triple point liquid. In the latter case, one expects to observe separate domains of square-root dependence in the shear viscosity associated with Bingham plastic-type behaviour (Evans 1982). It was conjectured that this is not the case and that, within the true error bars for the simulations, the results did satisfy the square-root dependence. Evidence confirming this conjecture is reported in this paper. The simulations reported in the earlier paper were performed on 108 molecules at five different strain /88/ IOP Publishing Ltd 1439
2 1440 P T Cummings and G P Morriss rates; in this paper, we present results for five additional strain rates with 108 molecules and results for six strain rates obtained using 256 molecules. Collectively, these results confirm that the square-root dependence of the shear viscosity is evident in the simulations and that the simulated system is in the liquid state at the experimental triple point. In Q 2, we briefly review the NEMD simulation method and in 0 3 we present results obtained at the triple point. 2. NEMD simulation method The sllod equations of motion induce steady-state planar Couette flow in a homogeneous system such that the streaming velocity has a non-zero component in the x direction, U,, which satisfies du,ldy = y where y is the constant strain rate. This corresponds to the linear velocity profile that distinguishes planar Couette flow (Bird et a1 1960). In order to do this, the positions r, and peculiar momentap, of the molecules in the simulation are obtained by solving the first-order isokinetic sllod equations of motion given by (Evans and Morriss 1984, Cummings and Morriss 1987) dr,ldt = p/m + r, Vu (20) dp,ldt = F, -pi - VU - Ap, (1) (2b) combined with the Lees-Edwards sliding brick boundary conditions (Lees and Edwards 1972). In equation (2), u=(yy, 0, 0) is the streaming velocity of the molecules, F, is the force on molecule i and A is given by A = (PI Ff -pip, Vu) p, * p,. I= I I, =I In this equation, indicates the full contraction of two second-order tensors and N is the number of particles (pseudoatoms in the current context) in the simulation cell. The formula for A is dictated by requiring that the temperature T be held constant, which is equivalent to demanding that E:; =,p: remain a constant of the motion. The yx component, Pyx, of the pressure tensor P is related to the velocity profile by the constitutive relation (3) where rj is the shear viscosity. The pressure tensor P, which is of second rank, is calculated from the expression N N N I= I r=l ]=I where m is the mass of each molecule, U, is the velocity of molecule i, V is the volume of the system, rf, = r, - r, and F,, is the force exerted on molecule i by molecule j and so
3 Liquid rubidium ut the triple point 1441 is given by the gradient of the pair potential q(r) between the molecules. Notice that the velocities in equation (5) have the streaming part U subtracted out and are thus peculiar velocities (Evans and Morriss 1984). The shear viscosity is therefore calculated from equation (4) using the time average of the yx component of P and the hydrostatic pressure P is the time average of the trace of P. The strain-rate-dependent q, P and configurational internal energy Uconf are then fitted to q = Til + VI Y I P= Pll+ Ply :? (7) win! = U,, + U,?3i? (8) where the subscript 0 indicates the zero strain rate limit of each quantity. Equations (6)-(8) are the leading-order strain rate dependencies derived theoretically (Kawasaki and Gunton 1973, Yamada and Kawasaki 1975, Ernst er ul 1978) and confirmed by extensive simulations (Evans 1979, 1981, 1983, Evans and Hanley 1980). The Newtonian viscosity of rubidium is given by q(,. The square-root dependence of the shear viscosity on strain rate, equation (6), is a consequence of the long-ranged t-? behaviour of the stress autocorrelation function at large time r. There has been some recent discussion as to whether the rrd? dependence (where d is the dimensionality of the system) is the true asymptotic longtime behaviour of the shear stress autocorrelation function (Dorfman and Kirkpatrick 1986, Kirkpatrick 1984, Kirkpatrick and Nieuwondt 1986, van Beijeren 1984). If this enhanced long-time tail behaviour is not the final asymptotic form then the small-y behaviour of the viscosity may deviate from its usual square-root dependence. Very recently Evans and Morriss (1987) have shown that at sufficiently small values of the shear rate, the viscosity becomes an analytic function, varying as y. This analytic behaviour is observed at shear rates that are much lower than those typically accessible in non-equilibrium molecular dynamics simulations. In practice, this means that the extrapolation to zero strain rate (required to obtain Newtonian viscosity) based on equation (6) may, in some cases, lead to an overestimate of the exact result. (6) 3. NEMD simulations at the triple point As in our previous paper (Cummings and Morriss 1987), all the simulations were performed on a Computer Signal Processing Inc (CSPI) bit word array processor using 108 or 256 rubidium ions using the Price potential (Price et a1 1970, Price 1971) for the interaction between rubidium pseudoatoms. The potential and force were assumed to be zero for r>3a and the potential shifted to make it continuous at the cut-off 3a. Other mechanical details of the simulation method used are described in P 2 of Cummings and Morriss (1987). In the earlier paper (Cummings and Morriss 1987) simulation results at the triple point (density p= g cm-, temperature T=318.9 K) were obtained for 108 rubidium ions at five reduced strain rates y* = yo(m/e) (equal to 0.1,O. 16, 0.36,0.64 and 1.0). In this paper, we report additional results for 108 molecules at strain rates y* of 0.01, 0.04, 0.25, 0.49 and 0.81 as well as results for 256 pseudoatom simulations at strain rates of 0.04, 0.09, 0.16, 0.36, 0.64 and 1.0. The 108 pseudoatom simulation results are tabulated in table 1 along with those from the earlier paper; the 256 pseudoatom results are given in table 2. All the shear viscosity results are plotted
4 ~ 1442 P T Cummings and G P Morriss Table 1. NEMD simulation results for liquid rubidium at the triple point (T= K and p= g ~m-~) using 108 rubidium pseudoatoms. Time steps y' ~(cp) LF"'/NE Pdls 80 OOO 80 OOO 160 OOO 100 OOO 240 OOO 100 OOO 400 OOO 360 OOO 360 OOO 240 OOO 1.o against y*"' in figure 1 and the pressure and internal energy results against Y"~'' in figure's 2 and 3 respectively. In tables 1 and 2, no numerical estimates of the errors in the shear viscosity are provided. This is in contrast to the previous paper (Cummings and Morriss 1987) in which the errors in the simulation quantities were estimated based on the standard deviation of block averages. As was pointed out in the previous paper, such error estimation is strongly dependent on the block size chosen. Rather than give numerical values of the error in the simulated quantities, we provide the following estimates as guides to the accuracy of the simulation results: the error in the shear viscosity varies from less than 1% for the high strain rate results (f' = 1) to 2-3% at the lowest strain rate; the error in the configurational internal energy and the pressure is 1% or less at all strain rates. The combined (108 and 256 pseudoatom) simulation results are fitted in a leastsquares sense to the strain rate equations (6)-(8). (In computing the least-squares fit, the results for 108 and 256 pseudoatom systems at strain rates 1.O, 0.64, 0.36, 0.16, 0.09 and 0.04 are averaged together to give a single value for 7, P and U at each of these strain rates.) The resulting fitted equations are ~= ~:'"' Pdl~ = ~"~" V'"'INE= t 0.418~*.~". (9) Table 2. NEMD simulation results for liquid rubidium at the triple point (T= K and p= g cm-') using 256 rubidium pseudoatoms. Time steps 20 OOO 40 OOO 4oOOO 40 OOO 40 OOO 80 OOO Y* 1.o tl (CP) LFOnfINe Pdle
5 Liquid rubidium at the triple point L Figure I. The shear viscosity q of liquid rubidium at the triple point as a function of the square root of the reduced strain rate y*. In all the figures, the circles (0) and triangles (A) represent the 108 and 256 pseudoatom simulation results respectively. The straight line is the least-squares fit of the simulation data to equation (6). Thus, the combined simulation results predict that the Newtonian viscosity of liquid rubidium at the triple point is cp. This is within 3.5% of the result of CP obtained using only the limited data from the previous paper (Cummings and Morriss 1987) and within 7% of the experimental result CP (Landolt-Bornstein 1969). It is our contention that the combined simulation results plotted in figure 1 exhibit the expected square-root dependence of the shear viscosity on strain rate, equation (6). There is greater regularity in the deviation of the 108 pseudoatom results from the square-root dependence (except at the lowest strain rate) than in the deviation of the 256 pseudoatom results. Nevertheless, all the simulation results cluster around the Figure 2. The dimensionless trace of the pressure tensor-see equation (5)-as a function of the cube of the square root of the strain rate y" for liquid rubidium at the triple point. The straight line is the least-squares fit of the simulation data to equation (7).
6 1444 P T Cummings and G P Morriss s a Figure3. The dimensionless configurational internal energy as a function of the cube of the square root of the strain rate y" for liquid rubidium at the triple point. The straight line is the least-squares fit of the simulation data to equation (8). square-root behaviour rather well, considering the large change in the viscosity over the range of simulated strain rates. One might be tempted to view the regularity in the deviation of the 108 pseudoatom results as being indicative of a dependence on strain rate that differs from equation (6). In particular, one might question two aspects of the asymptotic behaviour given in equation (6): first, is the large-y'l behaviour asymptotic in y":"'or is the power law different from 1/2? Second, is the small-y" result (q-y'") at all evident in the simulation? To answer these questions, the simulation results were leastsquares fitted to the following equations: q = qi, + q p l + q2y*2 (10) Equation (12) attempts to determine whether equation (6) is valid. If it is valid, one can expect that a will be near 1/2. Equation (1 1) is intended to determine whether the deviation from square-root behaviour in the simulation results is evidence of the transition to the classical small-y" behaviour. If equation (6) is valid, one can expect to see a small q? coefficient compared with the 7, coefficient (i.e., q2/ql<<l). Equation (10) is the most general fit of the simulation data that is considered here: again, validity of equation (6) would correspond to an a in this equation of approximately 112 and a small qrjq, ratio. The results of the least-squares fits are given in table 3 along with the squared error E' defined by I1 In this equation, the q,"""' are the simulation results for the shear viscosity (of which there are eleven when the 108 and 256 pseudoatom results at the same strain rate are combined), the y,* are the reduced strain rates at which the simulations were
7 Liquid rubidium at the triple point 1445 Table 3. Results obtained by least-squares fitting of simulation results to equations (10)- (12) and (6). Equation vi, 'I1 'I2 U E' (10) ( I.2 x 10. (11) X 10.' (12) W x 10 -' (6) I X 10 ' performed and fhci"! is the expression for the shear viscosity obtained from equations ( IO)-( 12) and (6) as appropriate. Simply on the grounds of flexibility, one can expect that equation (12) will exhibit the lowest squared error and equation (6) the highest. Only large differences between the squared errors can be taken as evidence that the asymptotic form of equation (6) is invalid. The resulting curves are plotted along with the simulation results in figure 4. It is evident that none of the curves fits the simulation results dramatically better than equation (6). In fact, equation (12) yields an exponent ( ) that is very close to 1/2. indicating that if the simulation data are dominated by a single power-law decay l a 0.2 Y - E I I I I , Figure4. The shear viscosity 7 of liquid rubidium at the triple point as a function of the reduced strain rate f. The full and dotted curves represent the least-squares fits of equations (IO) and (I 1) respectively in the upper figure. and of equations (12) and (6) respectively in the lower figure.
8 1446 P T Cummings and G P Morriss then the power is very close to 1/2. The equation that might be regarded as having the most theoretical legitimacy, equation (ll), fits the data well and moreover indicates that the correction to equation (6) is small. Perhaps more significant than any single fitted equation is that each equation extrapolates to a Newtonian viscosity in the range cp. This suggests strongly that the value obtained by extrapolating equation (6), 0.690cP, represents a good estimate of the zero strain rate viscosity implied by the simulation results. 4. Conclusions Additional simulation results at the triple point of liquid rubidium have been presented and analysed. It is clear from these results that the behaviour of the strainrate-dependent shear viscosity at the triple point is not anomalous (in this context, exhibiting Bingham plastic behaviour) and the square-root behaviour is that typically seen in non-equilibrium simulations at this range of values of strain rate. Based on the four functional forms used to extrapolate the simulation results to zero strain rate, the Price potential predicts a viscosity of liquid rubidium at the triple point of 0.69k0.04 cp. Acknowledgments PTC acknowledges the support of this research by the National Science Foundation (grant CBT ). The authors thank the National Science Foundation and the Australian Department and Science and Technology which, through the US/Australia Collaborative Research Program, provided travel funds which facilitated this research. The simulations reported in this paper were performed on a CSPI bit word array processor in the Center for Computer Aided Engineering at the University of Virginia. The authors are indebted to the Center for the provision of computational facilities in support of the array processor. The support of the National Science Foundation in providing matching funds (CPE ) for the purchase of the array processor to support research in statistical mechanics and chemical process dynamics is gratefully acknowledged. The authors are indebted to one of the referees who suggested fitting the shear viscosity results to functional forms other than equation (6), thus resulting in the discussion in the last three paragraphs of $3. References van Beijeren H 1984 Phys. Lett. 105A 191 Bird R B, Stewart W E and Lightfoot E N 1960 Transport Phenomena (New York: Wiley) Cummings P T and Morriss G P 1987 J. Phys. F: Mer. Phys Dorfman J R and Kirkpatrick T R 1986 Moleculur Dynumics Simuluiio,7 of SIuris/icul Mechunicul Sysrems ed. G Ciccotti and W G Hoover (Amsterdam: North-Holland) Ernst M H, Cichocki B, Dorfman J R, Sharma J and van Beijeren H 1978 J. Stat. Phys Evans D J 1979 Phys. Lett. 74A Phys. Rev. A Phys. Rev. A Molecular Dynamics Simulation of Statistical Mechanical Systems ed. G Ciccotti and W G Hoover (Amsterdam: North-Holland)
9 Liquid rubidium ut the triple point 1447 Evans D J and Hanley H J M 1980 Physica A Evans D J and Hoover W G 1986 Ann. Rev. Fluid Mech Evans D J and Morriss G P 1983 Chem. Phys Comput. Phys. Rep Phys. Rev. A Hoover W G and Ashurst W T 1975 Theor. Chem. Adu. Perspect. 1 1 Hoover W G, Ladd A J C and Moran B 1982 Phys. Rev. Lett Kawasaki K and Gunton J D 1973 Phys. Rev. A Kirkpatrick T R 1984 Phys. Rev. Len Kirkpatrick T R and Nieuwondt J 1986 Phys. Rev. A Landolt-Bornstein 1969 Zahlenwerte und Funktionem 5 Tiel (Berlin: Springer) Lees A W and Edwards S F 1972 J. Phys. C: Solid State Phys Morriss G P and Evans D E 1985 Mol. Phys Price D L 1971 Phys. Rev. A Price D L, Singwi K S and Tosi M P 1970 Phys. Rev. B Yamada T and Kawasaki K 1975 Prog. Theor. Phys
Non-equilibrium molecular dynamics calculation of the shear viscosity of liquid rubidium
J. Phys. F: Met. Phys. 17 (1987) 593-604. Printed in the UK Non-equilibrium molecular dynamics calculation of the shear viscosity of liquid rubidium P T Cummingsf and G P Morrisst + Department of Chemical
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 informationShear viscosity of model mixtures by nonequilibrium molecular dynamics. I. Argon-krypton mixtures
Shear viscosity of model mixtures by nonequilibrium molecular dynamics. I. Argon-krypton mixtures Song Hi Leea) and Peter T. Cummingsb) Department of ChemicaI Engineering, Thornton Hall, University of
More informationStatistical 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 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 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 informationThe viscosity-radius relationship from scaling arguments
The viscosity-radius relationship from scaling arguments D. E. Dunstan Department of Chemical and Biomolecular Engineering, University of Melbourne, VIC 3010, Australia. davided@unimelb.edu.au Abstract
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 informationNon-Equilibrium Molecular Dynamics Approach to the Rheology of Model Polymer Fluids
Fluid Phase Equilibria, 88 (1993) 99-113 Elsevier Science Publishers B.V., msterdam 99 Non-Equilibrium Molecular Dynamics pproach to the Rheology of Model Polymer Fluids J. W. Rudisill and P. T. Cummings
More informationNote on the KaplanYorke Dimension and Linear Transport Coefficients
Journal of Statistical Physics, Vol. 101, Nos. 12, 2000 Note on the KaplanYorke Dimension and Linear Transport Coefficients Denis J. Evans, 1 E. G. D. Cohen, 2 Debra J. Searles, 3 and F. Bonetto 4 Received
More informationarxiv: v1 [cond-mat.stat-mech] 26 Jan 2008
Long Time Tail of the Velocity Autocorrelation Function in a Two-Dimensional Moderately Dense Hard Disk Fluid Masaharu Isobe Graduate School of Engineering, Nagoya Institute of Technology, Nagoya 466-8555,
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 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 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 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 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 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 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 informationarxiv:nlin/ v1 [nlin.cd] 8 May 2001
Field dependent collision frequency of the two-dimensional driven random Lorentz gas arxiv:nlin/0105017v1 [nlin.cd] 8 May 2001 Christoph Dellago Department of Chemistry, University of Rochester, Rochester,
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 informationSupplementary Figure 1: Approach to the steady state. Energy vs. cycle for different
( = ) -6.9-6.94 (a) N = 8 < y.2 (T = 1).4 (T = 1).6 (T = 1).7 (T = 1) ( = ) -6.94 (b).8 (T = 1).9 (T = 1).12 (T = 1).14 (T = 1).8 (T =.466).9 (T =.466).12 (T =.466).14 (T =.466) > y n * 12 8 (c) N = 8
More informationCessation of Couette and Poiseuille flows of a Bingham plastic and finite stopping times
J. Non-Newtonian Fluid Mech. 29 2005) 7 27 Cessation of Couette and Poiseuille flows of a Bingham plastic and finite stopping times Maria Chatzimina a, Georgios C. Georgiou a,, Ioannis Argyropaidas b,
More informationInfluence of the Friedel oscillations on the structure of liquid rubidium
J. Phys. F: Met. Phys.. 12(1982)233-44 Printed in Great Britain Influence of the Friedel oscillations on the structure of liquid rubidium P T Cummingst and P A Egelstaff Department of Physics, University
More informationChemical and Biomolecular Engineering 150A Transport Processes Spring Semester 2017
Chemical and Biomolecular Engineering 150A Transport Processes Spring Semester 2017 Objective: Text: To introduce the basic concepts of fluid mechanics and heat transfer necessary for solution of engineering
More informationStress Overshoot of Polymer Solutions at High Rates of Shear
Stress Overshoot of Polymer Solutions at High Rates of Shear K. OSAKI, T. INOUE, T. ISOMURA Institute for Chemical Research, Kyoto University, Uji, Kyoto 611-0011, Japan Received 3 April 2000; revised
More informationNavier-Stokes Equation: Principle of Conservation of Momentum
Navier-tokes Equation: Principle of Conservation of Momentum R. hankar ubramanian Department of Chemical and Biomolecular Engineering Clarkson University Newton formulated the principle of conservation
More informationInstitute for Condensed Matter Theory and Beckman Institute, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
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,
More informationStudies on flow through and around a porous permeable sphere: II. Heat Transfer
Studies on flow through and around a porous permeable sphere: II. Heat Transfer A. K. Jain and S. Basu 1 Department of Chemical Engineering Indian Institute of Technology Delhi New Delhi 110016, India
More informationTowards a computational chemical potential for nonequilibrium steady-state systems
PHYSICAL REVIEW E VOLUME 60, NUMBER 5 NOVEMBER 1999 Towards a computational chemical potential for nonequilibrium steady-state systems András Baranyai Department of Theoretical Chemistry, Eötvös University,
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 informationElements of Continuum Elasticity. David M. Parks Mechanics and Materials II February 25, 2004
Elements of Continuum Elasticity David M. Parks Mechanics and Materials II 2.002 February 25, 2004 Solid Mechanics in 3 Dimensions: stress/equilibrium, strain/displacement, and intro to linear elastic
More informationAccelerating energy release prior to large events in simulated earthquake cycles: implications for earthquake forecasting
Accelerating energy release prior to large events in simulated earthquake cycles: implications for earthquake forecasting Peter Mora and David Place QUAKES, Department of Earth Sciences, The University
More informationChapter 5. The Differential Forms of the Fundamental Laws
Chapter 5 The Differential Forms of the Fundamental Laws 1 5.1 Introduction Two primary methods in deriving the differential forms of fundamental laws: Gauss s Theorem: Allows area integrals of the equations
More informationLaminar Forced Convection Heat Transfer from Two Heated Square Cylinders in a Bingham Plastic Fluid
Laminar Forced Convection Heat Transfer from Two Heated Square Cylinders in a Bingham Plastic Fluid E. Tejaswini 1*, B. Sreenivasulu 2, B. Srinivas 3 1,2,3 Gayatri Vidya Parishad College of Engineering
More informationViscoelastic Flows in Abrupt Contraction-Expansions
Viscoelastic Flows in Abrupt Contraction-Expansions I. Fluid Rheology extension. In this note (I of IV) we summarize the rheological properties of the test fluid in shear and The viscoelastic fluid consists
More informationInfluence of strain on transport in dense Lennard-Jones systems
Influence of strain on transport in dense Lennard-Jones systems Janka Petravic Citation: The Journal of Chemical Physics 120, 7041 (2004); doi: 10.1063/1.1683073 View online: http://dx.doi.org/10.1063/1.1683073
More informationPerformance evaluation of different model mixers by numerical simulation
Journal of Food Engineering 71 (2005) 295 303 www.elsevier.com/locate/jfoodeng Performance evaluation of different model mixers by numerical simulation Chenxu Yu, Sundaram Gunasekaran * Food and Bioprocess
More informationLong-time behavior and different shear regimes in quenched binary mixtures
PHYSICAL REVIEW E 75, 011501 2007 Long-time behavior and different shear regimes in quenched binary mixtures G. Gonnella* Dipartimento di Fisica, Università di Bari, and INFN, Sezione di Bari, Via Amendola
More informationChapter 9: Differential Analysis
9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control
More informationEvan Mitsoulis. School of Mining Engineering and Metallurgy. National Technical University of Athens. R.R. Huilgol
Cessation of Couette and Poiseuille ows of a Bingham plastic and nite stopping times Maria Chatzimina, Georgios C. Georgiou Department of Mathematics and Statistics, University of Cyprus P.O. Box 2537,
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 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 information2 Lyapunov Exponent exponent number. Lyapunov spectrum for colour conductivity of 8 WCA disks.
3 2 Lyapunov Exponent 1-1 -2-3 5 1 15 exponent number Lyapunov spectrum for colour conductivity of 8 WCA disks. 2 3 Lyapunov exponent 2 1-1 -2-3 5 1 15 exponent number Lyapunov spectrum for shear flow
More informationLecture 6 non-equilibrium statistical mechanics (part 2) 1 Derivations of hydrodynamic behaviour
Lecture 6 non-equilibrium statistical mechanics (part 2) I will outline two areas in which there has been a lot of work in the past 2-3 decades 1 Derivations of hydrodynamic behaviour 1.1 Independent random
More informationChapter 9: Differential Analysis of Fluid Flow
of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known
More informationUvA-DARE (Digital Academic Repository)
UvA-DARE (Digital Academic Repository) Poiseuille flow to measure the viscosity of particle model fluids. Backer, J.A.; Lowe, C.P.; Hoefsloot, H.C.J.; Iedema, P.D. Published in: Journal of Chemical Physics
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 informationCENG 501 Examination Problem: Estimation of Viscosity with a Falling - Cylinder Viscometer
CENG 501 Examination Problem: Estimation of Viscosity with a Falling - Cylinder Viscometer You are assigned to design a fallingcylinder viscometer to measure the viscosity of Newtonian liquids. A schematic
More informationUnsteady Boundary Layer Flow and Symmetry Analysis of a Carreau Fluid
Mathematics and Statistics: Open Access Received: Dec, 5, Accepted: Jan 5, 6, Published: Jan 8, 6 Math Stat, Volume, Issue http://crescopublications.org/pdf/msoa/msoa--.pdf Article Number: MSOA-- Research
More informationCessation of annular Poiseuille flows of Bingham plastics
J. Non-Newtonian Fluid Mech. 142 2007) 135 142 Cessation of annular Poiseuille flows of Bingham plastics Maria Chatzimina a, Christos Xenophontos a, Georgios C. Georgiou a,, Ioannis Argyropaidas b, Evan
More informationDifferential relations for fluid flow
Differential relations for fluid flow In this approach, we apply basic conservation laws to an infinitesimally small control volume. The differential approach provides point by point details of a flow
More informationComputational Fluid Dynamical Simulations of Droplet Flow and Merger in Microfluidic Channels
Computational Fluid Dynamical Simulations of Droplet Flow and Merger in Microfluidic Channels Katrina J. Donovan 1, Shelly Gulati 2, Xize Niu 3, Bradley M. Stone 1 and Andrew J. demello 4 1. Department
More informationFrom the last time, we ended with an expression for the energy equation. u = ρg u + (τ u) q (9.1)
Lecture 9 9. Administration None. 9. Continuation of energy equation From the last time, we ended with an expression for the energy equation ρ D (e + ) u = ρg u + (τ u) q (9.) Where ρg u changes in potential
More informationEvaluation of the CPY and PYX approximations for short ranged anisotropie potentials
MOLECULAR PHYSICS, 1983, VOL. 50, NO. 5, 1133-1140 Evaluation of the CPY and PYX approximations for short ranged anisotropie potentials by P. T. CUMMINGS t Departments of Mechanical Engineering and Chemistry,
More informationMolecular dynamics simulation of limiting conductances for LiCl, NaBr, and CsBr in supercritical water
JOURNAL OF CHEMICAL PHYSICS VOLUME 112, NUMBER 2 8 JANUARY 2000 Molecular dynamics simulation of limiting conductances for LiCl, NaBr, and CsBr in supercritical water S. H. Lee Department of Chemistry,
More informationExperimental test of an expression for the decay of an autocorrelation function
Supplemental Material for: Experimental test of an expression for the decay of an autocorrelation function Journal: Physical Review Letters Authors: Zach Haralson and J. Goree I. Experimental setup The
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 informationPharmaceutics I. Unit 6 Rheology of suspensions
Pharmaceutics I اينالديصيدلينيات 1 Unit 6 Rheology of suspensions 1 Rheology, the science of the flow or deformation of matter (liquid or soft solid) under the effect of an applied force. It addresses
More informationRock Rheology GEOL 5700 Physics and Chemistry of the Solid Earth
Rock Rheology GEOL 5700 Physics and Chemistry of the Solid Earth References: Turcotte and Schubert, Geodynamics, Sections 2.1,-2.4, 2.7, 3.1-3.8, 6.1, 6.2, 6.8, 7.1-7.4. Jaeger and Cook, Fundamentals of
More informationME 262 BASIC FLUID MECHANICS Assistant Professor Neslihan Semerci Lecture 4. (Buoyancy and Viscosity of water)
ME 262 BASIC FLUID MECHANICS Assistant Professor Neslihan Semerci Lecture 4 (Buoyancy and Viscosity of water) 16. BUOYANCY Whenever an object is floating in a fluid or when it is completely submerged in
More informationChapter 1 Introduction
Chapter 1 Introduction This thesis is concerned with the behaviour of polymers in flow. Both polymers in solutions and polymer melts will be discussed. The field of research that studies the flow behaviour
More informationStress, Strain, and Viscosity. San Andreas Fault Palmdale
Stress, Strain, and Viscosity San Andreas Fault Palmdale Solids and Liquids Solid Behavior: Liquid Behavior: - elastic - fluid - rebound - no rebound - retain original shape - shape changes - small deformations
More informationEffective Scraping in a Scraped Surface Heat Exchanger: Some Fluid Flow Analysis
ICEF 9 003 Effective Scraping in a Scraped Surface Heat Exchanger: Some Fluid Flow Analysis D. L. Pyle (1), K.-H. Sun (1), M. E. M. Lee (), C. P. Please (), A. D. Fitt (), S. K. Wilson (3), B. R. Duffy
More informationThe existence of Burnett coefficients in the periodic Lorentz gas
The existence of Burnett coefficients in the periodic Lorentz gas N. I. Chernov and C. P. Dettmann September 14, 2006 Abstract The linear super-burnett coefficient gives corrections to the diffusion equation
More informationDynamics of sheared gases
Computer Physics Communications 121 122 (1999) 225 230 www.elsevier.nl/locate/cpc Dynamics of sheared gases Patricio Cordero a,1, Dino Risso b a Departamento de Física, Facultad de Ciencias Físicas y Matemáticas,
More informationLecture 3. Properties of Fluids 11/01/2017. There are thermodynamic properties of fluids like:
11/01/2017 Lecture 3 Properties of Fluids There are thermodynamic properties of fluids like: Pressure, p (N/m 2 ) or [ML -1 T -2 ], Density, ρ (kg/m 3 ) or [ML -3 ], Specific weight, γ = ρg (N/m 3 ) or
More informationDetailed Outline, M E 521: Foundations of Fluid Mechanics I
Detailed Outline, M E 521: Foundations of Fluid Mechanics I I. Introduction and Review A. Notation 1. Vectors 2. Second-order tensors 3. Volume vs. velocity 4. Del operator B. Chapter 1: Review of Basic
More informationD.R. Rector, M.L. Stewart and A.P. Poloski Pacific Northwest National Laboratory P.O. Box 999, Richland, WA
Modeling of Sediment Bed Behavior for Critical Velocity in Horizontal Piping 9263 D.R. Rector, M.L. Stewart and A.P. Poloski Pacific Northwest National Laboratory P.O. Box 999, Richland, WA ABSTRACT A
More informationLQSAlamos National Laboratory
DEVELOPMENT AND APPLICATION OF ADVANCED ONEPOINT TURBULENCE MODELS Charles G. Speziale January 1996 Technical Report No. AM96002 DISCLAIMER This report was prepared as an account of work sponsored by an
More informationOn the oscillatory tube flow of healthy human blood
On the oscillatory tube flow of healthy human blood Miguel A. Moyers-Gonzalez Department of Mathematics and Statistics, University of Canterbury, New Zealand Robert G. Owens Département de mathématiques
More informationImproved model of nonaffine strain measure
Improved model of nonaffine strain measure S. T. Milner a) ExxonMobil Research & Engineering, Route 22 East, Annandale, New Jersey 08801 (Received 27 December 2000; final revision received 3 April 2001)
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 informationCONTRIBUTION TO EXTRUDATE SWELL FROM THE VELOCITY FACTOR IN NON- ISOTHERMAL EXTRUSION
Second International Conference on CFD in the Minerals and Process Industries CSIRO, Melbourne, Australia 6-8 December 1999 CONTRIBUTION TO EXTRUDATE SWELL FROM THE VELOCITY FACTOR IN NON- ISOTHERMAL EXTRUSION
More informationFinite-size analysis via the critical energy -subspace method in the Ising models
Materials Science-Poland, Vol. 23, No. 4, 25 Finite-size analysis via the critical energy -subspace method in the Ising models A. C. MAAKIS *, I. A. HADJIAGAPIOU, S. S. MARTINOS, N. G. FYTAS Faculty of
More informationDepartment of Materials and Environmental Chemistry Arrhenius Laboratory, Stockholm University Stockholm, Sweden *
CMST ) 9 50 07) DOI:0.9/cmst.06.0000066 Minimal Energy Dissipation Rate and Director Orientation Relative to External Dissipative Fields such as Temperature and Velocity Gradients in Nematic and Cholestric
More informationInternational Journal of Pure and Applied Mathematics Volume 58 No ,
International Journal of Pure and Applied Mathematics Volume 58 No. 2 2010, 195-208 A NOTE ON THE LINEARIZED FINITE THEORY OF ELASTICITY Maria Luisa Tonon Department of Mathematics University of Turin
More informationC subjected to large deformation rates exhibit
D. K. DOEKSEN, R. J. J. JONGSCHAAP, and H. KAMPHUIS Rheology Group Department of Applied Physics Twente University of Technology Ensc hed e The Netherlands A kinetic model for the rheological behavior
More informationElastic and Plastic Behavior of Model Solids
TECHNISCHE MECHANIK, Band 22, Heft 2, (2002), 79-88 Manuskripteingang: 18. Januar 2002 Elastic and Plastic Behavior of Model Solids S. Hess, M. Kröger A short ranged attractive (SHRAT) potential is employed
More informationLinear Constitutive Relations in Isotropic Finite Viscoelasticity
Journal of Elasticity 55: 73 77, 1999. 1999 Kluwer Academic Publishers. Printed in the Netherlands. 73 Linear Constitutive Relations in Isotropic Finite Viscoelasticity R.C. BATRA and JANG-HORNG YU Department
More informationCombined Stresses and Mohr s Circle. General Case of Combined Stresses. General Case of Combined Stresses con t. Two-dimensional stress condition
Combined Stresses and Mohr s Circle Material in this lecture was taken from chapter 4 of General Case of Combined Stresses Two-dimensional stress condition General Case of Combined Stresses con t The normal
More informationPoiseuille Advection of Chemical Reaction Fronts
Utah State University DigitalCommons@USU All Physics Faculty Publications Physics 8- Poiseuille Advection of Chemical Reaction Fronts Boyd F. Edwards Utah State University Follow this and additional works
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 informationSurface stress and relaxation in metals
J. Phys.: Condens. Matter 12 (2000) 5541 5550. Printed in the UK PII: S0953-8984(00)11386-4 Surface stress and relaxation in metals P M Marcus, Xianghong Qian and Wolfgang Hübner IBM Research Center, Yorktown
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 informationCOMPUTATIONAL STUDY OF PARTICLE/LIQUID FLOWS IN CURVED/COILED MEMBRANE SYSTEMS
COMPUTATIONAL STUDY OF PARTICLE/LIQUID FLOWS IN CURVED/COILED MEMBRANE SYSTEMS Prashant Tiwari 1, Steven P. Antal 1,2, Michael Z. Podowski 1,2 * 1 Department of Mechanical, Aerospace and Nuclear Engineering,
More informationThe Mathematical Analysis for Peristaltic Flow of Hyperbolic Tangent Fluid in a Curved Channel
Commun. Theor. Phys. 59 213 729 736 Vol. 59, No. 6, June 15, 213 The Mathematical Analysis for Peristaltic Flow of Hyperbolic Tangent Fluid in a Curved Channel S. Nadeem and E.N. Maraj Department of Mathematics,
More informationElec Eng 3BA3: Structure of Biological Materials
Elec Eng 3BA3: Structure of Biological Materials Page 1 of 12 Day Class Instructor: Dr. I. C. BRUCE Duration of Examination: 3 Hours McMaster University Final Examination December 5, 2008 This examination
More informationMath Review: Vectors and Tensors for Rheological Applications
Math Review: Vectors and Tensors for Rheological Applications Presented by Randy H. Ewoldt University of Illinois at Urbana-Champaign U. of Minnesota Rheological Measurements Short Course June 2016 MathReview-1
More informationTHE SUBORDINATION OF THE THREE- DIMENSIONAL FLOW INSTALLATION IN THE CONVERGING CHANNEL ON RHEOLOGICAL CHARACTERISTICS OF POLYMER STREAM
International Journal of Civil Engineering and Technology (IJCIET) Volume 9, Issue 13, December 2018, pp. 949-956, Article ID: IJCIET_09_13_095 Available online at http://www.iaeme.com/ijciet/issues.asp?jtype=ijciet&vtype=9&itype=13
More informationThe Kawasaki Identity and the Fluctuation Theorem
Chapter 6 The Kawasaki Identity and the Fluctuation Theorem This chapter describes the Kawasaki function, exp( Ω t ), and shows that the Kawasaki function follows an identity when the Fluctuation Theorem
More informationNonlinear stability of steady flow of Giesekus viscoelastic fluid
Nonlinear stability of steady flow of Giesekus viscoelastic fluid Mark Dostalík, V. Průša, K. Tůma August 9, 2018 Faculty of Mathematics and Physics, Charles University Table of contents 1. Introduction
More informationRheological Properties
Rheological Properties Shear Stress the force experienced by the liquid is proportional to the area of the upper plate and is defined as the shear stress, σ = F/A (Units = Nm -2 or Pascals [Pa]) Shear
More informationChapter 1 Fluid Characteristics
Chapter 1 Fluid Characteristics 1.1 Introduction 1.1.1 Phases Solid increasing increasing spacing and intermolecular liquid latitude of cohesive Fluid gas (vapor) molecular force plasma motion 1.1.2 Fluidity
More informationWe may have a general idea that a solid is hard and a fluid is soft. This is not satisfactory from
Chapter 1. Introduction 1.1 Some Characteristics of Fluids We may have a general idea that a solid is hard and a fluid is soft. This is not satisfactory from scientific or engineering point of view. In
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 informationNonequilibrium molecular dynamics via Gauss's principle of least constraint
PHYSICAL REVIEW A VOLUME 28, NUMBER 2 AUGUST 1983 Nonequilibrium molecular dynamics via Gauss's principle of least constraint Denis J. Evans Research School o/chemistry, Australian National University,
More informationComposite angle ply laminates and netting analysis
10.1098/rspa.2002.1066 FirstCite e-publishing Composite angle ply laminates and netting analysis By J. T. Evans and A. G. Gibson School of Mechanical and Systems Engineering, University of Newcastle upon
More informationFALLING FILM FLOW ALONG VERTICAL PLATE WITH TEMPERATURE DEPENDENT PROPERTIES
Proceedings of the International Conference on Mechanical Engineering 2 (ICME2) 8-2 December 2, Dhaka, Bangladesh ICME-TH-6 FALLING FILM FLOW ALONG VERTICAL PLATE WITH TEMPERATURE DEPENDENT PROPERTIES
More informationDilatancy Transition in a Granular Model. David Aristoff and Charles Radin * Mathematics Department, University of Texas, Austin, TX 78712
Dilatancy Transition in a Granular Model by David Aristoff and Charles Radin * Mathematics Department, University of Texas, Austin, TX 78712 Abstract We introduce a model of granular matter and use a stress
More information