仿真預處理的目標是根據(jù)分析創(chuàng)建所需的網(wǎng)格。在生成同時解析幾何和物理特征的網(wǎng)格時，我們的目標是提高計算效率?；诜抡骖A期，我們可以針對流動特征不明顯的特定網(wǎng)格區(qū)域進行網(wǎng)格增強。在下圖中，我們在車輛周圍添加了細化區(qū)域，因為我們預計這部分會有大量的物理特征，尤其是在尾流區(qū)域。這一過程需要運用大量的專業(yè)領域知識，并且非常依賴用戶輸入的數(shù)據(jù)。不過，對流動物理特性極少的區(qū)域進行過度細化會增加計算成本和仿真時間，這是不可取的。本文將詳細介紹 Cadence Fidelity Pointwise 中的自動自適應網(wǎng)格細化功能，該功能可管理數(shù)值誤差并遵循用戶定義的邊界，同時解析各種應用中的所有流動特征。

1.網(wǎng)格自適應有哪些要求？

網(wǎng)格應滿足以下條件：

符合底層 CAD 幾何圖形

細化符合局部各向異性邊界層

細化較大的和細微的流動特征

在自適應網(wǎng)格細化過程中，我們希望保留邊界層和近壁面的物理特性，使單元大小平滑漸變，確保求解器正常收斂。在自適應網(wǎng)格細化過程中，有必要定義一個自適應傳感器，以便描述需要額外細化的區(qū)域。在飛機外部研究中，馬赫數(shù)是一個合適的自適應變量，而在渦輪機械應用中，可以將速度大小用作自適應傳感器。

2.使用 Fidelity Pointwise 進行自適應網(wǎng)格細化

在 Fidelity Pointwise 中創(chuàng)建基礎網(wǎng)格，以啟動自適應網(wǎng)格細化。然后，在該網(wǎng)格上運行解。在每個邊緣對傳感器進行評估，如果傳感器在特定位置超過了預定閾值，就會標記為自適應細化。根據(jù)點的位置和該區(qū)域新的目標單元大小創(chuàng)建點云。然后將點云整合到 Fidelity Pointwise 中的基礎網(wǎng)格中，創(chuàng)建更新后的網(wǎng)格。反復進行這一過程，直到解達到與網(wǎng)格無關的狀態(tài)。

Fidelity Pointwise 中的網(wǎng)格自適應周期

3.應用測試用例

沖擊射流

沖擊射流的邊界條件

在上圖中，冷射流降落到熱板上，邊界條件如圖所示。

基準網(wǎng)格（左），初始網(wǎng)格（右）

本例的目的是比較完整結構化網(wǎng)格和調(diào)整后的網(wǎng)格。在上圖中，左邊的圖片顯示的是基線基礎六面體網(wǎng)格，右邊的圖片顯示的是將根據(jù)速度大小進行調(diào)整的初始非結構化網(wǎng)格。

速度大小等值線（左）和自適應邊緣縮放（右）

右圖顯示了當前邊緣通道所需的縮放比例。自適應過程的重點是射流和板之間的區(qū)域。根據(jù)解創(chuàng)建點云，有 25% 的節(jié)點超過了閾值，因此將進行相應的自適應細化。

初始點云，有 25% 的節(jié)點被標記為自適應節(jié)點

在第五個周期中，約有 70% 的節(jié)點被標記為自適應節(jié)點，而在最后一個周期中，94% 的節(jié)點被標記為自適應節(jié)點。一旦約 90% 的節(jié)點被標記為自適應節(jié)點，就可以認為已經(jīng)收斂。此時可以結束迭代。

第 5 個周期（左）和第 9 個周期（右）中經(jīng)過自適應細化的網(wǎng)格

檢查網(wǎng)格統(tǒng)計數(shù)據(jù)后發(fā)現(xiàn)，與高度細化的六邊形網(wǎng)格相比，修改后的網(wǎng)格節(jié)點和元素更少。放大沖擊區(qū)域后，明顯發(fā)現(xiàn)初始網(wǎng)格并沒有準確捕捉到數(shù)據(jù)。然而，網(wǎng)格在經(jīng)過每個周期后都有所改進，最終更接近實驗數(shù)據(jù)。

努塞爾數(shù) (Nu) 在經(jīng)過每個自適應周期后都更接近實驗 Nu 值

軸流式渦輪機葉

亞琛渦輪機

本測試案例中是一臺亞琛渦輪機，有 41 個葉片，轉速為 3500 RPM。入口和出口流動條件如下表所示：

自適應細化的初始網(wǎng)格（左），自適應細化的最終網(wǎng)格 - 葉片區(qū)域的切面（右）

這里再次使用速度大小作為自適應變量。在經(jīng)過自適應細化的網(wǎng)格中，沖擊波清晰可見。此外，經(jīng)過自適應細化的最終網(wǎng)格準確地反映了二次渦和沖擊波的存在。

經(jīng)過自適應細化的網(wǎng)格捕捉到了二次渦流

DrivAer 模型

自適應網(wǎng)格細化也可用于汽車應用，本例中我們使用 DrivAer 模型作為測試案例。將速度大小作為自適應變量。DrivAer 模型的 RANS 仿真使用了 SST 兩方程湍流模型。經(jīng)過自適應的網(wǎng)格和尾流區(qū)域的流線如下所示，顯示出良好的匹配度和對渦流的準確捕捉。

DrivAer 模型經(jīng)過自適應的網(wǎng)格

外部空氣動力學

DLR F6 模型

該測試案例是第二屆 AIAA 阻力預測研討會上的 DLR F6 模型。馬赫數(shù)為 0.75，迎角為 1°。這里的自適應變量是馬赫數(shù)。通過自適應網(wǎng)格細化，機翼頂部的沖擊波清晰可見。初始表面壓力和經(jīng)過自適應細化后的表面壓力如下圖所示。經(jīng)過每個自適應周期后，沖擊波都會變得更加清晰。

初始網(wǎng)格（左）和經(jīng)過自適應細化的網(wǎng)格（右）的表面壓力

觀察升力和阻力系數(shù)可以發(fā)現(xiàn)，每經(jīng)過一個周期，細化程度都會變得更加精確。這種網(wǎng)格自適應方法可以輕松集成到任何工作流程中。

雖然一開始需要花費一些精力進行設置，但一旦設置妥當，剩下的流程就可以自動完成。在自適應周期中，因為每次都會返回到基線基礎網(wǎng)格，所以拓撲結構將保持一致。

適應周期期間的升力和阻力系數(shù)

4.高階網(wǎng)格彎曲

高階網(wǎng)格彎曲是一項新興技術，對于在計算流體力學（CFD）求解器領域中使用有限元方法（FEM）的設計人員來說，這項技術將成為他們的得力助手。與有限體積法和有限差分法等傳統(tǒng) CFD 方法相比，有限元技術能以較少的元素提高準確度。這要歸功于有限元方法在線性元素的邊、面和內(nèi)部引入了額外的頂點（新的自由度）。對于與彎曲幾何體相鄰的元素，這些新自由度必須位于幾何體上，從而改變原始線性元素的形狀。當網(wǎng)格中有一組元素向粘性邊界聚類時，這一過程會更加困難。內(nèi)部元素的邊和面也必須根據(jù)邊界元素的曲率進行彎曲，以避免元素反轉。

對 Cadence Fidelity Pointwise 所使用的 WCN 平滑處理進行研究，允許混合階網(wǎng)格求解幾何曲率。在高度彎曲的幾何體附近，元素最高可升階至四次多項式（四次方），而在遠離彎曲幾何體的地方，元素仍然保持線性。網(wǎng)格平滑方法使用代價函數(shù)來強制保持所需的元素形狀和每個元素的正雅可比。當元素在幾何體附近彎曲時，粘性網(wǎng)格間距保持不變。研究顯示了復雜 3D 配置的效果。