【文章內(nèi)容簡介】 一個Glyph Types范圍，Arrow提供一個清晰的向量平面圖給Glyph Types。在一個用戶可以操作其為最佳效果的面板，每個圖形都有圖形控制選項。剩下的Properties面板主要包括圖像的Scale Mode菜單，最常見的Scale Mode選項為：Vector，其中圖像長度正比于向量量級；Off,其中每個圖像都是一樣的長度。Set Scale Factor 參數(shù)控制圖像的基本長度。向量默認在單元頂點繪制，但通常希望在單元中心繪制數(shù)據(jù)。首先對案例模塊應(yīng)用Cell Centers ?lter，然后對產(chǎn)生的單元中心數(shù)據(jù)應(yīng)用Glyph ?lter。 流線繪制同樣的，在ParaView中繼續(xù)后處理之前，要使之前所描述的向量繪制等模塊退出，現(xiàn)在希望繪制速度流線。在Pipeline ，在Filter菜單選擇Stream Tracer，點擊應(yīng)用。指定Seed points（原點）沿著Line Source，其運行垂直于幾何中心，也就是說，從(, 0, )到 (, , )，在本教程中的圖像應(yīng)用：point Resolution為21；Max Propagation為Length ；初始Step Length 為 Cell Length ；以及Integration Direction BOTH。默認參數(shù)使用RungeKutta 2 Integrator Type 。點擊應(yīng)用產(chǎn)生軌跡，然后從Filter中選擇Tube來產(chǎn)生高質(zhì)量的流線圖。在本圖中，選擇：Nun. Sides 6。Radius 。Radius factor 。 增加網(wǎng)格分辨率 在每個方向上通過一個2因數(shù)增加網(wǎng)格分辨率，粗網(wǎng)格得出的結(jié)果可以映射到細網(wǎng)格上，作為問題的初始條件。然后將細網(wǎng)格的解與粗網(wǎng)格的解作對比。 用已存在的案例創(chuàng)建一個新的案例 現(xiàn)在希望根據(jù)cavity創(chuàng)建一個名為cavityFine的新案例，用戶需要克隆cavity案例，并且編輯必須的文件。首先，在于cavity相同的目錄下創(chuàng)建一個新的案例文件，如： cd $FOAM RUN/tutorials/inpressible/icoFoam mkdir cavityFine 然后從cavity案例中拷貝基本的文件到cavityFine中，然后進入cavityFine案例： cp r cavity/constant cavityFine cp r cavity/system cavityFine cd cavityFine 創(chuàng)建細網(wǎng)格 希望用blockMesh增加網(wǎng)格單元數(shù)。打開 blockMeshDict 文件，編輯block的指定信息。塊的指定在blocks關(guān)鍵詞下的列表中。；在該階段只需知道緊隨hex的首先是塊頂點的列表，然后是是每個方向上許多單元數(shù)的列表。最初在cavity案例中設(shè)為（20 20 1），現(xiàn)在改為（40 40 1），得到新的細網(wǎng)格。 將粗網(wǎng)格結(jié)果映射到細網(wǎng)格 mapFields 應(yīng)用程序?qū)⑴c一個給定幾何結(jié)構(gòu)相關(guān)的場映射到另外一個幾何結(jié)構(gòu)的相應(yīng)場中，在我們的例子中，流場認為是連續(xù)的，由于源場及目的場的幾何結(jié)構(gòu)及邊界類型或條件都是恒定的。在該例子中執(zhí)行 mapFields時使用 consistent 命令行。 從目標案例controlDict 中startFrom/startTime 指定的 時間目錄中讀取mapFields映射的流場數(shù)據(jù)，也就是說，結(jié)果被映射到的地方。在本例中，希望從cavity案例粗網(wǎng)格的最終結(jié)果映射到cavityFine案例中的細網(wǎng)格，因此，在controlDict文件中設(shè)置 ，startFrom設(shè)為 startTime。 現(xiàn)在案例已準備好運行mapFields,輸入mapFields help 快速顯示mapFields要求源案例目錄作為自變量。使用 consistent選項，所以應(yīng)用程在cavityFine目錄執(zhí)行： mapFields ../cavity consistent 應(yīng)用程序?qū)⑦\行并在終端輸出：Source: .. cavity Target: . cavityFine Create databases as time Source time: Target time: Create meshes Source mesh size: 400 Target mesh size: 1600 Consistently creating and mapping fields for time Open? Liddriven cavity flow U33 interpolating p interpolating U End 控制調(diào)整 為保持Courant數(shù)小于1，由于所有單元的尺寸都評分了所以時間步長必須評分，所以 controlDict 。流場數(shù)據(jù)在固定數(shù)目的時間步長間隔時輸出。此處演示如何在固定時間間隔指定數(shù)據(jù)輸出。 在controlDict 中的writeControl 關(guān)鍵詞下，在用runTime輸入產(chǎn)生的輸出結(jié)果間指定一個固定數(shù)量的運行時間，而不是通過timeStep 輸入固定書目的時間步來輸出。在該案例中，因此應(yīng)設(shè)置writeInterval 為runTime。最終，由于案例是在粗網(wǎng)格結(jié)果上開始的，只需要運行很短的時間就能達到合理的收斂為穩(wěn)態(tài)。所以endTime 。 運行代碼作為后臺過程用戶應(yīng)該體驗運行icoFoam作為后臺程序，重新定向終端輸出到log文件，者在之后可以看到。在 cavityFine 目錄下，用戶執(zhí)行： icoFoam log amp。 cat log 在細網(wǎng)格上繪制向量 用戶可以在ParaView中同時打開多個案例，這是由于本質(zhì)上每個新的案例就是Pipeline Browser 中一個新的模塊。當在ParaView中打開一個新的案例時有一個次要的不便，因為有一個前提，所選數(shù)據(jù)是一個帶有擴展名的文件。但是在OpenFOAM中，每個案例存儲在指定目錄結(jié)構(gòu)的大批沒有擴展名的文件中，paraFoam程序自動運行的結(jié)果，就是產(chǎn)生一個帶有擴展名 .OpenFOAM 的空白文件，因此。 然而，如果用戶希望在ParaView中直接打開另一個案例，就需要創(chuàng)建這樣一個空白文件。例如，為加載cavityFine 案例，通過輸入下列命令來創(chuàng)建該文件: cd $FOAM RUN/tutorials/inpressible/icoFoam touch cavityFine/ 現(xiàn)在通過在File菜單選擇Open， ，cavityFine 案例可以加載到ParaView中。用戶現(xiàn)在可以在ParaView中由細網(wǎng)格繪制向量圖。通過同時使兩個案例的glyph激活，圖像可以與cavity案例相比較。 繪制圖像用戶也許希望通過抽取一些速度標量及沿著穿過區(qū)域的軸線繪制2維圖像來使結(jié)果具體化。 OpenFOAM 對這種數(shù)據(jù)處理很在行。有無數(shù)的應(yīng)用程序可以做專門的數(shù)據(jù)處理，一些簡單的計算包括在一個單個的應(yīng)用程序foamCalc 中。作為一個應(yīng)用程序，它是獨特的，因為它是這樣的： foamCalc calcType fieldName1 ... fieldNameN 在calcType中指定計算器操作；在輸出時，采用下列操作：addSubtract。 randomise。 div。 ponents。 mag。 magGrad。 magSqr。 interpolate 。用戶可以通過故意調(diào)用一個不存在的操作來獲得calcType 列表，所以foamCalc 產(chǎn)生一個錯誤信息并且列出可用的類型，例如： foamCalc xxxx Selecting calcType xxxx unknown calcType type xxxx, constructor not in hash table Valid calcType selections are: 8 ( randomise magSqr magGrad addSubtract div mag interpolate ponents ) ponents 及 mag calcTypes 提供有用的速度標量。當在案例如cavity上運行“foamCalc ponents U ”時，它從每個時間文件中讀取速度矢量場，在相應(yīng)的時間文件中，輸出標量場Ux，Uy，Uz分別代表速度的x,y,z分量。類似的，“foamCalc mag mag U” 輸出一個標量場magU 到每個時間文件中，代表速度的量級。 用戶可以在 cavity及 cavityFine 案例上用ponents calcType 運行foamCalc 。例如，對于cavity案例，用戶應(yīng)該進入cavity目錄并執(zhí)行foamCalc如下： cd /home/ying/RUN/tutorials/inpressible/icoFoam/cavity foamCalc ponents U 單獨的分量可以在ParaView 中繪制成圖。這是很快速的，方便的，并且在標志及格式上有合理的好的控制，所以打印的輸出是相當好的標準品。但是為發(fā)表圖像，用戶也許偏好輸出未加工的數(shù)據(jù)，然后用專門的圖像工具來繪制，比如gnuplot 或 Grace/xmgr 。為此，我們推薦使用sample應(yīng)用程序。 在開始繪制之前，用戶需要加載新產(chǎn)生的Ux, Uy 及 Uz 場到ParaView 中。為此， 中，需要檢查Properties 面板頂部的Update GUI 按鈕。點擊應(yīng)用使新的場加載到ParaView 中，將出現(xiàn)在Vol Field Status 窗口。確保已選擇新的場且改變得到應(yīng)用，也就是說如果需要再次點擊Apply。同樣，如果在 Region Status 面板上選擇了boundary regions ，邊界上的數(shù)據(jù)插入錯誤。因此用戶應(yīng)該取消Region Status 面板中的邊界，也就是說movingWall, fixedWall 及frontAndBack ，然后應(yīng)用這些改變。 現(xiàn)在為了在 ParaView 中顯示圖像，用戶應(yīng)該選擇所感興趣的模塊， ，從FilterData Analysis 菜單應(yīng)用Plot Over Line 。這在已存在的3DView窗口旁邊打開了一個新的XY Plot窗口，創(chuàng)建了一個ProbeLine 模塊，在其中用戶可以指定Properties面板中線條的終點。在這個例子中，用戶應(yīng)該使線條在區(qū)域中心垂直向上 ，也就是說，在Point1 和 Point2 文本框中，從(, 0, ) 到(, , ) ，設(shè)置Resolution 為100. 點擊應(yīng)用，在XYPlot窗口出現(xiàn)圖像，在Display面板，選擇 網(wǎng)格分級在 cavityGrade /constant/polyMesh /blockMeshDict 文件中，制定了分級等級，blockMesh根據(jù)blockMeshDict 文件劃分網(wǎng)格。convertToMeters 。 vertices ( (0 0 0) ( 0 0) (1 0 0) (0 0) ( 0) (1 0)