[Paper Review] An Indefinite Convection-Diffusion Operator With Real Spectrum

John Weir|ArXiv.org|Nov 8, 2007

Advanced Mathematical Modeling in Engineering1 references3 citations

TL;DR

This paper proves that the eigenvalues of an indefinite convection-diffusion operator $ H $, modeling liquid film dynamics in a rotating cylinder, are all real for $ 0 < \varepsilon < 2 $, using analytic continuation of eigenfunctions to the unit disk and transforming the problem into a self-adjoint differential equation on $[0,1]$ via Heun's equation. The key result confirms a conjecture by Benilov et al. and extends Davies' earlier spectral results by establishing full reality of the spectrum.

ABSTRACT

We confirm rigorously the conjecture, based on numerical and asymptotic evidence, that all the eigenvalues of a certain non-self-adjoint operator are real.

Motivation & Objective

• To resolve the conjecture that the spectrum of the indefinite convection-diffusion operator $ H $ is real, as posed in Benilov, O’Brien, and Sazonov (2007).
• To establish the reality of all eigenvalues of $ H $, building on Chugunova and Pelinovsky’s numerical evidence and Davies’ prior results on eigenvalue existence.
• To demonstrate that the operator $ -iH $, equivalent to a tridiagonal matrix $ A $, has only real eigenvalues by analyzing its action on $ l^2(\mathbf{Z}_+) $.
• To show that eigenfunctions of $ A_+ $ can be analytically continued to the unit disk, transforming the eigenvalue problem into a self-adjoint differential equation on $[0,1]$.

Proposed method

• The eigenvalue problem $ -iHf = \lambda f $ is transformed into a discrete tridiagonal matrix equation $ Av = \lambda v $ via Fourier series expansion, with $ A $ acting on $ l^2(\mathbf{Z}) $.
• The analysis is reduced to $ A_+ $, the restriction of $ A $ to $ l^2(\mathbf{Z}_+) $, due to unitary equivalence between $ A_+ $ and $ A_- $.
• The generating function $ u(z) = \sum_{k=1}^\infty v_k z^k $ is constructed and shown to satisfy a second-order Fuchsian differential equation (Heun’s equation) with regular singularities at $ 0, 1, -1, \infty $.
• Analytic continuation of $ u(z) $ to the unit disk is used to show that $ u $ is analytic on an open set containing $[0,1]$, enabling spectral analysis on $[0,1]$.
• The differential equation is rewritten in Sturm-Liouville form $ -(pu')' + qu = \mu w u $ on $[0,1]$ with $ \mu = 2\lambda/\varepsilon $, where $ p > 0 $, $ w > 0 $, and $ w $ has a pole at 0.
• Self-adjointness of the operator on $[0,1]$ with respect to the weight $ w $ is used to show $ \mu \in \mathbb{R} $, hence $ \lambda \in \mathbb{R} $.

Experimental results

Research questions

• RQ1 Are all eigenvalues of the indefinite convection-diffusion operator $ H $ real for $ 0 < \varepsilon < 2 $?
• RQ2 Can the eigenfunctions of $ H $ be analytically continued to the unit disk to enable spectral analysis via self-adjoint differential equations?
• RQ3 Does the transformation of the discrete eigenvalue problem into a Heun-type differential equation on $[0,1]$ lead to a self-adjoint operator with real spectrum?
• RQ4 Is the spectrum of $ -iH $ equal to its set of eigenvalues, and are all eigenvalues real?
• RQ5 Can Davies’ bound on Fourier coefficients be used to ensure absolute convergence of the generating function and analytic continuation?

Key findings

• All eigenvalues of the operator $ A_+ $, and hence of $ A $, are real for $ 0 < \varepsilon < 2 $ and $ 1/\varepsilon \notin \mathbb{Z} $.
• The eigenfunction generating function $ u(z) = \sum_{k=1}^\infty v_k z^k $ is absolutely convergent for $ |z| \leq 1 $, due to Davies’ bound on Fourier coefficients.
• The eigenfunction $ u(z) $ extends analytically to an open set containing $[0,1]$, allowing the use of differential equation techniques on $[0,1]$.
• The differential equation (8) is shown to be self-adjoint with respect to the weight $ w(x) = x^{-1}(1-x)^{1/\varepsilon}(x+1)^{-1/\varepsilon} $ on $[0,1]$, ensuring real $ \mu $.
• Since $ \mu = 2\lambda/\varepsilon $ and $ \mu \in \mathbb{R} $, it follows that $ \lambda \in \mathbb{R} $, proving the reality of all eigenvalues.
• The spectrum of $ -iH $ is purely discrete and equal to its set of eigenvalues, all of which are real, confirming the conjecture in [1].