BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

William McLean

William McLean contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NAmath.APNumerical Analysismath.PR

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
4topics
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

9 published item(s)

preprint2022arXiv

Numerical Evaluation of Mittag-Leffler Functions

The Mittag-Leffler function is computed via a quadrature approximation of a contour integral representation. We compare results for parabolic and hyperbolic contours, and give special attention to evaluation on the real line. The main point of difference with respect to similar approaches from the literature is the way that poles in the integrand are handled. Rational approximation of the Mittag-Leffler function on the negative real axis is also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Implementation of high-order, discontinuous Galerkin time stepping for fractional diffusion problems

The discontinuous Galerkin dG method provides a robust and flexible technique for the time integration of fractional diffusion problems. However, a practical implementation uses coefficients defined by integrals that are not easily evaluated. We describe specialised quadrature techniques that efficiently maintain the overall accuracy of the dG method. In addition, we observe in numerical experiments that known superconvergence properties of dG time stepping for classical diffusion problems carry over in a modified form to the fractional-order setting.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Uniform stability for a spatially-discrete, subdiffusive Fokker-Planck equation

We prove stability estimates for the spatially discrete, Galerkin solution of a fractional Fokker-Planck equation, improving on previous results in several respects. Our main goal is to establish that the stability constants are bounded uniformly in the fractional diffusion exponent $α\in(0,1]$. In addition, we account for the presence of an inhomogeneous term and show a stability estimate for the gradient of the Galerkin solution. As a by-product, the proofs of error bounds for a standard finite element approximation are simplified.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A second-order scheme with nonuniform time steps for a linear reaction-sudiffusion problem

Stability and convergence of a time-weighted discrete scheme with nonuniform time steps are established for linear reaction-subdiffusion equations. The Caupto derivative is approximated at an offset point by using linear and quadratic polynomial interpolation. Our analysis relies on two tools: a discrete fractional Grönwall inequality and the global consistency analysis. The new consistency analysis makes use of an interpolation error formula for quadratic polynomials, which leads to a convolution-type bound for the local truncation error. To exploit these two tools, some theoretical properties of the discrete kernels in the numerical Caputo formula are crucial and we investigate them intensively in the nonuniform setting. Taking the initial singularity of the solution into account, we obtain a sharp error estimate on nonuniform time meshes. The fully discrete scheme generates a second-order accurate solution on the graded mesh provided a proper grading parameter is employed. An example is presented to show the sharpness of our analysis.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Regularity theory for time-fractional advection-diffusion-reaction equations

We investigate the behavior of the time derivatives of the solution to a linear time-fractional, advection-diffusion-reaction equation, allowing space- and time-dependent coefficients as well as initial data that may have low regularity. Our focus is on proving estimates that are needed for the error analysis of numerical methods. The nonlocal nature of the fractional derivative creates substantial difficulties compared with the case of a classical parabolic PDE. In our analysis, we rely on novel energy methods in combination with a fractional Gronwall inequality and certain properties of fractional integrals.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Well-posedness of time-fractional, advection-diffusion-reaction equations

We establish the well-posedness of an initial-boundary value problem for a general class of time-fractional, advection-diffusion-reaction equations, allowing space- and time-dependent coefficients as well as initial data that may have low regularity. Our analysis relies on novel energy methods in combination with a fractional Gronwall inequality and properties of fractional integrals.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

Wider contours and adaptive contours

Contour integrals in the complex plane are the basis of effective numerical methods for computing matrix functions, such as the matrix exponential and the Mittag-Leffler function. These methods provide successful ways to solve partial differential equations, such as convection--diffusion models. Part of the success of these methods comes from exploiting the freedom to choose the contour, by appealing to Cauchy&#39;s theorem. However, the pseudospectra of non-normal matrices or operators present a challenge for these methods: if the contour is too close to regions where the norm of the resolvent matrix is large, then the accuracy suffers. Important applications that involve non-normal matrices or operators include the Black--Scholes equation of finance, and Fokker--Planck equations for stochastic models arising in biology. Consequently, it is crucial to choose the contour carefully. As a remedy, we discuss choosing a contour that is wider than it might otherwise have been for a normal matrix or operator. We also suggest a semi-analytic approach to adapting the contour, in the form of a parabolic bound that is derived by estimating the field of values. To demonstrate the utility of the approaches that we advocate, we study three models in biology: a monomolecular reaction, a bimolecular reaction and a trimolecular reaction. Modelling and simulation of these reactions is done within the framework of Markov processes. We also consider non-Markov generalisations that have Mittag-Leffler waiting times instead of the usual exponential waiting times of a Markov process.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Time-stepping error bounds for fractional diffusion problems with non-smooth initial data

We apply the piecewise constant, discontinuous Galerkin method to discretize a fractional diffusion equation with respect to time. Using Laplace transform techniques, we show that the method is first order accurate at the \$n\$th time level \$t_n\$, but the error bound includes a factor \$t_n^{-1}\$ if we assume no smoothness of the initial data. We also show that for smoother initial data the growth in the error bound as \$t_n\$ decreases is milder, and in some cases absent altogether. Our error bounds generalize known results for the classical heat equation and are illustrated for a model problem.

0 citations0 reviewsNo code linkedOpen access