相关论文: General relativistic null-cone evolutions with a h…
It has now become customary in the field of numerical relativity to couple high order finite difference schemes to mesh refinement algorithms. To this end, different modifications to the standard Berger-Oliger adaptive mesh refinement…
Gridless methods show great superiority in line spectral estimation. These methods need to solve an atomic $l_0$ norm (i.e., the continuous analog of $l_0$ norm) minimization problem to estimate frequencies and model order. Since this…
This paper is concerned with the construction of high order schemes on irregular grids for balance laws, including a discussion of an a-posteriori error indicator based on the numerical entropy production. We also impose well-balancing on…
We present compact semi-implicit finite difference schemes on structured grids for numerical solutions of the advection by an external velocity and by a speed in normal direction that are applicable in level set methods. The most involved…
We introduce a nonlinear structure preserving high-order scheme for anisotropic advection-diffusion equations. This scheme, based on Hybrid High-Order methods, can handle general meshes. It also has an entropy structure, and preserves the…
We use the canonical formalism developed together with David Robinson to st= udy the Einstein equations on a null surface. Coordinate and gauge conditions = are introduced to fix the triad and the coordinates on the null surface. Toget= her…
In this study, we explore multiple higher-order pole solutions in spinor Bose--Einstein condensates. These solutions are associated with different pairs of higher-order poles of the transmission coefficient in the inverse scattering…
A high-order accurate implicit-mesh discontinuous Galerkin framework for wave propagation in single-phase and bi-phase solids is presented. The framework belongs to the embedded-boundary techniques and its novelty regards the spatial…
A high-order accurate adjoint-based optimization framework is presented for unsteady multiphysics problems. The fully discrete adjoint solver relies on the high-order, linearly stable, partitioned solver introduced in [1], where different…
This is the first of a series of papers describing a numerical implementation of the conformally rescaled Einstein equation, an implementation designed to calculate asymptotically flat spacetimes, especially spacetimes containing black…
Exponential integrators based on contour integral representations lead to powerful numerical solvers for a variety of ODEs, PDEs, and other time-evolution equations. They are embarrassingly parallelizable and lead to global-in-time…
In this paper, a novel high-order, mass and energy-conserving scheme is proposed for the regularized logarithmic Schr\"{o}dinger equation(RLogSE). Based on the idea of the supplementary variable method (SVM), we firstly reformulate the…
We focus here on a class of fourth-order parabolic equations that can be written as a system of second-order equations by introducing an auxiliary variable. We design a novel second-order fully discrete mixed finite element method to…
We present a high-order method that provides numerical integration on volumes, surfaces, and lines defined implicitly by two smooth intersecting level sets. To approximate the integrals, the method maps quadrature rules defined on…
We study higher form Proca equations on Einstein manifolds with boundary data along conformal infinity. We solve these Laplace-type boundary problems formally, and to all orders, by constructing an operator which projects arbitrary forms to…
We present and compare third- as well as fifth-order accurate finite difference schemes for the numerical solution of the compressible ideal MHD equations in multiple spatial dimensions. The selected methods lean on four different…
We develop a high-order, explicit method for acoustic scattering in three space dimensions based on a combined-field time-domain integral equation. The spatial discretization, of Nystr\"om type, uses Gaussian quadrature on panels combined…
We present a fully grid-based approach for solving Hartree-Fock and all-electron Kohn-Sham equations based on low-rank approximation of three-dimensional electron orbitals. Due to the low-rank structure the total complexity of the algorithm…
We present a new computational framework (LEO), that enables us to carry out the very first large-scale, high-resolution computations in the context of the characteristic approach in numerical relativity. At the analytic level, our approach…
In this paper, we consider a non-linear fourth-order evolution equation of Cahn-Hilliard-type on evolving surfaces with prescribed velocity, where the non-linear terms are only assumed to have locally Lipschitz derivatives. High-order…