ADER-WENO Finite Volume Schemes with Space-Time Adaptive Mesh Refinement

We present the first high order one-step ADER-WENO finite volume scheme with Adaptive Mesh Refinement (AMR) in multiple space dimensions. High order spatial accuracy is obtained through a WENO reconstruction, while a high order one-step time discretization is achieved using a local space-time discontinuous Galerkin predictor method. Due to the one-step nature of the underlying scheme, the resulting algorithm is particularly well suited for an AMR strategy on space-time adaptive meshes, i.e.with time-accurate local time stepping. The AMR property has been implemented 'cell-by-cell', with a standard tree-type algorithm, while the scheme has been parallelized via the Message Passing Interface (MPI) paradigm. The new scheme has been tested over a wide range of examples for nonlinear systems of hyperbolic conservation laws, including the classical Euler equations of compressible gas dynamics and the equations of magnetohydrodynamics (MHD). High order in space and time have been confirmed via a numerical convergence study and a detailed analysis of the computational speed-up with respect to highly refined uniform meshes is also presented. We also show test problems where the presented high order AMR scheme behaves clearly better than traditional second order AMR methods. The proposed scheme that combines for the first time high order ADER methods with space--time adaptive grids in two and three space dimensions is likely to become a useful tool in several fields of computational physics, applied mathematics and mechanics.Comment: With updated bibliography informatio

Numerical Evidence that the Perturbation Expansion for a Non-Hermitian PT\mathcal{PT}-Symmetric Hamiltonian is Stieltjes

Recently, several studies of non-Hermitian Hamiltonians having $\mathcal{PT}$ symmetry have been conducted. Most striking about these complex Hamiltonians is how closely their properties resemble those of conventional Hermitian Hamiltonians. This paper presents further evidence of the similarity of these Hamiltonians to Hermitian Hamiltonians by examining the summation of the divergent weak-coupling perturbation series for the ground-state energy of the $\mathcal{PT}$-symmetric Hamiltonian $H=p^2+{1/4}x^2+i\lambda x^3$ recently studied by Bender and Dunne. For this purpose the first 193 (nonzero) coefficients of the Rayleigh-Schr\"odinger perturbation series in powers of $\lambda^2$ for the ground-state energy were calculated. Pad\'e-summation and Pad\'e-prediction techniques recently described by Weniger are applied to this perturbation series. The qualitative features of the results obtained in this way are indistinguishable from those obtained in the case of the perturbation series for the quartic anharmonic oscillator, which is known to be a Stieltjes series.Comment: 20 pages, 0 figure

The Olympic medals ranks, lexicographic ordering and numerical infinities

Several ways used to rank countries with respect to medals won during Olympic Games are discussed. In particular, it is shown that the unofficial rank used by the Olympic Committee is the only rank that does not allow one to use a numerical counter for ranking – this rank uses the lexicographic ordering to rank countries: one gold medal is more precious than any number of silver medals and one silver medal is more precious than any number of bronze medals. How can we quantify what do these words, more precious, mean? Can we introduce a counter that for any possible number of medals would allow us to compute a numerical rank of a country using the number of gold, silver, and bronze medals in such a way that the higher resulting number would put the country in the higher position in the rank? Here we show that it is impossible to solve this problem using the positional numeral system with any finite base. Then we demonstrate that this problem can be easily solved by applying numerical computations with recently developed actual infinite numbers. These computations can be done on a new kind of a computer – the recently patented Infinity Computer. Its working software prototype is described briefly and examples of computations are given. It is shown that the new way of counting can be used in all situations where the lexicographic ordering is required

The exact (up to infinitesimals) infinite perimeter of the Koch snowflake and its finite area

The Koch snowflake is one of the first fractals that were mathematically described. It is interesting because it has an infinite perimeter in the limit but its limit area is finite. In this paper, a recently proposed computational methodology allowing one to execute numerical computations with infinities and infinitesimals is applied to study the Koch snowflake at infinity. Numerical computations with actual infinite and infinitesimal numbers can be executed on the Infinity Computer being a new supercomputer patented in USA and EU. It is revealed in the paper that at infinity the snowflake is not unique, i.e., different snowflakes can be distinguished for different infinite numbers of steps executed during the process of their generation. It is then shown that for any given infinite number n of steps it becomes possible to calculate the exact infinite number, Nn, of sides of the snowflake, the exact infinitesimal length, Ln, of each side and the exact infinite perimeter, Pn, of the Koch snowflake as the result of multiplication of the infinite Nn by the infinitesimal Ln. It is established that for different infinite n and k the infinite perimeters Pn and Pk are also different and the difference can be infinite. It is shown that the finite areas An and Ak of the snowflakes can be also calculated exactly (up to infinitesimals) for different infinite n and k and the difference An − Ak results to be infinitesimal. Finally, snowflakes constructed starting from different initial conditions are also studied and their quantitative characteristics at infinity are computed