24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

lwzhao1990

木虫 (小有名气)

应助: 5 (幼儿园)
金币: 2057.2
散金: 610
红花: 3
帖子: 146
在线: 31.1小时
虫号: 660539
注册: 2008-11-23
性别: GG
专业: 工程热物理相关交叉领域

[交流] 【求助】冷凝的UDF程序已有7人参与

哪位仁兄有冷凝的UDF程序？上传一下，谢谢了！

回复此楼

» 猜你喜欢

导师想让我从独立一作变成了共一第一已经有9人回复
博士读完未来一定会好吗已经有23人回复
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有11人回复
读博已经有4人回复
JMPT 期刊投稿流程已经有4人回复
心脉受损已经有5人回复
Springer期刊投稿求助已经有4人回复
小论文投稿已经有3人回复
申请2026年博士已经有6人回复

1楼2009-10-29 09:13:02

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

1791887345

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 1
帖子: 3
在线: 1.5小时
虫号: 2436687
注册: 2013-04-25
专业: 分离过程

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

你的问题解决了吗？我也是类似的问题呢

赞一下

回复此楼

16楼2015-10-13 21:38:50

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

查看全部 16 个回答

tacaro

至尊木虫 (知名作家)

小木虫大能猫骑士勋章

仿真EPI: 1
应助: 21 (小学生)
贵宾: 0.06
金币: 32972.5
散金: 1128
红花: 17
沙发: 2
帖子: 6205
在线: 1291.9小时
虫号: 663213
注册: 2008-11-27
专业: 武术理论与方法

★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
nono2009(金币+3,VIP+0):今天表现不错！ 10-29 12:01

给你一份网上流传比较多的计算冷凝蒸发的udf

#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;
}

希望能有点启发