Related papers: Efficient discretization of the Laplacian on compl…
We introduce a novel discretization of the Monge-Ampere operator, simultaneously consistent and degenerate elliptic, hence accurate and robust in applications. These properties are achieved by exploiting the arithmetic structure of the…
The construction of stable, conservative, and accurate volume dissipation is extended to discretizations that possess a generalized summation-by-parts (SBP) property within a tensor-product framework. The dissipation operators can be…
We study the two-dimensional surface quasi-geostrophic equation on a bounded domain with a smooth boundary. Motivated by the three-dimensional incompressible Navier-Stokes equations and previous results in the entire space $\mathbb R^2$, we…
Developments in dynamical systems theory provides new support for the discretisation of \pde{}s and other microscale systems. By systematically resolving subgrid microscale dynamics the new approach constructs asymptotically accurate,…
This note is about promoting singularity subtraction as a helpful tool in the discretization of singular integral operators on curved surfaces. Singular and nearly singular kernels are expanded in series whose terms are integrated on…
This work focuses on developing high-order energy-stable schemes for wave-dominated problems in closed domains using staggered finite-difference summation-by-parts (SBP FD) operators. We extend the previously presented uniform staggered…
Smoothed particle hydrodynamics (SPH) has been extensively used to model high and low Reynolds number flows, free surface flows and collapse of dams, study pore-scale flow and dispersion, elasticity, and thermal problems. In different…
We present a general purpose method for solving partial differential equations on a closed surface, based on a technique for discretizing the surface introduced by Wenjun Ying and Wei-Cheng Wang [J. Comput. Phys. 252 (2013), pp. 606-624]…
An efficient $hp$-multigrid scheme is presented for local discontinuous Galerkin (LDG) discretizations of elliptic problems, formulated around the idea of separately coarsening the underlying discrete gradient and divergence operators. We…
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 an MPC-protected distributed computation, although the use of MPC assures data privacy during computation, sensitive information may still be inferred by curious MPC participants from the computation output. This can be observed, for…
A quasi-second order scheme is developed to obtain approximate solutions of the shallow water equationswith bathymetry. The scheme is based on a staggered finite volume scheme for the space discretization:the scalar unknowns are located in…
We analyze two types of summation-by-parts finite difference operators for approximating the second derivative with variable coefficient. The first type uses ghost points, while the second type does not use any ghost points. A previously…
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…
In this paper, we consider finite difference approximations of the second order wave equation. We use finite difference operators satisfying the summation-by-parts property to discretize the equation in space. Boundary conditions and grid…
A new and efficient neural-network and finite-difference hybrid method is developed for solving Poisson equation in a regular domain with jump discontinuities on embedded irregular interfaces. Since the solution has low regularity across…
In the present paper we describe a class of algorithms for the solution of Laplace's equation on polygonal domains with Neumann boundary conditions. It is well known that in such cases the solutions have singularities near the corners which…
This paper establishes strong convergence rates for the spatial finite element discretization of a two-dimensional stochastic Navier--Stokes system with transport noise and no-slip boundary conditions on a convex polygonal domain. The main…
By employing non-equispaced grid points near boundaries, boundary-optimized upwind finite-difference operators of orders up to nine are developed. The boundary closures are constructed within a diagonal-norm summation-by-parts (SBP)…
The paper describes a sparse direct solver for the linear systems that arise from the discretization of an elliptic PDE on a two dimensional domain. The scheme decomposes the domain into thin subdomains, or ``slabs'' and uses a two-level…