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

24小时热门版块排行榜    

微信图片_20250911175443_873_253.png
小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » fluent udf DEFINE_ADJUST求两相界面力能load但不能计算
51/1返回列表
查看: 1698  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭