Related papers: Arbitrary precision computation of hydrodynamic st…
A new implementation of the Roe scheme for solving two-layer shallow-water equations is presented in this paper. The proposed A-Roe scheme is based on the analytical solution to the characteristic quartic of the flux matrix, which is an…
Two asymptotically exact a posteriori error estimates are proposed for eigenvalues by the nonconforming Crouzeix--Raviart and enriched Crouzeix-- Raviart elements. The main challenge in the design of such error estimators comes from the…
In this paper, we compute the H- and Z-eigenvalues of even order symmetric tensors by using the adaptive cubic regularization algorithm.
We present rigorous error estimates towards a first-order unconditionally energy stable scheme designed for 3D hydrodynamic Q-tensor model of nematic liquid crystals. This scheme combines the scalar auxiliary variable (SAV), stabilization…
We propose a second-order accurate method to estimate the eigenvectors of extremely large matrices thereby addressing a problem of relevance to statisticians working in the analysis of very large datasets. More specifically, we show that…
We present a class of arbitrarily high order fully explicit kinetic numerical methods in compressible fluid dynamics, both in time and space, which include the relaxation schemes by S. Jin and Z. Xin. These methods can use CFL number larger…
The manuscript presents a new technique for computing the exponential of skew-Hermitian operators. Principal advantages of the proposed method include: stability even for large time-steps, the possibility to parallelize in time over many…
Zeitlin's model is a discretisation of the 2-D Euler equations that preserves the underlying geometric structure. This feature makes it suitable for studying the qualitative behaviour of the dynamics. Here, we utilise Arnold's geometric…
The QZ algorithm computes the Schur form of a matrix pencil. It is an iterative algorithm and at some point, it must decide that an eigenvalue has converged and move on with another one. Choosing a criterion that makes this decision is…
In the last decade matrix polynomials have been investigated with the primary focus on adequate linearizations and good scaling techniques for computing their eigenvalues and eigenvectors. In this article we propose a new method for…
We solve by Chebyshev spectral collocation some genuinely nonlinear Liouville-Bratu-Gelfand type, 1D and a 2D boundary value problems. The problems are formulated on the square domain $[-1, 1]\times[-1, 1]$ and the boundary condition…
For large-scale eigenvalue problems requiring many mutually orthogonal eigenvectors, traditional numerical methods suffer substantial computational and communication costs with limited parallel scalability, primarily due to explicit…
Hyperbolic conservation laws with stiff source terms appear in the study of a variety of physical systems. Early work showed that the use of formally second-order accurate semi-implicit methods could lead to a substantial loss of accuracy,…
We demonstrate that it is possible to compute wave function normalization constants for a class of Schr\"odinger type equations by an algorithm which scales linearly (in the number of eigenfunction evaluations) with the desired precision P…
The selective frequency damping method was applied to a bent flow. The method was used in an adaptive formulation. The most dangerous frequency was determined by solving an eigenvalue problem. It was found that one of the patterns,…
The existing doubling algorithms have been proven efficient for several important nonlinear matrix equations arising from real-world engineering applications. In a nutshell, the algorithms iteratively compute a basis matrix, in one of the…
A novel method to estimate unsteady aerodynamic force coefficients from pointwise velocity measurements is presented. The methodology is based on a resolvent-based reduced-order model which requires the mean flow to obtain physical flow…
In this article we prove, choosing an appropriately weighted $L_2$-energy equivalent to the classical energy, that the plane Couette and Poiseuille flows are nonlinearly stable with respect to streamwise perturbations for any Reynolds…
We establish sufficiency conditions in order to achieve approximation to discrete eigenvalues of self-adjoint operators in the second-order projection method suggested recently by Levitin and Shargorodsky, [math.SP/0212087]. We find…
In this study linear and nonlinear higher order singularly perturbed problems are examined by a numerical approach, the differential quadrature method. Here, the main idea is using Chebyshev polynomials to acquire the weighting coefficient…