Problems with discontinuous material boundaries have fields which suffer from jump discontinuities
This means the divergence of the Maxwell stress tensor is impulsive at material interfaces
Evaluating net force/torque using body forces requires both volumetric and surface integrals in these cases
This is a nontrivial excercise

In the meantime, the next step will be to add a special calculation method for cylindrical airgaps