double getCriticalPressure(TimeStep step) { return particleDensity * getCriticalVelocitySquared(step); }