Calculating Fx / Fy on surface without integration

This would involve adding a normal vector to your pressure, so you can use code from Obtain velocity normal to boundary - #2 by dokken and multiply by pressure. Simply change the space to be a VectorFunctionSpace over DG-1 (as the pressure in normal direction only makes sense wrt. a given cell).