## 1 Introduction

Most of the classical linear partial differential equations arising in mathematical physics can be written in a common operator form. It has been shown in

[4] that this form is an evolutionary problem, given by(1.1) |

where stands for the derivative with respect to time, and are bounded linear operators on some Hilbert space ,

is an unbounded skew-selfadjoint operator on

and is a given source term.We are interested in a unique solution of above equation. The theory presented in [4] deals vanishing initial conditions at , but any non-zero initial condition can be incorporated as data into . For this purpose let and define the weighted -function space

The space is a Hilbert space endowed with the natural inner product given by

for all , where is the inner product of . The associated weighted -function spaces are denoted by for . Now the solution theory of [4] states: If there exists a and a such that

then for all right hand sides exists a unique solution .

In [3] this class of problems was investigated numerically using a discontinuous Galerkin approach for the discretisation in time. Here we want to apply a continuous approach, namely the continuous Galerkin-Petrov method [7, 2, 5, 1]. Like in [3] we will consider a special subclass of the general class of evolutionary problems, although the approach of the time-discretisation is suitable in the general case.

For the precise class of problems let be the spatial dimension and be partitioned into measurable, disjoint sets and . Denoting by

the characteristic function of a domain

we define the linear operatorsAdditionally let the unbounded skew-selfadjoint operator A be given by

Here denotes the homogeneous Dirichlet boundary conditions. The Hilbert space can now be specified to

The resulting evolutionary problem is now of mixed type. More precisely, on we get an equation of elliptic type, on the equations becomes parabolic while on the problem is hyperbolic.

The solution theory of [4] does not only give for a final time and the existence of a unique solution . It holds in addition and therefore, by the Sobolev embedding theorem, is continuous in time. If the right-hand side is smoother, say , and and are in the range of , then the continuity in time of the solution is everywhere and we have, see [6]

In the following we want to utilise this additional regularity. Section 2 gives the precise formulation of the method considered. Section 3 then deals with the existence of discrete solutions, Section 4

presents error estimates and finally Section

5 gives some numerical examples and conclusions.## 2 Numerical method

Let us start wit a semi-discretisation in time. The semi-discrete variational form of (1.1) uses a decomposition of into disjoint intervals of length for and the piecewise polynomial-in-time spaces

for the trial and the test functions, where is the space of polynomials of degree up to on with values in . Let us localise the scalar product in to the time intervals by

Then the variational formulation using the continuous Galerkin-Petrov method reads:

Find such that for all and

(2.1) | ||||

where denotes the weighted -projection into and is the initial value.

In order to compute an approximation to the solution of problem (1.1), we need also to discretise the spatial domain and the function space . Let be discretised into by a regular simplicial mesh that resolves the sets and , and let be the maximal diameter of the cells of . Furthermore, let denote a polynomial degree. Then the discrete spaces are given by

where the spatial spaces are

Here is the space of polynomials of degree up to on the cell and is the Raviart-Thomas-space, defined by

The fully discrete method then reads similarly to (2.1):

Find such that for all and

(2.2) | ||||

where and is the initial value.

## 3 Existence of discrete solution

Let us start by choosing , the -projection of into the test space and denote by

a localised version of .

###### Lemma 3.1.

It holds for all

###### Proof.

Let us consider any interval . Then it holds

where the skew-symmetry of and the definition of was used.

For the first term we apply integration by parts and obtain due to the exponential weight

By the -orthogonality it follows

and therefore

With the general existence assumption and we obtain

(3.1) |

Summing over the intervals the statement follows. ∎

With Lemma 3.1 we have control over the -projection of the solution and of the final value of . In order to obtain control of the full solution we consider two approaches.

1) has a trivial null space

In this case we employ the norm equivalence in discrete spaces

with constants independent of and and we obtain

Applying a Young inequality we have

Therefore, in the case of with only trivial null-space the linear system is solvable and we also have continuous dependence of the solution on the data.

2) may have a non-trivial null space

Here we can apply another norm-equivalence to show the solvability of the discrete systems.
Using we have the norm equivalence

with constants independent of and . Using a similar approach as before we obtain

With a Young inequality follows

Therefore, also in the general case the linear system is solvable and the solution on each time interval depends continuously on the data and the initial value of the time interval.

## 4 Error-estimation

### 4.1 Semi discretisation in time

Note that we do have the Galerkin orthogonality

(4.1) |

for the solution of (1.1) for and of (2.1). We now want to estimate the error . Let

be the interpolation operator fulfilling locally for all

andWe decompose the error into , where

###### Lemma 4.1.

It holds for any

(4.2) |

###### Proof.

In the following estimates we will bound the terms of interest by a norm of the interpolation error given by

###### Lemma 4.2.

It holds for any

###### Proof.

Thus we have control over the weighted -projection of the error and its final value.

###### Corollary 4.3.

In the case of having only a trivial null-space we have immediate control of the full norm due to

for some implying

In the case of having a non-trivial null-space the norm equivalence considered in the existence proof of Section 3 do only provide a bound exponentially increasing in .

### 4.2 Full discretisation

Let be the maximal time-step length. Similarly to [3, Section 4] we can now proceed and estimate the error of the full discretisation. Here the only difference is the time-discretisation method.

###### Theorem 4.4.

###### Proof.

This is a direct consequence of the results in [3, Section 4] combined with the estimates of the previous section. ∎

## 5 Numerical examples

We consider two examples with unknown solutions. Simulations with known smooth solutions were also made and the theoretical orders were observed. The two following examples show a more realistic behaviour in the case of changing type systems. All computations were done in the finite-element framework , see github.com/SOFE-Developers/SOFE.

### 5.1 1+1d example

Let us consider as first example one spatial dimension and combine a hyperbolic and an elliptic region. To be more precise, let , , and . As final time we set . The problem is now given by

(5.1) |

with homogeneous Dirichlet-conditions for the first component of , the initial condition and a right-hand side , where is the characteristic function of the non-negative time line and

Thus, is continuous and it holds for . Therefore, the solution theory of [4] gives the existence of a unique solution that is continuous in time. Figure 1

shows plots of the components of the solution in the domain. Note that the first component has a kink along – it is continuous but not differentiable in . As mesh we use an equidistant mesh of cells in and cells in . In order to calculate the errors we use a reference solution instead of the unknown solution . The reference solution is computed on an mesh with polynomial degrees and . Table 1

cGP-method | dG-method | |||

, | ||||

256 | 8.890e-04 | 1.808e-04 | ||

512 | 3.136e-04 | 1.50 | 7.751e-05 | 1.22 |

1024 | 1.380e-04 | 1.18 | 3.496e-05 | 1.15 |

2048 | 6.739e-05 | 1.03 | 1.580e-05 | 1.15 |

, | ||||

256 | 1.187e-04 | 6.058e-05 | ||

512 | 5.489e-05 | 1.11 | 2.642e-05 | 1.20 |

1024 | 2.492e-05 | 1.14 | 1.137e-05 | 1.22 |

2048 | 1.114e-05 | 1.16 | 4.669e-06 | 1.28 |

shows the results for different values of and and polynomial degrees and .
We coupled as the theory gives for smooth the convergence order
if and are proportional. We observe for the continuous Galerkin-Petrov method
in the first column only a convergence rate between 1 and 2. Increasing the polynomial
degree reduces the error, but does not influence the rate much. A reason for this behaviour
could be that is not smooth enough for the error estimates to hold.

For comparison we also computed approximations with the discontinuous Galerkin method
from [3]. The errors given in the remaining columns show a similar behaviour
with convergence rates between 1 and 2. Nevertheless, the errors are smaller
for the discontinuous approach.

### 5.2 1+2d example

As second example we consider the last example of [3]. Let , , and The problem is given by

(5.2) |

where

Figure 2

shows some snapshots of the first component of the solution , approximated by a numerical simulation. Again we use equidistant meshes with cells in each dimension of space and cells in . As reference solution replacing the unknown exact solution we use an approximation calculated with and , resp. Table 2

cGP-method | dG-method | |||

, | ||||

16 | 1.961e-02 | 7.821e-03 | ||

32 | 9.199e-03 | 1.09 | 3.018e-03 | 1.37 |

64 | 3.751e-03 | 1.29 | 8.813e-04 | 1.78 |

96 | 1.324e-03 | 1.50 | 2.920e-04 | 1.59 |

, | ||||

16 | 5.499e-03 | 2.790e-03 | ||

32 | 1.435e-03 | 1.94 | 6.385e-04 | 2.13 |

64 | 4.430e-04 | 1.70 | 2.248e-04 | 1.51 |

shows the results. Similarly to the previous example we do not achieve the optimal convergence order for both methods. And again the discontinuous Galerkin method has smaller errors.

## Conclusions

The continuous solution of an evolutionary system with continuous right hand side can be approximated by several methods. Here we investigated the continuous Galerkin-Petrov method, that has optimal convergence order for smooth solutions. We have proved this for the -projection of the error into the test space. For operators

having only a trivial null space it also follows for the full error. The benefit of the continuous method compared to the discontinuous Galerkin method is the continuity that implies a non-dissipative behaviour. In our examples with unknown solutions, that are probably not smooth, the discontinuous Galerkin method is slightly better. Furthermore, these examples show that an increase of the polynomial degree in space over 2 and in time over 1 gives no huge benefit. This is different for smooth solutions – here both methods achieve the theoretical high convergence orders.

## References

- [1] N. Ahmed and G. Matthies. Higher order continuous Galerkin-Petrov time stepping schemes for transient convection-diffusion-reaction equations. ESAIM Math. Model. Numer. Anal., 49(5):1429–1450, 2015.
- [2] A. K. Aziz and P. Monk. Continuous finite elements in space and time for the heat equation. Math. Comp, 52(186):255–274, 1989.
- [3] S. Franz, S. Trostorff, and M. Waurick. Numerical methods for changing type systems. IMAJNA, 39(2):1009–1038, 2019.
- [4] R. Picard. A structural observation for linear material laws in classical mathematical physics. Math. Methods Appl. Sci., 32(14):1768–1803, 2009.
- [5] F. Schieweck. -stable discontinuous Galerkin-Petrov time discretization of higher order. J. Numer. Math., 18(1):25–57, 2010.
- [6] S. Trostorff and M. Waurick. On higher index differential-algebraic equations in infinite dimensions. In A. Böttcher, D. Potts, P. Stollmann, and D. Wenzel, editors, The Diversity and Beauty of Applied Operator Theory, pages 477–486, Cham, 2018. Springer International Publishing.
- [7] R. Winther. A stable finite element method for initial-boundary value problems for first-order hyperbolic systems. Math. Comp., 36(153):65–86, 1981.