24小时热门版块排行榜    

Znn3bq.jpeg
查看: 887  |  回复: 1

w420049436

新虫 (小有名气)

[求助] 做水冷凝结果不对,请教大神UDF哪里出错了? 已有1人参与

我现在正在学习水冷凝的UDF,然后上网找了例子模拟,发现结果出错,但是不知道错在哪里。恳求各位大神指正
模拟的参数如下:水蒸汽以400度进入管道,管道壁面温度为283度,水蒸汽的速度为0.1m/s。管径为0.05m,长度为3m。
#include "udf.h"
#include "sg_mphase.h"
#define T_SAT 373
#define LAT_HT 2.257e6



DEFINE_SOURCE(vap_src, cell, pri_th, dS, eqn)
{

  Thread *mix_th, *sec_th;
  real m_dot_v;

  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_v = 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_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.;

       }
   
  return m_dot_v;

}

DEFINE_SOURCE(lip_src, cell, sec_th, dS, eqn)

{

  Thread * mix_th, *pri_th;
  real m_dot_l;

  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_l = 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_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.1*C_R(cell, sec_th)*
      fabs(C_T(cell, sec_th) - T_SAT)/T_SAT;

       }
   
  return m_dot_l;

}

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哪里出错了?
这个是计算66s后的气相体积分数,根本看不到有水滴行程


做水冷凝结果不对,请教大神UDF哪里出错了?-1
这个是计算66s后的温度
回复此楼

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

博士科研新路历程

» 猜你喜欢

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

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

xiaobingwill

新虫 (初入文坛)

【答案】应助回帖

水蒸气速度很大,你把速度调小就应该能看到水珠了
一个人的战斗也要坚强
2楼2014-10-14 01:30:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 w420049436 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿西北工业大学289 085602 +24 yang婷 2026-04-10 24/1200 2026-04-10 10:25 by jiajinhpu
[考研] 0835学硕299求调剂 08大类可接受 +9 useryy 2026-04-03 9/450 2026-04-10 10:09 by 314126402
[考研] 083200 305分 求二轮调剂 不接受跨专业 +8 Claireyyyy 2026-04-09 8/400 2026-04-10 10:05 by luosha500
[考研] 289 分105500药学专硕求调剂(找B区学校) +5 白云123456789 2026-04-09 7/350 2026-04-09 21:03 by 白云123456789
[考研] 一志愿华工085600 331分 +6 天下ww 2026-04-09 6/300 2026-04-09 18:59 by l_paradox
[考研] 材料专硕初试分332一志愿西北工业大学, +12 故人?? 2026-04-09 12/600 2026-04-09 18:34 by Ccclqqq
[考研] 314求调剂 +16 xhhdjdjsjks 2026-04-09 17/850 2026-04-09 18:26 by 喝醉酒N
[考研] 296求调剂 +3 汪!?! 2026-04-08 3/150 2026-04-08 22:00 by zhouyuwinner
[考研] 求调剂,现在还能填的 +3 上岸小莹加油 2026-04-08 3/150 2026-04-08 14:30 by zhq0425
[考研] 265求调剂 +19 小木虫085600 2026-04-06 21/1050 2026-04-08 10:38 by 逆水乘风
[考研] 求调剂 +15 熊二想上岸 2026-04-06 15/750 2026-04-08 04:53 by 无际的草原
[考研] 277工科求调剂 +10 1915668 2026-04-04 11/550 2026-04-06 23:53 by 心媛意北
[考研] 考研调剂 +7 15615482637 2026-04-04 7/350 2026-04-06 22:56 by chenzhimin
[考研] 085500机械专硕初试288求调剂 +3 GZJguo666- 2026-04-05 3/150 2026-04-05 18:06 by jkddd
[考研] 调剂 +8 熊二想上岸 2026-04-04 8/400 2026-04-05 05:27 by houyaoxu
[考研] 11408,335分,本科211,求调剂,可转专业 +5 鳄梨大鳄鱼 2026-04-03 5/250 2026-04-04 22:49 by chongya
[考研] 323求调剂 +8 李佳乐1 2026-04-04 8/400 2026-04-04 22:26 by hemengdong
[考研] 考研调剂 +5 四川王涛 2026-04-04 5/250 2026-04-04 22:18 by 啵啵啵0119
[考研] 278求调剂 +14 范婷娜 2026-04-04 15/750 2026-04-04 22:15 by lqwchd
[考研] 305求调剂 +3 77Qi 2026-04-03 3/150 2026-04-03 23:01 by qzxyhcsy
信息提示
请填处理意见