Spectral-diagonalization-based matrix exponential integration for efficient and stable solutions of full-Bloch equations in surface NMR

Surface nuclear magnetic resonance (SNMR) is a geophysical extension of nuclear magnetic resonance (NMR) that enables non-invasive mapping of subsurface hydrogeological properties by measuring the relaxation response of groundwater hydrogen nuclei. Accurately modeling the transient spin dynamics in...

Celý popis

Uloženo v:
Podrobná bibliografie
Vydáno v:Computers & geosciences Ročník 207; s. 106073
Hlavní autoři: Lin, Tingting, Wang, Qingyue, Jiang, Chuandong, Ren, Chunpeng, Wang, Yunzhi, Wang, Liang
Médium: Journal Article
Jazyk:angličtina
Vydáno: Elsevier Ltd 01.02.2026
Témata:
Bloch equations
Matrix exponential integration
SNMR
Spectral diagonalization
Spectral diagonalization
Bloch equations
Matrix exponential integration
SNMR
ISSN:0098-3004
On-line přístup:Získat plný text
Tagy: Přidat tag
Žádné tagy, Buďte první, kdo vytvoří štítek k tomuto záznamu!
Popis
Shrnutí:Surface nuclear magnetic resonance (SNMR) is a geophysical extension of nuclear magnetic resonance (NMR) that enables non-invasive mapping of subsurface hydrogeological properties by measuring the relaxation response of groundwater hydrogen nuclei. Accurately modeling the transient spin dynamics in SNMR requires solving the full-Bloch equations under Earth’s geomagnetic field, where magnetic field inhomogeneities, multicomponent relaxation, and nonlinear pulsed excitations introduce significant mathematical and computational challenges. We present a spectral-diagonalization-based matrix exponential integration (SD-MEI) algorithm for efficient and stable solutions of full-Bloch equations in SNMR. Conventional explicit numerical methods exhibit cumulative discretization errors and escalating computational costs due to step-size dependence and finite precision limitations. SD-MEI integrates spectral diagonalization with matrix exponential operations, replacing iterative computations with a single eigendecomposition of the system matrix. This approach achieves parameter-robust computational complexity while maintaining numerical stability across broad B1 field strengths (10−10 T to 10−5 T) and relaxation times (10 ms to 1000 ms). Validated for steady-state free precession (SSFP) dynamics in heterogeneous geomagnetic environments, the method enables high-accuracy modeling of transient magnetization evolution with large time steps. The framework advances SNMR efficient forward modeling and inversion while optimizing protocols by resolving critical limitations in existing numerical and analytical approaches. •Modeling of Magnetization Dynamics: Accurately describes magnetization vectors over broad excitation fields.•Control of Computational Complexity: Keeps computation stable, improving efficiency by 20–1000x.•Improved Numerical Stability: Solves convergence issues in stiff and highly nonlinear systems.•Temporal Flexibility: Computes magnetization at any time without affecting overall accuracy.
ISSN:0098-3004
DOI:10.1016/j.cageo.2025.106073