摘要:\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E爲了對抗這種過度平滑,科技工作者們提出了混合的粒子-網格方法,然而這類方法存在着計算速度慢,計算結果差強人意的缺陷:往往由於粒子的非連續性,流體的湍流運動會在粒子間拉出空隙,最終在視覺上產生噪聲:\u003C\u002Fp\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cdiv class=\"pgc-img\"\u003E\u003Cimg src=\"http:\u002F\u002Fp1.pstatp.com\u002Flarge\u002Fpgc-image\u002F2b5ce83541c54587ace8bfaebc975788\" img_width=\"455\" img_height=\"256\" alt=\"北京電影學院發了一篇滿是數學公式的計算機頂會論文\" inline=\"0\"\u003E\u003Cp class=\"pgc-img-caption\"\u003E\u003C\u002Fp\u003E\u003C\u002Fdiv\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E我們提出的算法既能最大可能地保持流體的湍流細節,又能維持流體場的連續性,有效地解決特效解算中這一老大難問題。\u003C\u002Fp\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cdiv class=\"pgc-img\"\u003E\u003Cimg src=\"http:\u002F\u002Fp1.pstatp.com\u002Flarge\u002Fpgc-image\u002Fa12b70ab46af4afcae2bccbfbdd8d2c9\" img_width=\"455\" img_height=\"256\" alt=\"北京電影學院發了一篇滿是數學公式的計算機頂會論文\" inline=\"0\"\u003E\u003Cp class=\"pgc-img-caption\"\u003E\u003C\u002Fp\u003E\u003C\u002Fdiv\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E這是一個由傳統算法計算得到的煙霧動畫,由於數值粘性的存在,煙霧的形態被過度地模糊了。

"\u003Cdiv\u003E\u003Cdiv class=\"pgc-img\"\u003E\u003Cimg src=\"http:\u002F\u002Fp1.pstatp.com\u002Flarge\u002Fpgc-image\u002Fa32f94fcfa71468b8656dc16d633583f\" img_width=\"640\" img_height=\"86\" alt=\"北京電影學院發了一篇滿是數學公式的計算機頂會論文\" inline=\"0\"\u003E\u003Cp class=\"pgc-img-caption\"\u003E\u003C\u002Fp\u003E\u003C\u002Fdiv\u003E\u003Cp class=\"ql-align-justify\"\u003E來源:知乎專欄\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E本文約\u003Cstrong\u003E2400字,\u003C\u002Fstrong\u003E建議閱讀\u003Cstrong\u003E5分鐘。\u003C\u002Fstrong\u003E\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E這篇文章將介紹視覺特效在現代影視製作中穩定現象的產生方法。\u003C\u002Fp\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cdiv class=\"pgc-img\"\u003E\u003Cimg src=\"http:\u002F\u002Fp1.pstatp.com\u002Flarge\u002Fpgc-image\u002F703c3bb1addc40ac9785d9336b5c44bf\" img_width=\"1080\" img_height=\"594\" alt=\"北京電影學院發了一篇滿是數學公式的計算機頂會論文\" inline=\"0\"\u003E\u003Cp class=\"pgc-img-caption\"\u003E\u003C\u002Fp\u003E\u003C\u002Fdiv\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E大家都知道,視覺特效在現代的影視製作中出現十分頻繁。\u003C\u002Fp\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cdiv class=\"pgc-img\"\u003E\u003Cimg src=\"http:\u002F\u002Fp3.pstatp.com\u002Flarge\u002Fpgc-image\u002Ff6b52660b1654631bb2b4491717dbdc4\" img_width=\"455\" img_height=\"256\" alt=\"北京電影學院發了一篇滿是數學公式的計算機頂會論文\" inline=\"0\"\u003E\u003Cp class=\"pgc-img-caption\"\u003E\u003C\u002Fp\u003E\u003C\u002Fdiv\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E而諸如洪水、煙霧、爆炸等特效計算的背後,實際上是用計算機程序在求解已有百年曆史的“納維-斯托克斯方程”\u003C\u002Fp\u003E\u003Cdiv class=\"pgc-img\"\u003E\u003Cimg src=\"http:\u002F\u002Fp1.pstatp.com\u002Flarge\u002Fpgc-image\u002Fc9fc6d3d42ef4d48a8527996d634a73e\" img_width=\"480\" img_height=\"125\" alt=\"北京電影學院發了一篇滿是數學公式的計算機頂會論文\" inline=\"0\"\u003E\u003Cp class=\"pgc-img-caption\"\u003E\u003C\u002Fp\u003E\u003C\u002Fdiv\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E這個方程,對於做流體動力學的讀者一定不陌生,數十年來科學家們爲了計算機翼升力,已將其研究了百千萬遍。\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E然而基於影視製作的特別需求,影視科技工作者們對這個方程的求解提出了新的需求。\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E我們需要能夠處理更大的時間步長以及不損失精度細節!!這兩者本來是魚與熊掌不可兼得,這就等於是在說:\u003C\u002Fp\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cdiv class=\"pgc-img\"\u003E\u003Cimg src=\"http:\u002F\u002Fp9.pstatp.com\u002Flarge\u002Fpgc-image\u002Ffcfe42572403494884b1cc27ede13b2d\" img_width=\"417\" img_height=\"417\" alt=\"北京電影學院發了一篇滿是數學公式的計算機頂會論文\" inline=\"0\"\u003E\u003Cp class=\"pgc-img-caption\"\u003E\u003C\u002Fp\u003E\u003C\u002Fdiv\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E於是伴隨着特效應用的發展,十數年來,計算機圖形學的科技工作者爲了戰勝這個問題不斷地添磚加瓦,這裏面比如就有大家耳熟能詳的科研學府:斯坦福的科學家就曾經提出過使用渦旋力來增加流體的細節。\u003C\u002Fp\u003E\u003Cdiv class=\"pgc-img\"\u003E\u003Cimg src=\"http:\u002F\u002Fp3.pstatp.com\u002Flarge\u002Fpgc-image\u002F61610b7682aa4653894ac5f6a989ede8\" img_width=\"525\" img_height=\"134\" alt=\"北京電影學院發了一篇滿是數學公式的計算機頂會論文\" inline=\"0\"\u003E\u003Cp class=\"pgc-img-caption\"\u003E\u003C\u002Fp\u003E\u003C\u002Fdiv\u003E\u003Cp class=\"ql-align-justify\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cdiv class=\"pgc-img\"\u003E\u003Cimg src=\"http:\u002F\u002Fp1.pstatp.com\u002Flarge\u002Fpgc-image\u002F5d044279dd804b8daae94232ce82cd23\" img_width=\"720\" img_height=\"212\" alt=\"北京電影學院發了一篇滿是數學公式的計算機頂會論文\" inline=\"0\"\u003E\u003Cp class=\"pgc-img-caption\"\u003E\u003C\u002Fp\u003E\u003C\u002Fdiv\u003E\u003Cp class=\"ql-align-justify\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E他們論文研究的算法,甚至直接和工業光魔合作,使用在當年星球大戰特效鏡頭的製作上。\u003C\u002Fp\u003E\u003Cdiv class=\"pgc-img\"\u003E\u003Cimg src=\"http:\u002F\u002Fp1.pstatp.com\u002Flarge\u002Fpgc-image\u002Ff5baa4afcd0144aebe1ae25065f665d9\" img_width=\"327\" img_height=\"343\" alt=\"北京電影學院發了一篇滿是數學公式的計算機頂會論文\" inline=\"0\"\u003E\u003Cp class=\"pgc-img-caption\"\u003E\u003C\u002Fp\u003E\u003C\u002Fdiv\u003E\u003Cp class=\"ql-align-justify\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E還有很多科研機構也是不斷地在算法上推陳出新,甚至像迪士尼這樣家裏有礦的影視工作室還親自上陣,不斷地革新着美利堅的影視製作技術,帶動起了好萊塢產、研結合的影視工具生態。\u003C\u002Fp\u003E\u003Cdiv class=\"pgc-img\"\u003E\u003Cimg src=\"http:\u002F\u002Fp1.pstatp.com\u002Flarge\u002Fpgc-image\u002F963cde29a0f0485dbc6fce567a0c4234\" img_width=\"720\" img_height=\"240\" alt=\"北京電影學院發了一篇滿是數學公式的計算機頂會論文\" inline=\"0\"\u003E\u003Cp class=\"pgc-img-caption\"\u003E\u003C\u002Fp\u003E\u003C\u002Fdiv\u003E\u003Cp class=\"ql-align-justify\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E而不管是由工業光魔發表的用來補充煙火高精度細節的算法:\u003C\u002Fp\u003E\u003Cdiv class=\"pgc-img\"\u003E\u003Cimg src=\"http:\u002F\u002Fp3.pstatp.com\u002Flarge\u002Fpgc-image\u002F4f8c1ef807e54e4b9553cdab8a1c10e1\" img_width=\"720\" img_height=\"300\" alt=\"北京電影學院發了一篇滿是數學公式的計算機頂會論文\" inline=\"0\"\u003E\u003Cp class=\"pgc-img-caption\"\u003E\u003C\u002Fp\u003E\u003C\u002Fdiv\u003E\u003Cp class=\"ql-align-justify\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cdiv class=\"pgc-img\"\u003E\u003Cimg src=\"http:\u002F\u002Fp3.pstatp.com\u002Flarge\u002Fpgc-image\u002Fefe7cf89c73d4ad58b20b7e02b5bce53\" img_width=\"423\" img_height=\"338\" alt=\"北京電影學院發了一篇滿是數學公式的計算機頂會論文\" inline=\"0\"\u003E\u003Cp class=\"pgc-img-caption\"\u003E\u003C\u002Fp\u003E\u003C\u002Fdiv\u003E\u003Cp class=\"ql-align-justify\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E還是來源於迪士尼研究院-ETH Zurich的小波湍流論文:\u003C\u002Fp\u003E\u003Cdiv class=\"pgc-img\"\u003E\u003Cimg src=\"http:\u002F\u002Fp3.pstatp.com\u002Flarge\u002Fpgc-image\u002Fa2a9594edf7c4db39e76ac958c914802\" img_width=\"577\" img_height=\"292\" alt=\"北京電影學院發了一篇滿是數學公式的計算機頂會論文\" inline=\"0\"\u003E\u003Cp class=\"pgc-img-caption\"\u003E\u003C\u002Fp\u003E\u003C\u002Fdiv\u003E\u003Cp class=\"ql-align-justify\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cdiv class=\"pgc-img\"\u003E\u003Cimg src=\"http:\u002F\u002Fp9.pstatp.com\u002Flarge\u002Fpgc-image\u002F6ef3424752b943cbb67ffe4d7efc7aba\" img_width=\"279\" img_height=\"344\" alt=\"北京電影學院發了一篇滿是數學公式的計算機頂會論文\" inline=\"0\"\u003E\u003Cp class=\"pgc-img-caption\"\u003E\u003C\u002Fp\u003E\u003C\u002Fdiv\u003E\u003Cp class=\"ql-align-justify\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E都是不斷在嘗試攻克一個令人頭疼的問題:如何在使用高效穩定的流體特效仿真算法前提下,儘可能地保證流體運動的細節,這最終落在了對於N-S方程中的對流部分的求解上。\u003C\u002Fp\u003E\u003Cp class=\"ql-align-center\"\u003E \u003C\u002Fp\u003E\u003Cdiv class=\"pgc-img\"\u003E\u003Cimg src=\"http:\u002F\u002Fp1.pstatp.com\u002Flarge\u002Fpgc-image\u002Fca1a0edd42ee46f8af21ab0a7f35aef4\" img_width=\"484\" img_height=\"78\" alt=\"北京電影學院發了一篇滿是數學公式的計算機頂會論文\" inline=\"0\"\u003E\u003Cp class=\"pgc-img-caption\"\u003E\u003C\u002Fp\u003E\u003C\u002Fdiv\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E在圖形學特別是影視特效製作軟件中,爲了追求大時間步長的穩定性,對流方程的求解多采用半拉格朗日方法或其變種。這類方法雖然高效穩定,但有它最大的問題:\u003Cstrong\u003E數值粘性\u003C\u002Fstrong\u003E。\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E所謂數值粘性,顧名思義,就是物理系統中本不存在的粘性,是被人類設計的數值算法在求解方程的過程中帶入系統的,與真實的物理事實相左。我們先領略一下“\u003Cstrong\u003E數值粘性\u003C\u002Fstrong\u003E” 會帶來的問題:\u003C\u002Fp\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cdiv class=\"pgc-img\"\u003E\u003Cimg src=\"http:\u002F\u002Fp9.pstatp.com\u002Flarge\u002Fpgc-image\u002Ffb5ed072eed04718b39ca14fa6924a76\" img_width=\"256\" img_height=\"256\" alt=\"北京電影學院發了一篇滿是數學公式的計算機頂會論文\" inline=\"0\"\u003E\u003Cp class=\"pgc-img-caption\"\u003E\u003C\u002Fp\u003E\u003C\u002Fdiv\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cp class=\"ql-align-center\"\u003E比如我要沿着紅線所示的速度和方向“對流”白色的方塊\u003C\u002Fp\u003E\u003Cdiv class=\"pgc-img\"\u003E\u003Cimg src=\"http:\u002F\u002Fp1.pstatp.com\u002Flarge\u002Fpgc-image\u002F2dca40817ed349dba2682bf4efdbd4c1\" img_width=\"256\" img_height=\"256\" alt=\"北京電影學院發了一篇滿是數學公式的計算機頂會論文\" inline=\"0\"\u003E\u003Cp class=\"pgc-img-caption\"\u003E\u003C\u002Fp\u003E\u003C\u002Fdiv\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cp class=\"ql-align-center\"\u003E傳統的對流算法由於數值粘性會使物質模糊掉\u003C\u002Fp\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cdiv class=\"pgc-img\"\u003E\u003Cimg src=\"http:\u002F\u002Fp1.pstatp.com\u002Flarge\u002Fpgc-image\u002F1b8241bb258d4496ab3122a6171bb1c5\" img_width=\"256\" img_height=\"256\" alt=\"北京電影學院發了一篇滿是數學公式的計算機頂會論文\" inline=\"0\"\u003E\u003Cp class=\"pgc-img-caption\"\u003E\u003C\u002Fp\u003E\u003C\u002Fdiv\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cp class=\"ql-align-center\"\u003E我們的算法盡大可能地避免了數值粘性, 保留了清晰地邊界\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E流體的運動更爲複雜,計算機算法保持清晰(守恆)的能力,就顯得尤爲重要。\u003C\u002Fp\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cdiv class=\"pgc-img\"\u003E\u003Cimg src=\"http:\u002F\u002Fp1.pstatp.com\u002Flarge\u002Fpgc-image\u002Fa12b70ab46af4afcae2bccbfbdd8d2c9\" img_width=\"455\" img_height=\"256\" alt=\"北京電影學院發了一篇滿是數學公式的計算機頂會論文\" inline=\"0\"\u003E\u003Cp class=\"pgc-img-caption\"\u003E\u003C\u002Fp\u003E\u003C\u002Fdiv\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E這是一個由傳統算法計算得到的煙霧動畫,由於數值粘性的存在,煙霧的形態被過度地模糊了。\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E爲了對抗這種過度平滑,科技工作者們提出了混合的粒子-網格方法,然而這類方法存在着計算速度慢,計算結果差強人意的缺陷:往往由於粒子的非連續性,流體的湍流運動會在粒子間拉出空隙,最終在視覺上產生噪聲:\u003C\u002Fp\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cdiv class=\"pgc-img\"\u003E\u003Cimg src=\"http:\u002F\u002Fp1.pstatp.com\u002Flarge\u002Fpgc-image\u002F2b5ce83541c54587ace8bfaebc975788\" img_width=\"455\" img_height=\"256\" alt=\"北京電影學院發了一篇滿是數學公式的計算機頂會論文\" inline=\"0\"\u003E\u003Cp class=\"pgc-img-caption\"\u003E\u003C\u002Fp\u003E\u003C\u002Fdiv\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E我們提出的算法既能最大可能地保持流體的湍流細節,又能維持流體場的連續性,有效地解決特效解算中這一老大難問題。\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E把三個方法的結果並排放在一起對比如下:\u003C\u002Fp\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cdiv class=\"pgc-img\"\u003E\u003Cimg src=\"http:\u002F\u002Fp1.pstatp.com\u002Flarge\u002Fpgc-image\u002F7082508a17c94406b28865b08e72ee7b\" img_width=\"216\" img_height=\"386\" alt=\"北京電影學院發了一篇滿是數學公式的計算機頂會論文\" inline=\"0\"\u003E\u003Cp class=\"pgc-img-caption\"\u003E\u003C\u002Fp\u003E\u003C\u002Fdiv\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cp class=\"ql-align-center\"\u003E過度平滑的傳統方法\u003C\u002Fp\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cdiv class=\"pgc-img\"\u003E\u003Cimg src=\"http:\u002F\u002Fp3.pstatp.com\u002Flarge\u002Fpgc-image\u002F9db3ccfda1e84fd0a45439e90784bd22\" img_width=\"216\" img_height=\"387\" alt=\"北京電影學院發了一篇滿是數學公式的計算機頂會論文\" inline=\"0\"\u003E\u003Cp class=\"pgc-img-caption\"\u003E\u003C\u002Fp\u003E\u003C\u002Fdiv\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cp class=\"ql-align-center\"\u003E計算耗時,且噪聲過強的粒子方法\u003C\u002Fp\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cdiv class=\"pgc-img\"\u003E\u003Cimg src=\"http:\u002F\u002Fp1.pstatp.com\u002Flarge\u002Fpgc-image\u002Fe1a16cce13944784907c4377be1a78b0\" img_width=\"216\" img_height=\"387\" alt=\"北京電影學院發了一篇滿是數學公式的計算機頂會論文\" inline=\"0\"\u003E\u003Cp class=\"pgc-img-caption\"\u003E\u003C\u002Fp\u003E\u003C\u002Fdiv\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cp class=\"ql-align-center\"\u003E計算高效,即不失細節又平滑連續的我們的方法\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E爲了以極高的清晰度來對流體中的物質場,我們給出了以下觀察:\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E傳統的對流算法直接作用於流體空間中的“物理量”,“數值粘性”於是會累加在“物理量”上,並最終造成“物理量”的模糊。\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E如果能構造出一個映射,使得流體場中的每一個空間點在任何時候都能知道自己初始時刻的位置並去初始時刻訪問物理量的信息,我們就能夠避免在對流的過程中累加“數值粘性”。\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E這個哲學像極了愛情:\u003Cstrong\u003E一切若只如初見,是多麼的美好\u003C\u002Fstrong\u003E。\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E這個方法從求解的根本思路上跳出了傳統對流方程求解方式:\u003C\u002Fp\u003E\u003Cul\u003E\u003Cli class=\"ql-align-justify\"\u003E以往的方法對流“物理量”\u003C\u002Fli\u003E\u003Cli class=\"ql-align-justify\"\u003E我們的方法對流“空間”本身\u003C\u002Fli\u003E\u003C\u002Ful\u003E\u003Cp class=\"ql-align-justify\"\u003E爲此,我們定義了一個映射,這個映射的目的是將空間中的點映射回它初始時刻的位置 \u003C\u002Fp\u003E\u003Cdiv class=\"pgc-img\"\u003E\u003Cimg src=\"http:\u002F\u002Fp9.pstatp.com\u002Flarge\u002Fpgc-image\u002Ff76570c22460493e8cb01a44a19588d7\" img_width=\"265\" img_height=\"36\" alt=\"北京電影學院發了一篇滿是數學公式的計算機頂會論文\" inline=\"0\"\u003E\u003Cp class=\"pgc-img-caption\"\u003E\u003C\u002Fp\u003E\u003C\u002Fdiv\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cdiv class=\"pgc-img\"\u003E\u003Cimg src=\"http:\u002F\u002Fp3.pstatp.com\u002Flarge\u002Fpgc-image\u002Fc42be6e17110455489895c7ac4effd2b\" img_width=\"256\" img_height=\"256\" alt=\"北京電影學院發了一篇滿是數學公式的計算機頂會論文\" inline=\"0\"\u003E\u003Cp class=\"pgc-img-caption\"\u003E\u003C\u002Fp\u003E\u003C\u002Fdiv\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cp class=\"ql-align-center\"\u003E因爲流體運動而變化的空間映射的可視化\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E這個初始位置,滿足如下的數學關係,稱之爲後向映射。 \u003C\u002Fp\u003E\u003Cdiv class=\"pgc-img\"\u003E\u003Cimg src=\"http:\u002F\u002Fp9.pstatp.com\u002Flarge\u002Fpgc-image\u002F8aa47729221c4006a9c17176569d5191\" img_width=\"525\" img_height=\"61\" alt=\"北京電影學院發了一篇滿是數學公式的計算機頂會論文\" inline=\"0\"\u003E\u003Cp class=\"pgc-img-caption\"\u003E\u003C\u002Fp\u003E\u003C\u002Fdiv\u003E\u003Cp class=\"ql-align-justify\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E在這個基礎上,我們就有機會將原本的對流方程拆爲空間的對流和物理量的重構兩個部分。\u003C\u002Fp\u003E\u003Cdiv class=\"pgc-img\"\u003E\u003Cimg src=\"http:\u002F\u002Fp3.pstatp.com\u002Flarge\u002Fpgc-image\u002F324c2bc70d3a48229551ced6d4154f50\" img_width=\"720\" img_height=\"223\" alt=\"北京電影學院發了一篇滿是數學公式的計算機頂會論文\" inline=\"0\"\u003E\u003Cp class=\"pgc-img-caption\"\u003E\u003C\u002Fp\u003E\u003C\u002Fdiv\u003E\u003Cp class=\"ql-align-justify\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E爲了便於重構物理量時對於積分項的計算,我們提出了整個算法中最至關重要的部分:前向映射。並給出了前向映射的演化方程:\u003C\u002Fp\u003E\u003Cdiv class=\"pgc-img\"\u003E\u003Cimg src=\"http:\u002F\u002Fp1.pstatp.com\u002Flarge\u002Fpgc-image\u002Ff96a970523af4d8d897e172fc33ecf51\" img_width=\"568\" img_height=\"123\" alt=\"北京電影學院發了一篇滿是數學公式的計算機頂會論文\" inline=\"0\"\u003E\u003Cp class=\"pgc-img-caption\"\u003E\u003C\u002Fp\u003E\u003C\u002Fdiv\u003E\u003Cp class=\"ql-align-justify\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E有了前向映射,積分的過程變爲了演化累加的過程(大大減少了計算量)\u003C\u002Fp\u003E\u003Cdiv class=\"pgc-img\"\u003E\u003Cimg src=\"http:\u002F\u002Fp1.pstatp.com\u002Flarge\u002Fpgc-image\u002Fdcfca4a9de2b4770bc969ae788bc95e4\" img_width=\"598\" img_height=\"127\" alt=\"北京電影學院發了一篇滿是數學公式的計算機頂會論文\" inline=\"0\"\u003E\u003Cp class=\"pgc-img-caption\"\u003E\u003C\u002Fp\u003E\u003C\u002Fdiv\u003E\u003Cp class=\"ql-align-justify\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E好了,小編保證,不會再有數學公式了!\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E我們將我們發明的這樣層級前-後向映射來求解對流方程的形式稱之爲“n層級雙向跡徑映射對流求解法”(Bi-Directional Characteristic Mapping of Convective Quantaties of Level n BiMocq)。\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E雖然不打算在這裏展開討論,這樣的做法使得在保持流體物理量守恆的同時還保證計算極爲高效以及並行友好,大大地區別於混合的粒子-網格算法。我們粗暴實現的GPU代碼可以對比相應的CPU的併發式代碼加速50倍。\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E那麼,BiMocq好處都有啥?\u003C\u002Fp\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cdiv class=\"pgc-img\"\u003E\u003Cimg src=\"http:\u002F\u002Fp1.pstatp.com\u002Flarge\u002Fpgc-image\u002F27d2da9b8b4c4dcdb3814491ac4afb25\" img_width=\"378\" img_height=\"339\" alt=\"北京電影學院發了一篇滿是數學公式的計算機頂會論文\" inline=\"0\"\u003E\u003Cp class=\"pgc-img-caption\"\u003E\u003C\u002Fp\u003E\u003C\u002Fdiv\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E我們首先是對比了一個數學上有比較確定答案的問題,“泰勒渦旋”:兩個初始狀態如下的渦旋,在無粘的環境中,能量守恆的數值積分應該可以將其分開。\u003C\u002Fp\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cdiv class=\"pgc-img\"\u003E\u003Cimg src=\"http:\u002F\u002Fp1.pstatp.com\u002Flarge\u002Fpgc-image\u002Fe7767c3ffbe347edaae7531c56cd60d1\" img_width=\"460\" img_height=\"398\" alt=\"北京電影學院發了一篇滿是數學公式的計算機頂會論文\" inline=\"0\"\u003E\u003Cp class=\"pgc-img-caption\"\u003E\u003C\u002Fp\u003E\u003C\u002Fdiv\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cdiv class=\"pgc-img\"\u003E\u003Cimg src=\"http:\u002F\u002Fp1.pstatp.com\u002Flarge\u002Fpgc-image\u002F6f1f85e58b2448479cc5dcab665a83d7\" img_width=\"460\" img_height=\"398\" alt=\"北京電影學院發了一篇滿是數學公式的計算機頂會論文\" inline=\"0\"\u003E\u003Cp class=\"pgc-img-caption\"\u003E\u003C\u002Fp\u003E\u003C\u002Fdiv\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E對此,我們實現並且對比了幾乎所有主流使用的,先前發表於SIGGRAPH論文的算法,並驗證了我們的算法在能量守恆上的優越性。\u003C\u002Fp\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cdiv class=\"pgc-img\"\u003E\u003Cimg src=\"http:\u002F\u002Fp1.pstatp.com\u002Flarge\u002Fpgc-image\u002F876b631a1d4d4b988ad861d042daf641\" img_width=\"455\" img_height=\"256\" alt=\"北京電影學院發了一篇滿是數學公式的計算機頂會論文\" inline=\"0\"\u003E\u003Cp class=\"pgc-img-caption\"\u003E\u003C\u002Fp\u003E\u003C\u002Fdiv\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cp class=\"ql-align-center\"\u003E對比所有的SemiLagrangian類對流算法,BIMocqn 勝\u003C\u002Fp\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cdiv class=\"pgc-img\"\u003E\u003Cimg src=\"http:\u002F\u002Fp1.pstatp.com\u002Flarge\u002Fpgc-image\u002Fb3d871f1b8124e839de1b0e12bc58e28\" img_width=\"455\" img_height=\"256\" alt=\"北京電影學院發了一篇滿是數學公式的計算機頂會論文\" inline=\"0\"\u003E\u003Cp class=\"pgc-img-caption\"\u003E\u003C\u002Fp\u003E\u003C\u002Fdiv\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cp class=\"ql-align-center\"\u003EBIMocqn特寫\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E以及難以用Eulerian方法解算的leapfrogging vortex:\u003C\u002Fp\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cdiv class=\"pgc-img\"\u003E\u003Cimg src=\"http:\u002F\u002Fp9.pstatp.com\u002Flarge\u002Fpgc-image\u002Fc0280ae701e54dab825ec218e33f155e\" img_width=\"455\" img_height=\"256\" alt=\"北京電影學院發了一篇滿是數學公式的計算機頂會論文\" inline=\"0\"\u003E\u003Cp class=\"pgc-img-caption\"\u003E\u003C\u002Fp\u003E\u003C\u002Fdiv\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E更進一步地,用一個簡單的三維煙霧的例子來說明一下能量和質量守恆的好處:\u003C\u002Fp\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cdiv class=\"pgc-img\"\u003E\u003Cimg src=\"http:\u002F\u002Fp1.pstatp.com\u002Flarge\u002Fpgc-image\u002Fc1aa6e7a7f5249e2aac4e9b3588b6b8d\" img_width=\"540\" img_height=\"305\" alt=\"北京電影學院發了一篇滿是數學公式的計算機頂會論文\" inline=\"0\"\u003E\u003Cp class=\"pgc-img-caption\"\u003E\u003C\u002Fp\u003E\u003C\u002Fdiv\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E該例子中,左邊的是目前特效計算中常用的算法右邊的是我們的新算法。由於數值粘性,以往的算法無法保證煙霧運動時的清晰形狀導致了其在之後的效果計算中損失了煙霧的,湍流細節,那些更能體現煙霧真實感的捲曲細節。\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E大自然是最好的藝術家,我們想做的,僅僅是用程序在計算機中再現它的美。爲此,我們用我們的數值算法在計算機中再現瞭如下場景,渦環相撞實驗。\u003C\u002Fp\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cdiv class=\"pgc-img\"\u003E\u003Cimg src=\"http:\u002F\u002Fp1.pstatp.com\u002Flarge\u002Fpgc-image\u002F53f67a3cd5c44b93874936c460c8ac54\" img_width=\"455\" img_height=\"256\" alt=\"北京電影學院發了一篇滿是數學公式的計算機頂會論文\" inline=\"0\"\u003E\u003Cp class=\"pgc-img-caption\"\u003E\u003C\u002Fp\u003E\u003C\u002Fdiv\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E兩個相對運動的渦環的實拍視頻,在碰撞後會不斷往直徑方向擴大並最終由於任何微小的擾動在邊緣處射出花瓣狀的小渦。\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E這樣的現象以往的數值算法存在數值粘性就根本不可能重現出來。\u003C\u002Fp\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cdiv class=\"pgc-img\"\u003E\u003Cimg src=\"http:\u002F\u002Fp1.pstatp.com\u002Flarge\u002Fpgc-image\u002Fd72382aba6d44546b15b1690a5252b66\" img_width=\"238\" img_height=\"203\" alt=\"北京電影學院發了一篇滿是數學公式的計算機頂會論文\" inline=\"0\"\u003E\u003Cp class=\"pgc-img-caption\"\u003E\u003C\u002Fp\u003E\u003C\u002Fdiv\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E我們的算法成功地在計算機中再現了這個現象應該也是該領域第一次用高效的純歐拉網格算法完整地再現了這個過程。\u003C\u002Fp\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cdiv class=\"pgc-img\"\u003E\u003Cimg src=\"http:\u002F\u002Fp3.pstatp.com\u002Flarge\u002Fpgc-image\u002F061256882b10423cb0dddb99c76436ff\" img_width=\"383\" img_height=\"187\" alt=\"北京電影學院發了一篇滿是數學公式的計算機頂會論文\" inline=\"0\"\u003E\u003Cp class=\"pgc-img-caption\"\u003E\u003C\u002Fp\u003E\u003C\u002Fdiv\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E不知道有沒有細心的朋友注意到過這樣一種自然現象。\u003C\u002Fp\u003E\u003Cdiv class=\"pgc-img\"\u003E\u003Cimg src=\"http:\u002F\u002Fp3.pstatp.com\u002Flarge\u002Fpgc-image\u002F541b84e254764f0d9e76332f5ecaa80b\" img_width=\"526\" img_height=\"355\" alt=\"北京電影學院發了一篇滿是數學公式的計算機頂會論文\" inline=\"0\"\u003E\u003Cp class=\"pgc-img-caption\"\u003E\u003C\u002Fp\u003E\u003C\u002Fdiv\u003E\u003Cp class=\"ql-align-justify\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cdiv class=\"pgc-img\"\u003E\u003Cimg src=\"http:\u002F\u002Fp3.pstatp.com\u002Flarge\u002Fpgc-image\u002Ff92a0c466d5e458c84788fbb431ee0fa\" img_width=\"533\" img_height=\"484\" alt=\"北京電影學院發了一篇滿是數學公式的計算機頂會論文\" inline=\"0\"\u003E\u003Cp class=\"pgc-img-caption\"\u003E\u003C\u002Fp\u003E\u003C\u002Fdiv\u003E\u003Cp class=\"ql-align-justify\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E這種現象被稱爲開爾文-亥姆霍茲-瑞利-泰勒不穩定性(Kelvin–Helmholtz-Rayleigh-Taylor instability),這種現象通常由於兩種不同密度的流體在交界面處由於擾動而產生的不穩定現象。\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E同樣,我們使用我們的算法對這個現象進行了近似的模擬。我們可以看到由於密度本身的對流會對整個流場產生影響,對比以往的數值算法,在同等分辨率下由於我們的算法對密度對流有着更準確的估計,因此我們也能更準確地模擬這種不穩定性帶來的渦結構。\u003C\u002Fp\u003E\u003Cdiv class=\"pgc-img\"\u003E\u003Cimg src=\"http:\u002F\u002Fp1.pstatp.com\u002Flarge\u002Fpgc-image\u002F80af56573cd94a609777b04d82ad5474\" img_width=\"488\" img_height=\"335\" alt=\"北京電影學院發了一篇滿是數學公式的計算機頂會論文\" inline=\"0\"\u003E\u003Cp class=\"pgc-img-caption\"\u003E\u003C\u002Fp\u003E\u003C\u002Fdiv\u003E\u003Cp class=\"ql-align-justify\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E作爲展示這個算法在更通用的特效環境中的使用我們也增加了使用我們提出的新算法生成的爆炸以及汽車漂移特效的案例。\u003C\u002Fp\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cdiv class=\"pgc-img\"\u003E\u003Cimg src=\"http:\u002F\u002Fp1.pstatp.com\u002Flarge\u002Fpgc-image\u002F7021c61f18b443abb7c91a856ce857b8\" img_width=\"455\" img_height=\"256\" alt=\"北京電影學院發了一篇滿是數學公式的計算機頂會論文\" inline=\"0\"\u003E\u003Cp class=\"pgc-img-caption\"\u003E\u003C\u002Fp\u003E\u003C\u002Fdiv\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E\u003Cbr\u003E\u003C\u002Fp\u003E\u003Cdiv class=\"pgc-img\"\u003E\u003Cimg src=\"http:\u002F\u002Fp1.pstatp.com\u002Flarge\u002Fpgc-image\u002Fdf15325cf2f947989ec9b718cedc5689\" img_width=\"600\" img_height=\"278\" alt=\"北京電影學院發了一篇滿是數學公式的計算機頂會論文\" inline=\"0\"\u003E\u003Cp class=\"pgc-img-caption\"\u003E\u003C\u002Fp\u003E\u003C\u002Fdiv\u003E\u003Cp class=\"ql-align-justify\"\u003E原文鏈接:\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003Ehttps:\u002F\u002Fzhuanlan.zhihu.com\u002Fp\u002F72916159\u003C\u002Fp\u003E\u003Cp class=\"ql-align-right\"\u003E編輯:於騰凱\u003C\u002Fp\u003E\u003Cp class=\"ql-align-right\"\u003E校對:楊學俊\u003C\u002Fp\u003E\u003Cp class=\"ql-align-center\"\u003E\u003Cstrong\u003E— 完 —\u003C\u002Fstrong\u003E\u003C\u002Fp\u003E\u003Cp class=\"ql-align-justify\"\u003E關注清華-青島數據科學研究院官方微信公衆平臺“\u003Cstrong\u003ETHU數據派\u003C\u002Fstrong\u003E”及姊妹號“\u003Cstrong\u003E數據派THU\u003C\u002Fstrong\u003E”獲取更多講座福利及優質內容。\u003C\u002Fp\u003E\u003C\u002Fdiv\u003E"'.slice(6, -6), groupId: '6717827421392863751
相關文章