epagoge
Numerics note / 28 July 2026

The 2D Boussinesq corner blowup profile

Four convergence tests passed. The root was false.

Justin Hill, Independent Research. An independent computation of the self-similar profile underlying the Luo and Hou scenario for 3D axisymmetric Euler with boundary, solved as a root problem with an explicit sparse Jacobian rather than by marching or neural approximation.

Read the note (PDF) Code, logs, fields
New study · 13 August 2026
The viscous response is linear at leading order
Paper three: rho_bulk marches 1.243 to 1.029 toward a ν → 0 intercept of 1.0035 ± 0.0099, no finite viscous anomaly, three wall laws, and flux observables whose inviscid null is zero by theorem. Three of our own results killed en route.
Read paper three →
Study · 2 August 2026
Viscosity inverts the corner-flow blowup geometry
The direction-regularity exponent this note's profile made measurable, calibrated at +1.00 inviscid: window-matched, +0.59 becomes −1.2 deep in the viscous collapse. Below the type-I exclusion threshold, grid-certified.
Read the study →
The transferable result · one artifact · one functional

Deflated multistart at a frozen unstable exponent returned a state that satisfied every check the campaign had been applying. Under grid refinement its exponent moved away from its target with growing steps. It has no continuum limit.

5.4×10−14
Converged
Newton residual at machine precision.
1.70
Distinct
Relative distance from the ground state under deflation.
1.3×10−5
Parameter-stable
Motion across a decade of wedge truncation, where the true profile moves 2.72×10−3.
6172×
Morphologically coherent
Its own angular signature, against the ground branch's variation.
+0.88
Free residual
Against −0.00106 for the true profile. Three orders of magnitude, on the one functional the solve is not answerable to.

The gauge closure fixes two corner constants and leaves the corner identity c_l = 2Θxx / Wx unused. Reporting that leftover beside every converged solve costs one line of output and separates the cases the other four tests could not. Residual, distinctness, parameter stability and morphological coherence are jointly insufficient to certify a self-similar profile.

The exponent

Chen and Hou solve this profile by adaptive-mesh march (arXiv:2210.07191, arXiv:2305.05660). Wang and collaborators solve it with a physics-informed network refined by Gauss-Newton (arXiv:2509.14185). They share no discretization and agree to 3.4×10−7, which is why their value is treated here as reference rather than competitor.

Scaling exponent · five-rung extrapolation · model-class spread as the bar
−0.34240 ± 4.4×10−5 reference −0.34240009

Three extrapolation families fit the ladder with residuals five points cannot separate. Their spread sets the bar, added in quadrature with the measured corner-degree systematic. Quoting the deepest three rungs alone would report agreement at 2.5×10−6. That number is not quoted, because choosing the subset that tightens the bar is the same error as choosing the model family that does.

Independence is not total, and the note states the exception in its abstract. Two scalar corner constants are read off the reference profile and imposed as gauge targets. Freeing one moves the exponent by 9.7×10−5, about twice the bar. The claim is independence of method, not of data.

AxisWorst effect on the exponentStanding
Angular resolution, 36 to 481.16×10−8closed
Seed provenance, reference to three analytic9.65×10−8closed
Far-field truncation, 25 to 321.45×10−10closed · grids verified to differ
Corner degree, 24 to 281.94×10−5carried in the bar
Corner degree, 16 to 246.98×10−4not converged at 16
Wedge truncationextrapolated, layer analysedcarried in the bar

What went wrong on the way

An adversarial verification pass re-sourced every number in the draft against its log. It found three things worth publishing alongside the result.

Corrected
A mislabelled abscissa
A rung run at 2.5×10−5 printed as 3×10−5 by a rounding format, then fitted at the printed value. Correcting it tightened the spread and moved the quoted centre onto the reference.
Narrowed
An independence claim
The draft said the method shares no machinery with prior work. The solver hardcodes two of their corner constants. The claim now reads independence of method, and the cost of the shared data is measured.
Deleted
A model comparison that cannot exist
The draft called three fits AICc-indistinguishable. AICc is not computable on five points with three parameters. The statistic is gone and the spread is quoted directly.

Four retractions made during the work are recorded in the note rather than removed, including a spectral mode built on unconverged eigenvalues and a corner-angle derivative withdrawn when its slope proved formulation-dependent. The claims ledger grades every statement in the paper against its evidence, and no claim appears in the paper above its grade there.

Still open

The unstable branches remain unconfirmed by any method other than the network that found them. The route attempted here produced the artifact above. Treating the exponent as a Newton unknown on the untruncated problem is the next attempt, and it is not built.

A geometric reading of the published branch gaps puts the family's accumulation point at −1/2, the classical Leray exponent, and predicts a fourth branch near −0.4649. The network's own empirical law gives −0.46557 for the same branch. Both are extrapolations over the same four numbers, and nobody has measured the fourth reliably yet.

Share card

The image social platforms show for this page. 1200×630.

Four convergence tests passed, one free residual failed