摘  要： 提出一種簡單有效的、節約計算開銷的方法，用于制作實現三維場景中水域動畫效果。該方法的實現過程為利用二維模擬水流或波浪數據延拓到三維空間中，實現水沿其主導傳播方向產生波動動畫，其中利用了基于物理方程和粒子系統的方法。在動畫過程中，這些二維水流幀畫面被隨機地沿橫向方向取樣，在動畫時間和由用戶自定義噪聲的參數設置下取得的切片（稱為時間片）結構中進行。在重構三維水流畫面中，使用了帶平滑濾波器的幾何技術，以彌補因取樣所造成的切片損傷。使用隨機布朗運動噪聲數據補償三維畫面數據，逼真地顯示了漣漪波和怒濤波的動畫。

　　關鍵詞： 漣漪波；怒濤波；隨機延拓；仿真

0 引言

　　近些年來，計算機動畫技術表現得極為突出，取得了豐碩的成果。借助計算機圖形技術對各種自然現象的二維、三維動畫自動生成的研究方法不斷涌現，例如計算機模擬生成的爆炸、水流、波浪等動畫效果。這些自然現象的計算機生成方法給動畫制作帶來了極大的方便，為二維和三維動畫場景的制作縮短了周期、降低了成本。

　　對于創作大面積水域動畫，如海面的波浪飛濺及波浪洶涌澎湃的動畫時，因其過于復雜，計算量大，這在計算機動畫領域是一個極具挑戰的課題。本文提出一種簡單高效、節約計算量的方法來實現大面積海浪、波濤的動畫生成。其原理為利用水波動的二維仿真數據隨機延拓成三維水流波動效果。其實現過程如圖1所示。

1 相關研究文獻分析

　　計算機模擬水流動畫技術已經取得了長足發展。參考文獻[1]利用小振幅波理論和細胞自動機方法構建水波動畫的算法；參考文獻[2]使用參數方法描述水流動畫；參考文獻[3]基于水面的正弦關系建模描述大型海浪動畫及畫面渲染技術，在處理更復雜的水動力學有局限性；為獲得更為真實的水波動畫效果與飛濺效果，參考文獻[4]使用粒子系統結合高維場模型技術模擬水流動畫；而參考文獻[5]描述通過添加粒子產生噴霧和泡沫的水動畫的視覺效果；在參考文獻[6]中，表現一個全三維波浪粉碎效果的動畫是通過控制用戶定義的二維切片庫進行的；參考文獻[7]提出了基于粒子系統模擬的自適應采樣算法，著眼于幾何形狀復雜區域的計算，同時減少了流體深處的粒子數。

　　本文的方法由參考文獻[8]的思想改良而成，有關要點如下：

　?。?）引入新的隨機內容，提供了更好的可控性；

　　（2）隨機布朗運動（FBM）噪聲的使用，低/高通濾波器（LPF/HPF）和參數β（相對于波前的形狀/復雜性）的效果；

　?。?）多邊形生成渲染，包括表面粒子檢測、低通濾波器的切片取樣。

2 二維數據仿真方法

　　2.1 模擬仿真方法

　　一般來說，水動力學可以由如下微分方程描述：

　　其中，u是速度，t是時間，ρ為密度，p是壓力，ν是粘度，f是外力。

　　MPS方法[9]是通過粒子相互作用來解決式（1）和（2）的。對式（2）應用壓力映射后，MPS方法離散化并轉化成粒子相互作用的方程。粒子間的所有作用僅限于給定的截止半徑re內。兩粒子間相互作用的重量被定義為：

　　其中，r分別是兩粒子i和j在位置ri和rj間的距離。

　　在MPS方法中，流體的局部密度是用粒子數密度（PND）來表示的。粒子i的PND被定義為：

　　在自由表面的外側區域中不存在粒子時，自由表面上方的PND值減小。這樣，水粒子i滿足的條件為：

　　ni<δn0（6）

　　其中，n0是PND的標準值，δ是閾值，它們被認為是在自由表面上。

　　2.2 水粒子分類

　　在三維中重建流水表面，把水粒子根據它們的PND分類。這一簡單的分類的優勢在于，可以無需特殊處理或計算，就能得到水表面和飛濺水花粒子數的合理近似值。由以下條件根據式（6）進行分類：

　?。?）飛濺的水花粒子ni<δ1n0；

　?。?）水表面粒子δ1n0≤ni≤δ2n0；

　?。?）水中粒子ni>δ2n0。

　　δ1和δ2為用戶自定的閾值。

　　2.3 二維水波仿真

　　利用二維MPS方法對漣漪波和怒濤波兩類不同的水波進行模擬，其算法如下[10-11]：

　?。?）初始化模擬環境；

　?。?）for each時間步長；

　?。?）添加粒子；

　　（4）計算并賦予粒子作用力；

　?。?）平流粒子數；

　?。?）設置相鄰粒子和計算粒子數（PND）；

　?。?）提取水表面或飛濺水花粒子數；

　?。?）刪除粒子；

　?。?）end for。

　　程序實驗中使用了一些重要的仿真參數，各參數的值如表1所示。

3 三維動畫仿真

　　對上述的二維水波動畫方法，以隨機方式延拓成三維空間中的水波。對于三維表面重建，使用二維高斯濾波器。

　　3.1 隨機延拓方法

　　很容易想象，如果一個單一的二維波面在三維中被復制，所產生的表面在橫向方向上看起來絕對是均勻的，如圖2所示。

　　切片采樣過程以下列隨機公式取得：

　　其中，S為仿真動畫的幀；i為自然數，是動畫時間；k取值0，1，…，K-1，是動畫切片；N是噪聲函數；F為噪聲波動范圍參數；A為控制動畫速度的參數；中括號“[]”表示為運算后取最接近的整數。對于更換噪聲函數N，采用隨機布朗運動（fBm或也稱為1/fβ）的噪聲。fBm噪音完全由它的功率頻譜控制，與頻率f成反比，參數β為加權系數。

　　根據式（7）和式（8），幀S為在某個動畫時刻i的第k個切片，是根據隨機式（7）和在約束條件下式（8）下所取得的。如圖1所示。一旦所有的切片都填滿了由隨機取樣的幀，這些切片被線性采樣以構造三維水流波形。

　　3.2 曲面近似

　　為了處理三維表面重建，使用帶高斯濾波器的高度場技術。這能夠重建近似的三維水面，而補救線狀偽影所造成的二維切片的采樣。

　　在XZ平面上創建大小為N×M基礎網格如下：

　　為了取得更好的結果，在兩個相鄰的切片和截面積（3Δx，2Δz）高斯濾波參數σ=2.0之間插入一個額外的切片。

3.3 三維漣漪波

　　對于模擬漣漪波時，通過在不同方向合成幾類波的方式略有不同。用這種方法，消除了沿橫向方向的波形外觀，它產生在單一方向上的運動。為了最終的h（P）值在每個小區中心P上的基礎網格IXZ逼近，由式（10）合成臨近點在每個方向上的高度。

　　其中，ω是用戶自定義重量值，hd是相鄰點Q在d方向上的高度值。

4 實驗效果

　　根據前述的方法和相應程序算法，在Intel（R）Core（TM）主頻為2.83 GHz，內存為4 GB的環境下，利用Delphi結合OpenGL三維引擎庫進行編程實現。漣漪波和怒濤波在不同的布朗隨機運動噪聲參數下取得了不同動畫效果，如圖3、圖4所示。

5 結論

　　本文提出了利用二維水面波動曲線延拓形成三維的水波動畫效果，利用各類技術補償二維切片損傷重構三維曲面，形成逼真的海水波濤及漣漪動畫效果，在實驗過程中較同類方法實現的三維海浪波濤動畫有相當的效率。節約使用大量的粒子數，為后期動畫效果的渲染節約了系統資源。此研究為動畫制作中的大型水域動畫場景的計算機生成提供了技術支持，節約了大量類似開發的成本并縮短制作周期。此項技術方法將繼續向GPU領域邁進，以期在今后的研究中能更好地以快速渲染自動生成大型水域中的波濤洶涌、水浪渦旋等動畫效果。

參考文獻

　　[1] 楊懷平，胡事民，孫家廣.一種實現水波動畫的新算法[J].計算機學報，2002，25（6）：612-617.

　　[2] FOURNIER A， REEVES T. A simple model of ocean waves[J]. Proceedings of SIGGRAPH， 1986（20）：75-84.

　　[3] DARLES E， CRESPIN B， GAZANFARPOUR D. A survey of ocean simulation and rendering techniques in computer graphics[J]. Computer Graphics Forum，2011，30（1）：43-60.

　　[4] MAES M M， FUJIMOTO T， CHIBA N. Efficient animation of water flow on irregular terrains[C]. Proceedings of International Conference on Computer Graphics and Interactive Techniques， GRAPHITE′06， 2006：107-115.

　　[5] TAKAHASHI T， FUJII H， KUNIMATSU A， et al. Realistic animation of fluid with splash and foam[J]. Computer Graphics Forum， 2003，22（3）：391-400.

　　[6] ADAMS B， PAULY M， KEISER R， et al. Adaptively sampled particle fluids[J]. ACM Transactions on Graphics，2007，26（3）：48.

　　[7] HONG W， HOUSE D， KEYSER J. Adaptive partcles for incompressible fluid simulation[J]. The Visual Computer，2008，24（7）：535-543.

　　[8] FUJIMOTO T， MIYAUCHI S， SUZUKI T， et al. Noise-based animation of waving phenomena[C]. Proceedings of   IWAIT2005， IEICE Technical Report. Image Engineering， 2005，104（545）：93-98.

　　[9] KOSHIZUKA S， NOBE A， OKA Y. Numerical analysis of breaking waves using the MPS method[J]. International Journal for Numerical Methods in Fluids， 1998（26）：751-769.

　　[10] 李海，汪繼文.基于粒子系統的位置可控有聲煙花模擬[J].微型機與應用，2013，32（9）：32-34.

　　[11] 袁雪霞，王繼州.導彈煙霧的快速模擬[J].微型機與應用，2013，32（10）：32-34.