Re: [閒聊] 千禧年七大數學問題-納維斯托克斯

看板Physics作者 (harry901)時間7年前 (2017/06/26 03:18), 7年前編輯推噓2(2014)
留言16則, 5人參與, 最新討論串4/4 (看更多)
這邊 來討論介紹一下 NS eq在數學上目前的問題 首先把NS eq叫出來 他是長這樣的 → DV → 2→ ρ--- = -▽P + ρg + μ▽ V 慣性力 = 壓力 + 物體力 + 黏性力 ....eq(1) Dt 對於非黏流體 最後一項消去 則變成 → DV → ρ--- = -▽P + ρg + .....eq(2) Dt 這個東西叫做Euler eq. Euler eq. 看上面就知道跟NS eq只差了一個Laplacian 學過偏微分方程的都知道 這個差異可以說天差地別 簡單的PDE 隨便猜個通解丟進去翻雲覆雨 就可以得解 不然就是power series solution暴力 然而只因這樣 人類已經頭痛到唉唉叫了 原始的NS還是多了laplccian這個算子在裡面 更何況Eq(1)是3D的differential close form 在PDE的研究中 1D linear 叫做入門 2D linear叫做專業 至於3D還非linear的這種鬼東西根本就是有病 可想而知為什麼會有困難了 難在NS的表述是Nonlinear 但是去解他變成是一大工程 前文有介紹過NS的重要性 NS相當於 流體力學的F=ma 你只要假設合理 比如 1. 小明參加游泳比賽 初始時 在水中蹬牆給予500N的力量 水的黏滯係數已知 試問小明不出力可以游多遠? 這個是1D問題 丟進去NS eq 輕輕鬆鬆得解 2. 小明是一顆球體 重500N 今天掉入水中 小明的密度已知 水的什麼都知道 試問 小明掉若水中的過程是如何? 這個是3D問題 丟進去NS eq 鄉民算不出來 但是電腦該該叫算得出來 3. 小明是一顆球體 重500N 今天掉入水中 小明的密度已知 水的什麼都知道 試問 小明掉若水中的1s內 水波的形狀是? 這個是3D問題 丟進去NS eq 鄉民算不出來 但是電腦算了一星期然後燒掉 大致上是這樣的概念 因為NS是一個目前公認描速流體受力運動最基本的方程式 雖然計算流體力學(CFD)領域 已經做了很多合理的簡化 就好像在固體力學裡面ansis等現行軟體 鄉民不懂得統統丟進去 他都可以告訴你接下來會發生什麼事情 但是如果再回頭看看eq(1) 然後想想你活的維度是3d + 1t時間維度 那個東西展開來的PDE要去解 根本就是諾貝爾+等級了 CFD是可以解決很多現象 但是回歸到這個方程式 為什麼數學界那麼重視? 如果要討論數學跟物理之間的關係 這部份有好多好多可以討論的 (請其他鄉民補充 就算是wiki上來也沒關係 畢竟這部份很冷門但卻很有趣) 因為這個方程式跟我們生活非常密切 就這樣而已! 一台超級電腦花10年解出一個汽車外型流體問題 試著想想 這樣值得嗎? 一台超級電腦 CPU的設計因為散熱問題 丟給其他CPU去解我要怎麼散熱? (PS. 電腦CPU散熱其實是商業上最重要的根本問題 所以很多專利/專業甚麼的 這東西又回到流體力學/熱力學的問題了) 回到問題原點 數學界為什麼那麼重視?(以下言論請數學人不要生氣 物理人單純直覺) 如果回到歷史 可以發現數學跟物理兩者是相輔相成 物理先於數學而有數學領域的發展 反之 這兩種情況 履見不鮮 因此數學物理本一家 這點我個人是同意的 最明顯的case:牛頓被蘋果打到之後弄了微積分 雖然萊布尼茲喊冤 但這幾個世紀的數學物理語言都是那個長長的S表示積分 就知道數學物理本一體了 牛頓說:我當初就是假設 就是猜測 F 正比於 ma 這種狀況阿 (然後搞得後人要死要活) 隨著科技的發展 電腦的普及 物理界跟數學界 針對F=ma是否ok 經過實驗與驗證 這點不用懷疑 但是如果今天有個東西 物理人說 他是F=(ma)^2 呢? 以目前的人去觀測 1N的力推動1kg質量的物體 加速度? 數學跟物理早有定論 但若以什麼都未知的情況去想像如果發生這種歧見呢? 物理人會說:我實驗給你看 1N的時候去推小明然後怎樣 數學人會說:我證明給你看F=ma跟 F=(ma)^2不是等價的 就在這兩個爭議點 兩者取得平衡 達成共識 所以F=ma是OK的 回到原始ns-eq 小弟有說了 大家可以理解為流體的F=ma 物理人說:即使是這樣 統統丟進去給電腦處理就好拉 數學人說:你如何確認電腦是正確可信的? 當然以上兩句是敷衍了 各領域的研究以及嚴謹度 都是很嚴格的 如果"抽象且具體"的的解釋NS eq目前遇到的困難 物理人:唷 對於宇宙中的流體 統統丟進去電腦就好拉 數學人:幹!你先證明宇宙的流體都合適再來跟我說嘴 ps1. 你先證明地球上流體都ok再來談 ps2. 你先證明地太陽系上流體都ok再來談 物理人:窩~ 你看F=xa嘛 然後我就給他F^T=ma^T什麼的 不就好了嗎? 數學人:幹!那個F沙洨的^T 那不是數學語言阿 你在亂扯殺洨 你他媽的給我回去 重學數學 物理人:你他娘的感批評我的能力!? 去算算看小明到底會不會死再來跟我說 於是 這個爭議點 就在於數學家對於解之嚴格性的吹毛求疵 如同機械系/電機系 對於電阻單位 0.002歐姆就好了 數學系偏要說:幹你怎麼知道那個0.002之後 沒有其他歐姆? 就這樣 數學跟物理 常常在這邊打架 NS eq.是否在這邊打架呢? 請神人回應 (抱歉累了 先休息) -- ※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 114.36.62.7 ※ 文章網址: https://www.ptt.cc/bbs/Physics/M.1498418307.A.042.html ※ 編輯: harry901 (114.36.62.7), 06/26/2017 03:30:55 ※ 編輯: harry901 (114.36.62.7), 06/26/2017 03:36:08 ※ 編輯: harry901 (114.36.62.7), 06/26/2017 03:38:50

06/26 06:43, , 1F
但我想,數學家不是對已經解出來的case有懷疑(即使是電腦
06/26 06:43, 1F

06/26 06:51, , 2F
解的也一樣,解的穩定性夠、誤差不大的話,數學上也是接受
06/26 06:51, 2F

06/26 06:52, , 3F
的)。數學家擔心的是那些目前連數值都不知道應該怎麼跑的
06/26 06:52, 3F

06/26 06:53, , 4F
case,所以希望能有人證明出只要給了起始條件就能解出平滑
06/26 06:53, 4F

06/26 06:55, , 5F
解,最好能找到一套方法把現在解不動的NS eq.們一次解決。
06/26 06:55, 5F

06/26 06:56, , 6F
這樣還可以直接算出各種紊流模型,搞不好連分類定理都有了
06/26 06:56, 6F

06/26 06:56, , 7F
。以上目前都是夢話。
06/26 06:56, 7F

06/26 08:13, , 8F
ansys
06/26 08:13, 8F

06/26 09:48, , 9F
我只覺得流體力學好難...
06/26 09:48, 9F

06/28 02:53, , 10F
算得出來的時候還可以大混過去 算不出來的時候呢?
06/28 02:53, 10F

06/28 02:54, , 11F
有時候一樣的東西參數改一下數值軟體就爆了
06/28 02:54, 11F

06/28 02:55, , 12F
這時候是軟體出bug 還是問題根本ill-posed?
06/28 02:55, 12F

06/28 02:56, , 13F
如果問題是ill-posed 那debug一百年也是白費力氣
06/28 02:56, 13F

06/28 02:57, , 14F
但如果數學上能證明解存在 那至少你可以放心debug
06/28 02:57, 14F

06/28 02:59, , 15F
所以不要嫌數學家龜毛 應該把它當成worst case的保證
06/28 02:59, 15F

07/04 15:38, , 16F
ansys不僅僅在固力強大了我猜
07/04 15:38, 16F
文章代碼(AID): #1PK0o312 (Physics)
文章代碼(AID): #1PK0o312 (Physics)