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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » 水沸腾UDF总是报错
111/2返回列表12下一页
查看: 2788  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

fengyezhi

新虫 (小有名气)

[求助] 水沸腾UDF总是报错

水沸腾UDF导入后显示：line 66:parse error,好不容易数到第66行，还是没发现错误....还请大家帮帮忙，看看到底哪里出错了，谢啦
#include "udf.h"          /*包括常规宏*/
#include "sg_mphase.h"    /*包括体积分数宏*/
#define T_SAT 373         /*定义饱和温度*/
#define LAT_HT 1.e3       /*定义水蒸汽潜热*/

DEFINE_SOURCE(liq_src, cell, pri_th, dS, eqn)  /*定义液相源项*/
{
  Thread *mix_th, *sec_th;   /*混合相、第二相 定义计算区域指针*/
  real m_dot_l;             /*蒸发冷凝速率  定义液相质量转移*/
  mix_th = THREAD_SUPER_THREAD(pri_th);   /*指向混合区的主相即液相的指针*/
  sec_th = THREAD_SUB_THREAD(mix_th, 1); /* 指向单相控制区的气相的指针*/
  if(C_T(cell, mix_th)>=T_SAT){
    m_dot_l = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)*
                fabs(C_T(cell, pri_th) - T_SAT)/T_SAT;  /*如果液相单元的温度高于蒸发温度，液相向气相的质量转移*/
   dS[eqn] = -0.1*C_R(cell, pri_th)*
      fabs(C_T(cell, pri_th) - T_SAT)/T_SAT; /*定义源项对质量转移偏导*/
                               }
  else {
    m_dot_l = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*
                fabs(T_SAT-C_T(cell,mix_th))/T_SAT;
    dS[eqn] = 0.; /*于是气相向液相转移，所以液相的质量源项对质量转移的偏导为零*/
       }
   
  return m_dot_l;
}

DEFINE_SOURCE(vap_src, cell, sec_th, dS, eqn) /*定义气相源项*/
{
  Thread * mix_th, *pri_th;    /*混合相、第一相*/
  real m_dot_v;
  mix_th = THREAD_SUPER_THREAD(sec_th);   /*指向混合区的第二相即气相的指针*/
  pri_th = THREAD_SUB_THREAD(mix_th, 0);  /*指向单相控制区的液相的指针，液相为主相*/
  if(C_T(cell, mix_th)>=T_SAT){
    m_dot_v = 0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)*
                fabs(C_T(cell, mix_th) - T_SAT)/T_SAT; /*如果混合区单元的温度高于蒸发温度，液相向气相的质量质量转移*/
    dS[eqn] = 0.;    /*由于是液相向气相转移，所以气相的质量源项对来自液相的质量转移的偏导为零*/
                               }
  else {
    m_dot_v = -0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*
                fabs(T_SAT-C_T(cell,mix_th))/T_SAT;
    dS[eqn] = -0.1*C_R(cell, sec_th)*
      fabs(C_T(cell, sec_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;
  pri_th = THREAD_SUB_THREAD(mix_th, 0);   /*指向混合区的液相的指针*/
  sec_th = THREAD_SUB_THREAD(mix_th, 1);   /*指向混合区的气相的指针*/
  if(C_T(cell, mix_th)>=T_SAT){
    m_dot = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)*
                fabs(C_T(cell, pri_th) - T_SAT)/T_SAT;  /*如果混合区的单元温度高于蒸发温度。质量转移由液相向气相转移，吸热，质量转移量前有负号*/
    dS[eqn] = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)/T_SAT;
                               }
  else {
    m_dot = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*
                fabs(T_SAT-C_T(cell,mix_th))/T_SAT;
    dS[eqn] = -0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)/T_SAT;}  /*相反，气相向液相转移则放热*/
  return LAT_HT*m_dot;    /*能量大小气化潜热与质量转移率相乘得能量源项*/
回复此楼

» 收录本帖的淘帖专辑推荐

多相流UDF 我的收藏

» 猜你喜欢
1楼 2011-05-16 16:01:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

47425804

木虫 (正式写手)

【答案】应助回帖

★ ★
xiegangmai(金币+2): 谢谢应助！ 2011-05-16 22:30:34
fengyezhi(金币+5): 2011-05-18 08:31:37
如果你的程序完全是这么写的话，那么就是最后面少了一个大括号~~~
赞一下(1人) 回复此楼
2楼2011-05-16 17:34:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

贾晓林

新虫 (初入文坛)

【答案】应助回帖

★
fengyezhi(金币+5): 2011-05-18 08:29:38
xiegangmai(金币+1): 谢谢应助！ 2011-05-19 21:29:42
是不是你的首相和第二项设置错误或者你的相名称没改过来？我没发现错误呵呵
赞一下(1人) 回复此楼
3楼2011-05-17 09:17:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fengyezhi

新虫 (小有名气)
引用回帖:
Originally posted by 贾晓林 at 2011-05-17 09:17:59:
是不是你的首相和第二项设置错误或者你的相名称没改过来？我没发现错误呵呵

我用多相流模型算的，liquid和vapor两相，但是每次算的时候liquid残差都很大，后来就算不下去了....我的模型里有传热
赞一下 回复此楼
4楼2011-05-18 08:31:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fengyezhi

新虫 (小有名气)
引用回帖:
Originally posted by 47425804 at 2011-05-16 17:34:17:
如果你的程序完全是这么写的话，那么就是最后面少了一个大括号~~~

这个语法错误我已经改过来了，谢谢你呵呵。但是还是算不下去，我用mixture模型，每次引入liq质量源项时就会报错，计算不下去。只引入vap源项和能量源项就可以算，不知道为什么。我的有传热
赞一下 回复此楼
5楼2011-05-18 08:33:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

47425804

木虫 (正式写手)
★ ★
xiegangmai(金币+2): 鼓励讨论交流！ 2011-05-19 21:29:24
引用回帖:
Originally posted by fengyezhi at 2011-05-18 08:33:22:
这个语法错误我已经改过来了，谢谢你呵呵。但是还是算不下去，我用mixture模型，每次引入liq质量源项时就会报错，计算不下去。只引入vap源项和能量源项就可以算，不知道为什么。我的有传热

那就是编译顺利通过啦~~
可否详述下报错信息。
另外，有相变的情况下残差就是容易很大，建议适当减小时间步长。
赞一下(1人) 回复此楼
6楼2011-05-19 09:54:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

倚天屠龙love

金虫 (小有名气)
建议楼主看下韩占忠09年出的那本关于Fluent的书的139页，里面那个UDF程序和你的一样，你对比一下吧~
另外，我想请教一下楼主这个UDF能用于VOF模型么？楼主现在的沸腾换热的计算进行的怎么样了，我现在也在搞这个，希望能和楼主交流一下！
赞一下 回复此楼
7楼2011-10-28 19:28:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yanxiujuan

新虫 (初入文坛)
楼主的问题解决了吗？我也是做气液相变换热的，刚开始没多久，麻烦问一下m_dot = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)* fabs(C_T(cell, pri_th) - T_SAT)/T_SAT是怎么导出来的呀？谢谢啦
赞一下 回复此楼
8楼2012-10-17 18:52:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小小鸟110

铜虫 (正式写手)
#define LAT_HT 1.e3       /*定义水蒸汽潜热*/ 请问这一句中的潜热，为什么是1000呢，怎么查啊
赞一下 回复此楼
追求心智的成熟
9楼2014-03-25 20:42:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gearboy78

禁虫 (小有名气)
本帖内容被屏蔽
10楼2014-03-26 08:25:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 fengyezhi 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？ +3 gouxfjh 2026-08-28 3/150 2026-08-29 01:01 by jurkat.1640
[基金申请] 国自然评审意见 +11 wangmingqi 2026-08-28 15/750 2026-08-29 00:17 by fangyl2005
[基金申请] 为什么资助数各大高校都创新高，自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国际合作可查了，中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 国合现在查不到了吗？ +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
信息提示
关闭
请填处理意见
关闭