What is the difference between edge finite elements and nodal finite elements?

This for use in modeling the eddy current problem in classical electromagnetism. I am attempting to convert MATLAB code written for Electrical Impedance Tomography (EIT) into a code for Magnetic Induction Tomography (MIT). EIT involves electrodes with applied alternating currents and requires nodal FEM. Whereas MIT involves coils and induction of eddy currents and requires edge FEM.

The code for both EIT and MIT involves solving the forward and inverse problems in a simulation for tomographic image reconstruction of conductivity maps.

The MATLAB code for conversion is found at EIDORS:

  • 1
    $\begingroup$ Have you checked chapter 3 from the FEnICs book? $\endgroup$ – nicoguaro Feb 6 '17 at 22:53
  • $\begingroup$ No I haven't I will check it out. $\endgroup$ – Brendan Darrer Feb 7 '17 at 9:50

The paper: I. Anjam, J. Valdman, "Fast MATLAB assembly of FEM matrices in 2D and 3D: Edge elements", Applied Mathematics and Computation, 267, (2015), 252–263; states the following. http://www.sciencedirect.com/science/article/pii/S0096300315004191

"A finite element discretization is done in terms of edge elements, typically Raviart–Thomas elements [12] for H(div) problems and Nédélec elements [9] for H(curl) problems. Edge element basis functions are not defined on the nodes of 2D triangular or 3D tetrahedral meshes, but on edges and faces. Edge elements provide only partial continuity over element boundaries: continuity of normal vector component for H(div) problems and continuity of tangential vector component for H(curl) problems."

They give a figure for this representation.

enter image description here

Edge elements are derived from nodal shape functions:

See 4.1.3 Shape Function

"In Edge FEM on a tetrahedral mesh, a vector field is represented using a basis vector function" (1):

$N_{ij} = L_i\bigtriangledown L_j - L_j\bigtriangledown L_i$ ____________ (1)

The linear shape function of the node $i$ has the form,

$L_i (x,y,z) = a_i + b_i x + c_i y + d_i z$

enter image description here

(1) L. Ma, “Magnetic Induction Tomography for Non-destructive Evaluation and Process Tomography,” University of Bath, 2014, pp.12.


Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.