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

24小时热门版块排行榜    

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

jesseliuxl

金虫 (小有名气)

[求助] 介个UDF有啥问题

做的液液传质，用VOF+species transport模拟液液不相容体系界面之间的传质，把传质模型简单化，界面之间的传质量等于传质系数*浓度差，在两相作用里的mass transfer里加的UDF，迭代计算之后在display—counters—species里没有发现有组分传到另一相中，大家看看问题在哪？


#include "udf.h"
#include "sg_mphase.h"
#define kl 0.00002

DEFINE_MASS_TRANSFER

(oil_water_source,cell,thread,from_index,from_species_index,to_index,to_species_index)
{
  real m_l;
   
  Thread *oil = THREAD_SUB_THREAD(thread,from_index);
  Thread *water = THREAD_SUB_THREAD(thread,to_index);
  if (C_VOF(cell,oil)==1)
{
  m_l=0;

}
else if ((C_VOF(cell,oil)>0)&&(C_VOF(cell,oil)<1))

{
  m_l=kl*(C_YI(cell,oil,0)-C_YI(cell,water,0));

  }
return m_l;
}
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
生命不息折腾不止
1楼2011-06-09 08:44:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

NaOHWTO

金虫 (小有名气)

【答案】应助回帖

★ ★
jesseliuxl(金币+3): 2011-06-11 10:20:50
臭水沟(金币+2): 谢谢应助~~ 2011-06-11 12:50:41
else 那里是不是有点啰嗦了，C_VOF(c,t)本来的范围就是0到1之间。如果你是为了寻找两相间的界面的话，直接采用C_VOF(c,t)=0.5来判断。
if(C_VOF(c,t)>0.5)
m_l=0;
else
m_l=kl*(C_YI(cell,oil,0)-C_YI(cell,water,0));

VOF模型的volume fraction用几何重构方法吧。
另外
m_l=kl*(C_YI(cell,oil,0)-C_YI(cell,water,0));
应该不是用零吧，设0应该表示为混合相。
看你设置什么是主相了，如果water是主相，
m_l=kl*(C_YI(cell,oil,2)-C_YI(cell,water,1));
修改试试！
赞一下(3人) 回复此楼
6楼2011-06-11 09:04:57
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

47425804

木虫 (正式写手)

【答案】应助回帖

★
jesseliuxl(金币+3): 2011-06-10 09:57:12
臭水沟(金币+1): 谢谢应助~~ 2011-06-11 12:50:10
UDF写得虽然不是很严谨，但应该不会出现没有组分传递吧。
估计是模型设置上有些问题。
赞一下(1人) 回复此楼
2楼2011-06-09 14:56:01
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jesseliuxl

金虫 (小有名气)
引用回帖:
Originally posted by 47425804 at 2011-06-09 14:56:01:
UDF写得虽然不是很严谨，但应该不会出现没有组分传递吧。
估计是模型设置上有些问题。

模型上应该没有问题啊
赞一下 回复此楼
生命不息折腾不止
3楼2011-06-10 09:58:04
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tb701

金虫 (著名写手)

【答案】应助回帖

jesseliuxl(金币+3): 2011-06-11 08:23:22
m_l=kl*(C_YI(cell,oil,0)-C_YI(cell,water,0));
C-YI(c,t,i)其中i为species index .
赞一下 回复此楼
4楼2011-06-10 21:42:42
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tb701

金虫 (著名写手)

【答案】应助回帖

★
臭水沟(金币+1): 谢谢应助~~ 2011-06-10 22:10:49
jesseliuxl(金币+3): 2011-06-11 08:23:28
((C_VOF(cell,oil)>0，改成((C_VOF(cell,oil)>0.01试试
赞一下(1人) 回复此楼
5楼2011-06-10 21:44:34
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jesseliuxl

金虫 (小有名气)
我再算算试试，谢谢各位
赞一下 回复此楼
生命不息折腾不止
7楼2011-06-11 10:19:49
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lihesong

金虫 (小有名气)
楼主 你的问题解决没有？
一直以来vof模拟扩散就存在一些问题。例如，vof模型假设两相之间是不能穿插渗透，但是如果考虑扩散的话，是不是和这个前提矛盾。
第二，如果通过vof的值来表征扩散程度的话，就有一个问题——初始化的时候，界面两侧的网格被定义了一个体积分数初值，定义初值的目的是为了使体积分数在界面处连续。但是做数值计算时，初值的选取极大地影响了解的结果。不知楼主在这个问题上考虑没有！希望多交流交流。
赞一下 回复此楼
8楼2012-09-14 16:45:58
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luowenxu

木虫 (小有名气)
楼主 这个传质问题解决没有？udf有没有错！
赞一下 回复此楼
若想人前显贵，必先人后受罪
9楼2012-11-30 14:19:41
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jdzlplq

铁虫 (小有名气)
请问LZ，你是怎么定义oil属性的？是在material中自定义的吗？那么，自定义的属性的准确性如何，能贴出你油的属性吗？
最近在做制冷剂与润滑油的混合问题，润滑油是混合物，很多性质不是准确值，请问lz是怎么解决的？
赞一下 回复此楼
10楼2013-05-15 11:01:00
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jesseliuxl 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭