As discussed with @bauer and @holzer, we should try applying the force in moment space. This should be equivalent to the Guo force model, but save significant FLOPs.
To support this effort, the automatic Chapman-Enskog analysis should be extended to support forces.