Related papers: Numerical simulations with a first order BSSN form…
We extend the first order dissipative relativistic hydrodynamics model of Bemfica-Disconzi-Noronha- Kovtun (BDNK) in order to include the charge number current in full first order expansion with out-of-equilibrium contribution proportional…
In this paper we present a new first-order hyperbolic reformulation of the Cahn-Hilliard equation. The model is obtained from the combination of augmented Lagrangian techniques proposed earlier by the authors of this paper, with a classical…
We study a second order hyperbolic initial-boundary value partial differential equation with memory, that results in an integro-differential equation with a convolution kernel. The kernel is assumed to be either smooth or no worse than…
We develop a hybrid spatial discretization for the wave equation in second order form, based on high-order accurate finite difference methods and discontinuous Galerkin methods. The hybridization combines computational efficiency of finite…
This paper considers the discontinuous Galerkin (DG) methods for solving the Vlasov-Maxwell (VM) system, a fundamental model for collisionless magnetized plasma. The DG methods provide accurate numerical description with conservation and…
This paper is devoted to the study of Bingham flow with variable density. We propose a local bi-viscosity regularization of the stress tensor based on a Huber smoothing step. Next, our computational approach is based on a second-order,…
A new finite element method with discontinuous approximation is introduced for solving second order elliptic problem. Since this method combines the features of both conforming finite element method and discontinuous Galerkin (DG) method,…
Discontinuous Galerkin (DG) methods for solving elliptic equations are gaining popularity in the computational physics community for their high-order spectral convergence and their potential for parallelization on computing clusters.…
A novel explicit and implicit Kinetic Streamlined-Upwind Petrov Galerkin (KSUPG) scheme is presented for hyperbolic equations such as Burgers equation and compressible Euler equations. The proposed scheme performs better than the original…
We propose a Discontinuous Galerkin (DG) scheme for the numerical solution of the Hydrostatic Stokes equations in Oceanography. This new scheme is based on the introduction of the symmetric interior penalty (SIP) technique for the…
We propose a numerical method to solve general hyperbolic systems in any space dimension using forward Euler time stepping and continuous finite elements on non-uniform grids. The properties of the method are based on the introduction of an…
We show that the 3+1 vacuum Einstein field equations in Ashtekar's variables constitutes a first order symmetric hyperbolic system for arbitrary but fixed lapse and shift fields, by suitable adding to the system terms proportional to the…
A numerical solution scheme for the Einstein field equations based on generalized harmonic coordinates is described, focusing on details not provided before in the literature and that are of particular relevance to the binary black hole…
A hyperbolic integro-differential equation is considered, as a model problem, where the convolution kernel is assumed to be either smooth or no worse than weakly singular. Well-posedness of the problem is studied in the context of semigroup…
Discontinuous Galerkin (DG) discretizations with exact representation of the geometry and local polynomial degree adaptivity are revisited. Hybridization techniques are employed to reduce the computational cost of DG approximations and…
In this paper, numerical solutions of singularly perturbed boundary value problems are given by using variants of finite element method. Both Galerkin and subdomain Galerkin method based on quadratic B-spline functions are applied over the…
We extend the notion of numerical stability of finite difference approximations to include hyperbolic systems that are first order in time and second order in space, such as those that appear in Numerical Relativity. By analyzing the symbol…
We study gravitational collapse of a spherically symmetric scalar field in Einstein-aether theory (general relativity coupled to a dynamical unit timelike vector field). The initial value formulation is developed, and numerical simulations…
Due to added numerical stabilization (diffusion), the stationary states of numerical methods for hyperbolic problems need not be consistent discretizations of those of the PDEs. A closely related phenomenon is the lack of consistency of…
Standard discontinuous Galerkin methods, based on piecewise polynomials of degree $ \qq=0,1$, are considered for temporal semi-discretization for second order hyperbolic equations. The main goal of this paper is to present a simple and…