第296章 恭喜你
他拽過一張空白的草稿紙，筆尖落下去的速度快得像是在追趕什麼東西。
先用顧辛流型把飛行包線三組方程的子流形分別描述。
推進子流形、結構子流形、熱控子流形，三個子流形各自獨立，通過辛幾何的約化程序把各自的物理量映射到位形空間的同一個點上。
三個子流形的交集就是耦合面，而這個耦合面在顧辛流型上的原像，恰好是一個拉格朗日子流形。
在辛幾何里，拉格朗日子流形的曲率有一個非常好的性質，那就是它可以被分解為兩個部分，一部分來自子流形本身的內在幾何，而另一部分則來自子流形嵌入外圍空間的方式。
這個分解是辛幾何里一個經典的結果，叫做Darboux-Weinstein分解。
內在部分對應的是每個子系統自己的物理規律，這部分在之前的論文中里已經處理過了。
嵌入部分對應的就是兩個子系統在接口上相互作用的幾何效應，這就是曲率奇異的源頭。
他開始繼續往下推。
顧辛流型上的分層篩法可以把嵌入部分的曲率從整個流形的曲率里單獨拆出來。
拆出來之後，鞍點圓法沿著最速下降曲線去積分，繞開所有振盪項，積分路徑不再被動地被網格束縛，而是主動去選一條誤差最小的路。
這整套操作，和鄒楊在報告裡寫的降維映射是完全一致的思路，只是鄒楊他們用的是工程語言，而肖宿用的是幾何語言。
說穿了就是把時間維度變成一個新的坐標軸，把飛行包線從「時間演化的三維空間」映射成「四維的幾何空間」。
在這個四維空間裡，從亞音速到高超音速那段最難算的區間不再是時間的函數，而是一條靜態的曲線，曲線的曲率就是飛行包線的耦合剛性。
肖宿在草稿紙上寫滿了推導，從顧辛流型的約化程序一路推到鞍點圓法的積分路徑構造，中間還補了兩個引理。寫到最後一個等號畫出來的時候，他停了筆，盯著那個等號右邊的東西看了好一會兒。
那個表達式簡潔得不像話。
耦合剛性區映射到幾何空間後，複雜度從O(eⁿ)降到了O(n³log n)，比鄒楊估算的O(n³)還多了一個對數因子，但那個對數因子不是壞事，它是分層篩法在噪聲里拆主項時產生的額外精度餘量，相當於給後續的工程實現白送了一個誤差緩衝。
肖宿把這張草稿紙放到一邊，重新拿起了另一張。
既然飛行包線的模板有了，那接下來就是把這個模板套回到NS方程上了。
NS方程的法旗方向曲率，本質上就是慣性子流形和粘性子流形在交接面上的嵌入曲率。
把慣性項和粘性項分別用顧辛流型描述，兩個子流形的交接面就是渦拉伸發生的區域，而這個交接面的嵌入曲率就是法旗方向的曲率發散。
一旦思路通了，其他的就不是問題了。
肖宿逐漸開始越寫越快。
他直接用飛行包線推導里的Darboux-Weinstein分解，把NS方程流形上的曲率拆成了內在部分和嵌入部分。
內在部分對應的是純慣性運動和純粘性耗散，這部分在《Burgers方程與一維類比》里已經解決過了。
嵌入部分對應的就是慣性項和粘性項的相互作用，也就是渦拉伸。
然後就輪到分層篩法上場了。
渦拉伸在流場上並不是均勻分布的，它集中在某些特定的區域，這些區域恰好就是嵌入曲率最大的地方。
用分層篩法把這些區域從整個流場上篩出來，再用鞍點圓法沿著嵌入曲率的最速下降方向去積分，就可以把法旗方向的曲率發散精確地拆解成一列收斂的幾何級數。
每一級對應渦拉伸的一個尺度層級，從大到小依次排列，越往後的層級能量越小，收斂速度越快。
這比單純的「處理奇性」要深刻得多。
以前的數學家們面對NS方程的奇性，要麼試圖繞開它，要麼試圖用更高的數值精度去硬扛它，但肖宿的做法不是繞，也不是扛，而是直接把奇性本身變成了解的一部分。
奇性不再是需要被消除的東西，而是渦拉伸在幾何上的指紋。
每一個曲率發散的位置、每一個發散的模式、每一個發散的強度，都在描述渦拉伸是怎麼發生的、在哪裡發生的、以多大的能量發生的。
也就是說，他不僅解決了NS方程的正則性問題，還順便把渦拉伸的完整機制也解出來了。
……
第二天一早，幾乎是天剛亮肖宿就醒了，他的夢裡幾乎都是關於NS方程證明的公式。
桌面上攤著十幾張寫滿的草稿紙，從飛行包線的耦合面曲率到NS方程的法旗方向曲率分解，整條推導鏈條完整地鋪在上面，中間沒有一個斷點。
洗漱完後，肖宿又把最後那篇論文的草稿重新打開，光標還停在第三節那個位置。
他先把十幾張草稿紙依次排開，重新梳理了一下思緒，才開始輸入。
鍵盤聲又快又勻得響了起來，之後幾乎沒停過。
第三節的「法旗方向曲率分析」很快補齊了，後面又加了一節「基於顧辛流型嵌入曲率分解的渦拉伸機制」，再往後是「分層篩法在NS方程奇性結構中的應用」和「鞍點圓法的最優積分路徑構造」。
寫到最後一節「結論與展望」的時候，陽光已經從窗簾縫隙里漏進來了，在地板上畫出一道長長的光帶。
肖宿敲完最後一個句號，把論文從頭到尾掃了一遍。
全文四十三頁，不算長，但每一頁的信息密度都很大。
屏幕上的公式密密麻麻，從顧辛流型的度規張量 gμν在法旗分解下的曲率約束 Ric(X,Y)≤C⋅∥∇ω∥L2Ric(X,Y)≤C⋅∥∇ω∥L2 ，到利用耗散勢梯度推導出的渦量上界估計 ∥ω(t)∥L∞≤C1eC2t∥ω(t)∥L∞≤C1eC2t ，從弗洛爾同調群的邊界算子 ∂:F(L)→F(L)∂:F(L)→F(L) 在渦拉伸環面中的拓撲不變量計算，到鞍點圓法沿最速下降路徑積分時產生的誤差緩衝項 O(n3logn)O(n3logn) ……
每一個符號都嚴絲合縫，沒有一字多餘。
肖宿的目光緩緩掃過這些公式，往日裡冰冷的符號此刻仿佛有了溫度，那些曾經令人望而生畏的湍流奇性，那些在草稿紙上反覆塗改的推導，此刻都化作了一幅清晰而壯麗的幾何圖景。
從此，湍流將不再是混亂無序的代名詞，而是渦拉伸在四維幾何空間裡編織出的精美紋路，每一處曲率的起伏，每一次級數的收斂，都是自然規律最優雅的表達。
他的指尖輕輕拂過鍵盤，眼底褪去了往日的平靜，難得的泛起一絲難以察覺的震撼。
原來，湍流最本真、最極致之的美是這樣的。
原來，所有的複雜與晦澀的背後，是這樣的。
有時候連肖宿也會覺得震驚，那麼複雜混亂的現實世界，竟然可以用這樣的一行行優美的式子展示出來……
最後，他把論文保存，關掉了LaTeX，又輕車熟路的打開了arXiv的投稿頁面。
上傳。
填標題。
填摘要。
選分類。
提交。
屏幕上跳出「Submission Received」的確認頁面時，肖宿才感覺到眼睛有些發酸。
他揉了揉眼，往椅背上一靠，腦子裡那根緊繃的弦終於鬆了一點。
NS方程，徹底結束了。
肖宿靠在椅子上眯了一小會兒，等再睜開眼的時候已經快八點了。
他站起來正準備去倒杯水，門外就傳來了腳步聲，緊接著是肖宇的大嗓門。
「三哥！吃飯了！」
肖宇的聲音隔著門板傳進來，帶著十四五歲少年特有的那股子中氣，又響又脆，像是怕肖宿聽不見似的。
「嗯。」
肖宿應了一聲，換了件乾淨的T恤，開門下樓。
先用顧辛流型把飛行包線三組方程的子流形分別描述。
推進子流形、結構子流形、熱控子流形，三個子流形各自獨立，通過辛幾何的約化程序把各自的物理量映射到位形空間的同一個點上。
三個子流形的交集就是耦合面，而這個耦合面在顧辛流型上的原像，恰好是一個拉格朗日子流形。
在辛幾何里，拉格朗日子流形的曲率有一個非常好的性質，那就是它可以被分解為兩個部分，一部分來自子流形本身的內在幾何，而另一部分則來自子流形嵌入外圍空間的方式。
這個分解是辛幾何里一個經典的結果，叫做Darboux-Weinstein分解。
內在部分對應的是每個子系統自己的物理規律，這部分在之前的論文中里已經處理過了。
嵌入部分對應的就是兩個子系統在接口上相互作用的幾何效應，這就是曲率奇異的源頭。
他開始繼續往下推。
顧辛流型上的分層篩法可以把嵌入部分的曲率從整個流形的曲率里單獨拆出來。
拆出來之後，鞍點圓法沿著最速下降曲線去積分，繞開所有振盪項，積分路徑不再被動地被網格束縛，而是主動去選一條誤差最小的路。
這整套操作，和鄒楊在報告裡寫的降維映射是完全一致的思路，只是鄒楊他們用的是工程語言，而肖宿用的是幾何語言。
說穿了就是把時間維度變成一個新的坐標軸，把飛行包線從「時間演化的三維空間」映射成「四維的幾何空間」。
在這個四維空間裡，從亞音速到高超音速那段最難算的區間不再是時間的函數，而是一條靜態的曲線，曲線的曲率就是飛行包線的耦合剛性。
肖宿在草稿紙上寫滿了推導，從顧辛流型的約化程序一路推到鞍點圓法的積分路徑構造，中間還補了兩個引理。寫到最後一個等號畫出來的時候，他停了筆，盯著那個等號右邊的東西看了好一會兒。
那個表達式簡潔得不像話。
耦合剛性區映射到幾何空間後，複雜度從O(eⁿ)降到了O(n³log n)，比鄒楊估算的O(n³)還多了一個對數因子，但那個對數因子不是壞事，它是分層篩法在噪聲里拆主項時產生的額外精度餘量，相當於給後續的工程實現白送了一個誤差緩衝。
肖宿把這張草稿紙放到一邊，重新拿起了另一張。
既然飛行包線的模板有了，那接下來就是把這個模板套回到NS方程上了。
NS方程的法旗方向曲率，本質上就是慣性子流形和粘性子流形在交接面上的嵌入曲率。
把慣性項和粘性項分別用顧辛流型描述，兩個子流形的交接面就是渦拉伸發生的區域，而這個交接面的嵌入曲率就是法旗方向的曲率發散。
一旦思路通了，其他的就不是問題了。
肖宿逐漸開始越寫越快。
他直接用飛行包線推導里的Darboux-Weinstein分解，把NS方程流形上的曲率拆成了內在部分和嵌入部分。
內在部分對應的是純慣性運動和純粘性耗散，這部分在《Burgers方程與一維類比》里已經解決過了。
嵌入部分對應的就是慣性項和粘性項的相互作用，也就是渦拉伸。
然後就輪到分層篩法上場了。
渦拉伸在流場上並不是均勻分布的，它集中在某些特定的區域，這些區域恰好就是嵌入曲率最大的地方。
用分層篩法把這些區域從整個流場上篩出來，再用鞍點圓法沿著嵌入曲率的最速下降方向去積分，就可以把法旗方向的曲率發散精確地拆解成一列收斂的幾何級數。
每一級對應渦拉伸的一個尺度層級，從大到小依次排列，越往後的層級能量越小，收斂速度越快。
這比單純的「處理奇性」要深刻得多。
以前的數學家們面對NS方程的奇性，要麼試圖繞開它，要麼試圖用更高的數值精度去硬扛它，但肖宿的做法不是繞，也不是扛，而是直接把奇性本身變成了解的一部分。
奇性不再是需要被消除的東西，而是渦拉伸在幾何上的指紋。
每一個曲率發散的位置、每一個發散的模式、每一個發散的強度，都在描述渦拉伸是怎麼發生的、在哪裡發生的、以多大的能量發生的。
也就是說，他不僅解決了NS方程的正則性問題，還順便把渦拉伸的完整機制也解出來了。
……
第二天一早，幾乎是天剛亮肖宿就醒了，他的夢裡幾乎都是關於NS方程證明的公式。
桌面上攤著十幾張寫滿的草稿紙，從飛行包線的耦合面曲率到NS方程的法旗方向曲率分解，整條推導鏈條完整地鋪在上面，中間沒有一個斷點。
洗漱完後，肖宿又把最後那篇論文的草稿重新打開，光標還停在第三節那個位置。
他先把十幾張草稿紙依次排開，重新梳理了一下思緒，才開始輸入。
鍵盤聲又快又勻得響了起來，之後幾乎沒停過。
第三節的「法旗方向曲率分析」很快補齊了，後面又加了一節「基於顧辛流型嵌入曲率分解的渦拉伸機制」，再往後是「分層篩法在NS方程奇性結構中的應用」和「鞍點圓法的最優積分路徑構造」。
寫到最後一節「結論與展望」的時候，陽光已經從窗簾縫隙里漏進來了，在地板上畫出一道長長的光帶。
肖宿敲完最後一個句號，把論文從頭到尾掃了一遍。
全文四十三頁，不算長，但每一頁的信息密度都很大。
屏幕上的公式密密麻麻，從顧辛流型的度規張量 gμν在法旗分解下的曲率約束 Ric(X,Y)≤C⋅∥∇ω∥L2Ric(X,Y)≤C⋅∥∇ω∥L2 ，到利用耗散勢梯度推導出的渦量上界估計 ∥ω(t)∥L∞≤C1eC2t∥ω(t)∥L∞≤C1eC2t ，從弗洛爾同調群的邊界算子 ∂:F(L)→F(L)∂:F(L)→F(L) 在渦拉伸環面中的拓撲不變量計算，到鞍點圓法沿最速下降路徑積分時產生的誤差緩衝項 O(n3logn)O(n3logn) ……
每一個符號都嚴絲合縫，沒有一字多餘。
肖宿的目光緩緩掃過這些公式，往日裡冰冷的符號此刻仿佛有了溫度，那些曾經令人望而生畏的湍流奇性，那些在草稿紙上反覆塗改的推導，此刻都化作了一幅清晰而壯麗的幾何圖景。
從此，湍流將不再是混亂無序的代名詞，而是渦拉伸在四維幾何空間裡編織出的精美紋路，每一處曲率的起伏，每一次級數的收斂，都是自然規律最優雅的表達。
他的指尖輕輕拂過鍵盤，眼底褪去了往日的平靜，難得的泛起一絲難以察覺的震撼。
原來，湍流最本真、最極致之的美是這樣的。
原來，所有的複雜與晦澀的背後，是這樣的。
有時候連肖宿也會覺得震驚，那麼複雜混亂的現實世界，竟然可以用這樣的一行行優美的式子展示出來……
最後，他把論文保存，關掉了LaTeX，又輕車熟路的打開了arXiv的投稿頁面。
上傳。
填標題。
填摘要。
選分類。
提交。
屏幕上跳出「Submission Received」的確認頁面時，肖宿才感覺到眼睛有些發酸。
他揉了揉眼，往椅背上一靠，腦子裡那根緊繃的弦終於鬆了一點。
NS方程，徹底結束了。
肖宿靠在椅子上眯了一小會兒，等再睜開眼的時候已經快八點了。
他站起來正準備去倒杯水，門外就傳來了腳步聲，緊接著是肖宇的大嗓門。
「三哥！吃飯了！」
肖宇的聲音隔著門板傳進來，帶著十四五歲少年特有的那股子中氣，又響又脆，像是怕肖宿聽不見似的。
「嗯。」
肖宿應了一聲，換了件乾淨的T恤，開門下樓。