Researcher profile

M. C. Pugh

M. C. Pugh contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 15 - UnverifiedVerification L1Unclaimed author
3works
0followers
7topics
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

3 published item(s)

preprint2020arXiv

A Study of the Numerical Stability of an ImEx Scheme with Application to the Poisson-Nernst-Planck Equations

The Poisson-Nernst-Planck equations with generalized Frumkin-Butler-Volmer boundary conditions (PNP-FBV) describe ion transport with Faradaic reactions and have applications in a wide variety of fields. Using an adaptive time-stepper based on a second-order variable step-size, semi-implicit, backward differentiation formula (VSSBDF2), we observe that when the underlying dynamics is one that would have the solutions converge to a steady state solution, the adaptive time-stepper produces solutions that "nearly" converge to the steady state and that, simultaneously, the time-step sizes stabilize at a limiting size $dt_\infty$. Linearizing the SBDF2 scheme about the steady state solution, we demonstrate that the linearized scheme is conditionally stable and that this is the cause of the adaptive time-stepper's behaviour. Mesh-refinement, as well as a study of the eigenvectors corresponding to the critical eigenvalues, demonstrate that the conditional stability is not due to a time-step constraint caused by high-frequency contributions. We study the stability domain of the linearized scheme and find that it can have corners as well as jump discontinuities.

preprint2020arXiv

Adaptive Time-stepping Schemes for the Solution of the Poisson-Nernst-Planck Equations

The Poisson-Nernst-Planck equations with generalized Frumkin-Butler-Volmer boundary conditions (PNP-FBV) describe ion transport with Faradaic reactions, and have applications in a number of fields. In this article, we develop an adaptive time-stepping scheme for the solution of the PNP-FBV equations based on two time-stepping methods: a fully implicit (BDF2) method, and an implicit-explicit (SBDF2) method. We present simulations under both current and voltage boundary conditions and demonstrate the ability to simulate a large range of parameters, including any value of the singular perturbation parameter $ε$. When the underlying dynamics is one that would have the solutions converge to a steady-state solution, we observe that the adaptive time-stepper based on the SBDF2 method produces solutions that ``nearly'' converge to the steady state and that, simultaneously, the time-step sizes stabilize to a limiting size $dt_\infty$. In the companion to this article \cite{YPD_Part2}, we linearize the SBDF2 scheme about the steady-state solution and demonstrate that the linearized scheme is conditionally stable. This conditional stability is the cause of the adaptive time-stepper's behaviour. While the adaptive time-stepper based on the fully-implicit (BDF2) method is not subject to such time-step constraints, the required nonlinear solve yields run times that are significantly longer.

preprint2010arXiv

Research Announcement: Finite-time Blow Up and Long-wave Unstable Thin Film Equations

We study short--time existence, long--time existence, finite speed of propagation, and finite--time blow--up of nonnegative solutions for long-wave unstable thin film equations $h_t = -a_0(h^n h_{xxx})_x - a_1(h^m h_x)_x$ with $n>0$, $a_0 > 0$, and $a_1 >0$. The existence and finite speed of propagation results extend those of [Comm Pure Appl Math 51:625--661, 1998]. For $0<n<2$ we prove the existence of a nonnegative, compactly--supported, strong solution on the line that blows up in finite time. The construction requires that the initial data be nonnegative, compactly supported in $\R^1$, be in $H^1(\R^1)$, and have negative energy. The blow-up is proven for a large range of $(n,m)$ exponents and extends the results of [Indiana Univ Math J 49:1323--1366, 2000].