Weighted p-Laplacian Equation: Analysis & Applications

Updated 29 January 2026

Weighted p-Laplacian Equation is a nonlinear differential operator that incorporates spatial weight functions, extending the classical p-Laplacian to study degenerate elliptic and parabolic PDEs.
The operator leverages weighted Sobolev spaces and variational methods to establish existence, uniqueness, and regularity of weak solutions under Aₚ-type conditions.
It plays a central role in spectral theory, geometric analysis, nonlinear evolution problems, and finite element approximations for complex PDE models.

The weighted $p$ -Laplacian equation generalizes the nonlinear $p$ -Laplacian operator by incorporating a spatially varying weight function, appearing in both divergence-form second-order elliptic and parabolic PDEs. In its most generic PDE formulation, it reads

$-\mathrm{div}\left(w(x)|\nabla u|^{p-2}\nabla u\right) = f(x, u, \nabla u)$

on domains in $\mathbb{R}^n$ or on Riemannian manifolds, for a weight $w: \Omega \to (0,\infty)$ , $p>1$ , and suitably regular nonlinearities $f$ . The weighted $p$ -Laplacian forms the foundational quasilinear model in the theory of degenerate and singular elliptic PDEs, has deep interplay with weighted Sobolev and Orlicz spaces, and is central to spectral theory, nonlinear potential theory, calculus of variations, stochastic games, and geometric analysis. Its study involves tools from nonlinear analysis, variational methods, geometric measure theory, and (in the manifold setting) comparison geometry.

1. Definitions, Operator Structure, and Function Spaces

The weighted $p$ -Laplacian operator, often denoted $L_{p,w}$ or $-\mathrm{div}(w|\nabla u|^{p-2}\nabla u)$ in Euclidean domains, generalizes the classical divergence-form elliptic equation by the inclusion of a positive, measurable weight $w(x)$ . In smooth metric-measure spaces $(M,g,e^{-f}dV)$ , the weighted $p$ -Laplacian is given as

$\Delta_{p,f} u = e^{f}\,\mathrm{div}\!\bigl(e^{-f}|\nabla u|^{p-2}\nabla u\bigr)$

where $f$ is a smooth potential (Bakry–Émery theory) (Li et al., 2019, Dai et al., 2020).

Weighted Sobolev spaces $W^{1,p}(\Omega,w)$ consist of functions $u$ for which both $u \in L^p(\Omega)$ and $\nabla u \in L^p(\Omega, w \, dx)$ . The $A_p$ Muckenhoupt weights (those with $w \in L^1_{\mathrm{loc}}(\Omega)$ , $w>0$ a.e., $w^{-\frac{1}{p-1}}\in L^1_{\mathrm{loc}}(\Omega)$ ) guarantee reflexivity, rich embedding, and interpolation properties (1803.02014, Otarola et al., 2024). For variable-exponent problems, double-weighted variable exponent spaces $W^{1,q(.),p(.)}(\Omega, \theta_0, \theta)$ admit distinct weights for the function and its gradient (Unal et al., 2019).

Table: Formal definitions for weighted $p$ -Laplacian operators in selected settings

Setting	Operator Form	Reference
$\mathbb{R}^n$	$-\mathrm{div}(w(x)\|\nabla u\|^{p-2}\nabla u)$	(1803.02014)
Weighted Sobolev space	$W^{1,p}(\Omega,w)$ , norm via $L^p$ and $w\|\nabla u\|^p$	(Otarola et al., 2024)
Metric-measure (Bakry–Émery)	$\Delta_{p,f} u = e^{f}\mathrm{div}(e^{-f}\|\nabla u\|^{p-2}\nabla u)$	(Li et al., 2019)

2. Existence, Uniqueness, and Qualitative Analysis

Weighted $p$ -Laplacian boundary value problems exhibit fine existence, multiplicity, and regularity theory, with variational methods central in both the elliptic and parabolic cases. For Dirichlet or obstacle problems, solutions correspond to critical points of energy functionals over weighted Sobolev spaces: $\mathcal{J}(u) = \int_{\Omega} w(x)\frac{1}{p}|\nabla u|^p - \int_{\Omega} f(x)u,$ with homogeneous Dirichlet or obstacle constraints. Existence follows from coercivity, convexity, and lower semicontinuity provided the weights satisfy $A_p$ -type conditions and the nonlinearity is subcritical or fulfills a weak superlinearity property (Kumar et al., 2022, Unal et al., 2019, Otarola et al., 2024).

Regularity and qualitative properties (positivity, comparison principles, boundedness) are established using Picone identities, Hardy inequalities, maximum principles, and weighted Sobolev embeddings. Uniqueness of positive principal eigenfunctions, simplicity of the first eigenvalue, and Liouville-type non-existence results in unbounded domains have been developed via weighted test function techniques and Hardy/Poincaré inequalities (1803.02014, Drábek et al., 2018).

Problems with variable exponents $p(x)$ , double-phase structures, or singular/degenerate weights require refined compactness and embedding arguments to ensure weak solution existence (Unal et al., 2019, Repovš et al., 2021).

3. Spectral Theory and Variational Characterization

The spectral theory for the weighted $p$ -Laplacian is rich, with principal eigenvalues and spectral gaps having central applications in analysis and geometry. The principal (first nontrivial) eigenvalue in Dirichlet or mixed settings is variationally defined by

$\lambda_{1,p,w} = \inf_{u \neq 0} \frac{\int_\Omega w(x)|\nabla u|^p}{\int_\Omega |u|^p}$

(adjusted for constraints and weights), with corresponding weak eigenfunctions (Li et al., 2019, Drábek et al., 2018, Chen et al., 2015). In Bakry–Émery manifolds, sharp lower bounds for the first nonzero Neumann eigenvalue are given in terms of the weighted Ricci curvature $\mathrm{Ric}_f = \mathrm{Ric} + \nabla^2f$ and the manifold diameter, leading to extensions of the Lichnerowicz–Obata and Zhong–Yang theorems (Li et al., 2019).

Three equivalent variational principles for the principal eigenvalue exist: Rayleigh quotient, single/double integral operator forms, and a differential (quotient) formula, underpinning both qualitative and quantitative eigenvalue estimates (Chen et al., 2015). Spectral simplicity and capacity theory approaches provide further insight, especially when dealing with sign-changing weights or non-standard domains (Drábek et al., 2018, 1803.02014).

4. Nonlinear Evolution, Long-Time Behavior, and Random Media

The (weighted) $p$ -Laplacian evolution equation models nonlinear diffusion with density-dependent diffusivity: $u_t = \mathrm{div}(w(x)|\nabla u|^{p-2}\nabla u)$ with various boundary conditions. When $w \in A_p$ , this initial-value problem admits a unique strong/entropy solution in $L^1$ or $W^{1,p}(w)$ , with explicit $L^q$ decay rates

$\|u(t) - \overline{u}\|_{L^q} \leq C t^{-1/p} \|u_0 - \overline{u}\|_{L^2}^{2/p}$

and, for small $p<2$ , possible finite-time extinction. The weighted mass is always conserved under Neumann boundary conditions (Nerlich, 2016, Nerlich, 2017).

In random media—i.e., when the weight $w$ or diffusion coefficient is itself a stochastic process—the mild and strong solution frameworks extend via nonlinear semigroup theory; mass conservation, $L^1$ contraction, and long-time asymptotics hold pathwise (Nerlich, 2017). The deterministic and random settings share core functional-analytic techniques, with suitable measurability and integrability assumptions compensating for the lack of uniform ellipticity.

5. Geometric, Game-Theoretic, and Inverse Problem Perspectives

Geometry: In smooth metric-measure spaces, the weighted $p$ -Laplacian is fundamental for Bakry–Émery curvature-dimension inequalities and for the analysis of nonlinear PDEs on manifolds with weighted Ricci curvature lower bounds. Gradient comparison theorems, sharp spectral bounds, and corresponding rigidity or Liouville theorems are available under natural geometric constraints (Dai et al., 2020, Li et al., 2019).

Game-Theoretic Formulation: Recent developments link the viscosity solutions of weighted $p$ -Laplace equations to value functions of two-player stochastic games generalizing "tug-of-war," now with step probabilities and drifts determined by the weight's magnitude and gradient. The solution is characterized as the uniform limit (as step size $\varepsilon \to 0$ ) of a dynamic-programming value function involving weighted random walks and strategic moves—a framework extending Manfredi-Parviainen-Rossi's theory to the weighted case. The $p=\infty$ limiting equation involves a weighted infinity-Laplace operator and connects to drift-dominated tug-of-war (Aihara, 2024).

Inverse Problems: For equations of the form

$\nabla \cdot (\sigma(x)\nabla u + a(x)|\nabla u|^{p-2}\nabla u) = 0$

the recovery of the coefficients $(\sigma, a)$ from Dirichlet-to-Neumann data is achievable by a nonlinear linearization approach. The problem is well-posed in weighted spaces under ellipticity, and uniqueness/stability are established leveraging second-order Fréchet derivatives and CGO solutions analogously to Calderón-type results (Cârstea et al., 2020).

6. Numerical Analysis and Finite Element Approximation

Finite element discretization of the weighted $p$ -Laplace operator leverages the theory of weighted Orlicz–Sobolev spaces. With $A_1$ -class weights, the discrete variational problem is set in

$W^{1,p}_0(\omega, \Omega),\quad \|u\|_{W^{1,p}_0(\omega)} = \Bigl( \int_\Omega |u|^p + \omega|\nabla u|^p \, dx \Bigr)^{1/p}$

for strongly monotone and nonlinear operators. Céa-type best approximation and quasi-norm error bounds are derived using the $V$ -mapping $V(\nabla u) = |\nabla u|^{(p-2)/2}\nabla u$ and shifted $N$ -function techniques. Scott--Zhang and positivity-preserving interpolants yield a priori error estimates, including for the obstacle problem variant (Otarola et al., 2024).

Weighted $p(x)$ -Laplacian problems, variable exponent cases, and double-phase operators generalize the theory. Existence of weak solutions in double-weighted, variable-exponent Sobolev spaces is established for data in matching Muckenhoupt-type classes, leveraging compact embeddings and coercivity of the variational functional (Unal et al., 2019). Multiplicity results for systems and parametric families rely on critical point theory, Nehari manifold methods, and mountain-pass arguments, suitably adapted to the weighted and degenerate context (Repovš et al., 2021, Kumar et al., 2022).