We use the elliptic reconstruction technique in combination with a duality approach to prove apos... more 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({\Omega})) 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({\Omega}))-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.
We address the error control of Galerkin discretization (in space) of linear second order hyperbo... more 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
We present a continuous finite element method for some examples of fully nonlinear elliptic equat... more 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\`ere equation and Pucci's equation.
We present global existence results for solutions of reaction-diffusion systems on evolving domai... more We present global existence results for solutions of reaction-diffusion systems on evolving domains. Global existence results for a class of reaction-diffusion systems on fixed domains are extended to the same systems posed on spatially linear isotropically evolving domains. The results hold without any assumptions on the sign of the growth rate. The analysis is valid for many systems that commonly
Lecture Notes in Computational Science and Engineering, 2012
The 91st London Mathematical Society Durham Symposium took place from July 5th to 15th 2010, with... more The 91st London Mathematical Society Durham Symposium took place from July 5th to 15th 2010, with more than 100 international participants attending. The Symposium focused on Numerical Analysis of Multiscale Problems and this book contains 10 invited articles from ...
We use the elliptic reconstruction technique in combination with a duality approach to prove apos... more 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({\Omega})) 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({\Omega}))-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.
We address the error control of Galerkin discretization (in space) of linear second order hyperbo... more 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
We present a continuous finite element method for some examples of fully nonlinear elliptic equat... more 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\`ere equation and Pucci's equation.
We present global existence results for solutions of reaction-diffusion systems on evolving domai... more We present global existence results for solutions of reaction-diffusion systems on evolving domains. Global existence results for a class of reaction-diffusion systems on fixed domains are extended to the same systems posed on spatially linear isotropically evolving domains. The results hold without any assumptions on the sign of the growth rate. The analysis is valid for many systems that commonly
Lecture Notes in Computational Science and Engineering, 2012
The 91st London Mathematical Society Durham Symposium took place from July 5th to 15th 2010, with... more The 91st London Mathematical Society Durham Symposium took place from July 5th to 15th 2010, with more than 100 international participants attending. The Symposium focused on Numerical Analysis of Multiscale Problems and this book contains 10 invited articles from ...
Uploads
Papers by Omar Lakkis