24小时热门版块排行榜    

查看: 3181  |  回复: 8
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

linyi7445962

新虫 (初入文坛)

[求助] 有关水蒸汽冷凝的UDF 已有1人参与

一个球体模型,往里通入混合工质 其中包含水蒸汽 采用了多组分传输来确定混合工质的组成 需要考虑混合工质中水蒸气的冷凝情况 在网上查了下冷凝的UDF 都是多相流的 所以加入多相流模型 首相为混合工质 次相为水(液体) 在边界条件的fluid设置中的mixture选项中框选source terms 在engery中导入能量项 在首相中也框选source terms 并且在水蒸气项中导入vap_src 同理在次相中导入liq_src 设置好进行迭代 但是出现几种错误 1 能迭代但没产生质量转移 2不能进行迭代 提示错误 好像发散之类的 因为刚开始学习UDF 对于这方面不太理解 对于THREAD_SUPER_THREAD(pri_th)以及THREAD_SUB_THREAD(mix_th,1) 括号中的参数如何确定 以及指针到底指向哪里不是很清楚 在FLUENT HELP中看了下 也不是太明白 希望能帮忙解答下 或者有这方面经验的能交流下  以下是采用的UDF
由于刚进论坛 金币数太少 悬赏金少请见谅
#include "udf.h"
#include "sg_mphase.h"
#define T_SAT 373
#define LAT_HT 1.e3

DEFINE_SOURCE(vap_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,sec_th)*C_R(cell,sec_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT;

   dS[eqn]=0;
                               }
  else {
    m_dot_l=-0.1*C_VOF(cell,pri_th)*C_R(cell,pri_th)*fabs(T_SAT-C_T(cell,mix_th))/T_SAT;

    dS[eqn]=-0.1*C_R(cell,pri_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT;

       }
   
  return m_dot_l;

}

DEFINE_SOURCE(liq_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,sec_th)*C_R(cell,sec_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT;
    dS[eqn]=-0.1*C_R(cell,sec_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT;
                               }
  else {
    m_dot_v=0.1*C_VOF(cell,pri_th)*C_R(cell,pri_th)*fabs(T_SAT-C_T(cell,mix_th))/T_SAT;

    dS[eqn]=0;

       }
   
  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,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_VOF(cell,pri_th)*C_R(cell,pri_th)*fabs(T_SAT-C_T(cell,mix_th))/T_SAT;

    dS[eqn]=0.1*C_VOF(cell,pri_th)*C_R(cell,pri_th)/T_SAT;}

  return LAT_HT*m_dot;
}
回复此楼

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

Fluent中UDF

» 猜你喜欢

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

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

837729744

铜虫 (小有名气)

楼主问题解决了吗?这个水蒸气冷凝的udf,用到组分模型下,可以解决质量转移问题?应该还得改吧。
4楼2013-08-31 11:52:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

linyi7445962

新虫 (初入文坛)

对于质量转移问题解决 但是能量源项的问题在于导入模型之后进行迭代比不加UDF时 流体温度场下降速度很快 就像能量源项产生的是负能量 使温度下降很快,希望能解释enrg_src这块的错误在哪
2楼2012-05-31 14:59:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

linyi7445962

新虫 (初入文坛)

我的首相的ID是2 次相的ID是3 是否有关 有关应如何修改
3楼2012-05-31 16:55:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meiguimimi

银虫 (小有名气)

【答案】应助回帖

else {
    m_dot=0.1*C_VOF(cell,pri_th)*C_R(cell,pri_th)*fabs(T_SAT-C_T(cell,mix_th))/T_SAT;

    dS[eqn]=0.1*C_VOF(cell,pri_th)*C_R(cell,pri_th)/T_SAT;}

  return LAT_HT*m_dot;

在对质量的偏导数求导时,dS[eqn]项有一个负号
5楼2015-01-19 13:35:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 301求调剂 +6 yy要上岸呀 2026-03-17 6/300 2026-03-17 23:58 by 星空星月
[考研] 299求调剂 +4 △小透明* 2026-03-17 4/200 2026-03-17 20:09 by peike
[考研] 277调剂 +5 自由煎饼果子 2026-03-16 6/300 2026-03-17 19:26 by 李leezz
[考研] 328求调剂,英语六级551,有科研经历 +3 生物工程调剂 2026-03-16 8/400 2026-03-17 19:03 by Wangjingyue
[考研] 326求调剂 +5 上岸的小葡 2026-03-15 6/300 2026-03-17 17:26 by ruiyingmiao
[考研] 求调剂,总分315,考的生物医药,一志愿湖南师范大学。调剂到任何专业都可以 +4 小丁想进步 2026-03-11 5/250 2026-03-17 16:05 by 外星文明
[考研] 一志愿苏州大学材料工程(085601)专硕有科研经历三项国奖两个实用型专利一项省级立项 +6 大火山小火山 2026-03-16 8/400 2026-03-17 15:05 by 无懈可击111
[考研] 302求调剂 +4 小贾同学123 2026-03-15 8/400 2026-03-17 10:33 by 小贾同学123
[考研] 321求调剂 +5 大米饭! 2026-03-15 5/250 2026-03-16 16:33 by houyaoxu
[考研] 070303 总分349求调剂 +3 LJY9966 2026-03-15 5/250 2026-03-16 14:24 by xwxstudy
[考研] 求老师收留调剂 +4 jiang姜66 2026-03-14 5/250 2026-03-15 20:11 by Winj1e
[考博] 东华理工大学化材专业26届硕士博士申请 +6 zlingli 2026-03-13 6/300 2026-03-15 20:00 by ryzcf
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
[考研] 中科大材料与化工319求调剂 +3 孟鑫材料 2026-03-14 3/150 2026-03-14 20:10 by ms629
[考研] [0860]321分求调剂,ab区皆可 +4 宝贵热 2026-03-13 4/200 2026-03-13 22:01 by 星空星月
[考研] 一志愿西南交大,材料专硕317求调剂 +5 lx8568 2026-03-11 5/250 2026-03-13 21:43 by peike
[考研] 329求调剂 +3 miaodesi 2026-03-12 4/200 2026-03-13 20:53 by 18595523086
[考研] 材料工程调剂 +4 咪咪空空 2026-03-11 4/200 2026-03-13 19:57 by JourneyLucky
[考研] 工科调剂 +4 Jiang191123! 2026-03-11 4/200 2026-03-13 15:15 by Miko19
[考研] 求调剂 +3 程雨杭 2026-03-12 3/150 2026-03-13 15:06 by JourneyLucky
信息提示
请填处理意见