24小时热门版块排行榜    

查看: 1398  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料与化工一志愿南昌大学327求调剂推荐 +5 Ncdx123456 2026-03-13 6/300 2026-03-15 23:39 by lovewei0727
[考研] 材料工程专硕274一志愿211求调剂 +5 薛云鹏 2026-03-15 5/250 2026-03-15 20:38 by Logic2024
[考博] 东华理工大学化材专业26届硕士博士申请 +6 zlingli 2026-03-13 6/300 2026-03-15 20:00 by ryzcf
[考研] 294求调剂 +3 Zys010410@ 2026-03-13 4/200 2026-03-15 10:59 by zhq0425
[考研] 本科南京大学一志愿川大药学327 +3 麦田耕者 2026-03-14 3/150 2026-03-14 20:04 by 外星文明
[基金申请] 现在如何回避去年的某一个专家,不知道名字 +3 zk200107 2026-03-12 6/300 2026-03-14 17:13 by zk200107
[教师之家] 焦虑 +5 水冰月月野兔 2026-03-13 7/350 2026-03-14 15:14 by 农药害害
[考研] 290求调剂 +4 @将就将就看 2026-03-10 8/400 2026-03-14 14:23 by 千千运气
[考研] 材料080500调剂求收留 +3 一颗meteor 2026-03-13 3/150 2026-03-14 10:54 by peike
[考研] 307求调剂 +7 超级伊昂大王 2026-03-10 7/350 2026-03-14 00:49 by JourneyLucky
[考研] 311求调剂 +8 zchqwer 2026-03-10 8/400 2026-03-14 00:01 by JourneyLucky
[考研] 337一志愿华南理工0805材料求调剂 +7 mysdl 2026-03-11 9/450 2026-03-13 22:43 by JourneyLucky
[考研] 材料工程调剂 +9 咪咪空空 2026-03-12 9/450 2026-03-13 22:05 by 星空星月
[考研] 26调剂/材料科学与工程/总分295/求收留 +9 2026调剂侠 2026-03-12 9/450 2026-03-13 20:46 by 18595523086
[考研] 293求调剂 +3 世界首富 2026-03-11 3/150 2026-03-13 16:27 by JourneyLucky
[考研] 求调剂 +7 18880831720 2026-03-11 7/350 2026-03-13 16:10 by JourneyLucky
[考研] 070303一志愿西北大学学硕310找调剂 +3 d如愿上岸 2026-03-12 5/250 2026-03-13 10:56 by houyaoxu
[考研] 268求调剂 +4 好运连绵不绝 2026-03-12 4/200 2026-03-13 10:45 by hyswxzs
[考博] 读博申请 +5 感dd 2026-03-10 7/350 2026-03-11 17:02 by QGZDSYS
[考研] 279求调剂 +3 莫xiao 2026-03-10 4/200 2026-03-11 08:06 by 斩魂滴兔子!
信息提示
请填处理意见