| 18 | 1/1 | 返回列表 |
| 查看: 2667 | 回復(fù): 17 | |||
[交流]
是誤差嗎????【FLUENT的一個(gè)問題】
|
|||
|
FLUENT的材料庫中,可看到air的密度為1.225kg/m3, 現(xiàn)在我在一個(gè)圓柱管模型中,采用air作為流體,初始化后,管內(nèi)全為air; 執(zhí)行以下UDF,應(yīng)可算出管內(nèi)air的平均密度,實(shí)際算出的值為1.2353kg/m3。 這兩個(gè)值不一樣,是誤差嗎?好象誤差挺大。。。。 采用雙精度測試了一下,結(jié)果為1.225kg/m3。這個(gè)數(shù)值倒是準(zhǔn)確。 但是用混合物作測試值又不對了,請見下面一樓: [ 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 ] |
» 搶金幣啦!回帖就可以得到:
+4/228
+1/81
+1/81
+1/35
+1/27
+1/25
+1/18
+1/15
+1/13
+1/8
+1/6
+1/6
+1/5
+1/5
+1/4
+1/4
+1/3
+1/2
+1/1
+1/1
新蟲 (職業(yè)作家)
|
FLUENT中,air的平均分子量為28.966,利用組分輸運(yùn)模型, 我自定義一種分子量為28.966的氣體與air按摩爾分?jǐn)?shù)為1:9混合,理應(yīng)密度也為1.225kg/m3。采用雙精度,按上述UDF測試其管內(nèi)平均密度值為1.17667,UDM記錄的管內(nèi)每一個(gè)CELL的密度值都為1.176673870E+00; 與1.225相差挺有點(diǎn)大。! 再自定義一種分子量為30的氣體與air按摩爾分?jǐn)?shù)為1:9混合,理應(yīng)密度>1.225kg/m3。但按上述UDF測試其管內(nèi)平均密度值為1.18087,UDM記錄的管內(nèi)每一個(gè)CELL的密度值都為1.180874250E+00 ; 與1.225相差也很大!! 有誰能告訴,上述差值確實(shí)是由誤差導(dǎo)致嗎???還是我的測試有問題???上述UDF有問題嗎??? [ Last edited by kylafree on 2012-2-29 at 10:14 ] |
|
一個(gè)大氣壓按101325算的,在FLUENT材料面板中可以看到,air雖然是混合物,但FLUENT給出了它的平均分子量為28.966,密度為1.225kg/m3。 利用理想氣體狀態(tài)方程PV=nRT計(jì)算的話確實(shí)28.966的分子量與1.225kg/m3的密度剛好對應(yīng)上。 P=101325 , R=8.31451(理想氣體常數(shù)) , T=273.15+15 (15攝氏度) 取n=1 mol ,可算出此時(shí)的氣體摩爾體積為V=23.64L; 現(xiàn)在取1m3=1000L空氣,它的質(zhì)量為(1000/23.64)*28.966=1.225kg,所以空氣的密度為1.225kg/m3,對應(yīng)的這么好! 但又有一點(diǎn)兒問題,我又算了一下其它幾種氣體,C2H2,C2H6,NH3,由分子量計(jì)算的密度與FLUENT數(shù)據(jù)庫給的密度近似相等,有一定誤差;而H2O蒸汽差的比較遠(yuǎn)(這是不是由于15度時(shí)水的飽和蒸汽壓比較低的緣故?) [ 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 | 好運(yùn)連綿不絕 2026-03-12 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 0703化學(xué)336分求調(diào)劑 +4 | zbzihdhd 2026-03-15 | 5/250 |
|
|
[考研] 材料專碩326求調(diào)劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
[考研] 283求調(diào)劑 +3 | 聽風(fēng)就是雨; 2026-03-16 | 3/150 |
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
[考研] 0854控制工程 359求調(diào)劑 可跨專業(yè) +3 | 626776879 2026-03-14 | 9/450 |
|
|
[考研] 一志愿985,本科211,0817化學(xué)工程與技術(shù)319求調(diào)劑 +5 | Liwangman 2026-03-15 | 5/250 |
|
|
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +5 | d如愿上岸 2026-03-12 | 8/400 |
|
|
[考研] 材料與化工一志愿南昌大學(xué)327求調(diào)劑推薦 +7 | Ncdx123456 2026-03-13 | 8/400 |
|
|
[基金申請] 現(xiàn)在如何回避去年的某一個(gè)專家,不知道名字 +3 | zk200107 2026-03-12 | 6/300 |
|
|
[考研] 復(fù)試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 5/250 |
|
|
[考研] 297求調(diào)劑 +4 | 學(xué)海漂泊 2026-03-13 | 4/200 |
|
|
[考研] 266求調(diào)劑 +4 | 學(xué)員97LZgn 2026-03-13 | 4/200 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +6 | 邱gl 2026-03-12 | 7/350 |
|
|
[考研] 一志愿西南交大,材料專碩317求調(diào)劑 +5 | lx8568 2026-03-11 | 5/250 |
|
|
[考研] 0703化學(xué)一志愿211 總分320求調(diào)劑 +5 | 瑪卡巴卡啊哈 2026-03-11 | 5/250 |
|
|
[考研] 求調(diào)劑 +7 | 18880831720 2026-03-11 | 7/350 |
|
|
[考研] 314求調(diào)劑 +7 | 無懈可擊的巨人 2026-03-12 | 7/350 |
|
|
[考研] 274求調(diào)劑 +3 | S.H1 2026-03-12 | 3/150 |
|
|
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3 | d如愿上岸 2026-03-13 | 3/150 |
|