| 18 | 1/1 | 返回列表 |
| 查看: 2660 | 回復: 17 | |||
[交流]
是誤差嗎????【FLUENT的一個問題】
|
|||
|
FLUENT的材料庫中,可看到air的密度為1.225kg/m3, 現(xiàn)在我在一個圓柱管模型中,采用air作為流體,初始化后,管內(nèi)全為air; 執(zhí)行以下UDF,應可算出管內(nèi)air的平均密度,實際算出的值為1.2353kg/m3。 這兩個值不一樣,是誤差嗎?好象誤差挺大。。。! 采用雙精度測試了一下,結(jié)果為1.225kg/m3。這個數(shù)值倒是準確。 但是用混合物作測試值又不對了,請見下面一樓: [ Last edited by kylafree on 2012-2-28 at 19:56 ] #include "udf.h" DEFINE_ON_DEMAND(st2) { cell_t c; Thread *t = Lookup_Thread(Get_Domain(1), 8); real s=0.0; int n=0; begin_c_loop(c,t) { n=n+1; s=s+C_R(c,t); C_UDMI(c,t,0)=C_R(c,t); } end_c_loop(c,t) s=s/n; Message("\n s n =%g %d\n", s,n); } [ Last edited by kylafree on 2012-2-29 at 10:16 ] |
» 搶金幣啦!回帖就可以得到:
+1/85
+1/85
+1/79
+2/58
+1/48
+1/38
+1/34
+2/30
+1/28
+1/24
+1/18
+1/15
+1/12
+1/6
+1/5
+1/2
+1/1
+1/1
+1/1
+1/1
|
FLUENT中,air的平均分子量為28.966,利用組分輸運模型, 我自定義一種分子量為28.966的氣體與air按摩爾分數(shù)為1:9混合,理應密度也為1.225kg/m3。采用雙精度,按上述UDF測試其管內(nèi)平均密度值為1.17667,UDM記錄的管內(nèi)每一個CELL的密度值都為1.176673870E+00; 與1.225相差挺有點大。! 再自定義一種分子量為30的氣體與air按摩爾分數(shù)為1:9混合,理應密度>1.225kg/m3。但按上述UDF測試其管內(nèi)平均密度值為1.18087,UDM記錄的管內(nèi)每一個CELL的密度值都為1.180874250E+00 ; 與1.225相差也很大。! 有誰能告訴,上述差值確實是由誤差導致嗎???還是我的測試有問題???上述UDF有問題嗎??? [ Last edited by kylafree on 2012-2-29 at 10:14 ] |
|
一個大氣壓按101325算的,在FLUENT材料面板中可以看到,air雖然是混合物,但FLUENT給出了它的平均分子量為28.966,密度為1.225kg/m3。 利用理想氣體狀態(tài)方程PV=nRT計算的話確實28.966的分子量與1.225kg/m3的密度剛好對應上。 P=101325 , R=8.31451(理想氣體常數(shù)) , T=273.15+15 (15攝氏度) 取n=1 mol ,可算出此時的氣體摩爾體積為V=23.64L; 現(xiàn)在取1m3=1000L空氣,它的質(zhì)量為(1000/23.64)*28.966=1.225kg,所以空氣的密度為1.225kg/m3,對應的這么好! 但又有一點兒問題,我又算了一下其它幾種氣體,C2H2,C2H6,NH3,由分子量計算的密度與FLUENT數(shù)據(jù)庫給的密度近似相等,有一定誤差;而H2O蒸汽差的比較遠(這是不是由于15度時水的飽和蒸汽壓比較低的緣故?) [ Last edited by kylafree on 2012-2-29 at 16:52 ] |
|
按你說的求平均密度?諝夥肿恿繛28.966,自定義一種分子量為28.966的氣體,作測試,代碼如下,單精度,直接對密度平均,值為1.18422,密度*體積求和后再除總體積,值為1.18864。還是與1.225相差很多。 #include "udf.h" DEFINE_ON_DEMAND(st2) { cell_t c; Thread *t = Lookup_Thread(Get_Domain(1), 8); real m,s=0.0,q=0.0,r=0.0; int n=0; begin_c_loop(c,t) { n=n+1; s=s+C_R(c,t); C_UDMI(c,t,0)=C_R(c,t); r=r+C_VOLUME(c,t); q=q+C_VOLUME(c,t)*C_R(c,t); } end_c_loop(c,t) s=s/n; m=q/r; Message("\n s m n =%g %g %d\n", s,m,n); } |
| 18 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 268求調(diào)劑 +7 | 一定有學上- 2026-03-14 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 289求調(diào)劑 +6 | 步川酷紫123 2026-03-11 | 6/300 |
|
|
[考研] 285化工學碩求調(diào)劑(081700) +9 | 柴郡貓_ 2026-03-12 | 9/450 |
|
|
[考研] 材料專碩274一志愿陜西師范大學求調(diào)劑 +5 | 薛云鵬 2026-03-13 | 5/250 |
|
|
[考研] 機械專碩325,尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
[考研] 321求調(diào)劑 +5 | 大米飯! 2026-03-15 | 5/250 |
|
|
[考研] 085600調(diào)劑 +5 | 漾漾123sun 2026-03-12 | 6/300 |
|
|
[教師之家] 焦慮 +7 | 水冰月月野兔 2026-03-13 | 9/450 |
|
|
[考研] 326求調(diào)劑 +3 | mlpqaz03 2026-03-15 | 3/150 |
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|
|
[考研] 材料080500調(diào)劑求收留 +3 | 一顆meteor 2026-03-13 | 3/150 |
|
|
[考研] 337一志愿華南理工0805材料求調(diào)劑 +7 | mysdl 2026-03-11 | 9/450 |
|
|
[考研] 0703,333分求調(diào)劑 一志愿鄭州大學-物理化學 +3 | 李魔女斗篷 2026-03-11 | 3/150 |
|
|
[考研] 0856材料與化工301求調(diào)劑 +5 | 奕束光 2026-03-13 | 5/250 |
|
|
[考研] 四川大學085601材料工程專碩 初試294求調(diào)劑 +4 | 祝我們好在冬天 2026-03-11 | 4/200 |
|
|
[考研] 0703化學求調(diào)劑 +7 | 綠豆芹菜湯 2026-03-12 | 7/350 |
|
|
[考研] 290求調(diào)劑 +7 | ADT 2026-03-12 | 7/350 |
|
|
[考研] 工科調(diào)劑 +4 | Jiang191123! 2026-03-11 | 4/200 |
|
|
[考研] 化工學碩306求調(diào)劑 +9 | 42838695 2026-03-12 | 9/450 |
|
|
[考研] 333求調(diào)劑 +3 | 152697 2026-03-12 | 4/200 |
|