24小时热门版块排行榜    

查看: 2437  |  回复: 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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 281求调剂 +4 2026计算机_诚心 2026-03-01 4/200 2026-03-01 15:21 by chemisry
[考研] 304求调剂 +6 曼殊2266 2026-02-28 7/350 2026-03-01 15:14 by wjLi2017
[考研] 课题组接收材料类调剂研究生 +3 gaoxiaoniuma 2026-02-28 4/200 2026-03-01 14:30 by jjj三跨
[考研] 298求调剂 +9 人间唯你是清欢 2026-02-28 12/600 2026-03-01 14:23 by Ducount.Y
[考研] 材料学调剂 +8 提神豆沙包 2026-02-28 10/500 2026-03-01 14:14 by peike
[考研] 材料化工调剂 +8 今夏不夏 2026-03-01 9/450 2026-03-01 13:25 by yc258
[考研] 290求调剂 +8 材料专硕调剂; 2026-02-28 9/450 2026-03-01 12:46 by 闭眼看蓝天
[考研] 321求调剂一志愿东北林业大学材料与化工英二数二 +3 虫虫虫虫虫7 2026-03-01 6/300 2026-03-01 11:50 by gaoxiaoniuma
[考研] 调剂 +3 简木ChuFront 2026-02-28 3/150 2026-03-01 11:46 by 王伟要上岸啊
[考研] 高分子化学与物理调剂 +5 好好好1233 2026-02-28 9/450 2026-03-01 10:59 by fengyu211
[硕博家园] 2025届双非化工硕士毕业,申博 +3 更多的是 2026-02-27 4/200 2026-03-01 10:04 by ztg729
[考研] 272求调剂 +4 田智友 2026-02-28 4/200 2026-03-01 06:43 by 刘兵
[考研] 285求调剂 +6 满头大汗的学生 2026-02-28 6/300 2026-03-01 06:29 by Trying]
[考研] 材料调剂 +4 爱擦汗的可乐冰 2026-02-28 4/200 2026-03-01 00:38 by 猫猫球alter
[基金申请] 面上模板改不了页边距吧? +5 ieewxg 2026-02-25 6/300 2026-03-01 00:10 by addressing
[考研] 307求调剂 +4 73372112 2026-02-28 6/300 2026-03-01 00:04 by ll247
[考研] 304求调剂 +3 52hz~~ 2026-02-28 5/250 2026-03-01 00:00 by 52hz~~
[考博] 26申博 +4 想申博! 2026-02-26 4/200 2026-02-28 21:37 by limorning
[考研] 085600材料工程一志愿中科大总分312求调剂 +8 吃宵夜1 2026-02-28 10/500 2026-02-28 20:27 by L135790
[考研] 0856调剂 +3 刘梦微 2026-02-28 3/150 2026-02-28 13:22 by houyaoxu
信息提示
请填处理意见