Related papers: A fourth-order accurate compact difference scheme …
We present a method that gives highly accurate electrostatic potentials for systems where we have periodic boundary conditions in two spatial directions but free boundary conditions in the third direction. These boundary conditions are…
As part of our development of a computer code to perform 3D `constrained evolution' of Einstein's equations in 3+1 form, we discuss issues regarding the efficient solution of elliptic equations on domains containing holes (i.e., excised…
We propose a stabilized Nitsche-based cut finite element formulation for the Oseen problem in which the boundary of the domain is allowed to cut through the elements of an easy-to-generate background mesh. Our formulation is based on the…
A novel fourth-order finite difference formula coupling the Crank-Nicolson explicit linearized method is proposed to solve Riesz space fractional nonlinear reaction-diffusion equations in two dimensions. Theoretically, under the Lipschitz…
The Poisson-Boltzmann equation offers an efficient way to study electrostatics in molecular settings. Its numerical solution with the boundary element method is widely used, as the complicated molecular surface is accurately represented by…
In this paper, a novel immersed boundary method is developed, validated, and applied. Through devising a second-order three-step flow reconstruction scheme, the proposed method is able to enforce the Dirichlet, Neumann, Robin, and Cauchy…
This study aims to construct a stable, high-order compact finite difference method for solving Sobolev-type equations with Dirichlet boundary conditions in one-space dimension. Approximation of higher-order mixed derivatives in some…
In this paper we develop finite difference schemes for elliptic problems with piecewise continuous coefficients that have (possibly huge) jumps across fixed internal interfaces. In contrast with such problems involving one smooth…
This paper proposes and analyzes a finite difference method based on compact schemes for the Euler-Bernoulli beam equation with damping terms. The method achieves fourth-order accuracy in space and second-order accuracy in time, while…
In this paper, we study the biharmonic equation with Dirichlet boundary conditions in a polygonal domain. In particular, we propose a method that effectively decouples the fourth-order problem into a system of two Poison equations and one…
For elliptic interface problems with discontinuous coefficients, the maximum accuracy order for compact 9-point finite difference scheme in irregular points is three [7]. The discontinuous coefficients usually have abrupt jumps across the…
We introduce a refined immersed boundary (IB) methodology that is better-than-first-order accurate in practice, while preserving key properties of "continuous-forcing" IB approaches that retain a singular source term in the governing…
We present a polynomial multigrid method for nodal interior penalty and local discontinuous Galerkin formulations of the Poisson equation on Cartesian grids. For smoothing we propose two classes of overlapping Schwarz methods. The first…
The elliptic interface problems with discontinuous and high-contrast coefficients appear in many applications and often lead to huge condition numbers of the corresponding linear systems. Thus, it is highly desired to construct high order…
We introduce a technique that simplifies the problem of imposing jump conditions on interfaces that are not aligned with a computational grid in the context of the Correction Function Method (CFM). The CFM offers a general framework to…
A solver for the Poisson equation for 1D, 2D and 3D regular grids is presented. The solver applies the convolution theorem in order to efficiently solve the Poisson equation in spectral space over a rectangular computational domain.…
In this article, we have developed a higher order compact numerical method for variable coefficient parabolic problems with mixed derivatives. The finite difference scheme, presented here for two-dimensional domains, is based on fourth…
An immersed-boundary method for the incompressible Navier--Stokes equations is presented. It employs discrete forcing for a sharp discrimination of the solid-fluid interface, and achieves second-order accuracy, demonstrated in examples with…
In this paper, we present a parallel higher-order boundary integral method to solve the linear Poisson-Boltzmann (PB) equation. In our method, a well-posed boundary integral formulation is used to ensure the fast convergence of Krylov…
We consider the numerical solution of Poisson's equation on structured grids using geometric multigrid with nonstandard coarse grids and coarse level operators. We are motivated by the problem of developing high-order accurate numerical…