CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Fluent UDF and Scheme Programming

UDF for phase change

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree1Likes
  • 1 Post By N400

Reply
 
LinkBack Thread Tools Display Modes
Old   December 23, 2012, 21:46
Default UDF for phase change
  #1
New Member
 
nowbelle
Join Date: Dec 2012
Posts: 3
Rep Power: 4
N400 is on a distinguished road
I want to simulate the phase tranfer in the heat pipe, and a udf was compiled in FLUENT,which is as follow:

#include "udf.h"
#include "sg_mphase.h"
#define T_SAT 323
#define LAT_HT 2356e3
DEFINE_SOURCE(vap_src, cell, pri_th, dS, eqn)
{
Thread *mix_th, *sec_th;
real m_dot_v;
mix_th = THREAD_SUPER_THREAD(pri_th);
sec_th = THREAD_SUB_THREAD(mix_th, 1);
/*m_dot_v=0.;*/
if(C_T(cell, mix_th)>=T_SAT)
{
m_dot_v = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*fabs(C_T(cell, sec_th) - T_SAT)/T_SAT;

dS[eqn] = 0.;
}

if(C_T(cell, mix_th)<=T_SAT)
{
m_dot_v = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)*fabs(T_SAT-C_T(cell,pri_th))/T_SAT;
dS[eqn] = -0.1*C_R(cell, pri_th)*fabs(C_T(cell, pri_th) - T_SAT)/T_SAT;
}
return m_dot_v ;
}

DEFINE_SOURCE(liq_src, cell, sec_th, dS, eqn)
{
Thread * mix_th, *pri_th;
real m_dot_l;
mix_th = THREAD_SUPER_THREAD(sec_th);
pri_th = THREAD_SUB_THREAD(mix_th, 0);
/*m_dot_l=0.;*/
if(C_T(cell, mix_th)>=T_SAT)
{
m_dot_l = -0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*fabs(C_T(cell, sec_th) - T_SAT)/T_SAT;
dS[eqn] = -0.1*C_R(cell, sec_th)*fabs(C_T(cell, sec_th) - T_SAT)/T_SAT;
}
if(C_T(cell, mix_th)<=T_SAT)
{
m_dot_l = 0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)*fabs(T_SAT-C_T(cell,pri_th))/T_SAT;
dS[eqn] =0.;
}

return m_dot_l;
}

DEFINE_SOURCE(enrg_src, cell, mix_th, dS, eqn)
{
Thread *pri_th, *sec_th;
real m_dot;
pri_th = THREAD_SUB_THREAD(mix_th, 0);
sec_th = THREAD_SUB_THREAD(mix_th, 1);
if(C_T(cell, mix_th)>=T_SAT)
{
m_dot = -0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*fabs(C_T(cell, sec_th) - T_SAT)/T_SAT;
dS[eqn] = -0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*LAT_HT/T_SAT;
}

if(C_T(cell, mix_th)<=T_SAT)
{
m_dot = 0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)*fabs(T_SAT-C_T(cell,pri_th))/T_SAT;
dS[eqn] = 0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)*LAT_HT/T_SAT;
}
return LAT_HT*m_dot; }

Using the code above, I can simulate the evaporation, but the condensation did not appear, even the temperature of vapor was much lower than the saturation temperature.
Now I am so confused. Anyone can help me?
N400 is offline   Reply With Quote

Old   December 24, 2012, 00:35
Default
  #2
Senior Member
 
SSL
Join Date: Oct 2012
Posts: 227
Rep Power: 5
msaeedsadeghi is on a distinguished road
You may have not defined boundary conditions correctly. Please share files to check.
msaeedsadeghi is offline   Reply With Quote

Old   December 24, 2012, 01:09
Default
  #3
New Member
 
nowbelle
Join Date: Dec 2012
Posts: 3
Rep Power: 4
N400 is on a distinguished road
Quote:
Originally Posted by msaeedsadeghi View Post
You may have not defined boundary conditions correctly. Please share files to check.
My case and data files are bigger than the attachment allowing, may I send them to you via email?
My email address is nowbelle@yeah.net
Thank you !

PS:
Anyone else who is intersting in my problem or have the relevant experiences please contact me, we can communicate directly and quickly via emails! ^_^
sircorp likes this.

Last edited by N400; December 24, 2012 at 03:14.
N400 is offline   Reply With Quote

Old   December 24, 2012, 05:09
Default
  #4
Senior Member
 
SSL
Join Date: Oct 2012
Posts: 227
Rep Power: 5
msaeedsadeghi is on a distinguished road
Yes. Please email me.
msaeedsadeghi is offline   Reply With Quote

Old   December 24, 2012, 05:40
Default
  #5
New Member
 
nowbelle
Join Date: Dec 2012
Posts: 3
Rep Power: 4
N400 is on a distinguished road
Thanks!
I just sent you the files via email, and l am looking forward your reply!
N400 is offline   Reply With Quote

Old   April 19, 2013, 04:51
Default
  #6
New Member
 
Join Date: Apr 2013
Posts: 1
Rep Power: 0
Shengchao Hu is on a distinguished road
Hi,N400
I am writing the udf for the phase change of LNG,i have some questions to share with you.Could you send me your files ?
Thanks very much! My email adress:shengchaohu2012@163.com
Shengchao Hu is offline   Reply With Quote

Old   November 18, 2013, 02:33
Default
  #7
New Member
 
yangbin
Join Date: Nov 2013
Location: haerbin
Posts: 4
Rep Power: 3
yangbinsunshine is on a distinguished road
Hi,I am working on this udf recently.and I am confused about this sentence"m_dot_v = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*fabs(C_T(cell, sec_th) - T_SAT)/T_SAT".
First,in this sentence,the coefficent is 0.1,but i do not know why.
Second,what is the unit of "m_dot_v"?
Thank you very much,and any advice will be highly appreciated.
yangbinsunshine is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
UDF to change Rotation Speed in a MRF simulation Mike FLUENT 3 September 27, 2011 06:46
UDF Cmu change selçuk ataş FLUENT 1 September 29, 2008 14:16
UDF to change BC from Wall to Grid Line Aaron Smith FLUENT 0 February 21, 2006 10:14
Change Courant number from UDF HS FLUENT 0 December 1, 2005 06:16
UDF...UDF...UDF...UDF Luc SEMINEL FLUENT 0 November 25, 2002 05:03


All times are GMT -4. The time now is 01:13.