应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » 求程序查错，关于多组分多相流指针问题，解决的话送上所有金币
161/2返回列表12下一页
查看: 1388  |  回复: 15
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

malu430

金虫 (小有名气)

[求助] 求程序查错，关于多组分多相流指针问题，解决的话送上所有金币

程序目的：计算licl溶液的表面张力系数
Licl溶液的表面张力系数是温度和浓度的函数

有高手能帮忙看下下面的程序，为什么液相指针找的不对呢？
C_UDMI(cell, mixture_thread, 0)的值按理应该是在液相等于sigma_water*(1+a1*MF+a2*MF*T/647.3+a3*MF*pow(T/647.3,2)+a4*pow(MF,2)+a5*pow(MF,3))。在气相应该等于初始值0，因为我对udm赋值只在液相中循环的。

但是情况却不是这样：C_UDMI(cell, mixture_thread, 0)的值在液相没问题，但
在气相C_UDMI(cell, mixture_thread, 0)= sigma_water
感觉这个程序的结果是并没有在液相循环，而是在整个区域都循环了，由于调用的浓度MF=0，所以在气相算出的结果就变成了sigma_water。
这个是什么原因呢？

/*计算表面张力系数，公式来源为参考文献的第页*/


DEFINE_PROPERTY(sigma_licl,cell,mixture_thread)
{
        return C_UDMI(cell, mixture_thread, 0);
}

DEFINE_ADJUST(my_adjust,domain)
{
        int zone_ID=13;//fluid的zone ID是，从界面获得
        int phase_ID=1;//液相的ID是
        int comp_ID=0;//组分licl的ID是
        cell_t cell;
        Thread *mixture_thread = Lookup_Thread(Get_Domain(1), zone_ID);//因为adjust不传输mixture_thread这个指针变量,所以先通过domain找到mixture_thread
        Thread *thread_l = THREAD_SUB_THREAD(mixture_thread, phase_ID);//通过mixture_thread找到液相指针

        real sigma_water;//水的表面张力
        real T;//温度，K
        real MF;//licl的浓度
   
        real a1=2.757115;
        real a2=-12.011299;
        real a3=14.751818;

        real a4=2.443204;
        real a5=-3.147739;

        real a6=235.8;
        real a7=0.625;
        real a8=1.256;

     begin_c_loop(cell, thread_l) //Licl溶液的表面张力系数是温度和浓度的函数，所以只在液相循环

        {
          T=C_T(cell,thread_l);//得到温度
          sigma_water=a6*(1-a7*(1-T/647.3))*pow((1-T/647.3),a8)*0.001;//得到水的表面张力系数
          MF= C_YI(cell, thread_l, comp_ID);//得到licl溶液的浓度
          C_UDMI(cell, mixture_thread, 0) = sigma_water*(1+a1*MF+a2*MF*T/647.3+a3*MF*pow(T/647.3,2)+a4*pow(MF,2)+a5*pow(MF,3));//得到licl溶液的表面张力系数
        }
    end_c_loop(cell, thread_l)
}
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-10-26 18:03:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小熊猫进行时
★
xiegangmai: 金币-1, 专业版块交流贴，请勿灌水 2012-10-27 20:13:59
2楼2012-10-26 18:43:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kylafree

至尊木虫 (知名作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
malu430: 金币+10 2012-10-26 22:21:49
xiegangmai: 金币+2, 谢谢应助 2012-10-27 20:14:05
1.对你的情况把DEFINE_ADJUST换为DEFINE_EXECUTE_AT_END可能好些。
2.Thread *mixture_thread = Lookup_Thread(Get_Domain(1), zone_ID);//因为adjust不传输mixture_thread这个指针变量,所以先通过domain找到mixture_thread---///////我的理解中DEFINE_ADJUST(name,d)中的参数d就是mixture domain的指针。
3.int phase_ID=1;//液相的ID是---///////这里的ID=1???? 我猜测可能=3！
4.要找液相不妨试下如下方法
Domain *d = Get_Domain(3);   /* d为 secondary phase domain(液相)指针（上面第3条中的猜测） ,3来自phase panel. 1-mixture phase; 2-primary phase; 3-secondary phase*/
Thread *t;
  thread_loop_c(t,d)
  {
    begin_c_loop(c, t)     
    {     
      if(0.5<=C_VOF(c,t)&&C_VOF(c,t)<=1.0)    /*找到液相CELL，，，交界面处一个CELL中既有气相又有液相如何处理？！*/
      {
赞一下 回复此楼
3楼2012-10-26 21:40:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

malu430

金虫 (小有名气)
引用回帖:
3楼: Originally posted by kylafree at 2012-10-26 21:40:05
1.对你的情况把DEFINE_ADJUST换为DEFINE_EXECUTE_AT_END可能好些。
2.Thread *mixture_thread = Lookup_Thread(Get_Domain(1), zone_ID);//因为adjust不传输mixture_thread这个指针变量,所以先通过domain找到mixtu ...

谢谢kylafree的热心帮助，不过有几点还是不提清楚的
1. 不太明白为什么要在AT_END执行，能再具体解释下吗？
2. d是domain，不是指针吧
3.我刚刚试过这样
  Domain *domain_l=Get_Domain(3);
    thread_loop_c(thread_l,domain_l)
{
     begin_c_loop(cell, thread_l)
        {
         。。。。。。       
                 }
    end_c_loop(cell, thread_l)
}
也就是你提到的，但是还是不行，和之前的结果完全一样
4. 我也试过 if(0.5<=C_VOF(c,t)&&C_VOF(c,t)<=1.0)，但是计算结果很诡异。。。。但有一点要说的是，
表面张力系数不需要找界面再定义的，因为表面张力源项里是有一个体积分数梯度的，不是界面的话，这个梯度就为0。这样就实现了在界面加表面张力源项。
赞一下 回复此楼
4楼2012-10-26 22:29:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kylafree

至尊木虫 (知名作家)
引用回帖:
3楼: Originally posted by kylafree at 2012-10-26 21:40:05
1.对你的情况把DEFINE_ADJUST换为DEFINE_EXECUTE_AT_END可能好些。
2.Thread *mixture_thread = Lookup_Thread(Get_Domain(1), zone_ID);//因为adjust不传输mixture_thread这个指针变量,所以先通过domain找到mixtu ...

我错了
回复此楼
5楼2012-10-27 11:16:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kylafree

至尊木虫 (知名作家)
★ ★
xiegangmai: 金币+2, 辛苦了！ 2012-10-27 20:14:22
试试下面这种方式看看。
1）
Thread *mixture_thread = Lookup_Thread(Get_Domain(1), zone_ID);  这行中把Get_Domain(1)换为domain

2）
DEFINE_ADJUST(my_adjust,domain)
{
         int zone_ID=13;//fluid的zone ID是，从界面获得
         int comp_ID=0;//组分licl的ID是
         cell_t cell;
         Domain *d_l = Get_Domain(3);   
         Thread *thread_l = Lookup_Thread(d_l, zone_ID);

        real sigma_water;//水的表面张力
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

malu430
6楼2012-10-27 16:33:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

malu430

金虫 (小有名气)
送鲜花一朵
引用回帖:
6楼: Originally posted by kylafree at 2012-10-27 16:33:44
试试下面这种方式看看。
1）
Thread *mixture_thread = Lookup_Thread(Get_Domain(1), zone_ID);  这行中把Get_Domain(1)换为domain

2）
DEFINE_ADJUST(my_adjust,domain)
{
         int zone_ID=13;//fl ...

恩 谢谢 试试看 嘿嘿
赞一下 回复此楼
7楼2012-10-28 14:10:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

malu430

金虫 (小有名气)
引用回帖:
6楼: Originally posted by kylafree at 2012-10-27 16:33:44
试试下面这种方式看看。
1）
Thread *mixture_thread = Lookup_Thread(Get_Domain(1), zone_ID);  这行中把Get_Domain(1)换为domain

2）
DEFINE_ADJUST(my_adjust,domain)
{
         int zone_ID=13;//fl ...

结果还是一样啊
回复此楼
8楼2012-10-28 14:38:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fangdian

木虫 (正式写手)
可以的话将case发上来，这么感觉那个ID不是很清楚。。
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

malu430 malu430
9楼2012-10-28 15:46:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

malu430

金虫 (小有名气)
送鲜花一朵
引用回帖:
9楼: Originally posted by fangdian at 2012-10-28 15:46:02
可以的话将case发上来，这么感觉那个ID不是很清楚。。

发上来了
嘿嘿 谢谢
赞一下 回复此楼
10楼2012-10-28 17:49:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 malu430 的主题更新
161/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭