Related papers: General relativistic null-cone evolutions with a h…
Gradient schemes is a framework that enables the unified convergence analysis of many numerical methods for elliptic and parabolic partial differential equations: conforming and non-conforming Finite Element, Mixed Finite Element and Finite…
We revisit Winicour's affine-null metric initial value formulation of General Relativity, where the characteristic initial value formulation is set up with a null metric having two affine parameters. In comparison to past work, where the…
We present novel model reduction methods for rapid solution of parametrized nonlinear partial differential equations (PDEs) in real-time or many-query contexts. Our approach combines reduced basis (RB) space for rapidly convergent…
Although high-order Maxwell integral equation solvers provide significant advantages in terms of speed and accuracy over corresponding low-order integral methods, their performance significantly degrades in presence of non-smooth…
In this paper new innovative fourth order compact schemes for Robin and Neumann boundary conditions have been developed for boundary value problems of elliptic PDEs in two and three dimensions. Different from traditional finite difference…
The Lichnerowicz and Israel theorems are extended to higher order theories of gravity. In particular it is shown that Schwarzschild is the unique spherically symmetric, static, asymptotically flat, black-hole solution, provided the spatial…
In this paper, we design and analyze a Hybrid-High Order (HHO) approximation for a class of quasilinear elliptic problems of nonmonotone type. The proposed method has several advantages, for instance, it supports arbitrary order of…
In this paper, we propose high order numerical methods to solve a 2D advection diffusion equation, in the highly oscillatory regime. We use an integrator strategy that allows the construction of arbitrary high-order schemes {leading} to an…
In previous work, the numerical solution of the linearized gravitational field equations near space-like and null-infinity was discussed in the form of the spin-2 zero-rest-mass equation for the perturbations of the conformal Weyl…
We describe in this article a new code for evolving axisymmetric isolated systems in general relativity. Such systems are described by asymptotically flat space-times which have the property that they admit a conformal extension. We are…
This study proposes an intrusive projection-based model-order reduction framework for nonlinear problems with a polynomial structure, solved iteratively using a Newton solver in the reduced space. It is demonstrated that for the targeted…
A high-order Newton multigrid method is proposed for steady-state shallow water flows in open channels with regular and irregular geometries. The method integrates a finite volume discretization with third-order weighted essentially…
In this paper, we present an interior penalty discontinuous Galerkin finite element scheme for solving diffusion problems with strong anisotropy arising in magnetized plasmas for fusion applications. We demonstrate the accuracy produced by…
The space nonlocal Allen-Cahn equation is a famous example of fractional reaction-diffusion equations. It is also an extension of the classical Allen-Cahn equation, which is widely used in physics to describe the phenomenon of two-phase…
Higher-order exceptional points in the spectrum of non-Hermitian Hamiltonians describing open quantum or wave systems have a variety of potential applications in particular in optics and photonics. However, the experimental realization is…
Motivated by the problem of solving the Einstein equations, we discuss high order finite difference discretizations of first order in time, second order in space hyperbolic systems.Particular attention is paid to the case when first order…
We present simulations of the Einstein-Maxwell-Klein-Gordon system on compactified hyperboloidal slices. To the best of our knowledge, this is the first time that this setup is evolved with a common formulation like BSSN/Z4. Hyperboloidal…
The hybrid high-order method is a modern numerical framework for the approximation of elliptic PDEs. We present here an extension of the hybrid high-order method to meshes possessing curved edges/faces. Such an extension allows us to…
We introduce two multiscale numerical schemes for the time integration of weakly nonlinear Schr\"odinger equations, built upon the discretization of Picard iterates of the solution. These high-order schemes are designed to achieve high…
Motivated by the need for efficient and accurate simulation of the dynamics of the polar ice sheets, we design high-order finite element discretizations and scalable solvers for the solution of nonlinear incompressible Stokes equations. We…