Related papers: Convergence and round-off errors in a two-dimensio…
We consider filtered subspace iteration for approximating a cluster of eigenvalues (and its associated eigenspace) of a (possibly unbounded) selfadjoint operator in a Hilbert space. The algorithm is motivated by a quadrature approximation…
We present a simple technique for avoiding physically spurious eigenmodes that often occur in the solution of hydrodynamic stability problems by the Chebyshev collocation method. The method is demonstrated on the solution of the…
We develop and analyse residual-based a posteriori error estimates for the virtual element discretisation of a nonlinear stress-assisted diffusion problem in two and three dimensions. The model problem involves a two-way coupling between…
We consider second-order PDE problems set in unbounded domains and discretized by Lagrange finite elements on a finite mesh, thus introducing an artificial boundary in the discretization. Specifically, we consider the reaction diffusion…
This paper is devoted to study the Arnold-Winther mixed finite element method for two dimensional Stokes eigenvalue problems using the stress-velocity formulation. A priori error estimates for the eigenvalue and eigenfunction errors are…
We are concerned with the study of some classical spectral collocation methods as well as with the new software system Chebfun in computing high order eigenpairs of singular and regular Schrodinger eigenproblems. We want to highlight both…
This paper proposes a new hybrid high-order discretization for the biharmonic problem and the corresponding eigenvalue problem. The discrete ansatz space includes degrees of freedom in $n-2$ dimensional submanifolds (e.g., nodal values in…
The second moment method is a linear acceleration technique which couples the transport equation to a diffusion equation with transport-dependent additive closures. The resulting low-order diffusion equation can be discretized independent…
In this article, we prove the convergence of a semi-discrete numerical method applied to a general class of nonlocal nonlinear wave equations where the nonlocality is introduced through the convolution operator in space. The most important…
The randomized Arnoldi process has been used in large-scale scientific computing because it produces a well-conditioned basis for the Krylov subspace more quickly than the standard Arnoldi process. However, the resulting Hessenberg matrix…
This paper studies the $J$-method of [E. Jarlebring, S. Kvaal, W. Michiels. SIAM J. Sci. Comput. 36-4:A1978-A2001, 2014] for nonlinear eigenvector problems in a general Hilbert space framework. This is the basis for variational…
Understanding the orientation of geological structures is crucial for analyzing the complexity of the Earths' subsurface. For instance, information about geological structure orientation can be incorporated into local anisotropic…
We develop a variational multiscale proper orthogonal decomposition reduced-order model for turbulent incompressible Navier-Stokes equations. The error analysis of the full discretization of the model is presented. All error contributions…
We present an efficient method for computing dominant eigenvalues of large, nonsymmetric, diagonalizable matrices based on an adaptive block Lanczos algorithm combined with Chebyshev polynomial filtering. The proposed approach improves…
We obtain asymptotic formulae for the Steklov eigenvalues and eigenfunctions of curvilinear polygons in terms of their side lengths and angles. These formulae are quite precise: the errors tend to zero as the spectral parameter tends to…
A novel numerical approach to solving the shallow-water equations on the sphere using high-order numerical discretizations in both space and time is proposed. A space-time tensor formalism is used to express the equations of motion…
This study presents an efficient, accurate, effective and unconditionally stable time stepping scheme for the Darcy-Brinkman equations in double-diffusive convection. The stabilization within the proposed method uses the idea of stabilizing…
Spectral methods provide an elegant and efficient way of numerically solving differential equations of all kinds. For smooth problems, truncation error for spectral methods vanishes exponentially in the infinity norm and $L_2$-norm.…
A reaction-diffusion problem with a Caputo time derivative is considered. An integral discretization scheme on a graded mesh along with a decomposition of the exact solution is proposed. The truncation error estimate of the discretization…
The systematic errors inherent in the QCD sum rule approach to meson wavefunctions are examined in the context of QCD in 1+1 spacetime dimensions in the large N limit where the theory is exactly solvable. It is shown that the truncation of…