24小时热门版块排行榜    

查看: 1649  |  回复: 7

hihaoer

新虫 (小有名气)

[求助] 质量源项求助

一个分离提纯膜,简单的二维模型。这里为了简化定义进口只有一种气体氢气,膜有选择吸收氢气的能力,定义近壁温一个很薄的流体区域为源项,并画上网格。定义外界氢气的分压力为0。源项与压力有关,氢气不断减少,压力不断降低,每个边界网格的压力不一样,所以源项的值也不一样。
写了段源项的UDF,如果把源项改成个常数能算,但改成与压力有关的函数就算不动,老出错。在这困扰很久了,还希望大家帮忙看看,谢谢!
这是UDF:
#include "udf.h"
#define PRE_EXP 2.07e-4   /*指前因子*/
#define ACTIVE 17.7e+3  /*活化能*/
#define R 8.314          /*气体常数*/
#define T 623.15         /*膜工作温度*/
#define L 40e-6          /*膜工的厚度*/
DEFINE_SOURCE(permeate_source,c,t,ds,ean)
{
real source;
real P;           /*某点总压*/
P = C_P(c, t);           /*读取单元总压力*/
source = -(PRE_EXP * exp(-ACTIVE/(R*T)) * sqrt(P)) / L; /*膜的渗透流量*/
return source;
}

示意图.png

[ Last edited by hihaoer on 2012-9-9 at 10:44 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fangdian

木虫 (正式写手)

错误提示?

» 本帖已获得的红花(最新10朵)

2楼2012-09-09 14:26:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hihaoer

新虫 (小有名气)

送鲜花一朵
引用回帖:
2楼: Originally posted by fangdian at 2012-09-09 14:26:33
错误提示?

浮点错误。

错误提示.png

3楼2012-09-09 16:45:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hihaoer

新虫 (小有名气)

4楼2012-09-09 17:03:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

47425804

木虫 (正式写手)


C_P(c, t)获取的是表压,有可能为负值,不能开方。
需要再加上操作压力换算成绝压才行。

» 本帖已获得的红花(最新10朵)

5楼2012-09-09 19:32:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fangdian

木虫 (正式写手)

楼上说的有道理,你的case和data打不开,缺少:
Error: File PP.C doesn't exist!
Error Object: #f

» 本帖已获得的红花(最新10朵)

6楼2012-09-09 19:43:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hihaoer

新虫 (小有名气)

送鲜花一朵
引用回帖:
5楼: Originally posted by 47425804 at 2012-09-09 19:32:27
C_P(c, t)获取的是表压,有可能为负值,不能开方。
需要再加上操作压力换算成绝压才行。

确实是这样,我用的是速度进口,表压为负,前些天刚帮别人算了个多孔介质压降的,都是负数。。没想到。。加上操作压力就可以算了。太感谢了!
7楼2012-09-09 22:00:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hihaoer

新虫 (小有名气)

送鲜花一朵
引用回帖:
6楼: Originally posted by fangdian at 2012-09-09 19:43:13
楼上说的有道理,你的case和data打不开,缺少:
Error: File PP.C doesn't exist!
Error Object: #f

就是1L的UDF,谢谢你的回复。
8楼2012-09-09 22:00:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hihaoer 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 cqv1zKnzd7hr 2026-11-09 5/250 2026-11-11 07:08 by Z7qfCsdIiD1T
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 FChKX8meu1ds 2026-11-09 5/250 2026-11-11 04:29 by Z7qfCsdIiD1T
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 FChKX8meu1ds 2026-11-09 5/250 2026-11-11 04:23 by Z7qfCsdIiD1T
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 cqv1zKnzd7hr 2026-11-09 7/350 2026-11-11 04:21 by Z7qfCsdIiD1T
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 0TEMgxdT1eIw 2026-11-09 6/300 2026-11-11 04:10 by Z7qfCsdIiD1T
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 0TEMgxdT1eIw 2026-11-09 6/300 2026-11-11 04:07 by Z7qfCsdIiD1T
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 zSbvVX9JFl9J 2026-11-09 6/300 2026-11-11 03:59 by Z7qfCsdIiD1T
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 cqv1zKnzd7hr 2026-11-09 6/300 2026-11-11 01:21 by Z7qfCsdIiD1T
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 0TEMgxdT1eIw 2026-11-09 6/300 2026-11-11 01:13 by Z7qfCsdIiD1T
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 zSbvVX9JFl9J 2026-11-09 6/300 2026-11-11 01:02 by Z7qfCsdIiD1T
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 zSbvVX9JFl9J 2026-11-09 5/250 2026-11-11 00:59 by Z7qfCsdIiD1T
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 cqv1zKnzd7hr 2026-11-09 5/250 2026-11-10 20:45 by Z7qfCsdIiD1T
[考博] 27博士申请 +3 白白白ly 2026-11-09 3/150 2026-11-09 21:39 by szq0819
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +6 5sJFDKp2kwuI 2026-11-06 14/700 2026-11-09 11:52 by GwrQVVR6Oenx
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-08 23:01 by l48hFd8Nwae2
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 16:12 by 5sJFDKp2kwuI
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 ttH3f1oA6ace 2026-11-04 11/550 2026-11-06 16:10 by 5sJFDKp2kwuI
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 7/350 2026-11-06 02:13 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 02:04 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 hNpp4Y2CGvSC 2026-11-04 12/600 2026-11-06 01:53 by qL1cU46YIdK0
信息提示
请填处理意见