Weighted, Anisotropic & Hybrid Index Sets

• Weighted, anisotropic, and hybrid index sets are frameworks that characterize function space regularity near singularities to optimize PDE discretizations.
• They establish rigorous weighted Sobolev and analytic spaces by incorporating local geometric singularities, directional derivatives, and anisotropic scaling.
• These index sets guide optimal mesh design and error estimation in finite element and hp-methods, ensuring enhanced convergence rates on complex domains.

Weighted, anisotropic, and hybrid-regularity index sets constitute a foundational framework for the analysis of singularities and optimal discretizations in partial differential equations (PDEs) on non-smooth domains. These structures rigorously encode the interplay between localized geometric singularities and the directional regularity inherited by solutions. They underpin the design of optimal finite element and $hp$-methods, precise regularity results in function spaces, and the analysis of elliptic and parabolic equations with mixed or anisotropic scaling.

1. Weighted Sobolev and Analytic Spaces: Foundations

Let $\Omega\subset\mathbb{R}^3$ be a bounded polyhedral domain with sets of vertices ($V$), edges ($E$), and faces ($F$). The canonical distances to singular strata are

$$r_V(x) = \mathrm{dist}(x, \{\textrm{vertices}\}),\quad r_E(x) = \mathrm{dist}(x, \{\textrm{open edges}\}),\quad r_F(x) = \mathrm{dist}(x, \{\textrm{faces}\}).$$

Let $\beta = (\beta_V, \beta_E, \beta_F)$ be multi-exponents, typically constrained by the local geometry; e.g., for each vertex $V$ with interior angle $\omega_V$, $0\leq \beta_V<\pi/\omega_V$; analogously for each edge $E$, $0 \leq \beta_E < \pi/\alpha_E$ for dihedral angle $\alpha_E$.

The isotropically weighted Sobolev space is

$$H^{m,\beta}(\Omega) = \left\{ u\in L^2(\Omega): r_V(x)^{\beta_V-|\alpha|} r_E(x)^{\beta_E-|\alpha|} r_F(x)^{\beta_F-|\alpha|} \partial^\alpha u \in L^2(\Omega)\ \forall\,|\alpha|\leq m \right\},$$

equipped with the norm

$$\|u\|_{H^{m,\beta}(\Omega)} = \left( \sum_{|\alpha|\leq m} \int_\Omega r_V(x)^{2(\beta_V-|\alpha|)} r_E(x)^{2(\beta_E-|\alpha|)} r_F(x)^{2(\beta_F-|\alpha|)} |\partial^\alpha u|^2\,dx \right)^{1/2}.$$

This weighted structure naturally defines the domain of the Laplace operator on non-smooth domains and precisely describes regularity loss at geometric singularities (Bacuta et al., 2012, Costabel et al., 2010).

In analytic settings, analogous factorially-weighted seminorms and norms are constructed, e.g.,

$$|u|_{m,\beta;\omega} = \left(\sum_{|\alpha|=m} \|\prod_{c} r_c^{\beta_c+|\alpha|} \partial^{\alpha} u \|_{L^2(\omega)}^2\right)^{1/2},\qquad \|u\|_{A,\omega}^\beta = \sum_{m=0}^\infty \frac{1}{m!} |u|_{m,\beta;\omega},$$

encoding analytic regularity relative to edges and corners (Costabel et al., 2010).

2. Anisotropic Regularity and Multi-index Sets

Weighted regularity can be further refined by exploiting directional regularity. Distinguish mutually orthogonal differential operators

$\partial_e,\ \partial_f,\ \partial_i$

spanning, respectively, the edge direction, the face-transverse (but in-plane), and the fully transverse directions. A multi-index $$\alpha = (\alpha_e, \alpha_f, \alpha_i)\in\mathbb{N}^3$$ allows arbitrary derivatives to be indexed by the number taken in each direction.

Anisotropic regularity results state: for $f\in H^{m-1,\beta}(\Omega)$ and homogeneous boundary data, the solution $u$ to $-\Delta u = f$ on $\Omega$ admits

$$\|u\|_{H^{A,\beta}(\Omega)} = \left( \sum_{\alpha\in A} \int_\Omega r_V^{2(\beta_V-|\alpha|)} r_E^{2(\beta_E-|\alpha|)} r_F^{2(\beta_F-|\alpha|)} | \partial_e^{\alpha_e} \partial_f^{\alpha_f} \partial_i^{\alpha_i} u|^2 \right)^{1/2} \leq C \|f\|_{H^{m-1,\beta}(\Omega)},$$

where the index set $A$ encodes directional regularity constraints, e.g.,

$$A = \left\{ \alpha=(\alpha_e,\alpha_f,\alpha_i):\, 0\leq\alpha_e\leq a_e,\ 0\leq\alpha_f\leq a_f,\ 0\leq\alpha_i\leq a_i,\,\, \alpha_e+\alpha_f\leq b_{ef},\, \ldots\right\}$$

with the parameters determined by the domain geometry (Bacuta et al., 2012). The effect is that higher-order derivatives along edges may be controlled even if global Sobolev regularity is scarce.

In analytic regularity, this is mirrored by introducing weights (e.g., $(r_e/r_V)^{\beta_e+|\alpha_{\perp,e}|}$) and splitting multi-indices into edge-parallel and edge-transverse parts (Costabel et al., 2010).

In frequency space, as for vortex sheet problems, weighted anisotropic Sobolev spaces such as $H^{s,\sigma}_\gamma(\mathbb{R}^2)$ are defined with weights vanishing in singular directions, thus providing one-derivative gain away from critical rays and only tangential regularity otherwise. The index set in this setting is determined by the vanishing properties of the weight function in frequency variables, corresponding to anisotropic regularity dictated by the PDE symbol (Secchi, 2020).

3. Hybrid Index Sets and Flagged Spaces

Hybrid index sets unite isotropic and anisotropic viewpoints. The "hybrid" index set

$$H_{m,\mu} = \{ \alpha : |\alpha|\leq m \} \cup E_\mu, \text{ where } E_\mu = \{ \alpha=(\alpha_e,\alpha_f,\alpha_i): 0\leq\alpha_e\leq\mu_E,\,0\leq\alpha_f\leq\mu_F,\,0\leq\alpha_i\leq\mu_I \},$$

enables the enforcement of global $H^m$ control while exploiting additional regularity along principal directions (Bacuta et al., 2012).

In analytic settings, hybrid or "flagged" spaces $J^n_\beta(\Omega; V_0, E_0)$ require derivatives to satisfy different weights depending on the membership of a geometric singularity in subsets $V_0$ (special Dirichlet corners) and $E_0$ (edges with reduced weight). This construction interpolates between fully anisotropic and isotropic regularity regimes, matching the specific boundary condition and variational structure present in the problem (Costabel et al., 2010).

Hybrid index sets also appear in harmonic analysis and maximal regularity theory, as in weighted anisotropic mixed-norm function spaces. For example, in Lizorkin-Triebel spaces $F_{\vec{p},q}^{s,A}(\mathbb{R}^n, w; X)$, intersection representations across variable blocks partition the anisotropic scaling and integrability requirements, generalizing the Fubini property and enabling sharp boundary data descriptions in maximal regularity problems (Lindemulder, 2019).

4. Mesh Design, Discretization, and hp-FEM

The structure of weighted and anisotropic index sets directly informs optimal mesh refinement strategies for finite element and $hp$ methods on polyhedral and polygonal domains. For near-singularity mesh refinement:

• Near an edge $E$, set the mesh size $h_E \simeq h^{1/\mu_E}$, where $\mu_E$ is the highest controlled directional derivative in the edge direction.
• Near a face $F$, analogous relations $h_F\simeq h^{1/\mu_F}$ are used, with elements thin in the transverse direction.

This targeted grading ensures quasi-optimal finite element convergence rates: $$\|u-u_h\|_{H^1(\Omega)} \leq C h^{\min\{\mu_E,\mu_F,\mu_I\}} \|f\|_{L^2(\Omega)},$$ or equivalently in terms of degrees of freedom $N$, $C N^{-m/3} \|f\|_{H^{m-1}(\Omega)}$ for polynomial degree $m$ (Bacuta et al., 2012).

In the $hp$-FEM context, the analytic regularity classes $B_\beta^\ell(\Omega)$ are linked, via the index set correspondence, with exponential convergence rates provided that the polynomial degree $p_e$ is chosen to match the local anisotropic weight $\beta_e$, and mesh grading $h_e\sim p_e^{-1/(1+\beta_e)}$ compensates for singularity strength precisely (Costabel et al., 2010).

5. Index Set Structures in Function Space Theory

Weighted, anisotropic, and hybrid index sets have a parallel in harmonic analysis, especially in the context of weighted mixed-norm Besov and Lizorkin–Triebel spaces: $$B_{\vec{p},q}^{s,A}(\mathbb{R}^n, w; X),\qquad F_{\vec{p},q}^{s,A}(\mathbb{R}^n, w; X),$$ where $A=(A_1, ..., A_\ell)$ encodes the anisotropic dilation, $w$ is a product weight across variable blocks, and $\vec{p}$ specifies mixed $L^p$-integrability. These spaces admit intersection representations, e.g.,

$$F_{\vec{p},q}^{s,A}(\mathbb{R}^n, w; X) = \bigcap_{m=1}^L F_{\vec{p}_{J_m},q}^{s,A_{J_m}}(\mathbb{R}^{n_{J_m}},w_{J_m};X),$$

for any partition $\{J_1,\dots,J_L\}$ of variables, with full equivalence of norms (Lindemulder, 2019). This result refines the classical Fubini property, showing that the intersection holds even for distinct integrability exponents and general anisotropic scaling.

This structure is crucial in the analysis of parabolic and boundary value problems with mixed space-time regularity, where sharp data spaces are characterized as intersections of anisotropic function spaces tailored by the underlying index sets.

6. Synthesis and Applications

The combined apparatus of weighted, anisotropic, and hybrid index sets achieves:

• Rigorous description and control of singular behaviors near geometric singularities in domains.
• Precise regularity estimates (both Sobolev and analytic) necessary for optimal discretization and error estimation in finite element and $hp$ schemes.
• Flexible and sharp modeling of regularity in anisotropic settings, aligning higher regularity in favorable directions with weaker control in singular directions.
• Direct application to advanced harmonic analysis, intersection theorems, and boundary data spaces for PDEs on non-smooth domains.
• Optimal mesh design and degree distribution in $hp$-FEM, ensuring exponential or quasi-optimal convergence rates matched to solution regularity via the index set formalism.

The articulation of these index sets provides a unified framework for understanding, quantifying, and exploiting the complexity of regularity in PDEs and the function spaces governing their solutions (Bacuta et al., 2012, Costabel et al., 2010, Lindemulder, 2019, Secchi, 2020).