International Journal on Magnetic Particle Imaging IJMPI
Vol. 11 No. 1 Suppl 1 (2025): Int J Mag Part Imag
https://doi.org/10.18416/IJMPI.2025.2503011

Proceedings Articles

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

Main Article Content

Manfred Faldum (Institute for Geometry and Applied Mathematics, RWTH Aachen University), Markus Bachmayr (Institute for Geometry and Applied Mathematics, RWTH-Aachen University), Volkmar Schulz (1) Institute of Imaging and Computer Vision, RWTH-Aachen University 2) Fraunhofer Institute for Digital Medicine Medicine MEVIS, Aachen), Franziska Schrank (1) Institute of Imaging and Computer Vision, RWTH-Aachen University 2) Institute for Experimental Molecular Imaging, RWTH-Aachen University)

Copyright (c) 2025 Manfred Faldum, Markus Bachmayr, Volkmar Schulz, Franziska Schrank

Creative Commons License

This work is licensed under a Creative Commons Attribution 4.0 International License.

Abstract

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.

Article Details

References

[1] T. Kluth, P. Szwargulski, and T. Knopp. Towards accurate modeling of the multidimensional magnetic particle imaging physics. New Journal of Physics, 21(10):103032, 2019, doi:10.1088/1367-2630/ab4938.
[2] J. Weizenecker. The fokker–planck equation for coupled brown–néel-rotation. Physics in Medicine & Biology, 63(3):035004, 2018, doi:10.1088/1361-6560/aaa186.
[3] T. Kluth. Mathematical models for magnetic particle imaging. Inverse Problems, 34(8):083001, 2018, doi:10.1088/1361-6420/aac535.
[4] H. Albers, T. Kluth, and T. Knopp. A simulation framework for particle magnetization dynamics of large ensembles of single domain particles: Numerical treatment of brown/néel dynamics and parameter identification problems in magnetic particle imaging. ArXiv, abs/2010.07772, 2020.
[5] H. Egger and J. Schöberl. A hybrid mixed discontinuous Galerkin finite-element method for convection–diffusion problems. IMA Journal of Numerical Analysis, 30(4):1206–1234, 2009, doi:10.1093/imanum/drn083.
[6] M. Griebel. The combination technique for the sparse grid solution of pde’s on multiprocessor machines. Parallel Processing Letters, 02(01):61–70, 1992, doi:10.1142/S0129626492000180.
[7] M. Bachmayr. Low-rank tensor methods for partial differential equations. Acta Numerica, 32:1–121, 2023, doi:10.1017/S0962492922000125.