24小时热门版块排行榜    

查看: 1328  |  回复: 2
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

westlifexp

铜虫 (小有名气)

[求助] 欧拉模型--udf自定义粘度问题

先说下模型,气液固三相欧拉模型,其中我的固体相可以看做是液体,主相为空气,次相第一个为粘结剂,另一个为粉末。非稳态的,速度入口,压力出口。
因为材料属性中粘结剂和粉末的粘度随时间时刻发生变化,并且都是关于混合相的温度和速度的一个函数,通过提取混合相的温度和速度来对固体和液体的粘度进行计算。
   问题:主要是关于指针的选取,不是很清楚。
   我把两段程序贴出来,麻烦有哪位帮忙看看,所编的程序能否表达意思,谢谢。
   
   第一个是粘结剂,第二个是粉末的:
#include "udf.h"
#include

DEFINE_PROPERTY(binder_viscosity,cell,thread)

{

  real viscosity_0,viscosity,binder_viscosity,rate,temperture;

  real dudx,dudy,dudz,dvdx,dvdy,dvdz,dwdx,dwdy,dwdz;

  Thread *mixture_thread;

  Thread *thread_gas;

  Thread *thread_binder;

  Thread *thread_powder;

  thread_gas = THREAD_SUB_THREAD(mixture_thread, 2);      

  thread_binder = THREAD_SUB_THREAD(mixture_thread, 3);        

  thread_powder = THREAD_SUB_THREAD(mixture_thread, 4);
  
  temperture = C_T(cell,mixture_thread);

  dudx=C_DUDX(cell,mixture_thread);

  dudy=C_DUDY(cell,mixture_thread);

  dudz=C_DUDZ(cell,mixture_thread);

  dvdx=C_DVDX(cell,mixture_thread);

  dvdy=C_DVDY(cell,mixture_thread);

  dvdz=C_DVDZ(cell,mixture_thread);

  dwdx=C_DWDX(cell,mixture_thread);

  dwdy=C_DWDY(cell,mixture_thread);

  dwdz=C_DWDZ(cell,mixture_thread);

  rate =sqrt(0.5*(4*dudx*dudx+2*(dudy+dvdx)*(dudy+dvdx)+2*(dudz+dwdx)*
  
  (dudz+dwdx)+4*dvdy*dvdy+2*(dvdz+dwdy)*(dvdz+dwdy)+4*dwdz*dwdz));

  viscosity_0 = 9.539e-03*exp(5683/temperture);

  viscosity = viscosity_0/(1+4.648e-04*pow(viscosity_0*rate,0.6717));   

  binder_viscosity = 0.93 * viscosity;        

  return binder_viscosity;

}



#include "udf.h"
#include

DEFINE_PROPERTY(powder_viscosity,cell,thread)

{

  real viscosity_0,viscosity,powder_viscosity,rate,temperture;

  real dudx,dudy,dudz,dvdx,dvdy,dvdz,dwdx,dwdy,dwdz;

  Thread *mixture_thread;

  Thread *thread_gas;

  Thread *thread_binder;

  Thread *thread_powder;

  thread_gas = THREAD_SUB_THREAD(mixture_thread, 2);      

  thread_binder = THREAD_SUB_THREAD(mixture_thread, 3);        

  thread_powder = THREAD_SUB_THREAD(mixture_thread, 4);
  
  temperture = C_T(cell,mixture_thread);

  dudx=C_DUDX(cell,mixture_thread);

  dudy=C_DUDY(cell,mixture_thread);

  dudz=C_DUDZ(cell,mixture_thread);

  dvdx=C_DVDX(cell,mixture_thread);

  dvdy=C_DVDY(cell,mixture_thread);

  dvdz=C_DVDZ(cell,mixture_thread);

  dwdx=C_DWDX(cell,mixture_thread);

  dwdy=C_DWDY(cell,mixture_thread);

  dwdz=C_DWDZ(cell,mixture_thread);

  rate =sqrt(0.5*(4*dudx*dudx+2*(dudy+dvdx)*(dudy+dvdx)+2*(dudz+dwdx)*
  
  (dudz+dwdx)+4*dvdy*dvdy+2*(dvdz+dwdy)*(dvdz+dwdy)+4*dwdz*dwdz));

  viscosity_0 = 9.539e-03*exp(5683/temperture);

  viscosity = viscosity_0/(1+4.648e-04*pow(viscosity_0*rate,0.6717));   

  powder_viscosity = 1.33 * viscosity;        

  return powder_viscosity;

}

    导入fluent,编译成功了,但是在后来初始化的过程出现了问题。
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

fluent模拟粉末注塑成型
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

玲玲0328

新虫 (正式写手)

楼主解决问题了吗?我现在写UDF也遇到指针问题,如果楼主解决了可以请教一下吗r

发自小木虫IOS客户端
3楼2018-06-25 21:03:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +8 Kittylucky 2026-08-27 8/400 2026-08-27 18:33 by raiy123
[基金申请] 申请删除本帖 +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 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 怎么查啊 +5 huang1991js 2026-08-26 5/250 2026-08-26 16:09 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 宝贝虫子
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +4 yuleib84 2026-08-26 5/250 2026-08-26 13:40 by yuleib84
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
请填处理意见