Related papers: General relativistic null-cone evolutions with a h…
We introduce a new class of "filtered" schemes for some first order non-linear Hamilton-Jacobi-Bellman equations. The work follows recent ideas of Froese and Oberman (SIAM J. Numer. Anal., Vol 51, pp.423-444, 2013). The proposed schemes are…
A high-order quadrature algorithm is presented for computing integrals over curved surfaces and volumes whose geometry is implicitly defined by the level sets of (one or more) multivariate polynomials. The algorithm recasts the implicitly…
An approximation to General Relativity is presented that agrees with the Einstein field equations up to and including the fourth post-Newtonian (PN) order. This approximation is formulated in a fully constrained scheme: all involved…
A discretization scheme for variable coefficient elliptic PDEs in the plane is presented. The scheme is based on high-order Gaussian quadratures and is designed for problems with smooth solutions, such as scattering problems involving soft…
We present a higher order space-time unfitted finite element method for convection-diffusion problems on coupled (surface and bulk) domains. In that way, we combine a method suggested by Heimann, Lehrenfeld, Preu{\ss} (SIAM J. Sci. Comput.…
In the harmonic description of general relativity, the principle part of Einstein's equations reduces to 10 curved space wave equations for the componenets of the space-time metric. We present theorems regarding the stability of several…
The Generalized Riemann Problems (GRP) for nonlinear hyperbolic systems of balance laws in one space dimension are now well-known and can be formulated as follows: Given initial-data which are smooth on two sides of a discontinuity,…
This paper presents an efficient high-order sharp-interface method for solving the three-dimensional (3D) Poisson equation with Dirichlet boundary conditions on a nonuniform Cartesian grid with irregular domain boundaries. The new approach…
This paper considers the approximation of partial differential equations with a point collocation framework based on high-order local maximum-entropy schemes (HOLMES). In this approach, smooth basis functions are computed through an…
A straightforward algorithm for the symbolic computation of higher-order symmetries of nonlinear evolution equations and lattice equations is presented. The scaling properties of the evolution or lattice equations are used to determine the…
We present a new framework for statistical inference on Riemannian manifolds that achieves high-order accuracy, addressing the challenges posed by non-Euclidean parameter spaces frequently encountered in modern data science. Our approach…
Second-order superintegrable systems in dimensions two and three are essentially classified. With increasing dimension, however, the non-linear partial differential equations employed in current methods become unmanageable. Here we propose…
We study four-dimensional Einstein-Maxwell fields for which any higher-order corrections to the field equations effectively reduces to just a rescaling of the gravitational and the cosmological constant. These configurations are thus…
We linearize the Einstein equations when the metric is Bondi-Sachs, when the background is Schwarzschild or Minkowski, and when there is a matter source in the form of a thin shell whose density varies with time and angular position. By…
Perturbations of the linearized vacuum Einstein equations on a null cone in the Bondi-Sachs formulation of General Relativity can be derived from a single master function with spin weight two, which is related to the Weyl scalar \Psi_0, and…
A fourth-order finite volume embedded boundary (EB) method is presented for the unsteady Stokes equations. The algorithm represents complex geometries on a Cartesian grid using EB, employing a technique to mitigate the "small cut-cell"…
New implicit and implicit-explicit time-stepping methods for the wave equation in second-order form are described with application to two and three-dimensional problems discretized on overset grids. The implicit schemes are single step,…
This paper considers the problem for finding the $(\delta,\epsilon)$-Goldstein stationary point of Lipschitz continuous objective, which is a rich function class to cover a great number of important applications. We construct a zeroth-order…
We show the stability of the geometric optics approximation in general relativity by constructing a family $(g_\lambda)_{\lambda\in(0,1]}$ of high-frequency metrics solutions to the Einstein vacuum equations in 3+1 dimensions without any…
A third-order weighted essentially non-oscillatory compact least-squares scheme is developed for the finite volume method on structured curvilinear non-uniform grids. The proposed scheme features compact least-squares reconstruction with…