Related papers: Arbitrary precision computation of hydrodynamic st…
To provide mathematically rigorous eigenvalue bounds for the Steklov eigenvalue problem, an enhanced version of the eigenvalue estimation algorithm developed by the third author is proposed, which removes the requirements of the positive…
We revisit the numerical stability of the two-level orthogonal Arnoldi (TOAR) method for computing an orthonormal basis of a second--order Krylov subspace associated with two given matrices. We show that the computed basis is close (on…
Spectral methods are well suited for solving hydrodynamic problems in which the self-gravity of the flow needs to be considered. Because Poisson's equation is linear, the numerical solution for the gravitational potential for each…
A certified strategy for determining sharp intervals of enclosure for the eigenvalues of matrix differential operators with singular coefficients is examined. The strategy relies on computing the second order spectrum relative to subspaces…
We present a Chebyshev collocation method for linear ODE and DDE problems. We first give a posteriori estimates for the accuracy of the approximate solution of a scalar ODE initial value problem. Examples of the success of the estimate are…
Exact solutions to the equations of hydrodynamics provide valuable benchmark tests for numerical hydrodynamic codes and also provide useful insights into the nature of hydrodynamic flow. In this paper, we introduce two novel, closely…
We establish the first globally convergent algorithms for computing the Kreiss constant of a matrix to arbitrary accuracy. We propose three different iterations for continuous-time Kreiss constants and analogues for discrete-time Kreiss…
Local polynomial smoothing is a widespread technique in data analysis, and Savitzky-Golay (SG) filters are one of its most well-known realizations. In real settings, the effectiveness of SG filtering depends critically on proper tuning of…
This paper presents a new approach for the computation of eigenvalues of the generalized spheroidal wave equations. The novelty of the present method is in the use of the analytical derivatives of the eigenvalues to minimize losses in…
We consider a class of Fuchsian equations that, for instance, describes the evolution of compressible fluid flows on a cosmological spacetime. Using the method of lines, we introduce a numerical algorithm for the singular initial value…
The application of the Kudryashov method for finding exact solutions of the high order nonlinear evolution equations is considered. Some classes of solitary wave solutions for the families of nonlinear evolution equations of fifth, sixth…
This paper presents a more stable implementation and a highly accurate numerical tool for predicting flooding in urban areas. We started with the (linearised) well-posedness analysis by [1], where far-field boundary conditions were proposed…
The transient start-up flow solution with slip is a useful tool to verify Computational Fluid Dynamics (CFD) simulations. However, a highly accurate, open-source black box solution does not seem to be available. Our method provides a fast,…
In this chapter, we aim at presenting the basic techniques necessary to go beyond the widely accepted paradigm of second-order numerics. We specifically focus on finite-volume schemes for hyperbolic conservation laws occuring in fluid…
A new deterministic floating-point arithmetic called precision arithmetic is developed to track precision for arithmetic calculations. It uses a novel rounding scheme to avoid excessive rounding error propagation of conventional…
In the subcritical interval of the Reynolds number 4320\leq R\leq R_c\equiv 5772, the Navier--Stokes equations of the two--dimensional plane Poiseuille flow are approximated by a 22--dimensional Galerkin representation formed from…
An algorithm for computing eigenvalues and eigenfunctions of the angular spheroidal wave equation, based on a known but scarcely used method, is developed. By requiring the regularity of the wave function, represented by its series…
We present a systematic approach to deriving normal forms and related amplitude equations for flows and discrete dynamics on the center manifold of a dynamical system at local bifurcations and unfoldings of these. We derive a general,…
We develop a two-dimensional high-order numerical scheme that exactly preserves and captures the moving steady states of the shallow water equations with topography or Manning friction. The high-order accuracy relies on a suitable…
This paper presents entropy symmetrization and high-order accurate entropy stable schemes for the relativistic magnetohydrodynamic (RMHD) equations. It is shown that the conservative RMHD equations are not symmetrizable and do not admit a…