Shock-wave compression and Joule-Thomson expansion

Structurally stable atomistic one-dimensional shock waves 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, Txx≠Tyy and Maxwell's delayed response, with stress lagging strain rate 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.