We develop an optimization-based approach for additive decomposition and reconnection of algebraic problems arising from discretization of partial differential equations (PDEs). Application to a scalar advection-diffusion PDE illustrates the new approach. In particular, we obtain a robust iterative solver for advection-dominated problems using standard multi-level solvers for the Poisson equation.

In this paper we present a method to treat interface jump conditions for constant coefficients Poisson problems that allows the use of standard “black box” solvers, without compromising accuracy. The basic idea of the new approach is similar to the Ghost Fluid Method (GFM). The GFM relies on corrections applied on nodes located across the interface for discretization stencils that straddle the ...

Our research in theoretical biophysics [6] recently required us to solve a variant of the Poisson equation, which is a linear elliptic partial differential equation (PDE) in 3D. Ultimately, we want to make our method available as a applet, so we decided to implement a PDE solver in Java. Another reason for using Java was its promise of easy, portable multithreaded programming, because support f...

We study confined solutions of certain evolutionary partial differential equations (PDE) in 1 + 1 space–time. The PDE we study are Lie–Poisson Hamiltonian systems for quadratic Hamiltonians defined on the dual of the Lie algebra of vector fields on the real line. These systems are also Euler–Poincaré equations for geodesic motion on the diffeomorphism group in the sense of the Arnold program fo...

The spine of an object is an entity that can characterise the object’s topology and describes the object by a lower dimension. It has an intuitive appeal for supporting geometric modelling operations. The aim of this paper is to show how a spine for a PDE surface can be generated. For the purpose of the work presented here an analytic solution form for the chosen PDE is utilised. It is shown th...

Abstract Infinitesimal symmetries of a partial differential equation (PDE) can be defined algebraically as the solutions linearisation (Frechet derivative) holding on space to PDE, and they are well-known comprise linear having structure Lie algebra. Solutions adjoint PDE called adjoint-symmetries. Their algebraic for general systems is studied herein. This motivated by correspondence between v...

In this paper, we deal with the problem of computing the distance to a surface (a curve in 2D) and consider several distance function approximation methods which are based on solving partial differential equations (PDEs) and finding solutions to variational problems. In particular, we deal with distance function estimation methods related to the Poisson-like equations and generalized double-lay...

In this paper, we consider a PDE system arising in corrosion modelling. This system consists in two convection-diffusion equations on the densities of charge carriers and a Poisson equation on the electric potential. Boundary conditions are Robin boundary conditions. We discretize each equation by a finite volume scheme and we prove the convergence of the scheme towards a weak solution to the i...