張 南，金 生，艾叢芳

（大連理工大學(xué)海岸和近海工程國家重點(diǎn)實(shí)驗(yàn)室，大連116024）

基于WPF的三維流場可視化研究

張 南，金 生，艾叢芳

（大連理工大學(xué)海岸和近海工程國家重點(diǎn)實(shí)驗(yàn)室，大連116024）

通過對多種三維流場仿真算法的研究，以WPF圖形平臺為基礎(chǔ)，采用可視化編程語言開發(fā)了三維流場可視化系統(tǒng)。在三維地形、表面流場可視化的基礎(chǔ)上，提出了置于模型域的多剖面流場并疊加標(biāo)量場的三維流場呈現(xiàn)方式。同時(shí)采用示蹤球追蹤以及標(biāo)識點(diǎn)跡線描述的方式反映流場的結(jié)構(gòu)特性以及三維流場局部細(xì)節(jié)?？梢暬到y(tǒng)通過異步處理方式高效、流暢地呈現(xiàn)仿真流場。以三維潰壩模擬計(jì)算為例，驗(yàn)證了可視化平臺的實(shí)用性。

WPF；三維流場；剖面流場；跡線；示蹤球；動(dòng)畫

隨著三維數(shù)值算法的日益成熟，越來越多的實(shí)際工程問題采用三維計(jì)算方法來解決［1］。但是目前矢量場的可視化技術(shù)還遠(yuǎn)不能滿足實(shí)際的要求，這主要表現(xiàn)在缺乏有效的三維流場表示方法、三維交互技術(shù)和實(shí)時(shí)動(dòng)態(tài)顯示技術(shù)。因此，如何充分識別與分析流場信息仍然是流場表示方法和技術(shù)中具有挑戰(zhàn)性的研究課題［2-3］。龐大的數(shù)據(jù)量使得人們不得不尋求新的更為精密復(fù)雜的可視化算法和工具來分析三維流場。

基于OpenGL的三維地形可視化，是以往常見的對三維場景的實(shí)現(xiàn)方法形式［4-5］。另外一些商業(yè)化軟件（Tecplot、Mike21，Delft3D）雖然在呈現(xiàn)三維地形以及表面流場方面有一定優(yōu)勢，但在表現(xiàn)剖面流場方面，是將設(shè)定的剖面以二維圖形的方式在畫布上重新繪制，本質(zhì)上是二維方式。

本文從應(yīng)用WPF構(gòu)建三維場景的基本步驟出發(fā)，詳細(xì)描述了三維地形，表面流場，剖面流場，跡線及示蹤球以及流場動(dòng)畫的呈現(xiàn)方法，實(shí)現(xiàn)了歐拉法和拉格朗日法兩種方式的三維流場顯示，并且通過顏色材質(zhì)映射完成了對標(biāo)量場的模擬。

1基于WPF的三維場景建立

WPF是微軟新一代圖形系統(tǒng)。它是基于向量的呈現(xiàn)引擎，并且與分辨率無關(guān)，相對于DirextX是圖形硬件的底層接口，WPF提供的是上層抽象（底層圖形系統(tǒng)構(gòu)建于DirectX之上），它描述了你的場景（scene），只要硬件支持，它就會找出最好的呈現(xiàn)方法。

盡管WPF提供了方便的建模方法和三維圖形顯示類，但是構(gòu)建三維場景仍然需要復(fù)雜的工作。根據(jù)WPF三維建模需要的4個(gè)元素，將三維模型的建立歸納為如圖1步驟：

圖1 基于WPF的三維模型的建立步驟Fig.1 Process of establishing 3D model based on WPF

上述步驟主要是完成：（1）設(shè)置攝像機(jī)的位置和光線的方向和顏色等使三維場景看上去顏色、大小合適。（2）坐標(biāo)變換使模型的世界坐標(biāo)在屏幕上得以顯示。（3）建立三維模型。三維模型是通過三維曲面表現(xiàn)的，三維曲面通過網(wǎng)格化幾何圖形加材質(zhì)以及紋理映射到表面實(shí)現(xiàn)的。（4）在xaml文件中設(shè)計(jì)呈現(xiàn)三維場景的視口即viewport3D控件，用來顯示所有的三維空間的曲面。

2建立三維地形

地形也是以三維曲面方式表現(xiàn)的。

（1）網(wǎng)格化地形數(shù)據(jù)。描述地形的網(wǎng)格數(shù)據(jù)已經(jīng)在數(shù)值模擬前處理的時(shí)候生成，因此這部分不需要多做額外的工作，只需要讀取描述地形的網(wǎng)格數(shù)據(jù)即可，得到網(wǎng)格節(jié)點(diǎn)坐標(biāo)信息和連接關(guān)系信息。將節(jié)點(diǎn)和連接關(guān)系信息賦值給Positions和TriangleIndices屬性，WPF會自動(dòng)創(chuàng)建MeshGeometry3D對象，從而完成地形數(shù)據(jù)的網(wǎng)格化。

（2）建立三維地形模型。建立了描述三維曲面的網(wǎng)格或者基元，需要設(shè)置映射到圖形表面的材質(zhì)和紋理映射。通常選擇DiffuseMaterial材質(zhì)，它最能夠反映真實(shí)世界里平面向各個(gè)方面散射光的特性。常規(guī)情況都采用漸變顏色畫刷作為紋理映射到三維表面。顏色從藍(lán)過渡到紅，表示數(shù)值從小到大。

（3）對三維模型應(yīng)用變換矩陣。通常描述地形的網(wǎng)格數(shù)據(jù)采用世界坐標(biāo)系，因此必須采用適當(dāng)?shù)淖鴺?biāo)變換才能在屏幕上呈現(xiàn)出來。首先根據(jù)屏幕范圍以及需要顯示的整個(gè)模型來確定水平方向的縮放比例，為了圖形不失真，設(shè)置豎直方向的縮放比例與水平方向的縮放比例相同。WPF提供了Transform3DGroup對象進(jìn)行三維坐標(biāo)變換。它包括三種坐標(biāo)變換方式，即平移（TranslateTransfrom3D）、縮放（SacleTransfrom3D）以及旋轉(zhuǎn)（RotateTransform3D）。在設(shè)置應(yīng)用于地形的Model3DGroup對象的三維坐標(biāo)變換時(shí)，同樣需要設(shè)置應(yīng)用于攝像機(jī)的三維坐標(biāo)變換以到達(dá)預(yù)期的視覺效果。

（4）呈現(xiàn)三維地形。三維地形通過viewport3D類就可以呈現(xiàn)出來。只需要將MeshGeometry3D類轉(zhuǎn)換為ModelVisual3D即可。因?yàn)槿S的場景依賴于觀察者的觀察點(diǎn)。所以必須規(guī)定觀察點(diǎn)和光線。Camara類用來設(shè)置三維場景的觀察點(diǎn)。另外三維圖像里必須規(guī)定燈光效果使表面可見。更重要的一點(diǎn)，燈光確定了場景的哪一部分將被納入投影。Light對象用來創(chuàng)建各種燈光和陰影效果，仿照各種真實(shí)世界的燈光行為。

3表面流場與截面流場模擬

3.1表面流場模擬

單從水面的流動(dòng)也能夠認(rèn)識水流的部分三維特性，例如天然河道流有明顯的翻滾流、河道彎曲處有橫向流，這些現(xiàn)象都可通過表層水流運(yùn)動(dòng)而認(rèn)識。由于液體質(zhì)點(diǎn)運(yùn)動(dòng)的復(fù)雜性，在水力學(xué)中，用拉格朗日法研究液體運(yùn)動(dòng)將會遇到很多困難。同時(shí)，在解決實(shí)際問題時(shí)一般也不需要逐個(gè)弄清液體質(zhì)點(diǎn)的來龍去脈。因此，除個(gè)別問題（如波浪運(yùn)動(dòng)）外，在水力學(xué)中，普遍運(yùn)用歐拉法來分析和研究液體運(yùn)動(dòng)。目前大都以平面非結(jié)構(gòu)化、垂向分層方式構(gòu)建三維立體網(wǎng)格，因此表面流場的模擬可以擴(kuò)展成任意層的表面流場模擬。

為了方便實(shí)現(xiàn)表面流場的模擬，將繪制三維線段的方法封裝在ScreenSpaceLines3D類中，在繪制的時(shí)候給ScreenSpaceLines3D類的Points Thickness，Color屬性賦值，WPF自動(dòng)創(chuàng)建三維線段對象。值得注意的是將所有表示流場的矢量箭頭作為一個(gè)三維對象呈現(xiàn)在三維場景中，用以節(jié)省系統(tǒng)的開銷。

表面流場的繪制算法流程如下：

遍歷顯示層（ilayer）的所有網(wǎng)格節(jié)點(diǎn)以及節(jié)點(diǎn)上物理量的值。并找出該層的平面流場的最大流速作為繪制流速箭頭的比尺，以及各個(gè)網(wǎng)格節(jié)點(diǎn)所在層的x，y，z方向的流速u，v，w。然后根據(jù)該網(wǎng)格節(jié)點(diǎn)的流速值計(jì)算出描述流速的箭頭矢量，箭頭的長度表示流速大小，方向表示流速的方向。

流速矢量箭頭起點(diǎn)就是網(wǎng)格節(jié)點(diǎn)（x1，y1，z1），終點(diǎn)坐標(biāo)（x2，y2，z2）根據(jù)流速計(jì)算，如公式（1），其中s表示坐標(biāo)縮放比例。

循環(huán)計(jì)算該層所有網(wǎng)格節(jié)點(diǎn)的流場矢量箭頭（3條三維空間線段），利用全局坐標(biāo)變換，以WPF技術(shù)為基礎(chǔ)繪制。

若要設(shè)置流場的疏密，則需要設(shè)置分的份數(shù)（nx、ny），然后根據(jù)計(jì)算域的x，y方向的最大最小值確定標(biāo)識點(diǎn)的個(gè)數(shù)和坐標(biāo)，查詢離標(biāo)識點(diǎn)最近的網(wǎng)格節(jié)點(diǎn)作為描述流場的網(wǎng)格節(jié)點(diǎn)。獲取該網(wǎng)格節(jié)點(diǎn)上的物理量的值，計(jì)算流速箭頭，繪制在屏幕上。

如圖2所示，以三維潰壩的計(jì)算結(jié)果為例，表面三維流場的模擬。

圖2 表面三維流場的模擬Fig.2 3D flow field simulation of surface

3.2剖面流場模擬

在三維計(jì)算數(shù)據(jù)顯示中，剖面顯示是常用的手段，但是目前的商業(yè)化三維軟件（surfer、tecplot等）只能重新定義一個(gè)畫布以二維圖形的方式分析剖面數(shù)據(jù)，這樣不能直觀地顯示剖面在原計(jì)算域的位置，且不利于對多個(gè)剖面數(shù)據(jù)對比分析。因此，置于原計(jì)算域多剖面三維顯示的方法對現(xiàn)有三維剖面顯示進(jìn)行補(bǔ)充。

基于WPF呈現(xiàn)三維圖形的要素，截面流場繪制需要以下3個(gè)步驟。

（1）網(wǎng)格化剖面。

剖面是根據(jù)平面上設(shè)定的若干條剖線生成。剖線的起始點(diǎn)的坐標(biāo)可能在計(jì)算域內(nèi)或者計(jì)算域外，因此需要重置剖線的起始點(diǎn)坐標(biāo)。如果剖線兩個(gè)端點(diǎn)在計(jì)算域外，那么與計(jì)算域外邊界有兩個(gè)交點(diǎn)，視為剖線穿過計(jì)算域計(jì)入剖面線，并將剖線與計(jì)算域外邊界的交點(diǎn)作為剖線新的端點(diǎn)。如果剖線與計(jì)算域外邊界有一個(gè)交點(diǎn)或者沒有交點(diǎn)，視為剖線不穿過計(jì)算域則不計(jì)入剖面線。如果一個(gè)端點(diǎn)在計(jì)算域外，一個(gè)端點(diǎn)在計(jì)算域內(nèi)，那么將剖線與計(jì)算域的交點(diǎn)以及計(jì)算域內(nèi)的點(diǎn)作為剖線新的端點(diǎn)。如果剖線兩個(gè)端點(diǎn)在計(jì)算域內(nèi)則不做其他操作，直接計(jì)入剖面線。

由于計(jì)算網(wǎng)格采用平面非結(jié)構(gòu)化，垂向分層的方式劃分，因此剖面對計(jì)算域三維實(shí)體進(jìn)行切割時(shí)，剖面是由若干四邊形組成，這些四邊形分層排列，左右兩邊平行，上下兩邊不一定平行。

為了剖面上的四邊形更加均勻，不將剖線與平面網(wǎng)格的交點(diǎn)作為剖面網(wǎng)格化的網(wǎng)格節(jié)點(diǎn)，而是將剖線按網(wǎng)格尺度均為分成若干點(diǎn)n=dl/gl，其中dl是剖線長度，gl是網(wǎng)格尺度。經(jīng)過線性插值計(jì)算后得到的點(diǎn)作為剖線上的標(biāo)識點(diǎn)。這時(shí)剖面就由若干相對規(guī)則的四邊形組成。為了疊加剖面上的標(biāo)量場，也就是計(jì)算剖面上網(wǎng)格節(jié)點(diǎn)物理量的值，需要重置剖線上的標(biāo)識點(diǎn)。具體算法為，查找出離各個(gè)標(biāo)識點(diǎn)最近的網(wǎng)格節(jié)點(diǎn)，剔除最近點(diǎn)是同一個(gè)網(wǎng)格節(jié)點(diǎn)的點(diǎn)，剔除兩點(diǎn)間的距離小于網(wǎng)格尺度的點(diǎn)。重置后的剖線標(biāo)識點(diǎn)保存在哈希表（HashTable）中。其中key存放網(wǎng)格節(jié)點(diǎn)號，value存放該網(wǎng)格節(jié)點(diǎn)所在的單元（一個(gè)節(jié)點(diǎn)在一個(gè)或者多個(gè)單元里）。有了剖線上的標(biāo)識點(diǎn)坐標(biāo)，因?yàn)榇瓜蚴蔷鶆蚍謱?，根?jù)該時(shí)刻自由水面的高程值以及底高程值，很容易計(jì)算出剖面上各層標(biāo)識點(diǎn)的坐標(biāo)，也就是繪制剖面的網(wǎng)格節(jié)點(diǎn)。然后將每個(gè)小四邊形分解成兩個(gè)小三角形，也就是繪制剖面的連接關(guān)系。從而完成剖面的網(wǎng)格化。

（2）計(jì)算剖面上標(biāo)識點(diǎn)的標(biāo)量值。

上面已經(jīng)計(jì)算出離標(biāo)識點(diǎn)最近的網(wǎng)格節(jié)點(diǎn)所在的單元，用面積法判斷該標(biāo)識點(diǎn)具體在哪個(gè)三角形內(nèi)，以確定標(biāo)識點(diǎn)所在的單元號。用面積插值算法求出剖線上標(biāo)識點(diǎn)的物理量的值。

設(shè)網(wǎng)格頂點(diǎn)A、B、C，變量的值為v1、v2、v3，標(biāo)識點(diǎn)P，變量為v，則P點(diǎn)的變量值為

垂向上根據(jù)標(biāo)識點(diǎn)所在垂向?qū)拥奈恢帽壤?，插值得到其他層網(wǎng)格節(jié)點(diǎn)的物理量值。

（3）通過顏色的紋理映射疊加標(biāo)量場。

需要計(jì)算所要疊加的標(biāo)量（水位、水深、流速、壓強(qiáng)等）的最大、最小值，模擬開始前建立標(biāo)量場對應(yīng)的色譜表，根據(jù)各剖面節(jié)點(diǎn)位置處標(biāo)量值大小，通過插值得到對應(yīng)的顏色值，將顏色值以填充畫刷的顏色體現(xiàn)在截面紋理上以實(shí)現(xiàn)截面標(biāo)量場的顯示。

繪制剖面上的流速矢量場跟繪制表面流場的方法相似，也采用歐拉描述，所不同的是需要將流速向剖面投影，并計(jì)算出描述剖面流場的矢量箭頭。箭頭投影的方法參見文獻(xiàn)［6］。

如圖3所示，以三維潰壩的計(jì)算結(jié)果為例，采用多剖面的方式表現(xiàn)的三維流場。

圖3 多剖面流場模擬Fig.3 3D flow field simulation of multi?cross?sections

4跡線與示蹤球顯示

質(zhì)點(diǎn)的運(yùn)動(dòng)軌跡描述同樣可以幫助用戶更加方便的分析三維流場。其中跡線及示蹤球是比較常用的追蹤質(zhì)點(diǎn)運(yùn)動(dòng)軌跡的方式。在示蹤球三維模擬的過程中，在計(jì)算域內(nèi)布設(shè)一些示蹤球，追蹤不同時(shí)刻示蹤球所在的位置。在跡線的三維模擬過程中，追蹤一個(gè)剖面上若干標(biāo)識點(diǎn)在不同時(shí)刻所在的位置，并連接起來，形成描述標(biāo)識點(diǎn)運(yùn)動(dòng)軌跡的空間曲線。

質(zhì)點(diǎn)追蹤之前，需要在計(jì)算域內(nèi)布設(shè)標(biāo)識點(diǎn)，根據(jù)顯示的不同，進(jìn)行不同的初始化設(shè)置。

4.1初始化示蹤球

將能包含計(jì)算域的最小矩形區(qū)在x、y向等分為nx，ny份。根據(jù)示蹤球的坐標(biāo)，查找示蹤球所在的網(wǎng)格單元，并計(jì)算示蹤球所在位置初始時(shí)刻的流速值，算法同計(jì)算剖面流場標(biāo)識點(diǎn)的流速相同。示蹤球可以描述表層或者全部層的數(shù)據(jù)，經(jīng)實(shí)踐檢驗(yàn)最多顯示3層的數(shù)據(jù)，因?yàn)槿绻@示的層太多，界面上示蹤球就會重疊遮擋，反而顯示不清晰。

4.2初始化跡線

首先根據(jù)設(shè)定的剖線坐標(biāo)判斷剖面的有效性，根據(jù)劃分份數(shù)n，計(jì)算平面上標(biāo)識點(diǎn)的坐標(biāo)，然后根據(jù)數(shù)值模擬計(jì)算時(shí)分的層數(shù)nL，計(jì)算垂向上標(biāo)識點(diǎn)的坐標(biāo)，計(jì)算標(biāo)示點(diǎn)初始時(shí)刻的流速值。跡線及示蹤球的三維可視化效果如圖4所示。

圖4 跡線及示蹤球三維流場模擬Fig.4 3D flow field simulation of trace and tracer ball

4.3動(dòng)態(tài)生成算法

初始化完成后，在計(jì)算粒子經(jīng)過Δt時(shí)間后的位置的時(shí)候采用一階時(shí)間精度模擬。也就是說假設(shè)示蹤粒子以起始位置處的流速勻速運(yùn)動(dòng)。到下一個(gè)時(shí)刻，計(jì)算示蹤粒子新的流速值，并假設(shè)示蹤粒子又以新的流速勻速運(yùn)動(dòng)，循環(huán)往復(fù)，計(jì)算出流場的示蹤粒子的運(yùn)動(dòng)軌跡。

本質(zhì)上是認(rèn)為在足夠小的時(shí)間段內(nèi)示蹤粒子做勻速直線運(yùn)動(dòng)。把長時(shí)間段劃分為若干小時(shí)間段線性插值計(jì)算示蹤粒子新的位置的流速值。實(shí)際上是把示蹤粒子的曲線運(yùn)動(dòng)軌跡看作是若干個(gè)短距離直線運(yùn)動(dòng)組成［7-8］。采用時(shí)間的一階精度計(jì)算跡線的時(shí)候，保存的2個(gè)時(shí)刻的計(jì)算結(jié)果的時(shí)間間隔不能太長，示蹤粒子的每步軌跡也不能太大，否則將不足夠反應(yīng)流場的特性。如果采用高階精度，就需要多一次插值過程，計(jì)算增大很多，如果保存的時(shí)候時(shí)間間隔大，中間采用高階插值也沒有必要。Δt足夠小的話，采用一階精度可以滿足模擬的要求。

圖5 跡線及示蹤球的計(jì)算算法流程Fig.5 Algorithm flowchart of trace and tracer ball

一般情況跡線以及示蹤球的計(jì)算時(shí)間步長會比數(shù)值模擬計(jì)算的時(shí)間步長小一些。設(shè)跡線計(jì)算的當(dāng)前時(shí)刻為t，數(shù)值模擬計(jì)算的前一時(shí)刻為t0，下一時(shí)刻為t1，線性插值計(jì)算跡線上示蹤點(diǎn)的流速。算法如下。

跡線上的標(biāo)識點(diǎn)經(jīng)過Δt時(shí)間后的新的位置x+=u*Δt；y+=v*Δt；z+=w*Δt；根據(jù)新的位置確定所在的網(wǎng)格單元，如果該示蹤球已經(jīng)流出計(jì)算域，那么刪除該示蹤球。

如果計(jì)算的當(dāng)前時(shí)刻t+Δt＞t1，那么將原來的后一時(shí)刻的值賦給前一個(gè)時(shí)刻，重新讀取下一個(gè)時(shí)刻的流速值，直到t＞t總。算法的詳細(xì)流程如圖5所示。

5流場動(dòng)畫顯示

流場的動(dòng)畫演示，實(shí)際上就是按計(jì)算的時(shí)間順序呈現(xiàn)不同時(shí)刻的三維流場，標(biāo)量場和矢量的繪制方法以及步驟上面已經(jīng)介紹過，這里不再贅述。以WPF技術(shù)為基礎(chǔ)，流場的三維動(dòng)畫實(shí)現(xiàn)起來相對簡單，定義一個(gè)委托delegate用來實(shí)現(xiàn)繪制某一時(shí)刻三維流場，采用Dispatcher.BeginInvoke（）方法異步執(zhí)行該委托，這樣擺脫單個(gè)任務(wù)的牽制，提高三維圖形的呈現(xiàn)效率，能夠?qū)崿F(xiàn)在動(dòng)畫模擬的時(shí)候快速的響應(yīng)用戶的鼠標(biāo)指令（平移、縮放，旋轉(zhuǎn)）等。

6結(jié)論

從研究構(gòu)建三維流場的算法出發(fā)，應(yīng)用先進(jìn)的三維圖形平臺（WPF），實(shí)現(xiàn)了流場的歐拉法及拉格朗日法兩種表達(dá)方式。表面流場和剖面流場采用歐拉描述，以標(biāo)識點(diǎn)跡線和示蹤球方式追蹤粒子的運(yùn)動(dòng)軌跡來表現(xiàn)拉格朗日場。通過置于原計(jì)算域的多剖面的三維矢量場以及疊加標(biāo)量場的表現(xiàn)方法，直觀方便的對比分析剖面流場。同時(shí)按計(jì)算時(shí)間順序呈現(xiàn)三維流場，實(shí)現(xiàn)了實(shí)時(shí)、快速、平滑的流場矢量動(dòng)畫。為了更好的分析及表現(xiàn)三維流場，在構(gòu)建每種流場的時(shí)候都可以靈活地對流場各參數(shù)進(jìn)行設(shè)置。

［1］劉曉波，華祖林，何國建.計(jì)算流體力學(xué)的科學(xué)計(jì)算可視化研究進(jìn)展［J］.水動(dòng)力學(xué)研究與進(jìn)展，2004，19（1）：120-125. LIU X B，HUA Z L，HE G J.Advances in VISC of CFD［J］.Journal of Hydrodynamic，2004，19（1）：120-125.

［2］魯大營，朱登明，王兆其.三維流場的流線提取算法［J］.計(jì)算機(jī)輔助設(shè)計(jì)與圖形學(xué)學(xué)報(bào)，2013，25（5）：666-673. LU D Y，ZHU D M，WANG Z Q.Streamline Selection Algorithm for Three?Dimensional Flow Fields［J］.Journal of Computer?Aided Design&Computer Graphics，2013，25（5）：666-673.

［3］袁曉茹，張昕，肖何，等.可視化研究前沿及展望［J］.科研信息化技術(shù)與應(yīng)用，2011，2（4）：3-13. YUAN X R，ZHANG X，XIAO H，et al.Visualization Research Frontier and Prospect［J］.e?Science Technology&Application，2011，2（4）：3-13.

［4］吳杰，黃春生，范緒箕.基于OpenGL的CFD設(shè)計(jì)平臺中的流場可視化技術(shù)及其實(shí)現(xiàn)［J］.工程圖學(xué)學(xué)報(bào)，2004（2）：61-72. WU J，HUANG C S，F(xiàn)AN X Q.Visualization Technique of FlowFields on CFD Design Platform UsingOpenGL System［J］.Journal of Engineering Graphics，2004（2）：61-72.

［5］陳立華，梅亞東，王現(xiàn)勛，等.基于OpenGL三維河網(wǎng)地形與數(shù)據(jù)場的可視化［J］.武漢大學(xué)學(xué)報(bào)：工學(xué)版，2007，40（3）：34-37. CHEN L H，MEI Y D，WANG X X，et al.Visualization of 3D River Network Terrain and Data Field Based on OpenGL［J］.Engineer?ing Journal of Wuhan University，2007，40（3）：34-37.

［6］康蘇海，金生，劉洋.三維場景中截面流場繪制技術(shù)研究［J］.水道港口，2011，32（5）：361-366. KANG S H，JIN S，LIU Y.Technical Research on Cross?section Flow Drawing in 3?D scene［J］.Journal of Waterway and Harbor，2011，32（5）：361-366.

［7］康蘇海.三維流場立體仿真軟件研發(fā)［J］.水道港口，2010，32（5）：324-329.KANG S H.Software Development of 3?D Flow Simulation Using Anaglyph Technology［J］.Journal of Waterway and Harbor，2010，32（5）：324-329.

［8］KANG S H，JIN S.App1ication of video anag1yph maker for 3d flow simulation［J］.Journal of Hydrodynamics，2010，22（2）：289-294.

［9］Adam Nathan.WPF揭秘［M］.北京：人民郵電出版社，2008.

Study of three dimensional flow field visualization based on WPF

ZHANG Nan,JIN Sheng,AI Cong?fang
(State Key Laboratory of Coastal and Offshore Engineering,Dalian University of Technology,Dalian 116024,China)

According to the analysis of kinds of flow field simulation algorithm,3D flow field visualization sys?tem was developed based on WPF platform and visual programming language.On the basis of visualization of 3d to?pography and surface flow field,multiple cross?section flow placed into original computational domain was pro?posed.At the same time,the description of tracer sphere and path line was adopted to reflect the structure character?istics of the flow field characteristics and local details of the flow field.The 3D visualization system can present the simulation of flow field efficiently and smoothly by asynchronous update.On the basis of dam break computation,it is demonstrated that the visual platform can be viewed as a practical tool for 3D flow field analysis.

WPF;three dimensional flow field;cross?section flow;path line;tracer sphere;animation

TP 391.41

A

1005-8443（2016）04-0455-06

2016-01-22；

2016-04-11

張南（1983-），女，遼寧省人，博士研究生，主要從事數(shù)值模擬和流場仿真模擬研究。

Biography：ZHANG Nan（1983-），female,doctor student.