Source author record

Wm. G. Hoover

Wm. G. Hoover appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

20works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

20 published item(s)

preprint2020arXiv

Shockwaves and Local Hydrodynamics; Failure of the Navier-Stokes Equations

Shockwaves provide a useful and rewarding route to the nonequilibrium properties of simple fluids far from equilibrium. For simplicity, we study a strong shockwave in a dense two-dimensional fluid. Here, our study of nonlinear transport properties makes plain the connection between the observed local hydrodynamic variables (like the various gradients and fluxes) and the chosen recipes for defining (or "measuring") those variables. The range over which nonlocal hydrodynamic averages are computed turns out to be much more significant than are the other details of the averaging algorithms. The results show clearly the incompatibility of microscopic time-reversible cause-and-effect dynamics with macroscopic instantaneously-irreversible models like the Navier-Stokes equations.

preprint2015arXiv

What is liquid? Lyapunov instability reveals symmetry-breaking irreversibilities hidden within Hamilton's many-body equations of motion

Typical Hamiltonian liquids display exponential "Lyapunov instability", also called "sensitive dependence on initial conditions". Although Hamilton's equations are thoroughly time-reversible, the forward and backward Lyapunov instabilities can differ, qualitatively. In numerical work, the expected forward/backward pairing of Lyapunov exponents is also occasionally violated. To illustrate, we consider many-body inelastic collisions in two space dimensions. Two mirror-image colliding crystallites can either bounce, or not, giving rise to a single liquid drop, or to several smaller droplets, depending upon the initial kinetic energy and the interparticle forces. The difference between the forward and backward evolutionary instabilities of these problems can be correlated with dissipation and with the Second Law of Thermodynamics. Accordingly, these asymmetric stabilities of Hamilton's equations can provide an "Arrow of Time". We illustrate these facts for two small crystallites colliding so as to make a warm liquid. We use a specially-symmetrized form of Levesque and Verlet's bit-reversible Leapfrog integrator. We analyze trajectories over millions of collisions with several equally-spaced time reversals.

preprint2014arXiv

Time's Arrow for Shockwaves ; Bit-Reversible Lyapunov and "Covariant" Vectors ; Symmetry Breaking

Strong shockwaves generate entropy quickly and locally. The Newton-Hamilton equations of motion, which underly the dynamics, are perfectly time-reversible. How do they generate the irreversible shock entropy? What are the symptoms of this irreversibility? We investigate these questions using Levesque and Verlet's bit-reversible algorithm. In this way we can generate an entirely imaginary past consistent with the irreversibility observed in the present. We use Runge-Kutta integration to analyze the local Lyapunov instability of the forward and backward processes so as to identify those particles most intimately connected with the irreversibility described by the Second Law of Thermodynamics. Despite the perfect time symmetry of the particle trajectories, the fully-converged vectors associated with the largest Lyapunov exponents, forward and backward in time, are qualitatively different. The vectors display a time-symmetry breaking equivalent to Time's Arrow. That is, in autonomous Hamiltonian shockwaves the largest local Lyapunov exponents, forward and backward in time, are quite different.

preprint2013arXiv

Hamiltonian Thermostats Fail to Promote Heat Flow

Hamiltonian mechanics can be used to constrain temperature simultaneously with energy. We illustrate the interesting situations that develop when two different temperatures are imposed within a composite Hamiltonian system. The model systems we treat are "phi-4" chains, with quartic tethers and quadratic nearest-neighbor Hooke's-law interactions. This model is known to satisfy Fourier's law. Our prototypical problem sandwiches a Newtonian subsystem between hot and cold Hamiltonian reservoir regions. We have characterized four different Hamiltonian reservoir types. There is no tendency for any of these two-temperature Hamiltonian simulations to transfer heat from the hot to the cold degrees of freedom. Evidently steady heat flow simulations require energy sources and sinks, and are therefore incompatible with Hamiltonian mechanics.

preprint2013arXiv

Shockwave Compression and Joule-Thomson Expansion

Structurally-stable atomistic one-dimensional shockwaves have long been simulated by injecting fresh cool particles and extracting old hot particles at opposite ends of a simulation box. The resulting shock profiles demonstrate tensor temperature, with the longitudinal temperature exceeding the transverse, and Maxwell's delayed response, with stress lagging strainrate and heat flux lagging temperature gradient. Here this same geometry, supplemented by a short-ranged external "plug" field, is used to simulate steady Joule-Kelvin throttling flow of hot dense fluid through a porous plug, producing a dilute and cooler product fluid.

preprint2013arXiv

Time-Reversible Random Number Generators : Solution of Our Challenge by Federico Ricci-Tersenghi

Nearly all the evolution equations of physics are time-reversible, in the sense that a movie of the solution, played backwards, would obey exactly the same differential equations as the original forward solution. By way of contrast, stochastic approaches are typically not time-reversible, though they could be made so by the simple expedient of storing their underlying pseudorandom numbers in an array. Here we illustrate the notion of time-reversible random number generators. In Version 1 we offered a suitable reward for the first arXiv response furnishing a reversed version of an only slightly-more-complicated pseudorandom number generator. Here we include Professor Ricci-Tersenghi's prize-winning reversed version as described in his arXiv:1305.1805 contribution: "The Solution to the Challenge in `Time-Reversible Random Number Generators' by Wm. G. Hoover and Carol G. Hoover".

preprint2013arXiv

Time-Symmetry Breaking in Hamiltonian Mechanics

Hamiltonian trajectories are strictly time-reversible. Any time series of Hamiltonian coordinates {q} satisfying Hamilton's motion equations will likewise satisfy them when played "backwards", with the corresponding momenta changing signs : {+p} --> {-p}. Here we adopt Levesque and Verlet's precisely bit-reversible motion algorithm to ensure that the trajectory reversibility is exact, with the forward and backward sets of coordinates identical. Nevertheless, the associated instantaneous Lyapunov instability, or "sensitive dependence on initial conditions" of "chaotic" (or "Lyapunov unstable") bit-reversible coordinate trajectories can still exhibit an exponentially growing time-symmetry-breaking irreversibility. Surprisingly, the positive and negative exponents, as well as the forward and backward Lyapunov spectra, are usually not closely related, and so give four differing topological measures of "local" chaos. We have demonstrated this symmetry breaking for fluid shockwaves, for free expansions, and for chaotic molecular collisions. Here we illustrate and discuss this time-symmetry breaking for three statistical-mechanical systems, [1] a minimal (but still chaotic) one-body "cell model" with a four-dimensional phase space; [2] relatively small colliding crystallites, for which the whole Lyapunov spectrum is accessible; [3] a near-continuum inelastic collision of two larger 400-particle balls. In the last two of these pedagogical problems the two colliding bodies coalesce. The particles most prone to Lyapunov instability are dramatically different in the two time directions. Thus this Lyapunov-based symmetry breaking furnishes an interesting Arrow of Time.

preprint2013arXiv

Why Instantaneous Values of the "Covariant" Lyapunov Exponents Depend upon the Chosen State-Space Scale

We explore a simple example of a chaotic thermostated harmonic-oscillator system which exhibits qualitatively different local Lyapunov exponents for simple scale-model constant-volume transformations of its coordinate q and momentum p : { q,p } --> { (Q/s),(sP) } . The time-dependent thermostat variable zeta(t) is unchanged by such scaling. The original (q,p,zeta) motion and the scale-model (Q,P,zeta) version of the motion are physically identical. But both the local Gram-Schmidt Lyapunov exponents and the related local "covariant" exponents change with the change of scale. Thus this model furnishes a clearcut chaotic time-reversible example showing how and why both the local Lyapunov exponents and covariant exponents vary with the scale factor s .

preprint2012arXiv

Another Hamiltonian "Thermostat" - Comments on arXiv Contributions 1203.5968, 1204.4412, 1205.3478, and 1206.0188

Campisi, Zhan, Talkner, and Haenggi state, in promoting a new logarithmic computational thermostat [ arXiv 1203.5968 and 1204.4412 ], that (thermostated) Nose-Hoover mechanics is not Hamiltonian. First I point out that Dettmann clearly showed the Hamiltonian nature of Nose-Hoover mechanics. The trajectories {q(t)} generated by Dettmann's Hamiltonian are identical to those generated by Nose-Hoover mechanics. I also observe that when the (Hamiltonian) Campisi thermostat is applied to "nonequilibrium" heat transfer problems some very interesting, and somewhat paradoxical, phase portraits can result. See too Marc Meléndez' nice arXiv 1205.3478 and our joint work 1206.0188 .

preprint2012arXiv

Microscopic and Macroscopic Rayleigh-Benard Flows : Continuum and Particle Simulations, Turbulence, Fluctuations, Time Reversibility, and Lyapunov Instability

We discuss the irreversibility, nonlocality, and fluctuations, as well as the Lyapunov and hydrodynamic instabilities characterizing atomistic, smooth-particle, and finite-difference solutions of the two-dimensional Rayleigh-Bénard problem. To speed up the numerical analysis we control the time-dependence of the Rayleigh number so as to include many distinct flow morphologies in a single simulation. The relatively simple nature of these computational algorithms and the richness of the results they can yield make such studies and their interpretation particularly well suited to graduate-level research.

preprint2012arXiv

Steady Periodic Shear Flow is Stable in Two Space Dimensions . Nonequilibrium Molecular Dynamics vs Navier-Stokes-Fourier Stability Theory -- A Comment on two Arxiv Contributions

Dufty, Lee, Lutsko, Montanero, and Santos have carried out stability analyses of steady stationary shear flows. Their approach is based on the compressible and heat conducting Navier-Stokes-Fourier model. It predicts the unstable exponential growth of long-wavelength transverse perturbations for both two- and three-dimensional fluids. We point out that the patently-stable two-dimensional periodic shear flows studied earlier by Petravic, Posch, and ourselves contradict these predicted instabilities. The stable steady-state shear flows are based on nonequilibrium molecular dynamics with simple thermostats maintaining nonequilibrium stationary states in two space dimensions. The failure of the stability analyses remains unexplained.

preprint2011arXiv

Local Gram-Schmidt and Covariant Lyapunov Vectors and Exponents for Three Harmonic Oscillator Problems

We compare the Gram-Schmidt and covariant phase-space-basis-vector descriptions for three time-reversible harmonic oscillator problems, in two, three, and four phase-space dimensions respectively. The two-dimensional problem can be solved analytically. The three-dimensional and four-dimensional problems studied here are simultaneously chaotic, time-reversible, and dissipative. Our treatment is intended to be pedagogical, for publication in Communications in Nonlinear Science and Numerical Computation and for use in an updated version of our book on Time Reversibility, Computer Simulation, and Chaos. Comments are very welcome.

preprint2011arXiv

Maxwell and Cattaneo's Time-Delay Ideas Applied to Shockwaves and the Rayleigh-Benard Problem

We apply Maxwell and Cattaneo's relaxation approaches to the analysis of strong shockwaves in a two-dimensional viscous heat-conducting fluid. Good agreement results for reasonable values of Maxwell's relaxation times. Instability results if the viscous relaxation time is too large. These relaxation terms have negligible effects on slower-paced subsonic problems, as is shown here for two-roll and four-roll Rayleigh-Benard flows.

preprint2011arXiv

Nonequilibrium Fluctuations in a Gaussian Galton Board (or Periodic Lorentz Gas) Using Long Periodic Orbits

Predicting nonequilibrium fluctuations requires a knowledge of nonequilibrium distribution functions. Despite the distributions' fractal character some theoretical results, "Fluctuation Theorems", reminiscent of but distinct from, Gibbs' equilibrium statistical mechanics and the Central Limit Theorem, have been established away from equilibrium and applied to simple models. We summarize the simplest of these results for a Gaussian-thermostated Galton Board problem, a field-driven mass point moving through a periodic array of hard-disk scatterers. The billion-collision trillion-timestep data we analyze correspond to periodic orbits with up to 793,951,594 collisions and 447,064,397,614 timesteps.

preprint2010arXiv

Flexible Macroscopic Models for Dense-Fluid Shockwaves: Partitioning Heat and Work; Delaying Stress and Heat Flux; Two-Temperature Thermal Relaxation

Macroscopic models which distinguish the longitudinal and transverse temperatures can provide improved descriptions of the microscopic shock structures as revealed by molecular dynamics simulations. Additionally, we can include three relaxation times in the models, two based on Maxwell's viscoelasticity and its Cattaneo-equation analog for heat flow, and a third thermal, based on the Krook-Boltzmann equation. This approach can replicate the observed lags of stress (which lags behind the strain rate) and heat flux (which lags behind the temperature gradient), as well as the eventual equilibration of the two temperatures. For profile stability the time lags cannot be too large. By partitioning the longitudinal and transverse contributions of work and heat and including a tensor heat conductivity and bulk viscosity, all the qualitative microscopic features of strong simple-fluid shockwave structures can be reproduced.

preprint2010arXiv

Three Lectures: Nemd, Spam, and Shockwaves

We discuss three related subjects well suited to graduate research. The first, Nonequilibrium molecular dynamics or "NEMD", makes possible the simulation of atomistic systems driven by external fields, subject to dynamic constraints, and thermostated so as to yield stationary nonequilibrium states. The second subject, Smooth Particle Applied Mechanics or "SPAM", provides a particle method, resembling molecular dynamics, but designed to solve continuum problems. The numerical work is simplified because the SPAM particles obey ordinary, rather than partial, differential equations. The interpolation method used with SPAM is a powerful interpretive tool converting point particle variables to twice-differentiable field variables. This interpolation method is vital to the study and understanding of the third research topic we discuss, strong shockwaves in dense fluids. Such shockwaves exhibit stationary far-from-equilibrium states obtained with purely reversible Hamiltonian mechanics. The SPAM interpolation method, applied to this molecular dynamics problem, clearly demonstrates both the tensor character of kinetic temperature and the time-delayed response of stress and heat flux to the strain rate and temperature gradients. The dynamic Lyapunov instability of the shockwave problem can be analyzed in a variety of ways, both with and without symmetry in time. These three subjects suggest many topics suitable for graduate research in nonlinear nonequilibrium problems.

preprint2010arXiv

Well-Posed Two-Temperature Constitutive Equations for Stable Dense Fluid Shockwaves using Molecular Dynamics and Generalizations of Navier-Stokes-Fourier Continuum Mechanics

Guided by molecular dynamics simulations, we generalize the Navier-Stokes-Fourier constitutive equations and the continuum motion equations to include both transverse and longitudinal temperatures. To do so we partition the contributions of the heat transfer, the work done, and the heat flux vector between the longitudinal and transverse temperatures. With shockwave boundary conditions time-dependent solutions of these equations converge to give stationary shockwave profiles. The profiles include anisotropic temperature and can be fitted to molecular dynamics results, demonstrating the utility and simplicity of a two-temperature description of far-from-equilibrium states.

preprint2009arXiv

Nonlinear Stresses and Temperatures in Transient Adiabatic and Shear Flows via Nonequilibrium Molecular Dynamics -- Three Definitions of Temperature

We compare nonlinear stresses and temperatures for adiabatic shear flows, using up to 262,144 particles, with those from corresponding homogeneous and inhomogeneous flows. Two varieties of kinetic temperature tensors are compared to the configurational temperatures. This comparison leads to an improved form for the local and instantaneous smooth-particle averaged stream velocity and to a recognition of rotational contributions to the configurational temperature.

preprint2008arXiv

2D and 3D Dense-Fluid Shear Flows via Nonequilibrium Molecular Dynamics. Comparison of Time-and-Space-Averaged Tensor Temperature and Normal Stresses from Doll's, Sllod, and Boundary-Driven Shear Algorithms

Homogeneous shear flows (with constant strainrate du/dy) are generated with the Doll's and Sllod algorithms and compared to corresponding inhomogeneous boundary-driven flows. We use one-, two-, and three-dimensional smooth-particle weight functions for computing instantaneous spatial averages. The nonlinear stress differences are small, but significant, in both two and three space dimensions. In homogeneous systems the sign and magnitude of the shearplane stress difference, P(xx) - P(yy), depend on both the thermostat type and the chosen shearflow algorithm. The Doll's and Sllod algorithms predict opposite signs for this stress difference, with the Sllod approach definitely wrong, but somewhat closer to the (boundary-driven) truth. Neither of the homogeneous shear algorithms predicts the correct ordering of the kinetic temperatures, T(xx) > T(zz) > T(yy).

preprint2000arXiv

Entropy production and Lyapunov instability at the onset of turbulent convection

Computer simulations of a compressible fluid, convecting heat in two dimensions, suggest that, within a range of Rayleigh numbers, two distinctly different, but stable, time-dependent flow morphologies are possible. The simpler of the flows has two characteristic frequencies: the rotation frequency of the convecting rolls, and the vertical oscillation frequency of the rolls. Observables, such as the heat flux, have a simple-periodic (harmonic) time dependence. The more complex flow has at least one additional characteristic frequency -- the horizontal frequency of the cold, downward- and the warm, upward-flowing plumes. Observables of this latter flow have a broadband frequency distribution. The two flow morphologies, at the same Rayleigh number, have different rates of entropy production and different Lyapunov exponents. The simpler "harmonic" flow transports more heat (produces entropy at a greater rate), whereas the more complex "chaotic" flow has a larger maximum Lyapunov exponent (corresponding to a larger rate of phase-space information loss). A linear combination of these two rates is invariant for the two flow morphologies over the entire range of Rayleigh numbers for which the flows coexist, suggesting a relation between the two rates near the onset of convective turbulence.