流體中失效的歐拉方程

2020-12-03 木木聊科學

理想世界的流體方程

1757年,數學家歐拉(Leonhard Euler)發現了後來被稱為「歐拉方程」的流體方程,這些方程描述了流體隨時間的演化,就像牛頓的力學方程描述撞球在桌子上的運動一樣。

歐拉方程是一種理想化的對流體運動的數學描述,它們在一定的假設範圍內,模擬流體的運動。更確切地說,歐拉方程描述了流體中無窮小的粒子的瞬時運動。這個描述包括一個粒子的速度和它的渦量(即旋轉的速度和方向)。總的來說,這些信息匯聚成了一個「速度場」,描繪了流體在給定時刻的運動情況。歐拉方程從一個初始速度場開始,預測它在未來每一刻會發生的變化。

兩個多世紀以來,它們似乎做到了描述任何情況下的任何流體運動。然而多年來,一些數學家一直懷疑歐拉方程在某些特定的情況下會失效,因為歐拉方程並不是對真實世界流體的完全描述,它包括幾個非物理性的假設。例如,它們假設當流體的內流在流過彼此時,不會產生摩擦;再比如它們還假設流體是「不可壓縮的」,這意味著在歐拉方程的世界裡,流體是無法被壓縮到比它已經佔據的空間更小的空間裡的。

由不可壓縮的歐拉方程支配的理想流體的運動:u為流體的速度場,p為內壓強的力。

顯然,真正的流體內部是有摩擦的,歐拉方程的描述將流體的運動規定在了一個特定的理想化世界中。若要模擬更真實的流體運動,則需要使用納維-斯託克斯方程(NS方程)。許多數學家和其他研究人員想要知道,在科學上仍佔有非常崇高地位的歐拉方程,在無摩擦、不可壓縮的理想化世界中,是否總是能夠精確地描述流體的所有未來運動狀態。

終於,一個新的證明找到了會讓歐拉方程失效的特定條件。

尋找歐拉方程的缺陷

今年,加利福尼亞大學聖地牙哥分校(UCSD)的數學家Tarek Elgindi分別於4月和10月向arXiv提交了兩篇論文,這兩篇論文推翻了這組著名流體方程在幾個世紀以來的假設。

Elgindi證明了,在一組特定的情況下,歐拉方程會開始輸出無意義的東西。他所找到的特例令眾數學家們大吃一驚,因為這一「特例」是過去數學家們一直以為總能使方程有效的條件。但這並不代表歐拉方程從此將失去它在科學界的重要地位。

從理論上說,根據歐拉方程的運作原理,當你將當前狀態的值代入方程中之後,就能產生未來某一時刻的精確值;再將未來某一時刻的值代入方程,就能再次延伸預測。通常情況下,這個過程是一直有效的,它似乎能延伸到我們可預見的遙遠未來。但是,一旦在這個過程中歐拉方程開始產生一個無法繼續代入計算的值時,歐拉方程便失效了。

是什麼樣的值無法被繼續代入計算?它們一定是以某種非常不合理的方式放大了流體中的某個點的速度或渦量。這种放大非常極端,會將某一點的速度或渦量在有限的時間內放大到無窮大。一旦出現無窮大,方程就會崩潰,無法再繼續對未來狀態進行描述。

這些致命的極值被稱為「奇點」。因此當數學家在詢問「歐拉方程是否總是成立」時,他們實際上是在問:歐拉方程是否會在某些情況下產生奇點

多數數學家相信答案是肯定的,但他們從來沒有找到一個具體的實例。直到Elgindi的證明出現。他的結果雖然沒有表明歐拉方程會在某個確切條件下產生奇點,但這已經是迄今為止最接近這一目標的結果。

降低複雜度

為了實現這個目標,Elgindi考慮了一個簡化版的流體運動模型。在真實的三維流體中,任何粒子都有三個可以移動的軸,即x軸(左右)、y軸(上下)和z軸(前後),它們有很大的運動自由度;而且流體中的不同部分的粒子的運動不一定有任何密切聯繫。

在Elgindi的研究中,他簡化了歐拉方程需要處理的工作。他讓流體的運動關於z軸對稱,這種對稱在真實流體中雖然並不存在,但卻能使得對速度場的計算更加容易。他還限制了流體的運動範圍,流體中的粒子只可以沿z軸的方向,或朝著或遠離z軸運動,不能繞著z軸旋轉。這樣的設定基本上把問題簡化成一個二維問題。

最後,Elgindi對他輸入歐拉方程中的初始數據設定了一些額外規定。從某種意義上說,這些數據比描述真實世界流體的數值更加粗糙,更有可能形成奇點。

在現實世界裡,如果你從流體中的一個點移動非常小的距離到另一點,那麼第二個點的速度和第一個點的速度應該非常相似,它們的渦量應該也非常相似。具有這種特性的速度場被稱為是「平滑」的,也就是說,當你從一個點移動到下一個點時,速度場的值會連續平滑地變化,而不是快速變化。

但在Elgindi對流體的描述中情況卻不是這樣,他的數據中的渦量變化更大。雖然看起來Elgindi的簡化似乎與現實的流體行為偏離太多,但與許多其他數學家為研究歐拉公式所做過的簡化相比,這已經是非常溫和的設定。最終,Elgindi證明了在這些簡化過的情況下,歐拉方程已經開始產生非常意外的結果。

極值出現

在Elgindi的證明中,他設定的流體是沒有邊界的,就像是在空間漂浮的一個點。現在,我們用水箱中的水為例來理解他的證明。想像在水箱的兩端有兩個厚厚的水環,它們就像漩渦一樣在流體的主體內形成有組織的擾動。這種現象在自然界中確實存在。

現在,假設這兩個環朝著相對的方向移動。在前進的過程中,歐拉方程正常運行,計算出流體在每個時刻的速度場。但當環越靠越近時,方程就開始出現一些異常值。方程計算出的結果顯示,當兩個環越靠越近時,它們就以越來越大的強度相互吸引,導致環的中央被拉長了,看起來更像一對漏鬥。隨著它們的中心越靠越近,它們的速度也越來越快,最終相撞。

觀察相撞時的速度場,就能看到從未在歐拉方程的假設情況下所看到的東西——奇點。Elgindi證明了歐拉方程在相撞的點能計算出無窮大的渦量。

Elgindi的結果完全改變了數學家看待歐拉方程的方式。在此之前,數學家從來沒有證明過,在沒有邊界的情況下,歐拉方程只在短時間內有效,而不是永遠有效。在這場漫長的尋找歐拉方程中的「弱點」的拉鋸戰中,終於有一位數學家作出了突破。

來源:原理

相關焦點

  • 科學網—帶自由邊界歐拉方程的幾何分析與先驗估計
    帶自由邊界歐拉方程的幾何分析與先驗估計
  • 【遊戲流體力學基礎及Unity代碼(四)】用歐拉方程模擬無粘性染料之公式推導
    最後,1式等式左邊動量可以用下面這個式子表示: 單位時間內動量的變化就是 集齊2~5這四個式子,就召喚出了流體力學三大基本守恆方程之一——動量方程。世界萬物都遵循牛頓定律,流體也不例外。不過牛頓第一定律在流體力學中就是連續性方程,對不可壓縮流體而言,上篇我們已經推導過了,也就是速度的散度為0。牛頓第二定律則變成了動量方程。
  • 茶杯中的無窮大
    這項2013年的研究是應用數學家Thomas Y Hou(侯一釗)和Guo Luo(羅果)對流體進行的數值模擬,其研究的核心是由萊昂哈德·歐拉(Leonhard Euler)於1757年提出的歐拉方程。○ 歐拉方程:在流體力學中,它們描述了理想化的不可壓縮流體在沒有外力作用的情況下的運動,u為流體的速度,p為壓強。
  • 流體動力學NS方程的哲學缺陷
    在2014年,我與好幾個學友談論過流體力學中的NS方程。眼下有空,也就把自身近幾年的思考簡述如下。就NS方程的推導及其所反映的客觀現象而言,NS方程是對流體微元在瞬時意義上變形運動的描述。而在方程中的忽略只是表明方程本身的局限性,並不能被看成是正確的(物理上真實的)。因而,總可以通過把尺度變小來把適當尺度的小渦考察進來。得到的好處是應力(或壓力)更接近於局部實測量,但是,由此所付出的代價是在NS方程中不斷的補充新的應力(或壓力)項。方程趨向於無限的複雜化。最後乾脆就是把流體動力學運動看成是一個統計學現象,從而,在哲學上是對其原出發點NS方程的否定。
  • 茶杯中的無窮大 - 中科院物理所
    ○ 歐拉方程:在流體力學中,它們描述了理想化的不可壓縮流體在沒有外力作用的情況下的運動,u為流體的速度,p為壓強。| 公式來源:[1]歐拉方程被用來模擬流體隨時間的演化,比如當我們想要知道,如果一塊石頭被扔進平靜的池塘中,池塘中的水在五秒鐘後會怎樣運動時,就可以利用歐拉方程來找到答案。
  • 流體力學NS方程推導
    自然界中宏觀情況的流體運動畢竟佔據大多數,NS方程限定了自己的適用條件為宏觀運動,採用稍微專業一點難度術語是流體滿足連續介質假設。連續介質假設的意思就是說,我們在流場中隨意取出流體微團,這個流體微團在宏觀上是無窮小的,因此整個流場的物理量可以進行數學上的極限微分積分等運算;同時,這個流體微團在微觀上是無窮大的,微團中包含了大量分子,以至於可以進行分子層面的統計平均,獲得我們通常見到的流場變量。
  • 《歐拉方程及微分方程建模》思路與方法
    一、歐拉方程及其求解方法具有結構的變係數線性微分方程稱之為歐拉方程.令x=eu,則u=lnx,於是有將原歐拉方程中xky(k)全部用上式代入,則可以將原方程轉化為以y為函數,u為自變量的常係數線性微分方程
  • 變分法——歐拉-拉格朗日方程
    然後對這幾個積分的第二項進行分部積分變換:然後根據"變分法基本引理"(參見變分法基本引理)就可以導出歐拉-拉格朗日方程啦:需要注意的是,歐拉方程是泛函極值的必要條件,但不是充分條件,在處理實際泛函極值問題時,一般不去考慮充分條件,而是從實際問題的性質出發,間接地判斷泛函極值的存在性,直接使用歐拉方程求出極值曲線【往期精選】
  • 流體力學野史
    當時我想甩一個伯努利方程給他們看,可惜他們都不認識這個人,所以最後自來水流速不夠,鄉親們日常生活中也就只能忍受特別慢的農夫山泉了。 圖1:伯努利方程  數學小王子可能覺得伯努利方程只用到了初等數學,逼格不高。於是施展乾坤大挪移把伯努利方程從初等數學一下子拉到微積分層面。升級成伯努利方程2.0版本:歐拉公式。。。
  • 初中歷史-倫哈特·歐拉
    歐拉的天才使純數學和應用數學的每一個領域都得到了充實,他的數學物理成果有著無限廣闊的應用領域。  早在上一個世紀,艾薩克·牛頓就提出了力學的基本定律。歐拉特別擅長論證如何把這些定律運用到一些常見的物理現象中。例如,他把牛頓定律運用到流體運動,建立了流體力學方程。同樣他通過認真分析剛體的可能運動並應用牛頓定律建立了一個可以完全確定剛體運動的方程組。當然在實際中沒有物體是完全剛體。
  • 伯努利方程是動量方程還是能量方程
    伯努裡方程在一般流體力學的教科書中都有推導。例如,巴切勒經典的流體力學書(G. K.
  • 尋找「最好」(2)歐拉-拉格朗日方程
    歐拉-拉格朗日方程(Euler -Lagrange equation) 為變分法中的一條重要方程。
  • 納維-斯託克斯方程的來源
    1759年,瑞士數學家萊昂哈德·歐拉( Leonhard euler)推導出了套描述流體流動的方程。用數學的語言來說,流體是連續的非剛體物質。它可能是液體,比如流經船體的水,也可能是氣體,比如管道中流動的空氣。牛頓第二運動定律表明,作用在物體上的力等於物體的質量乘以加速度。歐拉在此基礎上對流體行為進行了研究,得出了一套方程。
  • 詳解伯努利方程的前世今生
    伯努利方程是流體力學最著名的方程之一。這個方程將流場中的壓強變化與速度變化聯繫在一起,因此在流場分析中成為不可或缺的分析工具。丹尼爾·伯努利在1738年發表的Hydrodynamica一書中最早關注到速度與壓強之間的相關關係。
  • 數學第一家族和「伯努利方程」
    1738年,於《流體力學》一書中,丹尼爾根據伯努利原理給出了等價的數學表述「伯努利方程」:p,代表流體中某點的壓強ν,代表該流體中某點的流速g,為重力加速度況且伯努利家族一向以擅長微積分聞名於世,歐拉是一個合格的助理,他自作主張把伯努利方程升級為伯努利方程2.0:再對比一下伯努利方程,有沒有感覺到一下子從中學到了碩士的水平。
  • 數學上最複雜的公式——納維斯託克斯方程,到底困難在哪兒?
    1775年,歐拉大神決定換個口味,去研究了一個與力學相關的數學領域——流體力學。他從最基本的無粘性流體的特性開始,仔細研究了無粘性流體的運動與動量變化的關係,於是寫成一本《流體運動的一般原理》。書裡留下了一個無粘性流體力學領域最重要的基礎方程。這本書也是流體力學的開山之作。
  • 第二章 流體靜力學
    歐拉平衡微分方程把質量力和表面力聯繫到了一起,宏觀上來理解,流體內部的壓力就是由於質量力的存在而引起的。例如,大海中的壓力,就是由於地球的重力作用引起的,並且和水深成正比。8. 重力作用下,靜止流體中個點的單位重量的流體的總勢能是相等的(勢能守恆)。*流體的勢能包括「位勢能」和「壓強勢能」。在流體中取任意兩點位置z1,z2,則有:
  • 用歐拉方程模擬無粘性染料之代碼實現
    然後繼續回顧我們討論的連續方程,也就是不可壓縮流體必須要滿足散度為0的條件,不僅要在時刻n滿足,也要在時刻n+1也滿足,寫成一維形式如下: 2式又被稱作壓力泊松方程(Pressure Possion Equation),或者說是雅可比迭代。泊松方程和第一章介紹的拉普拉斯方程的區別是,前者有源,後者無源。這裡的源就是密度乘以散度。
  • 關於流體的力學知識,珍藏版!
    但是,追蹤流動著的流體中這個質點或哪個質點是很麻煩的,實際上通常並不關心這個或哪個質點的命運,所以歐拉提出另一種方法,稱位歐拉方法。它和力學中慣用的方法不同,它不去考察流體中的某一質點的運動過程,而是研究每個時刻在空間各點流體的速度分布。這一方法現在被廣泛採用,包括我們下面的討論。3.
  • 拉格朗日力學:歐拉-拉格朗日方程的形象原理與描述
    在經典的牛頓物理學中,系統的拉格朗日是總動能減去總勢能,但在量子場論中,這種簡單的關係不再真實,並且每個時間點的拉格朗日方程是所有空間中所有領域的功能。因此我們要關注的不是任何一個特定理論中的拉格朗日方程,但拉格朗日如何用於預測系統的行為,這具有普遍的實踐和哲學意義。