24小时热门版块排行榜    

查看: 1399  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 290求调剂 +3 孔志浩 2026-03-12 8/400 2026-03-15 15:30 by 孔志浩
[考研] 085600材料与化工 求调剂 +11 enenenhui 2026-03-13 12/600 2026-03-14 21:48 by laoshidan
[基金申请] 现在如何回避去年的某一个专家,不知道名字 +3 zk200107 2026-03-12 6/300 2026-03-14 17:13 by zk200107
[考研] 330求调剂 +3 ?酱给调剂跪了 2026-03-13 3/150 2026-03-14 10:13 by JourneyLucky
[考研] 085600材料与化工 326 求调剂 +5 热爱生活ing 2026-03-09 5/250 2026-03-14 02:39 by JourneyLucky
[考研] 085600求调剂 +3 a邵星池 2026-03-09 3/150 2026-03-14 01:32 by JourneyLucky
[考研] 295复试调剂 +5 简木ChuFront 2026-03-09 5/250 2026-03-14 01:29 by JourneyLucky
[考研] 285 求调剂 资源与环境 一志愿北京化工大学 +3 未名考生 2026-03-10 3/150 2026-03-13 23:04 by JourneyLucky
[考研] 求材料调剂 +5 隔壁陈先生 2026-03-12 5/250 2026-03-13 22:03 by 星空星月
[考研] 26调剂/材料科学与工程/总分295/求收留 +9 2026调剂侠 2026-03-12 9/450 2026-03-13 20:46 by 18595523086
[考研] 311求调剂 +3 冬十三 2026-03-13 3/150 2026-03-13 20:41 by JourneyLucky
[硕博家园] 085600 260分求调剂 +3 天空还下雨么 2026-03-13 5/250 2026-03-13 18:46 by 天空还下雨么
[考研] 0703化学求调剂 +7 绿豆芹菜汤 2026-03-12 7/350 2026-03-13 17:25 by njzyff
[考研] 293求调剂 +3 世界首富 2026-03-11 3/150 2026-03-13 16:27 by JourneyLucky
[考研] 307求调剂 +5 超级伊昂大王 2026-03-12 5/250 2026-03-13 15:56 by 棒棒球手
[考研] 308求调剂 +3 是Lupa啊 2026-03-12 3/150 2026-03-13 14:30 by 求调剂zz
[论文投稿] 投稿问题 5+4 星光灿烂xt 2026-03-12 6/300 2026-03-13 14:17 by god_tian
[考研] 26考研求调剂 +5 丶宏Sir 2026-03-13 5/250 2026-03-13 13:05 by JourneyLucky
[考研] 08食品或轻工求调剂,本科发表3篇sci一区top论文,一志愿南师大食品科学与工程 +3 我是一个兵, 2026-03-10 3/150 2026-03-13 10:21 by Yuyi.
[考研] 327分求调剂086 +4 西红柿?小帅 2026-03-09 7/350 2026-03-10 14:47 by ruiyingmiao
信息提示
请填处理意见