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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » fluent udf DEFINE_ADJUST求两相界面力能load但不能计算
北京石油化工学院2026年研究生招生接收调剂公告
51/1返回列表
查看: 1763  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

guoyhchina

金虫 (小有名气)

[求助] fluent udf DEFINE_ADJUST求两相界面力能load但不能计算 已有1人参与

我用vof法计算气液两相流，想先找出两相界面，并求界面电密度，用DEFINE_ADJUST和DEFINE_ON_DEMAND都试过，能顺利load进去，但一hook这个界面电密度文件fluent就提示致命错误。下面是我的部分程序，请高手指导一下，先谢了。

#include "stdio.h"
#include "udf.h"
#include "sg.h"
#include "sg_mphase.h"
#include "math.h"
#include "flow.h"
#include "mem.h"
#include "metric.h"
#include "id.h"

#define EPS_1 4.25   /* 相对介电常数*/
#define EPS_0 8.8541878e-12  /* electrical permittivity 真空介电常数c^2/Nm^2*/


/* Define which user-defined scalars to use. */

enum                 /*枚举类型名*/
{
        phi,    /* 电势*/  

        N_REQUIRED_UDSI    /*存取用户自定义标量输运方程的数目*/
};


enum                 /*枚举类型名*/
{
        Ex,
        Ey,
        Ez,
        vfx,
        vfy,
        vfz,
        rhos,  /* 界面电密度*/
        N_REQUIRED_UDM    /*存取用户自定义标量输运方程的数目*/
};

DEFINE_ADJUST(store_VOF_gradient, d)  /*存储vof体积分数的梯度*/
{
        Thread *t;
        Thread *ppt;
        Thread **pt;
        cell_t c;
        int phase_domain_index=1;
        Domain *pDomain = DOMAIN_SUB_DOMAIN(d,phase_domain_index);
        Alloc_Storage_Vars(pDomain,SV_VOF_RG,SV_VOF_G,SV_NULL);
        Scalar_Reconstruction(pDomain, SV_VOF,-1,SV_VOF_RG,NULL);
        Scalar_Derivatives(pDomain,SV_VOF,-1,SV_VOF_G,SV_VOF_RG, Vof_Deriv_Accumulate);

        mp_thread_loop_c(t,d,pt)
        {
                if (FLUID_THREAD_P(t))
                {
                        ppt = pt[phase_domain_index];

                        begin_c_loop (c,t)
                        {
                                C_UDMI(c,t,vfx) = C_VOF_G(c,ppt)[0];
                                C_UDMI(c,t,vfy) = C_VOF_G(c,ppt)[1];
                                C_UDMI(c,t,vfz) = C_VOF_G(c,ppt)[2];
                        }
                        end_c_loop (c,t)
                }
        }

        Free_Storage_Vars(pDomain,SV_VOF_RG,SV_VOF_G,SV_NULL);
}


DEFINE_ADJUST(Calc_E_xyz, domain)  /*存储电场强度*/
{

        Thread *t;
        cell_t c;
        face_t f;
        Thread *t0;
        cell_t c0;
        real dr0[ND_ND],dr1[ND_ND], xf[ND_ND], dy;


&nbsp; &nbsp; &nbsp; &nbsp; if (n_uds < N_REQUIRED_UDSI)

&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; Internal_Error("not enough user-defined scalars allocated";


&nbsp; &nbsp; &nbsp; &nbsp; if(n_udm < N_REQUIRED_UDM)

&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; Error("Not enough udm allocated\n";


&nbsp; &nbsp; &nbsp; &nbsp; /* Do nothing if gradient isn't allocated yet. */

&nbsp; &nbsp; &nbsp; &nbsp; if (!Data_Valid_P()) return;

&nbsp; &nbsp; &nbsp; &nbsp; thread_loop_c(t,domain)
&nbsp; &nbsp; &nbsp; &nbsp; {
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; if (FLUID_THREAD_P(t))
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; if (NULL != THREAD_STORAGE(t,SV_UDS_I(phi)) && NULL != T_STORAGE_R_NV (t,SV_UDSI_G(phi)))
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; {
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; begin_c_loop(c,t)&nbsp; &nbsp;&nbsp; &nbsp; /*体循环*/
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; {
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; C_UDMI(c,t,Ex) = -1.*C_UDSI_G(c,t,phi)[0];&nbsp; &nbsp; /*Ex*/
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; C_UDMI(c,t,Ey) = -1.*C_UDSI_G(c,t,phi)[1];
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; C_UDMI(c,t,Ez) = -1.*C_UDSI_G(c,t,phi)[2];
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; }
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; end_c_loop(c,t)
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; /* Message("x electric strength: %g\n", C_UDMI(c,t,Ex)); */
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; }
&nbsp; &nbsp; &nbsp; &nbsp; }

}


DEFINE_ADJUST(calm_rhos, domain)&nbsp; &nbsp;/* 界面电密度，出问题的地方*/
{&nbsp; &nbsp; &nbsp; &nbsp;
&nbsp; &nbsp; &nbsp; &nbsp; Thread *t;
&nbsp; &nbsp; &nbsp; &nbsp; Thread *pri_th,*sec_th;
&nbsp; &nbsp; &nbsp; &nbsp; cell_t c;
&nbsp; &nbsp; &nbsp; &nbsp; real xc[ND_ND], EEp[3], EEs[3], nn[3], EE[3], rr;
&nbsp; &nbsp; &nbsp; &nbsp; /*&nbsp;&nbsp;Domain *domain = Get_Domain(1); Get domain pointer */
&nbsp; &nbsp; &nbsp; &nbsp; int ID = 1;
&nbsp; &nbsp; &nbsp; &nbsp; t =Lookup_Thread(domain, ID); /* mixture-level thread pointer */

&nbsp; &nbsp; &nbsp; &nbsp; pri_th = THREAD_SUB_THREAD(t,0);
&nbsp; &nbsp; &nbsp; &nbsp; sec_th = THREAD_SUB_THREAD(t,1);
&nbsp; &nbsp; &nbsp; &nbsp;
&nbsp; &nbsp; &nbsp; &nbsp; thread_loop_c(t,domain)
&nbsp; &nbsp; &nbsp; &nbsp; {
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; begin_c_loop(c,t)&nbsp; &nbsp;&nbsp; &nbsp; /*体循环*/
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; {
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; C_CENTROID(xc,c,t);
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; if (C_VOF(c,sec_th)>0. && C_VOF(c,sec_th)<1.)
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; {
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; EEs[0]=C_UDMI(c,sec_th,Ex);/* 次相电场强度 */
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; EEs[1]=C_UDMI(c,sec_th,Ey);
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; EEs[2]=C_UDMI(c,sec_th,Ez);
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; EEp[0]=C_UDMI(c,pri_th,Ex);/* 基相电场强度 */
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; EEp[1]=C_UDMI(c,pri_th,Ey);&nbsp; &nbsp; &nbsp; &nbsp;
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; EEp[2]=C_UDMI(c,pri_th,Ez);
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; nn[0]=C_UDMI(c,sec_th,vfx);/* vof梯度、法向矢量 */
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; nn[1]=C_UDMI(c,sec_th,vfy);
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; nn[2]=C_UDMI(c,sec_th,vfz);
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; NV_V(EEs,=,EEp);
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; rr=EPS_0*EPS_1;
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; NV_VS_VS(EE, =,EEp,*,EPS_0, -, EEs,*,rr);
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; C_UDMI(c,t,rhos) = NV_DOT(nn,EE);
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; }
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; else
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; {
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; C_UDMI(c,t,rhos) =0;
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; }
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; }
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; end_c_loop(c,t)
&nbsp; &nbsp; &nbsp; &nbsp; }
}
回复此楼

» 猜你喜欢
1楼 2017-02-26 11:36:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guoyhchina

金虫 (小有名气)
怎么全是乱码了，重新贴上程序。

#include "stdio.h"
#include "udf.h"
#include "sg.h"
#include "sg_mphase.h"
#include "math.h"
#include "flow.h"
#include "mem.h"
#include "metric.h"
#include "id.h"

#define EPS_1 4.25&nbsp;&nbsp; /* 相对介电常数*/
#define EPS_0 8.8541878e-12&nbsp; /* electrical permittivity 真空介电常数c^2/Nm^2*/

/* Define which user-defined scalars to use. */
enum&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; /*枚举类型名*/
{
&nbsp;phi,&nbsp;&nbsp;&nbsp; /* 电势*/&nbsp;
&nbsp;elrho,&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; /* 电密度*/
&nbsp;N_REQUIRED_UDSI&nbsp;&nbsp;&nbsp; /*存取用户自定义标量输运方程的数目*/
};

enum&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; /*枚举类型名*/
{
&nbsp;Ex,
&nbsp;Ey,
&nbsp;Ez,
&nbsp;vfx,
&nbsp;vfy,
&nbsp;vfz,
&nbsp;rhos,&nbsp; /* 界面电密度*/
&nbsp;N_REQUIRED_UDM&nbsp;&nbsp;&nbsp; /*存取用户自定义标量输运方程的数目*/
};
DEFINE_ADJUST(store_VOF_gradient, d)&nbsp; /*存储vof体积分数的梯度*/
{
&nbsp;Thread *t;
&nbsp;Thread *ppt;
&nbsp;Thread **pt;
&nbsp;cell_t c;
&nbsp;int phase_domain_index=1;
&nbsp;Domain *pDomain = DOMAIN_SUB_DOMAIN(d,phase_domain_index);
&nbsp;Alloc_Storage_Vars(pDomain,SV_VOF_RG,SV_VOF_G,SV_NULL);
&nbsp;Scalar_Reconstruction(pDomain, SV_VOF,-1,SV_VOF_RG,NULL);
&nbsp;Scalar_Derivatives(pDomain,SV_VOF,-1,SV_VOF_G,SV_VOF_RG, Vof_Deriv_Accumulate);
&nbsp;mp_thread_loop_c(t,d,pt)
&nbsp;{
&nbsp;&nbsp;if (FLUID_THREAD_P(t))
&nbsp;&nbsp;{
&nbsp;&nbsp;&nbsp;ppt = pt[phase_domain_index];
&nbsp;&nbsp;&nbsp;begin_c_loop (c,t)
&nbsp;&nbsp;&nbsp;{
&nbsp;&nbsp;&nbsp;&nbsp;C_UDMI(c,t,vfx) = C_VOF_G(c,ppt)[0];
&nbsp;&nbsp;&nbsp;&nbsp;C_UDMI(c,t,vfy) = C_VOF_G(c,ppt)[1];
&nbsp;&nbsp;&nbsp;&nbsp;C_UDMI(c,t,vfz) = C_VOF_G(c,ppt)[2];
&nbsp;&nbsp;&nbsp;}
&nbsp;&nbsp;&nbsp;end_c_loop (c,t)
&nbsp;&nbsp;}
&nbsp;}
&nbsp;Free_Storage_Vars(pDomain,SV_VOF_RG,SV_VOF_G,SV_NULL);
}

DEFINE_ADJUST(Calc_E_xyz, domain)&nbsp; /*存储电场强度*/
{
&nbsp;Thread *t;
&nbsp;cell_t c;
&nbsp;face_t f;
&nbsp;Thread *t0;
&nbsp;cell_t c0;
&nbsp;real dr0[ND_ND],dr1[ND_ND], xf[ND_ND], dy;

&nbsp;if (n_uds < N_REQUIRED_UDSI)
&nbsp;&nbsp;Internal_Error("not enough user-defined scalars allocated";

&nbsp;if(n_udm < N_REQUIRED_UDM)
&nbsp;&nbsp;Error("Not enough udm allocated\n";

&nbsp;/* Do nothing if gradient isn't allocated yet. */
&nbsp;if (!Data_Valid_P()) return;
&nbsp;thread_loop_c(t,domain)
&nbsp;{
&nbsp;&nbsp;if (FLUID_THREAD_P(t))
&nbsp;&nbsp;&nbsp;if (NULL != THREAD_STORAGE(t,SV_UDS_I(phi)) && NULL != T_STORAGE_R_NV (t,SV_UDSI_G(phi)))
&nbsp;&nbsp;&nbsp;{
&nbsp;&nbsp;&nbsp;&nbsp;begin_c_loop(c,t)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; /*体循环*/
&nbsp;&nbsp;&nbsp;&nbsp;{
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;C_UDMI(c,t,Ex) = -1.*C_UDSI_G(c,t,phi)[0];&nbsp;&nbsp;&nbsp; /*Ex*/
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;C_UDMI(c,t,Ey) = -1.*C_UDSI_G(c,t,phi)[1];
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;C_UDMI(c,t,Ez) = -1.*C_UDSI_G(c,t,phi)[2];
&nbsp;&nbsp;&nbsp;&nbsp;}
&nbsp;&nbsp;&nbsp;&nbsp;end_c_loop(c,t)
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;/* Message("x electric strength: %g\n", C_UDMI(c,t,Ex)); */
&nbsp;&nbsp;&nbsp;}
&nbsp;}
}

DEFINE_ADJUST(calm_rhos, domain)&nbsp;&nbsp; /* 界面电密度,出问题的地方，请帮忙指正*/
{&nbsp;
&nbsp;Thread *t;
&nbsp;Thread *pri_th,*sec_th;
&nbsp;cell_t c;
&nbsp;real xc[ND_ND], EEp[3], EEs[3], nn[3], EE[3], rr;
&nbsp;/*&nbsp; Domain *domain = Get_Domain(1); Get domain pointer */
&nbsp;int ID = 1;
&nbsp;t =Lookup_Thread(domain, ID); /* mixture-level thread pointer */
&nbsp;pri_th = THREAD_SUB_THREAD(t,0);
&nbsp;sec_th = THREAD_SUB_THREAD(t,1);
&nbsp;
&nbsp;thread_loop_c(t,domain)
&nbsp;{
&nbsp;&nbsp;begin_c_loop(c,t)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; /*体循环*/
&nbsp;&nbsp;{
&nbsp;&nbsp;&nbsp;C_CENTROID(xc,c,t);
&nbsp;&nbsp;&nbsp;if (C_VOF(c,sec_th)>0. && C_VOF(c,sec_th)<1.)
&nbsp;&nbsp;&nbsp;{
&nbsp;&nbsp;&nbsp;&nbsp;EEs[0]=C_UDMI(c,sec_th,Ex);/* 次相电场强度 */
&nbsp;&nbsp;&nbsp;&nbsp;EEs[1]=C_UDMI(c,sec_th,Ey);
&nbsp;&nbsp;&nbsp;&nbsp;EEs[2]=C_UDMI(c,sec_th,Ez);
&nbsp;&nbsp;&nbsp;&nbsp;EEp[0]=C_UDMI(c,pri_th,Ex);/* 基相电场强度 */
&nbsp;&nbsp;&nbsp;&nbsp;EEp[1]=C_UDMI(c,pri_th,Ey);&nbsp;
&nbsp;&nbsp;&nbsp;&nbsp;EEp[2]=C_UDMI(c,pri_th,Ez);
&nbsp;&nbsp;&nbsp;&nbsp;nn[0]=C_UDMI(c,sec_th,vfx);/* vof梯度、法向矢量 */
&nbsp;&nbsp;&nbsp;&nbsp;nn[1]=C_UDMI(c,sec_th,vfy);
&nbsp;&nbsp;&nbsp;&nbsp;nn[2]=C_UDMI(c,sec_th,vfz);
&nbsp;&nbsp;&nbsp;&nbsp;NV_V(EEs,=,EEp);
&nbsp;&nbsp;&nbsp;&nbsp;rr=EPS_0*EPS_1;
&nbsp;&nbsp;&nbsp;&nbsp;NV_VS_VS(EE, =,EEp,*,EPS_0, -, EEs,*,rr);
&nbsp;&nbsp;&nbsp;&nbsp;C_UDMI(c,t,rhos) = NV_DOT(nn,EE);
&nbsp;&nbsp;&nbsp;}
&nbsp;&nbsp;&nbsp;else
&nbsp;&nbsp;&nbsp;{
&nbsp;&nbsp;&nbsp;&nbsp;C_UDMI(c,t,rhos) =0;
&nbsp;&nbsp;&nbsp;}
&nbsp;&nbsp;}
&nbsp;&nbsp;end_c_loop(c,t)
&nbsp;}
}
赞一下 回复此楼
2楼2017-02-26 11:51:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guoyhchina

金虫 (小有名气)
还是乱码，附上附件，最后一个udf DEFINE_ADJUST(calm_rhos, domain)出错。
赞一下 回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : vof_interface_force.c
    • 2017-02-26 12:08:14, 3.23 K
3楼2017-02-26 12:08:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

想你睡不着

新虫 (正式写手)

【答案】应助回帖

你好，我想问下C_UDSI_G()后边的中括号是什么意思，代表什么？
赞一下 回复此楼
4楼2017-03-12 10:12:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卒死于猝

新虫 (初入文坛)
你好，我最近也被如何计算气液界面静电力的问题所困扰，想请问一下这个问题解决了么，能否帮忙指点一下呢，谢谢~
赞一下 回复此楼
5楼2018-07-25 17:48:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 guoyhchina 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 08工科275分求调剂 +13 AaAa7420 2026-03-31 13/650 2026-04-02 06:19 by 286640313
[考研] 301求调剂 +12 A_JiXing 2026-04-01 12/600 2026-04-01 22:47 by peike
[考研] 环境工程297分求调剂一志愿杭高院 +10 GENJIOW 2026-03-31 11/550 2026-04-01 21:38 by GENJIOW
[考研] 286分调剂 +18 Faune 2026-03-30 20/1000 2026-04-01 21:37 by 495374996
[考研] 材料科学与工程调剂 +16 深V宿舍吧 2026-03-30 17/850 2026-04-01 21:00 by 赖春艳
[考研] 261求调剂 +3 明仔· 2026-04-01 3/150 2026-04-01 20:52 by cq2548
[考研] 310分求调剂 +4 成功上岸wang 2026-04-01 4/200 2026-04-01 20:35 by liu823948201
[考研] 一志愿郑大材料工程290求调剂 +18 Youth_ 2026-03-30 18/900 2026-04-01 19:23 by 无际的草原
[考研] 【求调剂】085601材料工程专硕 | 总分272 | +10 脚滑的守法公民 2026-03-27 10/500 2026-04-01 17:23 by pies112
[考研] 311求调剂 +10 李芷新1 2026-03-31 10/500 2026-04-01 14:38 by chenqifeng666
[考研] 0855机械初试280求调剂 +3 kazenotori 2026-03-31 3/150 2026-04-01 10:08 by JourneyLucky
[考研] 生医工0831调剂求推荐 +4 小熊睿睿_s 2026-03-27 6/300 2026-03-31 17:12 by 记事本2026
[基金申请] 面上5B能上会吗？ +8 redcom 2026-03-29 8/400 2026-03-31 15:53 by niuailing
[考研] 286求调剂 +6 Faune 2026-03-30 6/300 2026-03-31 14:37 by jp9609
[考研] 272求调剂，接受跨专业调剂！ +3 闲鱼卢 2026-03-31 3/150 2026-03-31 13:00 by 替代品000
[考研] 085701环境工程求调剂 +11 多久上课 2026-03-27 12/600 2026-03-30 21:21 by 研究僧导导
[考研] 279求调剂 +12 j的立方 2026-03-29 12/600 2026-03-30 20:30 by dick_runner
[考研] 322求调剂 +10 宋明欣 2026-03-27 10/500 2026-03-30 18:47 by 544594351
[考研] 297求调剂 +17 田洪有 2026-03-26 18/900 2026-03-30 18:32 by nothing投稿中
[考研] 356求调剂 +3 gysy?s?a 2026-03-28 3/150 2026-03-29 00:33 by 544594351
信息提示
关闭
请填处理意见
关闭