Omar Lakkis

Emmanuil H. Georgoulis, Omar Lakkis, Juha M. Virtanen

We derive energy-norm a posteriori error bounds for an Euler time-stepping method combined with various spatial discontinuous Galerkin schemes for linear parabolic problems. For accessibility, we address first the spatially semidiscrete case, and then move to the fully discrete scheme by introducing the implicit Euler time-stepping. All results are presented in an abstract setting and then illustrated with particular applications. This enables the error bounds to hold for a variety of discontinuous Galerkin methods, provided that energy-norm a posteriori error bounds for the corresponding elliptic problem are available. To illustrate the method, we apply it to the interior penalty discontinuous Galerkin method, which requires the derivation of novel a posteriori error bounds. For the analysis of the time-dependent problems we use the elliptic reconstruction technique and we deal with the nonconforming part of the error by deriving appropriate computable a posteriori bounds for it.

Markos A. Katsoulakis, Georgios T. Kossioris, Omar Lakkis

We address the numerical discretization of the Allen-Cahn prob- lem with additive white noise in one-dimensional space. The discretization is conducted in two stages: (1) regularize the white noise and study the regularized problem, (2) approximate the regularized problem. We address (1) by introducing a piecewise constant random approximation of the white noise with respect to a space-time mesh. We analyze the regularized problem and study its relation to both the original problem and the deterministic Allen-Cahn problem. Step (2) is then performed leading to a practical Monte-Carlo method combined with a Finite Element-Implicit Euler scheme. The resulting numerical scheme is tested against theoretical benchmark results.

Omar Lakkis, Anotida Madzvamuse, Chandrasekhar Venkataraman

We present and analyse an implicit-explicit timestepping procedure with finite element spatial approximation for a semilinear reaction-diffusion systems on evolving domains arising from biological models, such as Schnakenberg's (1979). We employ a Lagrangian formulation of the model equations which permits the error analysis for parabolic equations on a fixed domain but introduces technical difficulties, foremost the space-time dependent conductivity and diffusion. We prove optimal-order error estimates in the $\Lp{\infty}(0,T;\Lp{2}(\W))$ and $\Lp{2}(0,T;\Hil{1}(\W))$ norms, and a pointwise stability result. We remark that these apply to Eulerian solutions. Details on the implementation of the Lagrangian and the Eulerian scheme are provided. We also report on a numerical experiment for an application to pattern formation on an evolving domain.

Alan Demlow, Omar Lakkis, Charalambos Makridakis

We derive a posteriori error estimates in the $L_\infty((0,T];L_\infty(Ω))$ norm for approximations of solutions to linear para bolic equations. Using the elliptic reconstruction technique introduced by Makridakis and Nochetto and heat kernel estimates for linear parabolic pr oblems, we first prove a posteriori bounds in the maximum norm for semidiscrete finite element approximations. We then establish a posteriori bounds for a fully discrete backward Euler finite element approximation. The elliptic reconstruction technique greatly simplifies our development by allow\ ing the straightforward combination of heat kernel estimates with existing elliptic maximum norm error estimators.

Emmanuil H. Georgoulis, Omar Lakkis, Charalambos Makridakis

We address the error control of Galerkin discretization (in space) of linear second order hyperbolic problems. More specifically, we derive a posteriori error bounds in the L\infty(L2)-norm for finite element methods for the linear wave equation, under minimal regularity assumptions. The theory is developed for both the space-discrete case, as well as for an implicit fully discrete scheme. The derivation of these bounds relies crucially on carefully constructed space- and time-reconstructions of the discrete numerical solutions, in conjunction with a technique introduced by Baker (1976, SIAM J. Numer. Anal., 13) in the context of a priori error analysis of Galerkin discretization of the wave problem in weaker-than-energy spatial norms.

Omar Lakkis, Tristan Pryer

We present a continuous finite element method for some examples of fully nonlinear elliptic equation. A key tool is the discretisation proposed in Lakkis & Pryer (2011, SISC) allowing us to work directly on the strong form of a linear PDE. An added benefit to making use of this discretisation method is that a recovered (finite element) Hessian is a biproduct of the solution process. We build on the linear basis and ultimately construct two different methodologies for the solution of second order fully nonlinear PDEs. Benchmark numerical results illustrate the convergence properties of the scheme for some test problems including the Monge-Ampère equation and Pucci's equation.

Omar Lakkis, Tristan Pryer

We propose a numerical method to approximate the solution of second order elliptic problems in nonvariational form. The method is of Galerkin type using conforming finite elements and applied directly to the nonvariational (nondivergence) form of a second order linear elliptic problem. The key tools are an appropriate concept of a 'finite element Hessian' and a Schur complement approach to solving the resulting linear algebra problem. The method is illustrated with computational experiments on three linear and one quasilinear PDE, all in nonvariational form.

Adrian Muntean, Omar Lakkis

We study a two-scale reaction-diffusion system with nonlinear reaction terms and a nonlinear transmission condition (remotely ressembling Henry's law) posed at air-liquid interfaces. We prove the rate of convergence of the two-scale Galerkin method proposed in Muntean & Neuss-Radu (2009) for approximating this system in the case when both the microstructure and macroscopic domain are two-dimensional. The main difficulty is created by the presence of a boundary nonlinear term entering the transmission condition. Besides using the particular two-scale structure of the system, the ingredients of the proof include two-scale interpolation-error estimates, an interpolation-trace inequality, and improved regularity estimates.

Omar Lakkis, Charalambos Makridakis, Tristan Pryer

We use the elliptic reconstruction technique in combination with a duality approach to prove aposteriori error estimates for fully discrete back- ward Euler scheme for linear parabolic equations. As an application, we com- bine our result with the residual based estimators from the aposteriori esti- mation for elliptic problems to derive space-error indicators and thus a fully practical version of the estimators bounding the error in the L \infty (0, T ; L2(Ω)) norm. These estimators, which are of optimal order, extend those introduced by Eriksson and Johnson (1991) by taking into account the error induced by the mesh changes and allowing for a more flexible use of the elliptic estima- tors. For comparison with previous results we derive also an energy-based aposteriori estimate for the L \infty (0, T ; L2(Ω))-error which simplifies a previous one given in Lakkis and Makridakis (2006). We then compare both estimators (duality vs. energy) in practical situations and draw conclusions.

Ellya Kawecki, Omar Lakkis, Tristan Pryer

We address the numerical solution via Galerkin type methods of the Monge-Ampère equation with transport boundary conditions arising in optimal mass transport, geometric optics and computational mesh or grid movement techniques. This fully nonlinear elliptic problem admits a linearisation via a Newton-Raphson iteration, which leads to an oblique derivative boundary value problem for elliptic equations in nondivergence form. We discretise these by employing the nonvariational finite element method, which lead to empirically observed optimal convergence rates, provided recovery techinques are used to approximate the gradient and the Hessian of the unknown functions. We provide extensive numerical testing to illustrate the strengths of our approach and the potential applications in optics and mesh movement.

Emmanuil H. Georgoulis, Omar Lakkis, Charalambos Makridakis et al.

We consider second order explicit and implicit two-step time-discrete schemes for wave-type equations. We derive optimal order aposteriori estimates controlling the time discretization error. Our analysis, has been motivated by the need to provide aposteriori estimates for the popular leap-frog method (also known as Verlet's method in molecular dynamics literature); it is extended, however, to general cosine-type second order methods. The estimators are based on a novel reconstruction of the time-dependent component of the approximation. Numerical experiments confirm similarity of convergence rates of the proposed estimators and of the theoretical convergence rate of the true error.

Omar Lakkis, Tristan Pryer

We derive energy-norm aposteriori error bounds, using gradient recovery (ZZ) estimators to control the spatial error, for fully discrete schemes for the linear heat equation. This appears to be the first completely rigorous derivation of ZZ estimators for fully discrete schemes for evolution problems, without any restrictive assumption on the timestep size. An essential tool for the analysis is the elliptic reconstruction technique. Our theoretical results are backed with extensive numerical experimentation aimed at (a) testing the practical sharpness and asymptotic behaviour of the error estimator against the error, and (b) deriving an adaptive method based on our estimators. An extra novelty provided is an implementation of a coarsening error "preindicator", with a complete implementation guide in ALBERTA.

Emmanuil H. Georgoulis, Omar Lakkis

We derive a posteriori error bounds for a quasilinear parabolic problem, which is approximated by the $hp$-version interior penalty discontinuous Galerkin method (IPDG). The error is measured in the energy norm. The theory is developed for the semidiscrete case for simplicity, allowing to focus on the challenges of a posteriori error control of IPDG space-discretizations of strictly monotone quasilinear parabolic problems. The a posteriori bounds are derived using the elliptic reconstruction framework, utilizing available a posteriori error bounds for the corresponding steady-state elliptic problem.