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

Matthew F. Causley

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

ResearcherUnclaimed source record
math.NA

Catalog footprint

What is connected

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

3 published item(s)

preprint2016arXiv

Method of lines transpose: High order L-stable O(N) schemes for parabolic equations using successive convolution

We present a new solver for nonlinear parabolic problems that is L-stable and achieves high order accuracy in space and time. The solver is built by first constructing a single-dimensional heat equation solver that uses fast O(N) convolution. This fundamental solver has arbitrary order of accuracy in space, and is based on the use of the Green's function to invert a modified Helmholtz equation. Higher orders of accuracy in time are then constructed through a novel technique known as successive convolution (or resolvent expansions). These resolvent expansions facilitate our proofs of stability and convergence, and permit us to construct schemes that have provable stiff decay. The multi-dimensional solver is built by repeated application of dimensionally split independent fundamental solvers. Finally, we solve nonlinear parabolic problems by using the integrating factor method, where we apply the basic scheme to invert linear terms (that look like a heat equation), and make use of Hermite-Birkhoff interpolants to integrate the remaining nonlinear terms. Our solver is applied to several linear and nonlinear equations including heat, Allen-Cahn, and the Fitzhugh-Nagumo system of equations in one and two dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Higher order A-stable schemes for the wave equation using a recursive convolution approach

In several recent works \cite{Causley2013a}, \cite{Causley2013}, we developed a new second order, A-stable approach to wave propagation problems based on the method of lines transpose (MOL$^T$) formulation combined with alternating direction implicit (ADI) schemes. In this work, we present several important modifications to our work, and thus obtain a family of wave solvers which are unconditionally stable, accurate of order 2P, and require $O(P^d N)$ operations per time step, where $N$ is the number of spatial points, and $d$ the number of spatial dimensions. We obtain these schemes by including higher derivatives of the solution, rather than increasing the number of time levels. The novel aspect of our approach is that the higher derivatives are constructed using successive applications of the convolution operator. We develop these schemes in one spatial dimension, and then extend the results to higher dimensions, by reformulating the ADI scheme to include recursive convolution. Thus, we retain a fast, unconditionally stable scheme, which does not suffer from the large dispersion errors characteristic to the ADI method. We demonstrate the utility of the method by applying it to a host of wave propagation problems. This method holds great promise for developing higher order, parallelizable algorithms for solving hyperbolic PDEs, and can also be extended to parabolic PDEs.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Method of Lines Transpose: A Fast Implicit Wave Propagator

As a follow up to \cite{Causley2013}, we provide a detailed description of the numerical implementation of an O(N), A-stable, second order accurate solution of the wave equation, constructed from semi-discrete boundary value problems. We improve on the previous algorithm by replacing the Lax-type correction used in \cite{Causley2013}, which was necessary for convergence when $Δt < Δx/c$, with a more accurate spatial quadrature, which we prove is convergent. We also demonstrate that the resulting solver remains fast even in the case of unstructured meshes, can incorporate domain decomposition, and allows for the implementation of Dirichlet, Neumann, periodic and outflow boundary conditions. Building upon results for the 1d formulation, we utilize alternate direction implicit (ADI) splitting to achieve a fast O(N) solver in higher spatial dimensions. Our solver is built upon line objects and, combined with the flexibility of the integral solver, allows us to solve problems on arbitrary spatial domains, by embedding the boundary in a regular Cartesian mesh. Our solver is designed to couple with particle codes, where scale separation is an issue. We therefore demonstrate the ability of our solver to take time steps well beyond that of the Courant-Friedrichs-Lewy (CFL) stability limit of explicit codes.

0 citations0 reviewsNo code linkedOpen access