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

How to deal latent heat during phase change through source term udf

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

Like Tree2Likes
  • 1 Post By yangbinsunshine
  • 1 Post By rampal

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 21, 2012, 10:00
Default How to deal latent heat during phase change through source term udf
  #1
New Member
 
deepak
Join Date: Sep 2011
Location: india
Posts: 8
Rep Power: 14
dkb_bg is on a distinguished road
Hi all
i have written udf for phase change from liquid to vapour in a narrow pipe.i think there is some problem with energy udf, since temperature is not decreasing but it should decrease.

DEFINE_SOURCE(enrg_src, cell, mix_th, dS, eqn)
{
Thread *pri_th;
Thread *sec_th;
real m_dot;
real P_SAT ;
real T_SAT;
P_SAT=exp(A+(B/(C_T(cell,mix_th)+C))+D*C_T(cell,mix_th)+E*log(C_T (cell,mix_th)))*1e5;
T_SAT= -P_SAT*P_SAT*2e-11+P_SAT*7e-5+254.1;

pri_th = THREAD_SUB_THREAD(mix_th, 0);
sec_th = THREAD_SUB_THREAD(mix_th, 1);

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

return LAT_HT*m_dot;

i am not very much clear with this udf. means due to phase change temp of liquid should drop.

plz enlighten me t will be grateful.
dkb_bg is offline   Reply With Quote

Old   November 14, 2013, 01:50
Default
  #2
New Member
 
yangbin
Join Date: Nov 2013
Location: haerbin
Posts: 4
Rep Power: 12
yangbinsunshine is on a distinguished road
Quote:
Originally Posted by dkb_bg View Post
Hi all
i have written udf for phase change from liquid to vapour in a narrow pipe.i think there is some problem with energy udf, since temperature is not decreasing but it should decrease.

DEFINE_SOURCE(enrg_src, cell, mix_th, dS, eqn)
{
Thread *pri_th;
Thread *sec_th;
real m_dot;
real P_SAT ;
real T_SAT;
P_SAT=exp(A+(B/(C_T(cell,mix_th)+C))+D*C_T(cell,mix_th)+E*log(C_T (cell,mix_th)))*1e5;
T_SAT= -P_SAT*P_SAT*2e-11+P_SAT*7e-5+254.1;

pri_th = THREAD_SUB_THREAD(mix_th, 0);
sec_th = THREAD_SUB_THREAD(mix_th, 1);

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

return LAT_HT*m_dot;

i am not very much clear with this udf. means due to phase change temp of liquid should drop.

plz enlighten me t will be grateful.
hi, i am very intrested ijn your problem.
i want to know how the coefficient 0.1 comes from?
is it suitable for every conditions?
SACHIN1639 likes this.
yangbinsunshine is offline   Reply With Quote

Old   November 18, 2013, 01:35
Default
  #3
New Member
 
yangbin
Join Date: Nov 2013
Location: haerbin
Posts: 4
Rep Power: 12
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

Old   January 8, 2016, 01:29
Default phase change problem inside heat pipe
  #4
Member
 
Ram Kumar Pal
Join Date: Apr 2015
Posts: 38
Rep Power: 11
rampal is on a distinguished road
Dear friend, I'm doing the simulation of Heat pipe in FLUENT as my project work. I have written udf for mass and energy source for phase change inside the device. But condensation is not happening in condenser zone. If you have completed this problem, please help me.
SHUJAN likes this.
rampal is offline   Reply With Quote

Old   April 12, 2018, 20:26
Default
  #5
New Member
 
England
Join Date: Apr 2018
Posts: 1
Rep Power: 0
chrischris is on a distinguished road
Quote:
Originally Posted by rampal View Post
Dear friend, I'm doing the simulation of Heat pipe in FLUENT as my project work. I have written udf for mass and energy source for phase change inside the device. But condensation is not happening in condenser zone. If you have completed this problem, please help me.
Same problem. How did you insert the udf ? Did you compile them ?
Could a VOF + EvaporationCondensation(Lee) model work without UDF ?
How about I just add an energy UDF on the water and vapor cell zone conditions ?
chrischris is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
friction forces icoFoam ofslcm OpenFOAM 3 April 7, 2012 10:57
About phase change heat and mass transfer Michael FLUENT 2 February 13, 2011 01:49
OpenFOAM on MinGW crosscompiler hosted on Linux allenzhao OpenFOAM Installation 127 January 30, 2009 19:08
udf with moving heat source Ryan FLUENT 0 April 10, 2003 19:13
help needed about phase change Yanhu Guo Main CFD Forum 4 January 23, 2001 23:16


All times are GMT -4. The time now is 03:07.