| 查看: 1356 | 回复: 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，编译成功了，但是在后来初始化的过程出现了问题。 |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教UDF中欧拉模型+反应
已经有20人回复
westlifexp
铜虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 72.9
- 散金: 3
- 帖子: 82
- 在线: 43.8小时
- 虫号: 1450881
- 注册: 2011-10-19
- 性别: GG
- 专业: 塑料
2楼2012-10-09 10:26:18
3楼2018-06-25 21:03:02
