Related papers: Revisiting high-order Taylor methods for astrodyna…
By extending Takens' embedding theorem (1981), Deyle and Sugihara (2011) provided a theoretical justification for using parallel measurement time series to reconstruct a system's attractor. Building on Takens' framework, Brunton et al.…
We introduce an extension to the AthenaK code for general-relativistic magnetohydrodynamics (GRMHD) in dynamical spacetimes using a 3+1 conservative Eulerian formulation. Like the fixed-spacetime GRMHD solver, we use standard finite-volume…
Analytic-bilinear approach for construction and study of integrable hierarchies, in particular, the KP hierarchy is discussed. It is based on the generalized Hirota identity. This approach allows to represent generalized hierarchies of…
The elucidation of many physical problems in science and engineering is subject to the accurate numerical modelling of complex wave propagation phenomena. Over the last decades, high-order numerical approximation for partial differential…
A globally converging numerical method to solve coupled sets of non-linear integral equations is presented. Such systems occur e.g. in the study of Dyson-Schwinger equations of Yang-Mills theory and QCD. The method is based on the knowledge…
This thesis deals with the investigation of a H(div)-conforming hybrid discontinuous Galerkin discretization for incompressible turbulent flows. The discretization method provides many physical and solving-oriented properties, which may be…
We present a new 1-D multi-physics simulation code with use cases intended for, but not limited to, hydrodynamic escapeproblems of planetary atmospheres and planetary accretion models. Our formulation treats an arbitrary number of species…
In this paper, we propose a first order energy stable linear semi-implicit method for solving the Allen-Cahn-Ohta-Kawasaki equation. By introducing a new nonlinear term in the Ohta-Kawasaki free energy functional, all the system forces in…
Second-order Newton-type algorithms that leverage the exact Hessian or its approximation are central to solve nonlinear optimization problems. However, their applications in solving large-scale nonconvex problems are hindered by three…
This paper continues the development of a discrete space-time concept that is recently used in the explanation of the cosmological constant. Instead of order estimation, a more theoretical treatment of the theory is introduced. Based on the…
An additive Runge-Kutta method is used for the time stepping, which integrates the linear stiff terms by an explicit singly diagonally implicit Runge-Kutta (ESDIRK) method and the nonlinear terms by an explicit Runge-Kutta (ERK) method. In…
In this article, we study the damped time-harmonic Galbrun's equation which models solar and stellar oscillations. We introduce and analyze hybrid discontinuous Galerkin discretizations (HDG) that are stable and optimally convergent for all…
The high-order numerical solution of the non-linear shallow water equations (and of hyperbolic systems in general) is susceptible to unphysical Gibbs oscillations that form in the proximity of strong gradients. The solution to this problem…
We present a new symplectic integrator designed for collisional gravitational $N$-body problems which makes use of Kepler solvers. The integrator is also reversible and conserves 9 integrals of motion of the $N$-body problem to machine…
Although the discovery of the chaotic motion of the inner planets in the solar system dates back to more than thirty years ago, the secular chaos of their orbits still dares more analytical analyses. Apart from the high-dimensional…
We study the integrability of a two-dimensional Hamiltonian system with a gyroscopic term and a non-homogeneous potential composed of two homogeneous components of different degrees. The model describes the motion of a particle in a plane…
Hirota's bilinear method ("direct method") has been very effective in constructing soliton solutions to many integrable equations. The construction of one- and two-soliton solutions is possible even for non-integrable bilinear equations,…
A new iterative solver is proposed to efficiently calculate the ground state electronic structure in Density Functional Theory calculations. This algorithm is particularly useful for simulating physical systems considered difficult to…
This research explores the development and application of the High-Order Dynamic Integration Method for solving integro-differential equations, with a specific focus on turbulent fluid dynamics. Traditional numerical methods, such as the…
Modeling interacting galaxies to reproduce observed systems is still a challenge due to the extended parameter space (among other problems). Orbit and basic galaxy parameters can be tackled by fast simulation techniques like the restricted…