Documentation

LeanPool.NavierStokesAndEuler.Euler.PacketForwardInitializedAllOrderBudget

The literal zero-history initialized packet supplies a complete drift-aware all-order correction budget, from fixed source data and explicit scalar frequency guards. No solution or energy estimate is assumed.

Actual finite-order correction budgets for the zero-history initialized packet. All coefficient and field bounds are supplied by the checked constructions.

Cutoff-independent background, derivative, drift and residual budgets for the actual zero-history correction data.

theorem EulerPacketTerminalDatum.forwardInitializedCorrection_background (M : EulerMeanPacketProvider.Data) {U : Type u_1} [NormedAddCommGroup U] [InnerProductSpace ℝ U] [CompleteSpace U] (D : EulerTransversePacketProvider.Data U) (hTime : M.T = D.T) (δ : ℝ) (hδ : 0 < δ) (ξ : U) (hs : tsupport EulerSpatialCutoffs.innerCutoff ⊆ D.support) (α : ℝ) (L : EulerTransversePacketForward.Budget D (Fin 4) 6) (NB : EulerTransversePacketJoin.NormalBudget D 6 L.R) (W : L.GradeGuards NB 1) (LM : EulerMeanPacketProvider.Budget M 6 L.R) (WM : LM.GradeGuards) (BC : EulerPacketCylinderField.CoefficientBudget (EulerPacketCylinderField.sourceCoefficientData period M D (EulerTransversePacketProvider.InitialData.zero period D) hTime)) (hRc : EulerParameterWordGevrey.sobolevCoefficientRadius (Fin 4) BC.Rc ≤ L.R) (hcost : BC.termCost ≤ L.R) (hδ1 : δ ≤ 1) (hα : 0 < α) (hR : wordRadius (Fin 4) δ ≤ L.R) (WP : L.GradeGuards NB (wordCost (Fin 4) 6 δ * ‖ξ‖)) (S : EulerPacketTimeProfile.Scales ↑(Set.Icc 0 M.T)) (hgrowth : EulerPacketCylinderField.timeProfileChange S.growth hTime = α • L.g) (Cagree : EulerPacketCylinderField.SourceCoefficientAgreement M D) (N : ℕ) (hN : 1 ≤ N) (k : ℝ) (hk : 4 ≤ k) (hbase : EulerPacketCoarseMajorant.tailBase L.R S.H0 BC.termCost N ≤ k ^ (1 / 100)) (s Q : ℕ) (hQ : Q + 6 ≤ s) (ρ : ℝ) (hρ : 0 < ρ) (hsmall : ρ * (4 * L.R) ≤ 1 / 2) (t : ↑(Set.Icc 0 D.T)) :
theorem EulerPacketTerminalDatum.forwardInitializedCorrection_background_derivative (M : EulerMeanPacketProvider.Data) {U : Type u_1} [NormedAddCommGroup U] [InnerProductSpace ℝ U] [CompleteSpace U] (D : EulerTransversePacketProvider.Data U) (hTime : M.T = D.T) (δ : ℝ) (hδ : 0 < δ) (ξ : U) (hs : tsupport EulerSpatialCutoffs.innerCutoff ⊆ D.support) (α : ℝ) (L : EulerTransversePacketForward.Budget D (Fin 4) 6) (NB : EulerTransversePacketJoin.NormalBudget D 6 L.R) (W : L.GradeGuards NB 1) (LM : EulerMeanPacketProvider.Budget M 6 L.R) (WM : LM.GradeGuards) (BC : EulerPacketCylinderField.CoefficientBudget (EulerPacketCylinderField.sourceCoefficientData period M D (EulerTransversePacketProvider.InitialData.zero period D) hTime)) (hRc : EulerParameterWordGevrey.sobolevCoefficientRadius (Fin 4) BC.Rc ≤ L.R) (hcost : BC.termCost ≤ L.R) (hδ1 : δ ≤ 1) (hα : 0 < α) (hR : wordRadius (Fin 4) δ ≤ L.R) (WP : L.GradeGuards NB (wordCost (Fin 4) 6 δ * ‖ξ‖)) (S : EulerPacketTimeProfile.Scales ↑(Set.Icc 0 M.T)) (hgrowth : EulerPacketCylinderField.timeProfileChange S.growth hTime = α • L.g) (Cagree : EulerPacketCylinderField.SourceCoefficientAgreement M D) (N : ℕ) (hN : 1 ≤ N) (k : ℝ) (hk : 4 ≤ k) (hbase : EulerPacketCoarseMajorant.tailBase L.R S.H0 BC.termCost N ≤ k ^ (1 / 100)) (s Q : ℕ) (hQ : Q + 6 ≤ s) (ρ : ℝ) (hρ : 0 < ρ) (hsmall : ρ * (4 * L.R) ≤ 1 / 2) (t : ↑(Set.Icc 0 D.T)) :
theorem EulerPacketTerminalDatum.forwardInitializedCorrection_drift (M : EulerMeanPacketProvider.Data) {U : Type u_1} [NormedAddCommGroup U] [InnerProductSpace ℝ U] [CompleteSpace U] (D : EulerTransversePacketProvider.Data U) (hTime : M.T = D.T) (δ : ℝ) (hδ : 0 < δ) (ξ : U) (hs : tsupport EulerSpatialCutoffs.innerCutoff ⊆ D.support) (α : ℝ) (L : EulerTransversePacketForward.Budget D (Fin 4) 6) (NB : EulerTransversePacketJoin.NormalBudget D 6 L.R) (W : L.GradeGuards NB 1) (LM : EulerMeanPacketProvider.Budget M 6 L.R) (WM : LM.GradeGuards) (BC : EulerPacketCylinderField.CoefficientBudget (EulerPacketCylinderField.sourceCoefficientData period M D (EulerTransversePacketProvider.InitialData.zero period D) hTime)) (hRc : EulerParameterWordGevrey.sobolevCoefficientRadius (Fin 4) BC.Rc ≤ L.R) (hcost : BC.termCost ≤ L.R) (hδ1 : δ ≤ 1) (hα : 0 < α) (hR : wordRadius (Fin 4) δ ≤ L.R) (WP : L.GradeGuards NB (wordCost (Fin 4) 6 δ * ‖ξ‖)) (S : EulerPacketTimeProfile.Scales ↑(Set.Icc 0 M.T)) (hgrowth : EulerPacketCylinderField.timeProfileChange S.growth hTime = α • L.g) (Cagree : EulerPacketCylinderField.SourceCoefficientAgreement M D) (N : ℕ) (hN : 1 ≤ N) (k : ℝ) (hk : 4 ≤ k) (hbase : EulerPacketCoarseMajorant.tailBase L.R S.H0 BC.termCost N ≤ k ^ (1 / 100)) (s Q : ℕ) (hQ : Q + 6 ≤ s) (ρ : ℝ) (hρ : 0 < ρ) (hsmall : ρ * (4 * L.R) ≤ 1 / 2) (t : ↑(Set.Icc 0 D.T)) :
theorem EulerPacketTerminalDatum.forwardInitializedCorrection_residual (M : EulerMeanPacketProvider.Data) {U : Type u_1} [NormedAddCommGroup U] [InnerProductSpace ℝ U] [CompleteSpace U] (D : EulerTransversePacketProvider.Data U) (hTime : M.T = D.T) (δ : ℝ) (hδ : 0 < δ) (ξ : U) (hs : tsupport EulerSpatialCutoffs.innerCutoff ⊆ D.support) (α : ℝ) (L : EulerTransversePacketForward.Budget D (Fin 4) 6) (NB : EulerTransversePacketJoin.NormalBudget D 6 L.R) (W : L.GradeGuards NB 1) (LM : EulerMeanPacketProvider.Budget M 6 L.R) (WM : LM.GradeGuards) (BC : EulerPacketCylinderField.CoefficientBudget (EulerPacketCylinderField.sourceCoefficientData period M D (EulerTransversePacketProvider.InitialData.zero period D) hTime)) (hRc : EulerParameterWordGevrey.sobolevCoefficientRadius (Fin 4) BC.Rc ≤ L.R) (hcost : BC.termCost ≤ L.R) (hδ1 : δ ≤ 1) (hα : 0 < α) (hR : wordRadius (Fin 4) δ ≤ L.R) (WP : L.GradeGuards NB (wordCost (Fin 4) 6 δ * ‖ξ‖)) (S : EulerPacketTimeProfile.Scales ↑(Set.Icc 0 M.T)) (hgrowth : EulerPacketCylinderField.timeProfileChange S.growth hTime = α • L.g) (Cagree : EulerPacketCylinderField.SourceCoefficientAgreement M D) (N : ℕ) (hN : 1 ≤ N) (k : ℝ) (hk : 4 ≤ k) (hbase : EulerPacketCoarseMajorant.tailBase L.R S.H0 BC.termCost N ≤ k ^ (1 / 100)) (X : ℝ) (hcoef : BC.multiplierCost ≤ k ^ (1 / 100)) (hX : 6 ≤ X) (hNX : X - 1 ≤ ↑N) (s Q : ℕ) (hQ : Q + 6 ≤ s) (ρ : ℝ) (hρ : 0 < ρ) (hsmall : ρ * (4 * L.R) ≤ 1 / 2) (t : ↑(Set.Icc 0 D.T)) :
EulerSobolevGevreyOperators.weightedNorm period 6 Q ρ (((forwardInitializedCorrectionData M D hTime δ hδ ξ hs α Cagree N hN k hk).residual.realization s) t) ≤ 2 * Real.exp (-(7 / 10) * X * Real.log k)

Forward initialized metric budget, constructed using sourceMetricBudgetOfFields.

Equations
  • One or more equations did not get rendered due to their size.
Instances For
    noncomputable def EulerPacketTerminalDatum.forwardInitializedSpatialBudget (M : EulerMeanPacketProvider.Data) {U : Type u_1} [NormedAddCommGroup U] [InnerProductSpace ℝ U] [CompleteSpace U] (D : EulerTransversePacketProvider.Data U) (hTime : M.T = D.T) (δ : ℝ) (hδ : 0 < δ) (ξ : U) (hs : tsupport EulerSpatialCutoffs.innerCutoff ⊆ D.support) (α : ℝ) (Cagree : EulerPacketCylinderField.SourceCoefficientAgreement M D) (N : ℕ) (hN : 1 ≤ N) (k : ℝ) (hk : 4 ≤ k) (L : EulerTransversePacketForward.Budget D (Fin 4) 6) (NB : EulerTransversePacketJoin.NormalBudget D 6 L.R) (W : L.GradeGuards NB 1) (LM : EulerMeanPacketProvider.Budget M 6 L.R) (WM : LM.GradeGuards) (BC : EulerPacketCylinderField.CoefficientBudget (EulerPacketCylinderField.sourceCoefficientData period M D (EulerTransversePacketProvider.InitialData.zero period D) hTime)) (hRc : EulerParameterWordGevrey.sobolevCoefficientRadius (Fin 4) BC.Rc ≤ L.R) (hcost : BC.termCost ≤ L.R) (hδ1 : δ ≤ 1) (hα : 0 < α) (hR : wordRadius (Fin 4) δ ≤ L.R) (WP : L.GradeGuards NB (wordCost (Fin 4) 6 δ * ‖ξ‖)) (S : EulerPacketTimeProfile.Scales ↑(Set.Icc 0 M.T)) (hgrowth : EulerPacketCylinderField.timeProfileChange S.growth hTime = α • L.g) (hbase : EulerPacketCoarseMajorant.tailBase L.R S.H0 BC.termCost N ≤ k ^ (1 / 100)) (X : ℝ) (hcoef : BC.multiplierCost ≤ k ^ (1 / 100)) (hX : 6 ≤ X) (hNX : X - 1 ≤ ↑N) (Kc : EulerPacketCorrectionCoefficients.CorrectionCoefficientBudget D period) (ρ : C(↑(Set.Icc 0 D.T), ℝ)) (hρ : ∀ (t : ↑(Set.Icc 0 D.T)), 0 < ρ t) (hpacket : ∀ (t : ↑(Set.Icc 0 D.T)), ρ t * (4 * L.R) ≤ 1 / 2) (hpressure : ∀ (t : ↑(Set.Icc 0 D.T)), 4 * Kc.M * (ρ t * Kc.Rc) ≤ 1) (q : ℕ) (hq : 6 ≤ q) :

    All four field estimates and all coefficient estimates are actual properties of the initialized source data at this finite Sobolev order.

    Equations
    • One or more equations did not get rendered due to their size.
    Instances For
      noncomputable def EulerPacketTerminalDatum.forwardInitializedDriftBudget (M : EulerMeanPacketProvider.Data) {U : Type u_1} [NormedAddCommGroup U] [InnerProductSpace ℝ U] [CompleteSpace U] (D : EulerTransversePacketProvider.Data U) (hTime : M.T = D.T) (δ : ℝ) (hδ : 0 < δ) (ξ : U) (hs : tsupport EulerSpatialCutoffs.innerCutoff ⊆ D.support) (α : ℝ) (Cagree : EulerPacketCylinderField.SourceCoefficientAgreement M D) (N : ℕ) (hN : 1 ≤ N) (k : ℝ) (hk : 4 ≤ k) (L : EulerTransversePacketForward.Budget D (Fin 4) 6) (NB : EulerTransversePacketJoin.NormalBudget D 6 L.R) (W : L.GradeGuards NB 1) (LM : EulerMeanPacketProvider.Budget M 6 L.R) (WM : LM.GradeGuards) (BC : EulerPacketCylinderField.CoefficientBudget (EulerPacketCylinderField.sourceCoefficientData period M D (EulerTransversePacketProvider.InitialData.zero period D) hTime)) (hRc : EulerParameterWordGevrey.sobolevCoefficientRadius (Fin 4) BC.Rc ≤ L.R) (hcost : BC.termCost ≤ L.R) (hδ1 : δ ≤ 1) (hα : 0 < α) (hR : wordRadius (Fin 4) δ ≤ L.R) (WP : L.GradeGuards NB (wordCost (Fin 4) 6 δ * ‖ξ‖)) (S : EulerPacketTimeProfile.Scales ↑(Set.Icc 0 M.T)) (hgrowth : EulerPacketCylinderField.timeProfileChange S.growth hTime = α • L.g) (hbase : EulerPacketCoarseMajorant.tailBase L.R S.H0 BC.termCost N ≤ k ^ (1 / 100)) (X : ℝ) (hcoef : BC.multiplierCost ≤ k ^ (1 / 100)) (hX : 6 ≤ X) (hNX : X - 1 ≤ ↑N) (Kc : EulerPacketCorrectionCoefficients.CorrectionCoefficientBudget D period) (ρ : C(↑(Set.Icc 0 D.T), ℝ)) (hρ : ∀ (t : ↑(Set.Icc 0 D.T)), 0 < ρ t) (hpacket : ∀ (t : ↑(Set.Icc 0 D.T)), ρ t * (4 * L.R) ≤ 1 / 2) (hpressure : ∀ (t : ↑(Set.Icc 0 D.T)), 4 * Kc.M * (ρ t * Kc.Rc) ≤ 1) (q : ℕ) (hq : 6 ≤ q) :

      The small drift envelope is kept separate from the full background.

      Equations
      • One or more equations did not get rendered due to their size.
      Instances For
        theorem EulerPacketTerminalDatum.forwardInitializedDriftBudget_growth (M : EulerMeanPacketProvider.Data) {U : Type u_1} [NormedAddCommGroup U] [InnerProductSpace ℝ U] [CompleteSpace U] (D : EulerTransversePacketProvider.Data U) (hTime : M.T = D.T) (δ : ℝ) (hδ : 0 < δ) (ξ : U) (hs : tsupport EulerSpatialCutoffs.innerCutoff ⊆ D.support) (α : ℝ) (Cagree : EulerPacketCylinderField.SourceCoefficientAgreement M D) (N : ℕ) (hN : 1 ≤ N) (k : ℝ) (hk : 4 ≤ k) (L : EulerTransversePacketForward.Budget D (Fin 4) 6) (NB : EulerTransversePacketJoin.NormalBudget D 6 L.R) (W : L.GradeGuards NB 1) (LM : EulerMeanPacketProvider.Budget M 6 L.R) (WM : LM.GradeGuards) (BC : EulerPacketCylinderField.CoefficientBudget (EulerPacketCylinderField.sourceCoefficientData period M D (EulerTransversePacketProvider.InitialData.zero period D) hTime)) (hRc : EulerParameterWordGevrey.sobolevCoefficientRadius (Fin 4) BC.Rc ≤ L.R) (hcost : BC.termCost ≤ L.R) (hδ1 : δ ≤ 1) (hα : 0 < α) (hR : wordRadius (Fin 4) δ ≤ L.R) (WP : L.GradeGuards NB (wordCost (Fin 4) 6 δ * ‖ξ‖)) (S : EulerPacketTimeProfile.Scales ↑(Set.Icc 0 M.T)) (hgrowth : EulerPacketCylinderField.timeProfileChange S.growth hTime = α • L.g) (hbase : EulerPacketCoarseMajorant.tailBase L.R S.H0 BC.termCost N ≤ k ^ (1 / 100)) (X : ℝ) (hcoef : BC.multiplierCost ≤ k ^ (1 / 100)) (hX : 6 ≤ X) (hNX : X - 1 ≤ ↑N) (Kc : EulerPacketCorrectionCoefficients.CorrectionCoefficientBudget D period) (ρ : C(↑(Set.Icc 0 D.T), ℝ)) (hρ : ∀ (t : ↑(Set.Icc 0 D.T)), 0 < ρ t) (hpacket : ∀ (t : ↑(Set.Icc 0 D.T)), ρ t * (4 * L.R) ≤ 1 / 2) (hpressure : ∀ (t : ↑(Set.Icc 0 D.T)), 4 * Kc.M * (ρ t * Kc.Rc) ≤ 1) (q : ℕ) (hq : 6 ≤ q) :
        EulerCorrectionEnergyMajorants.combinedConstant period (forwardInitializedDriftBudget M D hTime δ hδ ξ hs α Cagree N hN k hk L NB W LM WM BC hRc hcost hδ1 hα hR WP S hgrowth hbase X hcoef hX hNX Kc ρ hρ hpacket hpressure q hq).full (EulerAllOrderCorrectionData.Data.metricBudget period (forwardInitializedCorrectionData M D hTime δ hδ ξ hs α Cagree N hN k hk) ⋯ (forwardInitializedMetricBudget M D hTime δ hδ ξ hs α Cagree N hN k hk 0) (q + 1)) = EulerPacketCorrectionCoefficients.growthCoefficient D period Kc (2 * EulerPacketCorrectionConstants.velocity L.R S.H0 BC.multiplierCost) (12 * EulerPacketCorrectionConstants.velocity L.R S.H0 BC.multiplierCost * (4 * L.R))

        The zero-history initialized approximation satisfies the lifted divergence constraint for an actual volume-preserving source deformation.

        theorem EulerPacketTerminalDatum.forwardInitializedCorrectionData_divergence (M : EulerMeanPacketProvider.Data) {U : Type u_1} [NormedAddCommGroup U] [InnerProductSpace ℝ U] [CompleteSpace U] (D : EulerTransversePacketProvider.Data U) (hTime : M.T = D.T) (δ : ℝ) (hδ : 0 < δ) (ξ : U) (hs : tsupport EulerSpatialCutoffs.innerCutoff ⊆ D.support) (α : ℝ) (Cagree : EulerPacketCylinderField.SourceCoefficientAgreement M D) (N : ℕ) (hN : 1 ≤ N) (k : ℝ) (hk : 4 ≤ k) (Ξ : ↑(Set.Icc 0 D.T) → EulerSmoothLimit.Space → EulerSmoothLimit.Space) (hΞ : ∀ (t : ↑(Set.Icc 0 D.T)), ContDiff ℝ (↑⊤) (Ξ t)) (hF : ∀ (t : ↑(Set.Icc 0 D.T)) (x : EulerSmoothLimit.Space), fderiv ℝ (Ξ t) x = (D.F.field t) x) (hdet : ∀ (t : ↑(Set.Icc 0 D.T)) (x : EulerSmoothLimit.Space), Matrix.det (EulerPacketPiola.operatorMatrix ((D.F.field t) x)) = 1) (t : ↑(Set.Icc 0 D.T)) :
        noncomputable def EulerPacketTerminalDatum.forwardInitializedAllOrderBudget (M : EulerMeanPacketProvider.Data) {U : Type u_1} [NormedAddCommGroup U] [InnerProductSpace ℝ U] [CompleteSpace U] (D : EulerTransversePacketProvider.Data U) (hTime : M.T = D.T) (δ : ℝ) (hδ : 0 < δ) (ξ : U) (hs : tsupport EulerSpatialCutoffs.innerCutoff ⊆ D.support) (α : ℝ) (Cagree : EulerPacketCylinderField.SourceCoefficientAgreement M D) (L : EulerTransversePacketForward.Budget D (Fin 4) 6) (NB : EulerTransversePacketJoin.NormalBudget D 6 L.R) (W : L.GradeGuards NB 1) (LM : EulerMeanPacketProvider.Budget M 6 L.R) (WM : LM.GradeGuards) (BC : EulerPacketCylinderField.CoefficientBudget (EulerPacketCylinderField.sourceCoefficientData period M D (EulerTransversePacketProvider.InitialData.zero period D) hTime)) (hRc : EulerParameterWordGevrey.sobolevCoefficientRadius (Fin 4) BC.Rc ≤ L.R) (hcost : BC.termCost ≤ L.R) (hδ1 : δ ≤ 1) (hα : 0 < α) (hR : wordRadius (Fin 4) δ ≤ L.R) (WP : L.GradeGuards NB (wordCost (Fin 4) 6 δ * ‖ξ‖)) (S : EulerPacketTimeProfile.Scales ↑(Set.Icc 0 M.T)) (hgrowth : EulerPacketCylinderField.timeProfileChange S.growth hTime = α • L.g) (Kc : EulerPacketCorrectionCoefficients.CorrectionCoefficientBudget D period) (k : ℝ) (hk : 4 ≤ k) (hX : 64 ≤ EulerPacketSourceFrequency.expansion k) (hlog : 1 ≤ Real.log k) (htail : EulerPacketCoarseMajorant.tailPolynomialConstant L.R S.H0 BC.termCost ≤ EulerPacketSourceFrequency.smallPower k) (hcoefficient : BC.multiplierCost ≤ EulerPacketSourceFrequency.smallPower k) (hgrowthCost : 12 * EulerPacketCorrectionCoefficients.growthCoefficient D period Kc (2 * EulerPacketCorrectionConstants.velocity L.R S.H0 BC.multiplierCost) (12 * EulerPacketCorrectionConstants.velocity L.R S.H0 BC.multiplierCost * (4 * L.R)) * D.T ≤ EulerPacketSourceFrequency.smallPower k) (hdriftCost : 8 * EulerPacketCorrectionCoefficients.growthCoefficient D period Kc (2 * EulerPacketCorrectionConstants.velocity L.R S.H0 BC.multiplierCost) (12 * EulerPacketCorrectionConstants.velocity L.R S.H0 BC.multiplierCost * (4 * L.R)) * D.T * EulerPacketCorrectionConstants.drift L.R S.H0 BC.multiplierCost / EulerPacketCorrectionScalar.initialRadius L.R Kc.M Kc.Rc ≤ EulerPacketSourceFrequency.smallPower k) (herrorCost : 8 * EulerPacketCorrectionCoefficients.growthCoefficient D period Kc (2 * EulerPacketCorrectionConstants.velocity L.R S.H0 BC.multiplierCost) (12 * EulerPacketCorrectionConstants.velocity L.R S.H0 BC.multiplierCost * (4 * L.R)) * D.T / EulerPacketCorrectionScalar.initialRadius L.R Kc.M Kc.Rc ≤ EulerPacketSourceFrequency.smallPower k) (Ξ : ↑(Set.Icc 0 D.T) → EulerSmoothLimit.Space → EulerSmoothLimit.Space) (hΞ : ∀ (t : ↑(Set.Icc 0 D.T)), ContDiff ℝ (↑⊤) (Ξ t)) (hF : ∀ (t : ↑(Set.Icc 0 D.T)) (x : EulerSmoothLimit.Space), fderiv ℝ (Ξ t) x = (D.F.field t) x) (hdet : ∀ (t : ↑(Set.Icc 0 D.T)) (x : EulerSmoothLimit.Space), Matrix.det (EulerPacketPiola.operatorMatrix ((D.F.field t) x)) = 1) :

        Explicit scalar guards suffice because every analytic input to the all-order correction theorem is supplied by the constructed packet.

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