CFD Online Discussion Forums - heat generation inside battery

Code:
#include "udf.h"

#include "math.h"



DEFINE_ADJUST(cellheatgenerate,d)

{

        real SOC,ds,tm,Ri,tempu,volume,source,temp;

        Thread *t;

        cell_t c;

        tm = CURRENT_TIME;

        thread_loop_c(t,d)

        {        

                begin_c_loop(c,t)

                {                

                        SOC = 1.-((90.*tm)/5400.);



                        if ((SOC >= 0.) && (SOC <= 0.77))

                                ds = (99.88*SOC)-76.67;

                        else if ((SOC > 0.77) && (SOC <= 0.87))

                                ds = -30.;

                        else if ((SOC >0.87) && (SOC <=1.))

                                ds = -20.;



                        volume = C_VOLUME(c,t);

                        tempu = C_T(c,t);

                        if (tm == 0.)

                                temp = 291.15;

                        else if (tm > 0.)

                                temp = tempu;

                        

                        if (temp <= 298.15)

                                Ri = 0.000002258*pow(SOC,-0.3952);

                        else if ((temp > 298.15) && (temp <= 308.15))

                                Ri = 0.000001857*pow(SOC,-0.2787);

                        else if (temp > 308.15)

                                Ri = 0.000001659*pow(SOC,-0.1692);



                        temp = C_T(c,t);

                        source = ((Ri*8100.)-(90.*temp*ds/96.485))/volume;

                        C_UDMI(c,t,0) = source;                        

                }                

                end_c_loop(c,t)                

        }        

}



DEFINE_SOURCE(my_heat_source,c,t,dS,eqn)

{

        real source;

        source =  C_UDMI(c,t,0);        

        dS[eqn] = 0.0;

        return source;

}