Fortin–Soulie
Click here to read what the information on this page means.
Implementations
Symfem | "Fortin-Soulie" ↓ Show Symfem examples ↓↑ Hide Symfem examples ↑Before trying this example, you must install Symfem: pip3 install symfem This element can then be created with the following lines of Python: import symfem
# Create Fortin-Soulie order 2 on a triangle element = symfem.create_element("triangle", "Fortin-Soulie", 2) |
Examples
triangle
order 2
- \(R\) is the reference triangle. The following numbering of the subentities of the reference is used:
- \(\mathcal{V}\) is spanned by: \(1\), \(x\), \(x^{2}\), \(y\), \(x y\), \(y^{2}\)
- \(\mathcal{L}=\{l_0,...,l_{5}\}\)
- Functionals and basis functions:
\(\displaystyle l_{0}:v\mapsto\displaystyle\int_{e_{0}}(1 - s_{0})v\)
where \(e_{0}\) is the 0th edge;
and \(s_{0}\) is a parametrisation of \(e_{0}\).
\(\displaystyle \phi_{0} = \sqrt{2} \cdot \left(4 x^{2} - 2 x y - x - 2 y^{2} + 2 y - \frac{1}{3}\right)\)
This DOF is associated with edge 0 of the reference element.
\(\displaystyle l_{1}:v\mapsto\displaystyle\int_{e_{0}}(s_{0})v\)
where \(e_{0}\) is the 0th edge;
and \(s_{0}\) is a parametrisation of \(e_{0}\).
\(\displaystyle \phi_{1} = \sqrt{2} \cdot \left(12 x y - 3 x + 6 y^{2} - 6 y + 1\right)\)
This DOF is associated with edge 0 of the reference element.
\(\displaystyle l_{2}:v\mapsto\displaystyle\int_{e_{1}}(1 - s_{0})v\)
where \(e_{1}\) is the 1st edge;
and \(s_{0}\) is a parametrisation of \(e_{1}\).
\(\displaystyle \phi_{2} = 12 x^{2} + 24 x y - 18 x - 6 y + 4\)
This DOF is associated with edge 1 of the reference element.
\(\displaystyle l_{3}:v\mapsto\displaystyle\int_{e_{1}}(s_{0})v\)
where \(e_{1}\) is the 1st edge;
and \(s_{0}\) is a parametrisation of \(e_{1}\).
\(\displaystyle \phi_{3} = - 4 x^{2} - 28 x y + 10 x - 4 y^{2} + 10 y - \frac{8}{3}\)
This DOF is associated with edge 1 of the reference element.
\(\displaystyle l_{4}:v\mapsto\displaystyle\int_{e_{2}}(1 - s_{0})v\)
where \(e_{2}\) is the 2nd edge;
and \(s_{0}\) is a parametrisation of \(e_{2}\).
\(\displaystyle \phi_{4} = - 4 x^{2} - 4 x y + 4 x + 8 y^{2} - 8 y + \frac{4}{3}\)
This DOF is associated with edge 2 of the reference element.
\(\displaystyle l_{5}:v\mapsto v(\tfrac{1}{3},\tfrac{1}{3})\)
\(\displaystyle \phi_{5} = - 6 x^{2} - 6 x y + 6 x - 6 y^{2} + 6 y - 1\)
This DOF is associated with face 0 of the reference element.
References
- Fortin, Michel and Soulie, M. A nonāconforming piecewise quadratic finite element on triangles, International Journal for Numerical Methods in Engineering 19, 505–520, 1983. [DOI: 10.1002/nme.1620190405] [BibTeX]
DefElement stats
Element added | 06 March 2021 |
Element last updated | 10 February 2022 |