Related papers: Excision methods for high resolution shock capturi…
In this paper, a positivity-preserving fifth-order finite volume compact-WENO scheme is proposed for solving compressible Euler equations. As we know conservative compact finite volume schemes have high resolution properties while WENO…
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…
In this paper we present a full general relativistic one-dimensional hydro-code which incorporates a modern high-resolution shock-capturing algorithm, with an approximate Riemann solver, for the correct modelling of formation and…
This paper presents high-order numerical methods for solving boundary value problems associated with the Lane-Emden equation, which frequently arises in astrophysics and various nonlinear models. A major challenge in studying this equation…
We describe a high-order ADER-DG solver for the compressible Euler equations within the ExaHyPE framework. The implementation combines a high-order ADER-DG polynomial representation, a local space-time DG predictor, adaptive mesh…
Various methods of treating outer boundaries in numerical relativity are compared using a simple test problem: a Schwarzschild black hole with an outgoing gravitational wave perturbation. Numerical solutions computed using different…
In this paper, the general procedure to solve the General Relativistic Hydrodynamical(GRH) equations with Adaptive-Mesh Refinement (AMR) is presented. In order to achieve, the GRH equations are written in the conservation form to exploit…
We propose an open-boundary molecular dynamics method in which an atomistic system is in contact with an infinite particle reservoir at constant temperature, volume and chemical potential. In practice, following the Hamiltonian adaptive…
We construct a new Godunov type relativistic hydrodynamics code in Milne coordinates, using a Riemann solver based on the two-shock approximation which is stable under the existence of large shock waves. We check the correctness of the…
In the case of hyperbolic conservation laws, high-order methods, such as the classical DG method, experience the phenomenon of unwanted high-frequency oscillations in the vicinity of a shock. Shock-capturing methods such as artificial…
In this work, we provide a novel method to constrain the causal parameter space of a relativistic hydrodynamic system exclusively from its linear stability analysis at non-zero momenta. Our approach exploits the Lorentz-invariant stability…
Numerical schemes for the general relativistic hydrodynamic equations are discussed. The use of conservative algorithms based upon the characteristic structure of those equations, developed during the last decade building on ideas first…
Many of the recent numerical simulations of binary black holes in vacuum adopt the moving puncture approach. This successful approach avoids the need to impose numerical excision of the black hole interior and is easy to implement. Here we…
We introduce novel high order well-balanced finite volume methods for the full compressible Euler system with gravity source term. They require no a priori knowledge of the hydrostatic solution which is to be well-balanced and are not…
We present a new pseudo-spectral code for the simulation of evolution systems that are second order in space. We test this code by evolving a non-linear scalar wave equation. These non-linear waves can be stably evolved using very simple…
The reconstruction of water wave elevation from bottom pressure measurements is an important issue for coastal applications, but corresponds to a difficult mathematical problem. In this paper we present the derivation of a method which…
We present a set of inner boundary conditions for the numerical construction of dynamical black hole space-times, when employing a 3+1 constrained evolution scheme and an excision technique. These inner boundary conditions are heuristically…
This manuscript presents an adaptive high order discretization technique for elliptic boundary value problems. The technique is applied to an updated version of the Hierarchical Poincar\'e-Steklov (HPS) method. Roughly speaking, the HPS…
Many astrophysical systems can only be accurately modelled when the behaviour of their baryonic gas components is well understood. The residual distribution (RD) family of partial differential equation (PDE) solvers produce approximate…
We present a splitting method for the one-dimensional Saint-Venant-Exner equations used for describing the bed evolution in shallow water systems. We adapt the flux vector splitting approach of Toro and Vazquez-Cend\`on and identify one…