Computing the spectrum of non self-adjoint Sturm-Liouville problems with parameter dependent boundary conditions

This paper deals with the computation of the eigenvalues of non self-adjoint Sturm-Liouville problems with parameter dependent boundary conditions using the \textit{regularized sampling method}. A few numerical examples among which singular ones will be presented to illustrate the merit of the method and comparison made with the exact eigenvalues when they are available

Density-potential mappings in quantum dynamics

In a recent letter [Europhys. Lett. 95, 13001 (2011)] the question of whether the density of a time-dependent quantum system determines its external potential was reformulated as a fixed point problem. This idea was used to generalize the existence and uniqueness theorems underlying time-dependent density functional theory. In this work we extend this proof to allow for more general norms and provide a numerical implementation of the fixed-point iteration scheme. We focus on the one-dimensional case as it allows for a more in-depth analysis using singular Sturm-Liouville theory and at the same time provides an easy visualization of the numerical applications in space and time. We give an explicit relation between the boundary conditions on the density and the convergence properties of the fixed-point procedure via the spectral properties of the associated Sturm-Liouville operator. We show precisely under which conditions discrete and continuous spectra arise and give explicit examples. These conditions are then used to show that in the most physically relevant cases the fixed point procedure converges. This is further demonstrated with an example.Comment: 20 pages, 8 figures, 3 table

Incompressible viscous fluid flows in a thin spherical shell

Linearized stability of incompressible viscous fluid flows in a thin spherical shell is studied by using the two-dimensional Navier--Stokes equations on a sphere. The stationary flow on the sphere has two singularities (a sink and a source) at the North and South poles of the sphere. We prove analytically for the linearized Navier--Stokes equations that the stationary flow is asymptotically stable. When the spherical layer is truncated between two symmetrical rings, we study eigenvalues of the linearized equations numerically by using power series solutions and show that the stationary flow remains asymptotically stable for all Reynolds numbers.Comment: 28 pages, 10 figure

Eigenvalue enclosures and exclosures for non-self-adjoint problems in hydrodynamics

In this paper we present computer-assisted proofs of a number of results in theoretical fluid dynamics and in quantum mechanics. An algorithm based on interval arithmetic yields provably correct eigenvalue enclosures and exclosures for non-self-adjoint boundary eigenvalue problems, the eigenvalues of which are highly sensitive to perturbations. We apply the algorithm to: the Orr-Sommerfeld equation with Poiseuille profile to prove the existence of an eigenvalue in the classically unstable region for Reynolds number R=5772.221818; the Orr-Sommerfeld equation with Couette profile to prove upper bounds for the imaginary parts of all eigenvalues for fixed R and wave number α; the problem of natural oscillations of an incompressible inviscid fluid in the neighbourhood of an elliptical flow to obtain information about the unstable part of the spectrum off the imaginary axis; Squire's problem from hydrodynamics; and resonances of one-dimensional Schrödinger operators