24小时热门版块排行榜    

查看: 1684  |  回复: 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个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 柔引引进博士硕士。 +15 大发财树 2024-12-18 15/750 2024-12-22 09:20 by ChemMisty
[考博] 申博的科研计划书怎么写? +6 爱喝风的龙卷茶 2024-12-19 10/500 2024-12-22 02:03 by l4q8i@6n
[职场人生] 国企、民企行业不同,如何决择? +3 haihe623 2024-12-19 5/250 2024-12-21 21:42 by haihe623
[论文投稿] 挑数据是造假吗 +5 电话的建设 2024-12-21 5/250 2024-12-21 17:16 by 科研皮皮猪
[教师之家] 今年的科研绩效只奖励Q1了 +12 holypower 2024-12-20 17/850 2024-12-21 16:13 by willbuilder
[教师之家] 学校排名靠前发展好,最大的受益者是校领导还是普通老师? +6 akslis2024 2024-12-17 7/350 2024-12-21 11:22 by 凌晨一点393
[考博] 读博的欲望达到了巅峰 +14 小丸九 2024-12-17 14/700 2024-12-21 07:27 by Skywit520
[有机交流] 工艺需求,付费 +4 想当然灬 2024-12-17 6/300 2024-12-21 06:23 by a893069456
[考博] 有偿聘请兼职硕士 +10 大发财树 2024-12-19 10/500 2024-12-20 17:30 by mopsite
[论文投稿] 论文校稿 50+5 whale_full 2024-12-17 10/500 2024-12-20 16:39 by 北京莱茵润色
[考博] 东北石油大学博士招生两名 +4 13305437324 2024-12-15 4/200 2024-12-20 12:30 by 一路禅修
[考博] 申博找导师 +6 sepra 2024-12-15 6/300 2024-12-20 04:51 by 胖胖的大海
[教师之家] 某老师在组会上总是针对我带的学生怎么办? +15 lice_1987 2024-12-17 23/1150 2024-12-20 04:50 by 胖胖的大海
[论文投稿] Molecular Crystals and Liquid Crystals投稿8个月没有信息 +4 书剑如风 2024-12-17 12/600 2024-12-20 04:40 by 胖胖的大海
[论文投稿] 求推荐5分二区以上肿瘤接受生信结合实验的期刊。。。 +3 lice_1987 2024-12-17 6/300 2024-12-19 17:01 by Andy_124
[论文投稿] 毕业论文疑虑 20+6 风吹荷叶煞 2024-12-17 17/850 2024-12-18 17:35 by holypower
[教师之家] 我感觉学生水平很差,随随便便工资上万 +19 guilin2018 2024-12-15 23/1150 2024-12-18 08:57 by zxf984
[有机交流] 细丝状的晶体可以测单晶吗? +5 好气好气 2024-12-16 5/250 2024-12-17 22:53 by 请让我滚去学习
[公派出国] 意大利帕多瓦大学招收CSC博士生--光电催化材料在新能源及可持续绿色化学方向 +3 huangyuxing2 2024-12-16 3/150 2024-12-17 18:25 by zhizhan127
[硕博家园] 大连海事大学双一流学科招收2025年博士研究生 +3 steresewan 2024-12-15 3/150 2024-12-16 00:17 by luweidong401
信息提示
请填处理意见