 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 1204 | 回復(fù): 2 | ||
westlifexp銅蟲 (小有名氣)
|
[求助]
歐拉模型--udf自定義粘度問題
|
|
先說下模型，氣液固三相歐拉模型，其中我的固體相可以看做是液體，主相為空氣，次相第一個為粘結(jié)劑，另一個為粉末。非穩(wěn)態(tài)的，速度入口，壓力出口。 因為材料屬性中粘結(jié)劑和粉末的粘度隨時間時刻發(fā)生變化，并且都是關(guān)于混合相的溫度和速度的一個函數(shù)，通過提取混合相的溫度和速度來對固體和液體的粘度進行計算。 問題：主要是關(guān)于指針的選取，不是很清楚。 我把兩段程序貼出來，麻煩有哪位幫忙看看，所編的程序能否表達意思，謝謝。 第一個是粘結(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; } 導(dǎo)入fluent，編譯成功了，但是在后來初始化的過程出現(xiàn)了問題。 |
回復(fù)此樓» 猜你喜歡
求調(diào)劑
已經(jīng)有3人回復(fù)
-
265求調(diào)劑
已經(jīng)有4人回復(fù)
-
085700資源與環(huán)境308求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑
已經(jīng)有12人回復(fù)
-
286分人工智能專業(yè)請求調(diào)劑愿意跨考！
已經(jīng)有3人回復(fù)
-
329求調(diào)劑
已經(jīng)有5人回復(fù)
-
申請回稿延期一個月，編輯同意了。但系統(tǒng)上的時間沒變，給編輯又寫郵件了，沒回復(fù)
已經(jīng)有4人回復(fù)
-
材料學(xué)碩318求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑
已經(jīng)有6人回復(fù)
-
081700化工學(xué)碩調(diào)劑
已經(jīng)有3人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 請教UDF中歐拉模型+反應(yīng)
已經(jīng)有20人回復(fù)
westlifexp
銅蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 72.9
- 散金: 3
- 帖子: 82
- 在線: 43.8小時
- 蟲號: 1450881
- 注冊: 2011-10-19
- 性別: GG
- 專業(yè): 塑料
2樓2012-10-09 10:26:18
玲玲0328
新蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 6564.5
- 帖子: 711
- 在線: 30.4小時
- 蟲號: 2796028
- 注冊: 2013-11-12
- 專業(yè): 工程熱物理與能源利用
3樓2018-06-25 21:03:02
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 westlifexp 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +5 | 1孫悟空 2026-03-17 | 6/300 |
|
|---|---|---|---|---|
|
|
[考博] 申博26年 +3 | 八6八68 2026-03-19 | 3/150 |
|
|
|
[考研] 085601材料工程專碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
|
[考研] 材料與化工求調(diào)劑 +7 | 為學(xué)666 2026-03-16 | 7/350 |
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-18 | 3/150 |
|
|
|
[考研] 求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
|
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +4 | 生物工程調(diào)劑 2026-03-16 | 12/600 |
|
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +5 | 想上岸的鯉魚 2026-03-18 | 6/300 |
|
|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +7 | 步川酷紫123 2026-03-13 | 7/350 |
|
|
|
[考研] 311求調(diào)劑 +6 | 26研0 2026-03-15 | 6/300 |
|
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
|
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
|
[考研] 生物學(xué)071000 329分求調(diào)劑 +3 | 我愛生物生物愛?/a> 2026-03-17 | 3/150 |
|
|
|
[考研] 278求調(diào)劑 +5 | 煙火先于春 2026-03-17 | 5/250 |
|
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
|
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項國獎兩個實用型專利一項省級立項 +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
|
[考研] 274求調(diào)劑 +5 | 時間點 2026-03-13 | 5/250 |
|
|
|
[考研] 297一志愿上交085600求調(diào)劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|
|
|
[考研] 266求調(diào)劑 +4 | 學(xué)員97LZgn 2026-03-13 | 4/200 |
|