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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » 有关水蒸汽冷凝的UDF
91/1返回列表
查看: 3172  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-05-28 16:58:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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的回帖

837729744

铜虫 (小有名气)
楼主问题解决了吗？这个水蒸气冷凝的udf，用到组分模型下，可以解决质量转移问题？应该还得改吧。
赞一下 回复此楼
4楼2013-08-31 11:52:28
已阅   回复此楼   关注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的回帖

LRM2015

新虫 (初入文坛)
引用回帖:
5楼: Originally posted by meiguimimi at 2015-01-19 13:35:17
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=0.1*C_VOF(cell,pri_th)*C_R(cell,pri_th)/T_SAT;}

  return LAT_HT*m_dot;

在对质量的偏导数求 ...

前面的0.1是什么意思？
赞一下 回复此楼
6楼2015-12-31 20:29:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zjdycjl2012

新虫 (小有名气)
引用回帖:
6楼: Originally posted by LRM2015 at 2015-12-31 20:29:24
前面的0.1是什么意思？...

相变强度因子，是一个系数
赞一下 回复此楼
7楼2017-03-31 17:02:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zjdycjl2012

新虫 (小有名气)
用不用多组分模型会影响UDF吗？需要修改UDF吗？
赞一下 回复此楼
8楼2017-03-31 17:12:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

PeiyuChen

新虫 (初入文坛)
引用回帖:
7楼: Originally posted by zjdycjl2012 at 2017-03-31 17:02:53
相变强度因子，是一个系数...

楼主，问题解决了吗？我的模拟也是第二相液态水产生不了。

发自小木虫Android客户端
赞一下 回复此楼
9楼2017-04-15 14:24:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 linyi7445962 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 22408总分284求调剂 +3 InAspic 2026-03-13 3/150 2026-03-15 11:10 by zhq0425
[考研] 283求调剂 +6 小楼。 2026-03-12 9/450 2026-03-14 22:15 by cfxczx2015
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[考研] 一志愿哈工大材料324分求调剂 +5 闫旭东 2026-03-14 5/250 2026-03-14 14:53 by 木瓜膏
[考研] 337一志愿华南理工材料求调剂（有希望2吗？） +3 mysdl 2026-03-09 3/150 2026-03-14 02:53 by JourneyLucky
[考研] 301求调剂 +3 归零lbm 2026-03-09 3/150 2026-03-14 02:20 by JourneyLucky
[考研] 环境调剂 +6 晓看天暮看云 2026-03-09 6/300 2026-03-14 01:16 by JourneyLucky
[考研] 271求调剂 +10 生如夏花… 2026-03-11 10/500 2026-03-14 00:35 by 卖报员小雨
[考研] 318求调剂 +3 李新光 2026-03-10 3/150 2026-03-14 00:21 by JourneyLucky
[考研] 285化工学硕求调剂（081700） +6 柴郡猫_ 2026-03-12 6/300 2026-03-13 20:46 by hmn_wj
[考研] 【考研调剂求收留】 +3 Ceciilia 2026-03-11 3/150 2026-03-13 20:18 by JourneyLucky
[考研] 材料工程调剂 +4 咪咪空空 2026-03-11 4/200 2026-03-13 19:57 by JourneyLucky
[考研] 【0856】化学工程（085602）313 分，本科学科评估A类院校化学工程与工艺，诚求调剂 +7 小刘快快上岸 2026-03-11 7/350 2026-03-13 16:06 by ruiyingmiao
[考研] 295求调剂 +3 小匕仔汁 2026-03-12 3/150 2026-03-13 15:17 by vgtyfty
[论文投稿] 投稿问题 5+4 星光灿烂xt 2026-03-12 6/300 2026-03-13 14:17 by god_tian
[考研] 0857环境调剂 +5 熠熠_11 2026-03-10 5/250 2026-03-11 10:59 by wang_dand
[考研] 大连大学化学专业研究生调剂 +3 琪久. 2026-03-10 8/400 2026-03-11 10:02 by 琪久.
[基金申请] 提交后的基金本子，已让学校撤回了，可否换口子提交 +3 dut_pfx 2026-03-10 3/150 2026-03-11 08:38 by kudofaye
[考研] 调剂 +5 呵唔哦豁 2026-03-10 5/250 2026-03-10 22:00 by 28375m
[考研] 0703化学调剂 +3 三dd. 2026-03-10 3/150 2026-03-10 15:45 by peike
信息提示
关闭
请填处理意见
关闭