Efficient solvers for coupled Brown-Néel Fokker-Planck equations

In magnetic particle imaging (MPI), achieving efficient and accurate solutions to forward models is crucial for solving inverse problems. This work investigates the coupled Brown/Néel model, which leads to a convection-dominated Fokker-Planck equation defined over a higher-dimensional domain. To handle these challenges, we propose a joint angular-temporal discretization of this partial differential equation (PDE) combined with a reduced basis approach. Preliminary numerical results on simplified models demonstrate the efficacy of our approach, indicating its potential for broader applications in MPI.