Raviart–Thomas basis functions

From HandWiki

In applied mathematics, Raviart–Thomas basis functions are vector basis functions used in finite element and boundary element methods. They are regularly used as basis functions when working in electromagnetics. They are sometimes called Rao-Wilton-Glisson basis functions.[1] The space [math]\displaystyle{ \mathrm{RT}_q }[/math] spanned by the Raviart–Thomas basis functions of order [math]\displaystyle{ q }[/math] is the smallest polynomial space such that the divergence maps [math]\displaystyle{ \mathrm{RT}_q }[/math] onto [math]\displaystyle{ \mathrm{P}_q }[/math], the space of piecewise polynomials of order [math]\displaystyle{ q }[/math].[2]

Order 0 Raviart-Thomas Basis Functions in 2D

Raviart thomas labelled.png

In two-dimensional space, the lowest order Raviart Thomas space, [math]\displaystyle{ \mathrm{RT}_0 }[/math], has degrees of freedom on the edges of the elements of the finite element mesh. The [math]\displaystyle{ n }[/math]th edge has an associated basis function defined by[3]

[math]\displaystyle{ \mathbf{f}_n(\mathbf{r})=\left\{\begin{array}{ll} \frac{l_n}{2A_n^+}(\mathbf{r}-\mathbf{p}_+)\quad&\mathrm{if\ \mathbf{r}\in\ T_+}\\ -\frac{l_n}{2A_n^-}(\mathbf{r}-\mathbf{p}_-)\quad&\mathrm{if\ \mathbf{r}\in\ T_-}\\ \mathbf{0}\quad&\mathrm{otherwise} \end{array}\right. }[/math]

where [math]\displaystyle{ l_n }[/math] is the length of the edge, [math]\displaystyle{ T_+ }[/math] and [math]\displaystyle{ T_- }[/math] are the two triangles adjacent to the edge, [math]\displaystyle{ A_n^+ }[/math] and [math]\displaystyle{ A_n^- }[/math] are the areas of the triangles and [math]\displaystyle{ \mathbf{p}_+ }[/math] and [math]\displaystyle{ \mathbf{p}_- }[/math] are the opposite corners of the triangles.

Sometimes the basis functions are alternatively defined as

[math]\displaystyle{ \mathbf{f}_n(\mathbf{r})=\left\{\begin{array}{ll} \frac{1}{2A_n^+}(\mathbf{r}-\mathbf{p}_+)\quad&\mathrm{if\ \mathbf{r}\in\ T_+}\\ -\frac{1}{2A_n^-}(\mathbf{r}-\mathbf{p}_-)\quad&\mathrm{if\ \mathbf{r}\in\ T_-}\\ \mathbf{0}\quad&\mathrm{otherwise} \end{array}\right. }[/math]

with the length factor not included.

References

  1. Andriulli, Francesco P.; Cools; Bagci; Olyslager; Buffa; Christiansen; Michelssen (2008). "A Mulitiplicative Calderon Preconditioner for the Electric Field Integral Equation". IEEE Transactions on Antennas and Propagation 56 (8): 2398–2412. doi:10.1109/tap.2008.926788. Bibcode2008ITAP...56.2398A. 
  2. Logg, Anders, ed (2012). "Chapter 3. Common and unusual finite elements". Automated Solution of Differential Equations by the Finite Element Method. Lecture Notes in Computational Science and Engineering. 84. Berlin, Heidelberg: Springer Berlin Heidelberg. pp. 95–119. doi:10.1007/978-3-642-23099-8. ISBN 978-3-642-23098-1. http://link.springer.com/10.1007/978-3-642-23099-8. 
  3. Bahriawati, C.; Carstensen, C. (2005). "Three MATLAB Implementations Of The Lowest-Order Raviart-Thomas MFEM With a Posteriori Error Control". Computational Methods in Applied Mathematics 5 (4): 331–361. doi:10.2478/cmam-2005-0016. http://www2.mathematik.hu-berlin.de/~cc/cc_homepage/download/2005-BC_CC-Three_MATLAB_Implementations_Lowest-Order_Raviart-Thomas_MFEM.pdf. Retrieved 8 October 2015.