Eigenvalues of the Laplacian on perforated domains

This is based on a discussion I had with some colleagues, after attending the 9th Cascade RAIN, in OSU, Corvallis, OR. In particular, we were talking of the ongoing work of Shivam Patel (under the supervision of Dr. Jeffrey Ovall).

The "cheese grater" problem

(I invented this name). The idea is simple: consider the Laplace eigenvalue problem over Ωn=Ω(n;θ)\Omega_n = \Omega(n;\theta):

Δun=λnun, in Ωn,-\Delta u^n = \lambda^n u^n, \text{ in } \Omega^n,

and γnu=0\gamma^n u=0 on its boundary Ωn\partial \Omega^n. Here Δ\Delta is the Laplace operator, γn\gamma^n is a trace operator, θ\theta is a set of parameters that define the geometry of the domain, and nN.n\in\mathbb{N}.

A concrete example is the following:

In other words, we are removing small balls from the unit square[1].

From Shivam's presentation, at least in the Neumann case, it was observed that the eigenvalues somewhat stabilize as nn increases. Thus, the question is whether convergence can be established in this case, for different settings, and what is the dependence on the geometry of the domain.

Sierpinski carpet in the kitchen?

One concerns is quite natural, and pertains to the geometry of the domain. In what sense does the domain converge to a limit? More over, the limit geometry may not be a domain at all: not open, not simply connected, etc. It may even resemble a Sierpinski carpet, or a Cantor set[2].

So, my concern is the stability of the formulations and the spaces, as regularity theory stops working as nicely as we would like. Some Sobolev-like spaces can be defined; see Hu, J. (2003)., for example.

What happens in the kitchen stays in the kitchen

I had a quick though I have shared with some friends: suppose we can nest the domains ιn:Ωn+1Ωn\iota_n:\Omega_{n+1}\hookrightarrow\Omega_{n}. Then, naturally, we have a Sobolev-type contravariant functor V()V(\cdot), such that

V(Ωn)V(Ωn+1),V(\Omega_{n}) \hookrightarrow V(\Omega_{n+1}),

given by ffιnf\mapsto f\circ \iota_n.

The idea –despite the categorical language– is to use Courant-Fischer min-max principle to show monotonicity of the eigenvalues. The main issue being, of course, that a restriction operator may not preserve subspaces...

λkn=minVV(Ωn),dimV=kmaxuVΩnu2Ωnu2?λkn+1=minVV(Ωn+1),dimV=kmaxuVΩn+1u2Ωn+1u2.\lambda_k^n = \min_{V \subset V(\Omega_n), \dim V = k} \max_{u\in V} \frac{\int_{\Omega_n} |\nabla u|^2}{\int_{\Omega_n} |u|^2} \stackrel{?}{\geq} \lambda_k^{n+1} = \min_{V \subset V(\Omega_{n+1}), \dim V = k} \max_{u\in V} \frac{\int_{\Omega_{n+1}} |\nabla u|^2}{\int_{\Omega_{n+1}} |u|^2}.

[1] Let me know if this makes sense. Otherwise, the radius can be changed to be rn=12n+1r_n=\frac{1}{2n+1}, for example.
[2] The carpet and the Cantor set.

CC BY-SA 4.0 Gabriel Pinochet-Soto. Last modified: June 01, 2025. Website built with Franklin.jl and the Julia programming language.