BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Omar Lakkis

Omar Lakkis appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NAmath.APnlin.PSmath.DSNumerical Analysismath.PR

Catalog footprint

What is connected

13works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

13 published item(s)

preprint2025arXiv

Adaptive FEM with explicit time integration for the wave equation

Starting from a recent a posteriori error estimator for the finite element solution of the wave equation with explicit time-stepping [Grote, Lakkis, Santos, 2024], we devise a space-time adaptive strategy which includes both time evolving meshes and local time-stepping [Diaz, Grote, 2009] to overcome any overly stringent CFL stability restriction on the time-step due to local mesh refinement. Moreover, at each time-step the adaptive algorithm monitors the accuracy thanks to the error indicators and recomputes the current step on a refined mesh until the desired tolerance is met; meanwhile, the mesh is coarsened in regions of smaller errors. Leapfrog based local time-stepping is applied in all regions of local mesh refinement to incorporate adaptivity into fully explicit time integration with mesh change while retaining efficiency. Numerical results illustrate the optimal rate of convergence of the a posteriori error estimators on time evolving meshes.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Saint-Venant shallow water model for overland flows with precipitation and recharge

We propose a one-dimensional Saint-Venant (open channel) model overland flows including a water input--output source term modelling recharge via rainfall and infiltration (or exfiltration). We derive the model via asymptotic reduction from the two-dimensional Navier--Stokes equations under the shallow water assumption, with boundary conditions including recharge via ground infiltration and runoff. This new model recovers existing models as sepcial cases, and adds more scope by adding a water-mixing friction terms that depends on the rate of water recharge. We prospose a novel entropy function and its flux, that are useful in validating the model's conservation or dissipation properties. Based on this entropy function we propose a finite volume scheme extending a class of kinetic schemes and provide numerical comparisons with respect to the newly introduced mixing friction coefficient. We also provide a comparison with experimental data.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Implicit-explicit timestepping with finite element approximation of reaction-diffusion systems on evolving domains

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A finite element method for fully nonlinear elliptic problems

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Noise regularization and computations for the 1-dimensional stochastic Allen-Cahn problem

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A comparison of duality and energy aposteriori estimates for L?(0,T;L2(Ω)) in parabolic problems

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A finite element method for second order nonvariational elliptic problems

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A posteriori error bounds for discontinuous Galerkin methods for quasilinear parabolic problems

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A posteriori error control for discontinuous Galerkin methods for parabolic problems

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Global existence for semilinear reaction-diffusion systems on evolving domains

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 arise in the theory of pattern formation. We present numerical results illustrating our theoretical findings.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Gradient recovery in adaptive finite element methods for parabolic problems

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.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

A posteriori error estimates in the maximum norm for parabolic problems

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.

0 citations0 reviewsNo code linkedOpen access