| 查看: 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，编译成功了，但是在后来初始化的过程出现了问题。 |
回复此楼» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有17人回复
-
面上合作单位盖章
已经有5人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有8人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有7人回复
-
申请删除本帖
已经有7人回复
-
基金不中，共勉
已经有10人回复
-
国自然面上复盘~欢迎讨论
已经有11人回复
-
基金未中，这种答复是模板吗？
已经有6人回复
-
看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉……
已经有6人回复
-
怎么看青基中了没有啊
已经有5人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教UDF中欧拉模型+反应
已经有20人回复
3楼2018-06-25 21:03:02
10