  We develop a practical discrete model of hysteresis based on nonlinear play
and generalized play, for use in first-order conservation laws with
applications to adsorption-desorption hysteresis models. The model is easy to
calibrate from sparse data, and offers rich secondary curves. We compare it
with discrete regularized Preisach models. We also prove well-posedness and
numerical stability of the class of hysteresis operators involving all those
types, describe implementation and present numerical examples using
experimental data.
