Documentation

TauCeti.Analysis.Sobolev.Trace.Reflection.Basic

Even reflection across a hyperplane preserves weak differentiability #

Let H = normalHalfSpace a = {x | a < x.fst} in the Euclidean product ℝ × E, and let ρ (t, y) = (2a - t, y) be the reflection in its boundary hyperplane {a} × E. If w is weakly differentiable on H with weak gradient G, and w and G are locally integrable up to the boundary (we extend them by zero off H), then the even reflection

x ↦ w x + w (ρ x)

is weakly differentiable on the whole space, with weak gradient x ↦ G x + R (G (ρ x)), where R (s, z) = (-s, z) is the linear part of ρ. No boundary condition on w is needed: the reflection is what makes the two one-sided traces match. This is the analytic content of the extension operator W^{1,p}(H) → W^{1,p}(ℝ × E) constructed in TauCeti.Analysis.Sobolev.Trace.Extension.

The proof follows H. Brezis, Functional Analysis, Sobolev Spaces and Partial Differential Equations, Lemma 9.2. A test function φ on the whole space is tested against the even reflection; after the change of variables x ↦ ρ x this becomes the pairing of w on H with χ = φ ± φ ∘ ρ, where the sign is + for directions tangent to the boundary and - for the normal direction. Such a χ is not supported inside H, so it is multiplied by cutoffs η_n (x.fst) vanishing near the boundary. In tangential directions the cutoff commutes with differentiation. In the normal direction χ vanishes on the boundary, so |χ| ≤ M (x.fst - a) there, which exactly compensates the size O(n) of η_n' on its support of width O(1/n). Dominated convergence then passes to the limit.

Main declarations #

References #

H. Brezis, Functional Analysis, Sobolev Spaces and Partial Differential Equations, Lemma 9.2.

The reflection in the boundary hyperplane #

The linear reflection (t, y) ↦ (-t, y) of the Euclidean product ℝ × E in the hyperplane {0} × E, the kernel of the normal coordinate. It transforms the gradient of a function under TauCeti.normalReflection.

Instances For
    @[simp]

    The linear normal reflection negates the normal coordinate.

    @[simp]

    The linear normal reflection is its own inverse.

    The reflection R reverses the normal direction ν = (1, 0).

    The reflection R fixes every tangential direction τ y = (0, y).

    Pairing w - R w' with the normal direction ν = (1, 0) adds the normal components of w and w'.

    Pairing w - R w' with a tangential direction τ y = (0, y) subtracts the tangential components of w' from those of w.

    The reflection (t, y) ↦ (2a - t, y) of the Euclidean product ℝ × E in the hyperplane {a} × E, the boundary of TauCeti.normalHalfSpace a. Its linear part is TauCeti.normalLinearReflection.

    Instances For
      @[simp]

      The normal reflection in {a} × E replaces the normal coordinate t by 2a - t.

      The normal reflection is an involution.

      @[simp]

      The normal reflection is its own inverse.

      The fixed points of the normal reflection in {a} × E are the points of that hyperplane.

      The normal reflection in {a} × E exchanges the open half-spaces on either side.

      The normal reflection in {a} × E is its linear part TauCeti.normalLinearReflection followed by the translation by (2a, 0).

      The derivative of the normal reflection is its linear part.

      The normal reflection is smooth.

      theorem TauCeti.gradient_normalReflection_of_forall_eq_neg {E : Type u_1} [NormedAddCommGroup E] [InnerProductSpace ℝ E] [FiniteDimensional ℝ E] {a : ℝ} {U : Set (WithLp 2 (ℝ × E))} (hU : IsOpen U) {h : WithLp 2 (ℝ × E) → ℝ} (hodd : ∀ x ∈ U, h ((normalReflection E a) x) = -h x) {x : WithLp 2 (ℝ × E)} (hx : x ∈ U) (hd : DifferentiableAt ℝ h ((normalReflection E a) x)) :

      Gradients of odd functions. Let h be odd under the reflection ρ in {a} × E on an open set U, that is h (ρ x) = -h x for x ∈ U. If h is differentiable at ρ x, then its gradient transforms by minus the linear part R of ρ: ∇h (ρ x) = -R (∇h x). The normal component of ∇h is therefore even under ρ and its tangential components are odd.

      theorem TauCeti.mapsTo_normalReflection_ball {E : Type u_1} [NormedAddCommGroup E] [InnerProductSpace ℝ E] [FiniteDimensional ℝ E] {a : ℝ} {x₀ : WithLp 2 (ℝ × E)} (hx₀ : x₀.fst = a) (r : ℝ) :
      Set.MapsTo (⇑(normalReflection E a)) (Metric.ball x₀ r) (Metric.ball x₀ r)

      The normal reflection in {a} × E maps every ball centred on that hyperplane into itself.

      Change of variables by the normal reflection: integrals against Lebesgue measure are invariant under it.

      Precomposition with the normal reflection preserves local integrability.

      A function invariant under the normal reflection ρ on a measurable set U mapped into itself by ρ is integrable on U as soon as it is integrable on the half U ∩ H.

      theorem TauCeti.setIntegral_eq_two_smul_setIntegral_inter_normalHalfSpace {E : Type u_1} [NormedAddCommGroup E] [InnerProductSpace ℝ E] [FiniteDimensional ℝ E] [MeasurableSpace E] [BorelSpace E] {a : ℝ} {U : Set (WithLp 2 (ℝ × E))} (hUm : MeasurableSet U) (hU : Set.MapsTo (⇑(normalReflection E a)) U U) {F : Type u_2} [NormedAddCommGroup F] [NormedSpace ℝ F] {g : WithLp 2 (ℝ × E) → F} (hg : ∀ x ∈ U, g ((normalReflection E a) x) = g x) (hint : MeasureTheory.IntegrableOn g (U ∩ ↑(normalHalfSpace a)) MeasureTheory.volume) :
      ∫ (x : WithLp 2 (ℝ × E)) in U, g x = 2 • ∫ (x : WithLp 2 (ℝ × E)) in U ∩ ↑(normalHalfSpace a), g x

      Integrals of reflection-invariant functions. If U is a measurable set mapped into itself by the normal reflection ρ in {a} × E, for instance a ball centred on that hyperplane, and g is integrable on U ∩ H and invariant under ρ on U, then the integral of g over U is twice its integral over U ∩ H.

      theorem TauCeti.ae_clm_eq_of_forall_inner_normal_tangent {E : Type u_1} [NormedAddCommGroup E] [InnerProductSpace ℝ E] [FiniteDimensional ℝ E] [MeasurableSpace E] {Φ Ψ : WithLp 2 (ℝ × E) → WithLp 2 (ℝ × E) →L[ℝ] WithLp 2 (ℝ × E)} {μ : MeasureTheory.Measure (WithLp 2 (ℝ × E))} (h1 : ∀ (v : WithLp 2 (ℝ × E)), ∀ᵐ (x : WithLp 2 (ℝ × E)) ∂μ, inner ℝ ((Φ x) v) (WithLp.toLp 2 (1, 0)) = inner ℝ ((Ψ x) v) (WithLp.toLp 2 (1, 0))) (h2 : ∀ (y : E), ∀ᵐ (x : WithLp 2 (ℝ × E)) ∂μ, inner ℝ ((Φ x) (WithLp.toLp 2 (1, 0))) (WithLp.toLp 2 (0, y)) = inner ℝ ((Ψ x) (WithLp.toLp 2 (1, 0))) (WithLp.toLp 2 (0, y))) (h3 : ∀ (y z : E), ∀ᵐ (x : WithLp 2 (ℝ × E)) ∂μ, inner ℝ ((Φ x) (WithLp.toLp 2 (0, z))) (WithLp.toLp 2 (0, y)) = inner ℝ ((Ψ x) (WithLp.toLp 2 (0, z))) (WithLp.toLp 2 (0, y))) :
      ∀ᵐ (x : WithLp 2 (ℝ × E)) ∂μ, Φ x = Ψ x

      Two continuous linear endomorphisms agree as soon as their pairings agree with the normal direction on every vector and with the tangential directions on the normal and tangential directions.

      Cutoffs vanishing near the boundary hyperplane #

      noncomputable def TauCeti.normalCutoff {E : Type u_1} (a c : ℝ) (x : WithLp 2 (ℝ × E)) :

      The boundary-layer cutoff x ↦ σ(c (x.fst - a) - 1) of slope c in the normal coordinate, where σ is Real.smoothTransition. For 0 < c it vanishes when x.fst ≤ a + 1 / c and is one when x.fst ≥ a + 2 / c, so it is supported a positive distance inside the half-space TauCeti.normalHalfSpace a.

      Instances For
        theorem TauCeti.normalCutoff_def {E : Type u_1} (a c : ℝ) (x : WithLp 2 (ℝ × E)) :
        normalCutoff a c x = (c * (x.fst - a) - 1).smoothTransition

        The defining formula of the boundary-layer cutoff.

        theorem TauCeti.normalCutoff_nonneg {E : Type u_1} (a c : ℝ) (x : WithLp 2 (ℝ × E)) :

        The boundary-layer cutoff is nonnegative.

        theorem TauCeti.normalCutoff_le_one {E : Type u_1} (a c : ℝ) (x : WithLp 2 (ℝ × E)) :

        The boundary-layer cutoff is at most one.

        @[simp]
        theorem TauCeti.normalCutoff_eq_zero {E : Type u_1} {a c : ℝ} {x : WithLp 2 (ℝ × E)} (hx : c * (x.fst - a) ≤ 1) :
        normalCutoff a c x = 0

        The boundary-layer cutoff is zero where c (x.fst - a) ≤ 1.

        @[simp]
        theorem TauCeti.normalCutoff_eq_one {E : Type u_1} {a c : ℝ} {x : WithLp 2 (ℝ × E)} (hx : 2 ≤ c * (x.fst - a)) :
        normalCutoff a c x = 1

        The boundary-layer cutoff is one where c (x.fst - a) ≥ 2.

        The boundary-layer cutoff is smooth.

        The derivative of the boundary-layer cutoff is σ'(c (x.fst - a) - 1) c times the normal coordinate.

        theorem TauCeti.tsupport_normalCutoff_subset {E : Type u_1} [NormedAddCommGroup E] [NormedSpace ℝ E] (a : ℝ) {c : ℝ} (hc : 0 < c) :

        For a positive slope c, the boundary-layer cutoff has topological support inside the half-space TauCeti.normalHalfSpace a.

        theorem TauCeti.exists_normalCutoff_bound {E : Type u_1} [NormedAddCommGroup E] [InnerProductSpace ℝ E] [CompleteSpace E] (a : ℝ) {c : ℝ} (hc : 0 ≤ c) :
        ∃ (M : ℝ), 0 ≤ M ∧ (∀ (x : WithLp 2 (ℝ × E)), |normalCutoff a c x| ≤ M) ∧ ∀ (x : WithLp 2 (ℝ × E)), ‖gradient (normalCutoff a c) x‖ ≤ M

        The boundary-layer cutoff and its gradient are bounded by a common constant.

        Integration by parts up to the boundary #

        theorem TauCeti.tendsto_integral_normalCutoff_mul {E : Type u_1} [NormedAddCommGroup E] [InnerProductSpace ℝ E] [FiniteDimensional ℝ E] [MeasurableSpace E] [BorelSpace E] {a : ℝ} {F : WithLp 2 (ℝ × E) → ℝ} (hF : MeasureTheory.Integrable F MeasureTheory.volume) (hF0 : ∀ x ∉ normalHalfSpace a, F x = 0) :
        Filter.Tendsto (fun (n : ℕ) => ∫ (x : WithLp 2 (ℝ × E)), normalCutoff a (↑n + 1) x * F x) Filter.atTop (nhds (∫ (x : WithLp 2 (ℝ × E)), F x))

        Multiplying by the boundary cutoffs of slope n + 1 does not change, in the limit n → ∞, the integral of an integrable function vanishing off the half-space.

        theorem TauCeti.tendsto_integral_mul_fderiv_normalCutoff_apply {E : Type u_1} [NormedAddCommGroup E] [InnerProductSpace ℝ E] [FiniteDimensional ℝ E] [MeasurableSpace E] [BorelSpace E] {a : ℝ} {G : WithLp 2 (ℝ × E) → WithLp 2 (ℝ × E)} (hG : MeasureTheory.LocallyIntegrable G MeasureTheory.volume) {χ : WithLp 2 (ℝ × E) → ℝ} (hχ : ContDiff ℝ (↑⊤) χ) (hχc : HasCompactSupport χ) (hχa : ∀ (y : E), χ (WithLp.toLp 2 (a, y)) = 0) :
        Filter.Tendsto (fun (n : ℕ) => ∫ (x : WithLp 2 (ℝ × E)), χ x * (fderiv ℝ (normalCutoff a (↑n + 1)) x) (G x)) Filter.atTop (nhds 0)

        The boundary-layer error vanishes in the limit. Let G be a locally integrable vector field, and let χ be smooth with compact support and vanishing on the boundary hyperplane {a} × E. Then ∫ χ ⟪∇η_n, G⟫ → 0 as n → ∞, where η_n is the boundary cutoff of slope n + 1, even though ∇η_n grows like n on the boundary layer.

        Even reflection across the boundary of a half-space preserves weak differentiability. Let w have weak gradient G on the half-space {x | a < x.fst}, with w and G vanishing off the half-space and locally integrable on the whole space (so integrable up to the boundary on bounded sets). Then x ↦ w x + w (ρ x), where ρ is the reflection in the boundary hyperplane, is weakly differentiable on the whole space with weak gradient x ↦ G x + R (G (ρ x)), R the linear part of ρ. No boundary condition on w is assumed.

        Extension by zero in a tangential direction. Let w have weak derivative w' on the half-space {x | a < x.fst} in a direction v tangent to the boundary hyperplane, with w and w' vanishing off the half-space and locally integrable on the whole space. Then w' is a weak derivative of w in the direction v on the whole space. No boundary condition on w is needed: integration by parts in a tangential direction produces no boundary term.

        theorem TauCeti.HasWeakFDerivOn.hasWeakLineDerivOn_sub_comp_normalReflection {E : Type u_1} [NormedAddCommGroup E] [InnerProductSpace ℝ E] [FiniteDimensional ℝ E] [MeasurableSpace E] [BorelSpace E] {a : ℝ} {w : WithLp 2 (ℝ × E) → ℝ} {G : WithLp 2 (ℝ × E) → WithLp 2 (ℝ × E)} (hw : HasWeakFDerivOn MeasureTheory.volume (normalHalfSpace a) w fun (x : WithLp 2 (ℝ × E)) => (innerSL ℝ) (G x)) (hwl : MeasureTheory.LocallyIntegrable w MeasureTheory.volume) (hGl : MeasureTheory.LocallyIntegrable G MeasureTheory.volume) (hw0 : ∀ x ∉ normalHalfSpace a, w x = 0) (hG0 : ∀ x ∉ normalHalfSpace a, G x = 0) {d : WithLp 2 (ℝ × E)} (hd : d.fst = 0) :
        HasWeakLineDerivOn MeasureTheory.volume ⊤ (fun (x : WithLp 2 (ℝ × E)) => w x - w ((normalReflection E a) x)) (fun (x : WithLp 2 (ℝ × E)) => inner ℝ (G x - (normalLinearReflection E) (G ((normalReflection E a) x))) d) d

        Odd reflection in a tangential direction. Let w have weak gradient G on the half-space {x | a < x.fst}, with w and G vanishing off the half-space and locally integrable on the whole space. Then, in every direction d tangent to the boundary hyperplane, the odd reflection x ↦ w x - w (ρ x) is weakly differentiable on the whole space, with weak derivative x ↦ ⟪G x - R (G (ρ x)), d⟫. No boundary condition on w is assumed; in the normal direction the odd reflection is weakly differentiable only when w has zero trace.

        Precomposition with the reflection on the whole space. If w has weak gradient G on the whole space, then w ∘ ρ has weak gradient x ↦ R (G (ρ x)), where ρ is the reflection in {a} × E and R its linear part.