24小时热门版块排行榜    

查看: 1682  |  回复: 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个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 想找一位年轻的985博导 +5 艾弗森迪尔 2024-12-21 6/300 2024-12-21 22:11 by 艾弗森迪尔
[硕博家园] 四线城市公办高校无编or二线城市民办高校 +13 joshhartford 2024-12-20 29/1450 2024-12-21 20:52 by joshhartford
[论文投稿] MDPI投稿连续拒稿是什么原因?文章质量不好吗?还是其他原因? 2+8 多听多看多学 2024-12-17 17/850 2024-12-21 18:51 by 6681981668
[教师之家] 40多岁的教师危机感很强啊 (金币+5) +21 farmer1999 2024-12-16 28/1400 2024-12-21 16:07 by willbuilder
[有机交流] 装柱干法上样 +3 刘松垚 2024-12-20 3/150 2024-12-21 14:37 by hanjianwei
[找工作] 柔引引进博士硕士。 +14 大发财树 2024-12-18 14/700 2024-12-21 12:41 by xiaofengyue1
[考博] 求老师哦! +4 景小凯 2024-12-20 4/200 2024-12-21 12:09 by ujschem
[教师之家] 学校排名靠前发展好,最大的受益者是校领导还是普通老师? +6 akslis2024 2024-12-17 7/350 2024-12-21 11:22 by 凌晨一点393
[硕博家园] “认怂了”还是“认清了”:青年博士纷纷“逃离”科研的背后…… +13 柏舟0000 2024-12-16 14/700 2024-12-21 08:36 by steven_198377
[论文投稿] 论文校稿 50+5 whale_full 2024-12-17 10/500 2024-12-20 16:39 by 北京莱茵润色
[硕博家园] 柔性引进硕博人才 +6 大发财树 2024-12-20 6/300 2024-12-20 15:32 by owlmac
[考博] 申博找导师 +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 胖胖的大海
[找工作] 水电解制氢方向求职 +5 kyle1999 2024-12-19 8/400 2024-12-20 04:48 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
[论文投稿] 论文投稿 +4 @星尘@ 2024-12-17 6/300 2024-12-19 09:39 by YXLSYJK
[教师之家] 我感觉学生水平很差,随随便便工资上万 +19 guilin2018 2024-12-15 23/1150 2024-12-18 08:57 by zxf984
[基金申请] 网传董晨院士涉嫌24篇论文造假,本人回应:非恶意造假 +6 babu2015 2024-12-17 6/300 2024-12-18 08:46 by llhljsy
[硕博家园] 大连海事大学双一流学科招收2025年博士研究生 +3 steresewan 2024-12-15 3/150 2024-12-16 00:17 by luweidong401
信息提示
请填处理意见