24小时热门版块排行榜    

查看: 1526  |  回复: 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个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 希望今年自己国自然面上项目和老婆青年项目能中! +9 恐龙爸爸 2024-06-14 9/450 2024-06-20 08:04 by keyaner23
[找工作] 药学硕士找不到工作,打算去做科研助理了 +9 pom戴墨镜 2024-06-14 22/1100 2024-06-20 05:47 by 投必得科研顾问
[考博] 有机化学迷茫学生 +4 佛系摸鱼5 2024-06-18 5/250 2024-06-20 00:05 by 佛系摸鱼5
[基金申请] 江南大学到瑞士招聘,称取消非升即走,改预聘+长聘 +21 babu2015 2024-06-18 22/1100 2024-06-19 23:03 by feng6531
[海外博后] 墨尔本大学博后offer要不要接 +3 kyxblmm 2024-06-18 3/150 2024-06-19 22:39 by blake1111
[基金申请] 2024国社科通讯评审 +21 qsd10086 2024-06-13 43/2150 2024-06-19 20:57 by 05111401022
[考博] 34岁读博士晚吗 +39 emitdne 2024-06-13 39/1950 2024-06-19 19:03 by Dqtjs
[论文投稿] 审稿 +5 香瓜木香 2024-06-19 6/300 2024-06-19 17:44 by xli1984
[基金申请] 基金得中 +4 woaini0218 2024-06-18 4/200 2024-06-19 17:27 by mengzl
[基金申请] Nature 11日发文,中国著名学者们称造假迫不得已 +8 babu2015 2024-06-14 8/400 2024-06-19 15:25 by 风今25
[教师之家] 试用期辞职 +9 ZHONGWU_U 2024-06-18 14/700 2024-06-19 12:36 by ZHONGWU_U
[论文投稿] ACS AMI 返回审稿意见,一个大修,两个据稿,编辑给的修改重投 +5 智商已更新 2024-06-19 5/250 2024-06-19 12:35 by nono2009
[论文投稿] 审稿人含糊拒稿,还需要回复吗?如何回复? 20+4 BruceChum 2024-06-15 22/1100 2024-06-19 08:00 by kanyechris
[找工作] 杭电、天津科技、青农和宁波工程学院如何选? +12 味道很好啊 2024-06-13 24/1200 2024-06-18 23:09 by 药材
[教师之家] 请问事业编制和年薪制冲突吗? +12 ZHONGWU_U 2024-06-14 12/600 2024-06-18 19:31 by fangyl2005
[有机交流] 跑板能跑开,过柱过不纯怎么办 +4 小胡在努力 2024-06-18 6/300 2024-06-18 15:20 by long8811
[高分子] 烧瓶内合成聚酯 +3 大帝国乐 2024-06-17 7/350 2024-06-18 11:10 by 大帝国乐
[硕博家园] 博士毕业高校和就业的相关问题 +7 SCITOPPP 2024-06-14 11/550 2024-06-18 07:51 by yinxing1995
[论文投稿] 投稿后发现其他作者的邮箱填错了该怎么办呀 10+4 在飞的猪 2024-06-13 6/300 2024-06-14 04:45 by 小虫子咔咔
[基金申请] 博士后面上项目状态还是专家评审吗 10+9 Thatcheremu 2024-06-13 55/2750 2024-06-13 21:23 by 乌合麒麟
信息提示
请填处理意见