24小时热门版块排行榜    

查看: 2478  |  回复: 10
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

Mountisky

新虫 (初入文坛)

[交流] 【求助】UDF气体的密度 已有6人参与

需要用到真实气体的状态,需要自定义气体的密度,


#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;
}
这是我写的函数,可是这个用这个函数时总是出错,由于学学过C语言,不太懂。
有人说要循环,拜请高手指点下,看看这个函数哪里有问题。
回复此楼
孩子,你还年轻,努力啊!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

47425804

木虫 (正式写手)


★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
xiegangmai(金币+1): 鼓励讨论交流! 2011-03-15 12:00:10
引用回帖:
Originally posted by Mountisky at 2011-03-14 19:29:38:
气体的密度是不可能不引入压力的啊。。我的课题是气体压缩问题,压力是最主要参数的

哦,我明白你的意思,我做的课题气体压力变化不大,所以是作为不可压缩处理的,基于压力的分离式求解。这种情况下我记得不可用压力算密度,因为在密度中引入压力项就导致守恒方程格式有变化,fluent无法求解。

你做可压缩的流动,估计情况有不一样吧,我不大了解了,只是把我知道的情况都告诉你,呵呵~~~
7楼2011-03-15 08:54:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

47425804

木虫 (正式写手)


★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
xiegangmai(金币+2): 谢谢应助! 2011-03-14 18:35:20
C_P读出来的是表压或真空度,还要加上你设置的操作压力。
PS:mu_lam这个一般用来表示层流粘度的变量名,这里你拿来表示密度似乎。。。
当然,从程序上说是没问题的,只是影响可读性。密度习惯用rho来表示。

[ Last edited by 47425804 on 2011-3-14 at 15:51 ]
2楼2011-03-14 15:49:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Mountisky

新虫 (初入文坛)

引用回帖:
Originally posted by 47425804 at 2011-03-14 15:49:57:
C_P读出来的是表压或真空度,还要加上你设置的操作压力。
PS:mu_lam这个一般用来表示层流粘度的变量名,这里你拿来表示密度似乎。。。
当然,从程序上说是没问题的,只是影响可读性。密度习惯用rho来表示。

[ ...

#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;
}

是这个意思吗?
孩子,你还年轻,努力啊!
3楼2011-03-14 16:20:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

47425804

木虫 (正式写手)


★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
xiegangmai(金币+2): 谢谢应助! 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:注意标点问题,别用中文标点了。
#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;
}
4楼2011-03-14 16:44:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 298求调剂 +8 上岸6666@ 2026-03-20 8/400 2026-03-23 11:02 by laoshidan
[考研] 求调剂材料学硕080500,总分289分 5+3 @taotao 2026-03-19 21/1050 2026-03-23 10:17 by 冠c哥
[考研] 306求调剂 +7 chuanzhu川烛 2026-03-18 7/350 2026-03-23 09:22 by king123!
[考研] 328求调剂,英语六级551,有科研经历 +6 生物工程调剂 2026-03-17 10/500 2026-03-22 20:22 by edmund7
[考研] 一志愿上海交大生物与医药专硕324分,求调剂 +3 jiajunX 2026-03-22 3/150 2026-03-22 19:32 by brblmd
[考研] 材料 271求调剂 +5 展信悦_ 2026-03-21 5/250 2026-03-21 17:29 by 学员8dgXkO
[考研] 266求调剂 +3 哇呼哼呼哼 2026-03-20 3/150 2026-03-21 16:46 by barlinike
[考研] 材料与化工(0856)304求 B区 调剂 +3 邱gl 2026-03-21 3/150 2026-03-21 13:47 by lature00
[考研] 295求调剂 +4 一志愿京区211 2026-03-18 6/300 2026-03-20 23:41 by JourneyLucky
[考研] 294求调剂材料与化工专硕 +15 陌の森林 2026-03-18 15/750 2026-03-20 23:28 by JourneyLucky
[考研] 304求调剂 +7 司空. 2026-03-18 7/350 2026-03-20 23:08 by JourneyLucky
[考研] 一志愿南京理工大学085701资源与环境302分求调剂 +4 葵梓卫队 2026-03-18 6/300 2026-03-20 23:02 by JourneyLucky
[考研] 求调剂,一志愿:南京航空航天大学大学 ,080500材料科学与工程学硕,总分289分 +4 @taotao 2026-03-19 4/200 2026-03-20 22:14 by JourneyLucky
[考研] 290求调剂 +7 ^O^乜 2026-03-19 7/350 2026-03-20 21:43 by JourneyLucky
[考研] A区线材料学调剂 +5 周周无极 2026-03-20 5/250 2026-03-20 21:33 by laoshidan
[考研] 求调剂 +3 @taotao 2026-03-20 3/150 2026-03-20 19:35 by JourneyLucky
[考研] 261求B区调剂,科研经历丰富 +3 牛奶很忙 2026-03-20 4/200 2026-03-20 19:34 by JourneyLucky
[考研] 工科材料085601 279求调剂 +7 困于星晨 2026-03-17 9/450 2026-03-20 17:38 by 无懈可击111
[考研] 0856调剂,是学校就去 +8 sllhht 2026-03-19 9/450 2026-03-20 14:25 by 无懈可击111
[考研] 收复试调剂生 +4 雨后秋荷 2026-03-18 4/200 2026-03-18 14:16 by elevennnne
信息提示
请填处理意见