24小时热门版块排行榜    

查看: 586  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上合作单位盖章 +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 瞬息宇宙
[基金申请] 基金不中,共勉 +10 eulota 2026-08-26 10/500 2026-08-27 17:23 by lqllinqiaoli
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 系统查不到 +10 董八千 2026-08-26 10/500 2026-08-26 16:30 by Equinoxhua
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 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 宝贝虫子
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国合现在查不到了吗? +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
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
信息提示
请填处理意见