CFD前處理網格藝術 | CFD對計算網格的基本要求
利用數(shù)值計算方法得到的離散解是否比較滿意地逼近原偏微分方程組定解問的解,不僅取決于對原偏微分方程組所采用的離散化方法(即內點計算格式)及邊界條件的離散化方法(即邊界點計算格式),而且取決于離散點的分布情況。
另一方面,許多流體力學實際問題的邊界幾何形狀是非常復雜的,如戰(zhàn)斗機、運輸機全機構型。要得到高精度的數(shù)值解,邊界條件處理本身應保證適當?shù)挠嬎憔取6谶吔缣幚碇校行┪锢砹渴峭ㄟ^插值方法求得的。插值的精度直接影響邊界條件處理的精度,為此一般要求邊界附近的網格線盡可能與邊界正交,而且在物面邊界附近還需保證一定的網格節(jié)點密度,過稀的網格將導致計算精度的降低。
由此可知,對于數(shù)值求解偏微分方程(PDE)的定解問題而言,網格分布是十分重要的。在達到相同解的精度的前提下,合理的網格分布往往可以大大減少網格點的數(shù)目,從而大大節(jié)省所需要的計算機內存和計算時間。計算經驗表明,在某些問題中,不合適的網格分布有可能導致計算過程的不穩(wěn)定或不收斂。
CFD對計算網格的基本要求
網格質量是網格生成技術重點關注的研究領域。就結構網格而言,網格質量一般包括網格的光滑性、正交性、分布合理性等。對于非結構網格而言,網格的光滑性和分布合理性也是需要關注的重要方面,雖然不存在所謂的網格正交性,但也一般要求網格的形狀要盡量“正規(guī)”,即盡可能為正二角形正四面體。對于黏性流動計算問題,如在邊界層內采用“純”非結構網格(各異性四面體),其計算精度和離散效率均有不足,此時一般在邊界層內采用結構(六面體)或半結構(三棱柱)網格,因此網格的法向正交性也是需要關注的問題。以下進行簡要的分析說明。
01
網格光滑性要求
就結構網格而言,通常計算空間中的均勻網格對應于物理空間中的非均勻網格。為了得到高精度的計算結果,我們要求物理空間中的網格不均勻的變化是逐漸過渡的,而不是突然過渡的。
對于非結構網格而言,如果采用有限體積法,雖然不必進行從物理空間到計算空間的網格變化,但是由于現(xiàn)有的計算方法一般采用多項式分布來近似流場物理量的分布,非均勻的網格分布將導致重構多項式分布的“奇異”性,從而影響計算過程的穩(wěn)定性和計算結果的精確度。因此,網格分布的光滑性也是一個重要的質量指標。
注:網格的光滑性一般由相鄰單元的特征尺度之比來表征。
02
網格正交性要求
網格正交性指多維情況時結構網格的線(面)要盡可能正交。將此概念推廣至半結構的三棱柱網格,即要求棱柱的上下三角形面與其他面盡可能正交。其原因是正交性差的歪斜網格往往會帶來較大的計算誤差。
注:在物面邊界附近,網格的正交使得邊界條件的處理更加直截了當和精確,從而有利于提高計算精度。同時,對于非線性問題,一般采用迭代方法求解,實踐表明,網格的正交性對收斂歷程有較大影響,一般采用正交性好的網格計算收斂性更好,收斂速度更快;對于含激波的流場模擬,一般要求網格線與激波面正交,在波系結構非常復雜的情況下,要做到全場網格與激波正交是很不現(xiàn)實的,此時可以重點關注主激波的形狀。
03
網格分布要求
如前所述,為了保證計算精度同時又節(jié)省計算機內存和計算時間,必須在流場參數(shù)變化大的區(qū)域分布較密的網格。在某些問題中,在求解流場之前就大致知道流場的哪些區(qū)域流動參數(shù)變化大,如機翼后緣、翼身結合部等幾何曲率變化較大的區(qū)域。對于這類問題,我們可以事先在該區(qū)域布置較密的網格。但是,大多數(shù)情況下,我們事先并不清楚流場的哪些區(qū)域流動參數(shù)會發(fā)生劇烈變化,如超聲速情況下流場中會出現(xiàn)激波,而激波導致物理量的強間斷,為了得到清晰的流場結構,就必須在激波區(qū)域加密網格。然而事先我們并不知道激波的準確位置,為此必須根據流場解的結果自動的布置網格,這就是所謂的網格自適應技術。
圖文|子暃魚
排版|子暃魚
工程師必備
- 項目客服
- 培訓客服
- 平臺客服
TOP




















