Related papers: Convergence and round-off errors in a two-dimensio…
The role of round-off errors on the receptivity and instability of fluid flows are conclusively established for the first time using high accuracy simulations of the benchmark two-dimensional (2D) Taylor-Green vortex problem using double…
When using a finite difference method to solve an initial--boundary--value problem, the truncation error is often of lower order at a few grid points near boundaries than in the interior. Normal mode analysis is a powerful tool to analyze…
This paper focuses on the study of Sturm-Liouville eigenvalue problems. In the classical Chebyshev collocation method, the Sturm-Liouville problem is discretized to a generalized eigenvalue problem where the functions represent interpolants…
Global instability analysis of flows is often performed via time-stepping methods, based on the Arnoldi algorithm. When setting up these methods, several computational parameters must be chosen, which affect intrinsic errors of the…
This study presents a comprehensive spatial eigenanalysis of fully-discrete discontinuous spectral element methods, now generalizing previous spatial eigenanalysis that did not include time integration errors. The influence of discrete time…
The vertical modes of linearized equations of motion are widely used by the oceanographic community in numerous theoretical and observational contexts. However, the standard approach for solving the generalized eigenvalue problem using…
The Arnoldi-Tikhonov method is a well-established regularization technique for solving large-scale ill-posed linear inverse problems. This method leverages the Arnoldi decomposition to reduce computational complexity by projecting the…
Most of the literature on the solution of linear ill-posed operator equations, or their discretization, focuses only on the infinite-dimensional setting or only on the solution of the algebraic linear system of equations obtained by…
We show that by using higher order precision arithmetic, i.e., using floating point types with more significant bits than standard double precision numbers, one may accurately compute eigenvalues for non-normal matrices arising in…
In this work we present an a priori error analysis for solving the unsteady advection equation on cut cell meshes along a straight ramp in two dimensions. The space discretization uses a lowest order upwind-type discontinuous Galerkin…
The sum-of-squares method can give rigorous lower bounds on the energy of quantum Hamiltonians. Unfortunately, typically using this method requires solving a semidefinite program, which can be computationally expensive. Further, the…
Traditionally, finite differences and finite element methods have been by many regarded as the basic tools for obtaining numerical solutions in a variety of quantum mechanical problems emerging in atomic, nuclear and particle physics,…
We consider the forward problem of uncertainty quantification for the generalised Dirichlet eigenvalue problem for a coercive second order partial differential operator with random coefficients, motivated by problems in structural…
Spatial discretisation of geometrically complex computational domains often entails unstructured meshes of general topology for Computational Fluid Dynamics (CFD). Mesh skewness is then typically encountered causing severe deterioration of…
We establish guaranteed and practically computable a posteriori error bounds for source problems and eigenvalue problems involving linear Schr{\"o}dinger operators with atom-centered potentials discretized with linear combinations of atomic…
The efficient and reliable approximation of convection-dominated problems continues to remain a challenging task. To overcome the difficulties associated with the discretization of convection-dominated equations, stabilization techniques…
We calculate accurate eigenvalues and eigenfunctions of the Schr\"odinger equation for a two-dimensional quantum dipole. This model proved useful for the study of elastic effects of a single edge dislocation. We show that the Rayleigh-Ritz…
In the present paper we develop an approach to obtain sharp spectral asymptotics for Steklov type problems on planar domains with corners. Our main focus is on the two-dimensional sloshing problem, which is a mixed Steklov-Neumann boundary…
Two-dimensional (2D) flows are efficiently controlled with spanwise waviness, i.e. spanwise-periodic (SP) wall blowing/suction/deformation. We tackle the global linear stability of 2D flows subject to small-amplitude 3D SP control. Building…
In this work we consider a class of delay eigenvalue problems that admit a spectrum similar to that of a Hamiltonian matrix, in the sense that the spectrum is symmetric with respect to both the real and imaginary axis. More precisely, we…