Globalization of the positive slow profiles #
The radial coordinate in this file is the physical similarity radius R,
so that X = R² / 2. Smoothness at the axis is represented by genuine
smooth even functions of the signed radius. Moment corrections are supported
in a fixed positive-radius patch and are reflected with the appropriate parity.
All integrals and differential operators below are the actual ones. In particular the preceding radial source is retained when pressure is recomputed.
Field: an abbreviation for ℝ × ℝ → ℝ.
Instances For
Dr, given by ProfileHistories.radialPartial.
Instances For
De, given by ProfileHistories.parameterPartial.
Instances For
Regular algebra, bundling carrier, zero_mem, one_mem, add_mem and the required
compatibility proofs.
Equations
- One or more equations did not get rendered due to their size.
Instances For
Equations
- NavierStokes.GlobalSlowProfiles.instCoeFunEvenProfileField = { coe := fun (f : NavierStokes.GlobalSlowProfiles.EvenProfile S) => ↑f }
Bound constant, given by algebraMap ℝ (EvenProfile S) c.
Equations
Instances For
Radius squared, given by ⟨fun w => w.1 ^ 2 / 2, ⟨(contDiffOn_fst.pow 2).div_const 2, fun _ _ _ => by simp only [neg_sq]⟩⟩.
Instances For
Parameter, given by ⟨Prod.snd, ⟨contDiffOn_snd, fun _ _ _ => rfl⟩⟩.
Equations
Instances For
Inverse, given by ⟨fun w => (f w)⁻¹, ⟨f.smooth.inv hf, fun _ heta R => congrArg Inv.inv (f.even heta R)⟩⟩.
Instances For
Eta derivative as an element of EvenProfile S.
Equations
Instances For
The genuine X derivative, smoothly continued through R = 0.
Equations
- NavierStokes.GlobalSlowProfiles.xDerivative hS f = ⟨fun (w : ℝ × ℝ) => 2 * NavierStokes.BoundaryAxisJets.radialJet (↑f) 1 w.1 w.2, ⋯⟩
Instances For
Literal X-only cutoff. The same two radii will be used at every order.
Equations
- NavierStokes.GlobalSlowProfiles.coreCutoff inner stop w = 1 - NavierStokes.TransportPrimitive.cutoff inner stop (w.1 ^ 2 / 2)
Instances For
A normalized mass average. There is no division by the radius.
Equations
Instances For
Domain data, collecting isOpen, denominator.
- isOpen : IsOpen S
- denominator (eta : ℝ) : eta ∈ S → PositiveAxisSystem.ell h eta ≠ 0
Instances For
Inverse denominator, given by inverse (1 - boundConstant S (2 * h) * parameter S ^ 2) (fun w hw => d.denominator w.2 hw.2).
Equations
- One or more equations did not get rendered due to their size.
Instances For
The normalized divergence reconstruction V = X β.
Equations
- One or more equations did not get rendered due to their size.
Instances For
Time op as an element of EvenProfile S.
Equations
- One or more equations did not get rendered due to their size.
Instances For
Axial op as an element of EvenProfile S.
Equations
- One or more equations did not get rendered due to their size.
Instances For
Shifted axial as an element of ℕ → EvenProfile S | 0 => 0 | k + 1 => axialOp2 d (AxisSourceRegularity.slowOrder h k - 1) (beta k).
Equations
- NavierStokes.GlobalSlowProfiles.shiftedAxial d beta 0 = 0
- NavierStokes.GlobalSlowProfiles.shiftedAxial d beta k.succ = NavierStokes.GlobalSlowProfiles.axialOp2 d (NavierStokes.AxisSourceRegularity.slowOrder h k - 1) (beta k)
Instances For
The complete radial source divided by X, including both viscous terms.
Equations
- One or more equations did not get rendered due to their size.
Instances For
Previous omega div X as an element of ℕ → EvenProfile S | 0 => 0 | k + 1 => omegaDivX d u beta k.
Equations
- NavierStokes.GlobalSlowProfiles.previousOmegaDivX d u beta 0 = 0
- NavierStokes.GlobalSlowProfiles.previousOmegaDivX d u beta k.succ = NavierStokes.GlobalSlowProfiles.omegaDivX d u beta k
Instances For
Angular field, given by w.1 / C * phi w.
Equations
- NavierStokes.GlobalSlowProfiles.angularField C phi w = w.1 / C * ↑phi w
Instances For
Pressure source, given by boundConstant S (C ^ 2)⁻¹ * (∑ i ∈ Finset.range (n + 1), phi i * phi (n - i)) - boundConstant S (1 / 2) * omega.
Equations
- One or more equations did not get rendered due to their size.
Instances For
Pressure from source, given by radiusSquared S * massAverage hS source.
Equations
Instances For
Pull a local even axis function to the physical similarity radius.
Instances For
Localizing never evaluates derivatives outside the common hierarchy radius: the cutoff is already zero on a neighborhood of that boundary.
Equations
- One or more equations did not get rendered due to their size.
Instances For
Odd angular correction divided by R. The support gap proves
smoothness at zero instead of appealing to a formal cancellation.
Equations
Instances For
Exterior: an abbreviation for SlowStressSupport.exterior.
Instances For
One repair step. Only actual smooth profiles and their compact exterior are inputs; the five zero moments are proved by the explicit linear repair.
Coefficient data, collecting phi, axial, beta, pressure.
- phi : EvenProfile S
Phi of
Coefficient, of typeEvenProfile S. - axial : EvenProfile S
Axial of
Coefficient, of typeEvenProfile S. - beta : EvenProfile S
Beta of
Coefficient, of typeEvenProfile S. - pressure : EvenProfile S
Pressure field of
Coefficient, of typeEvenProfile S.
Instances For
The finite order-zero data and a sequence of localized seeds. The
constructor schemeFromHierarchy below supplies the seeds from one actual
local hierarchy and one common cutoff.
- domain : Domain S h
- lam : ℝ
Lam of
Scheme, of typeℝ. - a : ℝ
A of
Scheme, of typeℝ. - b : ℝ
B of
Scheme, of typeℝ. - B : ℝ
Bound parameter of
Scheme, of typeℝ. Amplitude of
Scheme, of typeℝ → ℝ.- amplitude_smooth : ContDiffOn ℝ (↑⊤) self.amplitude S
- base : Coefficient S
Base of
Scheme, of typeCoefficient S. - seedAxial : ℕ → EvenProfile S
Seed axial of
Scheme, of typeℕ → EvenProfile S. - seedPhi : ℕ → EvenProfile S
Seed phi of
Scheme, of typeℕ → EvenProfile S.
Instances For
Admissible data, collecting data, axial_exterior, beta_exterior, phi_exterior,
pressure_exterior, zero_data and their compatibility conditions.
- data : Coefficient S
Data of
Admissible, of typeCoefficient S.
Instances For
Previous axial, with branches according to hj : j < n.
Equations
Instances For
Previous phi, with branches according to hj : j < n.
Equations
Instances For
Previous beta, with branches according to hj : j < n.
Equations
- NavierStokes.GlobalSlowProfiles.previousBeta s n previous j = if hj : j < n then (previous j hj).data.beta else 0
Instances For
Previous source, given by previousOmegaDivX s.domain (previousAxial s n previous) (previousBeta s n previous) n.
Equations
- One or more equations did not get rendered due to their size.
Instances For
Step properties, constructed using out.data.beta.
Equations
- One or more equations did not get rendered due to their size.
Instances For
Step, given by Classical.choose (exists_step s hn previous).
Equations
- NavierStokes.GlobalSlowProfiles.step s hn previous = Classical.choose ⋯
Instances For
Recursion step used in global slow profiles.
Equations
- One or more equations did not get rendered due to their size.
- NavierStokes.GlobalSlowProfiles.recursionStep s n.succ previous = NavierStokes.GlobalSlowProfiles.step s ⋯ previous
Instances For
A single coherent infinite sequence, constructed by well-founded recursion; all coefficient indices used in the radial source are smaller.
Equations
Instances For
All five actual total rows hold at every positive order of the same infinite sequence. The source is recomputed from that sequence's lower fields.
Data required from the actual order-zero construction. No positive-order equations, moments, or extension properties are assumed here.
- fields : Coefficient S
Fields of
BaseData, of typeCoefficient S. Amplitude of
BaseData, of typeℝ → ℝ.- amplitude_smooth : ContDiffOn ℝ (↑⊤) self.amplitude S
Instances For
The common cutoff is chosen strictly inside the hierarchy's proved radius; the moment patch begins beyond that cutoff.
Equations
- One or more equations did not get rendered due to their size.
Instances For
The actual squared-radius profile. We use ordinary derivatives only at positive X; the axis jets are the right jets of this descent.
Instances For
A transparent interface supplied by the common cutoff constructor.
- parameter_open : IsOpen U
Instances For
Agreement of the finite order-zero input with the same local hierarchy. These are identities on an inner region, not positive-order hypotheses.
Instances For
An actual substitution in the mass integral, with no axis differentiability assumption on the squared-radius function.
Local extension, constructed using cutoffLift.
Equations
- NavierStokes.GlobalSlowProfiles.localExtension L n i = NavierStokes.GlobalSlowProfiles.cutoffLift ⋯ ⋯ ⋯ ⋯ ⋯ ⋯ (A.coefficients n i)
Instances For
As slow profiles, constructed using SlowResidualMatching.ofBeta.
Equations
- One or more equations did not get rendered due to their size.
Instances For
Component, given by ![q.phi, q.axial, q.beta, q.pressure] i.
Instances For
Every repaired coefficient agrees with the original local hierarchy on one common inner domain, including the axis.
Joint smoothness up to the axis is inherited from actual local equality.
All actual mixed right jets, not just the zero-order values, agree at the axis. The derivative is taken relative to the closed half-plane.
Germ locality transfers the two solved tangential equations. No agreement outside the fixed inner region is required.