2025-11-11T11:22:09.452968

The Stokes Eigenvalue Problem on balls and annuli in three dimensions: Solutions with Poloidal and Toroidal Fields

Rummler, Thäter
We consider the Stokes eigenvalue problem in open balls and open annuli in R3 with homogeneous Dirichlet boundary conditions. Using the frame of toroidal and poloidal fields we construct the othogonal decomposition of the Stokes eigenvalue problem in problems for toroidal and poloidal eigenfunctions. This provides the proof of the completeness of a system of explicitly calculated Stokes eigenfunctions given by one of the authors in 1999, [14].
academic

The Stokes Eigenvalue Problem on balls and annuli in three dimensions: Solutions with Poloidal and Toroidal Fields

Basic Information

  • Paper ID: 2408.06948
  • Title: The Stokes Eigenvalue Problem on balls and annuli in three dimensions: Solutions with Poloidal and Toroidal Fields
  • Authors: Bernd Rummler (Otto-von-Guericke-Universität Magdeburg), Gudrun Thäter (KIT, Karlsruhe Institute of Technology)
  • Classification: math.AP (Analysis of PDEs)
  • Submission Date: August 2024
  • Paper Link: https://arxiv.org/abs/2408.06948v2

Abstract

This paper investigates the Stokes eigenvalue problem with homogeneous Dirichlet boundary conditions in three-dimensional open balls and open annular regions. Utilizing the framework of toroidal and poloidal fields, the authors construct an orthogonal decomposition of the Stokes eigenvalue problem into toroidal and poloidal eigenfunction problems. This provides a proof for the completeness of the explicit system of Stokes eigenfunctions presented by one of the authors in 1999.

Research Background and Motivation

  1. Core Problem: Solving the Stokes eigenvalue problem in special regions with high geometric symmetry (open balls and open annular regions), which is of significant importance in fluid mechanics applications.
  2. Problem Significance:
    • Flow problems in open balls and open annular regions play important roles in applications
    • The high geometric symmetry enables obtaining more precise information compared to general domains
    • Facilitates analytical and numerical solutions of flow problems
  3. Limitations of Existing Methods:
    • Lack of rigorous proofs for the completeness of Stokes operator eigenfunction systems
    • Need for more systematic approaches to handle decomposition of divergence-free vector fields
  4. Research Motivation:
    • Establish a general method for decomposing divergence-free vector fields
    • Prove the completeness of previously computed real-valued eigenfunction systems
    • Provide theoretical foundations for Galerkin methods of the Navier-Stokes equations

Core Contributions

  1. Established Systematic Decomposition Theory: Proposed a complete theoretical framework for decomposing divergence-free vector fields into toroidal and poloidal fields
  2. Proved Completeness of Eigenfunction Systems: Provided rigorous completeness proofs for the Stokes eigenfunction system explicitly computed in 1999
  3. Constructed Orthogonal Decomposition: Orthogonally decomposed the Stokes eigenvalue problem into independent toroidal and poloidal eigenfunction problems
  4. Provided Explicit Solutions: Derived explicit expressions for Stokes eigenfunctions in open balls and open annular regions along with corresponding transcendental equations

Detailed Methodology

Problem Formulation

The paper studies the following Stokes eigenvalue problem:

Problem 1: Find solutions uu, λ\lambda, and pp (for σ:0σ<1\sigma: 0 \leq \sigma < 1) satisfying: Δu+p=λuand÷u=Tu=0inΩσ-\Delta u + \nabla p = \lambda u \quad \text{and} \quad \div u = \nabla^T \cdot u = 0 \quad \text{in} \quad \Omega_\sigma Boundary condition: u=0u = 0 on Ωσ\partial\Omega_\sigma

Where:

  • Ωo:={xR3:xE<1}\Omega_o := \{x \in \mathbb{R}^3 : \|x\|_E < 1\} is the open unit ball
  • Ωσ:={xR3:σ<xE<1}\Omega_\sigma := \{x \in \mathbb{R}^3 : \sigma < \|x\|_E < 1\} is the open annular ball

Core Theoretical Framework

1. Definition of Toroidal and Poloidal Fields

Definition 4: Toroidal and poloidal fields are defined as: t:=curl(ψx)=gradψ×xt := \text{curl}(\psi x) = \text{grad}\psi \times xp:=curl(curl(χx))p := \text{curl}(\text{curl}(\chi x))

where ψ\psi and χ\chi are called the toroidal and poloidal potentials, respectively.

2. Decomposition Theorem

Theorem 1 (Decomposition Theorem): Let uW22(Ωσ)u \in W_2^2(\Omega_\sigma) be a divergence-free vector field defined on Ωσ\Omega_\sigma satisfying the spherical mean value condition. Then there exist uniquely determined toroidal vector field tt and poloidal vector field pp such that u=t+pu = t + p (at least in the L2(Ωσ)L^2(\Omega_\sigma) sense).

3. Orthogonality Properties

The orthogonality of toroidal and poloidal fields in the L2(Ωσ)L^2(\Omega_\sigma) and W21(Ωσ)W_2^1(\Omega_\sigma) senses is established:

For any toroidal field tt and poloidal field pp: ΩσtsTpsdx=0\int_{\Omega_\sigma} t_s^T p_s dx = 0

Technical Innovations

  1. Systematic Application of Mie Representation: Systematically applies the divergence-free field decomposition method introduced by Gustav Mie in 1908 to the Stokes eigenvalue problem
  2. Clever Use of Spherical Harmonics: Utilizes eigenfunctions of the Laplace-Beltrami operator (spherical harmonics) as technical tools
  3. Uniqueness Guarantee: Ensures uniqueness of decomposition through vanishing mean value conditions on the sphere
  4. Derivation of Scalar Potential Equations: Transforms the vector eigenvalue problem into boundary value problems for scalar potentials

Main Results

Governing Equations for Scalar Potentials

Toroidal Eigenfunctions

Problem 3: For 0σ<10 \leq \sigma < 1, find solutions ψ\psi and λ\lambda satisfying: Δψ=λψinΩσ-\Delta\psi = \lambda\psi \quad \text{in} \quad \Omega_\sigma Boundary conditions: ψ=0\psi = 0 on Ωσ\partial\Omega_\sigma and satisfying the mean value condition

Poloidal Eigenfunctions

Problem 4: Find solutions χ\chi and λ\lambda satisfying: Δ2χ=λΔχinΩσ-\Delta^2\chi = \lambda\Delta\chi \quad \text{in} \quad \Omega_\sigma Boundary conditions: χ=χr=0\chi = \frac{\partial\chi}{\partial r} = 0 on Ωσ\partial\Omega_\sigma

Explicit Solutions in the Unit Ball

Theorem 3: For the unit ball Ωo\Omega_o, the complete eigenfunction system of the Stokes operator consists of:

  1. Toroidal Fields: tc,s(j),l,k:=curl(Jl+12(μl+12jr)rZlkx)t_{c,s}^{(j),l,k} := \text{curl}\left(\frac{J_{l+\frac{1}{2}}(\mu^j_{l+\frac{1}{2}}r)}{\sqrt{r}} Z_l^k x\right)

Eigenvalues: λ(j),l,k=(μl+12j)2\lambda^{(j),l,k} = (\mu^j_{l+\frac{1}{2}})^2, multiplicity (2l+1)(2l+1)

  1. Poloidal Fields: pc,s(j),l,k:=curlcurl(χ(j),l,kx)p_{c,s}^{(j),l,k} := \text{curl}\text{curl}(\chi^{(j),l,k} x)

where: χ(j),l,k=(Jl+12(μl+32jr)rμl+32jJl12(μl+32j)rl2l+1)Zlk\chi^{(j),l,k} = \left(\frac{J_{l+\frac{1}{2}}(\mu^j_{l+\frac{3}{2}}r)}{\sqrt{r}} - \frac{\mu^j_{l+\frac{3}{2}} J_{l-\frac{1}{2}}(\mu^j_{l+\frac{3}{2}}) r^l}{2l+1}\right) Z_l^k

Eigenvalues: λ(j),l,k=(μl+32j)2\lambda^{(j),l,k} = (\mu^j_{l+\frac{3}{2}})^2, multiplicity (2l+1)(2l+1)

Eigenvalue Equations in Annular Balls

For the annular ball Ωσ\Omega_\sigma, eigenvalues are determined by the roots of the following transcendental equations:

Toroidal Fields:

J_{l+\frac{1}{2}}(\mu\sigma) & J_{-l-\frac{1}{2}}(\mu\sigma) \\ J_{l+\frac{1}{2}}(\mu) & J_{-l-\frac{1}{2}}(\mu) \end{bmatrix} = 0$$ **Poloidal Fields**: $$\det\begin{bmatrix} J_{l-\frac{1}{2}}(\mu\sigma) & J_{-l+\frac{1}{2}}(\mu\sigma) & \sigma^{l-\frac{1}{2}} & 0 \\ J_{l-\frac{1}{2}}(\mu) & J_{-l+\frac{1}{2}}(\mu) & 1 & 0 \\ J_{l+\frac{3}{2}}(\mu\sigma) & J_{-l-\frac{3}{2}}(\mu\sigma) & 0 & \sigma^{-l-\frac{3}{2}} \\ J_{l+\frac{3}{2}}(\mu) & J_{-l-\frac{3}{2}}(\mu) & 0 & 1 \end{bmatrix} = 0$$ ## Theoretical Significance and Applications ### Properties of the Stokes Operator **Theorem 2**: The Stokes operator $A_\sigma$ is a positive self-adjoint operator whose inverse $A_\sigma^{-1}$ is injective, self-adjoint, and compact. **Lemma 3**: The Stokes operator possesses a pure point spectrum with positive eigenvalues of finite multiplicity, and the corresponding eigenfunctions are complete in the relevant function spaces. ### Numerical Application Value 1. **Galerkin Method**: Provides explicit basis functions for Galerkin methods of incompressible Navier-Stokes equations 2. **Numerical Turbulence Studies**: Can be used to construct Galerkin schemes for turbulence in arbitrary balls and annular regions 3. **Geophysical Applications**: Applications in geomagnetism, particularly in the analysis and interpretation of satellite observation data ## In-Depth Evaluation ### Strengths 1. **Theoretical Completeness**: Provides a complete theoretical framework from decomposition theorems to explicit solution construction 2. **Mathematical Rigor**: Employs modern functional analysis theory with rigorous and complete proofs 3. **Practical Value**: Possesses both theoretical significance and practical tools for numerical computation 4. **Historical Continuity**: Combines classical Mie decomposition theory with modern PDE theory ### Technical Depth 1. **Function Space Theory**: Skillfully applies Sobolev spaces and Hilbert space theory 2. **Operator Theory**: Deeply applies self-adjoint and compact operator theory 3. **Special Functions**: Cleverly utilizes properties of Bessel functions and spherical harmonics ### Limitations 1. **Geometric Restrictions**: Methods apply only to regions with special symmetry (balls and annuli) 2. **Boundary Conditions**: Primarily addresses homogeneous Dirichlet boundary conditions 3. **Dimensional Restrictions**: Specifically targets three-dimensional cases, though some results may generalize ### Impact Assessment 1. **Theoretical Contribution**: Provides new analytical tools for Stokes eigenvalue problems 2. **Computational Significance**: Establishes theoretical foundations for related numerical methods 3. **Interdisciplinary Value**: Possesses application potential in fluid mechanics, geophysics, and related fields ## Conclusions and Future Perspectives This paper successfully establishes a complete theory for the Stokes eigenvalue problem in spherical domains. Through decomposition into toroidal and poloidal fields, it not only proves the completeness of eigenfunction systems but also provides explicit construction methods. This work offers profound insights into the structure of three-dimensional divergence-free vector fields and establishes solid theoretical foundations for related numerical computations. Future research directions may include: 1. Generalization to other symmetric domains 2. Consideration of other types of boundary conditions 3. Verification through applications in practical fluid mechanics problems 4. Integration with modern computational fluid dynamics methods ## References The paper cites 18 important references covering multiple fields including elliptic equation theory, Stokes equations, and functional analysis, reflecting the academic depth and breadth of the work.