Documentation

TauCeti.Analysis.PDE.MaximumPrinciple.Drift

The weak maximum principle with bounded drift #

Let Ω ⊆ ℝⁿ be open with finite measure, and let

L u = -∂ᵢ(aⁱʲ ∂ⱼu) + bⁱ ∂ᵢu + c u

have measurable, uniformly elliptic a with constants 0 < λ ≤ Λ, bounded measurable drift ‖b‖ ≤ β and bounded measurable c ≥ 0. This file proves the weak maximum principle for weak subsolutions of L: if u ∈ H¹(Ω) satisfies a(u, v) ≤ 0 for every nonnegative v ∈ H¹₀(Ω), and u ≤ k on ∂Ω for a constant k ≥ 0, in the sense that (u - k)⁺ ∈ H¹₀(Ω), then u ≤ k almost everywhere in Ω.

No smallness of the drift is needed, and the energy form need not be coercive. By contrast, TauCeti.PDE.UniformlyEllipticOn.value_le_const_of_small_drift_of_poincare requires β P < λ for a Poincaré constant P, which makes the energy form coercive on H¹₀(Ω).

The argument is Gilbarg–Trudinger's. For l ≥ k put v = (u - l)⁺. Testing against v, ellipticity and c ≥ 0 give λ ‖∇v‖² ≤ β ∫ |∇v| v, and the integrand on the right lives on the set Γₗ = {u > l, ∇u ≠ 0}. Hölder's inequality on Γₗ and a Sobolev inequality ‖v‖_q ≤ S ‖∇v‖₂, q > 2, then give

λ ≤ β S |Γₗ|^{1/2 - 1/q}

whenever u > l on a set of positive measure. As l increases to the essential supremum of u, the sets Γₗ shrink to a subset of {u = sup u, ∇u ≠ 0} (or to the empty set if u is essentially unbounded), which is null because the weak gradient vanishes almost everywhere on level sets. So |Γₗ| → 0, which is impossible if the essential supremum exceeds k.

The Sobolev inequality is a hypothesis of TauCeti.PDE.UniformlyEllipticOn.value_le_const_of_forall_eLpNorm_le, which does not need Ω to have finite measure, since for l > k ≥ 0 the sets {u > l} do. On a set of finite measure the Sobolev inequality holds in every dimension (TauCeti.W1p.eLpNorm_value_le_mul_measure_rpow_mul_enorm_gradient_of_lt), which gives the unconditional form TauCeti.PDE.UniformlyEllipticOn.value_le_const_of_measure_ne_top.

The comparison principle follows, and with it uniqueness of weak solutions of the Dirichlet problem L u = f, u ∈ H¹₀(Ω). On a bounded domain, uniqueness feeds the Fredholm alternative TauCeti.PDE.fredholmAlternative_isWeakSolutionDirichlet_sub_const: although the energy form of L need not be coercive when the drift is large, the Dirichlet problem is uniquely solvable for every f ∈ L²(Ω) as soon as c ≥ 0.

Main declarations #

References #

theorem TauCeti.PDE.UniformlyEllipticOn.value_le_const_of_forall_eLpNorm_le {ι : Type u_1} [Fintype ι] {mu : MeasureTheory.Measure (EuclideanSpace ℝ ι)} [mu.IsAddHaarMeasure] {Omega : TopologicalSpace.Opens (EuclideanSpace ℝ ι)} {a : EuclideanSpace ℝ ι → Matrix ι ι ℝ} {b : EuclideanSpace ℝ ι → EuclideanSpace ℝ ι} {c : EuclideanSpace ℝ ι → ℝ} {lam Lam beta gamma : ℝ} [DecidableEq ι] (h : UniformlyEllipticOn (↑Omega) a lam Lam) (ha : MeasureTheory.AEStronglyMeasurable a (mu.restrict ↑Omega)) (hb : MeasureTheory.AEStronglyMeasurable b (mu.restrict ↑Omega)) (hc : MeasureTheory.AEStronglyMeasurable c (mu.restrict ↑Omega)) (hb_bound : ∀ x ∈ Omega, ‖b x‖ ≤ beta) (hc_bound : ∀ x ∈ Omega, ‖c x‖ ≤ gamma) (hc_nonneg : ∀ x ∈ Omega, 0 ≤ c x) {q : ENNReal} (hq : 2 < q) {S : NNReal} (hsob : ∀ v ∈ w1p0Submodule mu Omega 2, MeasureTheory.eLpNorm (↑↑(W1p.value v)) q (mu.restrict ↑Omega) ≤ ↑S * ‖W1p.gradient v‖ₑ) {k : ℝ} (hk : 0 ≤ k) {u : ↥(W1p mu Omega 2)} (hboundary : W1p.posPartAbove value_le_const_of_forall_eLpNorm_le._proof_1 hk u ∈ w1p0Submodule mu Omega 2) (hu : ∀ (v : ↥(W1p0 mu Omega 2)), (∀ᵐ (x : EuclideanSpace ℝ ι) ∂mu.restrict ↑Omega, 0 ≤ ↑↑(W1p.value ↑v) x) → energyFormH1 a b c u ↑v ≤ 0) :
∀ᵐ (x : EuclideanSpace ℝ ι) ∂mu.restrict ↑Omega, ↑↑(W1p.value u) x ≤ k

The weak maximum principle with bounded drift, under a Sobolev inequality. Let Ω ⊆ ℝⁿ be open, and suppose ‖w‖_q ≤ S ‖∇w‖₂ on H¹₀(Ω) for some exponent q > 2. Let L u = -∂ᵢ(aⁱʲ ∂ⱼu) + bⁱ ∂ᵢu + c u have measurable, uniformly elliptic a with constants λ, Λ, measurable b with ‖b‖ ≤ β and measurable c with 0 ≤ c ≤ γ. If u ∈ H¹(Ω) is a weak subsolution, a(u, v) ≤ 0 for every nonnegative v ∈ H¹₀(Ω), and (u - k)⁺ ∈ H¹₀(Ω) for some k ≥ 0, then u ≤ k almost everywhere. There is no smallness condition on β, and Ω need not have finite measure: the superlevel sets {u > l}, l > k ≥ 0, of u ∈ L² already do.

theorem TauCeti.PDE.UniformlyEllipticOn.value_le_const_of_measure_ne_top {ι : Type u_1} [Fintype ι] {mu : MeasureTheory.Measure (EuclideanSpace ℝ ι)} [mu.IsAddHaarMeasure] {Omega : TopologicalSpace.Opens (EuclideanSpace ℝ ι)} {a : EuclideanSpace ℝ ι → Matrix ι ι ℝ} {b : EuclideanSpace ℝ ι → EuclideanSpace ℝ ι} {c : EuclideanSpace ℝ ι → ℝ} {lam Lam beta gamma : ℝ} [DecidableEq ι] [Nonempty ι] (h : UniformlyEllipticOn (↑Omega) a lam Lam) (ha : MeasureTheory.AEStronglyMeasurable a (mu.restrict ↑Omega)) (hb : MeasureTheory.AEStronglyMeasurable b (mu.restrict ↑Omega)) (hc : MeasureTheory.AEStronglyMeasurable c (mu.restrict ↑Omega)) (hb_bound : ∀ x ∈ Omega, ‖b x‖ ≤ beta) (hc_bound : ∀ x ∈ Omega, ‖c x‖ ≤ gamma) (hc_nonneg : ∀ x ∈ Omega, 0 ≤ c x) (hOmega : mu ↑Omega ≠ ⊤) {k : ℝ} (hk : 0 ≤ k) {u : ↥(W1p mu Omega 2)} (hboundary : W1p.posPartAbove value_le_const_of_forall_eLpNorm_le._proof_1 hk u ∈ w1p0Submodule mu Omega 2) (hu : ∀ (v : ↥(W1p0 mu Omega 2)), (∀ᵐ (x : EuclideanSpace ℝ ι) ∂mu.restrict ↑Omega, 0 ≤ ↑↑(W1p.value ↑v) x) → energyFormH1 a b c u ↑v ≤ 0) :
∀ᵐ (x : EuclideanSpace ℝ ι) ∂mu.restrict ↑Omega, ↑↑(W1p.value u) x ≤ k

The weak maximum principle with bounded drift. Let Ω ⊆ ℝⁿ, n ≥ 1, be open with finite measure, and let L u = -∂ᵢ(aⁱʲ ∂ⱼu) + bⁱ ∂ᵢu + c u have measurable, uniformly elliptic a with constants λ, Λ, measurable b with ‖b‖ ≤ β and measurable c with 0 ≤ c ≤ γ. If u ∈ H¹(Ω) is a weak subsolution, a(u, v) ≤ 0 for every nonnegative v ∈ H¹₀(Ω), and (u - k)⁺ ∈ H¹₀(Ω) for some k ≥ 0, then u ≤ k almost everywhere. There is no smallness condition on β and no regularity condition on ∂Ω.

theorem TauCeti.PDE.UniformlyEllipticOn.value_le_of_forall_eLpNorm_le {ι : Type u_1} [Fintype ι] {mu : MeasureTheory.Measure (EuclideanSpace ℝ ι)} [mu.IsAddHaarMeasure] {Omega : TopologicalSpace.Opens (EuclideanSpace ℝ ι)} {a : EuclideanSpace ℝ ι → Matrix ι ι ℝ} {b : EuclideanSpace ℝ ι → EuclideanSpace ℝ ι} {c : EuclideanSpace ℝ ι → ℝ} {lam Lam beta gamma : ℝ} [DecidableEq ι] (h : UniformlyEllipticOn (↑Omega) a lam Lam) (ha : MeasureTheory.AEStronglyMeasurable a (mu.restrict ↑Omega)) (hb : MeasureTheory.AEStronglyMeasurable b (mu.restrict ↑Omega)) (hc : MeasureTheory.AEStronglyMeasurable c (mu.restrict ↑Omega)) (hb_bound : ∀ x ∈ Omega, ‖b x‖ ≤ beta) (hc_bound : ∀ x ∈ Omega, ‖c x‖ ≤ gamma) (hc_nonneg : ∀ x ∈ Omega, 0 ≤ c x) {q : ENNReal} (hq : 2 < q) {S : NNReal} (hsob : ∀ v ∈ w1p0Submodule mu Omega 2, MeasureTheory.eLpNorm (↑↑(W1p.value v)) q (mu.restrict ↑Omega) ≤ ↑S * ‖W1p.gradient v‖ₑ) {u v : ↥(W1p mu Omega 2)} (hboundary : W1p.posPart value_le_const_of_forall_eLpNorm_le._proof_1 (u - v) ∈ w1p0Submodule mu Omega 2) (huv : ∀ (w : ↥(W1p0 mu Omega 2)), (∀ᵐ (x : EuclideanSpace ℝ ι) ∂mu.restrict ↑Omega, 0 ≤ ↑↑(W1p.value ↑w) x) → energyFormH1 a b c u ↑w ≤ energyFormH1 a b c v ↑w) :
∀ᵐ (x : EuclideanSpace ℝ ι) ∂mu.restrict ↑Omega, ↑↑(W1p.value u) x ≤ ↑↑(W1p.value v) x

The comparison principle with bounded drift, under a Sobolev inequality. Under the hypotheses of TauCeti.PDE.UniformlyEllipticOn.value_le_const_of_forall_eLpNorm_le, if a(u, w) ≤ a(v, w) for every nonnegative w ∈ H¹₀(Ω), that is L u ≤ L v weakly, and (u - v)⁺ ∈ H¹₀(Ω), then u ≤ v almost everywhere.

theorem TauCeti.PDE.UniformlyEllipticOn.value_le_of_measure_ne_top {ι : Type u_1} [Fintype ι] {mu : MeasureTheory.Measure (EuclideanSpace ℝ ι)} [mu.IsAddHaarMeasure] {Omega : TopologicalSpace.Opens (EuclideanSpace ℝ ι)} {a : EuclideanSpace ℝ ι → Matrix ι ι ℝ} {b : EuclideanSpace ℝ ι → EuclideanSpace ℝ ι} {c : EuclideanSpace ℝ ι → ℝ} {lam Lam beta gamma : ℝ} [DecidableEq ι] [Nonempty ι] (h : UniformlyEllipticOn (↑Omega) a lam Lam) (ha : MeasureTheory.AEStronglyMeasurable a (mu.restrict ↑Omega)) (hb : MeasureTheory.AEStronglyMeasurable b (mu.restrict ↑Omega)) (hc : MeasureTheory.AEStronglyMeasurable c (mu.restrict ↑Omega)) (hb_bound : ∀ x ∈ Omega, ‖b x‖ ≤ beta) (hc_bound : ∀ x ∈ Omega, ‖c x‖ ≤ gamma) (hc_nonneg : ∀ x ∈ Omega, 0 ≤ c x) (hOmega : mu ↑Omega ≠ ⊤) {u v : ↥(W1p mu Omega 2)} (hboundary : W1p.posPart value_le_const_of_forall_eLpNorm_le._proof_1 (u - v) ∈ w1p0Submodule mu Omega 2) (huv : ∀ (w : ↥(W1p0 mu Omega 2)), (∀ᵐ (x : EuclideanSpace ℝ ι) ∂mu.restrict ↑Omega, 0 ≤ ↑↑(W1p.value ↑w) x) → energyFormH1 a b c u ↑w ≤ energyFormH1 a b c v ↑w) :
∀ᵐ (x : EuclideanSpace ℝ ι) ∂mu.restrict ↑Omega, ↑↑(W1p.value u) x ≤ ↑↑(W1p.value v) x

The comparison principle with bounded drift. Under the hypotheses of TauCeti.PDE.UniformlyEllipticOn.value_le_const_of_measure_ne_top, if a(u, w) ≤ a(v, w) for every nonnegative w ∈ H¹₀(Ω), that is L u ≤ L v weakly, and (u - v)⁺ ∈ H¹₀(Ω), then u ≤ v almost everywhere.

theorem TauCeti.PDE.UniformlyEllipticOn.eq_of_isWeakSolutionDirichlet_of_mass_nonneg_of_forall_eLpNorm_le {ι : Type u_1} [Fintype ι] {mu : MeasureTheory.Measure (EuclideanSpace ℝ ι)} [mu.IsAddHaarMeasure] {Omega : TopologicalSpace.Opens (EuclideanSpace ℝ ι)} {a : EuclideanSpace ℝ ι → Matrix ι ι ℝ} {b : EuclideanSpace ℝ ι → EuclideanSpace ℝ ι} {c : EuclideanSpace ℝ ι → ℝ} {lam Lam beta gamma : ℝ} [DecidableEq ι] (h : UniformlyEllipticOn (↑Omega) a lam Lam) (ha : MeasureTheory.AEStronglyMeasurable a (mu.restrict ↑Omega)) (hb : MeasureTheory.AEStronglyMeasurable b (mu.restrict ↑Omega)) (hc : MeasureTheory.AEStronglyMeasurable c (mu.restrict ↑Omega)) (hb_bound : ∀ x ∈ Omega, ‖b x‖ ≤ beta) (hc_bound : ∀ x ∈ Omega, ‖c x‖ ≤ gamma) (hc_nonneg : ∀ x ∈ Omega, 0 ≤ c x) {q : ENNReal} (hq : 2 < q) {S : NNReal} (hsob : ∀ v ∈ w1p0Submodule mu Omega 2, MeasureTheory.eLpNorm (↑↑(W1p.value v)) q (mu.restrict ↑Omega) ≤ ↑S * ‖W1p.gradient v‖ₑ) {f : ↥(MeasureTheory.Lp ℝ 2 (mu.restrict ↑Omega))} {u v : ↥(W1p0 mu Omega 2)} (hu : IsWeakSolutionDirichlet a b c f u) (hv : IsWeakSolutionDirichlet a b c f v) :
u = v

Uniqueness for the Dirichlet problem with bounded drift, under a Sobolev inequality. Under the hypotheses of TauCeti.PDE.UniformlyEllipticOn.value_le_const_of_forall_eLpNorm_le, two weak solutions u, v ∈ H¹₀(Ω) of L u = f are equal.

theorem TauCeti.PDE.UniformlyEllipticOn.eq_of_isWeakSolutionDirichlet_of_mass_nonneg_of_measure_ne_top {ι : Type u_1} [Fintype ι] {mu : MeasureTheory.Measure (EuclideanSpace ℝ ι)} [mu.IsAddHaarMeasure] {Omega : TopologicalSpace.Opens (EuclideanSpace ℝ ι)} {a : EuclideanSpace ℝ ι → Matrix ι ι ℝ} {b : EuclideanSpace ℝ ι → EuclideanSpace ℝ ι} {c : EuclideanSpace ℝ ι → ℝ} {lam Lam beta gamma : ℝ} [DecidableEq ι] [Nonempty ι] (h : UniformlyEllipticOn (↑Omega) a lam Lam) (ha : MeasureTheory.AEStronglyMeasurable a (mu.restrict ↑Omega)) (hb : MeasureTheory.AEStronglyMeasurable b (mu.restrict ↑Omega)) (hc : MeasureTheory.AEStronglyMeasurable c (mu.restrict ↑Omega)) (hb_bound : ∀ x ∈ Omega, ‖b x‖ ≤ beta) (hc_bound : ∀ x ∈ Omega, ‖c x‖ ≤ gamma) (hc_nonneg : ∀ x ∈ Omega, 0 ≤ c x) (hOmega : mu ↑Omega ≠ ⊤) {f : ↥(MeasureTheory.Lp ℝ 2 (mu.restrict ↑Omega))} {u v : ↥(W1p0 mu Omega 2)} (hu : IsWeakSolutionDirichlet a b c f u) (hv : IsWeakSolutionDirichlet a b c f v) :
u = v

Uniqueness for the Dirichlet problem with bounded drift. Under the hypotheses of TauCeti.PDE.UniformlyEllipticOn.value_le_const_of_measure_ne_top, two weak solutions u, v ∈ H¹₀(Ω) of L u = f are equal.

theorem TauCeti.PDE.UniformlyEllipticOn.existsUnique_isWeakSolutionDirichlet_of_mass_nonneg {ι : Type u_1} [Fintype ι] {mu : MeasureTheory.Measure (EuclideanSpace ℝ ι)} [mu.IsAddHaarMeasure] {Omega : TopologicalSpace.Opens (EuclideanSpace ℝ ι)} {a : EuclideanSpace ℝ ι → Matrix ι ι ℝ} {b : EuclideanSpace ℝ ι → EuclideanSpace ℝ ι} {c : EuclideanSpace ℝ ι → ℝ} {lam Lam beta gamma : ℝ} [DecidableEq ι] [Nonempty ι] (h : UniformlyEllipticOn (↑Omega) a lam Lam) (ha : MeasureTheory.AEStronglyMeasurable a (mu.restrict ↑Omega)) (hb : MeasureTheory.AEStronglyMeasurable b (mu.restrict ↑Omega)) (hc : MeasureTheory.AEStronglyMeasurable c (mu.restrict ↑Omega)) (hb_bound : ∀ x ∈ Omega, ‖b x‖ ≤ beta) (hc_bound : ∀ x ∈ Omega, ‖c x‖ ≤ gamma) (hc_nonneg : ∀ x ∈ Omega, 0 ≤ c x) (hOmega : Bornology.IsBounded ↑Omega) (f : ↥(MeasureTheory.Lp ℝ 2 (mu.restrict ↑Omega))) :
∃! u : ↥(W1p0 mu Omega 2), IsWeakSolutionDirichlet a b c f u

Existence and uniqueness for the Dirichlet problem with c ≥ 0. Let Ω ⊆ ℝⁿ, n ≥ 1, be a bounded open set, and let L u = -∂ᵢ(aⁱʲ ∂ⱼu) + bⁱ ∂ᵢu + c u have measurable, uniformly elliptic a, measurable b with ‖b‖ ≤ β and measurable c with 0 ≤ c ≤ γ. Then for every f ∈ L²(Ω) there is exactly one u ∈ H¹₀(Ω) with a(u, v) = ∫_Ω f v for all v ∈ H¹₀(Ω). Since β is arbitrary, the energy form of L need not be coercive.