[axt5488]

Hello Everyone,

I am working on to create a volumetric heat source (VHS) with a gaussian profile inside a simple cubic domain. I am working on Fluent 'Solidification & Melting' model. In order to create a moving volumetric heat source-VHS I am using DEFINE_SOURCE for my UDF. The equation that I am using for gaussian profile for VHS is attached [https://doi-org.ezproxy.uta.edu/10.1...0-019-04908-3]



However I have energy source equation as a function of time and position the energy source that I created did not move. Here you can see my UDF basically below

#include "udf.h"

#include "math.h"

#include "metric.h"

#include "sg.h"

#include "sg_mphase.h"

#include "mem.h"

#include "sg_mem.h"

#include "flow.h"

#include "unsteady.h"

#define A 0.24 /*Absorption coefficient*/

#define P 400 /*Laser Power*/

#define re 0.0001 /*radius on top*/

#define ri 0.0002 /*radius on bottom*/

#define zi 0.00005 /*cone length parameter lower z axis*/

#define ze 0.0001 /*cone length parameter upper z axis*/

#define v1 1.200 /*laser velocity*/

#define v2 0 /*laser velocity*/

#define pi 3.1415

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

{

real x[ND_ND];

C_CENTROID(x, c, t);

real x0 = 0.00061, y0 = 0, z0 = 0; /*initial position of the heat flux*/

real t0 = 1; /*the time that the heat flux is stationary*/

real time = CURRENT_TIME;

real H = ze - zi;

real r0;

r0 = ri; /*r0 = re + ((ri - re) / H) * (z - ze) z-ze equals to H then r0 equals to ri*/

real d;

real Q0;

double Q = 0.0;

if (time<= 1)

{

Q = 0.0;

}

else

{

d = M_PI * H * pow(re, 2) + (re * ri) + pow(ri, 2);

Q0 = 6 * A * P / d;

Q = Q0 * exp(-1 * pow(x[0] - x0 - v1 * time, 2.0) + pow(x[1] - y0, 2.0) / pow(r0, 2));

dS[eqn] = 0.0;

}

return Q;

}

The script compiled without any error but the heat source did not move. Do I need to add more comments for advective flux for scalar transport equation.
This is one of the topics previously discussed for welding/AM applications in the forum but it seems it is not solved yet. I hope this time we will figure it out.
Any comments are welcome and I will be more than happy if anybody in the forum knows how to UDF a moving heat source.

Thanks in advance

Updated by Ahmet Alptug
Metadata
Views 1
Last Activity less than a minute ago
Tags
Tags ·Add

[AlexanderZ]

coordinate of heat source location must be updated with time to make it move
as of know I can see that only magnitude of source is time dependent