無(wú)界區(qū)域上的流體力學(xué)數(shù)值模擬存在著本質(zhì)上的困難，近日由中國(guó)科學(xué)院力學(xué)研究所國(guó)家微重力實(shí)驗(yàn)室和美國(guó)達(dá)特茅斯學(xué)院物理與天文學(xué)系組成的聯(lián)合研究團(tuán)隊(duì)采用Hermite插值基函數(shù)，利用偽譜方法進(jìn)行了二維無(wú)界區(qū)域上的高精度流體力學(xué)數(shù)值模擬，相關(guān)成果在線發(fā)表于Journal of Computational Physics(Yin, JCP,2014)。

　　出現(xiàn)在科學(xué)與工程應(yīng)用中的大量問(wèn)題都可以歸結(jié)為無(wú)界區(qū)域上的偏微分方程的數(shù)值求解的問(wèn)題。如帶無(wú)限大彈性地基壩體的應(yīng)力分析，波（包括聲波、電磁波、彈性波、密度波等）在空間中的傳播都是典型的例子。在流體力學(xué)的研究中，許多問(wèn)題也都是在無(wú)界區(qū)域背景下進(jìn)行研究的，如液滴熱毛細(xì)遷移的研究中假設(shè)壁面無(wú)窮遠(yuǎn)，機(jī)翼的繞流問(wèn)題中假設(shè)的無(wú)窮遠(yuǎn)處的來(lái)流，大氣海洋數(shù)值模擬因計(jì)算規(guī)模的巨大只能采用無(wú)界區(qū)域假設(shè)等等。以往的無(wú)界數(shù)值模擬一般采用有界區(qū)域到無(wú)界區(qū)域的投影或人工邊界等方法，而直接對(duì)無(wú)界區(qū)域的模擬進(jìn)行的較少。

　　研究人員采用基底為無(wú)界區(qū)域的Hermite插值基函數(shù)對(duì)流體力學(xué)問(wèn)題進(jìn)行模擬，采用了具有指數(shù)精度的譜方法對(duì)控制方程進(jìn)行展開(kāi)。研究發(fā)現(xiàn)此類方法對(duì)時(shí)間方向上的離散要求較低，只需采用顯式離散格式就能達(dá)到其他方法隱式或半隱式格式所采用的較大時(shí)間步。對(duì)比采用簡(jiǎn)單擴(kuò)大區(qū)域方法進(jìn)行的數(shù)值模擬，新方法在計(jì)算網(wǎng)格較小的區(qū)域內(nèi)就能達(dá)到傳統(tǒng)方法在近百倍計(jì)算區(qū)域內(nèi)達(dá)到的效果。該研究對(duì)其它無(wú)界區(qū)域的物理問(wèn)題求解也將提供新的研究思路。

　　上述研究獲得了國(guó)家自然科學(xué)基金支持，由力學(xué)所國(guó)家微重力實(shí)驗(yàn)室的尹兆華副研究員和美國(guó)達(dá)特茅斯學(xué)院物理與天文學(xué)系的David C.Montgomery教授合作主持完成。

Hermite譜方法模擬（最后一行）和擴(kuò)大區(qū)域算法的Fourier譜方法模擬（前四行）的對(duì)比