A mathematical model of the critical state based on averaged fluxon motion has been implemented to solve for the current and field distributions inside a high temperature superconducting hysteresis machine. The machine consists of a rotor made from a solid cylindrical single domain HTS placed in a perpendicular rotating field. The solution technique uses the finite difference approximation for a two-dimensional domain, discretized in cylindrical polar co-ordinates. The torque generated or equivalently the hysteresis loss in such a machine has been investigated using the model. It was found that to maximize the efficiency, the field needs to penetrate the rotor such that B-0/mu(0)J(c)R = 0.56, where B-0 is the applied field amplitude, J(c) is the critical current density and R is the rotor radius. This corresponds to a penetration that is 27% greater than that which reaches the centre of the rotor. An examination of the torque density distributions across the rotor reveal that for situations where the field is less than optimal, a significant increase in the performance can be achieved by removing an inner cylinder from the rotor.