/***********************************************************************/
/* heat_conductivity.c */
/* UDF for specifying a temperature-dependent heat_conductivity property of gap gas*/
/***********************************************************************/
/***********************************************************************/
/* heat_source.c */
/* UDF for specifying a location-dependent heat_source of fuel */
/***********************************************************************/
DEFINE_SOURCE(cell_heat_source, cell, thread)
{
real x[ND_ND];
real heat_source;
real t,z;
t= RP_Get_Real("flow-time"
z = x[2];
heat_source = qmax * sin(pi*(z+0.1)/1.1) * sin (0.2 * pi * t);
return heat_source;
}