Documentation

LeanPool.NavierStokesAndEuler.Euler.AllOrderDriftCorrection

One all-order correction from the actual small-drift construction #

The finite solutions used here are Budget.solution from AllOrderDriftFinite. Their existence is proved by the drift-aware finite-Sobolev solver. The generic assembly proves compatibility from their literal equations and uniqueness; it does not require the coarse full-velocity shrinking-radius assumption.

The resulting single field retains the residual and target-error estimates at every external cutoff and satisfies the actual equation in every finite Sobolev order. Smoothness asserted here is spatial smoothness, with a jointly continuous representative and its genuine first time derivative.

Genuine pointwise time differentiation of the generically assembled correction.

noncomputable def EulerCorrectionAssembly.FiniteFamily.pointRawSource (period : ) [Fact (0 < period)] {T : } {hT : 0 < T} {A : EulerAllOrderCorrectionData.Data period T} (F : FiniteFamily period hT A) (t : (Set.Icc 0 T)) (x : EulerLiftedGradientSpace.LiftDomain period) :

The canonical pointwise nonlinear raw source of the actual common correction.

Equations
  • One or more equations did not get rendered due to their size.
Instances For
    noncomputable def EulerCorrectionAssembly.FiniteFamily.pointTimeDerivative (period : ) [Fact (0 < period)] {T : } {hT : 0 < T} {A : EulerAllOrderCorrectionData.Data period T} (F : FiniteFamily period hT A) (t : (Set.Icc 0 T)) (x : EulerLiftedGradientSpace.LiftDomain period) :

    The canonical actual time derivative, defined by bounded evaluation of the genuine continuous Sobolev source.

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

      The actual pointwise time derivative is jointly continuous in time and space.

      theorem EulerCorrectionAssembly.FiniteFamily.pointTimeDerivative_eq_pressure (period : ) [Fact (0 < period)] {T : } {hT : 0 < T} {A : EulerAllOrderCorrectionData.Data period T} (F : FiniteFamily period hT A) (t : (Set.Icc 0 T)) (x : EulerLiftedGradientSpace.LiftDomain period) :
      pointTimeDerivative period F t x = -pointRawSource period F t x - ((A.metric.coefficient t).coefficient x) (pointPressure period F t x)

      The pointwise time derivative is the literal raw-source and signed-pressure expression.

      theorem EulerCorrectionAssembly.FiniteFamily.pointField_hasDerivAt (period : ) [Fact (0 < period)] {T : } {hT : 0 < T} {A : EulerAllOrderCorrectionData.Data period T} (F : FiniteFamily period hT A) (x : EulerLiftedGradientSpace.LiftDomain period) (t : ) (ht : t Set.Ioo 0 T) :
      HasDerivAt (fun (r : ) => pointField period F (Set.projIcc 0 T r) x) (pointTimeDerivative period F t, x) t

      The canonical common field has its genuine pointwise first time derivative at every interior time.

      theorem EulerCorrectionAssembly.FiniteFamily.pointField_hasDerivAt_pressure (period : ) [Fact (0 < period)] {T : } {hT : 0 < T} {A : EulerAllOrderCorrectionData.Data period T} (F : FiniteFamily period hT A) (x : EulerLiftedGradientSpace.LiftDomain period) (t : ) (ht : t Set.Ioo 0 T) :
      HasDerivAt (fun (r : ) => pointField period F (Set.projIcc 0 T r) x) (-pointRawSource period F t, x - ((A.metric.coefficient t, ).coefficient x) (pointPressure period F t, x)) t

      The actual pointwise correction equation uses the reconstructed signed pressure and literal matrix multiplication.

      theorem EulerAllOrderDriftCorrection.Budget.solution_compatible (period : ) [Fact (0 < period)] {T : } {hT : 0 < T} {A : EulerAllOrderCorrectionData.Data period T} (B : Budget period hT A) (q : ) (hq : 6 q) :

      The drift-aware finite solutions coincide under the actual Sobolev restriction.

      noncomputable def EulerAllOrderDriftCorrection.Budget.commonPath (period : ) [Fact (0 < period)] {T : } {hT : 0 < T} {A : EulerAllOrderCorrectionData.Data period T} (B : Budget period hT A) :

      The common continuous L² correction constructed by the small-drift solver.

      Equations
      Instances For
        theorem EulerAllOrderDriftCorrection.Budget.solution_value_common (period : ) [Fact (0 < period)] {T : } {hT : 0 < T} {A : EulerAllOrderCorrectionData.Data period T} (B : Budget period hT A) (q : ) (hq : 6 q) (t : (Set.Icc 0 T)) :
        EulerCylinderSobolevSpace.value period ((solution period B q hq) t) = (commonPath period B) t

        Every actual drift-aware finite solution represents this same common field.

        theorem EulerAllOrderDriftCorrection.Budget.commonPath_initial (period : ) [Fact (0 < period)] {T : } {hT : 0 < T} {A : EulerAllOrderCorrectionData.Data period T} (B : Budget period hT A) :
        (commonPath period B) 0, = 0

        The common correction has zero initial trace.

        theorem EulerAllOrderDriftCorrection.Budget.commonPath_divergence (period : ) [Fact (0 < period)] {T : } {hT : 0 < T} {A : EulerAllOrderCorrectionData.Data period T} (B : Budget period hT A) (t : (Set.Icc 0 T)) :

        The common correction satisfies the genuine closed lifted divergence constraint.

        The common L² path satisfies the actual projected correction equation.

        noncomputable def EulerAllOrderDriftCorrection.Budget.commonJet (period : ) [Fact (0 < period)] {T : } {hT : 0 < T} {A : EulerAllOrderCorrectionData.Data period T} (B : Budget period hT A) (n : ) (t : (Set.Icc 0 T)) :

        An actual strong spatial jet of every order for the constructed common field.

        Equations
        • One or more equations did not get rendered due to their size.
        Instances For
          noncomputable def EulerAllOrderDriftCorrection.Budget.fieldTower (period : ) [Fact (0 < period)] {T : } {hT : 0 < T} {A : EulerAllOrderCorrectionData.Data period T} (B : Budget period hT A) :

          The common correction, with genuine continuous Sobolev realizations at all orders.

          Equations
          • One or more equations did not get rendered due to their size.
          Instances For
            theorem EulerAllOrderDriftCorrection.Budget.fieldTower_field (period : ) [Fact (0 < period)] {T : } {hT : 0 < T} {A : EulerAllOrderCorrectionData.Data period T} (B : Budget period hT A) :
            (fieldTower period B).field = commonPath period B

            The tower's underlying field is exactly the constructed common correction.

            theorem EulerAllOrderDriftCorrection.Budget.fieldTower_initial (period : ) [Fact (0 < period)] {T : } {hT : 0 < T} {A : EulerAllOrderCorrectionData.Data period T} (B : Budget period hT A) (q : ) :
            ((fieldTower period B).realization q) 0, = 0

            Every Sobolev realization of the common correction has zero initial data.

            theorem EulerAllOrderDriftCorrection.Budget.solution_eq_realization (period : ) [Fact (0 < period)] {T : } {hT : 0 < T} {A : EulerAllOrderCorrectionData.Data period T} (B : Budget period hT A) (q : ) (hq : 6 q) :
            solution period B q hq = (fieldTower period B).realization (q + 1)

            A finite drift-aware solution equals the common tower's realization at its order.

            theorem EulerAllOrderDriftCorrection.Budget.fieldTower_energy (period : ) [Fact (0 < period)] {T : } {hT : 0 < T} {A : EulerAllOrderCorrectionData.Data period T} (B : Budget period hT A) (P : ) (t : (Set.Icc 0 T)) :

            Both quantitative finite-solver bounds hold for the actual common realization. The radius and residual envelope are the drift-aware input budgets.

            theorem EulerAllOrderDriftCorrection.Budget.solution_sobolev_hasDerivAt (period : ) [Fact (0 < period)] {T : } {hT : 0 < T} {A : EulerAllOrderCorrectionData.Data period T} (B : Budget period hT A) (q : ) (hq : 6 q) (t : ) (ht : t Set.Ioo 0 T) :

            The constructed finite solution has its genuine time derivative in Hq.

            theorem EulerAllOrderDriftCorrection.Budget.fieldTower_hasDerivAt (period : ) [Fact (0 < period)] {T : } {hT : 0 < T} {A : EulerAllOrderCorrectionData.Data period T} (B : Budget period hT A) (q : ) (hq : 6 q) (t : ) (ht : t Set.Ioo 0 T) :

            Every finite Sobolev realization of the common correction satisfies the actual projected nonlinear evolution, not merely an equation for an unrelated finite solve.

            The common realization's actual derivative is the literal raw source plus signed pressure.

            noncomputable def EulerAllOrderDriftCorrection.Budget.pointField (period : ) [Fact (0 < period)] {T : } {hT : 0 < T} {A : EulerAllOrderCorrectionData.Data period T} (B : Budget period hT A) (t : (Set.Icc 0 T)) (x : EulerLiftedGradientSpace.LiftDomain period) :

            Canonical bounded evaluation fixes a pointwise representative of the correction.

            Equations
            Instances For
              theorem EulerAllOrderDriftCorrection.Budget.pointField_ae (period : ) [Fact (0 < period)] {T : } {hT : 0 < T} {A : EulerAllOrderCorrectionData.Data period T} (B : Budget period hT A) (t : (Set.Icc 0 T)) :

              The pointwise correction is an actual representative of the constructed L² field.

              theorem EulerAllOrderDriftCorrection.Budget.pointField_initial (period : ) [Fact (0 < period)] {T : } {hT : 0 < T} {A : EulerAllOrderCorrectionData.Data period T} (B : Budget period hT A) (x : EulerLiftedGradientSpace.LiftDomain period) :
              pointField period B 0, x = 0

              The canonical correction vanishes pointwise at the initial time.

              theorem EulerAllOrderDriftCorrection.Budget.pointField_joint_continuous (period : ) [Fact (0 < period)] {T : } {hT : 0 < T} {A : EulerAllOrderCorrectionData.Data period T} (B : Budget period hT A) :

              The canonical correction is jointly continuous in time and the cylinder point.

              theorem EulerAllOrderDriftCorrection.Budget.pointField_smooth (period : ) [Fact (0 < period)] {T : } {hT : 0 < T} {A : EulerAllOrderCorrectionData.Data period T} (B : Budget period hT A) (t : (Set.Icc 0 T)) (x : EulerLiftedGradientSpace.LiftDomain period) :

              The one common correction has a spatially smooth representative at every time.

              theorem EulerAllOrderDriftCorrection.Budget.pointField_divergence (period : ) [Fact (0 < period)] {T : } {hT : 0 < T} {A : EulerAllOrderCorrectionData.Data period T} (B : Budget period hT A) (t : (Set.Icc 0 T)) (x : EulerLiftedGradientSpace.LiftDomain period) :

              Its lifted divergence vanishes pointwise.

              theorem EulerAllOrderDriftCorrection.Budget.commonPath_odd (period : ) [Fact (0 < period)] {T : } {hT : 0 < T} {A : EulerAllOrderCorrectionData.Data period T} (B : Budget period hT A) (P : EulerCorrectionAssembly.ParityData period A) (t : (Set.Icc 0 T)) :
              -(EulerCylinderReflection.reflection period) ((commonPath period B) t) = (commonPath period B) t

              Odd input data give an odd common correction by the proved PDE uniqueness.

              theorem EulerAllOrderDriftCorrection.Budget.fieldTower_odd (period : ) [Fact (0 < period)] {T : } {hT : 0 < T} {A : EulerAllOrderCorrectionData.Data period T} (B : Budget period hT A) (P : EulerCorrectionAssembly.ParityData period A) (q : ) (t : (Set.Icc 0 T)) :

              The actual common Sobolev realizations retain the prescribed odd parity.

              theorem EulerAllOrderDriftCorrection.Budget.pointField_odd (period : ) [Fact (0 < period)] {T : } {hT : 0 < T} {A : EulerAllOrderCorrectionData.Data period T} (B : Budget period hT A) (P : EulerCorrectionAssembly.ParityData period A) (t : (Set.Icc 0 T)) (x : EulerLiftedGradientSpace.LiftDomain period) :
              pointField period B t (-x) = -pointField period B t x

              The canonical spatially smooth correction is pointwise odd for odd input data.

              noncomputable def EulerAllOrderDriftCorrection.Budget.pointTimeDerivative (period : ) [Fact (0 < period)] {T : } {hT : 0 < T} {A : EulerAllOrderCorrectionData.Data period T} (B : Budget period hT A) (t : (Set.Icc 0 T)) (x : EulerLiftedGradientSpace.LiftDomain period) :

              The actual pointwise time derivative, obtained from the constructed Sobolev source.

              Equations
              Instances For

                The correction's actual first time derivative is jointly continuous.

                The actual time derivative equals the raw source and the signed pressure term.

                theorem EulerAllOrderDriftCorrection.Budget.pointField_hasDerivAt (period : ) [Fact (0 < period)] {T : } {hT : 0 < T} {A : EulerAllOrderCorrectionData.Data period T} (B : Budget period hT A) (x : EulerLiftedGradientSpace.LiftDomain period) (t : ) (ht : t Set.Ioo 0 T) :
                HasDerivAt (fun (r : ) => pointField period B (Set.projIcc 0 T r) x) (pointTimeDerivative period B t, x) t

                The canonical common field has its genuine first time derivative at interior times.

                theorem EulerAllOrderDriftCorrection.Budget.pointField_hasDerivAt_pressure (period : ) [Fact (0 < period)] {T : } {hT : 0 < T} {A : EulerAllOrderCorrectionData.Data period T} (B : Budget period hT A) (x : EulerLiftedGradientSpace.LiftDomain period) (t : ) (ht : t Set.Ioo 0 T) :

                The pointwise correction equation contains the actual raw source and signed pressure.

                theorem EulerAllOrderDriftCorrection.exists_smooth_lifted_correction (period : ) [Fact (0 < period)] {T : } {hT : 0 < T} {A : EulerAllOrderCorrectionData.Data period T} (B : Budget period hT A) :
                ∃ (E : EulerAllOrderCorrectionData.FieldTower period T) (g : (Set.Icc 0 T)EulerLiftedGradientSpace.LiftDomain periodEulerLiftedGradientSpace.Vector3) (v : (Set.Icc 0 T)EulerLiftedGradientSpace.LiftDomain periodEulerLiftedGradientSpace.Vector3), E.field 0, = 0 (∀ (t : (Set.Icc 0 T)), E.field t EulerLiftedGradientSpace.divergenceFreeSpace period A.κ A.direction) Continuous (Function.uncurry g) (∀ (t : (Set.Icc 0 T)) (x : EulerLiftedGradientSpace.LiftDomain period), ContDiff (↑) (EulerMetricTransport.localFieldLift period (g t) x)) (∀ (t : (Set.Icc 0 T)), (E.field t) =ᵐ[EulerLiftedGradientSpace.liftMeasure period] g t) (∀ (t : (Set.Icc 0 T)) (x : EulerLiftedGradientSpace.LiftDomain period), i : Fin 3, (EulerTransportDerivatives.fieldDerivative period (EulerMetricTransport.coordinateDirection A.κ A.direction i) (g t) x).ofLp i = 0) Continuous (Function.uncurry v) (∀ (x : EulerLiftedGradientSpace.LiftDomain period) (t : ) (ht : t Set.Ioo 0 T), HasDerivAt (fun (r : ) => g (Set.projIcc 0 T r) x) (v t, x) t) (∀ (P : ) (t : (Set.Icc 0 T)), EulerGevreyMetricEstimate.energyNorm period P (B.radius t) ((EulerCorrectionEnergyData.MetricBudget.operatorPath period B.metric) t) ((E.realization (P + 6 + 1)) t) 2 * (B.spatial (P + 6) ).full.residual * Real.exp (3 * B.growthCoefficient * t) EulerGevreyMetricEstimate.energyNorm period P (B.radius t) ((EulerCorrectionEnergyData.MetricBudget.operatorPath period B.metric) t) ((E.realization (P + 6 + 1)) t) B.delta / 2) ∀ (q : ) (hq : 6 q) (t : ) (ht : t Set.Ioo 0 T), HasDerivAt (EulerVolterraConvolution.extendPath T (E.realization q)) ((EulerCorrectionOperators.CorrectionData.coefficients period (EulerAllOrderCorrectionData.Data.atOrder period A q) hq).apply t, ((E.realization (q + 1)) t, )) t

                Genuine coherent input bounds, with radius loss determined only by the actual transport drift, construct one smooth spatial correction with all-cutoff energy bounds and its actual finite-Sobolev evolution. Finite existence, compatibility, energy estimates and the correction equation are conclusions here.