# UDF for Temperature Dependent Energy Source

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

 LinkBack Thread Tools Display Modes
 December 10, 2012, 02:20 UDF for Temperature Dependent Energy Source #1 New Member   Mukesh Join Date: Sep 2012 Posts: 28 Rep Power: 7 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

 January 11, 2013, 05:28 #2 Senior Member   Join Date: Nov 2009 Posts: 125 Rep Power: 9 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

 January 18, 2013, 05:32 problem on UDF returned parameters #3 Senior Member   Join Date: Nov 2009 Posts: 125 Rep Power: 9 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

 January 18, 2013, 10:49 #4 New Member   Mukesh Join Date: Sep 2012 Posts: 28 Rep Power: 7 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

 January 18, 2013, 11:56 #5 Senior Member   Join Date: Nov 2009 Posts: 125 Rep Power: 9 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

January 19, 2013, 06:07
#6
New Member

Mukesh
Join Date: Sep 2012
Posts: 28
Rep Power: 7
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 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?

 January 19, 2013, 09:31 #7 Senior Member   Join Date: Nov 2009 Posts: 125 Rep Power: 9 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

 September 15, 2017, 10:27 #8 New Member   maziyar ghanikolahloo Join Date: Jan 2017 Posts: 5 Rep Power: 2 Dear machtech001, I want to write a UDF in which there is a constant term and also a negative term which is related to temperature. I am new to UDF how can I accomplish that? I want to use Average Temp of the solid zone and source term should be used just in solid zone. #include "udf.h" DEFINE_SOURCE(source_term,cell,thread,dS,eqn) { real source; real temp=C_T(cell,thread); int a=0.5546; int b=1.84*10^-4; int I=700; int SCR=1000; int A=373.83; source=261682243-(a*(1-b)*temp)*I*SCR*A return source }

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post mksca OpenFOAM Programming & Development 13 August 18, 2014 22:26 kevinz FLUENT 0 December 12, 2011 17:45 Muhammad Aamir Khan FLUENT 1 July 25, 2007 02:08 Nitesh Main CFD Forum 0 July 17, 2006 13:19 Zhong Lei Main CFD Forum 23 May 14, 1999 13:22

All times are GMT -4. The time now is 04:32.

 Contact Us - CFD Online - Privacy Statement - Top