 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 1203 | 回復: 2 | ||
westlifexp銅蟲 (小有名氣)
|
[求助]
歐拉模型--udf自定義粘度問題
|
|
先說下模型，氣液固三相歐拉模型，其中我的固體相可以看做是液體，主相為空氣，次相第一個為粘結(jié)劑，另一個為粉末。非穩(wěn)態(tài)的，速度入口，壓力出口。 因為材料屬性中粘結(jié)劑和粉末的粘度隨時間時刻發(fā)生變化，并且都是關于混合相的溫度和速度的一個函數(shù)，通過提取混合相的溫度和速度來對固體和液體的粘度進行計算。 問題：主要是關于指針的選取，不是很清楚。 我把兩段程序貼出來，麻煩有哪位幫忙看看，所編的程序能否表達意思，謝謝。 第一個是粘結(jié)劑，第二個是粉末的： #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，編譯成功了，但是在后來初始化的過程出現(xiàn)了問題。 |
回復此樓» 猜你喜歡
0856調(diào)劑，是學校就去
已經(jīng)有7人回復
-
0817 化學工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章
已經(jīng)有14人回復
-
申博26年
已經(jīng)有3人回復
-
一志愿吉林大學材料學碩321求調(diào)劑
已經(jīng)有10人回復
-
0703化學調(diào)劑 ，六級已過，有科研經(jīng)歷
已經(jīng)有12人回復
-
294求調(diào)劑材料與化工專碩
已經(jīng)有13人回復
-
288求調(diào)劑
已經(jīng)有14人回復
-
307求調(diào)劑
已經(jīng)有8人回復
-
一志愿南京理工大學085701資源與環(huán)境302分求調(diào)劑
已經(jīng)有5人回復
-
288求調(diào)劑，一志愿華南理工大學071005
已經(jīng)有5人回復
» 本主題相關商家推薦: (我也要在這里推廣)
» 本主題相關價值貼推薦，對您同樣有幫助:
- 請教UDF中歐拉模型+反應
已經(jīng)有20人回復
westlifexp
銅蟲 (小有名氣)
- 應助: 1 (幼兒園)
- 金幣: 72.9
- 散金: 3
- 帖子: 82
- 在線: 43.8小時
- 蟲號: 1450881
- 注冊: 2011-10-19
- 性別: GG
- 專業(yè): 塑料
2樓2012-10-09 10:26:18
玲玲0328
新蟲 (正式寫手)
- 應助: 0 (幼兒園)
- 金幣: 6564.5
- 帖子: 711
- 在線: 30.4小時
- 蟲號: 2796028
- 注冊: 2013-11-12
- 專業(yè): 工程熱物理與能源利用
3樓2018-06-25 21:03:02
相關版塊跳轉(zhuǎn)
我要訂閱樓主 westlifexp 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 288求調(diào)劑，一志愿華南理工大學071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 材料專碩274一志愿陜西師范大學求調(diào)劑 +8 | 薛云鵬 2026-03-13 | 8/400 |
|
|
|
[考研] 一志愿北京化工大學0703化學318分，有科研經(jīng)歷，求調(diào)劑 +3 | 一瓶苯甲酸 2026-03-14 | 3/150 |
|
|
|
[考研] 一志愿天津大學化學工藝專業(yè)（081702）315分求調(diào)劑 +11 | yangfz 2026-03-17 | 11/550 |
|
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-18 | 3/150 |
|
|
|
[考研] 0703化學調(diào)劑 +5 | pupcoco 2026-03-17 | 8/400 |
|
|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|
|
|
[考研] 330求調(diào)劑 +3 | 小材化本科 2026-03-18 | 3/150 |
|
|
|
[考研] 311求調(diào)劑 +4 | 冬十三 2026-03-18 | 4/200 |
|
|
|
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +3 | 生物工程調(diào)劑 2026-03-17 | 7/350 |
|
|
|
[考研] 0854可跨調(diào)劑，一作一項核心論文五項專利，省、國級證書40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
|
[考研] 0703化學求調(diào)劑 總分331 +3 | ZY-05 2026-03-13 | 3/150 |
|
|
|
[考研] 0703化學調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
|
[考研] 工科材料085601 279求調(diào)劑 +6 | 困于星晨 2026-03-17 | 6/300 |
|
|
|
[考研] 26考研求調(diào)劑 +6 | 丶宏Sir 2026-03-13 | 6/300 |
|
|
|
[考研] 274求調(diào)劑 +5 | 時間點 2026-03-13 | 5/250 |
|
|
|
[考研] 326求調(diào)劑 +4 | 諾貝爾化學獎覬?/a> 2026-03-15 | 7/350 |
|
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
|
[考研] 22408總分284求調(diào)劑 +3 | InAspic 2026-03-13 | 3/150 |
|
|
|
[考研] 本科南京大學一志愿川大藥學327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|