24小时热门版块排行榜    

查看: 2554  |  回复: 10

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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 9/450 2026-08-18 20:09 by 奶牛小黑
[基金申请] 朋友圈看到的 +4 wangzilk 2026-08-18 5/250 2026-08-18 19:25 by wangzilk
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +51 医学老男孩 2026-08-13 110/5500 2026-08-18 19:07 by icerianxia
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 时间戳变了,能看出什么问题? +8 基诺咪客 2026-08-17 11/550 2026-08-18 17:13 by 悠然熊猫
[基金申请] 重要消息,中午系统在维护 +9 yuleib84 2026-08-18 10/500 2026-08-18 16:52 by archvillain
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +4 Tide man 2026-08-14 5/250 2026-08-18 11:35 by Tide man
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] filecode=后面第一个是大写字母 +8 wangze12014 2026-08-14 10/500 2026-08-17 17:05 by xter9665
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
信息提示
请填处理意见