You should be use Heat flux input for moving heat source. Firstly you find heat flux by the help of moving componen

#include "udf.h"

DEFINE_PROFILE(Qfield, thread, position)

{

real coord[ND_ND];

real x,y,z; // coordinates

real xC=0.0075,yC=0.0075,zC=0.005; // initial position of heater

real D = 7e-5; // the heater diameter (m)

real V = 0.05; // the heater moving velocity (m/s)

real P = 5.0; // (W)

real i = 1;

real h = D*1;

real t; // time

real t0 = 1e-3; // the time that the heater is stationary

face_t f; // system parameter

//real RigionX = D*10;

real RigionX = 0.0075; // region (heating face) size (m)

real R,Rk,q0,r2;

R = D/2;

q0 = (2*P)/(3.14159268*(R*R));

t = RP_Get_Real("flow-time"); //

// heater path

if(t < t0)

{

xC = xC;

//yC = yC+(i*h);

//zC = zC;

yC = yC;

zC = zC;

}

else

{

t = t-t0;

while( (t*V) > 0.0075 )

{

t = t - (0.0075/V);

}

xC = xC+V*t;

yC = yC+V*t;

zC = zC;

}

//yC = yC+(i*h);

//xC = xC+(i*h);

//while( (t*V) > RigionX )

//while( (t*V) > RigionX )

// t = t - RigionX/V;

//i=i+1;

// Specifying the shape and position of the heater spot at a certain point in time

// iteration over the heating side

begin_f_loop(f,thread)

{

F_CENTROID(coord,f,thread);

x = coord[0];

y = coord[2];

z = coord[2];

r2 = sqrt(((x-xC)*(x-xC))+((y-yC)*(y-yC)));

F_PROFILE(f,thread,position) = 0.3*(((2*P)/(3.14159268*(R*R)))*exp((-2*(r2*r2))/(R*R)));

// F_PROFILE(f,thread,position) =(2*P/((3.1415926)*(R*R)))*exp(-2*(((r2*r2))/((R*R)))))

}

end_f_loop(f,thread)

}

follow. good question

Should I insert that code as a command block? and if so where would I put it?