Documentation

LeanPool.NavierStokesAndEuler.Euler.PacketExactPressureError

The canonical scalar pressure of the actual exact packet has the finite pressure's Hessian plus the Hessian of its actual correction. The normalization of the scalar potential does not affect this identity.

noncomputable def EulerPacketTerminalDatum.initializedExactPhysicalPressure (M : EulerMeanPacketProvider.Data) {U : Type u_1} [NormedAddCommGroup U] [InnerProductSpace U] [CompleteSpace U] (D : EulerTransversePacketProvider.Data U) (hTime : M.T = D.T) (τ : ) ( : 0 < τ) (hτT : τ < D.T) (B : EulerTransversePacketProvider.HistoryData (D.initial τ )) (δ : ) ( : 0 < δ) (ξ : U) (hs : tsupport EulerSpatialCutoffs.innerCutoffD.support) (α : ) (Cagree : EulerPacketCylinderField.SourceCoefficientAgreement M D) (N : ) (hN : 1 N) (k : ) (hk : 4 k) (Q : EulerAllOrderDriftCorrection.Budget period (initializedCorrectionData M D hTime τ hτT B δ ξ hs α Cagree N hN k hk)) (t : (Set.Icc 0 D.T)) (Y : EulerSmoothLimit.SpaceEulerSmoothLimit.Space) :

Initialized exact physical pressure, given by (initializedExactPacket M D hTime τ hτ hτT B δ hδ ξ hs α Cagree N hN k hk Q).graphPotential k t ∘ Y.

Equations
  • One or more equations did not get rendered due to their size.
Instances For
    theorem EulerPacketTerminalDatum.initializedExactPhysicalPressure_gradient (M : EulerMeanPacketProvider.Data) {U : Type u_1} [NormedAddCommGroup U] [InnerProductSpace U] [CompleteSpace U] (D : EulerTransversePacketProvider.Data U) (hTime : M.T = D.T) (τ : ) ( : 0 < τ) (hτT : τ < D.T) (B : EulerTransversePacketProvider.HistoryData (D.initial τ )) (δ : ) ( : 0 < δ) (ξ : U) (hs : tsupport EulerSpatialCutoffs.innerCutoffD.support) (α : ) (Cagree : EulerPacketCylinderField.SourceCoefficientAgreement M D) (N : ) (hN : 1 N) (k : ) (hk : 4 k) (Q : EulerAllOrderDriftCorrection.Budget period (initializedCorrectionData M D hTime τ hτT B δ ξ hs α Cagree N hN k hk)) (X Y : (Set.Icc 0 D.T)EulerSmoothLimit.SpaceEulerSmoothLimit.Space) (hX : ∀ (t : (Set.Icc 0 D.T)) (x : EulerSmoothLimit.Space), HasFDerivAt (X t) ((D.F.field t) x) x) (hXY : ∀ (t : (Set.Icc 0 D.T)) (x : EulerSmoothLimit.Space), X t (Y t x) = x) (hY : Continuous (Function.uncurry Y)) (t : (Set.Icc 0 D.T)) (x : EulerSmoothLimit.Space) :
    gradient (initializedExactPhysicalPressure M D hTime τ hτT B δ ξ hs α Cagree N hN k hk Q t (Y t)) x = gradient (fun (y : EulerSmoothLimit.Space) => initializedPressure M D τ hτT B δ ξ hs α N k⁻¹ (t, Y t y, k * inner D.m₀ (Y t y))) x + gradient (EulerAllOrderDriftCorrection.Budget.physicalPotential D period Q k Y t) x
    theorem EulerPacketTerminalDatum.initializedExactPhysicalPressure_hessian (M : EulerMeanPacketProvider.Data) {U : Type u_1} [NormedAddCommGroup U] [InnerProductSpace U] [CompleteSpace U] (D : EulerTransversePacketProvider.Data U) (hTime : M.T = D.T) (τ : ) ( : 0 < τ) (hτT : τ < D.T) (B : EulerTransversePacketProvider.HistoryData (D.initial τ )) (δ : ) ( : 0 < δ) (ξ : U) (hs : tsupport EulerSpatialCutoffs.innerCutoffD.support) (α : ) (Cagree : EulerPacketCylinderField.SourceCoefficientAgreement M D) (N : ) (hN : 1 N) (k : ) (hk : 4 k) (Q : EulerAllOrderDriftCorrection.Budget period (initializedCorrectionData M D hTime τ hτT B δ ξ hs α Cagree N hN k hk)) (X Y : (Set.Icc 0 D.T)EulerSmoothLimit.SpaceEulerSmoothLimit.Space) (hX : ∀ (t : (Set.Icc 0 D.T)) (x : EulerSmoothLimit.Space), HasFDerivAt (X t) ((D.F.field t) x) x) (hXY : ∀ (t : (Set.Icc 0 D.T)) (x : EulerSmoothLimit.Space), X t (Y t x) = x) (hY : Continuous (Function.uncurry Y)) (t : (Set.Icc 0 D.T)) (x : EulerSmoothLimit.Space) :
    fderiv (gradient (initializedExactPhysicalPressure M D hTime τ hτT B δ ξ hs α Cagree N hN k hk Q t (Y t))) x = fderiv (gradient fun (y : EulerSmoothLimit.Space) => initializedPressure M D τ hτT B δ ξ hs α N k⁻¹ (t, Y t y, k * inner D.m₀ (Y t y))) x + fderiv (gradient (EulerAllOrderDriftCorrection.Budget.physicalPotential D period Q k Y t)) x
    theorem EulerPacketTerminalDatum.initializedExactPhysicalPressure_hessian_error (M : EulerMeanPacketProvider.Data) {U : Type u_1} [NormedAddCommGroup U] [InnerProductSpace U] [CompleteSpace U] (D : EulerTransversePacketProvider.Data U) (hTime : M.T = D.T) (τ : ) ( : 0 < τ) (hτT : τ < D.T) (B : EulerTransversePacketProvider.HistoryData (D.initial τ )) (δ : ) ( : 0 < δ) (ξ : U) (hs : tsupport EulerSpatialCutoffs.innerCutoffD.support) (α : ) (Cagree : EulerPacketCylinderField.SourceCoefficientAgreement M D) (N : ) (hN : 1 N) (k : ) (hk : 4 k) (Q : EulerAllOrderDriftCorrection.Budget period (initializedCorrectionData M D hTime τ hτT B δ ξ hs α Cagree N hN k hk)) (X Y : (Set.Icc 0 D.T)EulerSmoothLimit.SpaceEulerSmoothLimit.Space) (hX : ∀ (t : (Set.Icc 0 D.T)) (x : EulerSmoothLimit.Space), HasFDerivAt (X t) ((D.F.field t) x) x) (hXY : ∀ (t : (Set.Icc 0 D.T)) (x : EulerSmoothLimit.Space), X t (Y t x) = x) (hY : Continuous (Function.uncurry Y)) (L : EulerTransversePacketJoin.Budget D τ hτT B (Fin 4) 6) (H : EulerTransversePacketPrimary.Budget L) (NB : EulerTransversePacketJoin.NormalBudget D 6 L.R) (W : L.GradeGuards NB) (LM : EulerMeanPacketProvider.Budget M 6 L.R) (WM : LM.GradeGuards) (BC : EulerPacketCylinderField.CoefficientBudget (EulerPacketCylinderField.joinedSourceCoefficientData period M D τ hτT B hTime)) (hRc : EulerParameterWordGevrey.sobolevCoefficientRadius (Fin 4) BC.Rc L.R) (hcost : BC.termCost L.R) (hδ1 : δ 1) ( : 0 < α) (hR : wordRadius (Fin 4) δ L.R) (WP : H.GradeGuards NB (wordCost (Fin 4) 6 δ * ξ)) (S : EulerPacketTimeProfile.Scales (Set.Icc 0 M.T)) (hgrowth : EulerPacketCylinderField.timeProfileChange S.growth hTime = α L.fullProfile) (hbase : EulerPacketCoarseMajorant.tailBase L.R S.H0 BC.termCost N k ^ (1 / 100)) (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)) (x : EulerSmoothLimit.Space) :
    fderiv (gradient (initializedExactPhysicalPressure M D hTime τ hτT B δ ξ hs α Cagree N hN k hk Q t (Y t))) x - (EulerPacketPrimaryPressure.coefficient τ hτT B ξ hs α t (Y t x) * deriv (EulerPeriodicProfile.profile δ) (k * inner D.m₀ (Y t x))) ((InnerProductSpace.rankOne ) ((D.normal.field t) (Y t x))) ((D.normal.field t) (Y t x)) initializedPressureHessianCost NB L.R S.H0 L.Rc L.C₀ / k + fderiv (gradient (EulerAllOrderDriftCorrection.Budget.physicalPotential D period Q k Y t)) x