A fast and accurate algorithm for the reinitialization of the signed distance function in two and three spatial dimensions is presented. The algorithm has computational complexity O(N logN) for the reinitialization of N grid points. The order of accuracy of the reinitialization is demonstrated to depend primarily on the interpolation algorithm used. Bicubic interpolation is demonstrated to resu...

We introduce a framework for the design of finite element methods for two-dimensional moving boundary problems with prescribed boundary evolution that have arbitrarily high order of accuracy, both in space and in time. At the core of our approach is the use of a universal mesh: a stationary background mesh containing the domain of interest for all times that adapts to the geometry of the immers...

In this paper we carry out the extension of the ADER approach to multidimensional non-linear systems of conservation laws. We implement non-linear schemes of up to fourth order of accuracy in both time and space. Numerical results for the compressible Euler equations illustrate the very high order of accuracy and non-oscillatory properties of the new schemes. Compared to the state-of-art finite...

We extend the technique of optimal time step (OTS) selection for finite difference (FD) schemes of time dependent PDEs to PDEs where the leading-order spatial derivative term has a spatially varying coefficient. The basic approach involves identifying a transformation of the domain that eliminates the spatial dependence of the coefficient for the leading-order term. This change of variables is ...

The commonly used one step methods and linear multi-step methods all have a global error that is of the same order as the local truncation error. In fact, this is true of the entire class of general linear methods. In practice, this means that the order of the method is typically defined solely by order conditions which are derived by studying the local truncation error. In this work we investi...

Spectral deferred corrections (SDC) is an iterative approach for constructing higherorder accurate numerical approximations of ordinary differential equations. SDC starts with an initial approximation of the solution defined at a set of Gaussian or spectral collocation nodes over a time interval and uses an iterative application of lower-order time discretizations applied to a correction equati...

In this paper, we analyze a multiscale method developed under the heterogeneous multiscale method (HMM) framework for numerical approximation of multiscale wave propagation problems in periodic media. In particular, we are interested in the long time O(ε−2) wave propagation, where ε represents the size of the microscopic variations in the media. In large time scales, the solutions of multiscale...

We give a further examination of the stabilized Residual Distribution schemes for the solution of the shallow water equations proposed in (Ricchiuto and Bollermann, J.Comp.Phys., available online 25 October 2008). Based on a non-linear variant of a Lax-Friedrichs scheme, the scheme is wellbalanced, able to handle dry areas and, for smooth regions of the solution, obtains second order of accurac...

A coupled system of two singularly perturbed linear reaction–diffusion two-point boundary value problems is examined. The leading term of each equation is multiplied by a small positive parameter, but these parameters may have different magnitudes. The solutions to the system have boundary layers that overlap and interact. The structure of these layers is analysed, and this leads to the constru...

The order of accuracy of the node-centered finite volume methods is analyzed, and the analysis is based on an exact derivation of the numerical errors in one dimension. The accuracy for various types of grids are considered. Numerical simulations and analysis are performed for both a hyperbolic and a eliptic case, and the results agree. The impact of weakly imposed boundary conditions is analyz...