Documentation

LeanPool.NavierStokesAndEuler.Euler.TransverseNormalResidual

Recovering the transverse normal residual #

A residual annihilated by the adjoint of a frame spanning m⊥ is exactly its normal component. Applying this elementary Hilbert-space fact to the proved projected coordinate equation gives the pressure coefficient in equation (11).

The literal scalar normal component, with the source denominator D_m.

Equations
Instances For
    theorem EulerTransverseNormalResidual.eq_normal_of_tangent_orthogonal {E : Type u_1} [NormedAddCommGroup E] [InnerProductSpace E] (m r : E) (hm : m 0) (hr : ∀ (x : E), inner m x = 0inner r x = 0) :

    Orthogonality to the tangent hyperplane determines the normal component.

    theorem EulerTransverseNormalResidual.eq_normal_of_adjoint_zero {E : Type u_1} [NormedAddCommGroup E] [InnerProductSpace E] {U : Type u_2} [NormedAddCommGroup U] [InnerProductSpace U] [CompleteSpace U] [CompleteSpace E] (Q : U →L[] E) (m r : E) (hm : m 0) (hRange : ∀ (η : E), inner m η = 0∃ (x : U), Q x = η) (hr : (ContinuousLinearMap.adjoint Q) r = 0) :

    A frame spanning the tangent hyperplane detects exactly the normal residual.

    theorem EulerTransverseNormalResidual.physical_velocity_balance {E : Type u_1} [NormedAddCommGroup E] [InnerProductSpace E] {U : Type u_2} [NormedAddCommGroup U] [InnerProductSpace U] [CompleteSpace U] [CompleteSpace E] (Q Q₁ : U →L[] E) (M : E →L[] E) (m : E) (hm : m 0) (hTangent : ∀ (x : U), inner m (Q x) = 0) (hRange : ∀ (η : E), inner m η = 0∃ (x : U), Q x = η) (hflow : Q₁ = M ∘SL Q) (v a : U) (f : E) (heq : (EulerTransverseGramInverse.gram Q) a = (ContinuousLinearMap.adjoint Q) (f - 2 Q₁ v)) :
    Q₁ v + Q a + M (Q v) + ((inner m f - 2 * inner m (M (Q v))) / m ^ 2) m = f

    Equation (10), together with Q_t = M Q, gives the exact source pressure coefficient and velocity equation (11). No pressure residual is assumed.