Numerical simulation of Bloch equations for dynamic magnetic resonance imaging. (English) Zbl 1378.78033

Summary: Magnetic resonance imaging (MRI) is a widely applied non-invasive imaging modality based on non-ionizing radiation which gives excellent images and soft tissue contrast of living tissues. We consider the modified Bloch problem as a model of MRI for flowing spins in an incompressible flow field. After establishing the well-posedness of the corresponding evolution problem, we analyze its spatial semi-discretization using discontinuous Galerkin methods. The high frequency time evolution requires a proper explicit and adaptive temporal discretization. The applicability of the approach is shown for basic examples.


78M10 Finite element, Galerkin and related methods applied to problems in optics and electromagnetic theory
78A25 Electromagnetic theory (general)
65N40 Method of lines for boundary value problems involving PDEs
65M60 Finite element, Rayleigh-Ritz and Galerkin methods for initial value and initial-boundary value problems involving PDEs
65L06 Multistep, Runge-Kutta and extrapolation methods for ordinary differential equations
92C55 Biomedical imaging and signal processing
78M20 Finite difference methods applied to problems in optics and electromagnetic theory


Full Text: DOI arXiv


[1] Benoit-Cattin, H.; Collewet, G.; Belaroussi, B.; Saint-Jalmes, H.; Odet, C., The SIMRI project: a versatile and interactive MRI simulator, J. Magn. Res., 173, 1, 97-115, (2005)
[2] Bittoun, J.; Taquin, J.; Sauzade, M., A computer algorithm for the simulation of any nuclear magnetic resonance (NMR) imaging method, Magn. Reson. Imaging, 2, 2, 113-120, (1984)
[3] Bloch, F., Nuclear induction, Phys. Rev., 70, 7-8, 460-474, (1946)
[4] Carr, H. Y., Steady-state free precession in nuclear magnetic resonance, Phys. Rev., 112, 5, 1693-1701, (1958)
[5] Deuflhard, P.; Bornemann, F., Scientific computing with ordinary differential equations, vol. 42, (2012), Springer Science & Business Media
[6] Di Pietro, D. A.; Ern, A., Mathematical aspects of discontinuous Galerkin methods, vol. 69, (2011), Springer Science & Business Media
[7] Ern, A.; Guermond, J.-L., Theory and practice of finite elements, vol. 159, (2004), Springer Science & Business Media
[8] J. Frahm, A. Haase, D. Matthaei, W. Hänicke, K.-D. Merboldt, Verfahren und Einrichtung zur schnellen Akquisition von Spinresonanzdaten für eine ortsaufgelöste Untersuchung eines Objekts, 1985.
[9] Haase, A.; Frahm, J.; Matthaei, D.; Hanicke, W.; Merboldt, K.-D., FLASH imaging. rapid NMR imaging using low flip-angle pulses, J. Magn. Res., 67, 2, 258-266, (1986)
[10] Hargreaves, B., Bloch equation simulator, (2004), published January 13
[11] Hazra, A., Numerical simulation of Bloch equations for dynamic magnetic resonance imaging, (2016), Institut für Numerische und Angewandte Mathematik, Ph.D. thesis · Zbl 1361.94003
[12] Hundsdorfer, W.; Verwer, J. G., Numerical solution of time-dependent advection-diffusion-reaction equations, (2003), Springer Science & Business Media · Zbl 1030.65100
[13] Joseph, A.; Kowallick, J. T.; Merboldt, K. D.; Voit, D.; Schaetz, S.; Zhang, S.; Sohns, J. M.; Lotz, J.; Frahm, J., Real-time flow MRI of the aorta at a resolution of 40 msec, J. Magn. Reson. Imaging, 40, 1, 206-213, (2014)
[14] Jou, L. D.; van Tyen, R.; Berger, S. A.; Saloner, D., Calculation of the magnetization distribution for fluid flow in curved vessels, Magn. Reson. Med., 35, 4, 577-584, (1996)
[15] Jurczuk, K.; Kretowski, M.; Bellanger, J. J.; Eliat, P. A.; Saint-Jalmes, H.; Bézy-Wendling, J., Computational modeling of MR flow imaging by the lattice Boltzmann method and Bloch equation, Magn. Reson. Imaging, 31, 7, 1163-1173, (2013)
[16] Lions, J. L., Perturbations singulières dans LES problèmes aux limites et en contrôle optimal, vol. 323, (2006), Springer · Zbl 0268.49001
[17] Lorthois, S.; Stroud-Rossman, J.; Berger, S.; Jou, L. D.; Saloner, D., Numerical simulation of magnetic resonance angiographies of an anatomically realistic stenotic carotid bifurcation, Ann. Biomed. Eng., 33, 3, 270-283, (2005)
[18] Nishimura, D. G., Principles of magnetic resonance imaging, (1996), Stanford University
[19] Roeloffs, V.; Voit, D.; Frahm, J., Spoiling without additional gradients: radial FLASH MRI with randomized radiofrequency phases, Magn. Reson. Med., 75, 5, 2094-2099, (2016)
[20] Shkarin, P.; Spencer, R. G.S., Time domain simulation of Fourier imaging by summation of isochromats, Int. J. Imaging Syst. Technol., 8, 419-426, (1997)
[21] Stejskal, E. O.; Tanner, J. E., Spin diffusion measurements: spin echoes in the presence of a time-dependant field gradient, J. Chem. Phys., 42, 1, 288-292, (1965)
[22] Stöcker, T.; Vahedipour, K.; Pflugfelder, D.; Shah, N. J., High-performance computing MRI simulations, Magn. Reson. Med., 64, 1, 186-193, (2010)
[23] Sun, H., Bloch equation simulator, (2012), published July 22
[24] Uecker, M.; Zhang, S.; Voit, D.; Karaus, A.; Merboldt, K.-D.; Frahm, J., Real-time MRI at a resolution of 20 ms, NMR Biomed., 23, 8, 986-994, (2010)
[25] Xanthis, C. G.; Venetis, I. E.; Chalkias, A. V.; Aletras, A. H., MRISIMUL: a GPU-based parallel approach to MRI simulations, IEEE Trans. Med. Imaging, 33, 3, 607-617, (2014)
[26] Yuan, C.; Gullberg, G. T.; Parker, D. L., The solution of Bloch equations for flowing spins during a selective pulse using a finite difference method, Med. Phys., 14, 6, 914-921, (1987)
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.