| 3 | 1/1 | 返回列表 |
| 查看: 1046 | 回復(fù): 2 | ||
akqinmian金蟲 (小有名氣)
|
[求助]
fluent單機多核并行版本變化物性參數(shù)UDF求助
|
|
求教各位大俠,小弟寫的隨溫度變化的材料參數(shù)和發(fā)熱源項的UDF ,由于遠(yuǎn)程使用別人的工作站(沒有裝VC),單機單核用解釋的方式導(dǎo)入沒有問題,單機多核情況下總是報錯,要不就是迭代開始就卡死,沒反應(yīng)。 求教各位,如果如何改寫為并行版本? 小弟計算機底子不行,用戶說明關(guān)于這部分的內(nèi)容反復(fù)讀也不太理解,幫助里也沒有這樣的例子。 #include "udf.h" DEFINE_PROPERTY(vis, cell, thread) { real temperature, nd; temperature = C_T(cell, thread); { if (temperature >=553.15 && temperature<=618.15) nd=-0.0000000000116472*pow(temperature,3)+0.0000000203453*pow(temperature,2)-0.0000122063*temperature+0.00259719; else if (temperature >618.15 && temperature<=632.15) nd=-0.000000000069519*pow(temperature,3)+0.000000127294*pow(temperature,2)-0.0000780956*temperature+0.0161298; else if (temperature >632.15 && temperature<=648.15) nd=-0.000000000454378*pow(temperature,3)+0.000000861208*pow(temperature,2)-0.00054462*temperature+0.114982; else if (temperature>648.15 && temperature<=656.55) nd=-0.0000000130651*pow(temperature,3)+0.0000254502*pow(temperature,2)-0.0165263*temperature+3.57745; else if (temperature>656.55 && temperature<=662.55) nd=-0.0000000110287*pow(temperature,3)+0.0000221444*pow(temperature,2)-0.01482*temperature+3.30587; else if (temperature>662.55 && temperature<=667.95) nd=-0.0000000039237*pow(temperature,3)+0.00000786908*pow(temperature,2)-0.00526079*temperature+1.17243; else if (temperature>667.95 && temperature<=674.35) nd=-0.000000000571417*pow(temperature,3)+0.00000116031*pow(temperature,2)-0.000785453*temperature+0.177281; else if (temperature>674.35 && temperature<=683.15) nd=-0.000000000132212*pow(temperature,3)+0.0000002727*pow(temperature,2)-0.000187513*temperature+0.0430127; else if (temperature>683.15 && temperature<=693.15) nd=-0.0000000000362821*pow(temperature,3)+0.0000000763507*pow(temperature,2)-0.0000535492*temperature+0.0125458; else if (temperature>693.15 && temperature<=765.15) nd=-0.00000000000256198*pow(temperature,3)+0.00000000579822*pow(temperature,2)-0.00000433921*temperature+0.00110355; else if (temperature>765.15 && temperature<=823.15) nd=-0.000000000000214632*pow(temperature,3)+0.000000000530378*pow(temperature,2)-0.000000398162*temperature+0.000120633; else nd=0.0; } return nd; #include "udf.h" DEFINE_SOURCE(cell_y_source1,c, t, dS, eqn) { real source; real x[ND_ND]; real z; C_CENTROID(x, c, t); z=x[2]; if(z>=0 && z<0.08) source = 1.96617E+08; else if(z>=0.08 && z<0.24) source = 1.94500E+08; else if(z>=0.24 && z<0.40) source = 1.89151E+08; else if(z>=0.40 && z<0.56) source = 2.10398E+08; else if(z>=0.56 && z<0.72) source = 2.43178E+08; else if(z>=0.72 && z<0.88) source = 3.19465E+08; else if(z>=0.88 && z<1.04) source = 3.93218E+08; else if(z>=1.04 && z<1.20) source = 4.65868E+08; else if(z>=1.20 && z<1.36) source = 5.30883E+08; else if(z>=1.36 && z<1.52) source = 5.63514E+08; else if(z>=1.52 && z<1.68) source = 5.75473E+08; else if(z>=1.68 && z<1.84) source = 5.54165E+08; else if(z>=1.84 && z<2.00) source = 5.08562E+08; else if(z>=2.00 && z<2.16) source = 4.70529E+08; else if(z>=2.16 && z<2.32) source = 4.20242E+08; else if(z>=2.32 && z<2.48) source = 3.56882E+08; else if(z>=2.48 && z<2.64) source = 3.11648E+08; else if(z>=2.64 && z<2.80) source = 2.77573E+08; else if(z>=2.80 && z<2.96) source = 2.51141E+08; else if(z>=2.96 && z<3.12) source = 2.29009E+08; else if(z>=3.12 && z<3.28) source = 2.16346E+08; else if(z>=3.28 && z<3.44) source = 2.07810E+08; else if(z>=3.44 && z<3.60) source = 2.13757E+08; else if(z>=3.60 && z<3.68) source = 2.23299E+08; else source =0.0; dS[eqn] = 0.0; return source; } |
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 280求調(diào)劑 +11 | 咕嚕曉曉 2026-03-18 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 考研調(diào)劑 +3 | 呼呼?~+123456 2026-03-21 | 3/150 |
|
|
[基金申請] 山東省面上項目限額評審 +3 | 石瑞0426 2026-03-19 | 3/150 |
|
|
[考研] 313求調(diào)劑 +4 | 肆叁貳壹22 2026-03-19 | 4/200 |
|
|
[考研] 296求調(diào)劑 +4 | www_q 2026-03-20 | 4/200 |
|
|
[考研] 機械專碩299求調(diào)劑至材料 +3 | kkcoco25 2026-03-16 | 4/200 |
|
|
[考研] 華東師范大學(xué)-071000生物學(xué)-293分-求調(diào)劑 +3 | 研究生何瑤明 2026-03-18 | 3/150 |
|
|
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +9 | Doleres 2026-03-19 | 9/450 |
|
|
[考研] 317求調(diào)劑 +5 | 申子申申 2026-03-19 | 9/450 |
|
|
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑,總分315(英一) +5 | sbdksD 2026-03-19 | 5/250 |
|
|
[考研] A區(qū)線材料學(xué)調(diào)劑 +5 | 周周無極 2026-03-20 | 5/250 |
|
|
[考研] 295材料求調(diào)劑,一志愿武漢理工085601專碩 +5 | Charlieyq 2026-03-19 | 5/250 |
|
|
[考研] 281求調(diào)劑(0805) +14 | 煙汐憶海 2026-03-16 | 25/1250 |
|
|
[考研] 求調(diào)劑 +3 | 暗涌afhb 2026-03-16 | 3/150 |
|
|
[考博] 申博26年 +3 | 八6八68 2026-03-19 | 3/150 |
|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
[考博] 26博士申請 +3 | 1042136743 2026-03-17 | 3/150 |
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
[考研] 一志愿南京大學(xué),080500材料科學(xué)與工程,調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|