A numerical meshless particle method in solving the magnetoencephalography forward problem

Abstract

In this paper, a numerical meshless particle method is presented in order to solve the magnetoencephalography\ud forward problem for analyzing the complex activation patterns in the human brain. The forward problem is\ud devoted to compute the scalp potential and magnetic field distribution generated by a set of current sources\ud representing the neural activity, and in this paper, it has been approached by means of the smoothed particle\ud hydrodynamics method suitably handled. The Poisson equation generated by the quasi-stationary Maxwell’s curl\ud equations, by assuming Neumann boundary conditions has been considered, and the current sources have been\ud simulated by current dipoles. The adopted meshless particle model has provided good results in agreement with\ud the analytical ones and by overcoming the drawback of the mesh generation. The numerical model has been validated,\ud at first, in computing the electric potential and the external magnetic field for a dipole plunged near the upper\ud surface of a homogeneous sphere simulating the human brain. Simulation results obtained by simulating two\ud concentric spheres with different conductivities are also reported. Moreover, in order to better assess the validity of\ud the proposed approach, a realistic human brain cortex model has been also simulated and compared with boundary\ud element method results. A satisfactory agreement has been reached

Toward a computational multiresolution analysis for radial point interpolation meshless method

“…The construction of functional in calculus of variations is similar in both FEM and RPIM based on the following inner product definition [10]. For two functions u and a, the inner product is defined as [10] hu; ai D Z u: N ad (1) where N a denotes the complex conjugate of a. Let us consider an operator equation as…”

“…A meshless method is a method used to establish the system of algebraic equations for the whole problem domain without the use of a predefined mesh for the domain discretization. Many meshless methods have found good applications in electromagnetics and shown very good potential to become powerful numerical tools [1,2]. The development of meshless methods can be traced back more than 70 years to the collocation methods.…”

Mixed basis function for radial point interpolation meshless method in electromagnetics

“…Meshless methods may also solve some other problems associated with the finite element method, such as locking and element distortion. [3][4][5][6][7] Some popular meshless techniques are the meshless local Petrov-Galerkin method (MLPG) and radial point interpolation meshless method (RPIM).…”

A Meshfree Solver for the MEG Forward Problem

“…Meshless methods have been previously proposed for solving the EEG [4] and the MEG [5] forward problem. However, they are domain methods, thus they may be outperformed by BEM from a computational efficiency standpoint.…”