# Discretised Euler-Lagrange equations

## The Morley element

Now that we have existence for our problem we would like to solve it numerically. For this we use finite element method. We will be using a combination of Morley elements (for ) and piecewise affine elements (for ). The Morley elements are non-conforming ones, which means that the discretised function space is not a subspace of the original space. In fact, the solution of the discretised problem may not even be continuous. However, the advantage of this approach is that we can use lower order polynomial basis functions in each triangle, which makes the implementation of the numerical scheme much easier.

The Morley element is used widely throughout different elliptic PDE problems (see e.g. [1], [2], [4]). However, there are no results known to authors where a combination of Morley elements and piecewise affine elements is used in the same system.

We approximate by a polyhedral domain . For a subset let denote the space of quadratic polynomials on . Let be a triangulation of the computational domain .

**Definition: **The Morley Element (see Figure 1) is the triple where is a triangular element, is the space of quadratic polynomials defined on , and is the set of nodal variables, where

Figure 1: The Morley element

More details about the Morley element can be found in [3] and [1]. From here we can define the discretised function set:

We note is the cartesian product of the Morley element space (for ) and the piecewise linear element space (for ).

Denote by the nodes of the triangulation that are located on the boundary of and let denote the midpoints of the edges located at the boundary of .

Let

....,

where , , and is the outward normal vector at , .

## The discretised Euler-Lagrange equations

Consider a symmetric, bilinear form defined by

The discretised Euler-Lagrange equations are:

Find such that

.

This problem has a unique solution if . We note that this is a weaker condition than , but we hope that this result could be improved. For the proof see Section 5.1.3 of the RSG report.

We do not prove convergence of the solution to the discretised problem to the exact solution. Some results regarding convergence of the Morley element method can be found in [4], [2], [5], [6] and [1], but the analysis of the convergence of the combination of Morley elements and piecewise affine elements exceeded the time-frame of the research study group. Instead we will present the numerical convergence results.

**References**

[1] M. Wang and J. Xu. The Morley element for fourth order elliptic equations in any dimensions. Numerische Mathematik, 2004.

[2] Shen, Q. A posteriori error estimates of the Morley element for the

fourth order elliptic eigenvalue problem

[3] Larson, M. G. and Bengzon, F. The Finite Element Method: Theory, Implementation and Applications

[4] Wang, L. H. Morley's element approximation to a fourth order variational

inequality with curvature obstacle

[5] Hu, J. and Shi, Z. and Xu, J. Convergence and optimality of the adaptive Morley element

method

[6] Li, M. and Guan, X. and Mao, S. New error estimates of the Morley element for the plate

bending problems