導(dǎo)語

偏微分方程（PDEs）是對自然物理定律最普適且簡潔的數(shù)學(xué)描述，能以緊湊的符號化形式捕捉豐富的多尺度物理現(xiàn)象。本文探討機器學(xué)習(xí)推動的PDE研究新范式，包括：（1）發(fā)現(xiàn)新PDE方程與粗粒度近似——從復(fù)雜自然與工程系統(tǒng)中挖掘隱含的控制方程；（2）學(xué)習(xí)有效坐標(biāo)系與降階模型——通過坐標(biāo)變換與維度壓縮提升PDE可分析性；（3）求解算子表示與傳統(tǒng)算法改進——用神經(jīng)網(wǎng)絡(luò)表示解算子并優(yōu)化數(shù)值方法。針對每個方向，我們總結(jié)了關(guān)鍵進展、現(xiàn)存挑戰(zhàn)與未來機遇。

關(guān)鍵詞：偏微分方程（PDEs）、機器學(xué)習(xí)、方程發(fā)現(xiàn)、降階模型、算子學(xué)習(xí)、數(shù)值加速、物理約束

計算模擬大講堂丨來源

論文題目：Promising directions of machine learning for partial differential equations 論文來源：https://www.nature.com/articles/s43588-024-00643-2

簡介

偏微分方程（PDEs）是描述自然規(guī)律與工程系統(tǒng)的核心數(shù)學(xué)工具，從流體力學(xué)的納維 - 斯托克斯方程到量子力學(xué)的薛定諤方程，從氣候模擬的熱傳導(dǎo)方程到生物領(lǐng)域的反應(yīng)擴散方程，PDEs 支撐著多學(xué)科的基礎(chǔ)研究與工程應(yīng)用。然而，傳統(tǒng) PDEs 研究長期面臨兩大瓶頸：一是強非線性、多尺度 PDEs 難以推導(dǎo)解析解，二是有限元、有限差分等數(shù)值方法計算成本極高，面對高維或復(fù)雜邊界問題時常常 “算力告急”。

由美國華盛頓大學(xué) Steven L. Brunton 與 J. Nathan Kutz 團隊發(fā)表于Nature Computational Science（2024 年 7 月，Volume 4，DOI: 10.1038/s43588-024-00643-2）的綜述論文《Promising directions of machine learning for partial differential equations》，系統(tǒng)性梳理了機器學(xué)習(xí)（ML）為 PDEs 研究帶來的突破。論文核心觀點明確：機器學(xué)習(xí)并非替代傳統(tǒng) PDEs 方法，而是通過 “數(shù)據(jù)驅(qū)動 + 物理約束” 的深度融合，從 “發(fā)現(xiàn)新 PDEs 與粗?；]合模型”“學(xué)習(xí)有效坐標(biāo)與降階模型”“優(yōu)化數(shù)值解法與解算子學(xué)習(xí)” 三大方向，解決傳統(tǒng)方法難以攻克的復(fù)雜難題。

研究基于流體、等離子體、神經(jīng)科學(xué)等跨領(lǐng)域數(shù)據(jù)驗證：ML 驅(qū)動的 PDEs 方法可將計算效率提升 1-2 個數(shù)量級（如湍流模擬速度加快 86 倍），同時在小數(shù)據(jù)、高噪聲場景下保持高精度，為 PDEs 研究開辟 “物理機理 + 數(shù)據(jù)智能” 雙驅(qū)動的新范式。

圖 1：機器學(xué)習(xí)推進 PDEs 研究的三大核心方向總覽

該圖以 “目標(biāo) - 方法 - 場景” 三維框架，清晰劃分 ML 在 PDEs 領(lǐng)域的三大應(yīng)用方向，是整篇論文的核心框架圖：

（a）方程發(fā)現(xiàn)（Automated equation discovery）：通過符號回歸技術(shù)，從神經(jīng)科學(xué)時空數(shù)據(jù)（如腦電信號的 u (t,x)）中學(xué)習(xí)可解釋的 PDEs 與粗?；枋觯鉀Q傳統(tǒng)手動推導(dǎo)無法覆蓋的復(fù)雜系統(tǒng)（如神經(jīng)元群體動力學(xué)）；

（b）坐標(biāo)學(xué)習(xí)（Uncovering coordinates）：利用自編碼器網(wǎng)絡(luò)，從圓柱繞流等高維流場仿真 / 實驗數(shù)據(jù)中學(xué)習(xí) “有效坐標(biāo)”，將高維數(shù)據(jù)壓縮至低維 latent 空間，該空間內(nèi)動力學(xué)演化更簡潔（如流場核心運動規(guī)律）；

（c）加速計算（Accelerating computations）：通過算子學(xué)習(xí)實現(xiàn)粗網(wǎng)格數(shù)據(jù)到高分辨率流場的超分辨預(yù)測（上圖流體流場對比），同時融合不同保真度模型（低成本低精度與高成本高精度），構(gòu)建 “低成本 + 高精度” 的混合模型（下圖成本 - 精度曲線）。

解讀：該圖直觀印證了論文 “ML 針對性解決 PDEs 三大痛點” 的核心邏輯 —— 針對 “未知系統(tǒng)無方程”“高維數(shù)據(jù)難分析”“數(shù)值計算效率低”，分別提供數(shù)據(jù)驅(qū)動發(fā)現(xiàn)、低維表征、高效計算的解決方案，為后續(xù)細(xì)分方向奠定理論框架。

圖 2：稀疏回歸從數(shù)據(jù)中發(fā)現(xiàn)納維 - 斯托克斯方程。以 2D 圓柱繞流為案例，完整展示 “從數(shù)據(jù)到 PDEs” 的發(fā)現(xiàn)流程，是 “方程發(fā)現(xiàn)” 方向的關(guān)鍵驗證：

（a）數(shù)據(jù)采集：獲取流場的 x/y 方向速度場（u、v）與渦量（ω，由速度場旋度計算）的快照數(shù)據(jù)；

（b）候選項構(gòu)建：計算數(shù)值偏導(dǎo)數(shù)（如 ω?、ω?、ω??），將狀態(tài)變量、偏導(dǎo)數(shù)及非線性組合（如 uω?、vω?）整合為候選項矩陣 Θ；

（c）稀疏回歸：通過帶 L0 正則的稀疏回歸（目標(biāo)函數(shù)為最小化預(yù)測誤差與系數(shù)稀疏性），篩選出非零系數(shù) ξ，剔除冗余候選項；

（d）PDEs 合成：最終得到渦量輸運方程 ω? + 0.9931uω? + 0.9910vω? = 0.0099ω?? + 0.0099ω??，其中雷諾數(shù)（Re）識別誤差僅 1%，與理論納維 - 斯托克斯方程高度一致。

解讀：證明 “數(shù)據(jù)驅(qū)動發(fā)現(xiàn) PDEs” 的可行性 —— 無需依賴第一性原理假設(shè)，僅通過數(shù)據(jù)與稀疏回歸即可還原經(jīng)典 PDEs，且保留物理可解釋性（如對流項、擴散項的系數(shù)符合理論預(yù)期）。這為神經(jīng)科學(xué)、流行病學(xué)等 “缺乏成熟理論模型” 的領(lǐng)域提供了新研究工具。

圖 3：稀疏貝葉斯方法構(gòu)建海洋模型的 LES 閉合項

聚焦 “粗粒化閉合模型” 方向，以大尺度海洋中尺度模擬為場景，展示 ML 解決多尺度 PDEs 難題的能力：

（a）輸入數(shù)據(jù)：高分辨率海洋流場模擬的速度場（u、v），從中診斷亞網(wǎng)格渦動量強迫項 S（傳統(tǒng)大渦模擬 LES 難以直接計算的核心量）；

（b）候選函數(shù)庫：基于流體力學(xué)機理，構(gòu)建含散度（σ=??u）、渦量（ζ=?×u）、剪切形變（D=?u/?y+?v/?x）等物理意義明確的候選函數(shù)庫；

（c）稀疏貝葉斯學(xué)習(xí)：采用相關(guān)向量機（RVM）進行迭代稀疏回歸，逐步剔除冗余函數(shù)，最終得到 S 的解析表達(dá)式（基于 resolved 變量的線性組合）；

（d）性能驗證：學(xué)習(xí)得到的閉合模型預(yù)測值（??）與真實值（S?）的均值、標(biāo)準(zhǔn)差高度吻合，且精度優(yōu)于深度學(xué)習(xí)黑箱模型。

解讀：解決 PDEs “多尺度耦合” 的傳統(tǒng)痛點 —— 大尺度系統(tǒng)（如海洋、氣候）的 PDEs 需考慮亞網(wǎng)格過程（如小尺度渦旋），傳統(tǒng)方法需手動設(shè)計閉合項，而 ML 可自動從數(shù)據(jù)中學(xué)習(xí)物理可解釋的閉合模型，且適配性更強（如不同海域的流場特征）。

圖 4：神經(jīng)網(wǎng)絡(luò)學(xué)習(xí) Koopman 線性化坐標(biāo)變換。以非線性伯格斯方程（Burgers’ equation）為對象，展示 “坐標(biāo)學(xué)習(xí)” 方向的核心技術(shù) —— 通過 ML 實現(xiàn)非線性 PDEs 的線性化：

（a）網(wǎng)絡(luò)架構(gòu)：包含編碼器 φ?、線性動力學(xué)算子 K、解碼器 φ_d 三部分；輸入為原始狀態(tài) u（服從非線性伯格斯方程），編碼器將 u 映射至低維 latent 空間變量 v；

（b）線性化過程：在 v 空間中，動力學(xué)演化呈線性（v???=Kv?，K 為線性算子），規(guī)避了原始空間的非線性復(fù)雜度；解碼器再將 v 映射回 u 空間，得到下一時刻的狀態(tài) u???；

（c）物理意義：該網(wǎng)絡(luò)自動學(xué)習(xí)到 Cole–Hopf 變換，成功將非線性伯格斯方程轉(zhuǎn)化為線性熱方程，實現(xiàn) “非線性問題線性化求解” 的突破。

解讀：突破 “非線性 PDEs 難分析” 的傳統(tǒng)瓶頸 ——Koopman 理論雖能將非線性系統(tǒng)線性化，但傳統(tǒng)方法難以找到坐標(biāo)變換；ML 通過神經(jīng)網(wǎng)絡(luò)自動學(xué)習(xí)變換關(guān)系，使復(fù)雜系統(tǒng)可通過成熟的線性工具（如譜分析）分析，大幅降低計算與理論推導(dǎo)成本。

圖 5：機器學(xué)習(xí)插值加速高分辨率流場計算。聚焦 “加速計算” 方向，以 2D 柯爾莫哥洛夫流為案例，驗證 ML 提升 PDEs 數(shù)值計算效率的能力：

（a）效率 - 精度對比：左側(cè)圖展示不同網(wǎng)格分辨率下的計算時間與精度（相關(guān)性 > 0.95 的時間），ML 插值方法（橙色虛線）比直接高分辨率模擬（藍(lán)色實線）快 86 倍，且精度損失極??；

（b）泛化性測試：右側(cè)上圖對比訓(xùn)練數(shù)據(jù)（強迫湍流）與驗證數(shù)據(jù)（更大域、衰減湍流、更高湍流強度）的流場，ML 模型在未見過的場景中仍保持高精度；

（c）模型架構(gòu)：右側(cè)下圖為 CNN-based 架構(gòu)，輸入為當(dāng)前時刻速度場 v (t)、外強迫 F (t)、對流通量 φ??，通過卷積層與插值模塊，輸出下一時刻高分辨率速度場 v (t+Δt)，實現(xiàn)端到端的粗→高分辨率映射。

解讀：印證 “ML 顛覆性提升計算效率” 的核心結(jié)論 —— 傳統(tǒng)高分辨率流場模擬需依賴超算（如 8192×8192 網(wǎng)格），而 ML 通過數(shù)據(jù)驅(qū)動的插值策略，在普通計算設(shè)備上即可實現(xiàn)高精度結(jié)果，為工程設(shè)計（如飛行器氣動模擬、氣候快速預(yù)測）降本提效。

線性代數(shù)：一名合格科研人的筑基課

在科研世界中，無論你研究的是人工智能、生物信息、網(wǎng)絡(luò)科學(xué)，還是物理與工程，幾乎所有復(fù)雜系統(tǒng)的建模與推理都指向同一種底層語言——線性代數(shù)。它不僅是計算公式的集合，更是一名科研人理解“結(jié)構(gòu)”、刻畫“變換”、判斷“穩(wěn)定性”、提取“信息”的基本思維框架。本課程以系統(tǒng)科學(xué)的視角重新解構(gòu)線性代數(shù)，帶你越過技巧、直達(dá)本質(zhì)，在跨學(xué)科的真實問題中建立起科研必備的數(shù)學(xué)基石。

集智學(xué)園聯(lián)合清華大學(xué)數(shù)學(xué)博士諸葛昌靖老師開設(shè)，課程將于12月20日開啟，現(xiàn)在加入可享早鳥價格。

詳情請見：

1.

2.

3.

4.

5.

6.