24小时热门版块排行榜

返回列表

westlifexp

铜虫 (小有名气)

应助: 1 (幼儿园)
金币: 72.9
散金: 3
帖子: 82
在线: 43.8小时
虫号: 1450881
注册: 2011-10-19
性别: GG
专业: 塑料

[求助] 欧拉模型--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，编译成功了，但是在后来初始化的过程出现了问题。

回复此楼