应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » 蒸发的程序出错，顺便请大神指导下C_YI的用法
41/1返回列表
查看: 1064  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

tjj_wei

铁虫 (初入文坛)

[求助] 蒸发的程序出错，顺便请大神指导下C_YI的用法

#include "udf.h"      //包括常规宏
#include "sg_mphase.h"  //包括体积分数
#include "mem.h"    //包括温度及组分质量分数
#include "sg.h"
#include "flow.h"


#define T_SAT  373  //定义蒸发温度100℃
#define LAT_HT 1.e3

DEFINE_SOURCE(liq_src,cell,sec_th,dS,eqn)
{
Thread *mix_th,*pri_th;   //定义计算区域指针
real m_dot_l;   //定义液相质量转移kg/(m2*s)
mix_th=THREAD_SUPER_THREAD(sec_th);
pri_th=THREAD_SUB_THREAD(mix_th,2);
if(C_T(cell,mix_th)>=T_SAT)
  {
    m_dot_l=-0.1*C_VOF(cell,sec_th)*C_R(cell,sec_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT;
    dS[eqn]=-0.1*C_VOF(cell,sec_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT;
  }
else
  {
     m_dot_l=0.1*C_YI(cell,pri_th,2)*C_VOF(cell,pri_th)*C_R(cell,pri_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT;
     dS[eqn]=0;
  }
  return  m_dot_l;
}

DEFINE_SOURCE(vap_src,cell,pri_th,dS,eqn)
{
Thread *mix_th,*sec_th;   //定义计算区域指针
real m_dot_v;   //定义气相质量转移kg/(m2*s)
mix_th=THREAD_SUPER_THREAD(pri_th);
sec_th=THREAD_SUB_THREAD(mix_th,3);
if(C_T(cell,mix_th)>=T_SAT)
  {
    m_dot_v=0.1*C_VOF(cell,sec_th)*C_R(cell,sec_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT;
    dS[eqn]=0;
  }
else
  {
     m_dot_v=-0.1*C_YI(cell,pri_th,2)*C_VOF(cell,pri_th)*C_R(cell,pri_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT;
     dS[eqn]=-0.1*C_YI(cell,pri_th,2)*C_R(cell,pri_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT;
  }
    return  m_dot_v;
}

DEFINE_SOURCE(enrg_src,cell,mix_th,dS,eqn)
{
Thread *pri_th,*sec_th;   //定义计算区域指针
real m_dot;   //定义液相质量转移kg/(m2*s)
pri_th=THREAD_SUPER_THREAD(mix_th,2);
sec_th=THREAD_SUB_THREAD(mix_th,3);
if(C_T(cell,mix_th)>=T_SAT)
  {
    m_dot=-0.1*C_VOF(cell,sec_th)*C_R(cell,sec_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT;
    dS[eqn]=-0.1*C_VOF(cell,sec_th)*C_R(cell,sec_th)/T_SAT;
  }
else
  {
     m_dot=0.1*C_YI(cell,pri_th,2)*C_VOF(cell,pri_th)*C_R(cell,pri_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT;
     dS[eqn]=0.1*C_YI(cell,pri_th,2)*C_VOF(cell,pri_th)*C_R(cell,pri_th)/T_SAT;
  }
  return  LAT_HT*m_dot;
}

DEFINE_SOURCE(h2o_src,cell,pri_th,dS,eqn)
{
Thread *mix_th,*sec_th;   //定义计算区域指针
real m_dot_vw;   //定义气相质量转移kg/(m2*s)
mix_th=THREAD_SUPER_THREAD(pri_th);
sec_th=THREAD_SUB_THREAD(mix_th,3);
if(C_T(cell,mix_th)>=T_SAT)
  {
    m_dot_vw=0.1*C_VOF(cell,sec_th)*C_R(cell,sec_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT;
    dS[eqn]=0;
  }
else
  {
     m_dot_vw=-0.1*C_YI(cell,pri_th,2)*C_VOF(cell,pri_th)*C_R(cell,pri_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT;
     dS[eqn]=-0.1*C_VOF(cell,pri_th)*C_R(cell,pri_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT;
  }
    return  m_dot_vw;
}
这是我的程序的内容，一初始化就出错，不知道是不是我的C_YI事业有问题，还有第一相是0还是查的ID号啊？请大家帮我看看
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-11-04 00:07:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小七工作室1

新虫 (著名写手)

【答案】应助回帖

感谢参与，应助指数 +1
这个是自带的程序  没有问题的
赞一下 回复此楼
2楼2013-11-04 10:23:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tjj_wei

铁虫 (初入文坛)
引用回帖:
2楼: Originally posted by 小七工作室1 at 2013-11-04 10:23:10
这个是自带的程序  没有问题的

额，这是我根据网上的程序改写的啊
怎么是自带的，，，，
赞一下 回复此楼
3楼2013-11-04 10:57:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

highway1

新虫 (初入文坛)
引用回帖:
3楼: Originally posted by tjj_wei at 2013-11-04 10:57:31
额，这是我根据网上的程序改写的啊
怎么是自带的，，，，...

fluent里面有这个模型，phase里面
赞一下 回复此楼
4楼2018-05-18 16:17:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 tjj_wei 的主题更新
41/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 312求调剂 +4 陌宸希 2026-03-16 5/250 2026-03-17 17:09 by ruiyingmiao
[考研] 085600材料与化工求调剂 +5 绪幸与子 2026-03-17 5/250 2026-03-17 16:40 by laoshidan
[考研] 085600材料与化工 +4 安全上岸！ 2026-03-16 4/200 2026-03-17 14:02 by 勇敢太监王公公
[考研] 302求调剂 +4 小贾同学123 2026-03-15 8/400 2026-03-17 10:33 by 小贾同学123
[考研] 机械专硕325，寻找调剂院校 +3 y9999 2026-03-15 5/250 2026-03-16 19:58 by y9999
[考研] 304求调剂 +5 素年祭语 2026-03-15 5/250 2026-03-16 17:00 by 我的船我的海
[考研] 311求调剂 +5 26研0 2026-03-15 5/250 2026-03-16 16:21 by a不易
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
[考研] 求老师收留调剂 +4 jiang姜66 2026-03-14 5/250 2026-03-15 20:11 by Winj1e
[考研] 22408总分284求调剂 +3 InAspic 2026-03-13 3/150 2026-03-15 11:10 by zhq0425
[考研] 080500，材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
[考研] 265求调剂 +4 威化饼07 2026-03-12 4/200 2026-03-14 17:23 by userper
[基金申请] 现在如何回避去年的某一个专家，不知道名字 +3 zk200107 2026-03-12 6/300 2026-03-14 17:13 by zk200107
[考研] 328求调剂 +3 5201314Lsy！ 2026-03-13 6/300 2026-03-14 15:31 by hyswxzs
[考研] 26调剂/材料/英一数二/总分289/已过A区线 +6 步川酷紫123 2026-03-13 6/300 2026-03-13 21:59 by 星空星月
[考研] 329求调剂 +3 miaodesi 2026-03-12 4/200 2026-03-13 20:53 by 18595523086
[硕博家园] 085600 260分求调剂 +3 天空还下雨么 2026-03-13 5/250 2026-03-13 18:46 by 天空还下雨么
[考研] 土木第一志愿276求调剂，科研和技能十分丰富，求新兴方向的导师收留 +3 土木小天才 2026-03-12 3/150 2026-03-13 15:01 by JourneyLucky
[考博] 福州大学杨黄浩课题组招收2026年专业学位博士研究生，2026.03.20截止 +3 Xiangyu_ou 2026-03-12 3/150 2026-03-13 09:36 by duanwu655
[考博] 26读博 +4 Rui135246 2026-03-12 10/500 2026-03-13 07:15 by gaobiao
信息提示
关闭
请填处理意见
关闭