Researcher profile

Martin Stynes

Martin Stynes contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 19 - UnverifiedVerification L1Unclaimed author
5works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

5 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.