Non-equilibrium molecular dynamics calculation of the shear viscosity of liquid rubidium
|
|
- Charla York
- 6 years ago
- Views:
Transcription
1 J. Phys. F: Met. Phys. 17 (1987) 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 Engineering, Thornton Hall, University of Virginia. Charlottesville, VA 22901, USA Research School of Chemistry. Australian National University, GPO Box 4, Canberra, ACT 260 I, Australia Received 18 April in final form 25 June 1986 Abstract. We report non-equilibrium molecular dynamics simulations of liquid rubidium using the isokinetic sllod algorithm for the dynamics and the Price potential for the interactions. The three state points simulated are the same three chosen by Bansal and Bruns who used a quite different non-equilibrium method. However. our results differ substantially from theirs. Specifically, we find that the Price potential consistently predicts better results for the shear viscosity. A more fundamental disagreement between our results and theirs is that our results underestimate the experimental data where theirs overestimate and vice versa. This suggests that the method used by Bansal and Bruns has systematic errors. 1. Introduction Over the past decade, 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 1984a, Evans 1985, Evans and Hoover 1986). The most promising recent developments include the derivation of synthetic techniques for calculating such properties as the shear viscosity, thermal conductivity and self-diffusion and 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 et a1 1982, Evans and Morriss 1983)). Finally, the formal status of these and other NEMD techniques has been clarified tremendously in recent years through linear response theory, so that previously raised objections to the orthodoxy of the NEMD methods have now been effectively answered (Evans and Morriss 1984a, Morriss and Evans 1985). Of all the NEMD techniques, methods for calculating the shear viscosity have the longest history and therefore have come under the greatest scrutiny. The methodsreviewed in $ 3 of Evans and Morriss (1984a)-amount to setting up a steady state planar Couette flow in a homogeneous system. The streaming velocity has a non-zero component in the x direction, U,, which satisfies du,/d>x= y where y is the constant strain rate. Equation (1) implies that there is a linear velocity profile in the x direction, the characteristic of planar Couette flow (Bird et a1 1960). The J'S /87/ $ IOP Publishing Ltd 593 (1)
2 594 P T Cummings and G P Morriss component, Pyx, of the pressure tensor Pis related to the velocity profile by the constitutive relation PYx=-~(d~,/d.v)=-~~ (2) where for Newtonian fluids q, the shear viscosity, is dependent on the state variables (pressure, temperature T and composition, for example). Non-Newtonian fluids are fluids for which q depends on y in addition to the dependence on the thermodynamic state. The pressure tensor P, which is of second rank, is calculated from the expression Y. Y Y I = I where m is the mass of each molecule, U, is the velocity of molecule i, Vis the volume of the system, r!, = r, - r, with rk being the position of the centre of mass of molecule k and F, is the force exerted by molecule i on moleculej and so is given by the gradient of the pair potential q(r) between the molecules. Notice that the velocities in equation (3) have the streaming part U subtracted out and are thus peculiar velocities (Evans and Morriss 1984a). The shear viscosity is therefore calculated from equation (2) using the time average of the jix-component of P and the hydrostatic pressure is the time average of the trace of P. In order to perform a molecular dynamics simulation of planar Couette flow, one necessary modification of the usual formalism is that the images of the central simulation cell located above and below it must move in time in opposite directions at a constant rate determined by the strain rate 7, and the x-component of the velocity of a molecule leaving either the top or the bottom of the central cell in the y direction must be modified to be consistent with the overall velocity profile (the so-called sliding brick boundary conditions). These boundary conditions were first introduced by Lees and Edwards (1972) and in their original NEMD algorithm the velocity profile was set up and maintained solely by these moving boundaries, the dynamics of the particles in the simulation cell being given by Newton s equations of motion, dri pi - dt m _- dpi dt - Fi where pi is the laboratory momentum of molecule i and Fi is the force on molecule i due to all the other molecules in the system. A convenient reformulation in terms of first-order equations of motion (Evans and Morriss 1984b) is the sllod algorithm given by dri Pi - +ri.vu dt m -- dpi - Fi -pi * VU dt combined with the Lees-Edwards sliding brick boundary conditions. As far as the trajectories r,(t) of the individual molecules are concerned, the sllod algorithm can be shown to be equivalent to the Lees-Edwards algorithm except at time f=o since elimination of pi from equations (5a) and (5b) leads to Newton s equations of motion. However, from equations (5) it is clear that the sllod equations embody the velocity profile
3 Shear riscositv of liquid rubidium 5 95 Table I. A comparison of predicted and experimental shear viscosities for liquid rubidium. The theoretical predictions of this work and Bansal and Bruns (1984) are obtained using two different non-equilibrium molecular dynamics techniques. 1 S O I a This work. Bansal and Bruns (1984). Landolt-Bornstein (1 969). (since the gradient of the streaming velocity is added to the peculiar velocity pi/m). The sllod algorithm is thus homogeneous and consistent with the boundary conditions. Evans and Morriss (1984a) review the development of the sllod algorithm (as well as the closely related Doll s tensor method) and the reasons for their utility in non-equilibrium simulations of planar Couette flow; we refer the reader to their paper for details. Other methods for calculating the shear viscosity include the Green-Kubo relation (which requires long equilibrium molecular dynamics simulations of large systems for accuracy equivalent to the synthetic sllod algorithm (Holian and Evans 1983)), the differential trajectory method (Ciccotti and Jacucci 1975, Ciccotti et a1 1976, 1979) and the method of Singer and co-workers (Singer et a1 1980). It has recently been shown (Paolini et a1 1986) that the differential trajectory method is more reliable when the perturbation is a delta function rather than a step function (as used by Singer et a1 1980). This last method was used by Bansal and Bruns (1984) to determine the shear viscosity of liquid Rb using the Price potential (Price et a1 1970, Price 1971) at three state points: T= K, p= g ~ m-~, T= K, p= g ~ m - and ~ T= K, p= g cm-3 where T is the absolute temperature and p is the density. Their results (shown in table 1) suggest that the Price potential is only moderately accurate in predicting the transport properties of liquid Rb. In fact, according to Bansal and Bruns. the error in the shear viscosity predicted by the Price potential ranges from 32 to 45%. Given the success of the Price potential in predicting thermodynamic and static and dynamic structural properties of liquid Rb (Rahman 1974a, b, Mountain 1978, Arlinghaus and Cummings 1987), this discrepancy seems somewhat large and prompted the reexamination of the predictions of the Price potential given in this paper. Using the best available technique-that is, Gaussian thermostatted sllod dynamics-we find the results given in table 1, which indicate that the Price potential results range in error from 3 to 11%, a significant improvement over that of Bansal and Bruns (1984). In 8 2 of this paper, we describe the method of calculation, the results obtained and their analysis in more detail. Section 3 contains our conclusions. 2. NEMD simulation method The potential used in the simulations reported in this paper is the Price potential calculated according to the prescription given by Price and co-workers (Price et a1 1970, Price 1971). As with all effective interionic potentials derived using pseudopotential theory (Harrison 1966), the Price potential is density dependent; thus the potential differs at each of the
4 596 P T Cutnmings and G P Morriss three state points simulated due to the different densities involved. The NEMD equations solved are the isokinetic sllod equations given by where,i is given by where 0 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 1 is dictated by requiring that the temperature T is held constant, which is equivalent to demanding that EL I p: remain a constant of the motion. All the simulations were performed using 108 rubidium ions. The calculations were carried out on a Computer Signal Processing Inc. (CSPI) bit word array processor in the Center for Computer Aided Engineering at the University of Virginia. Since the potential and the force are not given analytically, they were discretised onto a grid of 900 points equally spaced in r2 corresponding to the range O<r<3a where Y is the interparticle separation. Both the force and potential were obtained by linear interpolation of this grid. 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. Those simulation runs which start from a crystalline configuration begin with a face-centred cubic lattice which is equilibrated over time steps. At each state point, planar Couette flow was simulated at five values of the shear rate y given in dimensionless form by y* = ya(m/~) /~ where a is the diameter of an ion which can be defined to be the interionic separation at which the Price potential is zero so that the potential is positive for all r < U (= 4.4 8, for the densities simulated), m is the mass of the ions (= x g) and E is the depth of the Price potential (which is taken to be c/kb =400 K, kb being Boltzmann s constant, at the three state points simulated). In each simulation, the dimensionless time step used in the fourth-order predictor-corrector integration routine (Evans and Morriss 1984a) is Ar* = Ar(c/ma2)1/2 = which corresponds to a real time step of x s. The values of )I* along with the number of time steps per value of y* and the calculated values of the shear viscosity are given in tables 2 to 4. Notice that at each state point the shear viscosity reduces with Table 2. NEMD simulation results for liquid rubidium at a temperature of K and a density of g cm-j. Least-squares fitting of the data to equation (9) gives 7= y* I2 cp. The zero strain rate value of CP compares with the experimental value of CP and the Bansal and Bruns (1984) value of o 0.3 IO * i i i i * k * i k * 0.0 I 4.54 i 0.0 I
5 3.086 Shear riscositj' oj'liquid rubidium 597 Table 3. NEMD simulation results for liquid rubidium at a temperature of K and density g ~ m - ~ Least-squares, fitting of the data to equation (9) gives ~= cp. The zero strain rate value of CP compares with the experimental value of CP and the Bansal and Bruns (1984) value of o 0. I85 i ~ i I4 k 0.01 I z i i i I72 i i 0. I i T i i increasing shear rate-that is, rubidium exhibits shear thinning, and so is non-newtonian (Bird et a1 1960). in the range of strain rates used in the NEMD simulations. This is typical of the NEMD technique. Because of the small system sizes involved in these simulations, the strain rates must be quite large in order for there to be a discernible velocity profile across the simulation sample (so that a value of P,, can be determined with small fluctuations). For liquid rubidium, a reduced strain rate of y* = 1.0 corresponds to a strain rate in real units of y=4.482 x loii s-' which exceeds considerably strain rates currently accessible to experiment. (Note, however, that for many systems-particularly colloidal suspensions where the particles have a very high molecular mass-a reduced strain rate of y* = 1.O corresponds to strain rates in real units well within present laboratory capabilities.) Clearly then, the viscosity of interest which must be compared with experiment is the Newtonian viscosity obtained as the zero strain rate extrapolation of the simulation results for 7. Theoretical considerations (Kawasaki and Gunton 1973, Yamada and Kawasaki 1975, Ernst et a/ 1978) and extensive "ID simulation results (Evans 1979, 1981, Evans and Hanley 1980) indicate that to leading order v is related to the strain rate by where vo and vi are independent of 7. Consequently, we illustrate the simulation results for From these figures, it is clear the shear viscosity in figures 1 to 3 by plotting v against e. that the simulation results correlate very well with equation (9). To calculate vo and vl, the simulation results are least-squares fitted to equation (9) with the results given in tables 2 to 4. The zero-strain rate results are then reported in table 1. As an aside, the error estimates given in tables 2 to 4 and illustrated as error bars in figures 1 to 3 deserve some explanation. In the appendix, the individual simulation runs which were used to obtain the means and error estimates given in tables 2 to 4 are Table 4. hemd simulation results for liquid rubidium at a temperature of K and density g cm-3. Least-squares fitting of the data to equation (9) gives q= y*'" cp. The zero strain rate value of 38 CP compares with the experimental value of 55 CP and the Bansal and Bruns (1984) value of I.o 28 f ko i i l i i i f I 0. I34 i i i i i i 0.01
6 598 P T Cummings and G P Morriss Figure 1. The shear viscosity q of liquid rubidium at K and g cm-3 as a function of the square root of the reduced strain rate y*. The circles and the error bars represent the simulation results and the estimated error respectively (see text for explanation). The straight line is the least-squares fit of the experimental data to equation (9) and is given in table 2. described in detail. At most state points (defined by the density, temperature and strain rate), several simulation runs were performed, each of which yielded a value (or measurement ) for q, the pressure P and the configurational energy Uconf. The values reported in tables 2 to 4 are then obtained as the statistical mean of the measurements (weighted according to run length) and the errors represent one standard deviation. For Figure 2. The shear viscosity q of liquid rubidium at K and g cm-3 as a function of the square root of the reduced strain rate y*. The circles and the error bars represent the simulation results and the estimated error respectively (see text for explanation). The straight line is the least-squares fit of the experimental data to equation (9) and is given in table 3.
7 Shear viscosity of liquid rubidium o 112 Figure 3. The shear viscosity q of liquid rubidium at K and g cm-3 as a function of the square root of the reduced strain rate y*. The circles and the error bars represent the simulation results and the estimated error respectively (see text for explanation). The straight line is the least-squares fit of the experimental data to equation (9) and is given in table 4. those state points where only one long simulation run was performed, the errors have been conjectured, the conjecture being based in part on the block averages of 2000 timesteps performed throughout every simulation run. Thus, the reported values and error estimates in tables 2 to 4 represent block size weighted block averages and standard deviations. Recently, Straatsma et a1 (1986) have examined the formal correctness of estimating errors based on standard deviations obtained from block averaging and they conclude that the quality of the error estimate depends on the choice of the block size, which must be neither too small (as this increases correlation between measurements ) nor too large (as this reduces the number of measurements ). On the basis of their conclusions, it is possible that the errors reported in this paper may be artificially low. In view of this, we have included details of the individual simulation runs in the appendix for the reader to draw his or her own conclusions about the quality of our results. From the results in table 1, we conclude that the Price potential is quite accurate in predicting the shear viscosity of liquid rubidium at the three state points shown. It is interesting to note that at the high-density state point (which is very near the triple point), the Price potential prediction overestimates the experimental result while it underestimates it at the medium and lower density. Since the temperature is also varying at each state point, it is difficult to draw a firm conclusion about where the Price potential might yield exact viscosities, but it is clear that this region of exactitude of the Price potential is bounded by the temperatures and densities of the medium- and high-density state points simulated in this study. That the Price potential is most accurate near the triple point and less accurate at expanded (lower density) states is consistent with results obtained by Rahman (1974b) and Mountain (1978) for the structure factor of liquid rubidium. Using equilibrium molecular dynamics, Rahman found that the Price potential predicted the structure factor and inelastic neutron scattering data for rubidium at the triple point very well while Mountain found using Monte Carlo simulation that at expanded states the Price potential yielded less accurate predictions for the structure factor. The source of the inaccuracy of the Price
8 600 P T Crrtntiiitigs and G P Morriss potential at expanded states is not clear and will undoubtedly by the subject of further investigation. It is known (Kawasaki and Gunton 1973, Yamada and Kawasaki 1975, Ernst et al 1978; Evans 1979, 1981, Evans and Hanley 1980) that the pressure P depends on the strain rate according to the equation P=Po +Ply? Figure 4 shows the results obtained by plotting P against y*3'2 and clearly equation (10) is satisfied by the simulation results. (Strictly speaking, the results presented in figure 4 correspond to the trace of the pressure tensor P given in equation (3). The full expression for the pressure for liquid metals involves terms which take into account the electronic contribution and the density dependence of the Price potential.) According to Evans (1 983), the configurational part of the internal energy depends on the strain rate according to uconf = L'O + U] y3 2. The results obtained in this simulation are consistent with this dependence and can be fitted to equation (I 1) as indicated in figure 5. The NEMD results at the highest density state point, given in table 2, warrant further comment. Of the three state points simulated, the results at this state point show the most (10) (1 1) 'f T* 312 Figure 4. The dimensionless trace of the pressure tensor-see equation (3)-as a function of the cube of the square root of the strain rate Y* for liquid rubidium at the three state points T= K, p= g cm-3 (labelled I), T= K, p= g CITI-~ (labelled 11) and T= K.p= g ~ m ((labelled - ~ 111).
9 Shear riscositj- of liquid rubidium 60 1 deviation from the strain rate dependence for viscosity given in equations (9). In fact, plotting v against y* on a log-log plot suggests that these results fit a strain rate dependence of the form V=Vo + V I f (12) where a is slightly less than i. It has been observed that a failure to satisfy the strain rate dependence (9) for the shear viscosity over the whole range of simulated strain rates may indicate that there are several distinct regions where (9) is valid with different values of v0 and ql for each region corresponding to Bingham plastic behaviour (Evans 1982). This behaviour would indicate that at the density and temperature in question the simulated fluid would, at zero strain rate, be a crystalline solid rather than a liquid and that the imposition of a non-zero strain rate induces shear melting. However, there is no evidence in the present simulations nor in that of Rahman (1947b) that at zero strain rate the Price potential is predicting a solid rather than a liquid state. On the basis of the available information, we therefore conclude that the deviation from equation (9) is not due to shear melting effects. In view of the discussion above, we suspect that the errors in the results tabulated in table 2 represent underestimates and that equation (9) is satisfied to within the actual error of the simulation. To confirm this, we are in the process of performing further simulations at this state point using more values of y and evaluating more carefully the dependence of the results on system size. (To date, our investigation of system size dependence has been limited to relatively short calculations at this state point using 256 molecules for y* equal to 1.0 and which differ from the results reported in table 1 by less than 3%, suggesting that system size dependence of the results in table 2 is minimal.) 3. Conclusions We have demonstrated the accuracy of the Price potential in calculating the shear viscosity of liquid rubidium. The NEMD results presented in this paper represent errors of only 3-11% when compared with the corresponding experimental data. These results contrast
10 602 P T Cummings and G P Morriss with those obtained by Bansal and Bruns (1984) who found errors of 32-45% using a somewhat different NEMD simulation technique. Moreover, as is clear from table 1, Bansal and Bruns underestimate the experimental results where the results presented in this paper overestimate and tiice iwsa. Thus there is a strong discrepancy between the results presented here and those of Bansal and Bruns (1984) which must be attributed to the simulation methodology since the same pair potentials are used in both studies. (Since Bansal and Bruns also used the Price potential and the same size system (108 molecules), apart from the difference in simulation methodology the only other possible source of the discrepancy is the cut-off (the point beyond which the potential is assumed to be zero). In this paper, as noted in 0 2 a cut-off of 3a was used while Bansal and Bruns utilised a cut-off equal to half the box size which for the three densities g, g and g cm-3 correspond to cut-offs of 2.47, 2.57 and 2.67~ respectively. To probe the importance of the cut-off. we performed several short simulation runs with a cut-off of 2.5~ at the intermediate density state point. The results obtained for the shear viscosity lie within the error estimates given in table 3. This eliminates the difference in cut-off as a source of the discrepancy between the present results and those of Bansal and Bruns.) Given the success of the NEMD technique used in this paper on many other systems (both simple and molecular fluids), we are confident of the accuracy of the results presented in this paper. Accepting our results and considering the agreement of the simulation results presented in this paper with experimental data suggests that the Price potential is capable of predicting shear viscosities of liquid rubidium with quantitative accuracy. In future research, the accuracy of the Price potential for other transport and thermodynamic properties of liquid rubidium, as well as other liquid metals, will be examined. Acknowledgments PTC acknowledges the Camille and Henry Dreyfus Foundation for their support of this research through the award of a Grant for Newly Appointed Faculty in the Chemical Sciences. 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. Appendix Details of the individual simulation runs are provided in table Al. At each density and temperature, the simulations began with strain rate y* = 1.O and the molecules initially in a face-centred cubic crystalline configuration. For smaller values of the strain rate, a melted configuration obtained from an earlier simulation run was used as the initial configuration. At each value of the strain rate, a run of between and timesteps was performed and discarded in order to reach a steady state. The discarded runs are not recorded in table A 1.
11 Shear viscosity of liquid rubidium 603 Table A 1. Individual non-equilibrium molecular dynamics simulation runs performed to obtain the results given in tables 1 to 4. Timesteps T(K) p(gcm- ) y II (CP) Uconf/N& Pa3/& l o I.o I.O 1.o I.o 1.o o 1.o 1.o I I References Arlinghaus R T and Cummings P T I987 J. Phys. F: Met. Phys Bansal R and Bruns W 1984 J. Chern. Phys Bird R 8, Stewart W E and Lightfoot E N 1960 Transport Phenomena (New York: Wiley)
12 ~ 1984a 604 P T Curnmings arid G P Morriss Ciccotti G and Jacucci G 1975 Phw. Reo. Lett Ciccotti G, Jacucci G and McDonald I R 1976 Phj,s. Ret,. A J. Stat. Php Ernst M H, Cichocki B. Dorfman J R, Sharma J and van Beijeren H 1978 J. Stat. Phys Evans D J 1979 Ph.vs. Lett. A Phys.Rev.A Phys. Rev. A J. Chem. Phys Molecular Dj,namics Simulation ofstatistica1 Mechanical Systems presented at the International School of Physics Enrico Fermi 1985 to be published Evans D J and Hanley H J M 1980 Physica A Evans D J and Hoover W G 1986 Ann. Rer. Fluid Mech Evans D J and Morriss G P 1983 Chem. Phjs Compur. Phj,s. Rep. I b Php. Rec. A Harrison W A 1966 Pseudoporentials in the Theory of Metals (New York: Benjamin) Holian B L and Evans D J 1983 J. Chem. Phys Hoover W G and Ashurst W T 1975 Theoretical Chemistry, Advances and Perspectives Vol. I, ed. H Eying and D Henderson Hoover W G, Ladd A J C and Moran B 1982 Phrs. Rev. Lett Kawasaki K and Gunton J D 1973 Phys. Rev. A Landolt-Bornstein 1969 Zahltnwerte und Funktionem, 5. Tiel (Berlin: Springer) Lees A W and Edwards S F 1972 J. Phys. C: Solid State Php Morriss G P and Evans D E 1985 Mol. Phrs Mountain R D 1978 J. Phj,s. F:.Mer. Ph,rs Paolini G V, Ciccotti G and Massobrio C 1986 Phys. Rev. A to appear Price D L 1971 Phja. Rev. A Price D L, Singwi K S and Tosi M P 1970 Phys. Rev. B Rahman A 1974a Phys. Rev. Left b Ph.vs. Reis. A Singer K. Singer J V L and Fincham D 1980 Mol. Phys Straatsma T P, Berendsen H J C and Stam A J 1986 Mol. Php Yamada T and Kawasaki K 1975 Prog. Theor. Phys
Shear 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 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 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 informationGibbs ensemble simulation of phase equilibrium in the hard core two-yukawa fluid model for the Lennard-Jones fluid
MOLECULAR PHYSICS, 1989, VOL. 68, No. 3, 629-635 Gibbs ensemble simulation of phase equilibrium in the hard core two-yukawa fluid model for the Lennard-Jones fluid by E. N. RUDISILL and P. T. CUMMINGS
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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationComparison of the DSMC Method with an Exact Solution of the Boltzmann Equation
Comparison of the DSMC Method with an Exact Solution of the Boltzmann Equation J M Montanero and A Santos Departamento de Fisica, Universidad de Extremadura, Spain Abstract Results obtained from the Direct
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 informationThermal Conductivity of Complex Plasmas Using Novel Evan-Gillan Approach
Commun. Theor. Phys. 69 (218) 74 71 Vol. 69, No. 6, June 1, 218 Thermal Conductivity of Complex Plasmas Using Novel Evan-Gillan Approach Aamir Shahzad, 1,2, Syed Irfan Haider, 1 Muhammad Kashif, 1 Muhammad
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 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 informationTemperature and Pressure Controls
Ensembles Temperature and Pressure Controls 1. (E, V, N) microcanonical (constant energy) 2. (T, V, N) canonical, constant volume 3. (T, P N) constant pressure 4. (T, V, µ) grand canonical #2, 3 or 4 are
More informationTemperature and Pressure Controls
Ensembles Temperature and Pressure Controls 1. (E, V, N) microcanonical (constant energy) 2. (T, V, N) canonical, constant volume 3. (T, P N) constant pressure 4. (T, V, µ) grand canonical #2, 3 or 4 are
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 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 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 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 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 informationarxiv: v1 [cond-mat.stat-mech] 7 Mar 2019
Langevin thermostat for robust configurational and kinetic sampling Oded Farago, Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB EW, United Kingdom Department of Biomedical
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 informationBrownian Motion and Langevin Equations
1 Brownian Motion and Langevin Equations 1.1 Langevin Equation and the Fluctuation- Dissipation Theorem The theory of Brownian motion is perhaps the simplest approximate way to treat the dynamics of nonequilibrium
More informationThermodynamic Properties
Thermodynamic Properties (TP) Thermodynamic Properties Define and articulate some of the critical language and concepts of Thermodynamics Distinguish between the universe, system, surroundings, and boundary
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 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 informationON ALGORITHMS FOR BROWNIAN DYNAMICS COMPUTER SIMULATIONS
COMPUTATIONAL METHODS IN SCIENCE AND TECHNOLOGY 4,35-42 (1998) ON ALGORITHMS FOR BROWNIAN DYNAMICS COMPUTER SIMULATIONS ARKADIUSZ C. BRAŃKA Institute of Molecular Physics, Polish Academy of Sciences, Smoluchowskiego
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 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 informationvs. Chapter 4: Standard Flows Chapter 4: Standard Flows for Rheology shear elongation 2/1/2016 CM4650 Lectures 1-3: Intro, Mathematical Review
CM465 Lectures -3: Intro, Mathematical //6 Chapter 4: Standard Flows CM465 Polymer Rheology Michigan Tech Newtonian fluids: vs. non-newtonian fluids: How can we investigate non-newtonian behavior? CONSTANT
More informationSupplemental Material for Temperature-sensitive colloidal phase behavior induced by critical Casimir forces
Supplemental Material for Temperature-sensitive colloidal phase behavior induced by critical Casimir forces Minh Triet Dang, 1 Ana Vila Verde, 2 Van Duc Nguyen, 1 Peter G. Bolhuis, 3 and Peter Schall 1
More informationNonequilibrium temperature measurement in a thermal conduction process
PHYSICAL REVIEW E 95, 33 (7) Nonequilibrium temperature measurement in a thermal conduction process Puneet Kumar Patra * and Romesh C. Batra Department of Biomedical Engineering and Mechanics, M/C 9, Virginia
More informationPHONON DISPERSION CURVES OF LIQUID METALS
Journal of Non-Oxide Glasses Vol., No, 010, p. 91-96 PHONON DISPERSION CURVES O LIQUID METALS ADITYA M. VORA * Humanities and Social Science Department, S. T. B. S. College of Diploma Engineering, Opp.
More informationThe electric multipole expansion for a magic cube
INSTITUTE OF PHYSICS PUBLISHING Eur. J. Phys. 26 (2005) 809 813 EUROPEAN JOURNAL OF PHYSICS doi:10.1088/0143-0807/26/5/013 The electric multipole expansion for a magic cube Adam Rogers and Peter Loly Department
More informationEntropy production fluctuation theorem and the nonequilibrium work relation for free energy differences
PHYSICAL REVIEW E VOLUME 60, NUMBER 3 SEPTEMBER 1999 Entropy production fluctuation theorem and the nonequilibrium work relation for free energy differences Gavin E. Crooks* Department of Chemistry, University
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 informationPower Law of Molecular Weight of the Nucleation Rate of Folded Chain Crystals of Polyethylene
Macromolecules 2002, 35, 6985-6991 6985 Power Law of Molecular Weight of the Nucleation Rate of Folded Chain Crystals of Polyethylene Swapan K. Ghosh, Masamichi Hikosaka,*, Akihiko Toda, Shinichi Yamazaki,
More informationTHE ELECTRICAL AND THERMAL CONDUCTIVITIES OF MONOVALENT METALS. By P. G. KLEMENS* [Manuscript received November 5, 1953] Summary
THE ELECTRICAL AND THERMAL CONDUCTIVITIES OF MONOVALENT METALS By P. G. KLEMENS* [Manuscript received November 5, 1953] Summary The experimentally determined values of the high and low temperature electrical
More informationStructural properties of low-density liquid alkali metals
PRAMANA c Indian Academy of Sciences Vol. 65, No. 6 journal of December 2005 physics pp. 1085 1096 Structural properties of low-density liquid alkali metals A AKANDE 1, G A ADEBAYO 1,2 and O AKINLADE 2
More informationTopic2540 Newton-Laplace Equation
Topic540 Newton-Laplace Equation The Newton-Laplace Equation is the starting point for the determination of isentropic compressibilities of solutions [,] using the speed of sound u and density ρ; equation
More informationTake the measurement of a person's height as an example. Assuming that her height has been determined to be 5' 8", how accurate is our result?
Error Analysis Introduction The knowledge we have of the physical world is obtained by doing experiments and making measurements. It is important to understand how to express such data and how to analyze
More informationA tight-binding molecular dynamics study of phonon anharmonic effects in diamond and graphite
J. Phys.: Condens. Matter 9 (1997) 7071 7080. Printed in the UK PII: S0953-8984(97)83513-8 A tight-binding molecular dynamics study of phonon anharmonic effects in diamond and graphite G Kopidakis, C Z
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 informationStructural properties of low density liquid alkali metals
PRAMANA c Indian Academy of Sciences Vol. xx, No. x journal of xxxxxx 2005 physics pp. 1 12 Structural properties of low density liquid alkali metals A AKANDE 1, G A ADEBAYO 1,2 and O AKINLADE 2 1 The
More information150A Review Session 2/13/2014 Fluid Statics. Pressure acts in all directions, normal to the surrounding surfaces
Fluid Statics Pressure acts in all directions, normal to the surrounding surfaces or Whenever a pressure difference is the driving force, use gauge pressure o Bernoulli equation o Momentum balance with
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 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 informationThermodynamic behaviour of mixtures containing CO 2. A molecular simulation study
Thermodynamic behaviour of mixtures containing. A molecular simulation study V. Lachet, C. Nieto-Draghi, B. Creton (IFPEN) Å. Ervik, G. Skaugen, Ø. Wilhelmsen, M. Hammer (SINTEF) Introduction quality issues
More informationDetermining the shear viscosity of model liquids from molecular dynamics simulations Hess, Berk
University of Groningen Determining the shear viscosity of model liquids from molecular dynamics simulations Hess, Berk Published in: Journal of Chemical Physics DOI: 1.163/1.1421362 IMPORA OE: You are
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 informationTHE NATURE OF THERMODYNAMIC ENTROPY. 1 Introduction. James A. Putnam. 1.1 New Definitions for Mass and Force. Author of
THE NATURE OF THERMODYNAMIC ENTROPY James A. Putnam Author of http://newphysicstheory.com james@newphysicstheory.com Thermodynamic entropy is a pathway for transitioning from the mechanical world of fundamental
More informationMeasuring rheological properties using a slotted plate device
Korea-Australia Rheology Journal Vol. 19, No. 2, August 2007 pp. 75-80 Measuring rheological properties using a slotted plate device Daniel De Kee 1, Young Dae Kim* and Q. Dzuy Nguyen 2 Faculty of Applied
More informationDSMC Simulation of Binary Rarefied Gas Flows between Parallel Plates and Comparison to Other Methods
Simulation of Binary Rarefied Gas Flows between Parallel Plates and Comparison to Other Methods L. Szalmas Department of Mechanical Engineering, University of Thessaly, Pedion Areos, Volos 38334, Greece
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 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 informationA Mixed Finite Element Formulation for Solving Phase Change Problems with Convection
A Mixed Finite Element Formulation for Solving Phase Change Problems with Convection Youssef Belhamadia 1, Abdoulaye S. Kane 2, and André Fortin 3 1 University of Alberta, Campus Saint-Jean and Department
More informationSYMMETRY BREAKING PHENOMENA OF PURELY VISCOUS SHEAR-THINNING FLUID FLOW IN A LOCALLY CONSTRICTED CHANNEL
ISSN 1726-4529 Int j simul model 7 (2008) 4, 186-197 Original scientific paper SYMMETRY BREAKING PHENOMENA OF PURELY VISCOUS SHEAR-THINNING FLUID FLOW IN A LOCALLY CONSTRICTED CHANNEL Ternik, P. University
More informationJournal of Computational Physics
Journal of Computational Physics 231 (2012) 5306 5316 Contents lists available at SciVerse ScienceDirect Journal of Computational Physics journal homepage: www.elsevier.com/locate/jcp Nonequilibrium molecular
More informationPressure Dependent Study of the Solid-Solid Phase Change in 38-Atom Lennard-Jones Cluster
University of Rhode Island DigitalCommons@URI Chemistry Faculty Publications Chemistry 2005 Pressure Dependent Study of the Solid-Solid Phase Change in 38-Atom Lennard-Jones Cluster Dubravko Sabo University
More informationExercise: concepts from chapter 10
Reading:, Ch 10 1) The flow of magma with a viscosity as great as 10 10 Pa s, let alone that of rock with a viscosity of 10 20 Pa s, is difficult to comprehend because our common eperience is with s like
More informationAn explicit expression for finite-size corrections to the chemical potential
J. Phys.: Condens. Matter 1 (1989) 8659-8665. Printed in the UK An explicit expression for finite-size corrections to the chemical potential B Smitt and D Frenkelt t Koninklijke/Shell-Laboratorium, Amsterdam
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 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 informationPharmaceutics I صيدالنيات 1. Unit 6
Pharmaceutics I صيدالنيات 1 Unit 6 1 Rheology of suspensions Rheology, the study of flow, addresses the viscosity characteristics of powders, fluids, and semisolids. Materials are divided into two general
More informationNotes 4: Differential Form of the Conservation Equations
Low Speed Aerodynamics Notes 4: Differential Form of the Conservation Equations Deriving Conservation Equations From the Laws of Physics Physical Laws Fluids, being matter, must obey the laws of Physics.
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 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 informationMolecular Simulation Background
Molecular Simulation Background Why Simulation? 1. Predicting properties of (new) materials 2. Understanding phenomena on a molecular scale 3. Simulating known phenomena? Example: computing the melting
More informationNumerical study of the steady state fluctuation relations far from equilibrium
Numerical study of the steady state fluctuation relations far from equilibrium Stephen R. Williams, Debra J. Searles, and Denis J. Evans Citation: The Journal of Chemical Physics 124, 194102 (2006); doi:
More informationTime-Dependent Statistical Mechanics 1. Introduction
Time-Dependent Statistical Mechanics 1. Introduction c Hans C. Andersen Announcements September 24, 2009 Lecture 1 9/22/09 1 Topics of concern in the course We shall be concerned with the time dependent
More informationAvailable online at ScienceDirect. Procedia Engineering 90 (2014 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 9 (214 ) 599 64 1th International Conference on Mechanical Engineering, ICME 213 Validation criteria for DNS of turbulent heat
More informationAn Introduction to Two Phase Molecular Dynamics Simulation
An Introduction to Two Phase Molecular Dynamics Simulation David Keffer Department of Materials Science & Engineering University of Tennessee, Knoxville date begun: April 19, 2016 date last updated: April
More informationDynamic properties of Lennard-Jones fluids and liquid metals
PHYSICAL REVIEW E VOLUME 60, NUMBER 1 JULY 1999 Dynamic properties of Lennard-Jones fluids and liquid metals M. Canales Departament de Física i Enginyeria Nuclear, Universitat Politècnica de Catalunya,
More informationDepartment Mechanical Engineering IIT Kanpur. Instructors Dr. K. Muralidhar Dr. P.K. Panigrahi
NPTEL Online - IIT Bombay Course Name Optical Measurement Techniques in Thermal Science Department Mechanical Engineering IIT Kanpur Instructors Dr. K. Muralidhar Dr. P.K. Panigrahi Course Description
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 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 informationFUNDAMENTAL STUDY OF BINGHAM FLUID BY MEANS OF DAM-BREAK FLOW MODEL
Annual Journal of Hydraulic Engineering, JSCE, Vol.54, 2010, February FUNDAMENTAL STUDY OF BINGHAM FLUID BY MEANS OF DAM-BREAK FLOW MODEL How Tion PUAY1 and Takashi HOSODA2 1 Member of JSCE, Phd Student,
More informationAA210A Fundamentals of Compressible Flow. Chapter 1 - Introduction to fluid flow
AA210A Fundamentals of Compressible Flow Chapter 1 - Introduction to fluid flow 1 1.2 Conservation of mass Mass flux in the x-direction [ ρu ] = M L 3 L T = M L 2 T Momentum per unit volume Mass per unit
More informationA comparative study between dissipative particle dynamics and molecular dynamics for simple- and complex-geometry flows
THE JOURNAL OF CHEMICAL PHYSICS 123, 104107 2005 A comparative study between dissipative particle dynamics and molecular dynamics for simple- and complex-geometry flows Eric E. Keaveny, Igor V. Pivkin,
More informationWhat is Classical Molecular Dynamics?
What is Classical Molecular Dynamics? Simulation of explicit particles (atoms, ions,... ) Particles interact via relatively simple analytical potential functions Newton s equations of motion are integrated
More informationComputer Simulation of Shock Waves in Condensed Matter. Matthew R. Farrow 2 November 2007
Computer Simulation of Shock Waves in Condensed Matter Matthew R. Farrow 2 November 2007 Outline of talk Shock wave theory Results Conclusion Computer simulation of shock waves Shock Wave Theory Shock
More information