Related papers: Construction of second-order six-dimensional Hamil…
Black hole binaries on non-eccentric orbits form an important subclass of gravitational wave sources, but it is a non-trivial issue to construct numerical initial data with minimal initial eccentricity for numerical simulations. We compute…
In this paper, higher-order perturbation theory is applied and tailored to one-dimensional ring-shaped Bose-Hubbard systems. Spectral and geometrical properties are used to structurally simplify the contributions and reduce computational…
Motivated by the desire to understand the leading order nonlinear gravitational wave interactions around arbitrarily rapidly rotating Kerr black holes, we describe a numerical code designed to compute second order vacuum perturbations on…
In this paper, we consider numerical approximations for the viscous Cahn-Hilliard equation with hyperbolic relaxation. This type of equations processes energy-dissipative structure. The main challenge in solving such a diffusive system…
Quantum computation is one of the most promising new paradigms for the simulation of physical systems composed of electrons and atomic nuclei, with applications in chemistry, solid-state physics, materials science, and molecular biology.…
We present a new third-order, semi-discrete, central method for approximating solutions to multi-dimensional systems of hyperbolic conservation laws, convection-diffusion equations, and related problems. Our method is a high-order extension…
A novel routine to investigate the scalar fields in a cosmological context is discussed in the framework of the Hamiltonian formalism. Starting from the Einstein-Hilbert action coupled to a Lagrangian density that contains two components -…
We present an efficient second-order finite difference scheme for solving the 2D sine-Gordon equation, which can inherit the discrete energy conservation for the undamped model theoretically. Due to the semi-implicit treatment for the…
We present a novel positive kinetic scheme built on the efficient collide-and-stream algorithm of the lattice Boltzmann method (LBM) to address hyperbolic conservation laws. We focus on the compressible Euler equations with strong…
In this paper, we present a class of nonuniform time-stepping, high-order linear stabilized schemes that can preserve both the discrete energy stability and maximum-bound principle (MBP) for the time-fractional Allen-Cahn equation. To this…
In this paper we focus on the Cahn-Hilliard equation with dynamic boundary conditions, by adding two hyperbolic relaxation terms to the system. We verify that the energy of the total system is decreasing with time. By adding two…
Accretion onto black holes and compact stars brings material in a zone of strong gravitational and electromagnetic fields. We study dynamical properties of motion of electrically charged particles forming a highly diluted medium (a corona)…
We propose a class of weighted compact central (WCC) schemes for solving hyperbolic conservation laws. The linear version can be considered as a high-order extension of the central Lax-Friedrichs (LxF) scheme and the central conservation…
We present two strategies for designing passivity preserving higher order discretization methods for Maxwell's equations in nonlinear Kerr-type media. Both approaches are based on variational approximation schemes in space and time. This…
We present a novel methodology for constructing arbitrarily high-order structure-preserving methods tailored for damped Hamiltonian systems. This method combines the idea of exponential integrator and energy-preserving collocation methods,…
Inspired by the so-called Palatini formulation of General Relativity and of its modifications and extensions, we consider an analogous formulation of the dynamics of a self-interacting gauge field which is determined by non-linear extension…
Gyrocenter dynamics of charged particles plays a fundamental role in plasma physics. In particular, accuracy and conservation of energy are important features for correctly performing long-time simulations. For this purpose, we here propose…
This paper is the second in a series of two, and describes the current state of the art in modelling and prediction of chaotic time series. Sampled data from deterministic non-linear systems may look stochastic when analysed with linear…
A third order shock-capturing numerical scheme for three-dimensional special relativistic magnetohydrodynamics (3-D RMHD) is presented and validated against several numerical tests. The simple and efficient central scheme described in Paper…
We introduce a second-order numerical scheme for compressible atmospheric motions at small to planetary scales. The collocated finite volume method treats the advection of mass, momentum, and mass-weighted potential temperature in…