High-order Unfitted Discretizations for PDEs Coupled with Geometric Flow
Loading