24小时热门版块排行榜    

查看: 585  |  回复: 1

guolongma

新虫 (初入文坛)


[交流] 如何将一个实时计算得到的值有效作用于一个源项中

最近在多相流的模拟中遇到了一个难题,还请各位高手指点啊!

想达到的目的是:首先初始化一个h值(数值很小),再通过DEFINE_ADJUST宏实时计算得到一个h值,然后将其带入到最后的源项中,将其作为源项的一个判定条件和变量。

遇到的问题:得到的h值作用不到
1、用DEFINE_INIT初始化h值
2、如何将DEFINE_INIT初始化的值以及DEFINE_ADJUST宏实时得到的值带入到源项中
3、下面是UDF,h和Zmin定义为了全局变量
(1)初始化这样直接赋值行吗?
DEFINE_INIT(my_int,d)
{
h=1.0e-5
}
(2)这个是否有问题?编译可以通过
DEFINE_ADJUST(adjust_fraction, d)
{
    Thread *t;
    cell_t c;
    Thread **pt;  
    real x,y,z,e[ND_ND];
    d = Get_Domain(1);
    mp_thread_loop_c (t,d,pt)   
        {
                if (FLUID_THREAD_P(t))
                begin_c_loop (c,t)
                {
                    C_CENTROID(e,c,t);
                z=e[2];
             if(C_VOF(c,pt[0])>0.3 && C_VOF(c,pt[0])<0.9) //根据VOF数值判断界面.
                {
              if(Zmin              Zmin=z;
                }
                h=fabs(Zmin);
                }
                end_c_loop (c,t)
        }
}
4、想实现目的,还有没有其他的途径呢?

希望高手多多指点,散尽家财,只为求得一良方啊!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

得得315

禁虫 (正式写手)


guolongma(金币+1): 谢谢参与
本帖内容被屏蔽

2楼2016-06-07 16:57:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 guolongma 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 8/400 2026-08-28 08:26 by lfy8008
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 XLGfUIdOmEAO 2026-08-27 3/150 2026-08-28 00:49 by 3RAFyYXYSwGS
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +7 苏东坡二世 2026-08-23 7/350 2026-08-27 18:14 by 瞬息宇宙
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 基金不中,共勉 +10 eulota 2026-08-26 10/500 2026-08-27 17:23 by lqllinqiaoli
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
信息提示
请填处理意见