Documentation

LeanPool.CIVAxisymmetric.Closure.CirculationBound

The circulation bound eq:aniso:circulation:bound #

The last step of lem:aniso:annulus. On the regular annulus {R_- < |x| < R_+} × (t_0, 0) the velocity is bounded, so for R_* ∈ (R_-, R_+) the circulation Γ = r u_θ is bounded on the lateral boundary ∂B(R_*) × (t_0, 0) by 2 R_* M; and u is smooth at time t_0, so Γ is bounded on B(R_*) there. Between those two data sets the axis maximum principle CIV.axisMaximumPrinciple (lem:aniso:axis) fills in the interior: Γ obeys eq:aniso:axis:pde with drift (u_r, u_z), exponent k = −1, no zeroth-order term and source r f_θ, and it vanishes on the axis, so it is bounded on the closed half-disc by its parabolic-boundary values plus the accumulated source R_* M_f (t − t_0).

Two remarks on the shape of the argument. First, the maximum principle is applied with its final time s set to the time at which the bound is read, so the constant it produces never mentions an auxiliary s; the manuscript's s ↑ 0 is the monotone step t − t_0 ≤ −t_0 = |t_0|. Second, Γ is invariant under the rotations about the axis (circulation_rotZ_invariant) and every point of Vec3 is such a rotation of a meridional point with r ≥ 0 (exists_rotZ_meridional), so the half-disc bound transfers to the whole ball B(R_*) — axis and inner ball included, which is what eq:aniso:circulation:bound asserts.

The data of the circulation instantiation of lem:aniso:axis #

The radial drift u_r, read on the meridional plane.

Equations
Instances For

    The vertical drift u_z, read on the meridional plane.

    Equations
    Instances For

      The source r f_θ of the circulation equation, read on the meridional plane.

      Equations
      Instances For
        theorem CIV.abs_circulationSource_le (f : CKN.Foundation.Parabolic.ParabolicPoint → CKN.Foundation.Parabolic.Vec3) (Mf : ℝ) (hMf : ∀ z ∈ unitCylinder, |f z 1| ≤ Mf) (Rstar tη s : ℝ) (hRpos : 0 < Rstar) (hR1 : Rstar < 1) (hs0 : s < 0) (htη1 : -1 < tη) (p : (ℝ × ℝ) × ℝ) :
        p ∈ axisDomain Rstar tη s → |circulationSource f p| ≤ Rstar * Mf

        The source of the circulation equation obeys |r f_θ| ≤ R_* M_f on the half-disc of radius R_*, the uniform source bound CIV.axisMaximumPrinciple asks for.

        The circulation bound on the meridional half-disc #

        theorem CIV.abs_circulationProfile_le_of_nonneg (u : CKN.Foundation.Parabolic.ParabolicPoint → CKN.Foundation.Parabolic.Vec3) (pr : CKN.Foundation.Parabolic.ParabolicPoint → ℝ) (f : CKN.Foundation.Parabolic.ParabolicPoint → CKN.Foundation.Parabolic.Vec3) (hsol : IsClassicalSolutionOn u pr f unitCylinder) (haxi : IsAxisymmetricOn u unitCylinder) (hfaxi : IsAxisymmetricOn f unitCylinder) (Mf : ℝ) (hMf : ∀ z ∈ unitCylinder, |f z 1| ≤ Mf) (Rstar : ℝ) (hRpos : 0 < Rstar) (hR1 : Rstar < 1) (tη : ℝ) (htη1 : -1 < tη) (htη2 : tη < 0) (Mb : ℝ) (hbdry : ∀ (x : CKN.Foundation.Parabolic.Vec3), CKN.Foundation.Parabolic.vec3EuclideanNorm x = Rstar → ∀ t ∈ Set.Ioo tη 0, |circulation u (x, t)| ≤ Mb) (Mi : ℝ) (hinit : ∀ x ∈ CKN.Foundation.Parabolic.vec3Ball 0 Rstar, |circulation u (x, tη)| ≤ Mi) (t : ℝ) (htη_lt : tη < t) (ht0 : t < 0) (r z : ℝ) (hr0 : 0 ≤ r) (hrz : r ^ 2 + z ^ 2 < Rstar ^ 2) :
        |circulationProfile u ((r, z), t)| ≤ max Mb Mi + Rstar * Mf * (t - tη)

        The circulation bound on the closed meridional half-disc r ≥ 0, by one application of the axis maximum principle CIV.axisMaximumPrinciple to the circulation profile with drift (u_r, u_z), exponent k = −1, zeroth-order coefficient γ = 0 and source r f_θ. The bound is read at the final time s := t, so the constant max (M_b, M_i) + R_* M_f (t − t_η) does not depend on any auxiliary final time: this is the manuscript's "neither depends on s".

        From the half-disc to the whole ball #

        theorem CIV.abs_circulation_le_of_boundary_data (u : CKN.Foundation.Parabolic.ParabolicPoint → CKN.Foundation.Parabolic.Vec3) (pr : CKN.Foundation.Parabolic.ParabolicPoint → ℝ) (f : CKN.Foundation.Parabolic.ParabolicPoint → CKN.Foundation.Parabolic.Vec3) (hsol : IsClassicalSolutionOn u pr f unitCylinder) (haxi : IsAxisymmetricOn u unitCylinder) (hfaxi : IsAxisymmetricOn f unitCylinder) (Mf : ℝ) (hMf : ∀ z ∈ unitCylinder, |f z 1| ≤ Mf) (Rstar : ℝ) (hRpos : 0 < Rstar) (hR1 : Rstar < 1) (tη : ℝ) (htη1 : -1 < tη) (htη2 : tη < 0) (Mb : ℝ) (hbdry : ∀ (x : CKN.Foundation.Parabolic.Vec3), CKN.Foundation.Parabolic.vec3EuclideanNorm x = Rstar → ∀ t ∈ Set.Ioo tη 0, |circulation u (x, t)| ≤ Mb) (Mi : ℝ) (hinit : ∀ x ∈ CKN.Foundation.Parabolic.vec3Ball 0 Rstar, |circulation u (x, tη)| ≤ Mi) (t : ℝ) :
        t ∈ Set.Ioo tη 0 → ∀ x ∈ CKN.Foundation.Parabolic.vec3Ball 0 Rstar, |circulation u (x, t)| ≤ max Mb Mi + Rstar * Mf * -tη

        The circulation bound on the whole ball B(R_*), axis and inner ball included: the circulation is invariant under the rotations about the axis (circulation_rotZ_invariant), and every point of Vec3 is such a rotation of a meridional point with r ≥ 0 (exists_rotZ_meridional), so the half-disc bound transfers verbatim. The final time enters only through t − t_η ≤ −t_η, which is the manuscript's M_f |t_0| and its s ↑ 0.

        The initial data at time t_η #

        theorem CIV.exists_bound_circulation_slice (u : CKN.Foundation.Parabolic.ParabolicPoint → CKN.Foundation.Parabolic.Vec3) (hu1 : ContDiffOn ℝ 1 (fun (w : CKN.Foundation.Parabolic.Vec3 × ℝ) => u w) unitCylinder) (Rstar : ℝ) (hRpos : 0 < Rstar) (hR1 : Rstar < 1) (τ : ℝ) (hτ : τ ∈ Set.Ioo (-1) 0) :
        ∃ (Mi : ℝ), ∀ x ∈ CKN.Foundation.Parabolic.vec3Ball 0 Rstar, |circulation u (x, τ)| ≤ Mi

        At any time of (−1, 0) the circulation is bounded on B(R_*) for R_* < 1: u is smooth there, so Γ is continuous on the compact closure of the ball. This is the manuscript's "the former are finite since u is smooth at time t_0".

        The azimuthal force bound of eq:interior:force:c-two #

        The order-0 case of CIV.ForceC2Bounded (eq:interior:force:c-two): one constant bounds the azimuthal component of the force on the whole cylinder.

        eq:aniso:circulation:bound #

        theorem CIV.exists_circulation_bound_ball (u : CKN.Foundation.Parabolic.ParabolicPoint → CKN.Foundation.Parabolic.Vec3) (pr : CKN.Foundation.Parabolic.ParabolicPoint → ℝ) (f : CKN.Foundation.Parabolic.ParabolicPoint → CKN.Foundation.Parabolic.Vec3) (hsol : IsClassicalSolutionOn u pr f unitCylinder) (haxi : IsAxisymmetricOn u unitCylinder) (hfaxi : IsAxisymmetricOn f unitCylinder) (Mf : ℝ) (hMf : ∀ z ∈ unitCylinder, |f z 1| ≤ Mf) (Rstar : ℝ) (hR : 0 < Rstar ∧ Rstar < 1) (t₀ : ℝ) (ht₀ : t₀ ∈ Set.Ioo (-1) 0) (Mb : ℝ) (hbdry : ∀ (x : CKN.Foundation.Parabolic.Vec3), CKN.Foundation.Parabolic.vec3EuclideanNorm x = Rstar → ∀ t ∈ Set.Ioo t₀ 0, |circulation u (x, t)| ≤ Mb) :
        ∃ (CΓ : ℝ), ∀ t ∈ Set.Ioo t₀ 0, ∀ x ∈ CKN.Foundation.Parabolic.vec3Ball 0 Rstar, |circulation u (x, t)| ≤ CΓ

        The circulation bound eq:aniso:circulation:bound on B(R_*) × (t_0, 0), given a bound M_b for Γ on the lateral boundary ∂B(R_*) × (t_0, 0). The initial constant is produced here, from the smoothness of u at time t_0.

        The circulation bound eq:aniso:circulation:bound of lem:aniso:annulus: from the regular annulus chosen in its proof — a bound M for |u| on {R_- < |x| < R_+} × (t_0, 0) — the circulation obeys eq:aniso:circulation:bound on B(R_*) × (t_0, 0) for every R_* ∈ (R_-, R_+), axis and inner ball included. The lateral data is |Γ| ≤ 2 R_* M on ∂B(R_*) (circulation_hbdry_of_annulus_bound); the initial data comes from the smoothness of u at time t_0; the interior is filled in by the maximum principle CIV.axisMaximumPrinciple applied to Γ with k = −1, γ = 0 and source r f_θ.

        theorem CIV.exists_regular_annulus_circulation_bound_ball (q : ℝ) (u : CKN.Foundation.Parabolic.ParabolicPoint → CKN.Foundation.Parabolic.Vec3) (Du : CKN.Foundation.Parabolic.ParabolicPoint → Fin 3 → CKN.Foundation.Parabolic.Vec3) (pr : CKN.Foundation.Parabolic.ParabolicPoint → ℝ) (f : CKN.Foundation.Parabolic.ParabolicPoint → CKN.Foundation.Parabolic.Vec3) (hsw : CKN.IsSuitableWeakSolution (CKN.Foundation.Parabolic.vec3Ball 0 1) (Set.Ioo (-1) 0) q u Du pr f) (henergy : GlobalEnergyClass u Du pr) (hsol : IsClassicalSolutionOn u pr f unitCylinder) (haxi : IsAxisymmetricOn u unitCylinder) (hfaxi : IsAxisymmetricOn f unitCylinder) (hMf : ForceC2Bounded f) (R₁ R₀ : ℝ) (hR₁ : 0 ≤ R₁) (hR₁R₀ : R₁ < R₀) (hR₀ : R₀ < 1) :
        ∃ (Rlo : ℝ) (Rhi : ℝ) (t₀ : ℝ), R₁ < Rlo ∧ Rlo < Rhi ∧ Rhi < R₀ ∧ t₀ ∈ Set.Ioo (-1) 0 ∧ ∀ Rstar ∈ Set.Ioo Rlo Rhi, ∃ (CΓ : ℝ), ∀ t ∈ Set.Ioo t₀ 0, ∀ x ∈ CKN.Foundation.Parabolic.vec3Ball 0 Rstar, |circulation u (x, t)| ≤ CΓ

        lem:aniso:annulus's circulation clause in full: for 0 ≤ R₁ < R₀ < 1 there are radii R₁ < R_- < R_+ < R₀ and a time t_0 ∈ (−1, 0) such that, for every R_* ∈ (R_-, R_+), the circulation is bounded on the whole cylinder B(R_*) × (t_0, 0).

        The circulation bound of lem:aniso:annulus #

        The circulation bound eq:aniso:circulation:bound in the proof of lem:aniso:annulus. Fix R_* ∈ (R_-, R_+). For s ∈ (t_0, 0) the set closure B(R_*) × [t_0, s] is a compact subset of Q on which u is smooth, so Γ = r u_θ is continuous there (continuousOn_circulationProfile), smooth off the axis (hreg_circulationProfile), and vanishes on the axis since |Γ| ≤ r |u| (circulationProfile_axis_eq_zero). The circulation equation eq:aniso:circulation:pde is of the form eq:aniso:axis:pde with k = −1, γ = 0 and source r f_θ with |r f_θ| ≤ R_* M_f (hpde_circulationProfile, abs_circulationSource_le, from eq:interior:force:c-two). lem:aniso:axis, in the form CIV.axisMaximumPrinciple, bounds |Γ| on closure B(R_*) × [t_0, s] by its values at time t_0 and on ∂B(R_*) × [t_0, s] plus M_f |t_0|; the former are finite since u is smooth at time t_0 (exists_bound_circulation_slice), the latter at most 2 R_* times the bound on u on the annulus (circulation_hbdry_of_annulus_bound); neither depends on s, and s ↑ 0 gives eq:aniso:circulation:bound.