  We present a mixed finite element solver for the linearized R13 equations of
non-equilibrium gas dynamics. The Python implementation builds upon the
software tools provided by the FEniCS computing platform. We describe a new
tensorial approach utilizing the extension capabilities of FEniCS's Unified
Form Language (UFL) to define required differential operators for tensors above
second degree. The presented solver serves as an example for implementing
tensorial variational formulations in FEniCS, for which the documentation and
literature seem to be very sparse. Using the software abstraction levels
provided by the UFL allows an almost one-to-one correspondence between the
underlying mathematics and the resulting source code. Test cases support the
correctness of the proposed method using validation with exact solutions. To
justify the usage of extended gas flow models, we discuss typical application
cases involving rarefaction effects. We provide the documented and validated
solver publicly.
