Adjoint Solver實例<上>

2017-02-05  by:CAE仿真在線  來源:互聯(lián)網(wǎng)

本例來自于FLUENT文檔。

利用CFD進行流體計算,通常需要建立流體計算模型、設(shè)置邊界條件及初始條件、設(shè)置各種計算模型和控制參數(shù),從而輸出感興趣的物理量。在很多場合,往往需要關(guān)注輸出量對于輸入量的敏感性。通俗來說,就是輸入條件的改變對于輸出物理量的改變量有多大?FLUENT中的Adjoint Solver可以提供這些功能,讓使用者可以評價輸入變量對于輸出變量的影響程度。

本體以一個流體流經(jīng)圓柱的例子演示了如何產(chǎn)生敏感性數(shù)據(jù)、如何對數(shù)據(jù)進行后處理以及采用網(wǎng)格變形減小阻力。本例包含的內(nèi)容包括:

(1)如何加載adjoint solver

(2)如何選擇感興趣的觀察量

(3)如何進入求解器控制面板進行參數(shù)設(shè)置

(4)設(shè)置收斂標準以及顯示殘差曲線

(5)如何對計算數(shù)據(jù)進行后處理

(6)如何修改幾何形狀以減小阻力

本例要計算的模型為氣體流經(jīng)圓柱體,采用2D模型,雷諾數(shù)40,采用層流模型穩(wěn)態(tài)求解計算。最終目的為修改圓柱體外形以降低其阻力。在本例的后半部分通過采用升阻比最大化來優(yōu)化其氣動特性。

Adjoint solver是在模型求解之后進行的,因此本例從導入已計算完畢的casdat文件開始。

1、設(shè)置求解

Step 1:開啟FLUENT

采用2D、Double Precision啟動FLUENT。

選擇菜單項【File>Read>Case & Data…】讀取CASDat文件。

可以觀察網(wǎng)格及結(jié)果文件,如圖1所示。

1計算網(wǎng)格及速度場分布

Step 2:加載Adjoint Solver

Adjoint Solver為插件模型,需要在TUI中激活。

采用TUI命令 define/models/addon-module,出現(xiàn)如圖2所示的文本選擇命令。

2 TUI命令

輸入6選擇Adjoint Solver。

Step 3:定義Observable變量

本例定義的變量為圓柱所受的阻力。

點擊菜單【Adjoint>Observable】打開Adjoint Observables對話框,如圖3所示。

3 Adjont Observable對話框 4 Manage對話框

(1)點擊Manager…按鈕,進入Manager Adjoint Observables對話框,如圖4所示。

(2)點擊Create…按鈕,打開Create New Observable對話框,如圖5所示。

5 創(chuàng)建新的觀測器 6 Manage對話框

(3)選擇列表項中的force,更改變量名稱為force-drag,點擊OK按鈕確認操作并關(guān)閉對話框,同時返回至Manage Adjoint Observables對話框,如圖6所示。

(4)在圖6對話框中,選擇wall Zones列表項為wall,同時設(shè)置X Component1,Y Component0,表示檢測的力為x方向,即阻力。點擊OK按鈕關(guān)閉此對話框,返回至Adjoint Observables對話框,如圖7所示。選擇選項Minimize。

7 Adjont Observables對話框

(5)點擊圖7所示對話框中的Evaluate按鈕,則TUI窗口顯示所監(jiān)視的變量當前值,如圖8所示。當前監(jiān)測的阻力值為1337.8475N。

(6)點擊Close按鈕關(guān)閉對話框。

Step 4:定義Adjoint控制參數(shù)

選擇菜單【Adjoint>Contorls】,進入Adjoint Solution Contorls對話框,如圖8所示。

8 控制參數(shù)設(shè)置

(1)取消Use Stabilized Scheme選項,其他選項保持默認即可。

(2)點擊OK按鈕關(guān)閉對話框。

Step 5:設(shè)置監(jiān)視器

點擊菜單【Adjoint>Monitors】打開Adjoint Solver求解監(jiān)視器對話框,如圖9所示。

9 監(jiān)視器對話框

設(shè)置Adjoint continuityAdjoint velocity的殘差標準為1e-8,其他參數(shù)如圖所示設(shè)置。

Step 6:Adjoint計算

點擊菜單【Adjoint>Run Calculation】彈出如圖10所示對話框。

10 計算面板

(1)點擊按鈕Initialize進行初始化。

(2)設(shè)置Number of Iterations50

(3)點擊Calculate按鈕進行計算。

收斂殘差如圖11所示。

11 計算殘差

2、后處理

Step 1:邊界敏度輸出

點擊菜單【Adjoint>Reporting】進入Adjoint Reporting對話框,如圖12所示。

12 Adjoint Reporting對話框

選擇Boundary choice列表項inlet,在TUI窗口點擊按鈕Report輸出敏度信息。如圖13所示,可以看出,入口速度對于阻力的敏度為54.291264。

13 輸出信息

同樣可以Contours中查看各種敏度分布云圖。

Step 2:查看形狀敏度

點擊模型操作樹節(jié)點Graphics and Animations,選擇Vectors,點擊按鈕Set Up…,打開Vectors對話框,如圖13所示。

13 Vecotrs對話框

如圖所示進行設(shè)置。點擊按鈕Display,敏度矢量分布如圖14所示。

14 敏度分布

Step 3:修改形狀

(1)點擊菜單【Adjoint>Control-Volume Morphing】,彈出變形控制按鈕,如圖15所示。

15 變形控制定義

(2)在對話框中選擇邊界wall,點擊按鈕Get Bounds,下方參數(shù)框中將會列出wall邊界的幾何信息。

(3)點擊Larger Box數(shù)次,逐漸增大區(qū)域范圍至圖中所示。

(4)設(shè)置x pointy points參數(shù)值為20

點擊OK按鈕關(guān)閉對話框。

此時可以查看網(wǎng)格,如圖16所示。

16 網(wǎng)格顯示

(5)返回圖15所示對話框,設(shè)置Scale Factor參數(shù)值為0.1,點擊按鈕Update,此時激活其他按鈕。

(6)點擊按鈕Expected Change按鈕,阻力期望減少值顯示于TUI窗口。

(7)點擊按鈕Modify Mesh修改幾何形狀,點擊按鈕Accept確認設(shè)置。

此時重新查看網(wǎng)格,如圖17所示。

17 優(yōu)化后的網(wǎng)格形狀

(8)在優(yōu)化的網(wǎng)格基礎(chǔ)上重新進行計算。點擊Run Calculation節(jié)點,在面板中設(shè)置迭代次數(shù)100,點擊按鈕Calculate進行計算。收斂后查看阻力值,可以看到阻力值降低至1250.243N,相較之前的1337N有較大幅度的降低。

在氣動力學中,評價氣動特性常用升阻比,因此在下節(jié)中我們定義升阻比來對幾何形狀進行優(yōu)化。


開放分享:優(yōu)質(zhì)有限元技術(shù)文章,助你自學成才

相關(guān)標簽搜索:Adjoint Solver實例<上> Fluent培訓 Fluent流體培訓 Fluent軟件培訓 fluent技術(shù)教程 fluent在線視頻教程 fluent資料下載 fluent分析理論 fluent化學反應(yīng) fluent軟件下載 UDF編程代做 Fluent、CFX流體分析 HFSS電磁分析 

編輯
在線報名:
  • 客服在線請直接聯(lián)系我們的客服,您也可以通過下面的方式進行在線報名,我們會及時給您回復電話,謝謝!
驗證碼

全國服務(wù)熱線

1358-032-9919

廣州公司:
廣州市環(huán)市中路306號金鷹大廈3800
電話:13580329919
          135-8032-9919
培訓QQ咨詢:點擊咨詢 點擊咨詢
項目QQ咨詢:點擊咨詢
email:kf@1cae.com