是時(shí)候重新給業(yè)界帶來(lái)一些小小的震撼了……
……
在剛才瑪麗安·米爾扎哈尼發(fā)來(lái)的郵件中，提到了“反向利用人類證明龐加萊猜想的思路”。
但實(shí)際上，如果非要說(shuō)的話，水平集方法更貼近于反向的流形學(xué)習(xí)算法。
當(dāng)然，實(shí)際上并非如此，只是思路上可以這么概括罷了。
流形學(xué)習(xí)算法是把高維數(shù)據(jù)處理到低維，使人類能夠更容易理解。
而水平集方法則是把低維數(shù)據(jù)投射到高維，以便于計(jì)算機(jī)進(jìn)行運(yùn)算。
在數(shù)值計(jì)算領(lǐng)域，這算是一類在80年代末才被提出來(lái)的“新”算法，但因?yàn)閼?yīng)用范圍包括但不限于液體霧化、蒸發(fā)、燃燒、表面材料計(jì)算、圖像識(shí)別……總之跟流形學(xué)習(xí)一樣幾乎包羅萬(wàn)象，所以很快被推廣到了各個(gè)領(lǐng)域。
包括torch
ultiphysics，也在正式版當(dāng)中引入了幾種典型的水平集方法。
但問(wèn)題在于，目前的水平集方法，即便經(jīng)過(guò)改進(jìn)，也只能做到在特定區(qū)間內(nèi)“接近”守恒。
這是個(gè)很麻煩的問(wèn)題。
尤其是在它最主要的應(yīng)用領(lǐng)域——兩相流計(jì)算當(dāng)中更是如此。
因?yàn)樵诖蠖鄶?shù)情況下，兩相流問(wèn)題中涉及到的兩種流體都可以視為互不相溶的，這就意味著所采用的方法必須有效保持各相的質(zhì)量分別守恒。
所以，目前的水平集算法在長(zhǎng)軸距時(shí)間計(jì)算中，很容易出現(xiàn)劇烈震蕩乃至于發(fā)散——
而這本來(lái)應(yīng)該是水平集算法最擅長(zhǎng)的部分。
比如，在之前在設(shè)計(jì)渦扇10發(fā)動(dòng)機(jī)燃燒室的過(guò)程中，常浩南曾將流體體積方法和水平集方法結(jié)合，對(duì)航空燃油離開(kāi)霧化噴嘴之后的彌散行為進(jìn)行了計(jì)算。
結(jié)果么……
有用，但不多。
相比于壓氣機(jī)和渦輪部分相對(duì)“簡(jiǎn)單”的單相流體計(jì)算，項(xiàng)目團(tuán)隊(duì)在燃燒室的設(shè)計(jì)過(guò)程中獲得的指導(dǎo)最少，被迫采用了大量經(jīng)驗(yàn)結(jié)合中試的老辦法。
這也導(dǎo)致部件級(jí)測(cè)試中超過(guò)60的時(shí)間和經(jīng)費(fèi)都耗費(fèi)在了這上面。
好在，得益于小涵道比航發(fā)燃燒室本身的體積有限，擴(kuò)散過(guò)程的持續(xù)時(shí)間并不算特別長(zhǎng)，所以計(jì)算結(jié)果還不至于震蕩得太離譜，而且渦扇10所追求的性能指標(biāo)相對(duì)于其優(yōu)越的總體設(shè)計(jì)而實(shí)在是比較低，所以最后并沒(méi)有影響到什么。
但如果未來(lái)追求尺寸更大、數(shù)據(jù)更先進(jìn)的型號(hào)，比如9x那樣的龐然大物。
或者是另外一種應(yīng)用場(chǎng)景，火箭發(fā)動(dòng)機(jī)——
無(wú)論液體燃料還是固體燃料，由于火箭發(fā)動(dòng)機(jī)要自帶全部推進(jìn)劑，因此對(duì)燃燒-噴射過(guò)程的依賴程度遠(yuǎn)超航空發(fā)動(dòng)機(jī)。
那么這個(gè)守恒問(wèn)題還是需要解決的。
當(dāng)然，既然這篇論文的標(biāo)題叫做“一種……的方法”，那就說(shuō)明肯定不是從理論上標(biāo)本兼治。
只是在某種特定應(yīng)用場(chǎng)景下適用。
不過(guò)，即便如此，對(duì)于目前的水平集方法來(lái)說(shuō)也是一個(gè)巨大的進(jìn)步了。
……本文將提出一類處理帶發(fā)散自由速度場(chǎng)的二相流問(wèn)題的守恒性水平集方法，并在這一過(guò)程中開(kāi)發(fā)了不規(guī)則界面上rob邊界條件的高效處理方法，結(jié)合界面解析的相變求解方法來(lái)考慮液體形狀及內(nèi)部環(huán)流等對(duì)傳熱傳質(zhì)的影響……
首先假設(shè)任意區(qū)域Ω，以及其子區(qū)域Ω1和Ω2，且有Ω2=ΩΩ1。Γ為分割Ω的界面，構(gòu)造一個(gè)正則化函數(shù)Φ（通常為光滑的海維賽德函數(shù)）隱式表示Γ，使得當(dāng)經(jīng)過(guò)Γ時(shí)，Φ快速地從1變?yōu)?，而在子區(qū)域Ω1內(nèi)，則有Φ≈1，子區(qū)域Ω2內(nèi)，有Φ≈0，故而一般將Γ設(shè)為Φ的05水平集，在計(jì)算中，我們希望即使在有小擾動(dòng)存在的情況也可以保持界面的形狀……
……
時(shí)間，隨著常浩南雙手在鍵盤(pán)上的飛速敲擊而飛速流逝著。
不過(guò)在大約二十分鐘后，常浩南就停下了打字的動(dòng)作。&lt-->>;br>稍加思索之后，他抽出一張軟盤(pán)，把才寫(xiě)了個(gè)開(kāi)頭的文檔復(fù)制進(jìn)去，然后起身徑直前往了機(jī)房。
實(shí)際上，他開(kāi)發(fā)的這種方法在數(shù)學(xué)上并不復(fù)雜，實(shí)質(zhì)就是一種重新初始化過(guò)程，通過(guò)求解一個(gè)粘性項(xiàng)的穩(wěn)定解來(lái)修復(fù)水平集方程中的ψ值，以達(dá)到維持守恒性的目的。
而之所以過(guò)去從未有人涉及到，是因?yàn)槠湓诰唧w的計(jì)算方法上難以實(shí)現(xiàn)。
（請(qǐng)）
是時(shí)候重新給業(yè)界帶來(lái)一些小小的震撼了……
比如傳統(tǒng)的網(wǎng)格劃分方式，就很難適應(yīng)笛卡爾網(wǎng)格上高速變化的兩相流體。
但對(duì)于重生后首先就點(diǎn)出了數(shù)值計(jì)算專精，還全程主導(dǎo)了torch
ultiphysics軟件開(kāi)發(fā)和算法設(shè)計(jì)的常浩南來(lái)說(shuō)，這反而不是什么問(wèn)題。