Source author record

Martin Stynes

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

9works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

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

9 published item(s)

preprint2022arXiv

A posteriori error analysis for variable-coefficient multiterm time-fractional subdiffusion equations

An initial-boundary value problem of subdiffusion type is considered; the temporal component of the differential operator has the form $\sum_{i=1}^{\ell}q_i(t)\, D _t ^{α_i} u(x,t)$, where the $q_i$ are continuous functions, each $D _t ^{α_i}$ is a Caputo derivative, and the $α_i$ lie in $(0,1]$. Maximum/comparison principles for this problem are proved under weak hypotheses. A new positivity result for the multinomial Mittag-Leffler function is derived. A posteriori error bounds are obtained in $L_2(Ω)$ and $L_\infty(Ω)$, where the spatial domain $Ω$ lies in $\bR^d$ with $d\in\{1,2,3\}$. An adaptive algorithm based on this theory is tested extensively and shown to yield accurate numerical solutions on the meshes generated by the algorithm.

preprint2022arXiv

Optimal long-time decay rate of solutions of complete monotonicity-preserving schemes for nonlinear time-fractional evolutionary equations

The solution of the nonlinear initial-value problem $\mathcal{D}_{t}^αy(t)=-λy(t)^γ$ for $t>0$ with $y(0)>0$, where $\mathcal{D}_{t}^α$ is a Caputo derivative of order $α\in (0,1)$ and $λ, γ$ are positive parameters, is known to exhibit $O(t^{α/γ})$ decay as $t\to\infty$. No corresponding result for any discretisation of this problem has previously been proved. In the present paper it is shown that for the class of complete monotonicity-preserving ($\mathcal{CM}$-preserving) schemes (which includes the L1 and Grünwald-Letnikov schemes) on uniform meshes $\{t_n:=nh\}_{n=0}^\infty$, the discrete solution also has $O(t_{n}^{-α/γ})$ decay as $t_{n}\to\infty$. This result is then extended to $\mathcal{CM}$-preserving discretisations of certain time-fractional nonlinear subdiffusion problems such as the time-fractional porous media and $p$-Laplace equations. For the L1 scheme, the $O(t_{n}^{-α/γ})$ decay result is shown to remain valid on a very general class of nonuniform meshes. Our analysis uses a discrete comparison principle with discrete subsolutions and supersolutions that are carefully constructed to give tight bounds on the discrete solution. Numerical experiments are provided to confirm our theoretical analysis.

preprint2022arXiv

Supercloseness of the local discontinuous Galerkin method for a singularly perturbed convection-diffusion problem

A singularly perturbed convection-diffusion problem posed on the unit square in $\mathbb{R}^2$, whose solution has exponential boundary layers, is solved numerically using the local discontinuous Galerkin (LDG) method with piecewise polynomials of degree at most $k>0$ on three families of layer-adapted meshes: Shishkin-type, Bakhvalov-Shishkin-type and Bakhvalov-type.On Shishkin-type meshes this method is known to be no greater than $O(N^{-(k+1/2)})$ accurate in the energy norm induced by the bilinear form of the weak formulation, where $N$ mesh intervals are used in each coordinate direction. (Note: all bounds in this abstract are uniform in the singular perturbation parameter and neglect logarithmic factors that will appear in our detailed analysis.) A delicate argument is used in this paper to establish $O(N^{-(k+1)})$ energy-norm superconvergence on all three types of mesh for the difference between the LDG solution and a local Gauss-Radau projection of the exact solution into the finite element space. This supercloseness property implies a new $N^{-(k+1)}$ bound for the $L^2$ error between the LDG solution on each type of mesh and the exact solution of the problem; this bound is optimal (up to logarithmic factors). Numerical experiments confirm our theoretical results.

preprint2020arXiv

Why fractional derivatives with nonsingular kernels should not be used

In recent years, many papers discuss the theory and applications of new fractional-order derivatives that are constructed by replacing the singular kernel of the Caputo or Riemann-Liouville derivative by a non-singular (i.e., bounded) kernel. It will be shown here, through rigorous mathematical reasoning, that these non-singular kernel derivatives suffer from several drawbacks which should forbid their use. They fail to satisfy the fundamental theorem of fractional calculus since they do not admit the existence of a corresponding convolution integral of which the derivative is the left-inverse; and the value of the derivative at the initial time $t=0$ is always zero, which imposes an unnatural restriction on the differential equations and models where these derivatives can be used. For the particular cases of the so-called Caputo-Fabrizio and Atangana-Baleanu derivatives, it is shown that when this restriction holds the derivative can be simply expressed in terms of integer derivatives and standard Caputo fractional derivatives, thus demonstrating that these derivatives contain nothing new.

preprint2019arXiv

Existence, uniqueness and regularity of the solution of the time-fractional Fokker-Planck equation with general forcing

A time-fractional Fokker-Planck initial-boundary value problem is considered, with differential operator $u_t-\nabla\cdot(\partial_t^{1-α}κ_α\nabla u-\textbf{F}\partial_t^{1-α}u)$, where $0<α<1$. The forcing function $\textbf{F} = \textbf{F}(t,x)$, which is more difficult to analyse than the case $\textbf{F}=\textbf{F}(x)$ investigated previously by other authors. The spatial domain $Ω\subset\mathbb{R}^d$, where $d\ge 1$, has a smooth boundary. Existence, uniqueness and regularity of a mild solution $u$ is proved under the hypothesis that the initial data $u_0$ lies in $L^2(Ω)$. For $1/2<α<1$ and $u_0\in H^2(Ω)\cap H_0^1(Ω)$, it is shown that $u$ becomes a classical solution of the problem. Estimates of time derivatives of the classical solution are derived---these are known to be needed in numerical analyses of this problem.

preprint2016arXiv

Too much regularity may force too much uniqueness

Time-dependent fractional-derivative problems $D_t^δu + Au = f$ are considered, where $D_t^δ$ is a Caputo fractional derivative of order $δ\in (0,1)\cup (1,2)$ and~$A$ is a classical elliptic operator, and appropriate boundary and initial conditions are applied. The regularity of solutions to this class of problems is discussed, and it is shown that assuming more regularity than is generally true---as many researchers do---places a surprisingly severe restriction on the problem.

preprint2013arXiv

A finite difference method for a two-point boundary value problem with a Caputo fractional derivative

A two-point boundary value problem whose highest-order term is a Caputo fractional derivative of order $δ\in (1,2)$ is considered. Al-Refai's comparison principle is improved and modified to fit our problem. Sharp a priori bounds on derivatives of the solution $u$ of the boundary value problem are established, showing that $u''(x)$ may be unbounded at the interval endpoint $x=0$. These bounds and a discrete comparison principle are used to prove pointwise convergence of a finite difference method for the problem, where the convective term is discretized using simple upwinding to yield stability on coarse meshes for all values of $δ$. Numerical results are presented to illustrate the performance of the method.

preprint2013arXiv

Numerical methods for convection-diffusion problems or The 30 years war

Convection-diffusion problems arise in the modelling of many physical processes. Their typical solutions exhibit boundary and/or interior layers. Despite the linear nature of the differential operator, these problems pose still-unanswered questions to the numerical analyst. This talk will give a selective overview of numerical methods for the solution of convection-diffusion problems, while placing them in a historical context. It examines the principles that underpin the competing numerical techniques in this area and presents some recent developments.

preprint2013arXiv

Perturbed asymptotic expansions for interior-layer solutions of a semilinear reaction-diffusion problem with small diffusion

A semilinear reaction-diffusion two-point boundary value problem, whose second-order derivative is multiplied by a small positive parameter $\eps^2$, is considered. It can have multiple solutions. An asymptotic expansion is constructed for a solution that has an interior layer. Further properties are then established for a perturbation of this expansion. These are used in\cite{KoStMain} to obtain discrete sub-solutions and super-solutions for certain finite difference methods described there, and in this way yield convergence results for those methods.