Search: onr:"swepub:oai:DiVA.org:kth-234286" >
A partition of unit...
A partition of unity finite element method for computational diffusion MRI
-
- Nguyen, Van Dang, 1985- (author)
- KTH,Beräkningsvetenskap och beräkningsteknik (CST)
-
- Jansson, Johan (author)
- KTH,Beräkningsvetenskap och beräkningsteknik (CST)
-
- Hoffman, Johan, 1974- (author)
- KTH,Beräkningsvetenskap och beräkningsteknik (CST)
-
show more...
-
- Li, Jing-Rebecca (author)
- INRIA Saclay-Equipe DEFI, CMAP, Ecole Polytechnique Route de Saclay, 91128, Palaiseau Cedex, France
-
show less...
-
(creator_code:org_t)
- Elsevier, 2018
- 2018
- English.
-
In: Journal of Computational Physics. - : Elsevier. - 0021-9991 .- 1090-2716. ; 375, s. 271-290
- Related links:
-
https://bird.bcamath...
-
show more...
-
https://urn.kb.se/re...
-
https://doi.org/10.1...
-
show less...
Abstract
Subject headings
Close
- The Bloch–Torrey equation describes the evolution of the spin (usually water proton) magnetization under the influence of applied magnetic field gradients and is commonly used in numerical simulations for diffusion MRI and NMR. Microscopic heterogeneity inside the imaging voxel is modeled by interfaces inside the simulation domain, where a discontinuity in the magnetization across the interfaces is produced via a permeability coefficient on the interfaces. To avoid having to simulate on a computational domain that is the size of an entire imaging voxel, which is often much larger than the scale of the microscopic heterogeneity as well as the mean spin diffusion displacement, smaller representative volumes of the imaging medium can be used as the simulation domain. In this case, the exterior boundaries of a representative volume either must be far away from the initial positions of the spins or suitable boundary conditions must be found to allow the movement of spins across these exterior boundaries.Many approaches have been taken to solve the Bloch–Torrey equation but an efficient high-performance computing framework is still missing. In this paper, we present formulations of the interface as well as the exterior boundary conditions that are computationally efficient and suitable for arbitrary order finite elements and parallelization. In particular, the formulations are based on the partition of unity concept which allows for a discontinuous solution across interfaces conforming with the mesh with weak enforcement of real (in the case of interior interfaces) and artificial (in the case of exterior boundaries) permeability conditions as well as an operator splitting for the exterior boundary conditions. The method is straightforward to implement and it is available in FEniCS for moderate-scale simulations and in FEniCS-HPC for large-scale simulations. The order of accuracy of the resulting method is validated in numerical tests and a good scalability is shown for the parallel implementation. We show that the simulated dMRI signals offer good approximations to reference signals in cases where the latter are available and we performed simulations for a realistic model of a neuron to show that the method can be used for complex geometries.
Subject headings
- NATURVETENSKAP -- Matematik -- Beräkningsmatematik (hsv//swe)
- NATURAL SCIENCES -- Mathematics -- Computational Mathematics (hsv//eng)
Keyword
- Computational diffusion MRI
- Bloch–Torrey equation
- Partition of unity finite element method
- Interface conditions
- Weak pseudo-periodic conditions
- FEniCS/FEniCS-HPC
- Tillämpad matematik och beräkningsmatematik
- Applied and Computational Mathematics
- Biological Physics
- Biologisk fysik
- Computer Science
- Datalogi
Publication and Content Type
- ref (subject category)
- art (subject category)
Find in a library
To the university's database