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

Tristan Pryer

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

ResearcherUnclaimed source record
math.NAmath.APNumerical Analysis

Catalog footprint

What is connected

17works
3topics
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

17 published item(s)

preprint2026arXiv

A finite element method preserving the eigenvalue range of symmetric tensor fields

This paper presents a finite element method that preserves (at the degrees of freedom) the eigenvalue range of the solution of tensor-valued time-dependent convection--diffusion equations. Starting from a high-order spatial baseline discretisation (in this case, the CIP stabilised finite element method), our approach formulates the fully discrete problem as a variational inequality posed on a closed convex set of tensor-valued functions that respect the same eigenvalue bounds at their degrees of freedom. The numerical realisation of the scheme relies on the definition of a projection that, at each node, performs the diagonalisation of the tensor and then truncates the eigenvalues to lie within the prescribed bounds. The temporal discretisation is carried out using the implicit Euler method, and unconditional stability and optimal-order error estimates are proven for this choice. Numerical experiments confirm the theoretical findings and illustrate the method's ability to maintain eigenvalue constraints while accurately approximating solutions in the convection-dominated regime.

0 citations0 reviewsNo code linkedOpen access
preprint2024arXiv

Discretisation of an Oldroyd-B viscoelastic fluid flow using a Lie derivative formulation

In this article we present a numerical method for the Stokes flow of an Oldroyd-B fluid. The viscoelastic stress evolves according to a constitutive law formulated in terms of the upper convected time derivative. A finite difference method is used to discretise along fluid trajectories to approximate the advection and deformation terms of the upper convected derivative in a simple, cheap and cohesive manner, as well as ensuring that the discrete conformation tensor is positive definite. A full implementation with coupling to the fluid flow is presented, along with detailed discussion of the issues that arise with such schemes. We demonstrate the performance of this method with detailed numerical experiments in a lid-driven cavity setup. Numerical results are benchmarked against published data, and the method is shown to perform well in this challenging case.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Adaptive modelling of variably saturated seepage problems

In this article we present a goal-oriented adaptive finite element method for a class of subsurface flow problems in porous media, which exhibit seepage faces. We focus on a representative case of the steady state flows governed by a nonlinear Darcy-Buckingham law with physical constraints on subsurface-atmosphere boundaries. This leads to the formulation of the problem as a variational inequality. The solutions to this problem are investigated using an adaptive finite element method based on a dual-weighted a posteriori error estimate, derived with the aim of reducing error in a specific target quantity. The quantity of interest is chosen as volumetric water flux across the seepage face, and therefore depends on an a priori unknown free boundary. We apply our method to challenging numerical examples as well as specific case studies, from which this research originates, illustrating the major difficulties that arise in practical situations. We summarise extensive numerical results that clearly demonstrate the designed method produces rapid error reduction measured against the number of degrees of freedom.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A numerical implementation of the unified Fokas transform for evolution problems on a finite interval

We present the numerical solution of two-point boundary value problems for a third order linear PDE, representing a linear evolution in one space dimension. The difficulty of this problem is in the numerical imposition of the boundary conditions, and to our knowledge, no such computations exist. Instead of computing the evolution numerically, we evaluate the solution representation formula obtained by the unified transform of Fokas. This representation involves complex line integrals, but in order to evaluate these integrals numerically, it is necessary to deform the integration contours using appropriate deformation mappings. We formulate a strategy to implement effectively this deformation, which allows us to obtain accurate numerical results.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

An a posteriori analysis of some inconsistent, nonconforming Galerkin methods approximating elliptic problems

In this work we present an a posteriori analysis for classes of inconsistent, nonconforming schemes approximating elliptic problems. We show the estimates coincide with existing ones for interior penalty type discontinuous Galerkin approximations of the Laplacian and give new estimates for inconsistent discontinuous Galerkin approximation schemes of elliptic problems under quadrature approximation. We also examine the effect of inconsistencies on the a posteriori analysis of schemes applied to an unbalanced problem.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Discontinuous Galerkin methods for the $p$--biharmonic equation from a discrete variational perspective

We study discontinuous Galerkin approximations of the $p$--biharmonic equation from a variational perspective. We propose a discrete variational formulation of the problem based on a appropriate definition of a finite element Hessian and study convergence of the method (without rates) using a weak lower semicontinuity argument. We present numerical experiments aimed at testing the robustness of the method. We also note a superconvergence effect for some values of $p$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Noether type discrete conserved quantities arising from a finite element approximation of a variational problem

In this work we prove a weak Noether type theorem for a class of variational problems which include broken extremals. We then use this result to prove discrete Noether type conservation laws for certain classes of finite element discretisation of a model elliptic problem. In addition we study how well the finite element scheme satisfies the continuous conservation laws arising from the application of Noether's 1st Theorem (E. Noether 1918). We summarise extensive numerical tests, illustrating the conservativity of the discrete Noether law using the $p$--Laplacian as an example.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On the Numerical Approximation of $\infty$-Harmonic Mappings

Given a map $u : Ω\subseteq \mathbb{R}^n \longrightarrow \mathbb{R}^N$, the $\infty$-Laplacian is the system \[ \label{1} Δ_\infty u \, :=\, \Big(\text{D}u \otimes \text{D}u + |\text{D}u|^2 [\text{D}u]^\bot \! \otimes I \Big) : \text{D}^2 u\, = \, 0. \tag{1} \] \eqref{1} is the model system of vectorial Calculus of Variations in $L^\infty$ and arises as the "Euler-Lagrange" equation in relation to the supremal functional \[ \label{2} E_\infty(u,Ω)\, :=\, \| \text{D}u \|_{L^\infty(Ω)}. \tag{2} \] The scalar case of \eqref{1} has been introduced by Aronsson in the 1960s and by now is relatively classical and well understood. The general system \eqref{1} has been discovered and studied by the first author in a series of recent papers. Supremal functionals are fundamental for applications because they provide more realistic models as opposed to conventional integral models. Herein we provide numerical approximations of solutions to the Dirichlet problem when $n=2$ and $N=2,3$ for certain carefully selected boundary data on the unit square. Our experiments demonstrate interesting and unexpected phenomena occurring in $L^\infty$ and provide insights on the structure of general solutions and the natural separation to phases they present.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Energy consistent discontinuous Galerkin methods for a quasi-incompressible diffuse two phase flow model

We design consistent discontinuous Galerkin finite element schemes for the approximation of a quasi-incompressible two phase flow model of Allen-Cahn/Cahn-Hilliard/Navier-Stokes-Korteweg type which allows for phase transitions. We show that the scheme is mass conservative and monotonically energy dissipative. In this case the dissipation is isolated to discrete equivalents of those effects already causing dissipation on the continuous level, that is, there is no artificial numerical dissipation added into the scheme. In this sense the methods are consistent with the energy dissipation of the continuous PDE system.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Discontinuous Galerkin methods for nonvariational problems

We extend the finite element method introduced by Lakkis and Pryer [2011] to approximate the solution of second order elliptic problems in nonvariational form to incorporate the discontinuous Galerkin (DG) framework. This is done by viewing the NVFEM as a mixed method whereby the finite element Hessian is an auxiliary variable in the formulation. Representing the finite element Hessian in a discontinuous setting yields a linear system of the same size and having the same sparsity pattern of the compact DG methods for variational elliptic problems. Furthermore, the system matrix is very easy to assemble, Thus this approach greatly reduces the computational complexity of the discretisation compared to the continuous approach. We conduct a stability and consistency analysis making use of the unified framework set out in Arnold et. al. [2001]. We also give an apriori analysis of the method. The analysis applies to any consistent representation of the finite element Hessian, thus is applicable to the previous works making use of continuous Galerkin approximations.

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
preprint2012arXiv

Energy consistent DG methods for the Navier-Stokes-Korteweg system

We design consistent discontinuous Galerkin finite element schemes for the approximation of the Euler-Korteweg and the Navier-Stokes-Korteweg systems. We show that the scheme for the Euler-Korteweg system is energy and mass conservative and that the scheme for the Navier-Stokes-Korteweg system is mass conservative and monotonically energy dissipative. In this case the dissipation is isolated to viscous effects, that is, there is no numerical dissipation. In this sense the methods is consistent with the energy dissipation of the continuous PDE systems.

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

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