小白最近逛圖書館,發現最近關於Fluent的書是越來越多了,而且還發現這些關於Fluent教材中的案例都大同小異。小白接受小牛師兄的建議,找了一本結構比較鮮明的書照着上面的案例就練了起來。不過當練習的案例多了,小白卻發現了一個問題“在建立計算模型時,為什么有時候采用二維模型,而有時候又采用三維模型?二維模型和三維模型的差異在哪里?二維代替三維到底有沒有問題?”帶着這些問題,小白找到了小牛師兄。
“師兄,你說Fluent中計算的二維模型到底有沒有意義呢?”小白問。
“怎么會沒有意義呢,當然有意義啊,當你計算機性能不行的時候,你就會知道能算二維是一件多么有意義的事情了。”小牛師兄說。
“可是,現實世界中不可能存在二維的物體啊。”小白說。
“嗯,這倒是真的,現實世界的任何物體都是三維的,不存在沒有厚度的面,也不存在沒有截面積的線,更不存在沒有體積的點,所有的這一切都是近似的。但是近似並不表示一定就不行,實際上在很多時候利用這些近似,可以獲得不錯的結果。”小牛師兄說。
“那什么時候可以利用這種簡化呢?比如說用二維代替三維模型。”小白問。
“我們常說的以二維模型來近似模擬三維模型,只是針對某一維尺寸遠小於另外兩維的情況下。這時候需要注意的是,采用二維模型時無法考慮法向方向的物理信息,在進行二維幾何模型計算過程中,常常將面法向方向簡化為常數,即總是認為面幾何體的厚度為定值(通常默認為1 m)。在CFD中,二維模型一般都為平面模型。”小牛師兄說。
“我在做案例的時候,碰到過很奇怪的問題。比如說計算一個管道內部流動,有的案例采用三維幾何模型,如圖1所示;有的案例用半個三維幾何,如圖2所示;有的案例采用四分之一模型,如圖3所示;有的案例采用管道縱剖面幾何(二維矩形),如圖四所示;而有的案例又采用二維的半剖面幾何(矩形),如圖5所示。這個到底該怎么選呢?”小白問。
“通常來說,采用全三維模型自然是不存在任何問題的,對於比較簡單的幾何模型強烈建議使用全三維模型,因為其能提供所有的信息。在流場對稱的情況下可利用對稱性采用半模型或四分之一模型,當然嚴格意義上的流場對稱是不存在的,不過有時候這種簡化是可以認可的。比如圖2所示的模型,如果沿對稱面法向方向沒有外部力作用則可以使用,但如果有外部力(如重力等)沿法向方向,則不可使用。圖3也是一樣,對稱面法向法向不允許有外部力存在。至於圖4和圖5的簡化,則除了要求外部力方向不能沿面法向之外,還必須假設流場特性在任意方向截面分布一致。”小牛師兄說。
案例測試
分別測試以上五種模型的差異。
計算模型為圖1所示的直管道,截面直徑10mm,長度100mm。入口面流速 0.1 m/s,出口采用靜壓 0 Pa。流體介質為水。雷諾數約為1000,采用層流計算。
計算模型分別采用上述的五種模型。
計算結果
以下采用相同計算條件得到的速度雲圖分布(從左至右分別對應上圖1-5模型)。
分析:從圖中看出,五種不同模型計算得到的速度分布趨勢基本保持一致,但是利用模型4計算得到的最大速度明顯小於其他四種模型計算得到的結果。
根據理論計算結果,充分發展的核心最大速度應為平均速度的2倍,即 0.2 m/s。
改進計算
將計算區域延長,從當前的0.1m延長至0.4m重新計算。從上圖可以看出,采用軸對稱方式計算結果能夠與3D模型相吻合,因此改進計算只采用兩種平面模型。
下圖為延長至0.4m后的兩種模型計算結果,可以看到最大速度已經非常接近理論值 0.2m/s。而2D平面模型計算的最大值小於0.15m/s,距離理論值相去甚遠。
總結:對於類似直管流動問題,采用全3D模型,半3D模型、四分之一3D模型以及軸對稱二維模型都可接受,但是2D平面模型是不可接受的。
原因分析:2D平面模型實際上是當做拉伸體而不是旋轉體計算,拉伸厚度默認為1m。在計算過程中對於平面法向方向的考慮與真實的圓柱體不相符。
其他問題
實際上很多的模型都存在這樣的限制,在使用2D平面模型時需要格外注意,其只適用於厚度相對於其他兩個維度非常小的拉伸體的簡化,對於回轉體則不適用。
如下圖中的彎管模型及彎道模型。