Documentation

LeanPool.NavierStokesAndEuler.Euler.LinearDuhamelSobolevGevrey

The forward estimate in actual fixed-Sobolev external-word blocks #

The normalized Duhamel solution satisfies the constructed frozen equation. Its inverse at the base point is literally Id. Applying the finite base-order inverse estimate and then the external-word recurrence gives one factorial shift at the same input/output radius. Every constant is a fixed polynomial in the coefficient, data, propagator and time-length constants when q is fixed.

Quantitative bounds for the actual frozen forward equation #

The frozen coefficient has a polynomial tensor multiplier bound derived from the original coefficient and the H3 Green bound. Its fixed-Sobolev forcing block is bounded directly by the original initial/forcing blocks. No profile extremum, inverse amplitude, or raw weighted primitive is used.

A frozen bounded-operator equation for the actual forward solve #

At a chosen parameter x, Duhamel gives an exact equation with coefficient Id minus the fixed Green operator applied to the coefficient difference. Its coefficient at x is exactly Id. Thus the fixed-Sobolev inverse estimate can use the identity inverse; it never requires a norm for a profile-weighted raw time primitive.

@[instance_reducible]

Cache the standard NormedAddCommGroup (E →L[ℝ] E) instance to shorten typeclass synthesis.

Equations
Instances For
    @[instance_reducible]

    Cache the standard NormedSpace ℝ (E →L[ℝ] E) instance to shorten typeclass synthesis.

    Equations
    Instances For
      @[instance_reducible]

      Cache the standard NormedAddCommGroup C(Icc (0 : ℝ) T,E) instance to shorten typeclass synthesis.

      Equations
      Instances For
        @[instance_reducible]

        Cache the standard NormedSpace ℝ C(Icc (0 : ℝ) T,E) instance to shorten typeclass synthesis.

        Equations
        Instances For
          @[instance_reducible]

          Cache the standard NormedAddCommGroup C(Icc (0 : ℝ) T,E →L[ℝ] E) instance to shorten typeclass synthesis.

          Equations
          Instances For
            @[instance_reducible]

            Cache the standard NormedSpace ℝ C(Icc (0 : ℝ) T,E →L[ℝ] E) instance to shorten typeclass synthesis.

            Equations
            Instances For
              @[instance_reducible]

              Cache the standard NormedAddCommGroup (C(Icc (0 : ℝ) T,E) →L[ℝ] C(Icc (0 : ℝ) T,E)) instance to shorten typeclass synthesis.

              Equations
              Instances For
                @[instance_reducible]

                Cache the standard NormedSpace ℝ (C(Icc (0 : ℝ) T,E) →L[ℝ] C(Icc (0 : ℝ) T,E)) instance to shorten typeclass synthesis.

                Equations
                Instances For
                  noncomputable def EulerLinearDuhamel.frozenOperator {P : Type u_1} {E : Type u_2} [NormedAddCommGroup E] [NormedSpace E] [CompleteSpace E] (T : ) (hT : 0 T) (B : PC((Set.Icc 0 T), E →L[] E)) (U : (x : P) → Evolution T hT (B x)) (g : C((Set.Icc 0 T), )) (hg : ∀ (t : (Set.Icc 0 T)), 0 < g t) (x y : P) :
                  C((Set.Icc 0 T), E) →L[] C((Set.Icc 0 T), E)

                  The exact frozen coefficient, using the actual weighted Green operator.

                  Equations
                  • One or more equations did not get rendered due to their size.
                  Instances For
                    noncomputable def EulerLinearDuhamel.frozenForcing {P : Type u_1} {E : Type u_2} [NormedAddCommGroup E] [NormedSpace E] [CompleteSpace E] (T : ) (hT : 0 T) (B : PC((Set.Icc 0 T), E →L[] E)) (U : (x : P) → Evolution T hT (B x)) (g : C((Set.Icc 0 T), )) (hg : ∀ (t : (Set.Icc 0 T)), 0 < g t) (f : PC((Set.Icc 0 T), E)) (a₀ : PE) (x y : P) :
                    C((Set.Icc 0 T), E)

                    The exact transformed data under the fixed homogeneous and Green operators.

                    Equations
                    Instances For
                      theorem EulerLinearDuhamel.frozenOperator_self {P : Type u_1} {E : Type u_2} [NormedAddCommGroup E] [NormedSpace E] [CompleteSpace E] (T : ) (hT : 0 T) (B : PC((Set.Icc 0 T), E →L[] E)) (U : (x : P) → Evolution T hT (B x)) (g : C((Set.Icc 0 T), )) (hg : ∀ (t : (Set.Icc 0 T)), 0 < g t) (x : P) :

                      The frozen coefficient at its base parameter is the identity.

                      theorem EulerLinearDuhamel.frozenOperator_equation {P : Type u_1} {E : Type u_2} [NormedAddCommGroup E] [NormedSpace E] [CompleteSpace E] (T : ) (hT : 0 T) (B : PC((Set.Icc 0 T), E →L[] E)) (U : (x : P) → Evolution T hT (B x)) (g : C((Set.Icc 0 T), )) (hg : ∀ (t : (Set.Icc 0 T)), 0 < g t) (f : PC((Set.Icc 0 T), E)) (a₀ : PE) (x y : P) :
                      (frozenOperator T hT B U g hg x y) ((U y).weightedSolution g hg (f y) (a₀ y)) = frozenForcing T hT B U g hg f a₀ x y

                      The actual normalized Duhamel solution satisfies this bounded-operator equation.

                      theorem EulerLinearDuhamel.frozenOperator_contDiff {P : Type u_1} {E : Type u_2} [NormedAddCommGroup P] [NormedSpace P] [NormedAddCommGroup E] [NormedSpace E] [CompleteSpace E] (T : ) (hT : 0 T) (B : PC((Set.Icc 0 T), E →L[] E)) (U : (x : P) → Evolution T hT (B x)) (g : C((Set.Icc 0 T), )) (hg : ∀ (t : (Set.Icc 0 T)), 0 < g t) (hB : ContDiff (↑) B) (x : P) :
                      ContDiff (↑) (frozenOperator T hT B U g hg x)

                      The frozen coefficient is genuinely smooth in the translated coefficients.

                      theorem EulerLinearDuhamel.frozenForcing_contDiff {P : Type u_1} {E : Type u_2} [NormedAddCommGroup P] [NormedSpace P] [NormedAddCommGroup E] [NormedSpace E] [CompleteSpace E] (T : ) (hT : 0 T) (B : PC((Set.Icc 0 T), E →L[] E)) (U : (x : P) → Evolution T hT (B x)) (g : C((Set.Icc 0 T), )) (hg : ∀ (t : (Set.Icc 0 T)), 0 < g t) (f : PC((Set.Icc 0 T), E)) (a₀ : PE) (hf : ContDiff (↑) f) (ha₀ : ContDiff (↑) a₀) (x : P) :
                      ContDiff (↑) (frozenForcing T hT B U g hg f a₀ x)

                      The transformed data retain the actual parameter smoothness of the original data.

                      @[instance_reducible]

                      Cache the standard NormedAddCommGroup (E →L[ℝ] E) instance to shorten typeclass synthesis.

                      Equations
                      Instances For
                        @[instance_reducible]

                        Cache the standard NormedSpace ℝ (E →L[ℝ] E) instance to shorten typeclass synthesis.

                        Equations
                        Instances For
                          @[instance_reducible]

                          Cache the standard NormedAddCommGroup C(Icc (0 : ℝ) T,E) instance to shorten typeclass synthesis.

                          Equations
                          Instances For
                            @[instance_reducible]

                            Cache the standard NormedSpace ℝ C(Icc (0 : ℝ) T,E) instance to shorten typeclass synthesis.

                            Equations
                            Instances For
                              @[instance_reducible]

                              Cache the standard NormedAddCommGroup C(Icc (0 : ℝ) T,E →L[ℝ] E) instance to shorten typeclass synthesis.

                              Equations
                              Instances For
                                @[instance_reducible]

                                Cache the standard NormedSpace ℝ C(Icc (0 : ℝ) T,E →L[ℝ] E) instance to shorten typeclass synthesis.

                                Equations
                                Instances For
                                  @[instance_reducible]

                                  Cache the standard NormedAddCommGroup (C(Icc (0 : ℝ) T,E) →L[ℝ] C(Icc (0 : ℝ) T,E)) instance to shorten typeclass synthesis.

                                  Equations
                                  Instances For
                                    @[instance_reducible]

                                    Cache the standard NormedSpace ℝ (C(Icc (0 : ℝ) T,E) →L[ℝ] C(Icc (0 : ℝ) T,E)) instance to shorten typeclass synthesis.

                                    Equations
                                    Instances For

                                      The frozen coefficient amplitude is polynomial in the original coefficient and H3 constant.

                                      Equations
                                      Instances For
                                        theorem EulerLinearDuhamel.frozenOperator_bound {P : Type u_1} {E : Type u_2} [NormedAddCommGroup P] [NormedSpace P] [NormedAddCommGroup E] [NormedSpace E] [CompleteSpace E] (T : ) (hT : 0 T) (B : PC((Set.Icc 0 T), E →L[] E)) (U : (x : P) → Evolution T hT (B x)) (g : C((Set.Icc 0 T), )) (hg : ∀ (t : (Set.Icc 0 T)), 0 < g t) (hB : ContDiff (↑) B) (hg₀ : g 0, = 1) (C CB Rc : ) (hC : 0 C) (hCB : 0 CB) (hRc : 0 Rc) (hBb : ∀ (n : ) (y : P), iteratedFDeriv n B y CB * EulerGevrey.majorant Rc 0 n) (x : P) (hU : ∀ (t s : (Set.Icc 0 T)), s t(U x).propagator t s C * g t / g s) (n : ) (y : P) :

                                        Actual derivatives of the frozen coefficient have the stated polynomial bound.

                                        theorem EulerLinearDuhamel.frozenForcing_block_bound {P : Type u_1} {E : Type u_2} [NormedAddCommGroup P] [NormedSpace P] [NormedAddCommGroup E] [NormedSpace E] [CompleteSpace E] (T : ) (hT : 0 T) (B : PC((Set.Icc 0 T), E →L[] E)) (U : (x : P) → Evolution T hT (B x)) (g : C((Set.Icc 0 T), )) (hg : ∀ (t : (Set.Icc 0 T)), 0 < g t) {ι : Type u_3} [Fintype ι] (directions : ιP) (q : ) (f : PC((Set.Icc 0 T), E)) (a₀ : PE) (hf : ContDiff (↑) f) (ha₀ : ContDiff (↑) a₀) (hg₀ : g 0, = 1) (C A D R : ) (hC : 0 C) (x : P) (hU : ∀ (t s : (Set.Icc 0 T)), s t(U x).propagator t s C * g t / g s) (d : ) (hfa : ∀ (n : ), EulerParameterWordGevrey.block directions q f n x D * EulerGevrey.majorant R d n) (haa : ∀ (n : ), EulerParameterWordGevrey.block directions q a₀ n x A * EulerGevrey.majorant R d n) (n : ) :
                                        EulerParameterWordGevrey.block directions q (frozenForcing T hT B U g hg f a₀ x) n x (C * A + C * T * D) * EulerGevrey.majorant R d n

                                        The transformed right side preserves the original fixed-Sobolev external radius.

                                        The fixed-Sobolev inverse estimate needs bounds only at its base point #

                                        In particular a frozen Duhamel equation has identity as its base operator. The actual equation and smoothness hold as functions; every quantitative hypothesis, including invertibility, is needed only at the evaluation point.

                                        theorem EulerParameterWordGevrey.block_inverse_gevrey_at {P : Type u_1} {E : Type u_2} {ι : Type u_3} [NormedAddCommGroup P] [NormedSpace P] [NormedAddCommGroup E] [NormedSpace E] [Fintype ι] (directions : ιP) (q : ) (A : PE →L[] E) (u f : PE) (hA : ContDiff (↑) A) (hu : ContDiff (↑) u) (hf : ContDiff (↑) f) (heq : ∀ (y : P), (A y) (u y) = f y) (x : P) (inverse : E →L[] E) (hleft : ∀ (v : E), inverse ((A x) v) = v) (I B C D M Rc R : ) (_hC : 0 C) (_hD : 0 D) (hM : 1 M) (hMC : sobolevInverseCost I B q * C M) (hMD : sobolevInverseCost I B q * D M) (hRc : 0 Rc) (hR : 2 * M * (Rc + 1) R) (hinv : inverse I) (hbase : baseSize directions q A x B) (hcoeff : ∀ (j : ), coefficientBlock directions q A (j + 1) x C * (Rc ^ (j + 1) * (j + 1).factorial ^ 2)) (d : ) (hforce : ∀ (n : ), block directions q f n x D * EulerGevrey.majorant R d n) (n : ) :
                                        block directions q u n x EulerGevrey.majorant R (d + 1) n
                                        def EulerLinearDuhamel.forwardSobolevAmplitude (ι : Type u_1) [Fintype ι] (q : ) (T C CB Rc : ) :

                                        The once-enlarged coefficient amplitude for the frozen equation at fixed base order.

                                        Equations
                                        Instances For
                                          def EulerLinearDuhamel.forwardSobolevCost (ι : Type u_1) [Fintype ι] (q : ) (T C A D CB Rc : ) :

                                          A fixed polynomial cost for the source's forward Hq external-word estimate.

                                          Equations
                                          • One or more equations did not get rendered due to their size.
                                          Instances For
                                            @[instance_reducible]

                                            Cache the standard NormedAddCommGroup (E →L[ℝ] E) instance to shorten typeclass synthesis.

                                            Equations
                                            Instances For
                                              @[instance_reducible]

                                              Cache the standard NormedSpace ℝ (E →L[ℝ] E) instance to shorten typeclass synthesis.

                                              Equations
                                              Instances For
                                                @[instance_reducible]

                                                Cache the standard NormedAddCommGroup C(Icc (0 : ℝ) T,E) instance to shorten typeclass synthesis.

                                                Equations
                                                Instances For
                                                  @[instance_reducible]

                                                  Cache the standard NormedSpace ℝ C(Icc (0 : ℝ) T,E) instance to shorten typeclass synthesis.

                                                  Equations
                                                  Instances For
                                                    @[instance_reducible]

                                                    Cache the standard NormedAddCommGroup C(Icc (0 : ℝ) T,E →L[ℝ] E) instance to shorten typeclass synthesis.

                                                    Equations
                                                    Instances For
                                                      @[instance_reducible]

                                                      Cache the standard NormedSpace ℝ C(Icc (0 : ℝ) T,E →L[ℝ] E) instance to shorten typeclass synthesis.

                                                      Equations
                                                      Instances For
                                                        @[instance_reducible]

                                                        Cache the standard NormedAddCommGroup (C(Icc (0 : ℝ) T,E) →L[ℝ] C(Icc (0 : ℝ) T,E)) instance to shorten typeclass synthesis.

                                                        Equations
                                                        Instances For
                                                          @[instance_reducible]

                                                          Cache the standard NormedSpace ℝ (C(Icc (0 : ℝ) T,E) →L[ℝ] C(Icc (0 : ℝ) T,E)) instance to shorten typeclass synthesis.

                                                          Equations
                                                          Instances For
                                                            theorem EulerLinearDuhamel.weightedSolution_block_gevrey_at {P : Type u_1} {E : Type u_2} {ι : Type u_3} [NormedAddCommGroup P] [NormedSpace P] [NormedAddCommGroup E] [NormedSpace E] [CompleteSpace E] [Fintype ι] (directions : ιP) (hd : ∀ (i : ι), directions i 1) (q : ) (T : ) (hT : 0 T) (B : PC((Set.Icc 0 T), E →L[] E)) (U : (x : P) → Evolution T hT (B x)) (g : C((Set.Icc 0 T), )) (hg : ∀ (t : (Set.Icc 0 T)), 0 < g t) (f : PC((Set.Icc 0 T), E)) (a₀ : PE) (hB : ContDiff (↑) B) (hf : ContDiff (↑) f) (ha₀ : ContDiff (↑) a₀) (hg₀ : g 0, = 1) (C A D CB Rc R : ) (hC : 0 C) (hA : 0 A) (hD : 0 D) (hCB : 0 CB) (hRc : 0 Rc) (hR : 2 * forwardSobolevCost ι q T C A D CB Rc * (EulerParameterWordGevrey.sobolevCoefficientRadius ι Rc + 1) R) (hBb : ∀ (n : ) (y : P), iteratedFDeriv n B y CB * EulerGevrey.majorant Rc 0 n) (x : P) (hU : ∀ (t s : (Set.Icc 0 T)), s t(U x).propagator t s C * g t / g s) (d : ) (hforce : ∀ (n : ), EulerParameterWordGevrey.block directions q f n x D * EulerGevrey.majorant R d n) (hinitial : ∀ (n : ), EulerParameterWordGevrey.block directions q a₀ n x A * EulerGevrey.majorant R d n) (n : ) :
                                                            EulerParameterWordGevrey.block directions q (fun (y : P) => (U y).weightedSolution g hg (f y) (a₀ y)) n x EulerGevrey.majorant R (d + 1) n

                                                            The actual fixed-Hq block of the forward solution gains just one shift, with an unchanged radius and with H3 used only at the base parameter.