版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門(mén)版塊排行榜

返回列表

chang.z.s

銅蟲(chóng) (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 263.8
散金: 136
紅花: 2
帖子: 119
在線: 142.6小時(shí)
蟲(chóng)號(hào): 756128
注冊(cè): 2009-04-23
性別: GG
專業(yè): 可再生與替代能源利用中的

[求助] 用udf實(shí)現(xiàn)coupled邊界固體區(qū)域中溫度最值提取

物理問(wèn)題是想提取耦合邊界固體區(qū)域中心面上溫度最大值曲線和最小值取消（非穩(wěn)態(tài)）。monitor可實(shí)現(xiàn)監(jiān)控某一點(diǎn)或某面的平均值隨時(shí)間變化，但想得到的點(diǎn)比較多就考慮用udf，但一直沒(méi)成功過(guò)。
主要疑問(wèn)在于：
d=Get_Domain(1); 對(duì)單相流體獲得的是所有流體區(qū)域，邊界的固體區(qū)域如何獲得？
我在udf加了 t = Lookup_Thread(d,12); （12是固體壁面區(qū)域在boundary conditon 中的zone id）來(lái)獲得固體域的thread，但udf沒(méi)成功，是否是因?yàn)樵赿是流體域不能得到固體的thread。如果這樣的話，如何得到固體域呢？

回復(fù)此樓

» 猜你喜歡

一志愿中南化學(xué)337求調(diào)劑已經(jīng)有4人回復(fù)
281求調(diào)劑（0805）已經(jīng)有23人回復(fù)
265求調(diào)劑已經(jīng)有9人回復(fù)
279分求調(diào)劑一志愿211 已經(jīng)有10人回復(fù)
招收調(diào)劑碩士已經(jīng)有12人回復(fù)
296求調(diào)劑已經(jīng)有8人回復(fù)
招收博士1-2人已經(jīng)有3人回復(fù)
085700資源與環(huán)境308求調(diào)劑已經(jīng)有11人回復(fù)
274求調(diào)劑已經(jīng)有8人回復(fù)
290求調(diào)劑已經(jīng)有4人回復(fù)

1樓 2013-10-17 15:38:53

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

chang.z.s

銅蟲(chóng) (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 263.8
散金: 136
紅花: 2
帖子: 119
在線: 142.6小時(shí)
蟲(chóng)號(hào): 756128
注冊(cè): 2009-04-23
性別: GG
專業(yè): 可再生與替代能源利用中的

自頂一下把代碼貼出來(lái)，寫(xiě)的很挫請(qǐng)高手指教

DEFINE_EXECUTE_AT_END(excutend_twminmax)
{
Domain *d;
/* declare domain pointer since it is not passed as an argument to the DEFINE macro */
Thread *t;
cell_t c;
real temptf,ctime1;
int i;
FILE *twminmax;
real x[ND_ND];
real twmin0[100]={0}; real twmax0[100]={0}; real twmin[100]={0}; real twmax[100]={0}; //####grid is 500####only for this mesh;
real times0[100]={0}; real timeb0[100]={0}; real times[100]={0}; real timeb[100]={0};
i=0;
d=Get_Domain(1); // all fluid zone
t = Lookup_Thread(d,12);// ########pay attention to this part ,id of insert####### zone_ID 5
ctime1 = RP_Get_Real("flow-time"

;

for(i=0;i<250;i++)
{
thread_loop_c(t,d)
/* Compute max, min, volume-averaged temperature */
/* Loop over all cells */
{
begin_c_loop(c,t)
{
temptf = C_T(c,t);
C_CENTROID(x,c,t);

if ((fabs(x[1]-1.04)<0.02)&&(fabs(x[0]-(i+1)*0.02)<0.01))
{

if (ctime1<=1)
{
twmin=temptf; // temperature at current time
twmax=temptf;
times=ctime1;
timeb=ctime1;
}

if ((twmin0<twmin)&&(ctime1>1))// initial state skip these two
{
twmin=twmin0;
times=times0;
}

if ((twmax0>twmax)&&(ctime1>1))
{
twmax=twmax0;
timeb=timeb0;
}

}
twmin0=twmin; // 0 temperature at last time step, save the last time step data
twmax0=twmax;
times0=times;
timeb0=timeb;

}
end_c_loop(c,t)
}

}
/* write the data of twmin adn twmax */
twminmax=fopen("twminmax.txt","w";
for(i=0;i<250;i++)
{
fprintf(twminmax,"%g %g %g %g %g\n", ctime1, twmin, times, twmax, timeb);
}
Message("twmin[100];%g\n",twmin[100]);
fclose(twminmax);
}

贊一下

回復(fù)此樓

2樓2013-10-18 15:41:29

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

stryhood

至尊木蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 4 (幼兒園)
金幣: 15071.3
散金: 16
紅花: 1
帖子: 1397
在線: 385.9小時(shí)
蟲(chóng)號(hào): 1417580
注冊(cè): 2011-09-26
專業(yè): 可再生與替代能源利用中的

【答案】應(yīng)助回帖

★
臭水溝: 金幣+1, 謝謝交流~ 2013-10-25 20:30:17

loop thread和loop cell都用，你限定一下固體區(qū)域的坐標(biāo)，應(yīng)該能行

贊一下

回復(fù)此樓

3樓2013-10-21 08:58:34

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

chang.z.s

銅蟲(chóng) (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 263.8
散金: 136
紅花: 2
帖子: 119
在線: 142.6小時(shí)
蟲(chóng)號(hào): 756128
注冊(cè): 2009-04-23
性別: GG
專業(yè): 可再生與替代能源利用中的

引用回帖:

3樓: Originally posted by stryhood at 2013-10-21 08:58:34
loop thread和loop cell都用，你限定一下固體區(qū)域的坐標(biāo)，應(yīng)該能行

udf里是這樣處理的同時(shí)thread和cell的循環(huán)都有找固體區(qū)域的某個(gè)位置，
但最后輸出都是初值

thread_loop_c(t,d)
/* Compute max, min, volume-averaged temperature */
/* Loop over all cells */
{
begin_c_loop(c,t)

不確定是不是 if ((fabs(x[1]-1.04)<0.02)&&(fabs(x[0]-(i+1)*0.02)<0.01)) 判斷太嚴(yán)格，根本沒(méi)執(zhí)行，這種可能也不大。

贊一下

回復(fù)此樓

4樓2013-10-22 16:16:54

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

stryhood

至尊木蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 4 (幼兒園)
金幣: 15071.3
散金: 16
紅花: 1
帖子: 1397
在線: 385.9小時(shí)
蟲(chóng)號(hào): 1417580
注冊(cè): 2011-09-26
專業(yè): 可再生與替代能源利用中的

★
臭水溝: 金幣+1, 謝謝交流~ 2013-10-25 20:30:40

引用回帖:

4樓: Originally posted by chang.z.s at 2013-10-22 16:16:54
udf里是這樣處理的同時(shí)thread和cell的循環(huán)都有找固體區(qū)域的某個(gè)位置，
但最后輸出都是初值

thread_loop_c(t,d)
/* Compute max, min, volume-averaged temperature */
/* Loop over all cells */
{
...

thread和cell loop之后，應(yīng)該是全部循環(huán)的，看了你的代碼應(yīng)該是不用for循環(huán)的

贊一下

回復(fù)此樓

5樓2013-10-22 20:17:09

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

chang.z.s

銅蟲(chóng) (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 263.8
散金: 136
紅花: 2
帖子: 119
在線: 142.6小時(shí)
蟲(chóng)號(hào): 756128
注冊(cè): 2009-04-23
性別: GG
專業(yè): 可再生與替代能源利用中的

引用回帖:

5樓: Originally posted by stryhood at 2013-10-22 20:17:09
thread和cell loop之后，應(yīng)該是全部循環(huán)的，看了你的代碼應(yīng)該是不用for循環(huán)的...

for 是為了得到多個(gè)（250）點(diǎn)在時(shí)間進(jìn)程上的最大值和最小值，這個(gè)地方是比較低效，全域循環(huán)只找出一個(gè)點(diǎn)。

贊一下

回復(fù)此樓

6樓2013-10-25 15:52:22

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

stryhood

至尊木蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 4 (幼兒園)
金幣: 15071.3
散金: 16
紅花: 1
帖子: 1397
在線: 385.9小時(shí)
蟲(chóng)號(hào): 1417580
注冊(cè): 2011-09-26
專業(yè): 可再生與替代能源利用中的

引用回帖:

6樓: Originally posted by chang.z.s at 2013-10-25 15:52:22
for 是為了得到多個(gè)（250）點(diǎn)在時(shí)間進(jìn)程上的最大值和最小值，這個(gè)地方是比較低效，全域循環(huán)只找出一個(gè)點(diǎn)。...

不要用for循環(huán)，這個(gè)我有經(jīng)驗(yàn)，你一個(gè)一個(gè)的話能算，如果用for的話就不行，即使是for循環(huán)只有3個(gè)，原因我也搞不清楚。

贊一下

回復(fù)此樓

7樓2013-10-26 11:29:45

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zhaolei539

鐵蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 324.1
紅花: 1
帖子: 76
在線: 21.1小時(shí)
蟲(chóng)號(hào): 2870773
注冊(cè): 2013-12-14
性別: MM
專業(yè): 制造系統(tǒng)與自動(dòng)化

引用回帖:

7樓: Originally posted by stryhood at 2013-10-26 11:29:45
不要用for循環(huán)，這個(gè)我有經(jīng)驗(yàn)，你一個(gè)一個(gè)的話能算，如果用for的話就不行，即使是for循環(huán)只有3個(gè)，原因我也搞不清楚。...

你好！想請(qǐng)教提取某個(gè)邊界面的溫度，賦值給另一個(gè)邊界面
面上不同位置溫度不一樣不知道UDM怎么存儲(chǔ)的？我這樣寫(xiě)行不行？
#include "udf.h"
DEFINE_PROFILE(hnt_temperature,thread)
{
face_t f;
real tem_liquid;
real tem_solid;
Thread *thread_liquid=Lookup_Thread(domain,14); /*14是交界面液體側(cè)面的id*/;
Thread *thread_solid=Lookup_Thread(domain,15); /*15是交界面固體側(cè)面的id*/;

begin_f_loop(f,thread_liquid)
{
temp=F_T(f,thread_liquid);
F_UDMI(f,t,0)=temp_liquid;

}
  end_f_loop(f,thread_liquid)

begin_f_loop(f,thread_solid)
{
tem_solid=F_UDMI(f,t,0);
F_PROFILE(f,thread_solid)=tem_solid;

}
  end_f_loop(f,thread_solid)
}
麻煩您幫忙看看謝謝！

贊一下

回復(fù)此樓

努力！

8樓2014-06-14 10:33:45

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zhaolei539

鐵蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 324.1
紅花: 1
帖子: 76
在線: 21.1小時(shí)
蟲(chóng)號(hào): 2870773
注冊(cè): 2013-12-14
性別: MM
專業(yè): 制造系統(tǒng)與自動(dòng)化

引用回帖:

贊一下

回復(fù)此樓

努力！

9樓2014-06-14 15:29:12

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 chang.z.s 的主題更新

返回列表

普通表情龍兔虎貓

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 295復(fù)試調(diào)劑 +6	簡(jiǎn)木ChuFront 2026-03-19	6/300	2026-03-20 11:32 by lature00
[考研] 08工學(xué)調(diào)劑 +4	用戶573181 2026-03-20	4/200	2026-03-20 10:32 by 朗月清風(fēng)1
[考研] 0703化學(xué)調(diào)劑，六級(jí)已過(guò)，有科研經(jīng)歷 +12	曦熙兮 2026-03-15	12/600	2026-03-19 19:42 by maocaozhuxi
[考研] 321求調(diào)劑 +8	何潤(rùn)采123 2026-03-18	10/500	2026-03-19 16:46 by 何潤(rùn)采123
[考研] 0703化學(xué)調(diào)劑 +4	18889395102 2026-03-18	4/200	2026-03-19 16:13 by 30660438
[考研] 085600材料與化工求調(diào)劑 +6	緒幸與子 2026-03-17	6/300	2026-03-19 13:27 by houyaoxu
[考研] 354求調(diào)劑 +4	Tyoumou 2026-03-18	7/350	2026-03-18 21:45 by Tyoumou
[考研] 295求調(diào)劑 +3	一志愿京區(qū)211 2026-03-18	5/250	2026-03-18 17:03 by zhaoqian0518
[考研] 311求調(diào)劑 +11	冬十三 2026-03-15	12/600	2026-03-18 14:36 by 星空星月
[考研] 297求調(diào)劑 +8	戲精丹丹丹 2026-03-17	8/400	2026-03-18 14:30 by laoshidan
[考研] 085601求調(diào)劑 +4	Du.11 2026-03-16	4/200	2026-03-17 17:08 by ruiyingmiao
[考研] 332求調(diào)劑 +6	Zz版 2026-03-13	6/300	2026-03-17 17:03 by ruiyingmiao
[考研] 有沒(méi)有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3	TqlXswl 2026-03-16	7/350	2026-03-17 15:23 by TqlXswl
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4	Jy? 2026-03-16	4/200	2026-03-17 11:02 by gaoqiong
[考研] 藥學(xué)383 求調(diào)劑 +3	藥學(xué)chy 2026-03-15	4/200	2026-03-16 20:51 by 元子^0^
[考研] 326求調(diào)劑 +4	諾貝爾化學(xué)獎(jiǎng)覬?/a> 2026-03-15	7/350	2026-03-16 17:11 by 諾貝爾化學(xué)獎(jiǎng)覬?/a>
[考研] 22408總分284求調(diào)劑 +3	InAspic 2026-03-13	3/150	2026-03-15 11:10 by zhq0425
[考研] 中科大材料與化工319求調(diào)劑 +3	孟鑫材料 2026-03-14	3/150	2026-03-14 20:10 by ms629
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3	麥田耕者 2026-03-14	3/150	2026-03-14 20:04 by 外星文明
[考研] 復(fù)試調(diào)劑 +3	呼呼？~+123456 2026-03-14	3/150	2026-03-14 16:53 by WTUChen