Mountisky

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 5.5
帖子: 12
在線: 4.3小時(shí)
蟲號(hào): 998910
注冊(cè): 2010-04-17
性別: GG
專業(yè): 化工機(jī)械/氫能源設(shè)備

[交流] 【求助】UDF氣體的密度已有6人參與

需要用到真實(shí)氣體的狀態(tài)，需要自定義氣體的密度，

#include "udf.h"
DEFINE_PROPERTY(cell_density,cell,thread)
{
real T=C_T(cell,thread);
real P=C_P(cell,thread);
real R=8.314472;
real Mh=0.002;
real a=1.9155e-6;
real mu_lam;
mu_lam=P*Mh/(R*T+a*P*R);
return mu_lam;
}
這是我寫的函數(shù)，可是這個(gè)用這個(gè)函數(shù)時(shí)總是出錯(cuò)，由于學(xué)學(xué)過C語(yǔ)言，不太懂。
有人說要循環(huán)，拜請(qǐng)高手指點(diǎn)下，看看這個(gè)函數(shù)哪里有問題。

回復(fù)此樓

» 猜你喜歡

281求調(diào)劑（0805）已經(jīng)有9人回復(fù)
293求調(diào)劑已經(jīng)有14人回復(fù)
0703化學(xué)336分求調(diào)劑已經(jīng)有6人回復(fù)
304求調(diào)劑已經(jīng)有9人回復(fù)
環(huán)境領(lǐng)域全國(guó)重點(diǎn)實(shí)驗(yàn)室招收博士1-2名已經(jīng)有3人回復(fù)
材料專碩306英一數(shù)二已經(jīng)有10人回復(fù)
301求調(diào)劑已經(jīng)有6人回復(fù)
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑已經(jīng)有7人回復(fù)
302求調(diào)劑已經(jīng)有6人回復(fù)
26博士申請(qǐng) 已經(jīng)有3人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

是誤差嗎？？？？【FLUENT的一個(gè)問題】已經(jīng)有17人回復(fù)
udf編程通過，但是在FLUENT中運(yùn)行出錯(cuò)，麻煩udf高手幫忙看看已經(jīng)有4人回復(fù)
DPM模型中顆粒設(shè)定的問題已經(jīng)有17人回復(fù)
曳力函數(shù)的的UDF 已經(jīng)有15人回復(fù)
fluent在氣固流場(chǎng)中模擬兩種不同密度固體顆粒的分離已經(jīng)有9人回復(fù)
UDF應(yīng)用之后的初始化錯(cuò)誤已經(jīng)有6人回復(fù)
fluent 中自定義密度的問題已經(jīng)有9人回復(fù)
udf 重力方向變化已經(jīng)有17人回復(fù)
【求助】udf讀取混合物密度失敗、及平均流速讀�。� 已經(jīng)有8人回復(fù)
【求助】單機(jī)多核fluent并行udf 已經(jīng)有13人回復(fù)
【求助】密度函數(shù)的UDF 已經(jīng)有6人回復(fù)
【求助】fluent計(jì)算中，udf編譯通過，初始化出錯(cuò) 已經(jīng)有11人回復(fù)
【求助】混合氣體密度計(jì)算已經(jīng)有7人回復(fù)

孩子，你還年輕，努力��！

1樓 2011-03-14 15:18:04

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

回帖支持 ( 顯示支持度最高的前 50 名 )

47425804

木蟲 (正式寫手)

仿真EPI: 1
應(yīng)助: 24 (小學(xué)生)
金幣: 631.1
散金: 505
紅花: 27
帖子: 807
在線: 181.2小時(shí)
蟲號(hào): 905802
注冊(cè): 2009-11-18
專業(yè): 化學(xué)反應(yīng)工程

★ ★
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流
xiegangmai(金幣+1): 鼓勵(lì)討論交流！ 2011-03-14 18:35:56

誒。。。不過還有個(gè)問題~~~~
fluent里要自定義密度的話，好像只能與溫度有關(guān)，是不能引入壓力的。

贊一下(3人)

回復(fù)此樓

5樓2011-03-14 16:48:43

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

47425804

木蟲 (正式寫手)

仿真EPI: 1
應(yīng)助: 24 (小學(xué)生)
金幣: 631.1
散金: 505
紅花: 27
帖子: 807
在線: 181.2小時(shí)
蟲號(hào): 905802
注冊(cè): 2009-11-18
專業(yè): 化學(xué)反應(yīng)工程

★ ★
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流
xiegangmai(金幣+1): 鼓勵(lì)討論交流！ 2011-03-15 12:00:10

引用回帖:

Originally posted by Mountisky at 2011-03-14 19:29:38:
氣體的密度是不可能不引入壓力的啊。。我的課題是氣體壓縮問題，壓力是最主要參數(shù)的

哦，我明白你的意思，我做的課題氣體壓力變化不大，所以是作為不可壓縮處理的，基于壓力的分離式求解。這種情況下我記得不可用壓力算密度，因?yàn)樵诿芏戎幸雺毫︘?xiàng)就導(dǎo)致守恒方程格式有變化，fluent無(wú)法求解。

你做可壓縮的流動(dòng)，估計(jì)情況有不一樣吧，我不大了解了，只是把我知道的情況都告訴你，呵呵~~~

贊一下(3人)

回復(fù)此樓

7樓2011-03-15 08:54:01

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

普通回帖

47425804

木蟲 (正式寫手)

仿真EPI: 1
應(yīng)助: 24 (小學(xué)生)
金幣: 631.1
散金: 505
紅花: 27
帖子: 807
在線: 181.2小時(shí)
蟲號(hào): 905802
注冊(cè): 2009-11-18
專業(yè): 化學(xué)反應(yīng)工程

★ ★ ★
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流
xiegangmai(金幣+2): 謝謝應(yīng)助！ 2011-03-14 18:35:20

C_P讀出來(lái)的是表壓或真空度，還要加上你設(shè)置的操作壓力。
PS：mu_lam這個(gè)一般用來(lái)表示層流粘度的變量名，這里你拿來(lái)表示密度似乎。。。
當(dāng)然，從程序上說是沒問題的，只是影響可讀性。密度習(xí)慣用rho來(lái)表示。

[ Last edited by 47425804 on 2011-3-14 at 15:51 ]

贊一下(2人)

回復(fù)此樓

2樓2011-03-14 15:49:57

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

Mountisky

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 5.5
帖子: 12
在線: 4.3小時(shí)
蟲號(hào): 998910
注冊(cè): 2010-04-17
性別: GG
專業(yè): 化工機(jī)械/氫能源設(shè)備

引用回帖:

Originally posted by 47425804 at 2011-03-14 15:49:57:
C_P讀出來(lái)的是表壓或真空度，還要加上你設(shè)置的操作壓力。
PS：mu_lam這個(gè)一般用來(lái)表示層流粘度的變量名，這里你拿來(lái)表示密度似乎。。。
當(dāng)然，從程序上說是沒問題的，只是影響可讀性。密度習(xí)慣用rho來(lái)表示。

[ ...

#include "udf.h"
DEFINE_PROPERTY(cell_density,cell,thread)
{
real T=C_T(cell,thread);
real P=C_P(cell,thread);
real R=8.314472;
real Mh=0.002;
real a=1.9155e-6;
real rho;
real P0=P+101325；
rho=P0*Mh/(R*T+a*P0*R);
return rho;
}

是這個(gè)意思嗎？

贊一下

回復(fù)此樓

孩子，你還年輕，努力啊！

3樓2011-03-14 16:20:13

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

47425804

木蟲 (正式寫手)

仿真EPI: 1
應(yīng)助: 24 (小學(xué)生)
金幣: 631.1
散金: 505
紅花: 27
帖子: 807
在線: 181.2小時(shí)
蟲號(hào): 905802
注冊(cè): 2009-11-18
專業(yè): 化學(xué)反應(yīng)工程

★ ★ ★
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流
xiegangmai(金幣+2): 謝謝應(yīng)助！ 2011-03-14 18:35:40

引用回帖:

Originally posted by Mountisky at 2011-03-14 16:20:13:
#include "udf.h"
DEFINE_PROPERTY(cell_density,cell,thread)
{
real T=C_T(cell,thread);
real P=C_P(cell,thread);
real R=8.314472;
real Mh=0.002;
real a=1.9155e-6;
real rho;
r ...

差不多，你先試試看還有問題沒~~
PS：注意標(biāo)點(diǎn)問題，別用中文標(biāo)點(diǎn)了。
#include "udf.h"
DEFINE_PROPERTY(cell_density,cell,thread)
{
real T=C_T(cell,thread);
real P=C_P(cell,thread);
real R=8.314472;
real Mh=0.002;
real a=1.9155e-6;
real rho;
real P0=P+101325； //這個(gè)標(biāo)點(diǎn)
rho=P0*Mh/(R*T+a*P0*R);
return rho;
}

贊一下(2人)

回復(fù)此樓

4樓2011-03-14 16:44:56

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

Mountisky

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 5.5
帖子: 12
在線: 4.3小時(shí)
蟲號(hào): 998910
注冊(cè): 2010-04-17
性別: GG
專業(yè): 化工機(jī)械/氫能源設(shè)備

引用回帖:

Originally posted by 47425804 at 2011-03-14 16:48:43:
誒。。。不過還有個(gè)問題~~~~
fluent里要自定義密度的話，好像只能與溫度有關(guān)，是不能引入壓力的。

氣體的密度是不可能不引入壓力的啊。。我的課題是氣體壓縮問題，壓力是最主要參數(shù)的

贊一下

回復(fù)此樓

孩子，你還年輕，努力��！

6樓2011-03-14 19:29:38

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

liufangkan

木蟲 (著名寫手)

仿真EPI: 1
應(yīng)助: 25 (小學(xué)生)
金幣: 1693.1
散金: 1
紅花: 11
沙發(fā): 8
帖子: 1375
在線: 94.8小時(shí)
蟲號(hào): 1132792
注冊(cè): 2010-10-27
性別: GG
專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域

我是來(lái)學(xué)習(xí)的

回復(fù)此樓

有一種生活叫奮斗,努力就有希望!

8樓2011-03-15 09:48:31

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

cow092

jjj

9樓2013-03-31 14:44:34

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

6863523

金蟲 (正式寫手)

應(yīng)助: 3 (幼兒園)
金幣: 2471.5
散金: 208
紅花: 7
帖子: 837
在線: 264.1小時(shí)
蟲號(hào): 1396096
注冊(cè): 2011-09-09
性別: GG
專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域

引用回帖:

7樓: Originally posted by 47425804 at 2011-03-15 08:54:01
哦，我明白你的意思，我做的課題氣體壓力變化不大，所以是作為不可壓縮處理的，基于壓力的分離式求解。這種情況下我記得不可用壓力算密度，因?yàn)樵诿芏戎幸雺毫︘?xiàng)就導(dǎo)致守恒方程格式有變化，fluent無(wú)法求解。

...

正解，大牛

回復(fù)此樓

循序漸進(jìn)

10樓2014-05-30 15:27:41

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考博] 環(huán)境領(lǐng)域全國(guó)重點(diǎn)實(shí)驗(yàn)室招收博士1-2名 +3	QGZDSYS 2026-03-13	3/150	2026-03-18 00:13 by PY-Hzb
[考研] 303求調(diào)劑 +3	睿08 2026-03-17	5/250	2026-03-17 22:01 by 睿08
[考研] 085601專碩，總分342求調(diào)劑，地區(qū)不限 +4	share_joy 2026-03-16	4/200	2026-03-17 21:32 by hmn_wj
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3	1喜春8 2026-03-17	5/250	2026-03-17 17:18 by ber川cool子
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6	薛云鵬 2026-03-15	6/300	2026-03-17 11:05 by 學(xué)員h26Tkc
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4	Jy? 2026-03-16	4/200	2026-03-17 11:02 by gaoqiong
[考研] 275求調(diào)劑 +4	太陽(yáng)花天天開心 2026-03-16	4/200	2026-03-17 10:53 by 功夫瘋狂
[考研] 274求調(diào)劑 +5	時(shí)間點(diǎn) 2026-03-13	5/250	2026-03-17 07:34 by 熱情沙漠
[考研] 環(huán)境工程調(diào)劑 +6	大可digkids 2026-03-16	6/300	2026-03-16 17:16 by barlinike
[考研] 318求調(diào)劑 +3	Yanyali 2026-03-15	3/150	2026-03-16 16:41 by houyaoxu
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3	自由煎餅果子 2026-03-16	3/150	2026-03-16 14:10 by 運(yùn)氣yunqi
[考研] 材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大本人本科雙一流 +4	自由的_飛翔 2026-03-13	5/250	2026-03-14 19:39 by hmn_wj
[考研] 308求調(diào)劑 +5	是Lupa啊 2026-03-11	5/250	2026-03-13 22:13 by JourneyLucky
[考研] 0856材料與化工301求調(diào)劑 +5	奕束光 2026-03-13	5/250	2026-03-13 22:00 by 星空星月
[考研] 304求調(diào)劑 +7	7712b 2026-03-13	7/350	2026-03-13 21:42 by peike
[考研] 工科，求調(diào)劑 +3	我887 2026-03-11	3/150	2026-03-13 21:39 by JourneyLucky
[碩博家園] 085600 260分求調(diào)劑 +3	天空還下雨么 2026-03-13	5/250	2026-03-13 18:46 by 天空還下雨么
[考研] 材料專碩350 求調(diào)劑 +4	王金科 2026-03-12	4/200	2026-03-13 16:02 by ruiyingmiao
[考研] 295求調(diào)劑 +3	小匕仔汁 2026-03-12	3/150	2026-03-13 15:17 by vgtyfty
[考研] 工科調(diào)劑 +4	Jiang191123！ 2026-03-11	4/200	2026-03-13 15:15 by Miko19