Documentation

LeanPool.NavierStokesAndEuler.Euler.ParentStateGeometry

Geometric renewal between two actual Euler states. Their spatial regularity, particle velocity laws and fixed origin are supplied by the states themselves; only the quantitative packet expansion is an input.

The literal primary terms in source (20) initialize the next geometric parent frame. The forward coordinate is prescribed; the joined coordinate is the genuine stationary history trace and is proved nonzero.

The next source's center expansion follows from the actual physical velocity update. Odd particle displacements fix the origin, and the two literal velocity laws identify the source matrices there.

theorem EulerParentPacketFrames.Parent.position_zero_of_odd (G : Parent) (hodd : ∀ (t : (Set.Icc 0 G.T)), Function.Odd (G.displacement.field t)) (t : (Set.Icc 0 G.T)) :
G.position t 0 = 0
theorem EulerParentPacketFrames.Parent.center_update_of_odd (G N : Parent) (hTime : N.T = G.T) (u w : EulerSmoothLimit.SpaceEulerSmoothLimit.Space) (hu : ∀ (t : (Set.Icc 0 G.T)) (x : EulerSmoothLimit.Space), DifferentiableAt (u t) x) (hw : ∀ (t : (Set.Icc 0 N.T)) (x : EulerSmoothLimit.Space), DifferentiableAt (w t) x) (hGvelocity : ∀ (t : (Set.Icc 0 G.T)) (x : EulerSmoothLimit.Space), (G.velocity.field t) x = u (↑t) (G.position t x)) (hNvelocity : ∀ (t : (Set.Icc 0 N.T)) (x : EulerSmoothLimit.Space), (N.velocity.field t) x = u (↑t) (N.position t x) + w (↑t) (N.position t x)) (hGodd : ∀ (t : (Set.Icc 0 G.T)), Function.Odd (G.displacement.field t)) (hNodd : ∀ (t : (Set.Icc 0 N.T)), Function.Odd (N.displacement.field t)) (t : (Set.Icc 0 N.T)) :
(N.strain.field t) 0 = G.centerStrain t + fderiv (w t) 0
noncomputable def EulerParentPacketFrames.Parent.geometryFrameOfPhysicalUpdate (G N : Parent) (hTime : N.T = G.T) (u w : EulerSmoothLimit.SpaceEulerSmoothLimit.Space) (hu : ∀ (t : (Set.Icc 0 G.T)) (x : EulerSmoothLimit.Space), DifferentiableAt (u t) x) (hw : ∀ (t : (Set.Icc 0 N.T)) (x : EulerSmoothLimit.Space), DifferentiableAt (w t) x) (hGvelocity : ∀ (t : (Set.Icc 0 G.T)) (x : EulerSmoothLimit.Space), (G.velocity.field t) x = u (↑t) (G.position t x)) (hNvelocity : ∀ (t : (Set.Icc 0 N.T)) (x : EulerSmoothLimit.Space), (N.velocity.field t) x = u (↑t) (N.position t x) + w (↑t) (N.position t x)) (hGodd : ∀ (t : (Set.Icc 0 G.T)), Function.Odd (G.displacement.field t)) (hNodd : ∀ (t : (Set.Icc 0 N.T)), Function.Odd (N.displacement.field t)) {U : Type u_1} [NormedAddCommGroup U] [InnerProductSpace U] [CompleteSpace U] (m : EulerSmoothLimit.Space) (hm : m = 1) (R : U ≃ₗᵢ[] (EulerTransverseFrameCoordinates.referencePlane m)) (S : Set EulerSmoothLimit.Space) (hS : IsCompact S) {V : Type u_2} [NormedAddCommGroup V] [InnerProductSpace V] (mNew : EulerSmoothLimit.Space) (hmNew : mNew = 1) (RNew : V ≃ₗᵢ[] (EulerTransverseFrameCoordinates.referencePlane mNew)) (SNew : Set EulerSmoothLimit.Space) (hSNew : IsCompact SNew) (τ : ) ( : 0 τ) (η : U) ( : η 0) (c CM CH K error : ) (hCM : 0 CM) (hK : 1 K) (he : 0 error) (hMK : CM K) (hHK : CM ^ 2 + CH K ^ 2) (hM : tSet.Icc τ N.T, G.centerStrain t CM) (hH : tSet.Icc τ N.T, G.centerCurvature t CH) (hpacket : tSet.Icc τ N.T, fderiv (w t) 0 - c ((InnerProductSpace.rankOne ) (G.sourceVelocity m hm R S hS η t)) (G.sourceNormal m t) error) :
EulerPacketSourceGeometry.ParentFrame (N.transverseData mNew hmNew RNew SNew hSNew) τ

Concrete old/new-parent factory. The source matrix split is derived from the two physical velocity identities. The remaining quantitative inputs are the preceding packet's center shear error and parent low norm bounds, as used by the geometric induction.

Equations
  • One or more equations did not get rendered due to their size.
Instances For
    theorem EulerParentPacketFrames.Parent.joined_primary_center_term (G : Parent) {U : Type u_1} [NormedAddCommGroup U] [InnerProductSpace U] [CompleteSpace U] (m : EulerSmoothLimit.Space) (hm : m = 1) (R : U ≃ₗᵢ[] (EulerTransverseFrameCoordinates.referencePlane m)) (S : Set EulerSmoothLimit.Space) (hS : IsCompact S) (s : ) (hs : 0 < s) (hsT : s < G.T) (B : EulerTransversePacketProvider.HistoryData ((G.transverseData m hm R S hS).initial s hs )) (ξ : U) (hcut : tsupport EulerSpatialCutoffs.innerCutoffS) (δ : ) ( : 0 < δ) (α k : ) (t : (Set.Icc 0 G.T)) (Y : EulerSmoothLimit.SpaceEulerSmoothLimit.Space) (hY : Y 0 = 0) :
    (α * deriv (EulerPeriodicProfile.profile δ) (k * inner m (Y 0))) ((InnerProductSpace.rankOne ) (EulerPacketPrimaryFactorization.canonicalVelocity s hs hsT B ξ hcut (↑t) (Y 0))) (((G.transverseData m hm R S hS).normal.field t) (Y 0)) = (α / δ) ((InnerProductSpace.rankOne ) (G.sourceVelocity m hm R S hS (((B.coefficients.labelCoordinate 0) ξ) 0, ) t)) (G.sourceNormal m t)
    noncomputable def EulerParentPacketFrames.Parent.forwardGeometryFrame (G : Parent) {U : Type u_1} [NormedAddCommGroup U] [InnerProductSpace U] [CompleteSpace U] (m : EulerSmoothLimit.Space) (hm : m = 1) (R : U ≃ₗᵢ[] (EulerTransverseFrameCoordinates.referencePlane m)) (S : Set EulerSmoothLimit.Space) (hS : IsCompact S) (N : Parent) (hTime : N.T = G.T) (u w : EulerSmoothLimit.SpaceEulerSmoothLimit.Space) (hu : ∀ (t : (Set.Icc 0 G.T)) (x : EulerSmoothLimit.Space), DifferentiableAt (u t) x) (hw : ∀ (t : (Set.Icc 0 N.T)) (x : EulerSmoothLimit.Space), DifferentiableAt (w t) x) (hGvelocity : ∀ (t : (Set.Icc 0 G.T)) (x : EulerSmoothLimit.Space), (G.velocity.field t) x = u (↑t) (G.position t x)) (hNvelocity : ∀ (t : (Set.Icc 0 N.T)) (x : EulerSmoothLimit.Space), (N.velocity.field t) x = u (↑t) (N.position t x) + w (↑t) (N.position t x)) (hGodd : ∀ (t : (Set.Icc 0 G.T)), Function.Odd (G.displacement.field t)) (hNodd : ∀ (t : (Set.Icc 0 N.T)), Function.Odd (N.displacement.field t)) {V : Type u_2} [NormedAddCommGroup V] [InnerProductSpace V] (mNew : EulerSmoothLimit.Space) (hmNew : mNew = 1) (RNew : V ≃ₗᵢ[] (EulerTransverseFrameCoordinates.referencePlane mNew)) (SNew : Set EulerSmoothLimit.Space) (hSNew : IsCompact SNew) (τ : ) ( : 0 τ) (CM CH K error : ) (hCM : 0 CM) (hK : 1 K) (he : 0 error) (hMK : CM K) (hHK : CM ^ 2 + CH K ^ 2) (hM : tSet.Icc τ N.T, G.centerStrain t CM) (hH : tSet.Icc τ N.T, G.centerCurvature t CH) (Y : (Set.Icc 0 G.T)EulerSmoothLimit.SpaceEulerSmoothLimit.Space) (hY : ∀ (t : (Set.Icc 0 G.T)), Y t 0 = 0) (δ : ) ( : 0 < δ) (α k : ) (η : U) ( : η 0) (hsource20 : ∀ (t : (Set.Icc 0 G.T)), τ tfderiv (w t) 0 - (α * deriv (EulerPeriodicProfile.profile δ) (k * inner m (Y t 0))) ((InnerProductSpace.rankOne ) (EulerPacketForwardFactorization.canonicalVelocity (G.transverseData m hm R S hS) η (↑t) (Y t 0))) (((G.transverseData m hm R S hS).normal.field t) (Y t 0)) error) :
    EulerPacketSourceGeometry.ParentFrame (N.transverseData mNew hmNew RNew SNew hSNew) τ

    Forward geometry frame as an element of ParentFrame (N.transverseData mNew hmNew RNew SNew hSNew) τ.

    Equations
    • One or more equations did not get rendered due to their size.
    Instances For
      noncomputable def EulerParentPacketFrames.Parent.joinedGeometryFrame (G : Parent) {U : Type u_1} [NormedAddCommGroup U] [InnerProductSpace U] [CompleteSpace U] (m : EulerSmoothLimit.Space) (hm : m = 1) (R : U ≃ₗᵢ[] (EulerTransverseFrameCoordinates.referencePlane m)) (S : Set EulerSmoothLimit.Space) (hS : IsCompact S) (N : Parent) (hTime : N.T = G.T) (u w : EulerSmoothLimit.SpaceEulerSmoothLimit.Space) (hu : ∀ (t : (Set.Icc 0 G.T)) (x : EulerSmoothLimit.Space), DifferentiableAt (u t) x) (hw : ∀ (t : (Set.Icc 0 N.T)) (x : EulerSmoothLimit.Space), DifferentiableAt (w t) x) (hGvelocity : ∀ (t : (Set.Icc 0 G.T)) (x : EulerSmoothLimit.Space), (G.velocity.field t) x = u (↑t) (G.position t x)) (hNvelocity : ∀ (t : (Set.Icc 0 N.T)) (x : EulerSmoothLimit.Space), (N.velocity.field t) x = u (↑t) (N.position t x) + w (↑t) (N.position t x)) (hGodd : ∀ (t : (Set.Icc 0 G.T)), Function.Odd (G.displacement.field t)) (hNodd : ∀ (t : (Set.Icc 0 N.T)), Function.Odd (N.displacement.field t)) {V : Type u_2} [NormedAddCommGroup V] [InnerProductSpace V] (mNew : EulerSmoothLimit.Space) (hmNew : mNew = 1) (RNew : V ≃ₗᵢ[] (EulerTransverseFrameCoordinates.referencePlane mNew)) (SNew : Set EulerSmoothLimit.Space) (hSNew : IsCompact SNew) (τ : ) ( : 0 τ) (CM CH K error : ) (hCM : 0 CM) (hK : 1 K) (he : 0 error) (hMK : CM K) (hHK : CM ^ 2 + CH K ^ 2) (hM : tSet.Icc τ N.T, G.centerStrain t CM) (hH : tSet.Icc τ N.T, G.centerCurvature t CH) (Y : (Set.Icc 0 G.T)EulerSmoothLimit.SpaceEulerSmoothLimit.Space) (hY : ∀ (t : (Set.Icc 0 G.T)), Y t 0 = 0) (δ : ) ( : 0 < δ) (α k s : ) (hs : 0 < s) (hsT : s < G.T) (B : EulerTransversePacketProvider.HistoryData ((G.transverseData m hm R S hS).initial s hs )) (ξ : U) ( : ξ 0) (hcut : tsupport EulerSpatialCutoffs.innerCutoffS) (hsource20 : ∀ (t : (Set.Icc 0 G.T)), τ tfderiv (w t) 0 - (α * deriv (EulerPeriodicProfile.profile δ) (k * inner m (Y t 0))) ((InnerProductSpace.rankOne ) (EulerPacketPrimaryFactorization.canonicalVelocity s hs hsT B ξ hcut (↑t) (Y t 0))) (((G.transverseData m hm R S hS).normal.field t) (Y t 0)) error) :
      EulerPacketSourceGeometry.ParentFrame (N.transverseData mNew hmNew RNew SNew hSNew) τ

      Joined geometry frame as an element of ParentFrame (N.transverseData mNew hmNew RNew SNew hSNew) τ.

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

        Velocity increment, given by T.evolution.velocity (t,x)-S.evolution.velocity (t,x).

        Equations
        Instances For
          noncomputable def EulerParentPacketFrames.SmoothState.forwardRenewal {A N : Parent} (S : SmoothState A) (T : SmoothState N) (hTime : N.T = A.T) {U : Type u_1} [NormedAddCommGroup U] [InnerProductSpace U] [CompleteSpace U] (m : EulerSmoothLimit.Space) (hm : m = 1) (J : U ≃ₗᵢ[] (EulerTransverseFrameCoordinates.referencePlane m)) (support : Set EulerSmoothLimit.Space) (hSupport : IsCompact support) {V : Type u_2} [NormedAddCommGroup V] [InnerProductSpace V] (mNext : EulerSmoothLimit.Space) (hmNext : mNext = 1) (JNext : V ≃ₗᵢ[] (EulerTransverseFrameCoordinates.referencePlane mNext)) (supportNext : Set EulerSmoothLimit.Space) (hSupportNext : IsCompact supportNext) (τ : ) ( : 0 τ) (CM CH K error : ) (hCM : 0 CM) (hK : 1 K) (he : 0 error) (hMK : CM K) (hHK : CM ^ 2 + CH K ^ 2) (hM : tSet.Icc τ N.T, A.centerStrain t CM) (hH : tSet.Icc τ N.T, A.centerCurvature t CH) (δ : ) ( : 0 < δ) (α k : ) (ξ : U) ( : ξ 0) (hsource : ∀ (t : (Set.Icc 0 A.T)), τ tfderiv (S.velocityIncrement T t) 0 - (α * deriv (EulerPeriodicProfile.profile δ) (k * inner m (S.evolution.inverse.normalized t 0))) ((InnerProductSpace.rankOne ) (EulerPacketForwardFactorization.canonicalVelocity (A.transverseData m hm J support hSupport) ξ (↑t) (S.evolution.inverse.normalized t 0))) (((A.transverseData m hm J support hSupport).normal.field t) (S.evolution.inverse.normalized t 0)) error) :
          EulerPacketSourceGeometry.ParentFrame (N.transverseData mNext hmNext JNext supportNext hSupportNext) τ

          Forward renewal, constructed using A.forwardGeometryFrame.

          Equations
          • One or more equations did not get rendered due to their size.
          Instances For
            noncomputable def EulerParentPacketFrames.SmoothState.joinedRenewal {A N : Parent} (S : SmoothState A) (T : SmoothState N) (hTime : N.T = A.T) {U : Type u_1} [NormedAddCommGroup U] [InnerProductSpace U] [CompleteSpace U] (m : EulerSmoothLimit.Space) (hm : m = 1) (J : U ≃ₗᵢ[] (EulerTransverseFrameCoordinates.referencePlane m)) (support : Set EulerSmoothLimit.Space) (hSupport : IsCompact support) {V : Type u_2} [NormedAddCommGroup V] [InnerProductSpace V] (mNext : EulerSmoothLimit.Space) (hmNext : mNext = 1) (JNext : V ≃ₗᵢ[] (EulerTransverseFrameCoordinates.referencePlane mNext)) (supportNext : Set EulerSmoothLimit.Space) (hSupportNext : IsCompact supportNext) (τ : ) ( : 0 τ) (CM CH K error : ) (hCM : 0 CM) (hK : 1 K) (he : 0 error) (hMK : CM K) (hHK : CM ^ 2 + CH K ^ 2) (hM : tSet.Icc τ N.T, A.centerStrain t CM) (hH : tSet.Icc τ N.T, A.centerCurvature t CH) (δ : ) ( : 0 < δ) (α k s : ) (hs : 0 < s) (hsT : s < A.T) (H : EulerTransversePacketProvider.HistoryData ((A.transverseData m hm J support hSupport).initial s hs )) (ξ : U) ( : ξ 0) (hcut : tsupport EulerSpatialCutoffs.innerCutoffsupport) (hsource : ∀ (t : (Set.Icc 0 A.T)), τ tfderiv (S.velocityIncrement T t) 0 - (α * deriv (EulerPeriodicProfile.profile δ) (k * inner m (S.evolution.inverse.normalized t 0))) ((InnerProductSpace.rankOne ) (EulerPacketPrimaryFactorization.canonicalVelocity s hs hsT H ξ hcut (↑t) (S.evolution.inverse.normalized t 0))) (((A.transverseData m hm J support hSupport).normal.field t) (S.evolution.inverse.normalized t 0)) error) :
            EulerPacketSourceGeometry.ParentFrame (N.transverseData mNext hmNext JNext supportNext hSupportNext) τ

            Joined renewal, constructed using A.joinedGeometryFrame.

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