|
[Sponsors] |
March 7, 2011, 14:33 |
SOLVED:Udf function average temp inlet - outlet
|
#1 |
New Member
Join Date: Apr 2009
Location: Mexico
Posts: 7
Rep Power: 17 |
Hello,
I would appreciate if somebody take a look in the UDF and give me some advices. I am trying to get the average temperature at the outlet of a pipe, and after some manipulations use it in the pressure inlet of the next time step. It compiles and the simulation runs up to the first time step. It starts with reversed flow then just after the first time step prints Temperature limited to 1 in the domain 1. Tavg prints -nan. Thanks DEFINE_EXECUTE_AT_END(Temp_inlet) { Domain *d; real Temp_t; real area; real area_tot; real tavg; real A[ND_ND]; Thread *t; face_t f; d=Get_Domain(1); t=Lookup_Thread(d,11); begin_f_loop(f,t) { area = F_AREA(A,f,t); Temp_t += F_T(f,t); area_tot += area; tavg += Temp_t*area; } end_f_loop(f,t) tavg /= area_tot; tavg_tot1=tavg; printf("\n Tavg = %f\n",tavg_tot1); /*output -nan*/ } DEFINE_PROFILE(Temp_inlet,t,i) { face_t f; real time; real current_time; time=CURRENT_TIMESTEP; current_time=CURRENT_TIME; if (current_time >=10.) { begin_f_loop(f,t) { F_PROFILE(f,t,i)=((fabs(flow_tot)*(tavg_tot-tavg_tot1)*time)/15.)+tavg_tot1; } end_f_loop(f,t) } else { begin_f_loop(f,t) { F_PROFILE(f,t,i)=288.16; } end_f_loop(f,t) } } Last edited by erodv; March 7, 2011 at 19:56. Reason: SOLVED |
|
March 7, 2011, 15:08 |
|
#2 | |
Senior Member
|
Hi,
I think the problem is in your face loop, substitute this and check. hope that help you: Quote:
|
||
March 7, 2011, 19:54 |
Solved
|
#3 |
New Member
Join Date: Apr 2009
Location: Mexico
Posts: 7
Rep Power: 17 |
Thanks Amir, I did your proposed changes and it worked, this thread has been solved .
|
|
November 14, 2016, 09:29 |
|
#4 |
New Member
alper yıldırım
Join Date: Feb 2014
Posts: 10
Rep Power: 12 |
||
October 6, 2017, 11:16 |
inlet average temperature udf
|
#5 |
New Member
sina sadighi
Join Date: Jul 2017
Posts: 5
Rep Power: 9 |
hi
this is my code , when i put outlet id in it , it works . but when i put inlet id in it , it doesn't work . pls help what should i do . #include "udf.h" DEFINE_ADJUST(temp_ave, domain) { int ID = 9 ; FILE *output; Thread *thread=Lookup_Thread(domain, ID); real total_temp=0; real ave_T, T; int i=0; face_t face; begin_f_loop(face, thread) { T = F_T(face, thread); total_temp = total_temp + T; i = i+1; } end_f_loop(face, thread) ave_T = total_temp / i ; printf("inlet average of temp : %f\n" , ave_T); output = fopen("inlet average of temp.txt","a"); fprintf(output, "inlet average of temp: %f\n", ave_T); fclose(output); } |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
unsteady inlet velocity udf | jill | FLUENT | 8 | February 10, 2015 06:04 |
Inlet but not outlet | mireis | FLUENT | 1 | June 10, 2010 10:56 |
Error with Wmake | skabilan | OpenFOAM Installation | 3 | July 28, 2009 00:35 |
Question about inlet and outlet BC | babyface | Main CFD Forum | 1 | January 15, 2003 09:17 |
UDF in Fluent to Match Mass Flow at Pressure Outlet | Jonas Larsson | Main CFD Forum | 1 | April 29, 1999 10:44 |