[請益] 使用Fluent中k-ε模型計算紊流場

看板Mechanical (機械)作者 (who)時間13年前 (2012/06/20 16:17), 編輯推噓5(5013)
留言18則, 5人參與, 最新討論串1/1
小妹初次在這裡發文 請多指教 這是我計算的流場幾何http://imgur.com/thRWF
長25 高14 中央方塊為障礙物所在,邊長為1 流體為液態水 上邊界與下邊界均為no slip wall 左邊界為velocity-inlet 右邊界則是pressure outlet 雷諾數大小約為4萬 以這樣的條件,使用Fluent中的Laminar,以unsteady,PISO算則計算 可得出障礙物後方有渦流產生 但使用Fluent中的 標準k-εmodel,也以unsteady與PISO算則計算 就是沒辦法算出障礙物後方渦流(無法產生) 我的詳細計算設置: http://imgur.com/emwsd
http://imgur.com/SVcaX
http://imgur.com/13Usy
http://imgur.com/ZiUWw
不知道為何這樣的設置下,標準k-εmodel算不出障礙物後方渦流呢? 是否我有什麼地方設定有缺失或錯誤? 因事態緊急,要趕在7月初計算出來~"~ 就懇請版上各位高手多多賜教了 謝謝 -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 124.11.132.222 choose1:轉錄至看板 Master_D 06/20 16:22 choose1:轉錄至看板 ask 06/20 16:22

06/21 04:31, , 1F
你time step size調多少?
06/21 04:31, 1F

06/21 04:32, , 2F
障礙物後方的網格大小多少? 有考慮加密嗎
06/21 04:32, 2F
回k大:我的time step是0.05sec(不同model但計算條件均一致) 而且因為用standard k-εmodel計算的結果需要與另一個model相比較 因此網格是不能更動的@@

06/21 06:53, , 3F
為什麼第一行要加小妹這兩個字 (.▁.?)
06/21 06:53, 3F

06/21 06:53, , 4F
建議下次要詢問像這類的問題 要另外再附上妳計算後的
06/21 06:53, 4F

06/21 06:53, , 5F
結果 例如:壓力線圖、速度線圖......等
06/21 06:53, 5F

06/21 06:54, , 6F
看了一下 編寫的UDF正確嗎?我怎麼看到上面有錯誤顯示
06/21 06:54, 6F
回p大:因為我PO文的電腦不是跑模擬的那台電腦 所以UDF沒有放在資料夾裡~"~讀取時case&data會因找不到UDF而出現error http://imgur.com/JKI4W
0.05sec時 http://imgur.com/sJUUz
80.45sec時 由上面兩張stream function圖 可以看出用standard k-εmodel計算的流場根本沒變化orz PS.文章開頭只是想到這樣寫就這樣寫了~"~

06/21 09:04, , 7F
長,寬,邊長單位?Re值反推inlet velocity看有沒設錯...
06/21 09:04, 7F
回t大:長寬邊長單位目是m(公尺),velocity inlet確定是OK的 因為指導教授堅持要使用standard k-ε model算出渦流,所以仍需釐清此問題 我已經有嘗試過將使用其他model在相同設定下算出的渦流當成initial condition 再輸入於k-ε model中計算 結果原先的vortex消散,也沒有再產生新的vortex http://imgur.com/ZyfeL
(所有計算每個time step都是設為0.05sec) 在這些前述的條件之下,真的沒辦法以standard k-ε model計算出渦流現象了嗎? 但若是將提高雷諾數至40萬的話,standard k-ε model能不能算出障礙物後方的渦流呢? 就再請版上各位神人解惑了,謝謝 ※ 編輯: choose1 來自: 140.121.151.48 (06/21 14:40)

06/22 20:06, , 8F
婀~你這應該是要解一個 von Kármán vortex street吧
06/22 20:06, 8F

06/22 20:07, , 9F
你的結構網格可能還不夠密去抓到渦流 建議採用非結構
06/22 20:07, 9F

06/22 20:07, , 10F
你的雷諾數好像有算錯建議你在檢查一下
06/22 20:07, 10F

06/22 20:11, , 11F
von Kármán vortex 產生範圍 250 < Re < 2*10^5
06/22 20:11, 11F

06/22 20:13, , 12F
之前有跑過類似題目 搭配standard k-εmodel Re>400
06/22 20:13, 12F

06/22 20:14, , 13F
非結構網格 有抓到渦流 可以參考看看
06/22 20:14, 13F

06/23 01:12, , 14F
我覺得 0.05s 太大
06/23 01:12, 14F

06/23 01:13, , 15F
你沒有用 NITA 的話 那你的收斂標準是什麼
06/23 01:13, 15F

06/23 01:14, , 16F
還是迭代 20 次不管收斂標準就跑下一步了
06/23 01:14, 16F

06/24 03:50, , 17F
你先用個0.001s下去跑看看 不然就0.0001s
06/24 03:50, 17F

06/24 03:51, , 18F
紊流模式下 你單一網格如果比渦旋還大 那也會算不出來
06/24 03:51, 18F
文章代碼(AID): #1FuOR_lQ (Mechanical)
文章代碼(AID): #1FuOR_lQ (Mechanical)