24小时热门版块排行榜    

查看: 1530  |  回复: 10
【悬赏金币】回答本帖问题,作者jesseliuxl将赠送您 8 个金币

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;
}
回复此楼

» 猜你喜欢

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

生命不息折腾不止
已阅   回复此楼   关注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));
修改试试!
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写得虽然不是很严谨,但应该不会出现没有组分传递吧。
估计是模型设置上有些问题。
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试试
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 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 博后面上今天下午会公布吗?大家有无消息? +10 地理学1995 2024-06-24 14/700 2024-06-27 11:24 by jackiddd
[有机交流] 脱酯基合成步骤 5+4 小辉灰@anhui 2024-06-25 8/400 2024-06-27 11:17 by 火龙骑士
[有机交流] 大佬们,打的核磁氢谱与chemdraw预测的有偏差 +7 来了个 2024-06-26 7/350 2024-06-27 11:11 by wubiansiya
[硕博家园] 考博英语 +6 198新用户 2024-06-25 11/550 2024-06-27 10:13 by Bletilla
[基金申请] 国基在研影响申请结果吗 +13 WOWO159357 2024-06-26 21/1050 2024-06-27 09:44 by WOWO159357
[硕博家园] 回家两天,不想呆了 +7 368ghnf 2024-06-22 13/650 2024-06-26 19:07 by ZZZemmm
[有机交流] 过柱子,产品,杂质在是 pe:ea=100:1 也一起出来? +5 w256 2024-06-25 5/250 2024-06-26 09:31 by 小木木cc
[基金申请] 基金申请书名称有变化 +5 xuel2011 2024-06-25 7/350 2024-06-26 00:07 by 老虎当猫养
[基金申请] 这样的说辞是上会了吗 +12 学员d3zYCz 2024-06-24 15/750 2024-06-25 22:05 by 3115321
[基金申请] 博后面上到底什么时候出结果??? +7 爱学lsy 2024-06-24 7/350 2024-06-25 19:24 by 暴走的蒜泥
[基金申请] 能看出是否上会了吗 +10 articlefan 2024-06-23 15/750 2024-06-25 16:05 by 请慎重修改昵称
[基金申请] 2024安徽省哲社立项名单公示 +3 robin_work 2024-06-24 6/300 2024-06-25 14:37 by Pickfoot
[基金申请] 焦虑没有毛线用,默默前行是王道! +4 漠上藜梭 2024-06-24 8/400 2024-06-25 14:32 by 漠上藜梭
[教师之家] 复旦夏同学退学理由说明,哪儿可以下载? +7 苏东坡二世 2024-06-21 8/400 2024-06-25 12:42 by 药嘿1233
[教师之家] 有没有今年的影响因子? +6 jurkat.1640 2024-06-22 10/500 2024-06-24 22:48 by 潇湘之迷
[教师之家] 复旦夏同学提出高校成年人学生退学不应该让家长审核,大家认同吗? +10 苏东坡二世 2024-06-22 17/850 2024-06-24 16:52 by wanghuawei
[金属] 寻找钛合金热压缩代做 +4 liuyang358 2024-06-23 4/200 2024-06-24 13:50 by 搬砖狗不放弃
[基金申请] 说博后基金7月出的真打电话了吗? +12 antonysole 2024-06-24 14/700 2024-06-24 13:39 by sizhouyi
[基金申请] 国自然资助比率是不是要下降了?? +8 今晚推荐22 2024-06-21 12/600 2024-06-24 11:15 by Pickfoot
[催化] 求助原位红外 +3 jfdhj 2024-06-21 7/350 2024-06-24 10:52 by qvhm2609
信息提示
请填处理意见