Regularita přes strukturální zpětnou vazbu

Kde jsme skončili

Svitek 007 ustanovil tři věci. Za prvé: Millennium Prize se ptá na matematickou vlastnost konkrétní rovnice, ne na fyziku tekutin. Za druhé: ta rovnice obsahuje pouze pasivní disipaci (viskózní člen μ∇²v) bez strukturální zpětné vazby reagující na vírovitost. Za třetí: operátor ακ je kandidát na chybějící člen — energeticky konzistentní, strukturálně selektivní, spektrálně odlišný od viskozity.

Co jsme neudělali: nedefinovali jsme operátor formálně. Nepopsali jsme, jaké vlastnosti musí mít obecně, aby plnil stabilizační roli. Neformulovali jsme enstrofickou nerovnost explicitně. To dělá tento svitek.

Axiomatická definice: třída strukturálně stabilizačních operátorů

Místo toho, abychom prezentovali jeden konkrétní operátor, definujeme třídu. To je silnější pozice — netvrdíme 'tohle je ten správný člen,' tvrdíme 'jakýkoliv člen splňující tyto axiomy by měl stabilizační efekt.'

Nechť F: H¹(ℝ³) → L²(ℝ³) je operátor přiřazující rychlostnímu poli v korekční sílu f = F(v). Řekneme, že F je strukturálně stabilizační operátor, pokud splňuje čtyři axiomy:

Axiom 1 — Divergenční kompatibilita: ∇ · F(v) = 0 pro všechna v splňující ∇ · v = 0. Operátor nesmí narušit nestlačitelnost. Jakýkoliv korekční člen, který generuje divergentní složku, by porušil zachování hmoty. To je nepřekročitelná podmínka.

Axiom 2 — Energetická neutralita v rovnovážném stavu: ∫ F(v) · v dV = 0 pro v homogenní nebo v rovnovážné konfiguraci. Operátor nepřidává čistou energii do systému v rovnováze. Může redistribuovat energii mezi škálami nebo směry, ale celková energetická bilance zůstává zachována. V nehomogenním, nestabilním poli může být ∫ F · v dV ≠ 0 lokálně — ale celkový příspěvek musí být omezený a stabilizující (tedy záporný v oblastech nestability).

Axiom 3 — Subkritické škálování: Pod škálováním NS (v_λ(x,t) = λv(λx, λ²t)) se operátor transformuje nejvýše jako: ||F(v_λ)|| ≤ C λ^α ||F(v)|| kde α < 2. Viskózní člen μ∇²v škáluje přesně jako λ². Vortex stretching (ω·∇)v škáluje taky jako λ². Proto jsou standardní NS kritické — oba členy rostou stejně rychle a není jasné, který dominuje. Operátor F musí škálovat pomaleji než λ² (subkriticky), aby nepřidal další kritický člen, ale zároveň musí být dostatečně silný, aby v kritických oblastech dominoval nad vortex stretchingem.

Axiom 4 — Dominance v kritickém režimu: Existuje funkcionál Φ: H¹ → ℝ⁺ takový, že (a) Φ(v) ≥ 0 pro všechna v, (b) Φ(v) = 0 právě tehdy, když ω = ∇ × v = 0, (c) V oblastech, kde ||ω·S·ω|| → ∞ (intenzivní vortex stretching), platí: |⟨F(v), ω⟩| ≥ C_κ Φ(ω) kde C_κ > 0. To říká: tam, kde vortex stretching hrozí blowupem, tam operátor F generuje dostatečně silnou protisílu. Ne všude — jen tam, kde je to potřeba.

Čtyři axiomy definují třídu operátorů, které zachovávají nestlačitelnost, nepřidávají energii, škálují subkriticky a aktivně stabilizují v oblastech hrozícího blowupu.

Proč právě tyto čtyři axiomy

Každý axiom odpovídá jednomu ze známých problémů regularity:

Axiom 1 zaručuje, že rozšířená rovnice je stále rovnicí pro nestlačitelnou tekutinu. Bez něj bychom řešili jiný fyzikální problém.

Axiom 2 zaručuje, že operátor neřeší problém triviálně — nepřidává energii, která by 'přebila' vortex stretching zvenčí. To by byl podvod. Operátor musí stabilizovat strukturou, ne silou.

Axiom 3 reaguje přímo na Taovu supercriticality barrier. Tao ukázal, že energetická identita nestačí k důkazu regularity, protože NS jsou kritické (škálování viskozity a nelinearity je stejné). Přidaný operátor nesmí tuto kritikalitu zhoršit — musí ji zmírnit. Subkritické škálování to zajišťuje.

Axiom 4 je jádro celé konstrukce. Říká, že operátor není jen pasivní filtr (jako viskozita), ale aktivní stabilizátor, který se zapíná přesně v okamžiku, kdy hrozí blowup. Viskozita filtruje k² uniformně. Operátor F filtruje selektivně tam, kde je ω·S·ω velké.

Modifikovaná enstrofická nerovnost

Standardní enstrofická nerovnost pro 3D NS:

dΩ/dt = ∫ ω·S·ω dV − ν||∇ω||²

Problém: první člen (produkce) může dominovat nad druhým (disipace). Nerovnost nezaručuje, že Ω zůstane omezená.

S operátorem F splňujícím axiomy 1–4 dostáváme modifikovanou nerovnost:

dΩ/dt ≤ ∫ ω·S·ω dV − ν||∇ω||² − C_κ ∫ Φ(ω) dV

Třetí člen je nový. Je záporný (stabilizační) a aktivuje se přes funkcionál Φ, který je svázaný s intenzitou vírovitosti.

Klíčová otázka (otevřená): Existuje operátor F splňující axiomy 1–4 takový, že:

∫ ω·S·ω dV ≤ ν||∇ω||² + C_κ ∫ Φ(ω) dV

pro všechna t ≥ 0?

Pokud ano, pak Ω je omezená pro všechny časy, což implikuje regularitu modifikovaných rovnic.

Netvrdíme, že odpověď je 'ano.' Tvrdíme, že toto je správná otázka.

Modifikovaná enstrofická nerovnost s korekčním členem −C_κ ∫ Φ(ω) dV dává konkrétní otevřenou otázku: zda tento člen dokáže dominovat nad vortex stretchingem v kritickém režimu.

Škálovací analýza

NS rovnice mají přirozenou škálovací symetrii. Pokud v(x,t) je řešení, pak pro libovolné λ > 0 je v_λ(x,t) = λv(λx, λ²t) taky řešení (s přeškálovaným tlakem). Tato symetrie je klíčová pro problém regularity — blowup by musel respektovat škálování.

Jak se jednotlivé členy chovají pod škálováním v prostoru L²:

Konvektivní člen (v·∇)v: škáluje jako λ². Kritický.

Viskózní člen ν∇²v: škáluje jako λ². Kritický. (Proto je soutěž nerozhodná.)

Operátor F: škáluje jako λ^α, α < 2. Subkritický.

Na první pohled to vypadá jako nevýhoda — subkritický člen je 'slabší' než kritické členy. Ale to je přesně pointa. Operátor F nemusí dominovat všude. Musí dominovat jen v kritických oblastech (axiom 4), a tam ho lokální koncentrace enstrofie zesiluje přes funkcionál Φ. Na velkých škálách je zanedbatelný. Na malých škálách, kde hrozí blowup, se zapíná.

Toto chování je kvalitativně odlišné od viskozity. Viskozita je uniformně kritická — stejně silná všude. Operátor F je lokálně adaptivní — slabý kde není potřeba, silný kde je.

Operátor F škáluje subkriticky (λ^α, α < 2), čímž nezvyšuje kritikalitu NS, ale lokálně se zesiluje přes funkcionál Φ v oblastech intenzivní vírovitosti.

Napojení na Taovu bariéru

Tao (2016) dokázal: pro 'averaged NS' (kde je nelinearita nahrazena obecnějším operátorem zachovávajícím energetickou identitu) existuje finite-time blowup.

Co z toho plyne: jakýkoliv důkaz regularity skutečných NS musí využít specifickou strukturu nelinearity B(v,v) = (v·∇)v, ne jen její energetické vlastnosti.

Operátor F je přesně tento typ strukturální informace.

Jak? Tao ukazuje, že energetická identita ⟨B(v,v), v⟩ = 0 nestačí. Ale operátor F definovaný axiomy 1–4 reaguje na vnitřní strukturu B(v,v) — konkrétně na tensor deformační rychlosti S a jeho interakci s vírovitostí ω. To je informace, kterou energetická identita neobsahuje, ale Taovy 'averaged NS' ztrácejí.

Řečeno jinak: Tao říká 'potřebujete vědět víc než energii.' Axiomy 1–4 definují jaký typ informace potřebujete — strukturální zpětnou vazbu na ω·S·ω.

To neznamená, že máme důkaz. Znamená to, že máme strategii — a ta strategie přímo reaguje na nejsilnější negativní výsledek v oboru.

Axiomy 1–4 definují přesně typ strukturální informace, který Tao identifikoval jako nutný pro překonání supercriticality barrier — zpětnou vazbu na ω·S·ω.

Dvě interpretace (znovu, přesněji)

Ze svitku 007 zůstávají dvě interpretace operátoru F. V axiomatickém rámci je lze formulovat přesněji:

Interpretace A: F jako izolace existující struktury NS. Lawson et al. (2022) pozorovali samoregularizační anti-twist ve standardních NS. Pokud existuje operátor F splňující axiomy 1–4, který je důsledkem nelinearity B(v,v) — tedy pokud lze z (v·∇)v algebraicky izolovat člen s vlastnostmi axiomů — pak standardní NS už obsahují strukturální stabilizaci. Problém regularity by se pak redukoval na formální důkaz, že izolovaný člen skutečně dominuje v kritickém režimu. To by vedlo k důkazu regularity standardních NS.

Interpretace B: F jako nový fyzikální člen. Pokud F nelze izolovat z nelinearity standardních NS — pokud axiomy 1–4 vyžadují informaci, kterou B(v,v) neobsahuje — pak je F rozšíření NS o novou fyziku. Blowup ve standardních NS by pak byl reálný, ale fyzikálně nerelevantní, protože reálná tekutina obsahuje mechanismus, který standardní rovnice nezachycují. To by vedlo k protipříkladu (blowup) pro standardní NS, ale zároveň k lépe definovaným rovnicím pro reálné tekutiny.

Obě interpretace jsou produktivní a obě vedou k pokroku. Axiomatický rámec je nezávislý na tom, která interpretace je správná.

Spektrální stabilizační kritérium

V Fourierově prostoru lze formulovat specifické kritérium pro operátor F:

Nechť Ê(k,t) je energetické spektrum a Ω̂(k,t) je enstrofické spektrum. Operátor F splňuje spektrální stabilizační kritérium, pokud existuje vlnové číslo k* závislé na okamžitém stavu pole takové, že:

pro k < k*: příspěvek F ke spektru je zanedbatelný (nenarušuje inerciální rozsah),
pro k > k*: příspěvek F tlumí enstrofické spektrum rychleji než ν k² (aktivní stabilizace).

Hodnota k* není fixní — závisí na lokální konfiguraci pole. V laminárním proudění k* → ∞ (operátor je neaktivní). V oblasti intenzivního vortex stretchingu k* klesá a operátor se aktivuje na menších škálách.

Tohle je kvalitativně jiné než jakýkoliv LES model. LES filtruje na fixní škále. Operátor F filtruje na adaptivní škále definované lokální nestabilitou.

Návrh numerického testu

Axiomatický rámec dává konkrétní numericky testovatelnou predikci:

Test: DNS simulace standardních NS vs modifikovaných NS (se členem F) pro Taylor-Greenův vortex (dobře definované počáteční podmínky s intenzivním vortex stretchingem).

Metriky:

(a) Časový vývoj maximální enstrofie: standardní NS by měly vykazovat neomezený nebo rychle rostoucí trend; modifikované NS by měly vykazovat saturaci.

(b) Spektrální distribuce: v modifikovaných NS by měl být identifikovatelný zlom v enstrofickém spektru kolem adaptivního k*.

(c) Topologie vírovitých struktur: operátor F by měl generovat pozorovatelné odlišnosti ve struktuře vírovitých trubic — specificky potlačení extrémního natahování bez eliminace vírovitosti jako takové.

Referenční data: Buaria, Pumir & Bodenschatz (2020) poskytují DNS data pro enstrofickou produkci při Re_λ = 140–1300. Srovnání s těmito daty by ukázalo, zda modifikovaný model lépe odpovídá pozorované 'self-attenuation of extreme events.'

Tento test jsme neprovedli. Formulujeme ho jako otevřený experimentální návrh.

Co tvrdíme a co ne

Tvrdíme: Definovali jsme třídu strukturálně stabilizačních operátorů čtyřmi axiomy (divergenční kompatibilita, energetická neutralita, subkritické škálování, dominance v kritickém režimu). Tyto axiomy přímo reagují na Taovu supercriticality barrier a na Lawsonův samoregularizační efekt. Formulovali jsme modifikovanou enstrofickou nerovnost s korekčním členem a identifikovali klíčovou otevřenou otázku: zda korekční člen dominuje nad vortex stretchingem v kritickém režimu. Navrhli jsme konkrétní numerický test.

Netvrdíme: Nemáme důkaz regularity — ani standardních, ani modifikovaných NS. Netvrdíme, že axiomy 1–4 jsou splnitelné — to je samo o sobě otevřená otázka. Netvrdíme, že operátor ακ ze svitku 007 splňuje všechny čtyři axiomy — ověření vyžaduje formální analýzu. Netvrdíme, že toto řeší Millennium Prize.

Co nabízíme: Matematicky formulovanou strategii, jak přemýšlet o problému regularity z pohledu strukturální dynamiky pole. Ne odpověď, ale směr.

Otevřené otázky

Je třída operátorů definovaná axiomy 1–4 neprázdná? Existuje alespoň jeden konkrétní operátor, který všechny čtyři axiomy splňuje?

Pokud ano — je modifikovaná enstrofická nerovnost uzavřitelná? Lze ukázat, že korekční člen C_κ∫Φ(ω) dV skutečně dominuje nad ∫ω·S·ω dV v kritickém režimu?

Jaký je vztah k Cafarelli-Kohn-Nirenberg? Jejich parciální regularita ukazuje, že singulární množina standardních NS má parabolický Hausdorffův rozměr nejvýše 1. Operátor F by měl tuto množinu redukovat na prázdnou — ale jak to formálně ukázat?

Lze axiomy 1–4 odvodit z prvních principů? Z kinetické teorie (Boltzmannova rovnice), z molekulární dynamiky, nebo z nerovnovážné statistické mechaniky?

Co se stane v limitě ν → 0 (Eulerovy rovnice)? Pokud F ≠ 0, mají modifikované Eulerovy rovnice globálně hladká řešení? Pokud ano, je to silnější výsledek než regularita NS.

Kde se potkáváme s konvenční fyzikou

Zpětná vazba v dynamických systémech je známá z teorie regulace (Ljapunov) a samoorganizace (Prigogine). Navier-Stokesovy rovnice ji ale neobsahují. Standardní komunita hledá důkaz regularity pro existující rovnice. Co když je potřeba nejdřív opravit rovnici a pak hledat důkaz?