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

UDF for Temperature Dependent Energy Source

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

Reply
 
LinkBack Thread Tools Display Modes
Old   December 10, 2012, 02:20
Default UDF for Temperature Dependent Energy Source
  #1
New Member
 
Mukesh
Join Date: Sep 2012
Posts: 28
Rep Power: 4
er.mkumar is on a distinguished road
Hello All,

I am solving a problem wherein there is an energy source term with two parts:
1) Part dependent on Temperature [i.e. T*(Cp_g - Cp_s) ]
2) Part independent of Temperature

Can anybody suggest how to include temperature dependent source term in the UDF as we know Energy equation solves for Temperature.

Thanks a lot
Mukesh
er.mkumar is offline   Reply With Quote

Old   January 11, 2013, 05:28
Default
  #2
Senior Member
 
Join Date: Nov 2009
Posts: 125
Rep Power: 7
mactech001 is on a distinguished road
Hi Mukesh,

did you get any response or help in this please? i'm also having the same problem.

i was thinking using DEFINE_SOURCE and C_T macro, but not yet sure how to relate them.
__________________
Thank you for your kind attention.

Kind regards,
mactech001
Currently using: ANSYS v13
mactech001 is offline   Reply With Quote

Old   January 18, 2013, 05:32
Default problem on UDF returned parameters
  #3
Senior Member
 
Join Date: Nov 2009
Posts: 125
Rep Power: 7
mactech001 is on a distinguished road
Dear all,

i've the following UDF code that attempts to apply a temperature-dependent source term to a solid cell.

#include "udf.h"
#include "mem.h"

DEFINE_SOURCE(cell_x_source,cell,thread,dS,eqn)
{
real source;
int Ro=12e-3; /*12mOhm at 20C*/
int To=20; /*temperature 20C*/
int Iphase=20; /*Arms*/

source=Iphase*Iphase*Ro*((C_T(cell,thread)-273)/To)/C_VOLUME(cell,thread);
Message("C_T output: %g\n",C_T(cell,thread));
Message("source output: %g\n",source);
Message("C_VOLUME output: %g\n",C_VOLUME(cell,thread));

return source;
}


While running 100 iterations, i observe in the Message output that the C_VOLUME of the solid cell is changing rather than keeping at a constant value as expected. why is this happening please?

The final outcome of the run was:
Error: Divergence detected in AMG solver: temperature
Error Object: #f

Hope to receive any comments/ideas you may have.
__________________
Thank you for your kind attention.

Kind regards,
mactech001
Currently using: ANSYS v13
mactech001 is offline   Reply With Quote

Old   January 18, 2013, 10:49
Default
  #4
New Member
 
Mukesh
Join Date: Sep 2012
Posts: 28
Rep Power: 4
er.mkumar is on a distinguished road
Hello,

Though I'm not sure if we could use the MACRO C_T(c,t) for SOURCE term in Energy equation, still I suggest you remove that C_VOLUME(c,t) as we must define source term as per unit volume to correctly solve the Energy eqn in FLUENT.

Hope that works.

All The Best,
Mukesh
er.mkumar is offline   Reply With Quote

Old   January 18, 2013, 11:56
Default
  #5
Senior Member
 
Join Date: Nov 2009
Posts: 125
Rep Power: 7
mactech001 is on a distinguished road
Dear Mukesh,
thank you for your response.

The C_VOLUME in my equation is to obtain the volume of the solid cell, and use it to divide my energy/power in order to express the source in per volume. Note that in the source equation, i've divided by C_VOLUME.

My only query is: Why the C_VOLUME value is changing from one iteration to the next?
__________________
Thank you for your kind attention.

Kind regards,
mactech001
Currently using: ANSYS v13
mactech001 is offline   Reply With Quote

Old   January 19, 2013, 06:07
Default
  #6
New Member
 
Mukesh
Join Date: Sep 2012
Posts: 28
Rep Power: 4
er.mkumar is on a distinguished road
Could you please send in your energy equation so as to check the units. FLUENT solves equations with source in Unit/volume terms only. So please check the units.

Best of Luck......Mukesh

Quote:
Originally Posted by mactech001 View Post
Dear Mukesh,
thank you for your response.

The C_VOLUME in my equation is to obtain the volume of the solid cell, and use it to divide my energy/power in order to express the source in per volume. Note that in the source equation, i've divided by C_VOLUME.

My only query is: Why the C_VOLUME value is changing from one iteration to the next?
er.mkumar is offline   Reply With Quote

Old   January 19, 2013, 09:31
Default
  #7
Senior Member
 
Join Date: Nov 2009
Posts: 125
Rep Power: 7
mactech001 is on a distinguished road
Dear Mukesh,

In my application, my source term:

source=Iphase*Iphase*Ro*((C_T(cell,thread)-273)/To)/C_VOLUME(cell,thread);

where,
Iphase*Iphase*Ro*((C_T(cell,thread)-273)/To) is in units of Watts

C_VOLUME(cell,thread) is in units of cubic meters;
__________________
Thank you for your kind attention.

Kind regards,
mactech001
Currently using: ANSYS v13
mactech001 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
wmake compiling new solver mksca OpenFOAM Programming & Development 13 August 18, 2014 22:26
time dependent energy source term kevinz FLUENT 0 December 12, 2011 17:45
UDF for energy source term Muhammad Aamir Khan FLUENT 1 July 25, 2007 02:08
Time dependent Energy source Nitesh Main CFD Forum 0 July 17, 2006 13:19
Why FVM for high-Re flows? Zhong Lei Main CFD Forum 23 May 14, 1999 13:22


All times are GMT -4. The time now is 06:54.